SEARCH · USGS Science

Results for “Earth and Planetary Science Letters”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 325 records · Page 18Linked to original sources

USGS2015

Widespread occurrence of (per)chlorate in the Solar System

Perchlorate (ClO− 4 ) and chlorate (ClO− 3 ) are ubiquitous on Earth and ClO− 4 has also been found on Mars. These species can play important roles in geochemical processes such as oxidation of organic matter and as biological electron acceptors, and are also indicators of important photochemical reactions involving oxyanions; on Mars they could be relevant for human habitability both in terms of in situ resource utilization and potential human health effects. For the first time, we extracted, detected and quantified ClO− 4 and ClO− 3 in extraterrestrial, non-planetary samples: regolith and rock samples from the Moon, and two chondrite meteorites (Murchison and Fayetteville). Lunar samples were collected by astronauts during the Apollo program, and meteorite samples were recovered immediately after their fall. This fact, together with the heterogeneous distribution of ClO− 4 and ClO− 3 within some of the samples, and their relative abundance with respect to other soluble species (e.g., NO− 3 ) are consistent with an extraterrestrial origin of the oxychlorine species. Our results, combined with the previously reported widespread occurrence on Earth and Mars, indicate that ClO− 4 and ClO− 3 could be present throughout the Solar System.

Earth and Planetary Science Letters↗

USGS2016

Radiometric dating of marine-influenced coal using Re–Os geochronology

Coal deposits are integral to understanding the structural evolution and thermal history of sedimentary basins and correlating contemporeous estuarine and fluvial delatic strata with marine sections. While marine shales may readily lend themselves to Re–Os dating due to the dominance of hydrogenous Re and Os, the lack of a chronometer for near-shore sedimentary environments hampers basinwide correlations in absolute time. Here, we employ the Re–Os geochronometer, along with total organic carbon (TOC) and Rock–Eval data, to determine the timing and conditions of a marine incursion at the top of the Matewan coal bed, Kanawha Formation, Pottsville Group, West Virginia, USA. The observed range for hydrogen index (HI: 267–290 mg hydrocarbon/gram total organic carbon) for these coal samples suggests dominance of aliphatic hydrocarbons with low carbon (<C 19 ) chain length. Average Re ( 107.6±16.4 ng/g ) and Os ( 0.52±0.09 ng/g ) concentrations of the marine-influenced Matewan coal are higher by few orders of magnitude than published data for terrestrial coal. A Re–Os isochron for the Matewan coal provides an age of 325±14 Ma (Model 3; MSWD = 12; n=19 ; 2 σ ). This is the first Re–Os age derived from coal samples; the age overlaps a new composite Re–Os age of 317±2 Ma for the immediately overlying Betsie Shale Member. External precision for replicate Os analyses carried out for several Matewan coal samples shows a positive correlation with their HI. The HI, which is low in terrestrial organic matter, reflects the degree of marine influence. Thus, samples with the most profound marine influence also have the best analytical reproducibility. Equilibration of Os isotopes with seawater under marine conditions overwhelms variability inherited from terrestrial plant debris, decreasing scatter on the isochron. The 187 Re/ 188 Os ratios of the Matewan coal (&sim;3300–5135) are higher than most of those previously published for Phanerozoic black shale (mostly <2000). Mass balance calculations based on Re/TOC and Os/TOC ratios for the Matewan coal indicate that both Re and Os are primarily marine in origin, and their high 187 Re/ 188 Os ratios confirm efficient removal of both elements from a sulfidic water column into the coal. We show that Re–Os geochronology of marine-influenced coal can be a viable tool for constraining depositional ages.

West Virginia↗

USGS2015

Coral 13C/12C records of vertical seafloor displacement during megathrust earthquakes west of Sumatra

