Search

Search Results

• 1. [Article] A 3-D Tomographic Model of the P-wave Velocity Structure of the Central Cascadia Forearc

  Signals from the R/V Langseth’s tuned airgun array were recorded on an array of 33 EarthScope FlexArray seismometers deployed in the Oregon Coast Range from 43.5° to 45°N as part of the 2012 Ridge2Trench ...

  Citation × Citation Citation Abstract Copy Title: A 3-D Tomographic Model of the P-wave Velocity Structure of the Central Cascadia Forearc Author: Kenyon, Christopher Bradford Signals from the R/V Langseth’s tuned airgun array were recorded on an array of 33 EarthScope FlexArray seismometers deployed in the Oregon Coast Range from 43.5° to 45°N as part of the 2012 Ridge2Trench experiment to image the structure of the Juan de Fuca plate. This segment of the Juan de Fuca/North America plate boundary slipped in 2 moderate low-angle thrust earthquakes in 2004, and continuing seismicity has clustered around these events since that time. Previously acquired bathymetric, potential field, and seismic data indicate the presence of subducted seamounts in this region and hinted at large velocity variations of the overlying forearc crust immediately above the plate boundary. Approximately 18000 airgun shots along 4 profiles on the continental shelf and upper slope and a profile oblique to the coast that extends across the margin to the trench provided sources for this study. Data quality is excellent, with strong Pg and PmP arrivals observed at most stations for most shots. We present velocity models that merge the new data with constraints from onshore/offshore data acquired in 1989 and 1996 to improve the resolution and extend the velocity model for this segment of the Cascadia subduction zone, which appears to be a segment boundary in the paleoseismic record and currently exhibits anomalous locking behavior as indicated by geodetic data. Our approach to modeling involved starting with a 1D P-wave velocity model and using an inversion process to modify the model to beetter fit the observed first arrival times of the airgun shots. In a subsequent step, we used a 2D starting model that incorporated constraints from prior 2D experiments that extended farther seaward and landward of the 3D study region. Finally, a smoothed version of the model obtained from the 2D inversion was used as a starting model to see if the final model could be further refined. Our final preferred model used a 2D starting model and moderate smoothing. The RMS misfit of predicted travel times is reduced from 0.48s for the initial 2D model to 0.040s for the final preferred 3D model. This model includes interesting velocity structure including possible subducted seamounts on the subducting Juan de Fuca plate and a mafic dike related to Cape Perpetua volcanics within the upper plate. These results contribute to the broader goal of understanding the relationship between plate boundary behavior and crustal structure. Title: A 3-D Tomographic Model of the P-wave Velocity Structure of the Central Cascadia Forearc Author: Kenyon, Christopher Bradford Signals from the R/V Langseth’s tuned airgun array were recorded on an array of 33 EarthScope FlexArray seismometers deployed in the Oregon Coast Range from 43.5° to 45°N as part of the 2012 Ridge2Trench experiment to image the structure of the Juan de Fuca plate. This segment of the Juan de Fuca/North America plate boundary slipped in 2 moderate low-angle thrust earthquakes in 2004, and continuing seismicity has clustered around these events since that time. Previously acquired bathymetric, potential field, and seismic data indicate the presence of subducted seamounts in this region and hinted at large velocity variations of the overlying forearc crust immediately above the plate boundary. Approximately 18000 airgun shots along 4 profiles on the continental shelf and upper slope and a profile oblique to the coast that extends across the margin to the trench provided sources for this study. Data quality is excellent, with strong Pg and PmP arrivals observed at most stations for most shots. We present velocity models that merge the new data with constraints from onshore/offshore data acquired in 1989 and 1996 to improve the resolution and extend the velocity model for this segment of the Cascadia subduction zone, which appears to be a segment boundary in the paleoseismic record and currently exhibits anomalous locking behavior as indicated by geodetic data. Our approach to modeling involved starting with a 1D P-wave velocity model and using an inversion process to modify the model to beetter fit the observed first arrival times of the airgun shots. In a subsequent step, we used a 2D starting model that incorporated constraints from prior 2D experiments that extended farther seaward and landward of the 3D study region. Finally, a smoothed version of the model obtained from the 2D inversion was used as a starting model to see if the final model could be further refined. Our final preferred model used a 2D starting model and moderate smoothing. The RMS misfit of predicted travel times is reduced from 0.48s for the initial 2D model to 0.040s for the final preferred 3D model. This model includes interesting velocity structure including possible subducted seamounts on the subducting Juan de Fuca plate and a mafic dike related to Cape Perpetua volcanics within the upper plate. These results contribute to the broader goal of understanding the relationship between plate boundary behavior and crustal structure. Close

