ALBERT

Description: The high elevation of the Southern Puna Plateau, the widespread melting of its crust, the gap in intermediate depth seismicity and the recent eruptions of ignimbrite complexes can be explained by delamination of the lithospheric mantle beneath it. To test this hypothesis, an array consisting of 73 broad band and short period seismic stations was deployed in the region for a period of two years starting in 2007. We inverted the data using the two plane wave approach and obtained 1D and 3D Rayleigh wave phase velocities. Our dispersion curve shows that at short periods (〈70 s) the phase velocities are slightly higher than those of the Tibetan plateau and lower than those of the Anatolian plateau. At periods of 100-140 s we observe a low velocity zone that might be remnant hot asthenosphere below a flat slab (7-10 Ma). We estimate the average continental lithosphere thickness for the region to be between 100 and 130 km. Our three dimensional Rayleigh wave phase velocities show a high velocity anomaly at low frequencies (0.007, 0.008 and 0.009 Hz) slightly to the north of Cerro Galan. This would be consistent with the hypothesis of delamination in which a piece of lithosphere has detached and caused upwelling of hot asthenosphere which in turn caused widespread alkaline-collision related volcanism. This interpretation is also corroborated by our shear wave velocity model where a high velocity anomaly beneath the northern edge of Cerro Galan at 130 km depth is interpreted as the delaminated block on top of the subducting Nazca slab.

Description: Chlorine isotope compositions of high-pressure (~2.3 GPa) serpentinite, rodingite, and hydrothermally altered oceanic crust (AOC) differ significantly from high- and ultrahigh-pressure (〉3.2 GPa) metasedimentary rocks in the Aosta region, Italy. Texturally early serpentinites, rodingites, and AOC have bulk δ 37 Cl values indistinguishable from those of modern seafloor analogues (δ 37 Cl = -1.0 to +1.0‰). In contrast, serpentinites and AOC samples that recrystallized during exhumation have low δ 37 Cl values (-2.7 to -0.5‰); 37 Cl depletion correlates with progressive changes in bulk chemistry. HP/UHP metasediments have low δ 37 Cl values (median = -2.5‰) that differ statistically from modern marine sediments (median = -0.6‰). Cl in metasedimentary rocks is concentrated in texturally early minerals, indicating modification of seafloor compositions early in the subduction history. The data constrain fluid sources during both subduction- and exhumation-related phases of fluid-rock interaction: (1) Marine sediments at the top of the downgoing plate likely interacted with isotopically light pore fluids from the accretionary wedge in the early stages of subduction. (2) No pervasive interaction with externally derived fluid occurred during subsequent subduction to the maximum depths of burial. (3) Localized mixing between serpentinites and fluids released by previously isotopically modified metasediments occurred during exhumation in the subduction channel. Most samples, however, preserved protolith signatures during subduction to near-arc depths.

Description: We have measured Ni, Ca, and Mn in olivine phenocrysts from volcanoes in the Galápagos Archipelago to infer the mantle source lithologies. Results show that peridotite is the dominant source lithology for Fernandina, Floreana, Genovesa, Wolf Island, and Darwin Island. These volcanoes largely characterize the PLUME, WD, FLO and DUM Nd, Sr, and Pb isotopic endmembers of Harpp and White (2001). Volcan Wolf, Alcedo, Marchena, and Cerro Azul, also produced from the melting of peridotite sources, have isotopic compositions that can be defined by mixing of the 4 isotopic endmembers. Our analysis suggests that peridotite was present in the sources of the volcanoes covered in this study and therefore is the dominant source lithology of the Galápagos plume. Pyroxenite melting is generally focused in two isotopically distinct domains: Roca Redonda, Volcan Ecuador, and Sierra Negra in the enriched western part of the archipelago, and Santiago, Santa Cruz, and Santa Fe in the depleted east. One implication of this finding is that the Western and Eastern Pyroxenite Domains represent two separate bodies of recycled crust within the Galápagos mantle plume. Furthermore, both isotopically enriched and depleted domains of the archipelago were generated from mixtures of peridotite and pyroxenite. This suggests that there is no relationship between the source lithology of the Galápagos plume and its isotopic characteristics. The identification of peridotite source melting in volcanoes with isotopic characteristics that have been attributed to recycled crust points to the importance of mixing in OIB genesis, consistent with studies in the Canary Islands.

Description: Paleomagnetic analysis and radiocarbon dating of an expanded Holocene deep-sea sediment sequence recovered by Integrated Ocean Drilling Program (IODP) Expedition 303 from Labrador Sea Site U1305 (Lat: 57°28.5 N, Long. 48°31.8 W, water depth 3459 m) provides insights into mechanisms that drive both paleomagnetic secular variation (PSV) and magnetization acquisition in deep-sea sediments. Seventeen radiocarbon dates on planktonic foraminifera define postglacial (c. 8 ka) sedimentation rates as ranging from 35 to 〉 90 cm/kyr. Alternating field (AF) demagnetization of u-channel samples show that these homogeneous sediments preserve a strong, stable, and consistently well-defined component magnetization. Normalized remanence records pass reliability criteria for relative paleointensity (RPI) estimates. Assuming that the age of magnetization is most accurately defined by well dated PSV records with the highest sedimentation rates, allows us to estimate and correct for temporal offsets at Site U1305 interpreted to result from post-depositional remanence acquisition at a depth of ~ 20 cm. Comparisons indicate that the northern North Atlantic PSV and RPI records are more consistent with European than North American records, and the evolution of virtual geomagnetic poles (VGP) are temporally and longitudinally similar to global reconstructions, though with much larger latitudinal variation. The largest deviations from a geocentric axial dipole (GAD) are observed during times of the highest intensities, in contrast to the usual assumption. These observations are consistent with the idea that PSV in the North Atlantic and elsewhere during the Holocene results from temporal oscillations of high latitude flux concentrations at a few recurrent locations.

Description: The conditions permitting mantle serpentinization during continental rifting are explored within 2D thermotectonostratigraphic basin models, which track the rheological evolution of the continental crust, account for sediment blanketing effects, and allow for kinetically controlled mantle serpentinization processes. The basic idea is that the entire extending continental crust has to be brittle for crustal scale faulting and mantle serpentinization to occur [ Perez-Gussinye and Reston , 2001]. The isostatic and latent heat effects of the reaction are fully coupled to the structural and thermal solutions. A systematic parameter study shows that a critical stretching factor exists for which complete crustal embrittlement and serpentinization occurs. Increased sedimentation rates shift this critical stretching factor to higher values as sediment blanketing effects result in higher crustal temperatures. Sediment supply has therefore, through the temperature-dependence of the viscous flow laws, strong control on crustal strength and mantle serpentinization reactions are only likely when sedimentation rates are low and stretching factors high. In a case study for the Norwegian margin we test whether the inner lower crustal bodies (LCB) imaged beneath the Møre and Vøring margin could be serpentinized mantle. Multiple 2D transects have been reconstructed through the 3D data set by Scheck-Wenderoth and Maystrenko [2011]. We find that serpentinization reactions are possible and likely during the Jurassic rift phase. Predicted thicknesses and locations of partially serpentinized mantle rocks fit to information on LCBs from seismic and gravity data. We conclude that some of the inner LCBs beneath the Norwegian margin may be partially serpentinized mantle.

Description: The origin of the Bermuda swell and volcanism remains enigmatic. The lack of an associated time-progressive hotspot track and absence of present-day volcanic activity make it difficult to reconcile with a deep mantle plume model. We analyze shear wave splitting measurements to estimate mantle flow direction and receiver function stacks to place constraints on the mantle transition zone thermal structure. *KS phases exhibit well-resolved null arrivals (no splitting) beneath the swell over a range of back azimuths. We find that the 410 and 660 km discontinuities are 49 ± 5 km and 19 ± 5 km deeper than the global average, respectively, leading to a transition zone thickness that is 27 ± 4 km thinner than average. Together, an apparently isotropic upper mantle and a thinned mantle transition zone suggest that mantle flow is primarily vertical beneath the swell, consistent with the presence of hot, buoyant mantle at depth.

Description: We present an improved density model and a new structural map of the Neapolitan Yellow Tuff caldera, the active portion of the nested Campi Flegrei caldera. The model was built using a new 3D inversion of the available high-precision gravity data, and a new digital terrain and marine model. The inversion procedure, based on a variable-depth lumped assembling of the subsurface gravity distribution via cell aggregation, gives better defined insights into the internal caldera architecture, that well agree with the available geological, geophysical and geochemical data. The adopted 3D gravity method is highly efficient for characterizing the shallow caldera structure (down to 3 km depth) and defining features related to regional or volcano tectonic lineaments and dynamics. In particular, the resulting density distribution highlights a pronounced density low in correspondence of the central portion of the caldera with a detail not available till now. The joint interpretation of the available data, suggests a subsurface structural setting that supports a piecemeal collapse of the caldera, and allows the identification of its headwall. Positive gravity anomalies localize dense intrusions (presently covered by late volcanic deposits) along the caldera marginal faults, and the main structural lineaments both bordering the resurgent block and cutting the caldera floor. These results allow us to both refine the current geological-structural framework and propose a new structural map that highlights the caldera boundary and its internal setting. This map is useful to interpret the phenomena occurring during unrest, and to improve both short- and long-term volcanic hazards assessment.

Description: Catastrophic collapses of submarine volcanoes have the potential to generate major tsunami, threatening many coastal populations. Recognizing the difficulties surrounding anticipations of these events, quantitative assessment of collapse-prone regions based on detailed morphological, geological and geophysical mapping can still provide important information about the hazards associated with these collapses. Rumble III is one of the shallowest, and largest, submarine volcanoes found along the Kermadec arc, and is both volcanically and hydrothermally active. Previous surveys have delineated major collapse features at Rumble III; based on time-lapse bathymetry, dramatic changes in the volcano morphology have been shown to have occurred over the interval 2007 to 2009. Furthermore, this volcano is located just ˜300 km from the east coast of the North Island of New Zealand. Here, we present a geophysical model for Rumble III, that provides the locations and sizes of potential weak regions of this volcano. Shipborne and near-seafloor geological and geophysical data collected by the AUV Sentry are used to determine the subsurface distribution of weak and unstable volcanic rocks. The resulting model provides evidence for potentially unstable areas located in the Southeastern flank of this volcano which should be included in future hazard predictions.

Description: Broadband seismic experiments over the last two decades have produced dense data coverage across Tibet. Yet, the mechanism of the India-Asia lithospheric convergence beneath it remains a puzzle, with even its basic features debated and with very different end-member models advocated today. We measured highly accurate Rayleigh- and Love-wave phase-velocity curves in broad period ranges (up to 5-200 s) for a few tens of pairs and groups of stations across Tibet, combining, in each case, hundreds to thousands of inter-station measurements made with cross-correlation and waveform-inversion methods. Robust shear-velocity profiles were then determined by extensive series of non-linear inversions of the data, designed to constrain the depth-dependent ranges of isotropic-average shear speeds and radial anisotropy. Temperature anomalies in the upper mantle were estimated from shear velocities using accurate petro-physical relationships. Our results reveal strong heterogeneity in the upper mantle beneath Tibet. Very large high-velocity anomalies in the upper mantle are consistent with the presence of underthrust (beneath southwestern Tibet) and subducted (beneath central and eastern Tibet) Indian lithosphere. The lithosphere. In contrast to the Indian lithosphere, Tibetan lithosphere and asthenosphere display low to normal shear speeds; Tibetan lithosphere is thus warm and thin. Radial anisotropy in the upper mantle is weak in central and strong in northeastern Tibet, possibly reflecting asthenospheric flow above the subducting Indian lithospheric slab.

Description: The Pacific Northwest (PNW) has experienced voluminous intraplate volcanism over the past ~17 Ma, beginning with the Steens/Columbia River flood basalts and continuing with the still-ongoing volcanism in the High Lava Plains (HLP) and eastern Snake River Plain (SRP). Here we present two complementary datasets (SKS splitting and Rayleigh wave phase velocity anisotropy) that place constraints on the anisotropic structure of the upper mantle beneath the HLP and SRP regions. Beneath the HLP, SKS phases reveal dominantly E-W fast splitting directions and large (up to ~2.7 sec) delay times, with pronounced lateral variations in δ t . Lateral and depth variability in the strength of anisotropy beneath the HLP is also evident from Rayleigh wave dispersion. Beneath the SRP, SKS splitting delay times are much smaller (~0.5 sec) and surface wave observations suggest a region of upper mantle anisotropy (~50-150 km depth) with a geometry that deviates significantly from the generally plate motion parallel fast directions observed just outside of the SRP. Beneath the HLP, the geometry of the anomalously strong anisotropy is similar to the anisotropy in the deeper parts of the upper mantle, resulting in constructive interference and large SKS splitting delay times. Beneath the SRP, the geometry of the anomalous anisotropic region in the shallow mantle is different, resulting in destructive interference and reduced SKS splitting delay times. We discuss several possible explanations for these observations, including variations in olivine lattice preferred orientation (LPO) strength, transitions in olivine fabric type, and a contribution from aligned partial melt.

Description: The bulk composition of the silicate portion of the Earth (BSE) has long been assumed to be tied to chondrites, in which refractory, lithophile elements like Sm and Nd exist in chondritic relative abundances. However, the 142 Nd/ 144 Nd ratios of modern terrestrial samples are 18±5 ppm higher than the ordinary-chondrite reservoir, and this challenges the traditional BSE model. Here we investigate a hypothesis that this terrestrial 142 Nd excess is related to a Sm/Nd ratio 6% higher than chondritic. This Sm/Nd ratio yields a superchondritic 143 Nd/ 144 Nd (~0.5130) similar to that identified in the highest 3 He/ 4 He mantle reservoir, and we argue that this reservoir represents the BSE composition for lithophile elements. We develop a compositional model for BSE in which the elevated Sm/Nd requires a shift of 143 Nd/ 144 Nd from 0.51263 (chondritic) to 0.51300. The new BSE composition is depleted in highly incompatible elements, including K, relative to the chondrite-based BSE, and offers a solution the “missing” 40 Ar paradox. This BSE compositional model requires that 〉83% of the mantle is depleted to form continental crust. It also implies a ~30% reduction in BSE U, Th and K, and therefore in the current rate of radiogenic heating and, thus, a proportional increase in the heat flow delivered to surface by plate tectonics. We explore thermal history models including effects related to a newly recognized evolution in the style of plate tectonics over Earth history: The lower radiogenic heat production may delay the onset of core convection and dynamo action to as late as 3.5 Gyr.

Description: The internal geological structure of the Northeast German Basin (NEGB) is affected by intense salt diapirism and by the presence of several stratified aquifer complexes of regional relevance. The shallow Quaternary to late Tertiary freshwater aquifer is separated from the underlying Mesozoic saline aquifers by an embedded Tertiary clay enriched aquitard (Rupelian Aquitard). An important feature of this aquitard is that hydraulic connections between the upper and lower aquifers do exist in areas where the Rupelian Aquitard is missing (hydrogeological windows). Three-dimensional thermohaline numerical simulations are carried out to investigate the effects of such hydrogeological windows in the Rupelian Aquitard on the resulting groundwater, temperature and salinity distributions. Numerical results suggest that hydrogeological windows act as preferential domains of hydraulic interconnectivity between the different aquifers at depth, and enable vigorous heat and mass transport which causes a mixing of warm and saline groundwater with cold and less saline groundwater within both aquifers. In areas where the Rupelian Aquitard confines the Mesozoic aquifer, dissolved solutes from major salt structures are transported laterally giving rise to plumes of variable salinity content ranging from few hundreds of meters to several tens of kilometers. Furthermore, destabilizing thermal buoyancy forces may overwhelm counteracting stabilizing salinity induced forces offside of salt domes. This may result in buoyant upward groundwater flow transporting heat and mass to shallower levels within the same Mesozoic Aquifer.