The recent surge of megathrust earthquakes and tsunami disasters has highlighted the need for a comprehensive understanding of earthquake cycles along convergent plate boundaries. Space geodesy has been used to document recent crustal deformation patterns with unprecedented precision, however the production of long paleogeodetic records of vertical seafloor motion is still a major challenge. Here we show that carbon isotope ratios ( ) in the skeletons of massive Porites corals from west Sumatra record abrupt changes in light exposure resulting from coseismic seafloor displacements. Validation of the method is based on the coral response to uplift (and subsidence) produced by the March 2005 M w 8.6 Nias–Simeulue earthquake, and uplift further south around Sipora Island during a M&sim;8.4 megathrust earthquake in February 1797. At Nias, the average step-change in coral was 0.6±0.1&permil;/m for coseismic displacements of +1.8 m and −0.4 m in 2005. At Sipora, a distinct change in Porites microatoll growth morphology marks coseismic uplift of 0.7 m in 1797. In this shallow water setting, with a steep light attenuation gradient, the step-change in microatoll is 2.3&permil;/m , nearly four times greater than for the Nias Porites . Considering the natural variability in coral skeletal , we show that the lower detection limit of the method is around 0.2 m of vertical seafloor motion. Analysis of vertical displacement for well-documented earthquakes suggests this sensitivity equates to shallow events exceeding M w &sim;7.2 in central megathrust and back-arc thrust fault settings. Our findings indicate that the coral paleogeodesy technique could be applied to convergent tectonic margins throughout the tropical western Pacific and eastern Indian oceans, which host prolific coral reefs, and some of the world's greatest earthquake catastrophes. While our focus here is the link between coral , light exposure and coseismic crustal deformation, the same principles could be used to characterize interseismic strain during earthquake cycles over the last several millennia.

Earth and Planetary Science Letters↗

USGS2016

Lithospheric rheology constrained from twenty-five years of postseismic deformation following the 1989 Mw 6.9 Loma Prieta earthquake

The October 17, 1989 M w 6.9 Loma Prieta earthquake provides the first opportunity of probing the crustal and upper mantle rheology in the San Francisco Bay Area since the 1906 M w 7.9 San Francisco earthquake. Here we use geodetic observations including GPS and InSAR to characterize the Loma Prieta earthquake postseismic displacements from 1989 to 2013. Pre-earthquake deformation rates are constrained by nearly 20 yr of USGS trilateration measurements and removed from the postseismic measurements prior to the analysis. We observe GPS horizontal displacements at mean rates of 1–4 mm/yr toward Loma Prieta Mountain until 2000, and ∼2 mm/yr surface subsidence of the northern Santa Cruz Mountains between 1992 and 2002 shown by InSAR, which is not associated with the seasonal and longer-term hydrological deformation in the adjoining Santa Clara Valley. Previous work indicates afterslip dominated in the early (1989–1994) postseismic period, so we focus on modeling the postseismic viscoelastic relaxation constrained by the geodetic observations after 1994. The best fitting model shows an elastic 19-km-thick upper crust above an 11-km-thick viscoelastic lower crust with viscosity of ∼6 × 10 18 Pas, underlain by a viscous upper mantle with viscosity between 3 × 1018 and 2 × 10 19 Pas. The millimeter-scale postseismic deformation does not resolve the viscosity in the different layers very well, and the lower-crustal relaxation may be localized in a narrow shear zone. However, the inferred lithospheric rheology is consistent with previous estimates based on post-1906 San Francisco earthquake measurements along the San Andreas fault system. The viscoelastic relaxation may also contribute to the enduring increase of aseismic slip and repeating earthquake activity on the San Andreas fault near San Juan Bautista, which continued for at least a decade after the Loma Prieta event.

California↗

USGS2016

Corrigendum to “Widespread occurrence of (per)chlorate in the Solar System” [Earth Planet. Sci. Lett. 430 (2015) 470–476]

The authors regret that two sets of data (Atacama (Rao et al., 2010) and Mars Meteorite Range (Kounaves et al., 2014)) in Fig. 2 of our article were plotted in the wrong units. The correction does not change the relationship between ClO 3 − "> ClO3− and ClO 4 − "> ClO4− ; it only shifts the magnitude of the concentrations. The conclusions of the article are not affected. The corrected Fig. 2 appears below.