• 2. [Article] Evolution and inheritance of Cascadia sub-arc mantle reservoirs

  Inheritance from pre-existing mantle domains and fluid and melt contributions from active subduction together produce the geochemical signatures of mantle-derived arc basalts. In this context, this work ...

  Citation × Citation Citation Abstract Copy Title: Evolution and inheritance of Cascadia sub-arc mantle reservoirs Author: Bromley, Sarah Ashley Inheritance from pre-existing mantle domains and fluid and melt contributions from active subduction together produce the geochemical signatures of mantle-derived arc basalts. In this context, this work evaluates the evolution of Cascadia mantle sources by documenting the isotopic and compositional characteristics of primitive basalts along a transect across the Eocene-Oligocene Proto-Cascadia (EOPC) arc at ~44.5-45.5° N. Primitive EOPC flows, dikes, and sills are exposed across a ~300 km transect that includes the Oregon Coast Range in the Cascadia forearc, the Western Cascades, flanking the modern arc, and the John Day and Eastern Clarno formations east of the Cascades. Like the modern arc, EOPC was built upon accreted terranes of western North America and within the Columbia embayment, which is lithosphere of oceanic affinity that crops out as the Siletzia terrane in the forearc and extends beneath the arc to the backarc. Potential mantle source reservoirs for EOPC magmas include contributions from mantle domains related to pre-existing underlying terranes, distinct North America lithosphere, and depleted Pacific-like upper mantle. In addition, the geochemical characteristics of EOPC magmas have likely been overprinted by subduction processes. Major, trace element, and isotopic data from the EOPC reveal a heterogeneous mantle source that was variably influenced by subduction processes. In the forearc, the high field strength (HFSE) enriched basalts of the Oregon Coast Range represent low degree partial melts of a relatively enriched mantle source. Despite this enriched character, there is little evidence in the trace element or isotopic data to suggest that the OCR (Oregon Coast Range) samples have been strongly influenced by either a crustal or subduction component. Their distinctive ²⁰⁶Pb/²⁰⁴Pb enrichment (as compared to ²⁰⁷Pb/²⁰⁴Pb) distinguishes the forearc magmas from the arc and backarc magmas and also from a hypothetical Cascadia subduction component. Forearc EOPC magmas share a mantle source with the accreted Siletzia terrane, as evidenced by their shared ²⁰⁶Pb/²⁰⁴Pb enrichment. At the apparent arc axis, the Western Cascades produced a diversity of primitive magmas that are, for the most part, interpreted to result from higher degrees of partial melting of a less enriched source. Fluid fluxing appears to have facilitated mantle melting beneath the Western Cascades. Additionally, the mantle beneath the arc may be slightly influenced by the ²⁰⁶Pb/²⁰⁴Pb enriched source underlying the forearc. Though this effect is difficult to resolve, the Western Cascades samples appear to be slightly more variable in ²⁰⁶Pb/²⁰⁴Pb (and generally higher in ²⁰⁶Pb/²⁰⁴Pb at a given ²⁰⁷Pb/²⁰⁴Pb) than their backarc counterparts, suggesting such an influence. Both trace element and isotopic data suggest a significant subduction contribution to EOPC backarc primitive magmas. The EOPC backarc magmas appear to have originated from a heterogeneous mantle and are variably influenced by a subduction component. Though this subduction influence could be an inherited feature, the affinity between John Day and Eastern Clarno magmas and a modern Cascadia sediment source suggest that this is at least a relatively young feature. Taken together, these results demonstrate the heterogeneity of source and process across the Eocene-Oligocene arc. Apparent subduction contribution increases from the geographic forearc to backarc, perhaps indicating a wider arc than is typically envisioned at the time. The distinct ²⁰⁶Pb/²⁰⁴Pb enrichment associated with the Siletzia mantle and apparent in the forearc appears to wane in influence to the east, having at most a minor influence on the Western Cascade arc magmas. As with the Siletzia mantle domain in EOPC arc and backarc, mantle domains associated with the Western Cascades and John Day/Eastern Clarno magmas do not persist in the High Cascades or High Lava Plains. Beneath and behind the arc, mantle reservoirs appear to have been largely replaced since the Eocene-Oligocene. However, the High Lava Plains basalts appear to carry a paleo-enrichment signature that may be an inherited feature from advected mantle. If this enriched mantle was advected from beneath North America, it may explain the observation that High Lava Plains basalts are more enriched to the east. Title: Evolution and inheritance of Cascadia sub-arc mantle reservoirs Author: Bromley, Sarah Ashley Inheritance from pre-existing mantle domains and fluid and melt contributions from active subduction together produce the geochemical signatures of mantle-derived arc basalts. In this context, this work evaluates the evolution of Cascadia mantle sources by documenting the isotopic and compositional characteristics of primitive basalts along a transect across the Eocene-Oligocene Proto-Cascadia (EOPC) arc at ~44.5-45.5° N. Primitive EOPC flows, dikes, and sills are exposed across a ~300 km transect that includes the Oregon Coast Range in the Cascadia forearc, the Western Cascades, flanking the modern arc, and the John Day and Eastern Clarno formations east of the Cascades. Like the modern arc, EOPC was built upon accreted terranes of western North America and within the Columbia embayment, which is lithosphere of oceanic affinity that crops out as the Siletzia terrane in the forearc and extends beneath the arc to the backarc. Potential mantle source reservoirs for EOPC magmas include contributions from mantle domains related to pre-existing underlying terranes, distinct North America lithosphere, and depleted Pacific-like upper mantle. In addition, the geochemical characteristics of EOPC magmas have likely been overprinted by subduction processes. Major, trace element, and isotopic data from the EOPC reveal a heterogeneous mantle source that was variably influenced by subduction processes. In the forearc, the high field strength (HFSE) enriched basalts of the Oregon Coast Range represent low degree partial melts of a relatively enriched mantle source. Despite this enriched character, there is little evidence in the trace element or isotopic data to suggest that the OCR (Oregon Coast Range) samples have been strongly influenced by either a crustal or subduction component. Their distinctive ²⁰⁶Pb/²⁰⁴Pb enrichment (as compared to ²⁰⁷Pb/²⁰⁴Pb) distinguishes the forearc magmas from the arc and backarc magmas and also from a hypothetical Cascadia subduction component. Forearc EOPC magmas share a mantle source with the accreted Siletzia terrane, as evidenced by their shared ²⁰⁶Pb/²⁰⁴Pb enrichment. At the apparent arc axis, the Western Cascades produced a diversity of primitive magmas that are, for the most part, interpreted to result from higher degrees of partial melting of a less enriched source. Fluid fluxing appears to have facilitated mantle melting beneath the Western Cascades. Additionally, the mantle beneath the arc may be slightly influenced by the ²⁰⁶Pb/²⁰⁴Pb enriched source underlying the forearc. Though this effect is difficult to resolve, the Western Cascades samples appear to be slightly more variable in ²⁰⁶Pb/²⁰⁴Pb (and generally higher in ²⁰⁶Pb/²⁰⁴Pb at a given ²⁰⁷Pb/²⁰⁴Pb) than their backarc counterparts, suggesting such an influence. Both trace element and isotopic data suggest a significant subduction contribution to EOPC backarc primitive magmas. The EOPC backarc magmas appear to have originated from a heterogeneous mantle and are variably influenced by a subduction component. Though this subduction influence could be an inherited feature, the affinity between John Day and Eastern Clarno magmas and a modern Cascadia sediment source suggest that this is at least a relatively young feature. Taken together, these results demonstrate the heterogeneity of source and process across the Eocene-Oligocene arc. Apparent subduction contribution increases from the geographic forearc to backarc, perhaps indicating a wider arc than is typically envisioned at the time. The distinct ²⁰⁶Pb/²⁰⁴Pb enrichment associated with the Siletzia mantle and apparent in the forearc appears to wane in influence to the east, having at most a minor influence on the Western Cascade arc magmas. As with the Siletzia mantle domain in EOPC arc and backarc, mantle domains associated with the Western Cascades and John Day/Eastern Clarno magmas do not persist in the High Cascades or High Lava Plains. Beneath and behind the arc, mantle reservoirs appear to have been largely replaced since the Eocene-Oligocene. However, the High Lava Plains basalts appear to carry a paleo-enrichment signature that may be an inherited feature from advected mantle. If this enriched mantle was advected from beneath North America, it may explain the observation that High Lava Plains basalts are more enriched to the east. Close