Description: Long continuous seismic data recorded at five broadband seismic stations during 2006 at Campi Flegrei caldera have been analyzed. Introducing a coarse-grained method, we evaluate the time evolution of amplitude and polarization of the seismic noise in the frequency band common to Long-Period events. The series are modulated on tidal time scales: the root-mean square is basically dominated by solar contribution, while the azimuth of the polarization vector shows lunar diurnal and semidiurnal constituents. In addition, we find that in the frequency band common to Long-Period events the azimuths are polarized towards a specific area, suggesting that these persistent oscillations can be induced by the activity of the shallow geothermal reservoir.

Description: We report on the mineralogical assemblages found in the hyper-alkaline springs hosted on Liguria and Oman ophiolites based on exhaustive XRD and SEM analyses. In Liguria, hyper-alkaline springs produce a thin brownish calcite precipitate that covers the bedrock due to the concomitant atmospheric CO 2 uptake and neutralization of the hyper-alkaline waters. No brucite and portlandite minerals are observed. The discharge of alkaline waters in Oman ophiolite forms white-orange precipitates. Calcium carbonate minerals (calcite and/or aragonite) are the most abundant and ubiquitous precipitates and are produced by the same mechanism as in Liguria. This process is observed as a thin surface crust made of rhombohedral calcite. Morphological features of aragonite vary from needle-, bouquet-, dumbbell-, spheroidal-like habitus according to the origin of carbon, the temperature and the ionic composition of the hyper-alkaline springs, and the biochemical and organic compounds. Brucite is observed both at hyper-alkaline springs located at the thrust plane and at the paleo-Moho. The varying mixing proportions between the surface run-off waters and the hyper-alkaline ones control brucite precipitation. The Layered Double Hydroxide minerals occur solely in vicinity of hyper-alkaline springs emerging within the bedded gabbros. Finally, the dominant mineralogical associations we found in Oman (Ca-bearing carbonates and brucite) in a serpentinizing environment driven by the meteoric waters are surprisingly the same as those observed at the Lost City hydrothermal site in a totally marine environment.

Description: A computer program (PBUQ) that uses Monte Carlo simulations to propagate uncertainty through regression equations and the equation for the paleosol carbonate CO 2 paleobarometer is presented. PBUQ includes options for all of the common approaches to determining values for input variables and incorporates several recent advancements relevant to determining values for soil-respired CO 2 concentrations, δ 13 C values of respired CO 2 , δ 13 C values of atmospheric CO 2 and temperatures of soil carbonate formation. PBUQ is intended to improve confidence in paleoatmospheric CO 2 research by helping researchers draw statistically significant conclusions. PBUQ can also be used to attribute and partition error among various sources and thereby advance this technique. Sensitivity analysis indicates that S(z) is the largest source of uncertainty for most paleosols and that uncertainty is minimized for soils in which CO 2 is an evenly balanced mixture between soil-derived and atmospheric components. Evenly balanced mixtures are most likely for paleosols formed in deserts and for weakly-developed paleosols. Development of proxies for soil-respired CO 2 concentrations and δ 13 C values of soil-respired CO 2 specifically for such soils is perhaps the most crucial next step for improving this technique. Currently, calcic paleosols are best used to test the significance of trends and/or differences among time slices in paleoatmospheric CO 2 concentration. Application to quantifying Earth System Sensitivity will require large scale averaging of determinations from individual paleosols and/or reduced uncertainty associated with input variables.

Description: Long-lived detachment faults are now known to be important in tectonic evolution of slow-spreading mid-ocean ridges, and there is increasing evidence that fluid flow plays a critical role in development of detachment systems. Here we document a new manifestation of low-temperature hydrothermal venting associated with the detachment fault that formed Kane Megamullion ~3.3-2.1 m.y. ago in the western rift-valley wall of the Mid-Atlantic Ridge. Hydrothermal effects on the detachment surface include 1) cemented mounds of igneous rock and chalk debris containing hydrothermal Mn oxides and Fe oxyhydroxides, and 2) layered deposits of similar Fe-Mn minerals ± interbedded chalks. Mounds are roughly conical, ~1-10 meters high, and contain primarily basalts with lesser gabbro, serpentinite, and polymict breccia. The layered Fe-Mn-rich sediments are flat-bedded to contorted and locally are buckled into low-relief linear or polygonal ridges. We propose that the mounds formed where hydrothermal fluids discharged through the detachment hanging wall near the active fault trace. Hydrothermal precipitates cemented hanging-wall debris and welded it to the footwall, and this debris persisted as mounds as the footwall was exhumed and surrounding unconsolidated material sloughed off the sloping detachment surface. Some of the layered Fe-Mn-rich deposits may have precipitated from fluids discharging from the hanging-wall vents, but they also precipitated from low-temperature fluids venting from the exposed footwall through overlying chalks. Observed natural disturbance and abnormally thin hydrogenous Fe-Mn crusts on some contorted, hydrothermal Fe-Mn-rich chalks on ~2.7 Ma crust suggest diffuse venting that is geologically recent. Results of this study imply that there are significant fluid pathways through all parts of detachment systems and that low-temperature venting through fractured detachment footwalls may continue for several million years off-axis.

Description: P and S relative arrival time residuals from teleseismic earthquakes recorded on over 60 temporary AfricaArray broadband seismic stations deployed in Uganda, Tanzania and Zambia between 2007 and 2011 have been inverted, together with relative arrival time residuals from earthquakes recorded by previous deployments, for a tomographic image of mantle wave speed variations extending to a depth of 1200 km beneath eastern Africa. The image shows a low wave speed anomaly (LWA) well developed at shallow depths (100-200 km) beneath the Eastern and Western branches of the Cenozoic East African rift system and northwestern Zambia, and a fast wave speed anomaly at depths ≤ 350 km beneath the central and northern parts of the East African Plateau and the eastern and central parts of Zambia. At depths ≥350 km the LWA is most prominent under the central and southern parts of the East African Plateau and dips to the southwest beneath northern Zambia, extending to a depth of at least 900 km. The amplitude of the LWA is consistent with a ~150-300 K thermal perturbation, and its depth extent indicates that the African superplume, originally identified as a lower mantle anomaly, is likely a whole mantle structure. A superplume extending from the core-mantle boundary to the surface implies an origin for the Cenozoic extension, volcanism and plateau uplift in eastern Africa rooted in the dynamics of the lower mantle.

Description: In the Garibaldi Belt, the northern segment of the Cascade arc, basalts at Bridge River Cones, Salal Glacier, and Mt. Meager (BSM volcanic centers) are alkalic, atypical for an arc setting. Subduction signatures are negligible or absent from primitive alkalic basalts from Salal Glacier and Bridge River, while altered oceanic crust may have contributed a minimal amount of fluid at Mt. Meager. More evolved BSM basalts display trace element signatures considered typical of arc lavas, but this is a consequence of deep crustal assimilation rather than primary input from the subducted slab. Primary BSM basalts represent 3-8% melts that segregated from enriched garnet lherzolite at significantly higher temperatures and pressures (70-105 km) than calc-alkaline Cascade arc basalts. The BSM mantle source is significantly more incompatible element-enriched than the depleted mantle tapped by calc-alkaline Cascade arc basalts. The BSM basalts are also isotopically distinct from calc-alkaline Cascade arc basalts, more similar to MORB and intraplate basalts of the NE Pacific and NW North America. The relatively deep, hot, and geochemically distinct mantle source for BSM basalts is consistent with upwelling asthenosphere. The BSM volcanic centers are close to the projected trace of the Nootka fault, which forms the boundary between the subducting Juan de Fuca plate and the near-stagnant Explorer plate. A gap or attenuated zone between the plates may promote upwelling of enriched asthenosphere that undergoes low-degree decompression melting to generate alkalic basalts that are essentially free of slab input yet occur in an arc setting.

Description: The temporal evolution of the mantle melting processes in the Asal Rift is evaluated from the chemical composition of 56 new lava flows sampled along 10 km of the rift axis and 9 km off-axis (i.e., erupted within the last 620 ky). Petrological and primary geochemical results show that most of the samples of the inner floor of the Asal Rift are affected by plagioclase accumulation. Trace element ratios and major element compositions corrected for mineral accumulation and crystallization show a symmetric pattern relative to the rift axis and preserved a clear signal of mantle melting depth variations. While FeO, Fe 8.0 , Zr/Y and (Dy/Yb) N decrease from the rift shoulders to the rift axis, SiO 2 , Na/Ti, Lu/Hf increase and Na 2 O and Na 8.0 are constant across the rift. These variations are qualitatively consistent with shallow melting beneath the rift axis and deeper melting for off-axis lava flows. Na 8.0 and Fe 8.0 contents show that beneath the rift axis, melting paths are shallow, from 81 ± 4 km to 43 ± 5 km. These melting paths are consistent with adiabatic melting in normal-temperature fertile asthenosphere, beneath an extensively thinned mantle lithosphere. On the contrary, melting on the rift shoulders (from 107 ± 7 km to 67 ± 8 km) occurred beneath thicker lithosphere, requiring a mantle solidus temperature 100 ± 40 °C hotter. In this geodynamic environment, the calculated rate of lithospheric thinning appears to be 4.0 ± 2.0 cm yr -1 , a value close to the mean spreading rate (2.9 ± 0.2 cm yr -1 ) over the last 620 ky.

Description: Global Positioning System (GPS) and Differential Interferometric Synthetic Aperture Radar (DInSAR) data, collected from July 2007 to July 2008 on Mt. Etna, are analysed to define the dynamics preceding and accompanying the onset of the eruption on 13 May 2008. Short and long-term comparisons have been made on both GPS and radar data, covering similar time windows. Thanks to the availability of three GPS surveys the year before the eruption onset, an increase in the seawards movement of the NE flank of the volcano has been detected in the few months before the dike intrusion. The GPS ground deformation pattern also shows a slight inflation centred on the western side of the volcano in the pre-eruptive long-term comparison (from July 2007 to May 2008). The GPS has been integrated with DInSAR data by the SISTEM approach, to take advantage of the different methodologies and provide high spatial sampling of the 3D ground displacement pattern. We inverted the SISTEM results in order to model the pressure source causing the observed pre-eruptive inflation. The subsequent emplacement of the eruptive dike was imaged by two GPS surveys carried out on a dense network over the uppermost part of the volcano on May 6 and 13, i.e. a few days before and a few hours after the beginning of the eruption. We inverted this comparison to define the position, geometry and kinematics of the dike. The dike intrusion was also imaged by DInSAR data with temporal baselines of 2-3 months, which confirm strong displacements localized on the summit area, rapidly decreasing towards the middle flanks of the volcano, as detected by very short-term GPS data; furthermore, the comparison between DInSAR and GPS data highlighted the presence of a depressurizing source localized beneath the upper south-western area, acting just after the dike intrusion. Finally, the long period (one year) GPS and DInSAR data were integrated by SISTEM in order to finely depict the 3D ground deformation pattern with the highest spatial resolution. The long-period data allowed the complex kinematics of the volcano to be finely imaged and highlighting the interaction between flank dynamics and magma injection.

Description: A regional survey of alkaline springs in Oman and Ligurian ophiolites shows that the alkaline water compositions significantly vary from one ophiolite to the other and within the same ophiolite. The first order correlation between the Na (and K) and Cl concentrations points to fluid compositions only partly due to evaporation. The scatter around the evaporation line implies that Na and Cl may not be conservative during the alteration of the ultramafic rocks. Mg is almost entirely depleted at pH 〉 10.5 as a result of serpentine formation within the ultramafic body and of brucite (and minor hydrotalcite) precipitation at the springs. Ca accumulates in the high pH fluids and is consumed by Ca-carbonate formation at the springs, by mixing with river waters or by the CO2 supply from the atmosphere. Thermodynamic calculations show that brucite saturation is reached at pH values around 10.5 which triggers major changes in the water composition. The waters evolve from a quartz-saturated low pH continental environment to a brucite-dominated high pH serpentinizing system at low temperature. The highest water salinities are found in springs located along the basal thrust plane of the ophiolite. The highest Al concentrations are found in some springs located on the crustal side of the mantle/crust boundary. This poses the question of the hydrologic pathways and of the role of the mineralogical composition of the altered formations.

Description: [1]  This study examines the effect of gas hydrate formation on seismic wave velocities of fine-grained sediments. Synthesis of gas hydrates in fine-grained sediments has proved to be challenging, and how hydrate formation would affect the seismic wave velocities and stiffness of clay-rich sediments has not yet been fully understood. In this study, CO 2 hydrate was synthesized in remolded and partially water-saturated clayey silt sediments that were originally cored from a hydrate occurrence region in the Ulleung Basin, East Sea, offshore Korea. After achieving excess water conditions, compressional wave and shear wave velocities were measured for different hydrate saturations and under different vertical effective stresses. The results reveal that the compressional wave velocity V P and shear wave velocity V S increase and the stress-dependency of V P and V S decreases as the hydrate saturation S H increases from 0% to ~60%. In particular, the V S ‒S H trend lies between the grain-cementing model and the load-bearing model, suggesting that gas hydrate formation in clayey silt sediments causes weak cementation from a hydrate saturation less than ~28%. The weak cementation in fine-grained sediments can be explained by the breakage of hydrate bonds that are cementing grains during sediment compression and/or the innate weakness in bonding between hydrate crystals and fine mineral grains owing to the presence of unfrozen water films on clay mineral surfaces. In addition, it is found that at low S H the cementation effect on V P is masked by the high stiffness of pore-filling phases, but it becomes pronounced at S H greater than 47%.

Description: [1]  Mountain rivers play a key role in the delivery of particulate organic carbon (POC) to large river systems and the ocean. Due to the extent of its drainage area and runoff, the Amazon River is one of Earth's most important biogeochemical systems. However the source of POC eroded from the humid region of the Eastern Andes and the input of fossil POC from sedimentary rocks (POC fossil ) remains poorly constrained. Here we collected suspended sediments from the Kosñipata River during flood events to better characterise Andean POC, measuring the nitrogen to organic carbon ratio (N/C), stable carbon isotopes ( δ 13 C org ) and radiocarbon ( Δ 14 C org ). Δ 14 C org values ranged from -711‰ to -15‰ and significant linear trends between Δ 14 C org, N/C and δ 13 C org suggested that this reflects the mixing of POC fossil with very young organic matter ( Δ 14 C org  ~ 50‰) from the terrestrial biosphere (POC non-fossil ). Using N/C and Δ 14 C org in an end member mixing analysis, we quantify the fraction of POC fossil (to within 0.1) and find that it contributes a constant proportion of the suspended sediment mass (0.37 ± 0.03%) and up to 80% of total POC. In contrast, the relative contribution of POC non-fossil was variable, being most important during the rising limb and peak discharges of flood-events. The new data shed light on published measurements of ‘old’ POC (low Δ 14 C org ) in Andean-fed tributaries of the Amazon River, with their Δ 14 C org and δ 13 C org values consistent with variable addition of POC fossil . The findings suggest a greater persistence of Andean POC in the lowland Amazon than previously recognised.