Earth and Planetary Science Letters↗

USGS2016

A submarine landslide source for the devastating 1964 Chenega tsunami, southern Alaska

During the 1964 Great Alaska earthquake (M w 9.2), several fjords, straits, and bays throughout southern Alaska experienced significant tsunami runup of localized, but unexplained origin. Dangerous Passage is a glacimarine fjord in western Prince William Sound, which experienced a tsunami that devastated the village of Chenega where 23 of 75 inhabitants were lost – the highest relative loss of any community during the earthquake. Previous studies suggested the source of the devastating tsunami was either from a local submarine landslide of unknown origin or from coseismic tectonic displacement. Here we present new observations from high-resolution multibeam bathymetry and seismic reflection surveys conducted in the waters adjacent to the village of Chenega. The seabed morphology and substrate architecture reveal a large submarine landslide complex in water depths of 120–360 m. Analysis of bathymetric change between 1957 and 2014 indicates the upper 20–50 m ( &sim; 0.7 km 3 ) of glacimarine sediment was destabilized and evacuated from the steep face of a submerged moraine and an adjacent &sim; 21 km 2 perched sedimentary basin. Once mobilized, landslide debris poured over the steep, 130 m-high face of a deeper moraine and then blanketed the terminal basin ( &sim; 465 m water depth) in 11 ± 5 m of sediment. These results, combined with inverse tsunami travel-time modeling, suggest that earthquake- triggered submarine landslides generated the tsunami that struck the village of Chenega roughly 4 min after shaking began. Unlike other tsunamigenic landslides observed in and around Prince William Sound in 1964, the failures in Dangerous Passage are not linked to an active submarine delta. The requisite environmental conditions needed to generate large submarine landslides in glacimarine fjords around the world may be more common than previously thought.

Alaska↗

USGS2016

Sensitivity of Pliocene Arctic climate to orbital forcing, atmospheric CO 2 and sea ice albedo parameterisation

General circulation model (GCM) simulations of the mid-Pliocene Warm Period (mPWP, 3.264 to 3.025 Myr ago) do not reproduce the magnitude of Northern Hemisphere high latitude surface air and sea surface temperature (SAT and SST) warming that proxy data indicate. There is also large uncertainty regarding the state of sea ice cover in the mPWP. Evidence for both perennial and seasonal mPWP Arctic sea ice is found through analyses of marine sediments, whilst in a multi-model ensemble of mPWP climate simulations, half of the ensemble simulated ice-free summer Arctic conditions. Given the strong influence that sea ice exerts on high latitude temperatures, an understanding of the nature of mPWP Arctic sea ice would be highly beneficial. Using the HadCM3 GCM, this paper explores the impact of various combinations of potential mPWP orbital forcing, atmospheric CO 2 concentrations and minimum sea ice albedo on sea ice extent and high latitude warming. The focus is on the Northern Hemisphere, due to availability of proxy data, and the large data–model discrepancies in this region. Changes in orbital forcings are demonstrated to be sufficient to alter the Arctic sea ice simulated by HadCM3 from perennial to seasonal. However, this occurs only when atmospheric CO 2 concentrations exceed 300 ppm. Reduction of the minimum sea ice albedo from 0.5 to 0.2 is also sufficient to simulate seasonal sea ice, with any of the combinations of atmospheric CO 2 and orbital forcing. Compared to a mPWP control simulation, monthly mean increases north of 60°N of up to 4.2 °C (SST) and 9.8 °C (SAT) are simulated. With varying CO 2 , orbit and sea ice albedo values we are able to reproduce proxy temperature records that lean towards modest levels of high latitude warming, but other proxy data showing greater warming remain beyond the reach of our model. This highlights the importance of additional proxy records at high latitudes and ongoing efforts to compare proxy signals between sites.

Earth and Planetary Science Letters↗

USGS2016

Short-period volcanic gas precursors to phreatic eruptions: Insights from Poás Volcano, Costa Rica