• 3. [Article] Kinematic Constraints on Tremor and Slow Slip in Cascadia and Implications for Fault Properties

  Subduction zone fault processes range from tsunami-generating megathrust events to aseismic creep along the deeper portions of the fault. Episodic tremor and slow slip (ETS) represents the transition between ...

  Citation × Citation Citation Abstract Copy Title: Kinematic Constraints on Tremor and Slow Slip in Cascadia and Implications for Fault Properties Author: Krogstad, Randy Year: 2016 Subduction zone fault processes range from tsunami-generating megathrust events to aseismic creep along the deeper portions of the fault. Episodic tremor and slow slip (ETS) represents the transition between these two regimes, where slip occurs at semi-regular recurrence intervals of months-to-years. These events are also accompanied by low frequency earthquakes, referred to as tremor. The study of ETS in Cascadia has been made possible by the enhancement of large-scale seismic and geodetic networks. In this dissertation, I use a range of geodetic and seismic observations at sub-daily to decadal time scales to investigate the kinematic behavior of individual ETS events, as well as the long-term behavior of the ETS zone and its relationship with the updip seismogenic zone. In Cascadia, current seismic hazard maps use the ETS zone as the downdip limit of rupture during future megathrust events. In Chapter II, I utilize uplift rates derived from 80 years of leveling measurements to explore the possibility that long-term strain accumulation exists near the ETS zone. The uplift rates are consistent with a region of 10-20% locking on the updip side of the ETS zone. The lack of associated topography indicates that the accumulated strain must be released during the megathrust cycle. The correlation of tremor and slip in Cascadia suggests there is an inherent relationship between the two. In Chapter III, I develop a method for using tremor as a proxy for slip to assess the spatial relationship of tremor and slip. I compare predictions of tremor-derived slip models to results from static inversions of GPS offsets by modeling slip based on the density of tremor. These comparisons suggest that the correlation of tremor and slip is variable along strike and along dip. In Chapter IV, I explore how borehole strainmeters can improve our resolution of slip on the plate interface. I incorporate strainmeters into joint, time-dependent kinematic inversions with GPS data. The temporal resolution of strainmeters provides improved constraints when deriving time-dependent slip estimates during slow slip events, allowing us to better image the kinematics of slow slip. This dissertation includes previously published and unpublished material. Title: Kinematic Constraints on Tremor and Slow Slip in Cascadia and Implications for Fault Properties Author: Krogstad, Randy Year: 2016 Subduction zone fault processes range from tsunami-generating megathrust events to aseismic creep along the deeper portions of the fault. Episodic tremor and slow slip (ETS) represents the transition between these two regimes, where slip occurs at semi-regular recurrence intervals of months-to-years. These events are also accompanied by low frequency earthquakes, referred to as tremor. The study of ETS in Cascadia has been made possible by the enhancement of large-scale seismic and geodetic networks. In this dissertation, I use a range of geodetic and seismic observations at sub-daily to decadal time scales to investigate the kinematic behavior of individual ETS events, as well as the long-term behavior of the ETS zone and its relationship with the updip seismogenic zone. In Cascadia, current seismic hazard maps use the ETS zone as the downdip limit of rupture during future megathrust events. In Chapter II, I utilize uplift rates derived from 80 years of leveling measurements to explore the possibility that long-term strain accumulation exists near the ETS zone. The uplift rates are consistent with a region of 10-20% locking on the updip side of the ETS zone. The lack of associated topography indicates that the accumulated strain must be released during the megathrust cycle. The correlation of tremor and slip in Cascadia suggests there is an inherent relationship between the two. In Chapter III, I develop a method for using tremor as a proxy for slip to assess the spatial relationship of tremor and slip. I compare predictions of tremor-derived slip models to results from static inversions of GPS offsets by modeling slip based on the density of tremor. These comparisons suggest that the correlation of tremor and slip is variable along strike and along dip. In Chapter IV, I explore how borehole strainmeters can improve our resolution of slip on the plate interface. I incorporate strainmeters into joint, time-dependent kinematic inversions with GPS data. The temporal resolution of strainmeters provides improved constraints when deriving time-dependent slip estimates during slow slip events, allowing us to better image the kinematics of slow slip. This dissertation includes previously published and unpublished material. Close