Description: [1]  The D’ Entrecasteaux Island (DEI) gneiss domes are fault-bounded domes with ~2.5 km of relief exposing ultrahigh-pressure (UHP) and high-pressure (HP) metamorphic gneisses and migmatites exhumed in an Oligocene-Miocene arc-continent collision and subduction zone subject to Late Miocene to Recent continental extension. Multi-channel seismic (MCS) reflection data and well data show the Trobriand basin formed as a forearc basin caused by southward Miocene subduction at the Trobriand trench. Subduction slowed at ~8 Ma as the margin transitioned to an extensional tectonic environment. Since then, the Trobriand basin has subsided 1–2.5 km as a broad sag basin with few normal faults deforming the basin fill. South of the DEI, the Good enough rift basin developed after extension began (~8 Ma) as the hanging-wall of the north-dipping Owen-Stanley normal fault that bounds the basin's southern margin. The lack of upper crustal extension accompanying subsidence in the Trobriand and Good enough basins suggests depth-dependent lithospheric extension since 8 Ma has accompanied uplift of the DEI gneiss domes. Structural reconstructions of seismic profiles show 2.3 to 13.4 km of basin extension in the upper crust, while syn-rift basin subsidence values indicate at least 20.7 to 23.6 km of extension occurred in the entire crust since ~8 Ma. Results indicating thinning is preferentially accommodated in the lower crust surrounding the DEI are used to constrain a schematic model of uplift of the DEI domes involving vertical exhumation of buoyant, post-orogenic lower crust, far-field extension from slab rollback, and an inverted two-layer crustal density structure.

Description: [1]  The relationship between lithospheric evolution of eastern Eurasia and subduction of the Pacific plate has long been debated. However, the timing and implications of subduction on the tectonics of eastern China are not well constrained. Here, we present new zircon U-Pb ages and Hf isotopes, elemental and Sr-Nd-Pb isotopic data on Cretaceous volcanic rocks from the Ningwu basin, eastern China to further address this issue. Our age data reveal rapid eruption of the volcanic rocks within a short duration from 133 to 130 Ma. The rocks, mostly characterized by shoshonitic and high-K calc-alkaline signatures, display light rare earth element and Pb enrichment, Nb, Ta and Ti depletion, highly radiogenic Sr-Pb isotopic ratios and variable ε Hf (t) (+1.8 to −10), suggesting derivation from an enriched lithospheric mantle metasomatized by marine sediments. The early lavas (133.3 ± 1.1 Ma) show stronger subducted-related signatures than the late lavas (130.1 ± 1.0 Ma), which we interpret to reflect consumption of a significant volume of fusible subducted components in the early melting phase. The large ε Hf (t) variation of late lavas suggests greater involvement of asthenospheric melts and lower crust in their petrogenesis. The youngest age (130 Ma) appears to coincide with an inferred change in the direction of Pacific-Eurasia convergence, manifested as a change from extension to transpression in eastern China. The narrow window of eruption may signify a rapid change of the tectonic regime in the Early Cretaceous.