Volcanic eruptions involving interaction with water are amongst the most violent and unpredictable geologic phenomena on Earth. Phreatic eruptions are exceptionally difficult to forecast by traditional geophysical techniques. Here we report on short-term precursory variations in gas emissions related to phreatic blasts at Poás volcano, Costa Rica, as measured with an in situ multiple gas analyzer that was deployed at the edge of the erupting lake. Gas emitted from this hyper-acid crater lake approaches magmatic values of SO 2 /CO 2 1–6 days prior to eruption. The SO 2 flux derived from magmatic degassing through the lake is measureable by differential optical absorption spectrometry (sporadic campaign measurements), which allows us to constrain lake gas output and input for the major gas species during eruptive and non-eruptive periods. We can further calculate power supply to the hydrothermal system using volatile mass balance and thermodynamics, which indicates that the magmatic heat flux into the shallow hydrothermal system increases from &sim;27 MW during quiescence to &sim;59 MW during periods of phreatic events. These transient pulses of gas and heat from the deeper magmatic system generate both phreatic eruptions and the observed short-term changes in gas composition, because at high gas flux scrubbing of sulfur by the hydrothermal system is both kinetically and thermodynamically inhibited whereas CO 2 gas is always essentially inert in hyperacid conditions. Thus, the SO 2 /CO 2 of lake emissions approaches magmatic values as gas and power supply to the sub-limnic hydrothermal system increase, vaporizing fluids and priming the hydrothermal system for eruption. Our results suggest that high-frequency real-time gas monitoring could provide useful short-term eruptive precursors at volcanoes prone to phreatic explosions.

Poás Volcano↗

USGS2016

Slab melting and magma formation beneath the southern Cascade arc

The processes that drive magma formation beneath the Cascade arc and other warm-slab subduction zones have been debated because young oceanic crust is predicted to largely dehydrate beneath the forearc during subduction. In addition, geochemical variability along strike in the Cascades has led to contrasting interpretations about the role of volatiles in magma generation. Here, we focus on the Lassen segment of the Cascade arc, where previous work has demonstrated across-arc geochemical variations related to subduction enrichment, and H-isotope data suggest that H 2 O in basaltic magmas is derived from the final breakdown of chlorite in the mantle portion of the slab. We use naturally glassy, olivine-hosted melt inclusions (MI) from the tephra deposits of eight primitive ( MgO>7 wt% ) basaltic cinder cones to quantify the pre-eruptive volatile contents of mantle-derived melts in this region. The melt inclusions have B concentrations and isotope ratios that are similar to mid-ocean ridge basalt (MORB), suggesting extensive dehydration of the downgoing plate prior to reaching sub-arc depths and little input of slab-derived B into the mantle wedge. However, correlations of volatile and trace element ratios (H 2 O/Ce, Cl/Nb, Sr/Nd) in the melt inclusions demonstrate that geochemical variability is the result of variable addition of a hydrous subduction component to the mantle wedge. Furthermore, correlations between subduction component tracers and radiogenic isotope ratios show that the subduction component has less radiogenic Sr and Pb than the Lassen sub-arc mantle, which can be explained by melting of subducted Gorda MORB beneath the arc. Agreement between pMELTS melting models and melt inclusion volatile, major, and trace element data suggests that hydrous slab melt addition to the mantle wedge can produce the range in primitive compositions erupted in the Lassen region. Our results provide further evidence that chlorite-derived fluids from the mantle portion of the slab ( ∼7–9 km below the slab top) cause flux melting of the subducted oceanic crust, producing hydrous slab melts that migrate into the overlying mantle, where they react with peridotite to induce further melting.

Cascades↗

USGS2016

Bayesian estimation of magma supply, storage, and eruption rates using a multiphysical volcano model: Kīlauea Volcano, 2000–2012