• 4. [Article] Thermal effects of hydrothermal circulation and seamount subduction: Temperatures in the Nankai Trough Seismogenic Zone Experiment transect, Japan

  This is the publisher's version of record. The original submission is copyrighted by American Geophysical Union and can be found here: http://www.agu.org/

  Citation × Citation Citation Abstract Copy Title: Thermal effects of hydrothermal circulation and seamount subduction: Temperatures in the Nankai Trough Seismogenic Zone Experiment transect, Japan Author: Harris, R. N., Spinelli, G. A. This is the publisher's version of record. The original submission is copyrighted by American Geophysical Union and can be found here: http://www.agu.org/ Title: Thermal effects of hydrothermal circulation and seamount subduction: Temperatures in the Nankai Trough Seismogenic Zone Experiment transect, Japan Author: Harris, R. N., Spinelli, G. A. This is the publisher's version of record. The original submission is copyrighted by American Geophysical Union and can be found here: http://www.agu.org/ Close

• 5. [Article] Thermal effects of hydrothermal circulation and seamount subduction: Temperatures in the Nankai Trough Seismogenic Zone Experiment transect, Japan

  This is the publisher’s final pdf. The published article is copyrighted by the American Geophysical Union and can be found at: http://www.agu.org/journals/gc/.

  Citation × Citation Citation Abstract Copy Title: Thermal effects of hydrothermal circulation and seamount subduction: Temperatures in the Nankai Trough Seismogenic Zone Experiment transect, Japan Author: Spinelli, G. A., Harris, R. N. This is the publisher’s final pdf. The published article is copyrighted by the American Geophysical Union and can be found at: http://www.agu.org/journals/gc/. Title: Thermal effects of hydrothermal circulation and seamount subduction: Temperatures in the Nankai Trough Seismogenic Zone Experiment transect, Japan Author: Spinelli, G. A., Harris, R. N. This is the publisher’s final pdf. The published article is copyrighted by the American Geophysical Union and can be found at: http://www.agu.org/journals/gc/. Close