Description: Mid-ocean ridges magmatism is by and large considered to be mostly dry. Nevertheless, numerous works in the last decade have shown that a hydrous component is likely to be involved in ocean ridges magmas genesis and / or evolution. The petrology and geochemistry of peculiar coarse grained gabbros sampled in the upper part of the gabbroic sequence from the Northern Oman ophiolite (Wadi Rajmi) provide information on the origin and fate of hydrous melts in fast spreading oceanic settings. Uncommon crystallization sequences for oceanic settings (clinopyroxene crystallizing before plagioclase), extreme mineral compositions (plagioclase An% up to 99, and clinopyroxene Mg# up to 96), and the presence of magmatic amphibole, imply the presence of a high water activity during crystallization. Various petrological and geochemical constraints point to hydration resulting from the recycling of hydrothermal fluids. This recycling event may have occurred at the top of the axial magma chamber where assimilation of anatectic hydrous melts is recurrent along mid-ocean ridges, or close to segments ends where fresh magma intrudes previously hydrothermally altered crust. In ophiolitic settings, hydration and remelting of hydrothermally altered rocks producing hydrous melts may also occur during the obduction process. Although dry magmatism dominates oceanic magmatism, the dynamic behavior of fast spreading ocean ridge magma chambers has the potential to produce the observed hydrous melts (either in ophiolites or at spreading centers), which are thus part of the general mid-ocean ridges lineage.

Description: The Seismic Array HiKurangi Experiment (SAHKE) investigated the structure of the forearc and subduction plate boundary beneath the southern North Island along a 350 km transect. Tomographic inversion of first-arrival travel times was used to derive a 15-20 km deep P-wave image of the crust. The refracted phases and migrated reflection events image subducting slab geometry and crustal structure. In the west, Australian Plate Moho depth decreases westward across the Taranaki Fault system from 35 to ~28-30 km. In the east, subducted Pacific Plate oceanic crust is recognised to have a positive velocity gradient, but becomes less distinct beneath the Tararua Ranges, where the interface increases in dip at about 15 km depth from 〈5° to 〉15°. This bend in the subducted plate is associated with vertical clusters in seismicity, splay fault branching, and low-velocity high-attenuation material that we interpret to be an underplated subduction sedimentary channel. We infer that a step down in the decollément transfers slip on the plate interface at the top of a subduction channel to the oceanic crust and drives local uplift of the Tararua Ranges. Reflections from the Wairarapa Fault show that it is listric and soles into the top of underplated sediments, which in turn abut the Moho of the over-riding plate at ~32 km depth, near the downdip end of the strongly locked zone. The change in dip of the Hikurangi subduction interface is spatially correlated with the transition from geodetically determined locked to unlocked areas of the plate interface.

Description: [1]  The origin of crustal-scale silicic magmatism remains a matter of debate, and notable uncertainty exists concerning the physical mechanisms that drive ascent and emplacement of felsic magmas in upper crustal regions. A 2D numerical model demonstrates that injection of mantle-derived mafic magma into a partially molten hot zone in the lower crust can drive felsic magma ascent and intrusion into upper crustal levels. The injection of mafic magma induces over pressure in the reservoir, which increases crustal stresses and triggers development of brittle/plastic shear zones, and can drive significant surface uplift. The emerging topography causes a non-uniform over pressure distribution in the reservoir and can trigger felsic magma ascent along crustal shear zones. Based on systematic numerical experiments we investigate the influence of crustal strength and injection rate. The initial upper crustal strength controls the degree of crustal faulting and surface uplift and, therefore, whether felsic magma ascent can be initiated or not. The final upper crustal strength influences the depth and final style of felsic intrusion. The injection rate of mafic magma determines the time scale of overpressure growth and surface uplift stage. In contrast, the duration of the subsequent felsic ascent and intrusion emplacement stages remains nearly constant. Our results imply that mafic underplating and intrusion into the lower crust may not only be a prime control for the generation of felsic magmas in the lower crust, but may also be an important physical driving mechanism for felsic magma ascent and intrusion into upper crustal levels.

Description: The Paleocene-Eocene Thermal Maximum (PETM) is marked by a prominent negative carbon isotope excursion (CIE) of 3-5‰ that has a characteristic rapid onset, stable body, and recovery to near pre-CIE isotopic composition. Although the CIE is the major criterion for global correlation of the Paleocene-Eocene boundary, spatial variations in the position and shape of the CIE have not been systematically evaluated. We measured carbon isotope ratios of bulk organic matter (δ 13 C org ) and pedogenic carbonate (δ 13 C carb ) at six PETM sections across a 16 km transect in the SE Bighorn Basin, Wyoming. Bed tracing and high-resolution floral and faunal biostratigraphy allowed correlation of the sections independent of chemostratigraphy. The onset of the CIE in bulk organic matter at all six sections occurs within a single laterally extensive geosol. The magnitude of the CIE varies from 2.1-3.8‰. The absolute and relative stratigraphic thickness of the body of the CIE in bulk organic matter varies significantly across the field area and under-represents the thickness of the PETM body by 30%-80%. The variations cannot be explained by basinal position and instead suggest that δ 13 C org values were influenced by local factors such as reworking of older carbon. The stratigraphic thickness and shape of the CIE have been used to correlate sections, estimate timing of biotic and climatic changes relative to the presumed carbon isotope composition of the atmosphere, and calculate rates of environmental and biotic change. Localized controls on δ 13 C org values place these inferences in question by influencing the apparent shape and duration of the CIE.

Description: Estimates of the relative motion between the Hawaiian and Louisville hotspots have consequences for understanding the role and character of deep Pacific-mantle return flow. The relative motion between these primary hotspots can be inferred by comparing the age records for their seamount trails. We report 40 Ar/ 39 Ar ages for 18 lavas from 10 seamounts along the Hawaiian-Emperor Seamount Chain (HESC), showing that volcanism started in the sharp portion of the Hawaiian-Emperor Bend (HEB) at ≥47.5 Ma and continued for ≥5 Myr. The slope of the along-track distance from the currently active Hawaiian hotspot plotted versus age is constant (57±2 km/Myr) between ~57 and 25 Ma in the central ~1900 km of the seamount chain, including the HEB. This model predicts an age for the oldest Emperor Seamounts that matches published ages, implying that a linear age-distance relationship might extend back to at least 82 Ma. In contrast, Hawaiian age progression was much faster since at least ~15 Ma and possibly as early as ~27 Ma. Linear age-distance relations for the Hawaii-Emperor and Louisville seamount chains predict ~300 km overall hotspot relative motion between 80 and 47.5 Ma, in broad agreement with numerical models of plumes in a convecting mantle, and paleomagnetic data We show that a change in hotspot relative motion may also have occurred between ~55 Ma and ~50 Ma. We interpret this change in hotspot motion as evidence that the HEB reflects a combination of hotspot and plate motion changes driven by the same plate/mantle reorganization.

Description: ABSTRACT Axial Seamount, an active submarine volcano on the Juan de Fuca Ridge at 46°N, 130°W, erupted in January 1998 along 11 km of its upper south rift zone. We use ship-based multibeam sonar, high-resolution (1-m) bathymetry, sidescan sonar imagery, and submersible dive observations to map four separate 1998 lava flows that were fed from eleven eruptive fissures. These new mapping results give an eruption volume of 31 x 10 6 m 3 , 70% of which was in the northern-most flow, 23% in the southern-most flow, and 7% in two smaller flows in between. We introduce the concept of map-scale submarine lava flow morphology (observed at a scale of 100s of meters, as revealed by the high-resolution bathymetry), and an interpretive model in which two map-scale morphologies are produced by high effusion-rate eruptions: “inflated lobate flows” are formed near eruptive vents, and where they drain downslope more than 0.5-1.0 km, they transition to “inflated pillow flows.” These two morphologies are observed on the 1998 lava flows at Axial. A third map-scale flow morphology that was not produced during this eruption, “pillow mounds,” is formed by low effusion-rate eruptions in which pillow lava piles up directly over the eruptive vents. Axial Seamount erupted again in April 2011 and there are remarkable similarities between the 1998 and 2011 eruptions, particularly the locations of eruptive vents and lava flow morphologies. Because the 2011 eruption reused most of the same eruptive fissures, 58% of the area of the 1998 lava flows is now covered by 2011 lava.

Description: A global map of surface heat flow is presented on a 2° by 2° equal area grid. It is based on a global heat flow data set of over 38,000 measurements. The map consists of three components. Firstly, in regions of young ocean crust (〈67.7Ma) the model estimate uses a half-space conduction model based on the age of the oceanic crust, since it is well known that raw data measurements are frequently influenced by significant hydrothermal circulation. Secondly in other regions of data coverage the estimate is based on data measurements. At the map resolution these two categories (young ocean, data covered) cover 65% of Earth's surface. Thirdly, for all other regions the estimate is based on the assumption that there is a correlation between heat-flow and geology. This assumption is assessed and the correlation is found to provide a minor improvement over assuming that heat flow would be represented by the global average. The map is made available digitally.

Description: We present 36 new 40 Ar- 39 Ar incremental heating age determinations from the Caribbean Large Igneous Province (CLIP) providing evidence for extended periods of volcanic activity and suggest a new tectonomagmatic model for the province's timing and construction. These new 40 Ar- 39 Ar ages for the Curaçao Lava Formation (CLF) and Haiti's Dumisseau Formation show evidence for active CLIP volcanism from 94 to 63 Ma. No clear changes in geochemical character are evident over this period. The CLF has trace element signatures (e.g., Zr/Nb = 10-20) and flat rare earth element (REE) trends consistent with plume volcanism. The Dumisseau Formation also has plume-like geochemistry and steeper REE trends similar to ocean island basalts. Volcanism in the Dumisseau Formation appears to have largely ceased by 83 Ma while at Curaçao it continued until 63 Ma. A rapidly surfacing and melting plume head alone does not fit this age distribution. Instead, we propose that the residual Galapagos plume head, following initial ocean plateau construction, was advected eastward by asthenospheric flow induced by subducting oceanic lithosphere. Slab rollback at the Lesser Antilles and Central America subduction zones created an extensional regime within the Caribbean plate. Mixing of plume with upwelling asthenospheric mantle provided a source for intermittent melting and eruption through the original plateau over a ~30 Ma period.

Description: To study how an impacting plume modifies the mantle lithosphere, we analyzed the microstructures and crystal preferred orientations (CPO) of 29 peridotites and 37 pyroxenites that sample the mantle root of the Ontong Java Plateau (OJP) from 60 to 120 km depth. The peridotites show a strong compositional variability, but homogeneous coarse-granular to tabular microstructures, except for those equilibrated at the shallowest and deepest depths, which are porphyroclastic. All peridotites have clear olivine CPO, with dominant fiber-[010] patterns. Low intragranular misorientations and straight grain boundaries in olivine suggest that, above 100 km depth, annealing often followed deformation. Calculated density and P-wave velocities of the peridotites decrease weakly with depth. S-wave velocities decrease faster, resulting in increasing Vp/Vs ratio with depth. Calculated densities and seismic velocity profiles are consistent with those estimated for normal mantle compositions under a cold oceanic geotherm. Enrichment in pyroxenites may further increase seismic velocities. The calculated seismic properties cannot therefore explain the low S-waves velocities predicted by Rayleigh wave tomography and ScS data in the mantle beneath the OJP. Calculated P- and S-wave anisotropy is variable (2-12%). It is higher on average in the deeper section of the lithosphere. Because olivine has dominantly [010]-fiber CPO patterns, if foliations are horizontal, vertically propagating S-waves and Rayleigh waves will sample very weak anisotropy in the OJP mantle lithosphere. Moreover, if the orientation of the lineation changes with depth, the anisotropy-induced contrast in seismic properties might produce an intralithospheric reflector marking the stratification of the OJP mantle root.

Description: The quantification of heat and mass flow between deep reservoirs and the surface is important for understanding magmatic and hydrothermal systems. Here, we use high-resolution measurement of carbon dioxide flux (ϕCO 2 ) and heat flow at the surface to characterize the mass (CO 2 and steam) and heat released to the atmosphere from two magma-hydrothermal systems. Our soil gas and heat flow surveys at Rotokawa and White Island in the Taupō Volcanic Zone, New Zealand, include over 3,000 direct measurements of ϕCO 2 and soil temperature and 60 carbon isotopic values on soil gases. Carbon dioxide flux was separated into background and magmatic/hydrothermal populations based on the measured values and isotopic characterization. Total CO 2 emission rates (ΣCO 2 ) of 441 ± 84 t d -1 and 124 ± 18 t d -1 were calculated for Rotokawa (2.9 km 2 ) and for the crater floor at White Island (0.3 km 2 ), respectively. The total CO 2 emissions differ from previously published values by +386 t d -1 at Rotokawa and +25 t d -1 at White Island, demonstrating that earlier research underestimated emissions by 700% (Rotokawa) and 25% (White Island). These differences suggest that soil CO 2 emissions facilitate more robust estimates of the thermal energy and mass flux in geothermal systems than traditional approaches. Combining the magmatic/hydrothermal-sourced CO 2 emission (constrained using stable isotopes) with reservoir H 2 O:CO 2 mass ratios and the enthalpy of evaporation, the surface expression of thermal energy release for the Rotokawa hydrothermal system (226 MW t ) is 10 times greater than the White Island crater floor (22.5 MW t ).

Description: To better constrain the mechanical behavior of sediments accreted to accretionary prism, we conducted triaxial mechanical tests on natural samples from the Miura-Boso paleo-accretionary prism (Japan) in drained conditions with confining pressures up to 200 MPa as well as post-experiments P-wave velocity (V p ) measurements. During experiments, deformation is principally non-coaxial and accommodated by two successive modes of deformation, both associated with strain-hardening and velocity-strengthening behavior: (1) compaction-assisted shearing, distributed in a several mm-wide shear zone and (2) faulting, localized within a few tens of µm-wide, dilatant fault zone. Deformation is also associated with (1) a decrease in Young's modulus all over the tests, (2) anomalously low V p in the deformed samples compared to their porosity and (3) an increase in sensitivity of V p to effective pressure. We interpret this evolution of the poroelastic properties of the material as reflecting the progressive breakage of intergrain cement and the formation of microcracks along with macroscopic deformation. When applied to natural conditions, these results suggest that the deformation style (localized vs distributed) of shallow (z 〈 a few km) sediments is mainly controlled by the variations in stress/strain rate during the seismic cycle and is therefore independent of the porosity of sediments. Finally, we show that the effect of strain, through cement breakage and microcracks formation, may lower V p for effective pressure up to 40 MPa. As a consequence, the low V p anomalies observed in Nankai accretionary prisms by seismic imaging between 2 and 4km depth [ Kitajima and Saffer , 2012] could reflect sediment deformation rather than porosity anomalies. This article is protected by copyright. All rights reserved.

Description: This study presents a new thermodynamic model for the calculation of phase relations during the melting of anhydrous spinel lherzolite at pressures of 1–2.5 GPa. The model is based on the total energy minimization algorithm for calculating phase equilibria within multicomponent systems and the thermodynamic configuration of Ueki and Iwamori [2013]. The model is based on a SiO 2 – Al 2 O 3 – FeO–Fe 3 O 4 –MgO–CaO system that includes silicate melt, olivine, clinopyroxene, orthopyroxene, and spinel as possible phases. The molar Gibbs free energy of the melt phase is modeled quasi-empirically, and the thermodynamic parameters for silicate melt end-member components are calibrated with a polybaric calibration database. The temperatures and pressures used in this newly compiled calibration dataset are 1230–1600 ∘ C and 0.9–3 GPa, corresponding to the stability range of spinel lherzolite. The modeling undertaken during this study reproduces the general features of experimentally determined melting phase relations of spinel lherzolite at 1–2.5 GPa, including the solidus temperature, the melt composition, the chemical reaction during melting and the degree of melting. This new thermodynamic modeling also reproduces phase relations of various bulk compositions from fertile to deplete spinel lherzolite and can be used in the modeling of polybaric mantle melting within various natural settings. Comparing the results derived from this new modeling with those produced using previous models indicates that the new approach outlined here, involving a combination of total energy minimization and the direct calibration of melt thermodynamic parameters at pressure and temperature conditions corresponding to mantle melting with a relatively simple melt thermodynamic equation, can accurately model polybaric melting phase relations. This article is protected by copyright. All rights reserved.

Description: The volatile contents of olivine-hosted (Fo89–71) melt inclusion glasses in rapidly quenched mafic tephras from volcanic front volcanoes of the Central American Volcanic Arc (CAVA) in Guatemala, Nicaragua, and Costa Rica, were analyzed by secondary ion mass spectrometry (SIMS) in order to derive the minimum eruptive output of CO2, along with H2O, Cl, and S. Details of the analytical method are provided that establish melt inclusion CO2 analyses with the Cameca ims6f at the Helmholtz Centre Potsdam. The highest CO2 concentrations (up to 1800 μg/g) are observed in Nicaraguan samples, while melt inclusions from Guatemala and Costa Rica have CO2 contents between 50 and 500 μg/g. CO2 does not positively covary with sediment/slab fluid tracers such as Ba/La, Ba/Th, or U/La. Instead, the highest CO2 concentrations occur in the inclusions with the most depleted incompatible element compositions and low H2O, approaching the composition of mid-ocean ridge basalts (MORBs), whereas the most H2O-rich inclusions are relatively CO2-poor (

Description: We determined the focal mechanism solutions (FMS) of 191 crustal earthquakes as well as the stress tensor in the source area of the 2008 Iwate-Miyagi earthquake (2008 IMEQ, M7.2) that occurred in the central portion of northeast (NE) Japan. The FMS and the stress tensors were determined by using both 1-D and 3-D velocity models, which exhibit almost the same results. The differences caused by the use of 1-D and 3-D models can be neglected when compared with the differences due to the different methods, which indicates that the FMS and the stress tensor determined with a 1-D model are accurate enough to study the crustal stress field in the study region. The obtained P axis (σ1) trends WNW-ESE subhorizontally, and the T axis (σ3) is oriented subvertically in a NNE-SSW belt perpendicular to σ1. The σ1 orientation is consistent with the motion of the Pacific plate relative to NE Japan, which indicates that the plate boundary forces dominate the intraplate stress regime. Both temporal and spatial variations of the stress field in the IMEQ source area are detected, which may be induced by the stress rotation accompanying the main shock and its aftershocks. The seismogenic faults in the study area are estimated to be very weak, which argues against the concept of strong crust. The faults may be weakened by the high-temperature magma and the fluids in the lower crust and uppermost mantle that intrude upward into the shallower crust.

Description: The Integrated Ocean Drilling Program (IODP) Hole 1301A on the eastern flank of Juan de Fuca Ridge was used in the first long-term deployment of microbial enrichment flow cells using osmotically driven pumps in a subseafloor borehole. Three novel osmotically driven colonization systems with unidirectional flow were deployed in the borehole and incubated for 4 years to determine the microbial colonization preferences for 12 minerals and glasses present in igneous rocks. Following recovery of the colonization systems, we measured cell density on the minerals and glasses by fluorescent staining and direct counting and found some significant differences between mineral samples. We also determined the abundance of mesophilic and thermophilic culturable organotrophs grown on marine R2A medium and identified isolates by partial 16S or 18S rDNA sequencing. We found that nine distinct phylotypes of culturable mesophilic oligotrophs were present on the minerals and glasses and that eight of the nine can reduce nitrate and oxidize iron. Fe(II)-rich olivine minerals had the highest density of total countable cells and culturable organotrophic mesophiles, as well as the only culturable organotrophic thermophiles. These results suggest that olivine (a common igneous mineral) in seawater-recharged ocean crust is capable of supporting microbial communities, that iron oxidation and nitrate reduction may be important physiological characteristics of ocean crust microbes, and that heterogeneously distributed minerals in marine igneous rocks likely influence the distribution of microbial communities in the ocean crust.

Description: We present the first comprehensive study of mass wasting processes in the continental slope of a convergent margin of a subduction zone where tectonic processes are dominated by subduction erosion. We have used multibeam bathymetry along ∼1300 km of the Middle America Trench of the Central America Subduction Zone and deep-towed side-scan sonar data. We found abundant evidence of large-scale slope failures that were mostly previously unmapped. The features are classified into a variety of slope failure types, creating an inventory of 147 slope failure structures. Their type distribution and abundance define a segmentation of the continental slope in six sectors. The segmentation in slope stability processes does not appear to be related to slope preconditioning due to changes in physical properties of sediment, presence/absence of gas hydrates, or apparent changes in the hydrogeological system. The segmentation appears to be better explained by changes in slope preconditioning due to variations in tectonic processes. The region is an optimal setting to study how tectonic processes related to variations in intensity of subduction erosion and changes in relief of the underthrusting plate affect mass wasting processes of the continental slope. The largest slope failures occur offshore Costa Rica. There, subducting ridges and seamounts produce failures with up to hundreds of meters high headwalls, with detachment planes that penetrate deep into the continental margin, in some cases reaching the plate boundary. Offshore northern Costa Rica a smooth oceanic seafloor underthrusts the least disturbed continental slope. Offshore Nicaragua, the ocean plate is ornamented with smaller seamounts and horst and graben topography of variable intensity. Here mass wasting structures are numerous and comparatively smaller, but when combined, they affect a large part of the margin segment. Farther north, offshore El Salvador and Guatemala the downgoing plate has no large seamounts but well-defined horst and graben topography. Off El Salvador slope failure is least developed and mainly occurs in the uppermost continental slope at canyon walls. Off Guatemala mass wasting is abundant and possibly related to normal faulting across the slope. Collapse in the wake of subducting ocean plate topography is a likely failure trigger of slumps. Rapid oversteepening above subducting relief may trigger translational slides in the middle Nicaraguan upper Costa Rican slope. Earthquake shaking may be a trigger, but we interpret that slope failure rate is lower than recurrence time of large earthquakes in the region. Generally, our analysis indicates that the importance of mass wasting processes in the evolution of margins dominated by subduction erosion and its role in sediment dynamics may have been previously underestimated.

Description: In the past decade, major advancements in precision and accuracy of U-Pb geochronology, which stem from improved sample pretreatment and refined measurement techniques, have revealed previously unresolvable discrepancies among analyses from different laboratories. One solution to evaluating and resolving many of these discrepancies is the adoption of a common software platform that standardizes data-processing protocols, enabling robust interlaboratory comparisons. We present the results of a collaboration to develop cyber infrastructure for high-precision U-Pb geochronology based on analyzing accessory minerals by isotope dilution-thermal ionization mass spectrometry. This cyber infrastructure implements an architecture specifying the workflows of data acquisition, statistical filtering, analysis and interpretation, publication, community-based archiving, and the compilation and comparison of data from different laboratories. The backbone of the cyber infrastructure consists of two open-source software programs: Tripoli and U-Pb_Redux. Tripoli interfaces with commercially available mass spectrometers using standardized protocols, statistical filtering, and interactive visualizations to aid the analyst in preparing raw data for analysis in U-Pb_Redux. U-Pb_Redux implements the architecture by orchestrating the analyst's workflow with interactive visualizations and provides data reduction and uncertainty propagation that support data interpretations. Finally, U-Pb_Redux enables production of publication-ready graphics and data tables, the archiving of results, and the comparative compilation of archived results to support cooperative science.

Description: The late Holocene activity of a restraining bend of the northern strand of the North Anatolian Fault in Izmit Bay was investigated by a sedimentological, geochemical, and paleoecological analysis of sediment cores from Hersek coastal lagoon, NW Turkey. The sediment cores show a succession of sedimentary sequences composed of three units separated by gradual transitions. The first unit is composed of a thin layer of shell debris-rich sediment in abrupt contact with the underlying organic-rich deposits. This unit is overlain by a thick foraminifera-rich mud deposit, and the sequences are capped by an organic-rich mud unit. These sequences are interpreted as silting up, shallowing upward deposits, typical of a lagoon becoming isolated from the sea. We suggest that they represent the sedimentary signature of coseismic subsidence, which was caused by reverse slip at the Hersek bend, and tsunamis in Izmit Bay. Our radiocarbon-dated paleoseismological record indicates (1) the atypical collapse of the hanging wall during the 740 earthquake and (2) subsidence of the footwall during the 987, 1509, and 1719 earthquakes. This study contributes to the understanding of the dynamics of restraining bends, and it highlights the potential of coastal sediments for reconstructing past earthquakes and tsunamis in regions dominated by strike-slip deformations.

Description: Estimates of slab curvature from 133 profiles across 12 major subduction zones are obtained by fitting splines through earthquakes locations. Several geometric quantities are computed from the spline fits, including the minimum radius of curvature. We also introduce a quantity called the equivalent radius of curvature to precisely define the bending dissipation in either a viscous or a plastic slab. The minimum and equivalent radius of curvature are approximately equal in slabs with simple geometry (i.e., a single bend) but deviate substantially when the geometry is more complicated. A compilation of estimates shows that the equivalent radius of curvature is insensitive to the age of the plate at the trench. Mean values are 175 km or 191 km, depending on the rheology of the slab. When bending occurs by viscous deformation, the resulting dissipation can exceed the power available from subduction with reasonable choices for viscosity. A more consistent description of the energetics is obtained by assuming a plastic rheology.

Description: New geochemical and microstructural data constrain the origins and pathways of paleofluids during the growth of Sheep Mountain Anticline, Wyoming, United States. Oxygen, carbon, and strontium isotope studies were performed on prefolding and fold-related calcite veins and their sedimentary host rocks and combined to fluid inclusion microthermometry results. We show that most of the cements precipitated from Paleogene meteoric fluid. Stable isotopes and fluid inclusion homogenization temperatures further indicate that most veins were mineralized from upward moving fluids after these fluids were heated at depth (T 〉 110°C). This implies that fluids migrated along the basement thrust underlying the fold and/or at the base of the cover. Above the fault tip, the fluids circulated rapidly in the diffuse synfolding (and early folding) fracture network. The zone of preferential migration of the warm fluids is currently located in the backlimb of the fold, which supports some of the previously published structural interpretation of the subsurface. This study also highlights the potential of combined fracture analysis and geochemical analyses of paleofluid flows in fractures to constrain both the deformation history and the fluid flow during basement-involved shortening in Laramide-style forelands.

Description: Basalts from the Mid-Atlantic Ridge change progressively in composition with increasing distance from the Azores platform. Study of the Lucky Strike and Menez Gwen segments reveals much complexity in the gradient. Both segments contain only basalts enriched relative to normal mid-oceanic ridge basalt, but in two distinct groups. Moderately enriched basalts occur throughout the segments, with proximal Menez Gwen enriched relative to Lucky Strike. Highly enriched basalts occur at segment centers. Incompatible element ratios of the highly enriched basalts exceed those of the Azores platform, while isotopic compositions are less enriched. These observations can be explained by a low-degree melt of garnet-bearing Azores mantle added to mantle depleted by previous melt extraction. Melting this “metasomatized” mantle produces lavas that match the enriched samples. The Azores gradient cannot be explained by simple two-component mixing; rather, it reflects recent melt extraction and addition processes related to southward flow of the Azores plume. The Azores gradient also permits tests of segmentation models. Central supply models predict step functions in chemical compositions between segments. Within-segment gradients require vertical supply. Central supply is supported by robust central volcanoes, thicker crust at segment centers, and a step function in isotopes between the segments. The lava diversity at segment centers, however, requires batches of distinct magma that are preserved through melting and melt delivery. Within-segment gradients in moderately incompatible element ratios support a component of multiple supply. The data suggest partial homogenization of magma within a segment and preferential melt focusing to segment centers with some vertical transport.

Description: The nitrogen isotopic composition of organic matter contained within diatom frustules of different sediment size fractions from the Guaymas Basin, Gulf of California, suggests that multiple nutrient cycling processes can be examined from a single sediment sample. The diatom-bound δ15N values from the greater and less than 63 μm size fractions were analyzed from six intervals of core MD02–2517, each representing ∼30–40 years. The diatom-bound δ15N values are compared to bulk sedimentary δ15N values, weight percent biogenic opal, and weight percent organic carbon. Bulk sedimentary δ15N values range from 10‰ to 12.5‰ and diatom-bound δ15N values of the 63 μm fraction range from 1‰ to 7‰. The best explanation for the low diatom δ15N values in the 〉63 μm fraction is that the large diatoms grew at the base of the euphotic zone in excess nitrate and under low light conditions. Large diatom δ15N values correlate significantly with bulk sedimentary δ15N values (r = 0.68, 0.82), indicating that large diatoms, such as Thalassiothrix longissima, are an important component of export production fueled by a nonupwelling source of nitrogen. The 〉63 μm diatoms are estimated to contribute between 3% and 15% of the total nitrogen in the sediments based upon an assumed Si-to-N ratio of 1.2.

Description: Episodic tremor and slip (ETS) has been observed in many subduction zones, but its mechanical underpinnings as well as its potential for triggering damaging earthquakes have proven difficult to assess. Here we use a seismic array in Cascadia of unprecedented density to monitor seismicity around a moderate 16 day ETS episode. In the 4 months of data we examine, we observe five tiny earthquakes within the subducting slab during the episode and only one more in the same area, which was just before and nearby the next ETS burst. These earthquakes concentrate along the sides and updip edge of the ETS region, consistent with greater stress concentration there than near the middle and downdip edge of the tremor area. Most of the seismicity is below the megathrust, with a similar depth extent to the background intraslab seismicity. The pattern of earthquakes that we find suggests slow slip has a more continuous temporal and spatial pattern than the tremor loci, which notoriously appear in bursts, jumps, and streaks.

Description: Dust samples collected from subway platforms in Shanghai, China, have been examined using magnetic measurements and geochemical analysis. Our results indicate that the Shanghai subway platform dusts have extremely strong magnetic signatures. These results, combined with X-ray diffraction analysis and scanning and transmission electron microscope examinations, indicate that the magnetic mineralogy of the dust is dominated by iron scraps due to wheel-rail mechanical abrasion and spherules rich in magnetite from fossil fuel combustion. Although the magnetic particles are primarily micrometer sized, fine submicron magnetic grains are also evident in the dust. The underground platform dusts have a much higher iron flake abundance and magnetic susceptibility than those from aboveground platforms because the latter ones are diluted by inputs of magnetically weaker ambient aerosols with a higher proportion of magnetite spherules. Geochemical analysis indicates that underground platform dusts have elevated Fe and Mn, but lower Al and Ti contents relative to aboveground subway dust. This is consistent with the closed nature of underground platforms, which therefore reduces exposure to soil-derived dust. Since the adverse environmental effects of subway particles may be linked to higher contents of iron and other metals, our results demonstrate that magnetic measurements provide a novel and effective approach for characterizing iron mineralogy and grain size in subway dusts.

Description: The Cocos plate experiences extensional faulting as it bends into the Middle American Trench (MAT) west of Nicaragua, which may lead to hydration of the subducting mantle. To estimate the along strike variations of volatile input from the Cocos plate into the subduction zone, we gathered marine seismic refraction data with the R/V Marcus Langseth along a 396 km long trench parallel transect offshore of Nicaragua and Costa Rica. Our inversion of crustal and mantle seismic phases shows two notable features in the deep structure of the Cocos plate: (1) Normal oceanic crust of 6 km thickness from the East Pacific Rise (EPR) lies offshore Nicaragua, but offshore central Costa Rica we find oceanic crust from the northern flank of the Cocos Nazca (CN) spreading center with more complex seismic velocity structure and a thickness of 10 km. We attribute the unusual seismic structure offshore Costa Rica to the midplate volcanism in the vicinity of the Galápagos hot spot. (2) A decrease in Cocos plate mantle seismic velocities from ∼7.9 km/s offshore Nicoya Peninsula to ∼6.9 km/s offshore central Nicaragua correlates well with the northward increase in the degree of crustal faulting outboard of the MAT. The negative seismic velocity anomaly reaches a depth of ∼12 km beneath the Moho offshore Nicaragua, which suggests that larger amounts of water are stored deep in the subducting mantle lithosphere than previously thought. If most of the mantle low velocity zone can be interpreted as serpentinization, the amount of water stored in the Cocos plate offshore central Nicaragua may be about 2.5 times larger than offshore Nicoya Peninsula. Hydration of oceanic lithosphere at deep sea trenches may be the most important mechanism for the transfer of aqueous fluids to volcanic arcs and the deeper mantle.

Description: The Cenozoic Song Hong–Yinggehai Basin in the South China Sea contains a large volume of sediment that has been used in previous studies, together with regional geomorphology, to argue for the existence of a large palaeodrainage system that connected eastern Tibet with the South China Sea. To test this and to understand the significance of sediment volumes deposited in the Song Hong–Yinggehai Basin, this study compared erosion histories of source regions with sediment volumes deposited during the two main stages in basin evolution spanning active rifting and subsidence (30–15.5 Ma) and postrift sedimentation (15.5 Ma to present). The study of basin provenance by detrital zircon U-Pb dating revealed Hainan was an important and continuous source of sediment, and a bedrock thermochronological study quantified its overall contribution to basin sedimentation. Comparison between the accumulated mass of basin sediment and volumes of eroded bedrock, calculated from apatite thermochronometry across the modern Red River drainage in northern Vietnam as well as Hainan Island, accounted for the bulk of sediment deposited since 30 Ma. Consequently, if an expanded paleodrainage ever existed it must have predated the Oligocene.

Description: Quartz mylonites from the Tonale Fault Zone in the Alps (northern Italy) have been investigated by the Ti-in-quartz geothermometer (TitaniQ) in order to test its applicability to measure deformation temperatures. The eastern part of the Tonale Fault Zone was contact metamorphosed by the synkinematic intrusion of the Adamello pluton, forming an ∼800 m wide mylonitic shear zone, with a synkinematic temperature gradient from ∼280°C at the frictional-viscous transition to ∼700°C at the pluton contact as derived from metamorphic mineral assemblages. Deformation microstructures from quartz mylonite samples, systematically collected across the mylonitic shear zone, display the entire range of dynamic recrystallization in quartz, which comprise bulging recrystallization (BLG), subgrain rotation recrystallization (SGR), and grain boundary migration recrystallization (GBM). TitaniQ geothermometry yields the near-peak deformation temperature for quartz mylonites deformed at metamorphic temperatures above ∼540°C in the zone of GBM. However, for mylonites formed under lower temperatures in the zones of SGR and BLG, the preexisting Ti concentrations were not reset. It is suggested that this is due to the sluggish Ti volume diffusion rates below 500°C and the short duration of contact metamorphism and deformation. Even in the higher temperature samples the reequilibration of Ti-in-quartz content was achieved by grain boundary migration rather than by volume diffusion. Hence, our results show that GBM is crucial for the reequilibration of Ti-in-quartz, while quartz mylonites deformed by either BLG or SGR, which predominate in natural shear zones at greenschist facies metamorphic conditions, most likely yield inherited temperatures.

Description: High-precision U-Pb geochronology by isotope dilution-thermal ionization mass spectrometry is integral to a variety of Earth science disciplines, but its ultimate resolving power is quantified by the uncertainties of calculated U-Pb dates. As analytical techniques have advanced, formerly small sources of uncertainty are increasingly important, and thus previous simplifications for data reduction and uncertainty propagation are no longer valid. Although notable previous efforts have treated propagation of correlated uncertainties for the U-Pb system, the equations, uncertainties, and correlations have been limited in number and subject to simplification during propagation through intermediary calculations. We derive and present a transparent U-Pb data reduction algorithm that transforms raw isotopic data and measured or assumed laboratory parameters into the isotopic ratios and dates geochronologists interpret without making assumptions about the relative size of sample components. To propagate uncertainties and their correlations, we describe, in detail, a linear algebraic algorithm that incorporates all input uncertainties and correlations without limiting or simplifying covariance terms to propagate them though intermediate calculations. Finally, a weighted mean algorithm is presented that utilizes matrix elements from the uncertainty propagation algorithm to propagate random and systematic uncertainties for data comparison between other U-Pb labs and other geochronometers. The linear uncertainty propagation algorithms are verified with Monte Carlo simulations of several typical analyses. We propose that our algorithms be considered by the community for implementation to improve the collaborative science envisioned by the EARTHTIME initiative.

Description: Orogeny in the geologically young northern Apennines is explained by the eastward retreat of a convergence zone in which Adriatic lithosphere subducts to the west. It is unclear, however, whether all of the lithosphere or only the lower portion has subducted since 15 Ma, when the convergence zone lay near Corsica. We combine teleseismic P and S wave arrival time data from the RETREAT seismic network (2003–2006) and surrounding permanent stations to estimate tomographic images of the upper mantle structure beneath the northern Apennines. We image a vertically oriented slab that that extends to only ∼300 km depth. Our slab termination is shallower than previous studies but is confirmed by resolution tests. Furthermore, our images resolve the southern edge of the northern Apennines slab at ∼43°N, with no deep continuity with any slab segment to the south, as earlier proposed. Our results suggest that only a 300–400 km strip of lithosphere has subducted since 15 Ma beneath the northern Tyrrhenian Sea, a length roughly equivalent to the distance from the present-day west coast of Corsica to the crest of the northern Apennines. Although not a definitive indicator, the largely aseismic vertical slab configuration and its limited extent, coupled with other volcanic and geophysical indicators, suggests the delamination scenario of Bird (1979), in which only the lower portion of the continental lithosphere subducts.

Description: Over the past several decades, contrasting models have been proposed for the physical and chemical processes responsible for the uplift and long-term stability of the Colorado Plateau (CP) and crustal thinning beneath the Basin and Range Province (BRP) in the southwestern United States. Here we provide new constraints on the models by modeling gravity anomalies and by systematically analyzing over 15,500 P-to-S receiver functions recorded at 72 USArray and other broadband seismic stations on the southwestern CP and the southern BRP. Our results reveal that the BRP is characterized by a thin crust (28.2 ± 0.5 km), a mean Vp/Vs of 1.761 ± 0.014 and a mean amplitude (R) of P-to-S converted wave (relative to that of the direct P wave) of 0.181 ± 0.014 that are similar to a typical continental crust, consistent with the model that the thin crust was the consequence of lithospheric stretching during the Cenozoic. The CP is characterized by the thickest crust (42.3 ± 0.8 km), largest Vp/Vs (1.825 ± 0.009) and smallest R (0.105 ± 0.007) values in the study area. In addition, many stations on the CP exhibit a clear arrival before the P-to-S converted phase from the Moho, corresponding to a lower crustal layer of about 12 km thick with a mafic composition. We hypothesize that the lower crustal layer, which has an anomalously large density as revealed by gravity modeling and high velocities in seismic refraction lines, contributed to the long-term stability and preuplift low elevation of the Colorado Plateau.

Description: We examined the effects of high pressure on thermal conductivity in core samples from the slope–apron facies and the upper part of the accretionary prism at site C0001 of the NanTroSEIZE drilling program and in other samples of five terrestrial rock types. Thermal conductivity clearly increased with increasing pressure for both wet (water saturated) and dry samples. We determined the rate of thermal conductivity change of the NanTroSEIZE sediments to be 0.014 Wm−1K−1/MPa when pressure was increased, and 0.01 Wm−1K−1/MPa when pressure was decreased. Using the rate determined for decreasing pressure, we estimated that thermal conductivities measured at atmospheric pressure rather than at in situ pressure may be underestimated by 7% for a core sample from around 1 km depth and by 20% for a core sample from around 3 km depth. In general, the rate of thermal conductivity change with pressure showed a positive correlation with porosity. However, the relationship of the rate of thermal conductivity change to porosity is also dependent on the fabric, mineral composition, and pore structure of the sediments and rocks. Furthermore, for two sandstones we tested, the effect of pressure on thermal conductivity for dry samples was greater than that for wet samples.

Description: A continuous section from extrusive lavas, through sheeted dikes, and uppermost gabbros recovered from Integrated Ocean Drilling Program Hole 1256D provides important information regarding magma plumbing systems beneath superfast spreading ridges. Petrological examination demonstrates that a model of fractional crystallization from a magma of composition similar to one of the more primitive gabbros in a shallow (∼50–100 MPa) melt lens reasonably explains mineral and whole rock compositions of many lavas and dikes. Elevated concentrations of trace elements in some rocks appear to have resulted from mixing between primitive magma and highly evolved magma. About half of the dike samples have more evolved Fe-rich compositions than the extrusive lava samples. Magma densities of the Fe-rich dikes are a little higher (∼30 kg/m3) than those of lavas, suggesting that these dike magmas would not reach the surface. Mineralogical investigations reveal that both lavas and dikes contain oscillatory zoned plagioclase xenocrysts, implying magma mixing caused by successive episodes of fractionation and magma replenishment in the melt lens. The plagioclase xenocrysts contain high-Anorthite sections [An: 100 × Ca/(Ca+Na) in mole percent] whose compositions are not in equilibrium with host liquids. The high-An sections were likely crystallized when primitive magmas with high CaO/Na2O were injected in the melt lens. Since the oscillatory zoned plagioclase generally forms crystal clots, they were probably accumulated in a mush zone. The petrographical examination favors a model suggesting that injection of primitive magma into the melt lens broke the mush zone and pushed out the oscillatory zoned plagioclase.

Description: Analyses of radiogenic neodymium (Nd), strontium (Sr), and lead (Pb) isotope compositions of clay-sized detrital sediments allow detailed tracing of source areas of sediment supply and present and past transport of particles by water masses in the eastern Indian Ocean. Isotope signatures in surface sediments range from −21.5 ($\varepsilon$Nd), 0.8299 (87Sr/86Sr), and 19.89 (206Pb/204Pb) off northwest Australia to +0.7 ($\varepsilon$Nd), 0.7069 (87Sr/86Sr), and 17.44 (206Pb/204Pb) southwest of Java. The radiogenic isotope signatures primarily reflect petrographic characteristics of the surrounding continental bedrocks but are also influenced by weathering-induced grain size effects of Pb and Sr isotope systems with superimposed features that are caused by current transport of clay-sized particles, as evidenced off Australia where a peculiar isotopic signature characterizes sediments underlying the southward flowing Leeuwin Current and the northward flowing West Australian Current (WAC). Gravity core FR10/95-GC17 off west Australia recorded a major isotopic change from Last Glacial Maximum values of −10 ($\varepsilon$Nd), 0.745 (87Sr/86Sr), and 18.8 (206Pb/204Pb) to Holocene values of −22 ($\varepsilon$Nd), 0.8 (87Sr/86Sr), and 19.3 (206Pb/204Pb), which documents major climatically driven changes of the WAC and in local riverine particle supply from Australia during the past 20 kyr. In contrast, gravity core FR10/95-GC5 located below the present-day pathway of the Indonesian throughflow (ITF) shows a much smaller isotopic variability, indicating a relatively stable ITF hydrography over most of the past 92 kyr. Only the surface sediments differ significantly in their isotopic composition, indicating substantial changes in erosional sources attributed to a change of the current regime during the past 5 kyr.

Description: The Amphitheater Mountains and southern central Alaska Range expose a thick sequence of Triassic Nikolai basalts that is underlain by several mafic-ultramafic complexes, the largest and best exposed being the Fish Lake and Tangle (FL-T) mafic-ultramafic sills that flank the Amphitheater Mountains synform. Three-dimensional (3-D) modeling of gravity and magnetic data reveals details of the structure of the Amphitheater Mountains, such as the orientation and thickness of Nikolai basalts, and the geometry of the FL-T intrusions. The 3-D model (50 × 70 km) includes the full geographic extent of the FL-T complexes and consists of 11 layers. Layer surfaces and properties (density and magnetic susceptibility) were modified by forward and inverse methods to reduce differences between the observed and calculated gravity and magnetic grids. The model suggests that the outcropping FL-T sills are apparently connected and traceable at depth and reveals variations in thickness, shape, and orientation of the ultramafic bodies that may identify paths of magma flow. The model shows that a significant volume (2000 km3) of ultramafic material occurs in the subsurface, gradually thickening and plunging westward to depths exceeding 4 km. This deep ultramafic material is interpreted as the top of a keel or root system that supplied magma to the Nikolai lavas and controlled emplacement of related magmatic intrusions. The presence of this deep, keel-like structure, and asymmetry of the synform, supports a sag basin model for development of the Amphitheater Mountains structure and reveals that the feeders to the Nikolai are much more extensive than previously known.

Description: Multiple lines of evidence exist for a range of sediment mass movement processes within the shallow megasplay fault zone (MSFZ) area and the adjacent slope basin in the outer fore arc of the Nankai subduction zone, Japan. Diagnostic features observed in three-dimensional reflection seismic data and in cores of the Integrated Ocean Drilling Program (IODP) document a multifarious mass movement history spanning ∼2.87 million years. Various modes and scales of sediment remobilization can be related to the different morphotectonic settings in which they occurred. From this evidence, we decipher the tectonic control on slumping and mass transport deposition in the Nankai fore arc. Three periods of intensified mass wasting coincided with pulses of enhanced activity on the splay fault: (1) an initial phase of juvenile out-of-sequence thrusting ∼1.95 to 1.7 Ma, (2) a reactivation phase between ∼1.55 and 1.24 Ma, and (3) at about 1 Ma, during a phase of uplift of the fore-arc high and motion along the MSFZ. We suggest that slope oversteepening, extensional stress regimes, and lateral transmission of fluid overpressures may have preconditioned the slope sediments to fail. Individual mass-wasting events may have been triggered by dynamic loading from earthquake waves and/or transient pulses of pore pressure along the splay fault. Overall, our results provide insights into the complicated interplay between tectonic and submarine mass movement processes. We demonstrate that detailed knowledge about the spatial and temporal distribution of submarine mass movements can be integrated into a holistic reconstruction of tectonostratigraphic evolution of accretionary margins.

Description: Experimental studies indicate that crystal-bearing magma exhibits non-Newtonian behavior at high strain rates and solid fractions. We use a zero-dimensional (0-D) inversion model to reevaluate rheological parameters and shear heating effects from laboratory data on crystal-bearing magma. The results indicate non-Newtonian behavior with power law coefficients of up to n = 13.5. It has been speculated that finite strain effects, shear heating, power law melt rheology, or plasticity are responsible for this non-Newtonian behavior. We use 2-D direct numerical crystal-scale simulations to study the relative importance of these mechanisms. These simulations demonstrate that shear heating has little effect on aggregate (bulk) rheologies. Finite strain effects result in both strain weakening and hardening, but the resulting power law coefficient is modest (maximum n = 1.3). For simulations with spherical crystals the strain weakening and hardening behavior is related to rearrangement of crystals rather than strain rate related weakening. Finite strain effects were insignificant in a numerical simulation with naturally shaped crystals. Strain partitioning into the melt phase may induce microscopic stresses that are adequate to provoke a nonlinear viscous response in the melt. Large differential stresses and low effective stresses revealed by the simulations are sufficient to cause crystals to fail plastically. Numerical experiments that account for plastic failure show large power law coefficients (n ≈ 50 in some simulations). We conclude that this effect is the dominant cause of the strong nonlinear viscous response of crystal-bearing magmas observed in laboratory experiments.

Description: Postspreading volcanism has built large seamounts and volcanic ridges along the short axes of a highly segmented part of the East Pacific Rise crest that ceased spreading at the end of the middle Miocene, offshore Baja California Sur, Mexico. Lava samples from Rosa Seamount, the largest volcano, are predominantly alkalic basalts, mugearites, and benmoreites. This lavas series was generated through fractional crystallization and is compositionally similar to the moderately alkalic lava series in many oceanic islands. Samples from volcanic ridges at three adjacent failed spreading axes include mildly alkalic, transitional, and tholeiitic basalts and differentiated trachyandesites and andesite. The subtle but distinct petrologic and isotopic differences among the four sites may be due to differences in the degree of partial melting of a common, heterogeneous source. Postspreading lavas from these four abandoned axes off Baja California Sur together with those from other fossil spreading axes and from seamount volcanoes that grew on the East Pacific Rise flanks define a compositional continuum ranging from normal mid-ocean ridge basalt (NMORB)-like to ocean island basalt (OIB)-like. We propose that the compositional spectrum of these intraplate volcanic lavas is due to different degrees of partial melting of the compositionally heterogeneous suboceanic mantle in the eastern Pacific. A large degree of partial melting of this heterogeneous mantle during vigorous mantle upwelling at an active spreading center produces NMORB melts, whereas a lesser degree of partial melting during weak mantle upwelling following cessation of spreading produces OIB-like melts. The latter melts have a low (

Description: Here we present high-precision (2σ = ±3 ppm) 87Sr/86Sr measurements in coastal waters, together with salinity, to evaluate water mass mixing and the influence of submarine groundwater discharge (SGD) in coastal waters and marginal seas. Nonhomogeneous Sr isotopic variations in water columns were documented in the Southern Okinawa Trough (SOT), South China Sea, and Kao-ping Canyon (KPC), where seawater 87Sr/86Sr varied up to 70 ppm. Seawater Sr isotopic composition changes only slightly in the upper 200 m of the SOT but was detectable and highly correlated with salinity, indicating a mixing between radiogenic North Pacific Tropical Water (high 87Sr/86Sr and high salinity) at 100–150 m and a less radiogenic component with low 87Sr/86Sr and low salinity at ∼200 m. Vertical profiles of seawater 87Sr/86Sr along the KPC show significant variations, suggesting dynamic mixing affected by continental inputs (i.e., river runoff and SGD) in this region. These results highlight the potential use of seawater Sr isotopes as a powerful tracer for determining mixing ratios and the dynamic mixing of oceanic water masses, especially in coastal and marginal seas.

Description: Gravity-derived crustal thickness models were calculated for the North Atlantic Ocean between 76°N and the Chain Fracture Zone and calibrated using seismically determined crustal thickness. About 7% of the ocean crust is 7 km thick and is interpreted to have been affected by excess magmatism. Thin crust probably reflects reduced melt production from relatively cold or refractory mantle at scales of up to hundreds of kilometers along the spreading axis. By far the most prominent thick crust anomaly is associated with Iceland and adjacent areas, which accounts for 57% of total crustal volume in excess of 7 km. Much smaller anomalies include the Azores (8%), Cape Verde Islands (6%), Canary Islands (5%), Madeira (

Description: The central North Island of New Zealand is characterized by extension, high heat flow, and active volcanism caused by the Pacific plate subducting beneath the Australian plate off the east coast of the North Island. The Taupo Volcanic Zone (TVZ) is the current locus of rifting and active volcanism. The shear velocity structure derived here from short-period surface wave analysis fills a gap in the geophysical information obtained in previous active and passive source studies. Using ambient noise correlation techniques, we reprocess data from four three-component temporary seismic arrays that were deployed between 2001 and 2005 in the central North Island and data from permanent seismic stations operational at the same time. Special attention has been paid to the correction of timing errors and incomplete instrument response information. Low shear velocities in the top 15–20 km of the crust inferred from Rayleigh wave dispersion analysis coincide with the presumed source regions of rhyolitic magma in the TVZ and are correlated with conductive bodies inferred from a 3-D magnetotelluric survey. Based on the observed changes in shear velocity we use empirical relations to infer an average percentage of partial melt in the upper crust of

Description: We developed a Keyhole Markup Language (KML) generator for converting geochemical and isotopic data sets for rocks and sediments stored in database systems into KML. Our program allows users to visualize geochemical or isotopic data easily in Google Earth. The generator accepts data files produced by the database systems PetDB, SedDB, GEOROC, and GANSEKI. The data are plotted three-dimensionally as a bar graph on the surface of the virtual Earth at the sampling site. This type of visual presentation, including information on sample localities, directly shows the distribution of isotopic or compositional anomalies of specific samples on the Earth's surface. We provide a Web application for the generator, so anyone can set the parameters for visualization over the Internet. With other KML generators we developed earlier, geochemical data can be overlain on a seismic tomographic model. This overlay image can provide information on the origin of samples in the tomographic model.

Description: Many studies have suggested that continental weathering inputs have controlled the dissolved oceanic budget of hafnium (Hf). However, whether the offset of seawater Nd-Hf isotope compositions from the terrestrial array can be fully generated by incongruent weathering of continental rocks (the zircon effect) is still not well constrained. In recent years, an increasing amount of combined U-Pb ages and Hf-isotopic compositions of riverine detrital zircons have been published. Here a new model of the Nd-Hf isotopic compositions of the weathered zircon-free part of the upper continental crust is presented, which is based on published Hf isotope compositions and formation ages of modern riverine detrital zircons combined with Nd isotopic compositions of rocks from the upper continental crust. Our model results indicate that the Nd-Hf isotopic composition of the weathered zircon-free part of the upper continental crust is not consistent with the seawater isotopic compositions. This suggests that the elevated seawater Hf isotope compositions for given Nd isotope compositions cannot be fully explained by incongruent zircon weathering of the continents, which is also supported by a recent study demonstrating incongruent weathering of other minerals than zircon.

Description: The use of phase equilibrium calculations to compute physical properties of rocks has become commonplace in geophysical modeling. Typically, the phase equilibrium calculations are used to construct two-dimensional tables of rock properties as a function of pressure and temperature. We document a computer program that can be used to assemble a three-dimensional table that accounts for compositional variations from such two-dimensional tables. A user-selected interpolation scheme is used to recover data from the tables as continuous function of its independent variables. We illustrate the utility of the program using two examples. The first explores the effect of water content of an average continental crust composition. As water content is the primary parameter controlling the amount of melt and/or low-density fluid present in the crust, this model provides a basis for assessing the impact of fluids on the geophysical properties of the crust. In the second example, we model the mantle composition as a mixture between harzburgitic and basaltic end-members. We show that for purposes of seismic velocity calculations, the continuum model is well approximated by interpolation from a table in which the compositional variable is sampled at 10% intervals. The tables for both examples are provided electronically as ready-to-use geophysical tools.

Description: Resistivity at the bit tools typically provide images of wellbore breakouts only a few minutes after the hole is drilled. In certain cases images are taken tens of minutes to days after drilling of the borehole. The sonic caliper can also image borehole geometry. We present four examples comparing imaging a few minutes after drilling to imaging from about 30 min to 3 days after drilling. In all cases the borehole breakouts widen with time. The tendency to widen with time is most pronounced within a few hundred meters below the seafloor (mbsf), but may occur at depths greater than 600 mbsf. In one example the widening may be due to reduced borehole fluid pressure that would enhance borehole failure. In the three other cases, significant decreases in fluid pressure during temporal evolution of breakouts are unlikely. The latter examples may be explained by time-dependent failure of porous sediments that are in an overconsolidated state due to drilling of the borehole. This time-dependent failure could be a consequence of dilational deformation, decrease of pore fluid pressure, and maintenance of sediment strength until migrating pore fluids weaken shear surfaces and allow spallation into the borehole. Breakout orientations, and thus estimates of stress orientations, remain consistent during widening in all four cases. In vertical boreholes, breakouts wider than those initially estimated by resistivity imaging would result in higher estimates of horizontal stress magnitudes. Because the vertical overburden stress is fixed, higher estimated horizontal stresses would favor strike-slip or thrust faulting over normal faulting.

Description: A new U-Pb zircon dating protocol for small (10–50 μm) zircons has been developed using an automated searching method to locate zircon grains in a polished rock mount. The scanning electron microscope-energy-dispersive X ray spectrum-based automated searching method can routinely find in situ zircon grains larger than 5 μm across. A selection of these grains was ablated using a 10 μm laser spot and analyzed in an inductively coupled plasma-quadrupole mass spectrometer (ICP-QMS). The technique has lower precision (∼6% uncertainty at 95% confidence on individual spot analyses) than typical laser ablation ICP-MS (∼2%), secondary ion mass spectrometry (

Description: Trace metals are known to be essential elements in marine ecosystems. Radiogenic isotopes of neodymium (Nd) have been used as tracers in many recent oceanic trace metal studies, although, among rare earth elements, cerium (Ce) isotopes might be an interesting complementary tracer for particle reactive and lithogenic metals such as manganese. This study determined the 138Ce/142Ce ratios in surface waters of the Pacific Ocean and its surrounding marginal seas: the Sulu Sea, the South China Sea, the East China Sea, and the South Australian Basin. The 138Ce/142Ce and 143Nd/144Nd data are discussed in terms of the sources of rare earth elements and elemental fractionation between Ce and Nd in the marine environment. In the Western North Pacific Central Water, East China Sea, and South China Sea, isotopic compositions of Ce ($\varepsilon$Ce = +0.7 to 1.4) are most affected by radiogenic Ce of continental origin. In contrast, less radiogenic isotopic compositions of Ce ($\varepsilon$Ce = −0.4 to +0.3) in the Pacific Equatorial Water were observed locally near volcanic islands such as New Guinea Island, suggesting the influence of mantle-derived Ce. Compared with Nd, the isotopic composition of Ce showed a heterogeneous distribution in a given surface water mass, reflecting the importance of local sources. Variations of isotopic compositions and concentrations of Ce in the western Equatorial Pacific and the East China Sea suggest that lithogenic Ce is supplied and scavenged by particle-dissolved interaction near the margins. Radiogenic Ce in the Western North Pacific Central Water, which is more continental-like than Nd isotopes, suggests direct input by atmospheric dust into the North Pacific Ocean. The isotopic distribution of Ce is sensitive to aeolian supply to the surface waters of the open ocean. This unique feature indicates that the 138Ce/142Ce ratio can be a useful chemical tracer for lithogenic trace elements such as iron and manganese, which have short oceanic residence time.

Description: Oceanic plates vary in temperature, topography, and sediment load as they enter subduction zones. These variations persist along the subduction interface causing perturbations in coupling and earthquake rupture processes. We explore the effects of variable subducting plate structure on microseismicity rupture characteristics along the Nicoya Peninsula, Costa Rica. The subducting Cocos Plate has low relief along the northern and central portion of the peninsula, with seamounts present at the southern tip of the peninsula. We compute apparent stresses for 94 ML 2.5–4.2 earthquakes along the plate interface using waveform coda and find along-strike variations that mimic bathymetric variability. Median stress values are higher (3.2 MPa) in the smooth northern region, with lower values in the central (2.1 MPa) and southern (0.7 MPa) segments. Higher apparent stresses along a zone of little suspected subduction topography imply increased coupling or higher friction along the interface. These results agree with geodetic and other seismic studies that suggest variable plate coupling along the Nicoya Peninsula.

Description: Identification of the mineral remains of magnetotactic bacteria (MTB), known as magnetofossils, is of particular interest because their occurrence can be used for environmental and climatic reconstructions. Single-domain magnetite particles, which are biomineralized in the cell body of MTB, have characteristic properties that can be used to detect their fossil remains. Acquisition of anhysteretic and isothermal remanent magnetization (ARM and IRM), first-order reversal curve (FORC) diagrams, and ferromagnetic resonance (FMR) spectra were used to detect the magnetic mineral inventory in Holocene lake sediments. A comparative analysis in terms of the discriminatory power of these methods is presented. The FORC diagrams contain two distinct features: a sharp horizontal ridge centered on the horizontal axis Bc and a feature with symmetric spread along the vertical Bb axis. The coercivity spectra derived from the central ridge coincides with that derived from ARM and IRM acquisition curves and is compatible with the presence of noninteracting linear chains of single-domain magnetite. The second feature on FORC diagrams is indicative of interacting particles in clusters. In the FMR spectra from bulk sediment, two populations are separated empirically based on the FORC information. An asymmetric signal is taken to describe the population, which contains single-domain particles in clusters. Empirical spectral separation of this contribution results in FMR spectra that are similar to those of intact MTB, which strongly suggests that a fraction of linear magnetosome chains is present. Combination of FMR and FORC results demonstrates the strong potential of these methods for identifying magnetofossils, based on alignment and interaction patterns of magnetic particles.

Description: The West Antarctic ice sheet (WAIS) has long been considered vulnerable to rapid retreat and today parts are rapidly losing ice. Projection of future change in WAIS is, however, hampered by our poor understanding of past changes, especially during interglacial periods that could be analogs for the future, but which undoubtedly provide an opportunity for testing predictive models. We consider how ice-loss would open seaways across WAIS; these would likely alter Southern Ocean circulation and climate, and would broadly define the de-glacial state, but they may also have left evidence of their existence in the coastal seas they once connected. We show the most likely routes for such seaways, and that a direct seaway between Weddell and Ross seas, which did not pass through the Amundsen Sea sector, is unlikely. Continued ice-loss at present rates would open seaways between Amundsen and Weddell seas (A-W), and Amundsen and Bellingshausen seas (A-B), in around one thousand years. This timescale indicates potential future vulnerability, but also suggests seaways may have opened in recent interglacial periods. We attempt to test this hypothesis using contemporary bryozoan species assemblages around Antarctica, concluding that anomalously high similarity in assemblages in the Weddell and Amundsen seas supports recent migration through A-W. Other authors have suggested opening of seaways last occurred during Marine Isotope Stage 7a (209 ka BP), but we conclude that opening could have occurred in MIS 5e (100 ka BP) when Antarctica was warmer than present and likely contributed to global sea levels higher than today.

Description: Detailed heat flow surveys were carried out in the central part of the Nankai Trough southeast of the Kii Peninsula (off Kumano) for investigation of the thermal structure of the subducting plate interface. At stations in the Kumano Trough (forearc basin) and its vicinity, long-term monitoring of temperature profiles in surface sediments was conducted because bottom water temperature variations (BTV) significantly disturb subbottom sediment temperatures. Heat flow values were successfully determined at seven stations by removing the influence of BTV from temperature records for 300 to 400 days. The surface heat flow data were combined with estimates from depths of methane hydrate bottom simulating reflectors to construct a heat flow profile across the subduction zone. Heat flow decreases from 90–110 mW/m2 on the floor of the Nankai Trough to 50–60 mW/m2 at around 30 km from the deformation front, while it is rather uniform, 40–60 mW/m2, in the Kumano Trough. The values measured on the Nankai Trough floor are concordant with the value estimated from the age of the subducting Philippine Sea plate, about 20 m.y., taking into account the effect of sedimentation. The obtained heat flow profile was used to constrain thermal models of the subduction zone. The subsurface thermal structure was calculated using a two-dimensional, steady state model, in which the frictional heating along the plate interface and the radioactive heat production are treated as unknown parameters. Comparison of the calculated surface heat flow in the Kumano Trough with the observed data indicates that the effective coefficient of friction is small, about 0.1 or less, and thus the shear stress on the plate interface is very low in this subduction zone.

Description: A record of soil CO2 flux data on the lower flanks of Mt. Etna, measured from 1989 to 2008, is discussed in the framework of multidisciplinary observations, including seismicity, ground deformations, and flank instability. A huge increase in soil CO2 fluxes appears to be related to the dynamics of magma ascent in the upper portions of the volcano (0–3 km) and the intrusion of dykes along the southern rift, as mainly observed before the 1991–93 eruption. Even after the 1991–93 eruption, the recharge of the shallow/central reservoir was accompanied by a long-term increase in soil CO2 degassing mainly in the southwestern area. The 2001 eruption marked dramatic changes in the areal distribution of seismicity, the deformation pattern, and the soil CO2 degassing. Indeed, while the soil CO2 degassing showed background values in the southwestern area, it progressively increased in the eastern sector and along the Pernicana fault. This has been related to the marked sliding of the eastern flank since the 2002–03 eruption and the associated seismicity. This study provides evidence that the extent of soil CO2 degassing on the lower flanks of Mt. Etna is controlled by (1) the volume of involved magma, (2) the intrusion of dykes in the upper parts of the volcano, and (3) fault movement and seismicity. This implies that different degassing structures must be monitored simultaneously when attempting to understand the behavior of the volcano as a whole.

Description: Analysis of sediment samples in the visible to mid infrared (IR) region requires small amounts of sample material and enables rapid and cost efficient geochemical analysis of mineral and organic sediment components. Here we use geochemical properties (total organic and inorganic carbon, biogenic silica, total nitrogen) from the ICDP deep drilling project PASADO to compare three different IR spectroscopy techniques: Diffuse Reflectance Fourier Transform IR Spectrometry (DRIFTS), Attenuated Total Reflectance Fourier Transform IR Spectroscopy (ATR-FTIRS) and Visible Near IR Spectroscopy (VNIRS). ATR-FTIRS and VNIRS are more rapid techniques compared to DRIFTS. Results show that calibration models developed using DRIFTS are most robust (correlation coefficient: R = 0.92 for TIC, R = 0.84 for BSi, R = 0.97 for TOC, R = 0.95 for TN). However, good statistical performance was also obtained by using ATR-FTIRS and VNIRS. When time and costs are limiting factors, these tools may be given preference for rapid biogeochemical screening.

Description: In geochronology, isotopic ages are determined from the ratio of parent and daughter nuclide concentrations in minerals. For dating of geological material using the K-Ar system, the simultaneous determination of 40Ar and 40K concentrations on the same aliquot is not possible. Therefore, a widely used variant, the 40Ar/39Ar technique, involves the production of 39Ar from 39K by neutron bombardment and the reliance on indirect natural calibrators of the neutron flux, referred to as “mineral standards.” Many mineral standards still in use rely on decades-old determinations of 40Ar concentrations; resulting uncertainties, both systematic and analytical, impede the determination of higher accuracy ages using the K-Ar decay system. We discuss the theoretical approach and technical design of a gas delivery system which emits metrologically traceable amounts of 40Ar and will allow for the sensitivity calibration of noble gas mass spectrometers. The design of this system is based on a rigorous assessment of the uncertainty budget and detailed tests of a prototype system. A number of obstacles and proposed resolutions are discussed along with the selection of components and their integration into a pipette system.

Description: Variations in gas emissions of open-vent volcanoes are investigated using a model of magma convection in narrow conduits. Laboratory experiments with both vertical and inclined conduits and dimensional analysis show that for Grashof numbers lower than 100 the volumetric rate of magma ascent is a simple function of equivalent conduit radius, density difference between the magmas, and viscosity of the degassed magma that descends back to the reservoir. The rate of magma ascent depends on the flux coefficient, estimated as 0.1 and 0.2 for vertical and inclined conduits, respectively. The equivalent radius parameter accounts for the dimensions of the conduit(s) regardless of its geometry, thus extending the treatment by previous models that used flow in pipes. The volume flow rate of convection increases with higher density difference and conduit size, but is also highly influenced by the large variations in viscosity of the degassed magma as volatile content and crystallinity change. The model presented here can be used to constrain the degassing and ascent rates of volatile-rich magma when combined with petrologic data on magmatic volatile content. Application of the model to Villarrica volcano (Chile) reveals that the background degassing levels observed (∼3 kg s−1 SO2) are associated with convective ascent of a relatively degassed magma (0.04 wt% S, ∼0.5 wt% H2O), while episodes of higher SO2 emissions (measurements up to 15 kg s−1) can be explained by the ascent of magma with higher volatile content (up to 0.09 wt% S, ∼1.5 wt% H2O).

Description: We analyze focal mechanisms of shallow-intermediate earthquakes in a NW-SE transect along the western Betic Cordillera and Alboran Sea, and deep earthquakes located in the central Betics to constrain the state of stress at the Gibraltar Arc slow convergence area. Shallow earthquakes ( 620 km) show very similar focal mechanisms, fitting the general slab behavior of resistance to further descent at the 660 km discontinuity. Seismicity features evidence the present-day stress distribution in a context of transition from subduction to continental collision.

Description: Hematite impact “bombs” are one of the most striking (and enigmatic) features of the large Araguainha impact structure in central Brazil. They have both porous or massive textures, elongated shapes from 5 to 50 cm in diameter, and botryoidal textures that suggest hydrothermal origin. Some authors have considered these objects as a possible analog of hematite nodules found in Mars, and consequently related to a hydrothermal system. Here we report rock magnetic measurements, X-ray diffraction and Mössbauer spectra for both massive and porous samples for a detailed description of the hematite. Room temperature magnetic measurements, including hysteresis loops, back-field and saturation magnetization acquisition, FORC, as well as X-ray diffraction and Mössbauer experiments are compatible with both massive and porous types being almost pure hematite. Room temperature FORCs after heating in a He atmosphere show two peaks; the original high-coercivity peak of hematite and a low-coercivity one (probably maghemite) at the Bc and Bb origin, thus indicating significant modification of the magnetic mineralogy of the material during thermal treatment in reducing conditions. However, conditioning in an oxidizing environment (heating in air) seemed to block generation of this low coercivity material in subsequent heating in a reducing (Ar) atmosphere. Therefore, we conclude that this material was not heated greatly in its generation. This would not be likely for impact-ejected bombs, so origin from post-impact hydrothermal activity seems likely.

Description: The Afar depression is an ideal locale to study the role of extension and magmatism as rifting progresses to seafloor spreading. Here we present receiver function results from new and legacy experiments. Crustal thickness ranges from ∼45 km beneath the highlands to ∼16 km beneath an incipient oceanic spreading center in northern Afar. The crust beneath Afar has a thickness of 20–26 km outside the currently active rift segments and thins northward. It is bounded by thick crust beneath the highlands of the western plateau (∼40 km) and southeastern plateau (∼35 km). The western plateau shows VP/VS ranging between 1.7–1.9, suggesting a mafic altered crust, likely associated with Cenozoic flood basalts, or current magmatism. The southeastern plateau shows VP/VS more typical of silicic continental crust (∼1.78). For crustal thicknesses 2.0) can only be explained by significant amounts of magmatic intrusions in the lower crust. This suggests that melt emplacement plays an important role in late stage rifting, and melt in the lower crust likely feeds magmatic activity. The crust between the location of the Miocene Red Sea rift axis and the current rift axis is thinner ( 2.0) than beneath the eastern part of Afar (〉26 km, VP/VS 〈 1.9). This suggests that the eastern region contains less partial melt, has undergone less stretching/extension and has preserved a more continental crustal signature than west of the current rift axis. The Red Sea rift axis appears to have migrated eastward through time to accommodate the migration of the Afar triple junction.

Description: Identical high precision U-Pb baddeleyite ages, together with paleomagnetic and geochemical data, on mafic dykes occurring over an area of 140,000 km2, define a Paleoproterozoic giant dyke swarm at ca. 2.367 Ga in the Dharwar craton, south India, referred to here as the Dharwar giant dyke swarm. All six U-Pb ages on these dykes are identical within error and suggest emplacement of this swarm within a geologically short time span of ∼5 Myr. A systematic southward progression in the trend of dykes from N48°E to N90°E, defines a fan angle of about 40° with convergence to a focal point about 300 km west of the present-day Dharwar craton boundary, resulting in a spectacular radiating dyke swarm extending across the entire eastern Dharwar craton. The large areal extent, radiating dyke pattern and short duration imply a mantle plume origin for the Dharwar giant dyke swarm. Despite their large areal distribution, all dykes in this swarm are geochemically coherent and have similar primitive mantle-normalized trace element patterns and rare earth element characteristics. Although the NE part of the swarm is magnetically overprinted, a remanence survives that has the same direction as primary magnetizations from dykes in the southern part of the swarm.

Description: Surface displacements solutions of elastic deformation around an inflating magma chamber generally assume that the associated internal overpressure is limited by the bedrock tensile strength. When considering stress equilibrium in the bedrock adjacent to a spherical or infinitely long cylinder, the gravity body force actually resists tensile failure, thus leading to a much larger pressure threshold. And when considering a Coulomb failure criterion, analytical and numerical models predict that shear failure develops instead of tensile failure. Here, three numerical codes are used to compare elasto-plastic solutions of surface displacements and patterns of failure in plane-strain. Shear failure propagates independently from the surface downward, then from the chamber walls upwards, and finally the two plasticized domains connect. Another test with internal underpressure (simulating source deflation) fits standard solutions from tunneling engineering. The effect of pore fluid pressures is also explored. In case of lithostatic fluid pore pressure in the bedrock, the gravity effect cancels out, and tensile failure is enabled for an overpressure close to the tensile strength. Coupled hydromechanical models in undrained conditions indicate that the initial bedrock porosity modifies the evolution of fluid pressure, volumetric strain and effective normal stress, and consequently also the pressure threshold for the onset of failure. We show that a bedrock of low porosity is more prone to fail than a bedrock of high porosity. In summary, our elasto-plastic and hydromechanical models illustrate the contexts for either tensile or shear failure around magmatic bodies, at the same time complementing and delimiting predictions deduced from elasticity.

Description: In this paper we look at changes of the statistics of the stationary explosive process at a basaltic volcano as a proxy for departures from thermodynamic equilibrium conditions in the shallow plumbing system. Specifically, we investigate the explosion process of Stromboli volcano that occurred during 2002–2003, 2006–2007 and 2010–2011. The first two periods were characterized by eruptions with significant lava effusion and strong paroxysmal events, while the last one shows persistent explosive activity accompanied by minor episodes of lava flow. We use three-component seismic data recorded by broadband stations operating on the volcano and, for 2007 and 2010–2011 cases, strainmeter data from a Sacks-Evertson borehole dilatometer. For each time interval we study the explosive process by looking at the inter-occurrence times and at the amplitude distribution. Moreover, we analyze its waveforms, spectral content and polarization properties. In all three cases we find sharp increases of the explosion rate associated with swarms. Swarms are characterized by quasi-monochromatic seismic events with frequency peak close to about 3 Hz, higher amplitude than the usual explosions and variability coefficient of the inter-occurrence times close to 0.5. In correspondence to the swarms, we also observe negative variations in the strain signals, which indicate a depressurization in the shallow plumbing system. This depressurization emerges clearly from the data collected during 2010–2011, whereas it is less sharp for the 2007 episode, and has been estimated in about 105 Pa. From the polarization analysis we infer that this depressurization affects the upper 0.3–0.8 km of the plumbing system.

Description: Convergent plate boundaries accommodate intraplate displacement within a ∼100–1000 m thick shear zone. Marine geophysicists typically define this zone, the subduction channel (SC), as the sedimentary layer between the downgoing oceanic crust and the base of the upper plate. Geologists and modelers, instead, perceive the SC as a specific type of shear zone. The original theory of SCs was developed when the net accretion of marine sediments to the forearc was thought to typify a convergent margin. While erosive margins were briefly mentioned, their mechanics were not discussed in any detail. We now realize that subduction erosion is taking place at roughly half of the modern subduction margins. Here we review and revise the theory of erosive SCs (1) to unify this concept across disciplines, focusing on the meaning of the channel's boundaries; (2) to redefine the portions of the forearc included in the SC concept; and (3) to better idealize this dynamic system where material supply to the channel, fluid content, and the heterogeneity of deformation all influence the SC's upper and lower boundaries. Migration of the channel boundaries controls the downdip variation of tectonic mechanisms that shape the margin. Within the shallow,

Description: The abundances of 60 elements in 616 Ocean Floor Basaltic (OFB) glasses from the Abyssal Volcanic Glass Data File (AVGDF) of the Smithsonian Institution have been determined by laser-ablation (LA)-ICP-MS and electron microprobe analysis (EMPA). The elements analyzed include all 28 of the refractory lithophile elements, which provide the framework for establishing the geochemical behavior and source abundances of volatile, chalcophile and siderophile elements. In addition to the traditionally analyzed elements (rare earth elements (REE), high field strength elements (HFSE), large ion lithophile elements (LILE) and first row transition elements (FRTE)), we report analyses for lesser-analyzed elements (Li, Be, Ga, Ge, As, Se, Mo, Ag, Cd, In, Sn, Sb, W, Tl and Bi). The precision of the method for most elements is between 2 and 4%, one standard deviation, although ratios of elements determined simultaneously are more precise (e.g., REE, Zr/Hf). Subsets of 329 glasses were analyzed by electron microprobe for S and 154 glasses for Cl. The results define a representative trace element geochemistry of OFB, against which local variations resulting from differences in basalt petrogenesis in a range of tectonic settings or different styles of magmatic differentiation may be compared.

Description: Several factors influence the topography of the overriding plate in a subduction zone: mantle wedge flow caused by motion of the descending slab, lithosphere elasticity, and interaction between the solid lithosphere and the viscous asthenosphere. This paper presents the results of new subduction modeling that incorporates both elastic/elastoplastic deformation in the overriding plate and viscous deformation in the asthenosphere. The effects of fluid-solid coupling and interaction were considered using an arbitrary Lagrangian-Eulerian technique that solves for the continuous deformation of fluid flow in contact with the rigid plate. Modeling results show that thicker lithosphere, which is harder to deform, contributes to the development of broadly shallower back-arc basins. As lithospheric strength decreases, its ability to withstand outside force (fluid load) resulting from flow in the mantle wedge also declines, leading to the development of deep and narrow basins. Large-scale deformation of the lithosphere modifies both the mantle wedge geometry and the wedge flow pattern, which in turn increases the fluid load and thereby enhances topographic variation. Generally, the effect of interaction is enhanced with lithospheric deformation. The elastoplastic model also leads to a weaker lithosphere and a correspondingly larger deflection in the lithosphere. Similarly to previous studies, back-arc basins are deeper when mantle viscosity is higher. For a 50 km thick elastic lithosphere, the lithosphere-asthenosphere interaction is negligible when the asthenospheric viscosity is less than 4 × 1019 Pa s. Increasing the rate of subduction or decreasing the dip angle will also lead to the development of deeper basins.

Description: We make use of observations on orogenic strain accumulation and deformation partitioning in the Central Andes to explore the backarc strength evolution at the lithospheric scale. In plan view, the Altiplano-Puna plateaux experienced rapid initial increase of surface area undergoing active deformation during the Cenozoic. Beyond the maximum lateral extent reached around 10–15 Ma (40–50% of entire proto-Andes undergoing deformation) at 10–20% total strain, rapid localization initiated at the eastern flank of the Altiplano (Inter- and Subandean thrust belt) but not at the Puna latitude. Localization was associated with a significant increase in bulk shortening rate. Average fault slip rates equally increased by an order of magnitude following a protracted period of stable average rates. Estimates of strength evolution based on force balance calculations and critical wedge analysis suggest significant backarc weakening driving this change after the Middle Miocene. Strain accumulation led to localization and weakening with development of a detachment propagating through crust and upper mantle. We find that lithosphere-scale failure resulting from strain weakening beyond a critical strain threshold (c. 20%) and fault coalescence with formation of a weak detachment in shales (effective coefficient of friction 〈 0.1) plays a key role in the evolution of the Andes. Strain-related lithosphere weakening appears to dominate over the impact of external forcing mechanisms, such as variations of plate convergence, mantle-assisted processes, or erosion. Comparison of these orogen-scale observations with experimental rock rheology indicates substantial similarity of deformation behavior with similar weakening thresholds across a wide range of scales.

Description: Indonesia is arguably one of the tectonically most complex regions on Earth today due to its location at the junction of several major tectonic plates and its long history of collision and accretion. It is thus an ideal location to study the interaction between subducting plates and mantle convection. Seismic anisotropy can serve as a diagnostic tool for identifying various subsurface deformational processes, such as mantle flow, for example. Here, we present novel shear wave splitting results across the Indonesian region. Using three different shear phases (local S, SKS, and downgoing S) to improve spatial resolution of anisotropic fabrics allows us to distinguish several deformational features. For example, the block rotation history of Borneo is reflected in coast-parallel fast directions, which we attribute to fossil anisotropy. Furthermore, we are able to unravel the mantle flow pattern in the Sulawesi and Banda region: We detect toroidal flow around the Celebes Sea slab, oblique corner flow in the Banda wedge, and sub-slab mantle flow around the arcuate Banda slab. We present evidence for deep, sub-520 km anisotropy at the Java subduction zone. In the Sumatran backarc, we measure trench-perpendicular fast orientations, which we assume to be due to mantle flow beneath the overriding Eurasian plate. These observations will allow to test ideas of, for example, slab–mantle coupling in subduction regions.

Description: We present the interpretation of newly acquired high-quality industry-standard deep seismic reflection and swath bathymetry data to provide insight into the structural style and evolution of the Mentawai Fault Zone (MFZ). The MFZ lies along the boundary between the accretionary wedge and the proposed continental backstop. This zone exhibits arcuate ridges on the seafloor, convex toward the east. Beneath these ridges the structures developed as landward-vergent imbricated backthrusts in the inner part of the accretionary wedge and higher-angle backthrusts that deformed the forearc basin sediments. In the forearc high, anticlines were developed due to the seaward-vergent forearc high thrusts originating in the accretionary wedge. The imbricated backthrusts may have initiated during the Early-Middle Miocene contemporaneously with the slide and back-rotation of forearc high thrusts. In the Late Miocene, the higher-angle backthrusts were initiated. Continuous contraction induced the frontal higher-angle backthrusts and formed a fold-thrust belt toward the east during the Pliocene. The folds and thrusts were disturbed by diapirs and mud volcanoes. Backthrusting and fold-thrust belts developed in the MFZ may explain the compressional features observed at the boundary between the accretionary wedge and continental backstop along the southern Sumatra margin. The backthrusts along the MFZ are waning in activity and hence the risk of a large earthquake and associated tsunami at the present time should be small.

Description: This article introduces the software Melange, a 3D lattice-particle hybrid model. The software was specifically designed in order to simulate ductile visco-elasto-plastic deformation and can be used to study tectonic processes in the lithosphere from the micro to the macro scale. Melange is under an open source license. The code takes both relevant yield mechanisms for the deformation of lithospheric material into account: dynamic brittle failure and ductile creep, where ductile creep is modeled as viscoelasticity. The software considers effects of the local geology, of the inherent disorder of geomaterials, of rheological layering of the lithosphere and applies repulsion when the material fractures. Driving forces are the externally applied strain and the gravitational load. Melange applies an elastically isotropic regular 3D lattice with HCP geometry and next-nearest neighbor interactions. Young's modulus, viscosity, material density and system size can be freely chosen. Poisson's ratio is restricted to values

Description: [1]  Plagioclase and spinel lherzolite thermometry and barometry is applied to an extensive geochemical dataset of young (〈10.5 Ma) primitive basaltic lavas from across Oregon's High Lava Plains, California's Modoc Plateau, and the central - southern Cascades volcanic arc to calculate the depths and temperatures of mantle melting. This study focuses on basalts with low pre-eruptive H 2 O-contents that are little fractionated near-primary melts of mantle peridotite (i.e., basalts thought to be products of anhydrous decompression mantle melting). Calculated minimum depths of nominally anhydrous melt extraction are 40–58 km below Oregon's High Lava Plains, 41–51 km below the Modoc Plateau, and 37–60 km below the central and southern Cascades arc. The calculated depths are very close to Moho depths as determined from a number of regional geophysical studies and suggest that the geophysical Moho and lithosphere-asthenosphere boundary in this region are located in very close proximity to one another (within 5–10 km). The basalts originated at 1185–1383 °C and point to a generally warm mantle beneath this area, but not one hot enough to obviously require a plume contribution. Our results, combined with a range of other geologic, geophysical, and geochemical constraints, are consistent with a regional model whereby anhydrous mantle melting over the last 10.5 Ma in a modern convergent margin and back-arc was driven by subduction-induced corner flow in the mantle wedge, and to a lesser extent toroidal flow around the southern edge of the subducting Juan de Fuca and Gorda plates, and crustal extension-related upwelling of the shallow mantle.

Description: [1]  We determine the 3-D melt geometry of partially molten samples of dunite containing 1.6 and 3.6 volume % of basaltic melt that were held in a piston cylinder apparatus at upper mantle conditions for 430 hours. Our approach involves serial sectioning and high-resolution field emission SEM imaging. Resolution is such that melt pockets approaching ~30 nm in size were resolved while covering an area of ~300 by 230 micron. The principal result of this study is to show that thin layers (typically 100 nm or less in thickness) between adjacent grains observed in 2-D images persist with depth, and are therefore wetted two-grain boundaries. Melt geometries most closely resembling triple junction tubules of the isotropic equilibrium model occur at all three grain edges, but are small compared to larger pockets. The wetted grain boundaries at a dihedral angle 〉0° for this system are inferred to be due to slow expulsion of melt from dynamically reorganizing grain boundaries during steady state grain growth. The attenuation peak observed in forced torsional oscillation experiments on similar samples ( [18]) is likely related to the wetted grain boundaries. Grain growth, driven by surface energy reduction, occurs also at the larger grain sizes expected for the mantle. This suggests the presence of wetted grain boundaries and significant velocity reduction and attenuation in partially molten upper mantle, as observed for example in back-arc basins.

Description: [1]  Thellier-type experiments are a method used to estimate the intensity of the ancient geomagnetic field from samples carrying thermoremanent magnetization. The analysis of Thellier-type experimental data is conventionally done by manually interpreting data from each specimen individually. The main limitations of this approach are: (1) manual interpretation is highly subjective and can be biased by misleading concepts (2) the procedure is time consuming, and (3) unless the measurement data are published, the final results cannot be reproduced by readers. These issues compound when trying to combine together paleointensity data from a collection of studies. Here we address these problems by introducing the Thellier GUI : a comprehensive tool for interpreting Thellier type experimental data. The tool presents a graphical user interface, which allows manual interpretation of the data, but also includes two new interpretation tools: 1) Thellier Auto Interpreter : an automatic interpretation procedure based on a given set of experimental requirements, and 2) Consistency Test : a self-test for the consistency of the results assuming groups of samples that should have the same paleointensity values. We apply the new tools to data from two case studies. These demonstrate that interpretation of non-ideal Arai plots is non-unique and different selection criteria can lead to significantly different conclusions. Hence, we recommend adopting the automatic interpretation approach, as it allows a more objective interpretation, which can be easily repeated or revised by others. When the analysis is combined with a Consistency Test the credibility of the interpretations is enhanced. We also make the case that published paleointensity studies should include the measurement data (as supplementary files or as a contributions to the MagIC database) so that results based on a particular data set can be reproduced and assessed by others.

Description: [1]  The mean composition of mid-ocean ridge basalts (MORB) is determined using a global dataset of major elements, trace elements, and isotopes compiled from new and previously published data. A global catalog of 771 ridge segments, including their mean depth, length, and spreading rate, enables calculation of average compositions for each segment. Segment averages permit weighting by segment length and spreading rate and reduce the bias introduced by uneven sampling. A bootstrapping statistical technique provides rigorous error estimates. Based on the characteristics of the data, we suggest a revised nomenclature for MORB. “ all morb ” is the total composition of the crust apart from back-arc basins, N-MORB the most likely basalt composition encountered along the ridge 〉500 km from hot spots, and D-MORB the depleted end-member. all morb and N-MORB are substantially more enriched than early estimates of normal ridge basalts. The mean composition of back-arc spreading centers requires higher extents of melting and greater concentrations of fluid-mobile elements, reflecting the influence of water on back-arc petrogenesis. The average data permit a re-evaluation of several problems of global geochemistry. The K/U ratio reported here (12,340±840) is in accord with previous estimates, much lower than the estimate of Arevalo et al. (2009). The low Sm/Nd and 143 Nd/ 144 Nd ratio of all morb and N-MORB provide constraints on the hypothesis that Earth has a non-chondritic primitive mantle. Either Earth is chondritic in Sm/Nd and the hypothesis is incorrect, or MORB preferentially sample an enriched reservoir, requiring a large depleted reservoir in the deep mantle.

Description: [1]  Long-term changes of geomagnetic field intensity, including possible dependence on lengths of polarity intervals, provide fundamentally important information for understanding the geodynamo. A positive correlation between paleointensity and polarity interval length was previously suggested from an Oligocene (ca 23–34 Ma) relative paleointensity record at Deep Sea Drilling Program Site 522 in the Atlantic Ocean, which is the only continuous paleointensity data set published so far for this age interval. We have conducted a paleomagnetic study of Eocene to Oligocene sediments at three sites in the eastern equatorial Pacific Ocean. Our objectives include revisiting the issue of the paleointensity-polarity length correlation. Magnetic properties of the sediments meet the frequently used criteria for reliable relative paleointensity estimation. Although short-wavelength normalized remanence intensity fluctuations associated with polarity boundaries and possible geomagnetic excursions agree among the three sites, long-term changes are inconsistent. Apparent positive correlation between normalized intensity and polarity length was observed, but the normalized intensity has an obvious anti-correlation with the ratio of anhysteretic remanent magnetization (ARM) to isothermal remanent magnetization (IRM), which is mainly controlled by the relative abundance of biogenic and terrigenous magnetic minerals. Furthermore, the normalized intensity correlates with sedimentation rate. These facts indicate a lithological contamination on the normalized intensity records. The dependence on ARM/IRM and sedimentation rate is also evident at Site 522. It is inferred that variations in sedimentation rate and the relative abundance of biogenic magnetite on depositional remanent magnetization acquisition efficiency may not be well compensated by the normalization. It is therefore premature to conclude that stronger geomagnetic fields were recorded during longer polarity intervals from currently available normalized intensity records.

Description: Concentrations of cosmogenic iodine, 129 I, in the pore fluid of marine sediments often indicate that the pore fluid is much older than the host sediment, even when vertical flow due to sediment compaction is taken into account. Old pore fluid has been used in previous studies to argue for pervasive upward fluid flow and a deep methane source for hydrate deposits. Alternatively, old pore fluid age may reflect more complex flow patterns. We use a two-dimensional numerical transport model to account for the effects of topography and fractures on pore fluid pathlines when sediment permeability is anisotropic. We find that fluid focusing can cause significant lateral migration as well as regions where downward flow reverses direction and returns toward the seafloor. Longer pathlines can produce pore fluid ages much older than that expected with a one-dimensional compaction model. For steady-state models with geometry representative of Blake Ridge (USA), a well-studied hydrate province, we find pore fluid ages beneath regions of topography and within fractured zones that are up to 70 Ma old. Our results suggest that the measurements of 129 I/ 127 I reflect a mixture of new and old pore fluid. However, old pore fluid need not originate at great depths. Methane within pore fluids can travel laterally several kilometers, implying an extensive source region around the deposit. This type of focusing should aid hydrate formation beneath topographic highs. © 2013 American Geophysical Union. All rights reserved