Estimating rates of magma supply to the world's volcanoes remains one of the most fundamental aims of volcanology. Yet, supply rates can be difficult to estimate even at well-monitored volcanoes, in part because observations are noisy and are usually considered independently rather than as part of a holistic system. In this work we demonstrate a technique for probabilistically estimating time-variable rates of magma supply to a volcano through probabilistic constraint on storage and eruption rates. This approach utilizes Bayesian joint inversion of diverse datasets using predictions from a multiphysical volcano model, and independent prior information derived from previous geophysical, geochemical, and geological studies. The solution to the inverse problem takes the form of a probability density function which takes into account uncertainties in observations and prior information, and which we sample using a Markov chain Monte Carlo algorithm. Applying the technique to Kīlauea Volcano, we develop a model which relates magma flow rates with deformation of the volcano's surface, sulfur dioxide emission rates, lava flow field volumes, and composition of the volcano's basaltic magma. This model accounts for effects and processes mostly neglected in previous supply rate estimates at Kīlauea, including magma compressibility, loss of sulfur to the hydrothermal system, and potential magma storage in the volcano's deep rift zones. We jointly invert data and prior information to estimate rates of supply, storage, and eruption during three recent quasi-steady-state periods at the volcano. Results shed new light on the time-variability of magma supply to Kīlauea, which we find to have increased by 35–100% between 2001 and 2006 (from 0.11–0.17 to 0.18–0.28 km 3 /yr), before subsequently decreasing to 0.08–0.12 km 3 /yr by 2012. Changes in supply rate directly impact hazard at the volcano, and were largely responsible for an increase in eruption rate of 60–150% between 2001 and 2006, and subsequent decline by as much as 60% by 2012. We also demonstrate the occurrence of temporal changes in the proportion of Kīlauea's magma supply that is stored versus erupted, with the supply “surge” in 2006 associated with increased accumulation of magma at the summit. Finally, we are able to place some constraints on sulfur concentrations in Kīlauea magma and the scrubbing of sulfur by the volcano's hydrothermal system. Multiphysical, Bayesian constraint on magma flow rates may be used to monitor evolving volcanic hazard not just at Kīlauea but at other volcanoes around the world.

Hawai'i↗

USGS2016

Insights into methane dynamics from analysis of authigenic carbonates and chemosynthetic mussels at newly-discovered Atlantic Margin seeps

The recent discovery of active methane venting along the US northern and mid-Atlantic margin represents a new source of global methane not previously accounted for in carbon budgets from this region. However, uncertainty remains as to the origin and history of methane seepage along this tectonically inactive passive margin. Here we present the first isotopic analyses of authigenic carbonates and methanotrophic deep-sea mussels, Bathymodiolus sp., and the first direct constraints on the timing of past methane emission, based on samples collected at the upper slope Baltimore Canyon (&sim;385 m water depth) and deepwater Norfolk (&sim;1600 m) seep fields within the area of newly-discovered venting. The authigenic carbonates at both sites were dominated by aragonite, with an average signature of −47&permil; , a value consistent with microbially driven anaerobic oxidation of methane-rich fluids occurring at or near the sediment–water interface. Authigenic carbonate U and Sr isotope data further support the inference of carbonate precipitation from seawater-derived fluids rather than from formation fluids from deep aquifers. Carbonate stable and radiocarbon ( and ) isotope values from living Bathymodiolus sp. specimens are lighter than those of seawater dissolved inorganic carbon, highlighting the influence of fossil carbon from methane on carbonate precipitation. U–Th dates on authigenic carbonates suggest seepage at Baltimore Canyon between 14.7±0.6 ka to 15.7±1.6 ka , and at the Norfolk seep field between 1.0±0.7 ka to 3.3±1.3 ka , providing constraint on the longevity of methane efflux at these sites. The age of the brecciated authigenic carbonates and the occurrence of pockmarks at the Baltimore Canyon upper slope could suggest a link between sediment delivery during Pleistocene sea-level lowstand, accumulation of pore fluid overpressure from sediment compaction, and release of overpressure through subsequent venting. Calculations show that the Baltimore Canyon site probably has not been within the gas hydrate stability zone (GHSZ) in the past 20 ka, meaning that in-situ release of methane from dissociating gas hydrate cannot be sustaining the seep. We cannot rule out updip migration of methane from dissociation of gas hydrate that occurs farther down the slope as a source of the venting at Baltimore Canyon, but consider that the history of rapid sediment accumulation and overpressure may play a more important role in methane emissions at this site.