• 6. [Article] From Cinder Cones to Subduction Zones: Volatile Recycling and Magma Formation beneath the Southern Cascade Arc

  Volatiles (H2O, CO2, S, Cl) play a key role in magmatic processes at subduction zones. In this study, the dissolved volatile contents of olivine-hosted melt inclusions from cinder cones in the Lassen segment ...

  Citation × Citation Citation Abstract Copy Title: From Cinder Cones to Subduction Zones: Volatile Recycling and Magma Formation beneath the Southern Cascade Arc Author: Walowski, Kristina Year: 2015 Volatiles (H2O, CO2, S, Cl) play a key role in magmatic processes at subduction zones. In this study, the dissolved volatile contents of olivine-hosted melt inclusions from cinder cones in the Lassen segment of the Cascade arc are used to investigate dehydration of subducted oceanic lithosphere, magma formation in the sub-arc mantle wedge, and mafic magma storage and evolution in the crust. Relatively young, hot oceanic lithosphere subducts beneath the Cascade arc. The hydrogen-isotope and trace-element compositions of melt inclusions, when integrated with thermo-petrologic modeling, demonstrate that fluids in Cascade magmas are sourced from hydrated peridotite in the deep slab interior and that the oceanic crustal part of the slab extensively dehydrates beneath the forearc. In contrast to their slab-derived H, the melt inclusions have B concentrations and isotope ratios that are similar to mid-ocean ridge basalt (MORB), requiring little to no slab contribution of B, which is also consistent with extensive dehydration of the downgoing plate before it reaches sub-arc depths. Correlations of volatile and trace element ratios in the melt inclusions (H2O/Ce, Cl/Nb, Sr/Nd) demonstrate that geochemical variability in the magmas is the result of variable amounts of addition of a hydrous subduction component to the mantle wedge. Radiogenic isotope ratios require that the subduction component has less radiogenic Sr and Pb and more radiogenic Nd than the Lassen sub-arc mantle and is therefore likely to be a partial melt of subducted Gorda MORB. These results provide evidence that chlorite-derived fluids from the deep slab interior flux-melt the oceanic crust, producing hydrous slab melts that migrate into the overlying mantle, where they react with peridotite to induce further melting. The basaltic magmas that erupted at Cinder Cone near Mt. Lassen trapped melt inclusions during olivine crystallization at ~7-15 km depth. The melt inclusion compositions require that two different mantle-derived magmas were involved in the eruption, and temporal changes show that arrival of the two batches correlates with two explosive phases of activity. Both magmas experienced rapid crustal contamination before erupting, illustrating the complexities of cinder cone eruptions. This dissertation includes previously published and unpublished co-authored material. Title: From Cinder Cones to Subduction Zones: Volatile Recycling and Magma Formation beneath the Southern Cascade Arc Author: Walowski, Kristina Year: 2015 Volatiles (H2O, CO2, S, Cl) play a key role in magmatic processes at subduction zones. In this study, the dissolved volatile contents of olivine-hosted melt inclusions from cinder cones in the Lassen segment of the Cascade arc are used to investigate dehydration of subducted oceanic lithosphere, magma formation in the sub-arc mantle wedge, and mafic magma storage and evolution in the crust. Relatively young, hot oceanic lithosphere subducts beneath the Cascade arc. The hydrogen-isotope and trace-element compositions of melt inclusions, when integrated with thermo-petrologic modeling, demonstrate that fluids in Cascade magmas are sourced from hydrated peridotite in the deep slab interior and that the oceanic crustal part of the slab extensively dehydrates beneath the forearc. In contrast to their slab-derived H, the melt inclusions have B concentrations and isotope ratios that are similar to mid-ocean ridge basalt (MORB), requiring little to no slab contribution of B, which is also consistent with extensive dehydration of the downgoing plate before it reaches sub-arc depths. Correlations of volatile and trace element ratios in the melt inclusions (H2O/Ce, Cl/Nb, Sr/Nd) demonstrate that geochemical variability in the magmas is the result of variable amounts of addition of a hydrous subduction component to the mantle wedge. Radiogenic isotope ratios require that the subduction component has less radiogenic Sr and Pb and more radiogenic Nd than the Lassen sub-arc mantle and is therefore likely to be a partial melt of subducted Gorda MORB. These results provide evidence that chlorite-derived fluids from the deep slab interior flux-melt the oceanic crust, producing hydrous slab melts that migrate into the overlying mantle, where they react with peridotite to induce further melting. The basaltic magmas that erupted at Cinder Cone near Mt. Lassen trapped melt inclusions during olivine crystallization at ~7-15 km depth. The melt inclusion compositions require that two different mantle-derived magmas were involved in the eruption, and temporal changes show that arrival of the two batches correlates with two explosive phases of activity. Both magmas experienced rapid crustal contamination before erupting, illustrating the complexities of cinder cone eruptions. This dissertation includes previously published and unpublished co-authored material. Close

• 7. [Article] Volatiles in basaltic magmas from central Mexico: From subduction to eruption

  xvi, 167 p. ; ill. (some col.) A print copy of this thesis is available through the UO Libraries. Search the library catalog for the location and call number.

  Citation × Citation Citation Abstract Copy Title: Volatiles in basaltic magmas from central Mexico: From subduction to eruption Author: Johnson, Emily Renee Year: 2009, 2008 xvi, 167 p. ; ill. (some col.) A print copy of this thesis is available through the UO Libraries. Search the library catalog for the location and call number. Title: Volatiles in basaltic magmas from central Mexico: From subduction to eruption Author: Johnson, Emily Renee Year: 2009, 2008 xvi, 167 p. ; ill. (some col.) A print copy of this thesis is available through the UO Libraries. Search the library catalog for the location and call number. Close

• 8. [Article] Crustal Heat Flow and Thermal Modeling of the Hikurangi Trough, New Zealand

  Measuring surface heat flow at the Hikurangi Trough is key toward characterizing the local thermal regime and the influence of thermally sensitive processes such as fluid flow and slow slip. Marine heat ...

  Citation × Citation Citation Abstract Copy Title: Crustal Heat Flow and Thermal Modeling of the Hikurangi Trough, New Zealand Author: Antriasian, Anson M. Measuring surface heat flow at the Hikurangi Trough is key toward characterizing the local thermal regime and the influence of thermally sensitive processes such as fluid flow and slow slip. Marine heat flow data were collected during May– June 2015 in the northern Hikurangi Trough at sites seaward and landward of the deformation front, giving mean and one standard deviation values of 58 ± 16 and 46 ± 15 mW m⁻² respectively. A 3.5 m violin-bow style probe was used to measure thermal conductivity and temperature gradient of the seafloor, and heat flow is calculated using Fourier’s law. Two-dimensional steady-state, finite element modeling was undertaken to analyze the influence of subduction parameters on surface heat flow and temperature along the subduction thrust. Heat flow values seaward of the deformation front have been found to be sensitive to the structure of buried and outcropping extrusive crust, but the mean heat flow is within the range of expected values for lithosphere this age (120 Myr) according to the GDH1 cooling model. Landward of the deformation front, measured heat flow is in excess of preferred modeled values. I interpret these observations as evidence for fluid flow within a crustal aquifer system and through splay faults within the forearc. A heat flow discontinuity of 40 mW m⁻² landward of the deformation front supports the presence of compaction sourced fluid along the subduction thrust, which is consistent with previous evidence for excess pore pressure and the occurrence of shallow slow slip. Title: Crustal Heat Flow and Thermal Modeling of the Hikurangi Trough, New Zealand Author: Antriasian, Anson M. Measuring surface heat flow at the Hikurangi Trough is key toward characterizing the local thermal regime and the influence of thermally sensitive processes such as fluid flow and slow slip. Marine heat flow data were collected during May– June 2015 in the northern Hikurangi Trough at sites seaward and landward of the deformation front, giving mean and one standard deviation values of 58 ± 16 and 46 ± 15 mW m⁻² respectively. A 3.5 m violin-bow style probe was used to measure thermal conductivity and temperature gradient of the seafloor, and heat flow is calculated using Fourier’s law. Two-dimensional steady-state, finite element modeling was undertaken to analyze the influence of subduction parameters on surface heat flow and temperature along the subduction thrust. Heat flow values seaward of the deformation front have been found to be sensitive to the structure of buried and outcropping extrusive crust, but the mean heat flow is within the range of expected values for lithosphere this age (120 Myr) according to the GDH1 cooling model. Landward of the deformation front, measured heat flow is in excess of preferred modeled values. I interpret these observations as evidence for fluid flow within a crustal aquifer system and through splay faults within the forearc. A heat flow discontinuity of 40 mW m⁻² landward of the deformation front supports the presence of compaction sourced fluid along the subduction thrust, which is consistent with previous evidence for excess pore pressure and the occurrence of shallow slow slip. Close

• 9. [Article] Heat flow along the NanTroSEIZE transect: Results from IODP Expeditions 315 and 316 offshore the Kii Peninsula, Japan

  [1] We report detailed thermal measurements undertaken during Integrated Ocean Drilling Program (IODP) Expeditions 315 and 316 along the Nankai Trough Seismogenic Zone Experiment (NanTroSEIZE) transect offshore ...

  Citation × Citation Citation Abstract Copy Title: Heat flow along the NanTroSEIZE transect: Results from IODP Expeditions 315 and 316 offshore the Kii Peninsula, Japan Author: Schmidt‐Schierhorn, Friederike, Spinelli, Glenn, Harris, Robert N. [1] We report detailed thermal measurements undertaken during Integrated Ocean Drilling Program (IODP) Expeditions 315 and 316 along the Nankai Trough Seismogenic Zone Experiment (NanTroSEIZE) transect offshore the Kii Peninsula, Japan. Geothermal objectives included determining thermophysical rock properties of the cored material and characterizing the background thermal regime along this transect. New thermal conductivity measurements made with a divided bar are reported and supplement shipboard full‐ and half‐space thermal conductivity measurements for a total of 938 thermal conductivity measurements. Thermal conductivity varies between about 1.0 W m⁻¹ K⁻¹ near the seafloor to 1.6 W m⁻¹ K⁻¹ at a depth of 1 km. Thermal conductivity generally increases with depth and correlates with variations in porosity and lithology. Thermal gradients along the transect are characterized by 48 sediment temperature measurements from 6 sites. Thermal corrections for the effects of bathymetric relief and sedimentation improve the confidence with which the heat flow values can be interpreted. Heat flow generally decreases with landward distance from the deformation front and varies from 70 mW m⁻² just landward of the deformation front to 54 mW m⁻² at sites characterizing the outer fore‐arc high and to 57 mW m⁻² at the Kumano Basin Site. IODP heat flow measurements are significantly lower than nearby seafloor heat flow measurements. This difference is most likely due to variations in bottom water temperature that have a large effect on values of seafloor heat flow. Thus the heat flow of the Nankai accretionary prism is lower than previously thought. We present thermal models of subduction along this transect and explore the impact of the initial geotherm. Conductive plate cooling based on the age of subducting seafloor (20 Myr) under predicts the observed heat flow. We find a good fit to the data using a geotherm appropriate for 10 Myr seafloor. The extra heat is interpreted in terms of back‐arc thermal environments. Title: Heat flow along the NanTroSEIZE transect: Results from IODP Expeditions 315 and 316 offshore the Kii Peninsula, Japan Author: Schmidt‐Schierhorn, Friederike, Spinelli, Glenn, Harris, Robert N. [1] We report detailed thermal measurements undertaken during Integrated Ocean Drilling Program (IODP) Expeditions 315 and 316 along the Nankai Trough Seismogenic Zone Experiment (NanTroSEIZE) transect offshore the Kii Peninsula, Japan. Geothermal objectives included determining thermophysical rock properties of the cored material and characterizing the background thermal regime along this transect. New thermal conductivity measurements made with a divided bar are reported and supplement shipboard full‐ and half‐space thermal conductivity measurements for a total of 938 thermal conductivity measurements. Thermal conductivity varies between about 1.0 W m⁻¹ K⁻¹ near the seafloor to 1.6 W m⁻¹ K⁻¹ at a depth of 1 km. Thermal conductivity generally increases with depth and correlates with variations in porosity and lithology. Thermal gradients along the transect are characterized by 48 sediment temperature measurements from 6 sites. Thermal corrections for the effects of bathymetric relief and sedimentation improve the confidence with which the heat flow values can be interpreted. Heat flow generally decreases with landward distance from the deformation front and varies from 70 mW m⁻² just landward of the deformation front to 54 mW m⁻² at sites characterizing the outer fore‐arc high and to 57 mW m⁻² at the Kumano Basin Site. IODP heat flow measurements are significantly lower than nearby seafloor heat flow measurements. This difference is most likely due to variations in bottom water temperature that have a large effect on values of seafloor heat flow. Thus the heat flow of the Nankai accretionary prism is lower than previously thought. We present thermal models of subduction along this transect and explore the impact of the initial geotherm. Conductive plate cooling based on the age of subducting seafloor (20 Myr) under predicts the observed heat flow. We find a good fit to the data using a geotherm appropriate for 10 Myr seafloor. The extra heat is interpreted in terms of back‐arc thermal environments. Close