Baltimore Canyon seep fields, Norfolk seep fields↗

USGS2016

Large along-strike variations in the onset of Subandean exhumation: Implications for Central Andean orogenic growth

Plate tectonics drives mountain building in general, but the space-time pattern and style of deformation is influenced by how climate, geodynamics, and basement structure modify the orogenic wedge. Growth of the Subandean thrust belt, which lies at the boundary between the arid, high-elevation Central Andean Plateau and its humid, low-elevation eastern foreland, figures prominently into debates of orogenic wedge evolution. We integrate new apatite and zircon (U-Th)/He thermochronometer data with previously published apatite fission-track data from samples collected along four Subandean structural cross-sections in Bolivia between 15° and 20°S. We interpret cooling ages vs. structural depth to indicate the onset of Subandean exhumation and signify the forward propagation of deformation. We find that Subandean growth is diachronous south (11 ± 3 Ma) vs. north (6 ± 2 Ma) of the Bolivian orocline and that Subandean exhumation magnitudes vary by more than a factor of two. Similar north-south contrasts are present in foreland deposition, hinterland erosion, and paleoclimate; these observations both corroborate diachronous orogenic growth and illuminate potential propagation mechanisms. Of particular interest is an abrupt shift to cooler, more arid conditions in the Altiplano hinterland that is diachronous in southern Bolivia (16-13 Ma) vs. northern Bolivia (10-7 Ma) and precedes the timing of Subandean propagation in each region. Others have interpreted the paleoclimate shift to reflect either rapid surface uplift due to lithosphere removal or an abrupt change in climate dynamics once orographic threshold elevations were exceeded. These mechanisms are not mutually exclusive and both would drive forward propagation of the orogenic wedge by augmenting the hinterland backstop, either through surface uplift or spatially variable erosion. In summary, we suggest that diachronous Subandean exhumation was driven by piecemeal hinterland uplift, orography, and the outward propagation of deformation.

Andes↗

USGS2016

The timing of compositionally-zoned magma reservoirs and mafic 'priming' weeks before the 1912 Novarupta-Katmai rhyolite eruption

The June 6, 1912 eruption of more than 13 km3 of dense rock equivalent (DRE) magma at Novarupta vent, Alaska was the largest of the 20th century. It ejected >7 km3 of rhyolite, ~1.3 km3 of andesite and ~4.6 km3 of dacite. Early ideas about the origin of pyroclastic flows and magmatic differentiation (e.g., compositional zonation of reservoirs) were shaped by this eruption. Despite being well studied, the timing of events that led to the chemically and mineralogically zoned magma reservoir remain poorly known. Here we provide new insights using the textures and chemical compositions of plagioclase and orthopyroxene crystals and by reevaluating previous U-Th isotope data. Compositional zoning of the magma reservoir likely developed a few thousand years before the eruption by several additions of mafic magma below an extant silicic reservoir. Melt compositions calculated from Sr contents in plagioclase fill the compositional gap between 68 and 76% SiO2 in whole pumice clasts, consistent with uninterrupted crystal growth from a continuum of liquids. Thus, our findings support a general model in which large volumes of crystal-poor rhyolite are related to intermediate magmas through gradual separation of melt from crystal-rich mush. The rhyolite is incubated by, but not mixed with, episodic recharge pulses of mafic magma that interact thermochemically with the mush and intermediate magmas. Hot, Mg-, Ca-, and Al-rich mafic magma intruded into, and mixed with, deeper parts of the reservoir (andesite and dacite) multiple times. Modeling the relaxation of the Fe-Mg concentrations in orthopyroxene and Mg in plagioclase rims indicates that the final recharge event occurred just weeks prior to the eruption. Rapid addition of mass, volatiles, and heat from the recharge magma, perhaps aided by partial melting of cumulate mush below the andesite and dacite, pressurized the reservoir and likely propelled a ~10 km lateral dike that allowed the overlying rhyolite to reach the surface.

Earth and Planetary Science Letters↗

USGS2016

Interaction between climate, volcanism, and isostatic rebound in Southeast Alaska during the last deglaciation

Observations of enhanced volcanic frequency during the last deglaciation have led to the hypothesis that ice unloading in glaciated volcanic terrains can promote volcanism through decompression melting in the shallow mantle or a reduction in crustal magma storage time. However, a direct link between regional climate change, isostatic adjustment, and the initiation of volcanism remains to be demonstrated due to the difficulty of obtaining high-resolution well-dated records that capture short-term climate and volcanic variability traced to a particular source region. Here we present an exceptionally resolved record of 19 tephra layers paired with foraminiferal oxygen isotopes and alkenone paleotemperatures from marine sediment cores along the Southeast Alaska margin spanning the last deglacial transition. Major element compositions of the tephras indicate a predominant source from the nearby Mt. Edgecumbe Volcanic Field (MEVF). We constrain the timing of this regional eruptive sequence to 14.6–13.1 ka. The sudden increase in volcanic activity from the MEVF coincides with the onset of Bølling–Allerød interstadial warmth, the disappearance of ice-rafted detritus, and rapid vertical land motion associated with modeled regional isostatic rebound in response to glacier retreat. These data support the hypothesis that regional deglaciation can rapidly trigger volcanic activity. Rapid sea surface temperature fluctuations and an increase in local salinity (i.e., δ 18 O sw ) variability are associated with the interval of intense volcanic activity, consistent with a two-way interaction between climate and volcanism in which rapid volcanic response to ice unloading may in turn enhance short-term melting of the glaciers, plausibly via albedo effects on glacier ablation zones.

Earth and Planetary Science Letters↗

USGS2016

Joint analysis of geodetic and earthquake fault-plane solution data to constrain magmatic sources: A case study from Kīlauea Volcano

A joint analysis of geodetic and seismic datasets from Kīlauea Volcano during a period of magmatic unrest in 2006 demonstrates the effectiveness of this combination for testing and constraining models of magma dynamics for a complex, multi-source system. At the end of 2003, Kīlauea's summit began a four-year-long period of inflation due to a surge in magma supply to the volcano. In 2006, for the first time since 1982, Kīlauea's Southwest Rift Zone (SWRZ) also experienced inflation. To investigate the characteristics of active magma sources and the nature of their interactions with faults in the SWRZ during 2006, we integrate, through Coulomb stress modeling, contemporary geodetic data from InSAR and GPS with a new catalogue of double-couple fault-plane solutions for volcano-tectonic earthquakes. We define two periods of inflation during 2006 based on the rate of deformation measured in daily GPS data, spanning February to 15 March 2006 (Period 1) and 16 March to 30 September 2006 (Period 2). InSAR data for these two periods are inverted to determine the position, change in size, and shape of inflation sources in each period. Our new models are consistent with microseismic activity from each period. They suggest that, during Period 1, deformation in the SWRZ can be explained by pressurization of magma in a spherical reservoir beneath the south caldera, and that, during Period 2, magma was also aseismically intruded farther to the southwest into the SWRZ along a sub-horizontal plane. Our Coulomb stress analysis shows that the microseismicity recorded in the SWRZ is induced by overpressurization of the south caldera reservoir, and not by magma intrusion into the SWRZ. This study highlights the importance of a joint analysis of independent geophysical datasets to fully constrain the nature of magma accumulation.

Hawaii↗

USGS2016

Comparison of the mineral composition of the sediment found in two Mars dunefields: Ogygis Undae and Gale crater – three distinct endmembers identified