• 10. [Article] The geochemical evolution of the Aucanquilcha Volcanic Cluster : prolonged magmatism and its crustal consequences

  The interaction of magma with continental crust at convergent margins is fundamental to understanding if and how continents grow. Isotopic and elemental data constrain the progressive stages of development ...

  Citation × Citation Citation Abstract Copy Title: The geochemical evolution of the Aucanquilcha Volcanic Cluster : prolonged magmatism and its crustal consequences Author: Walker, Barry Alan The interaction of magma with continental crust at convergent margins is fundamental to understanding if and how continents grow. Isotopic and elemental data constrain the progressive stages of development of the magmatic underpinnings of the long-lived Aucanquilcha Volcanic Cluster (AVC), situated atop the thick continental crust of the central Andes in northern Chile. Whole rock data are used in conjunction with mineral compositions to infer processes that gave rise to eleven million years of intermediate, dominantly dacite, arc volcanism. A pulse of volcanic activity at the AVC between ~5 and 2 Ma is bracketed by more sluggish rates. We document chemical changes in the lavas that accompany this eruptive evolution. Trace element data suggest that crystal fractionation and magma mixing were the dominant mechanisms generating the diversity observed in the AVC whole rock data. Fractionation was dominant during early and waning stages of magmatism, and magma mixing was an important process during the high flux period. Peak thermal maturity of the AVC underpinnings coincided with the high magma flux and likely promoted open system processes during this time. Mineral compositions from zircon, amphibole, pyroxene, and Fe-Ti oxides confirm the importance of material recycling in the production of evolved AVC rocks. Various geothermometers were employed to calculate the pre-eruptive conditions of AVC magma using mineral compositions. Pressure estimates from amphibole and two-pyroxene barometry indicate crystallization depths of 1 – 5 kb and 4 – 6 kb, respectively. Temperature estimates from zircon, Fe-Ti oxides, amphiboles, and pyroxenes indicate temperatures ranging from ~700°C to 1100°C. Zircon temperatures are always the lowest (700°C - 950°C), and pyroxene temperatures are always the highest (1000°C - 1100°C), with Fe-Ti oxide and amphiboles temperatures falling in between. U-Pb ages from zircons and thermometry from individual samples evidence the thermal maturation and consolidation of the underpinnings below the AVC, presumably culminating in a large, crystal-rich mush zone where magmas were trapped and processed. It is in these middle to upper crustal zones where magmatic diversity is attenuated and giant, relatively homogeneous batholiths are formed. Isotopes of AVC lavas are similar to values observed from other central Andes volcanic centers. Lead isotopes are consistent with the AVC's location within a Pb isotope transition zone between the Antofalla and Arequipa basement terranes. Oxygen and Sr isotopic ratios are high and Nd isotopic ratios low with respect to a depleted mantle. Through time, ⁸⁷Sr/⁸⁶Sr values of AVC lavas progressively increase from lows of ~0.70507 to ~0.70579 (upper values of 0.70526 to 0.70680), and εNd values decrease from highs of -1.0 to -4.6 (lows of -1.6 to -7.3). Similarly, O isotopes (δ¹⁸O) show a slight increase in base level through time from lows of 6.5‰ to 7.0‰ (highs of 6.75‰ – 7.5‰). Dy/Yb and Sm/Yb ratios also increased systematically from highs of 2.11 to 3.45, and 2.76 to 6.67, respectively. Despite the temporal isotopic variation, there is little isotopic variation with indices of fractionation, suggesting this signal is the consequence of deep magmatic processing, here attributed to an expanding zone of melting, assimilation, storage, and homogenization (MASH) of mantle-derived magma in the deep crust. Upward expansion brought the MASH zone into contact with rocks that were increasingly evolved with respect to Sr and Nd isotopes, explaining the isotopic shifts. Downward expansion of the MASH zone enhanced garnet stability during basalt fractionation, explaining the increased Dy/Yb and Sm/Yb ratios. Mass balance calculations involving Sr, Nd, and O isotope modeling are consistent with a crustal component making up 10 - 30% of AVC lavas, implying that although the history of central Andean magmatism is replete with large scale crustal recycling, the current phase is largely a crust formation event. Title: The geochemical evolution of the Aucanquilcha Volcanic Cluster : prolonged magmatism and its crustal consequences Author: Walker, Barry Alan The interaction of magma with continental crust at convergent margins is fundamental to understanding if and how continents grow. Isotopic and elemental data constrain the progressive stages of development of the magmatic underpinnings of the long-lived Aucanquilcha Volcanic Cluster (AVC), situated atop the thick continental crust of the central Andes in northern Chile. Whole rock data are used in conjunction with mineral compositions to infer processes that gave rise to eleven million years of intermediate, dominantly dacite, arc volcanism. A pulse of volcanic activity at the AVC between ~5 and 2 Ma is bracketed by more sluggish rates. We document chemical changes in the lavas that accompany this eruptive evolution. Trace element data suggest that crystal fractionation and magma mixing were the dominant mechanisms generating the diversity observed in the AVC whole rock data. Fractionation was dominant during early and waning stages of magmatism, and magma mixing was an important process during the high flux period. Peak thermal maturity of the AVC underpinnings coincided with the high magma flux and likely promoted open system processes during this time. Mineral compositions from zircon, amphibole, pyroxene, and Fe-Ti oxides confirm the importance of material recycling in the production of evolved AVC rocks. Various geothermometers were employed to calculate the pre-eruptive conditions of AVC magma using mineral compositions. Pressure estimates from amphibole and two-pyroxene barometry indicate crystallization depths of 1 – 5 kb and 4 – 6 kb, respectively. Temperature estimates from zircon, Fe-Ti oxides, amphiboles, and pyroxenes indicate temperatures ranging from ~700°C to 1100°C. Zircon temperatures are always the lowest (700°C - 950°C), and pyroxene temperatures are always the highest (1000°C - 1100°C), with Fe-Ti oxide and amphiboles temperatures falling in between. U-Pb ages from zircons and thermometry from individual samples evidence the thermal maturation and consolidation of the underpinnings below the AVC, presumably culminating in a large, crystal-rich mush zone where magmas were trapped and processed. It is in these middle to upper crustal zones where magmatic diversity is attenuated and giant, relatively homogeneous batholiths are formed. Isotopes of AVC lavas are similar to values observed from other central Andes volcanic centers. Lead isotopes are consistent with the AVC's location within a Pb isotope transition zone between the Antofalla and Arequipa basement terranes. Oxygen and Sr isotopic ratios are high and Nd isotopic ratios low with respect to a depleted mantle. Through time, ⁸⁷Sr/⁸⁶Sr values of AVC lavas progressively increase from lows of ~0.70507 to ~0.70579 (upper values of 0.70526 to 0.70680), and εNd values decrease from highs of -1.0 to -4.6 (lows of -1.6 to -7.3). Similarly, O isotopes (δ¹⁸O) show a slight increase in base level through time from lows of 6.5‰ to 7.0‰ (highs of 6.75‰ – 7.5‰). Dy/Yb and Sm/Yb ratios also increased systematically from highs of 2.11 to 3.45, and 2.76 to 6.67, respectively. Despite the temporal isotopic variation, there is little isotopic variation with indices of fractionation, suggesting this signal is the consequence of deep magmatic processing, here attributed to an expanding zone of melting, assimilation, storage, and homogenization (MASH) of mantle-derived magma in the deep crust. Upward expansion brought the MASH zone into contact with rocks that were increasingly evolved with respect to Sr and Nd isotopes, explaining the isotopic shifts. Downward expansion of the MASH zone enhanced garnet stability during basalt fractionation, explaining the increased Dy/Yb and Sm/Yb ratios. Mass balance calculations involving Sr, Nd, and O isotope modeling are consistent with a crustal component making up 10 - 30% of AVC lavas, implying that although the history of central Andean magmatism is replete with large scale crustal recycling, the current phase is largely a crust formation event. Close