The composition of two dune fields, Ogygis Undae and the NE–SW trending dune field in Gale crater (the “Bagnold Dune Field” and “Western Dune Field”), were analyzed using thermal emission spectra from the Mars Global Surveyor (MGS) Thermal Emission Spectrometer (TES) and the Mars Odyssey Thermal Emission Imaging System (THEMIS). The Gale crater dune field was used as a baseline as other orbital compositional analyses have been conducted, and in situ sampling results will soon be available. Results from unmixing thermal emission spectra showed a spatial variation between feldspar mineral abundances and pyroxene mineral abundances in Ogygis Undae. Other datasets, including nighttime thermal inertia values, also showed variation throughout the dune field. One explanation proposed for this variation is a bimodal distribution of two sand populations. This distribution is seen in some terrestrial dune fields. The two dune fields varied in both mineral types present and in uniformity of composition. These differences point to different source lithologies and different distances travelled from source material. Examining these differences further will allow for a greater understanding of aeolian processes on Mars.

Earth and Planetary Science Letters↗

USGS2017

Lithospheric density structure beneath the Tarim basin and surroundings, northwestern China, from the joint inversion of gravity and topography

Intraplate strain generally focuses in discrete zones, but despite the profound impact of this partitioning on global tectonics, geodynamics, and seismic hazard, the processes by which deformation becomes localized are not well understood. Such heterogeneous intraplate strain is exemplified in central Asia, where the Indo-Eurasian collision has caused widespread deformation while the Tarim block has experienced minimal Cenozoic shortening. The apparent stability of Tarim may arise either because strain is dominantly accommodated by pre-existing faults in the continental suture zones that bound it—essentially discretizing Eurasia into microplates—or because the lithospheric-scale strength (i.e., viscosity) of the Tarim block is greater than its surroundings. Here, we jointly analyze seismic velocity, gravity, topography, and temperature to develop a 3-D density model of the crust and upper mantle in this region. The Tarim crust is characterized by high density, v s , v p , and v p /v s , consistent with a dominantly mafic composition and with the presence of an oceanic plateau beneath Tarim. Low-density but high-velocity mantle lithosphere beneath southern (southwestern) Tarim underlies a suite of Permian plume-related mafic intrusions and A-type granites sourced in previously depleted mantle lithosphere; we posit that this region was further depleted, dehydrated, and strengthened by Permian plume magmatism. The actively deforming western and southern margins of Tarim—the Tien Shan, Kunlun Shan, and Altyn Tagh fault—are underlain by buoyant upper mantle with low velocity; we hypothesize that this material has been hydrated by mantle-derived fluids that have preferentially migrated along Paleozoic continental sutures. Such hydrous material should be weak, and herein strain focuses there because of lithospheric-scale variations in rheology rather than the pre-existence of faults in the brittle crust. Thus this world-class example of strain partitioning arises not simply from the pre-existence of brittle faults but from the thermo-chemical and therefore rheological variations inherited from prior tectonism.

Earth and Planetary Science Letters↗

USGS2017

Pulsed strain release on the Altyn Tagh fault, northwest China

Earthquake recurrence models assume that major surface-rupturing earthquakes are followed by periods of reduced rupture probability as stress rebuilds. Although purely periodic, time- or slip-predictable rupture models are known to be oversimplifications, a paucity of long records of fault slip clouds understanding of fault behavior and earthquake recurrence over multiple ruptures. Here, we report a 16 kyr history of fault slip—including a pulse of accelerated slip from 6.4 to 6.0 ka—determined using a Monte Carlo analysis of well-dated offset landforms along the central Altyn Tagh strike-slip fault (ATF) in northwest China. This pulse punctuates a median rate of 8.1 +1.2 / −0.9 mm/a and likely resulted from either a flurry of temporally clustered ∼Mw 7.5 ground-rupturing earthquakes or a single large >Mw 8.2 earthquake. The clustered earthquake scenario implies rapid re-rupture of a fault reach >195 km long and indicates decoupled rates of elastic strain energy accumulation versus dissipation, conceptualized as a crustal stress battery. If the pulse reflects a single event, slip-magnitude scaling implies that it ruptured much of the ATF with slip similar to, or exceeding, the largest documented historical ruptures. Both scenarios indicate fault rupture behavior that deviates from classic time- or slip-predictable models.

Altyn Tagh fault↗