USGS Science⌕ Search

SEARCH · USGS Science

Results for “Evolution”

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

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

At least 1,783 records · Page 99Linked to original sources

Middle and Late Pleistocene glaciations in the southwestern Pamir and their effects on topography

Glacial chronologies provide insight into the evolution of paleo-landscapes, paleoclimate, topography, and the erosion processes that shape mountain ranges. In the Pamir of Central Asia, glacial morphologies and deposits indicate extensive past glaciations, whose timing and extent remain poorly constrained. Geomorphic data and 15 new 10 Be exposure ages from moraine boulders and roches moutonnées in the southwestern Pamir document multiple Pleistocene glacial stages. The oldest exposure ages, , underestimate the age of the earliest preserved glacial advance and imply that the modern relief of the southwestern Pamir (peaks at ∼5000–6000 m a.s.l.; valleys at ∼2000–3000 m a.s.l.) already existed in the late Middle Pleistocene. Younger exposure ages (∼40–80 ka, ∼30 ka) complement the existing Central Asian glacial chronology and reflect successively less extensive Late Pleistocene glaciations. The topography of the Pamir and the glacial chronologies suggest that, in the Middle Pleistocene, an ice cap or ice field occupied the eastern Pamir high-altitude plateau, whereas westward flowing valley glaciers incised the southwestern Pamir. Since the Late Pleistocene deglaciation, the rivers of the southwestern Pamir adjusted to the glacially shaped landscape. Localized rapid fluvial incision and drainage network reorganization reflect the transient nature of the deglaciated landscape.

Earth and Planetary Science Letters↗

Pliocene erosional pulse and glacier-landscape feedbacks in the western Alaska Range

Pliocene–Pleistocene glaciation modified the topography and erosion of most middle- and high-latitude mountain belts, because the evolution of catchment topography controls long-term glacier mass balance and erosion. Hence, characterizing how erosion rates change during repeated glaciations can help test hypothesized glacier erosion-landscape feedbacks across a range of settings. To better understand how glaciations and landscapes coevolve on geologic timescales, I quantify erosion rates in the glaciated western Alaska Range with low-temperature thermochronometric data and modeling. Zircon (U–Th)/He and apatite fission track data suggest mountain-building was underway by the early Miocene. In contrast, lower-temperature apatite (U–Th)/He age-elevation and grain age-kinetic data indicate that erosion accelerated coincident with regional Pliocene glaciation ca. 4 Ma. Furthermore, erosion rates calculated within an eroding half-space indicate slow erosion at a rate ≤0.3 km/m.y. before 4.2 Ma, an initial pulse of rapid erosion at a rate of 1.0–1.6 km/m.y. during 4.2–2.9 Ma, and more moderate erosion at a rate of 0.4–0.7 km/m.y. since 2.9 Ma. The initial erosion pulse suggests a significant transient landscape adjustment to the introduction of efficient glacial erosion. The subsequent decrease in Pleistocene erosion rates is consistent with a negative feedback between continuing glaciation and glacier size/erosivity: If glacial erosion outpaces rock uplift, glacier erosion decreases over time as topography, mass balance, valley gradients, and ice flux are reduced. These findings imply that in areas of moderate rock uplift rates, the onset of local Plio–Pleistocene glaciation may have been punctuated by an initial pulse of rapid landscape change, after which change became more gradual.

Alaska↗

Physicochemical models of effusive rhyolitic eruptions constrained with InSAR and DEM data: A case study of the 2011-2012 Cordon Caulle eruption

The 9 month long 2011-2012 eruption of Cordon Caulle (Southern Andes, Chile) is the best instrumentally recorded rhyolitic eruption to date and the first time that the effusion of a rhyolitic flow has been observed in detail. We use Interferometric Synthetic Aperture Radar (InSAR), with time-lapse DEMs and numerical models to study the dynamics of coupled magma reservoir deflation and lava effusion. InSAR recorded ~2.2-2.5 m of deflation after the first three days of the eruption, which can be modeled using a spheroidal magma reservoir at a depth of ∼5 km, ∼20 km long, and with a pressure drop of ~20-30 MPa. The source is elongated in the NW-SE direction and its large dimensions imply a large plumbing system spanning neighboring volcanoes and active throughout the eruption, with a slight change halfway through the effusive phase. TanDEM-X DEMs record the extrusion of both the rhyolitic lava flow and the intrusion of a shallow laccolith around the eruptive vent, with a total volume of ~1.2 km3 DRE. The laccolith was emplaced during the first month of the eruption, during both the eruption explosive and effusive stages. Both the reservoir pressure drop and the extruded volume time series follow quasi-exponential trends, and can be explained by a model that couples the reservoir pres- sure decrease, time- and pressured ependent variations in the magma properties inside of the reservoir, and conduit flow. This model predicts both the temporal evolution and amplitude of both time series, and a magma compressibility of ∼10^−10 Pa−1, half the compressibility of the magma of the sub-Plinian explosive phase. Further, we estimate that the reservoir contained 1-3 wt.% dissolved H2O at the onset of lava effusion, with no exsolved CO2 and H2O in the reservoir throughout the eruption. This is in accord with a magma that was significantly degassed after the explosive phase. These remaining volatiles might have been responsible for magma fragmentation, consistent with the hybrid explosive and effusive style observed during the waning of the eruption.

Cordón Caulle↗

The role of preexisting upper plate strike-slip faults during long-lived (ca. 30 Myr) oblique flat slab subduction, southern Alaska

Upper plates of subduction zones commonly respond to flat slab subduction by structural reactivation, magmatic arc disruption, and foreland basin inversion. However, the role of active strike-slip faults in focusing convergent deformation and magmatism in response to oblique flat slab subduction remains less clear. Here, we present new detrital apatite fission-track (dAFT) ages from 12 modern catchments in the eastern Alaska Range, Alaska, USA, to reveal how the dextral Denali fault system has facilitated bedrock exhumation and topographic growth during ca. 30 Ma-to-present oblique flat slab subduction of the Yakutat oceanic plateau. Additionally, a 940 ka ( 40 Ar/ 39 Ar whole rock) basalt flow is spatially associated with Cenozoic structures, locally reset AFT ages and provides the first evidence for Quaternary volcanism along the southern flank of the eastern Alaska Range. We integrate our new data with other thermochronologic, geochronologic, and regional geologic datasets to show that (1) most high topography regions in southern Alaska have undergone rapid bedrock cooling and exhumation since ca. 30 Ma; (2) elevated terrain and young cooling are spatially associated with long-lived active strike-slip fault systems; (3) topographic growth associated with strike-slip fault deformation led to local inversion of basin systems and drainage reorganization; (4) the onset of oblique oceanic plateau subduction is coeval with a southward shift in arc magmatism from one region of active strike-slip faulting to another above the northeastern edge of the flat slab; and (5) Quaternary volcanism marks the revival of magmatism in the eastern Alaska Range above the geophysically imaged northeastern edge of the flat slab. Our analysis of the post-30 Ma geologic evolution of southern Alaska demonstrates that strike-slip fault systems that were active at the time of slab flattening evolved into transpression zones that focused bedrock cooling, rock exhumation, and topographic growth.

Alaska↗

Can we accurately estimate sediment budgets on Mars?

Sediment budgets are fundamentally important for planetary science. However, only one primary method, based on remote sensing, is currently available for determining extraterrestrial sediment budgets. For determining sediment budgets on Earth, both in-situ and remote sensing methods are available. Despite the widespread use of the two methods, there has been surprisingly little research on how well the sediment budgets produced by these two approaches reconcile with one another, which highlights the lack of quantitative understanding of errors for sediment budgets measured with remote sensing in planetary research. Therefore, there is a general need to expand our knowledge of sediment budgets. Here we use a background review and analog case study of an aeolian dunefield in Grand Canyon, Earth to frame a path forward for addressing shortcomings of remote sensing sediment budgets on Mars. We estimate a 53% percent difference in the sediment budget determined with remote sensing relative to in-situ methods for a simple endmember scenario of a dunefield within a unimodal wind directional regime and no external sediment supply. However, when we incorporated key sources of uncertainty in remote sensing change detection following methods commonly used by geomorphologists on Earth, the estimates of sediment budget differences relative to the in-situ method spanned a much larger range, from 3% to 138%. Our case study also suggests that sediment budget errors could be much larger under more complex wind direction, sediment supply, and physiographic settings, and that variability in those landscape characteristics might be used to better estimate errors for dunefield sediment budgets. We conclude that by comparing sediment budgets derived from in-situ measurements of sediment fluxes and from remote sensing measurements at many more analog sites on Earth, the aeolian research community, and the geomorphology discipline, could gain an understanding of the errors of the remote sensing method, which is used by investigators on other planetary bodies such as Mars. This could improve the ability to quantify sediment budgets on Mars – and, in the future, other planetary environments where high-resolution topographic data are available – as well as directly improve our ability to interpret extraterrestrial landscape evolution related to climate, weather, and geologic history.

Earth and Planetary Science Letters↗

Geochemical evidence for evolving Proterozoic crustal thickness and orogenic styles in southwestern Laurentia

It has long been challenging for researchers to track the crustal thickness and mode(s) of crustal modification in ancient convergent margins, limiting evaluation of the tectonic styles and processes that modify continental crust during orogenesis. We present trace element igneous geochemical crustal thickness proxies that quantitatively track the crustal thickness evolution of the long-lived Proterozoic active margin in the southwestern U.S.A. We integrate these results with geobarometric data to constrain the mode of crustal modification. The data indicate a complex record of crustal thickness change in space and time and evolving orogenic styles. Geochemical proxies at 1.84–1.72 Ga are consistent with 20–40 km thick magmatic arcs that were locally thickened to ∼50 km during ∼1.75 Ga tectonism. During the Yavapai orogeny, 1.72–1.69 Ga, a ∼200-km-wide belt of 50-60 km thick crust extended from southern California to northern Colorado and was rapidly thinned and exhumed by ∼1.68 Ga. Crustal thickening and thinning during the Yavapai orogeny largely occurred by shortening and exhumation, respectively, in the upper 25 km of the crust. Subsequent 1.68–1.60 Ga tectonism involved crustal growth, local crustal thickening, and low-P, high-T metamorphism, consistent with extensional accretionary orogenesis. The 1.47–1.37 Ga Picuris orogeny was associated with 50–60 km thick crust across much of the Southwest and involved crustal shortening with ∼10 km of magmatic underplating. Advective heat from the emplacement of ferroan granites in the mid-crust likely contributed to elevated geothermal gradients and rheologically weakened the crust. Our results suggest evolving orogenic styles in the Southwest from 1.75–1.69 Ga short-lived crustal thickening associated with terrane accretion to 1.69–1.60 Ga largely extensional accretionary orogenesis, and regional, long-lived crustal thickening at 1.47–1.37 Ga involving extensive basaltic underplating. Contrasting with some recent hypotheses, our data document a complex middle Proterozoic record for the Southwest that was not orogenically quiescent or tectonically stagnant but involved complex mountain building styles.

Earth and Planetary Science Letters↗

Seismicity zoning at Coso geothermal field and stress changes from fluid production and migration

The Coso geothermal field is a major geothermal power production site in the western United States. It has been observed that low-magnitude seismic events ( M < 3.71) are unevenly distributed in three distinct zones, namely, nearfield (<3 km), midfield (3–6 km), and farfield (> 6 km) from the Coso geothermal plant. These zones exhibit distinct changes in earthquake location before and during geothermal production episodes that began in 1986. After 1986, the midfield region of the main flank experiences a significantly lower seismicity rate than the surrounding areas before production episodes. During 2014–2019, the farfield earthquakes cluster in the eastern and western parts of the greater Coso area, which is discernably different from how those pre-production earthquake events were distributed along the conjugate NW-SE and SW-NW trending structures across the main flank. Here, we analyze the stage of stress with finite-element-based poroelastic simulations to illustrate how the spatiotemporal evolution of the seismicity is associated with the pattern of stress perturbations caused by fluid migration amid the operations of geothermal power plants. Generally, ∼70% of co-production seismicity is found in zones of increased Coulomb stress between 2014 and 2019 at >99% confidence. Meanwhile, the midfield zone of seismic paucity overlaps with the zone of decreasing pore-fluid pressure. Overall, the results provide a physical explanation of how decadal geothermal operations at Coso have perturbed stress-field changes and contributed to the evolving characteristic seismic pattern, shedding insights into assessing the seismic hazard in other geothermal settings.

California↗

Dynamic drainage reorganization in Eastern Tibet: Insights from the Yangtze River first bend

The modern drainage network of eastern Tibet is widely believed to have developed through a series of river capture and flow reversal events; however, the timing and mechanisms driving this reorganization remain contentious. Among these events, the river capture that formed the First Bend of the Yangtze River (YFB) stands out as both iconic and particularly debated. Here we present sedimentary provenance data from the Late Miocene–Quaternary Dali Basin, located south of the YFB, which indicate that a southward-flowing Jinsha River (i.e., the present-day upper Yangtze River) sourced sediment to the Dali basin at ∼7.4–6.4 Ma in a drainage configuration different from that of today. Because this interval postdates the initial establishment of a near-modern Jinsha River system prior to the Miocene, our results imply at least two discrete fluvial reorganizations occurred at the YFB—one preceding ∼7.4 Ma and another following ∼6.4 Ma. By integrating these findings with landscape evolution modeling, we infer that the initiation of rapid uplift of the Yulong-Haba Mountains and the Diancang Shan may have been responsible for these drainage reorganizations. These results underscore that Cenozoic drainage systems on the eastern Tibetan Plateau have evolved dynamically on a short timescale of ∼10 5 –10 6 -year, rather than remaining in a long-term stationary configuration on ∼10 7 -year timescales.

eastern Tibetan Plateau, first bend of the Yangtze↗

Recycling radiogenic osmium by crustal foundering in subduction zones: Evidence from pyroxenite xenoliths in the north Andean arc

Removal of mafic-ultramafic lower crust (e.g., via delamination) is fundamental to making andesitic continental crust, yet direct evidence of this process remains elusive. A unique suite of garnet clinopyroxenite and hornblendite (arclogite) xenoliths from the northern volcanic zone of the Andes, erupted to the surface in the mid-Pleistocene, have bulk-rock geochemistry and osmium (Os) isotopic compositions that demonstrate their crustal affinity, yet have equilibration pressures and temperatures below the arc Mohorovičić discontinuity (~53 km) and some as deep as within the sub-arc mantle wedge (~105 km). Garnet websterites from the same xenolith suite, sourced from depths approaching the Wadati-Benioff zone (~140 km), have elemental and isotopic compositions indicative of a mantle origin, likely formed as products of peridotite-liquid reactions above the subducting slab. Variability in bulk-rock 143 Nd/ 144 Nd values and garnet oxygen isotope ratios ( 𝛿 18 O) for these samples is attributed to minor assimilation of subducted components and/or older crustal material, but assimilation alone cannot account for the highly radiogenic Os isotopic compositions. These results provide direct petrologic evidence for modern arclogite formation and foundering in the Andean orogen, the archetypal active continental subduction system. Elevated 187 Os/ 188 Os values of Mercaderes arclogites imply that recycling of gravitationally unstable arc cumulates during continental crust formation introduces highly radiogenic Os into the convective mantle, which has implications for mass fluxes across the crust-mantle boundary as well as the sources and evolution of mantle heterogeneity as seen in orogenic peridotite massifs and in the sources of oceanic basalts

Earth and Planetary Science Letters↗

Comprehensive characterization of Mauna Loa basalt rheology based on the 2022 eruption

Rheological data are crucial for the understanding of storage, transport, eruption, and emplacement of magma and lava. While viscosity estimates for lavas from Mauna Loa, the largest active volcano on Earth, exist, no direct rheological measurements have been done to date. This is especially surprising given the role Mauna Loa lavas and resulting landforms play in the interpretation of volcanic processes across our solar system and beyond. We present the first complete rheological characterization of Mauna Loa lava, based on the 2022 eruption. We constrain the melt viscosity and map the lava's rheology during crystallization at isothermal subliquidus conditions and at cooling at rates between 0.25 and 3.00 °C/minute. These experiments were performed at oxygen fugacities relevant to storage and eruption conditions at Mauna Loa (log fO 2 = -8.7). We integrate textural and thermal characteristics obtained from the experiments with those from water quenched natural samples (i.e. sampled from actively flowing lava and quenched within seconds to preserve the textural state of the lava while flowing) to reconstruct the evolution of lava rheology during emplacement. This integration suggests that the Mauna Loa lavas were emplaced at viscosities between 10 1.5 to 10 4.5 Pa s over the temperature range of 1150 – 1090 °C before entering a region of rapid viscosity increase and effective solidification caused by groundmass crystallization. We provide a detailed comparison of our data to rheological reconstructions of the 1984 eruption. The resulting rheological map can help guide physical property-based magma ascent and lava emplacement models, as well as the interpretation of flow morphologies and the conditions of formation of basaltic landforms on other planets.

Hawaii↗

Mechanisms of forest resilience

Ecosystems are dynamic systems with complex responses to environmental variation. In response to pervasive stressors of changing climate and disturbance regimes, many ecosystems are realigning rapidly across spatial scales, in many cases moving outside of their observed historical range of variation into alternative ecological states. In some cases, these new states are transitory and represent successional stages that may ultimately revert to the pre-disturbance condition; in other cases, alternative states are persistent and potentially self-reinforcing, especially under conditions of altered climate, disturbance regimes, and influences of non-native species. These reorganized states may appear novel, but reorganization is a characteristic ecosystem response to environmental variation that has been expressed and documented throughout the paleoecological record. Resilience, the ability of an ecosystem to recover or adapt following disturbance, is an emergent property that results from the expression of multiple mechanisms operating across levels of organism, population, and community. We outline a unifying framework of ecological resilience based on ecological mechanisms that lead to outcomes of persistence, recovery, and reorganization. Persistence is the ability of individuals to tolerate exposure to environmental stress, disturbance, or competitive interactions. As a direct expression of life history evolution and adaptation to environmental variation and stress, persistence is manifested most directly in survivorship and continued growth and reproduction of established individuals. When persistence has been overcome (e.g., following mortality from stress, disturbance, or both), populations must recover by reproduction. Recovery requires the establishment of new individuals from seed or other propagules following dispersal from the parent plant. When recovery fails to re-establish the pre-disturbance community, the ecosystem will assemble into a new state. Reorganization occurs along a gradient of magnitude, from changes in the relative dominance of species present in a community, to individual species replacements within an essentially intact community, to complete species turnover and shift to dominance by plants of different functional types, e.g. transition from forest to shrub or grass dominance. When this latter outcome is persistent and involves reinforcing mechanisms, the resulting state represents a vegetation type conversion (VTC), which in this framework represents an end member of reorganization processes. We explore reorganization in greater detail as this phase is increasingly observed but the least understood of the resilience responses. This resilience framework provides a direct and actionable basis for ecosystem management in a rapidly changing world, by targeting specific components of ecological response and managing for sustainable change.

Forest Ecology and Management↗

Organic nitrogen chemistry during low-grade metamorphism

Most of the organic nitrogen (Norg) on Earth is disseminated in crustal sediments and rocks in the form of fossil nitrogen-containing organic matter. The chemical speciation of fossil Norg within the overall molecular structure of organic matter changes with time and heating during burial. Progressive thermal evolution of organic matter involves phases of enhanced elimination of Norg and ultimately produces graphite containing only traces of nitrogen. Long-term chemical and thermal instability makes the chemical speciation of Norg a valuable tracer to constrain the history of sub-surface metamorphism and to shed light on the subsurface biogeochemical nitrogen cycle and its participating organic and inorganic nitrogen pools. This study documents the evolutionary path of Norg speciation, transformation and elimination before and during metamorphism and advocates the use of X-ray photoelectron spectroscopy (XPS) to monitor changes in Norg speciation as a diagnostic tool for organic metamorphism. Our multidisciplinary evidence from XPS, stable isotopes, traditional quantitative coal analyses, and other analytical approaches shows that at the metamorphic onset Norg is dominantly present as pyrrolic and pyridinic nitrogen. The relative abundance of nitrogen substituting for carbon in condensed, partially aromatic systems (where N is covalently bonded to three C atoms) increases exponentially with increasing metamorphic grade, at the expense of pyridinic and pyrrolic nitrogen. At the same time, much Norg is eliminated without significant nitrogen isotope fractionation. The apparent absence of Rayleigh-type nitrogen isotopic fractionation suggests that direct thermal loss of nitrogen from an organic matrix does not serve as a major pathway for Norg elimination. Instead, we propose that hot H, O-containing fluids or some of their components gradually penetrate into the carbonaceous matrix and eliminate Norg along a progressing reaction front, without causing nitrogen isotope fractionation in the residual Norg in the unreacted core of the carbonaceous matrix. Before the reaction front can reach the core, an increasing part of core Norg chemically stabilizes in the form of nitrogen atoms substituting for carbon in condensed, partially aromatic systems forming graphite-like structural domains with delocalized ??-electron systems (nitrogen atoms substituting for "graphitic" carbon in natural metamorphic organic matter). Thus, this nitrogen species with a conservative isotopic composition is the dominant form of residual nitrogen at higher metamorphic grade. ?? 2007 Elsevier Ltd. All rights reserved.

Geochimica et Cosmochimica Acta↗

Modeling aluminum-silicon chemistries and application to Australian acidic playa lakes as analogues for Mars

Recent Mars missions have stimulated considerable thinking about the surficial geochemical evolution of Mars. Among the major relevant findings are the presence in Meridiani Planum sediments of the mineral jarosite (a ferric sulfate salt) and related minerals that require formation from an acid-salt brine and oxidizing environment. Similar mineralogies have been observed in acidic saline lake sediments in Western Australia (WA), and these lakes have been proposed as analogues for acidic sedimentary environments on Mars. The prior version of the equilibrium chemical thermodynamic FREZCHEM model lacked Al and Si chemistries that are needed to appropriately model acidic aqueous geochemistries on Earth and Mars. The objectives of this work were to (1) add Al and Si chemistries to the FREZCHEM model, (2) extend these chemistries to low temperatures (<0 ??C), if possible, and (3) use the reformulated model to investigate parallels in the mineral precipitation behavior of acidic Australian lakes and hypothetical Martian brines. FREZCHEM is an equilibrium chemical thermodynamic model parameterized for concentrated electrolyte solutions using the Pitzer approach for the temperature range from <-70 to 25 ??C and the pressure range from 1 to 1000 bars. Aluminum chloride and sulfate mineral parameterizations were based on experimental data. Aluminum hydroxide and silicon mineral parameterizations were based on Gibbs free energy and enthalpy data. New aluminum and silicon parameterizations added 12 new aluminum/silicon minerals to this Na-K-Mg-Ca-Fe(II)-Fe(III)-Al-H-Cl-Br-SO 4 -NO 3 -OH-HCO 3 -CO 3 -CO 2 -O 2 -CH 4 -Si-H 2 O system that now contain 95 solid phases. There were similarities, differences, and uncertainties between Australian acidic, saline playa lakes and waters that likely led to the Burns formation salt accumulations on Mars. Both systems are similar in that they are dominated by (1) acidic, saline ground waters and sediments, (2) Ca and/or Mg sulfates, and (3) iron precipitates such as jarosite and hematite. Differences include: (1) the dominance of NaCl in many WA lakes, versus the dominance of Fe-Mg-Ca-SO 4 in Meridiani Planum, (2) excessively low K + concentrations in Meridiani Planum due to jarosite precipitation, (3) higher acid production in the presence of high iron concentrations in Meridiani Planum, and probably lower rates of acid neutralization and hence, higher acidities on Mars owing to colder temperatures, and (4) lateral salt patterns in WA lakes. The WA playa lakes display significant lateral variations in mineralogy and water chemistry over short distances, reflecting the interaction of acid ground waters with neutral to alkaline lake waters derived from ponded surface runoff. Meridiani Planum observations indicate that such lateral variations are much less pronounced, pointing to the dominant influence of ground water chemistry, vertical ground water movements, and aeolian processes on the Martian surface mineralogy. ?? 2009 Elsevier Ltd.

Geochimica et Cosmochimica Acta↗

Geochemistry of hydrothermal fluids from the PACMANUS, Northeast Pual and Vienna Woods hydrothermal fields, Manus Basin, Papua New Guinea

Processes controlling the composition of seafloor hydrothermal fluids in silicic back-arc or near-arc crustal settings remain poorly constrained despite growing evidence for extensive magmatic–hydrothermal activity in such environments. We conducted a survey of vent fluid compositions from two contrasting sites in the Manus back-arc basin, Papua New Guinea, to examine the influence of variations in host rock composition and magmatic inputs (both a function of arc proximity) on hydrothermal fluid chemistry. Fluid samples were collected from felsic-hosted hydrothermal vent fields located on Pual Ridge (PACMANUS and Northeast (NE) Pual) near the active New Britain Arc and a basalt-hosted vent field (Vienna Woods) located farther from the arc on the Manus Spreading Center. Vienna Woods fluids were characterized by relatively uniform endmember temperatures (273–285 °C) and major element compositions, low dissolved CO 2 concentrations (4.4 mmol/kg) and high measured pH (4.2–4.9 at 25 °C). Temperatures and compositions were highly variable at PACMANUS/NE Pual and a large, newly discovered vent area (Fenway) was observed to be vigorously venting boiling (358 °C) fluid. All PACMANUS fluids are characterized by negative &#x3B4; D H 2 O "> δDH 2 O values, in contrast to positive values at Vienna Woods, suggesting substantial magmatic water input to circulating fluids at Pual Ridge. Low measured pH (25 °C) values (∼2.6–2.7), high endmember CO 2 (up to 274 mmol/kg) and negative &#x3B4; 34 S H 2 S "> δ 34 S H2 S values (down to −2.7‰) in some vent fluids are also consistent with degassing of acid-volatile species from evolved magma. Dissolved CO 2 at PACMANUS is more enriched in 13 C (−4.1‰ to −2.3‰) than Vienna Woods (−5.2‰ to −5.7‰), suggesting a contribution of slab-derived carbon. The mobile elements (e.g. Li, K, Rb, Cs and B) are also greatly enriched in PACMANUS fluids reflecting increased abundances in the crust there relative to the Manus Spreading Center. Variations in alkali and dissolved gas abundances with Cl at PACMANUS and NE Pual suggest that phase separation has affected fluid chemistry despite the low temperatures of many vents. In further contrast to Vienna Woods, substantial modification of PACMANUS/NE Pual fluids has taken place as a result of seawater ingress into the upflow zone. Consistently high measured Mg concentrations as well as trends of increasingly non-conservative SO 4 behavior, decreasing endmember Ca/Cl and Sr/Cl ratios with increased Mg indicate extensive subsurface anhydrite deposition is occurring as a result of subsurface seawater entrainment. Decreased pH and endmember Fe/Mn ratios in higher Mg fluids indicate that the associated mixing/cooling gives rise to sulfide deposition and secondary acidity production. Several low temperature (⩽80 °C) fluids at PACMANUS/NE Pual also show evidence for anhydrite dissolution and water–rock interaction (fixation of B) subsequent to seawater entrainment. Hence, the evolution of fluid compositions at Pual Ridge reflects the cumulative effects of water/rock interaction, admixing and reaction of fluids exsolved from silicic magma, phase separation/segregation and seawater ingress into upflow zones.

Manus Basin↗

Seasonal dripwater Mg/Ca and Sr/Ca variations driven by cave ventilation: Implications for and modeling of speleothem paleoclimate records

A 4-year study in a central Texas cave quantifies multiple mechanisms that control dripwater composition and how these mechanisms vary at different drip sites. We monitored cave-air compositions, in situ calcite growth, dripwater composition and drip rate every 4–6 weeks. Three groups of drip sites are delineated (Groups 1–3) based on geochemical variations in dripwater composition. Quantitative modeling of mineral-solution reactions within the host carbonate rock and cave environments is used to identify mechanisms that can account for variations in dripwater compositions. The covariation of Mg/Ca (and Sr/Ca) and Sr isotopes is key in delineating whether Mg/Ca and Sr/Ca variations are dictated by water–rock interaction (i.e., calcite or dolomite recrystallization) or prior calcite precipitation (PCP). Group 1 dripwater compositions reflects a narrow range of the extent of water–rock interaction followed by varying amounts of prior calcite precipitation (PCP). Group 2 dripwater compositions are controlled by varying amounts of water–rock interaction with little to no PCP influence. Group 3 dripwater compositions are dictated by variable extents of both water–rock interaction and PCP. Group 1 drip sites show seasonal variations in dripwater Mg/Ca and Sr/Ca, whereas the other drip sites do not. In contrast to the findings of most previous dripwater Mg/Ca–Sr/Ca studies, these seasonal variations (at Group 1 drip sites) are independent of changes in water flux (i.e., rainfall and/or drip rate), and instead significantly correlate with changes in cave-air CO 2 concentrations. These results are consistent with lower cave-air CO 2 , related to cool season ventilation of the cave atmosphere, enhancing calcite precipitation and leading to dripwater geochemical evolution via PCP. Group 1 dripwater Mg/Ca and Sr/Ca seasonality and evidence for PCP as a mechanism that can account for that seasonality, have two implications for many other regions where seasonal ventilation of caves is likely: (1) speleothem trace-element records may provide seasonal signals, and (2) such records may be biased toward recording climate conditions during the season when calcite is depositing. Additionally, we use our results to construct a forward model that illustrates the types of speleothem Mg/Ca and Sr/Ca variations that would result from varying controls on dripwater compositions. The model provides a basis for interpreting paleo-dripwater controls from high frequency Mg/Ca and Sr/Ca variations for speleothems from caves at which long term monitoring studies are not feasible.

Geochimica et Cosmochimica Acta↗

Uranium isotopes and dissolved organic carbon in loess permafrost: Modeling the age of ancient ice

The residence time of ice in permafrost is an indicator of past climate history, and of the resilience and vulnerability of high-latitude ecosystems to global change. Development of geochemical indicators of ground-ice residence times in permafrost will advance understanding of the circumstances and evidence of permafrost formation, preservation, and thaw in response to climate warming and other disturbance. We used uranium isotopes to evaluate the residence time of segregated ground ice from ice-rich loess permafrost cores in central Alaska. Activity ratios of 234 U vs. 238 U ( 234 U/ 238 U) in water from thawed core sections ranged between 1.163 and 1.904 due to contact of ice and associated liquid water with mineral surfaces over time. Measured ( 234 U/ 238 U) values in ground ice showed an overall increase with depth in a series of five neighboring cores up to 21 m deep. This is consistent with increasing residence time of ice with depth as a result of accumulation of loess over time, as well as characteristic ice morphologies, high segregated ice content, and wedge ice, all of which support an interpretation of syngenetic permafrost formation associated with loess deposition. At the same time, stratigraphic evidence indicates some past sediment redistribution and possibly shallow thaw among cores, with local mixing of aged thaw waters. Using measures of surface area and a leaching experiment to determine U distribution, a geometric model of ( 234 U/ 238 U) evolution suggests mean ages of up to &sim;200 ky BP in the deepest core, with estimated uncertainties of up to an order of magnitude. Evidence of secondary coatings on loess grains with elevated ( 234 U/ 238 U) values and U concentrations suggests that refinement of the geometric model to account for weathering processes is needed to reduce uncertainty. We suggest that in this area of deep ice-rich loess permafrost, ice bodies have been preserved from the last glacial period (10&ndash;100 ky BP), despite subsequent fluctuations in climate, fire disturbance and vegetation. Radiocarbon ( 14 C) analysis of dissolved organic carbon (DOC) in thaw waters supports ages greater than &sim;40 ky BP below 10 m. DOC concentrations in thaw waters increased with depth to maxima of >1000 ppm, despite little change in ice content or cryostructures. These relations suggest time-dependent production of old DOC that will be released upon permafrost thaw at a rate that is mediated by sediment transport, among other factors.

Alaska↗

Global patterns and environmental controls of perchlorate and nitrate co-occurrence in arid and semi-arid environments

Natural perchlorate (ClO 4 − ) is of increasing interest due to its wide-spread occurrence on Earth and Mars, yet little information exists on the relative abundance of ClO 4 − compared to other major anions, its stability, or long-term variations in production that may impact the observed distributions. Our objectives were to evaluate the occurrence and fate of ClO 4 − in groundwater and soils/caliche in arid and semi-arid environments (southwestern United States, southern Africa, United Arab Emirates, China, Antarctica, and Chile) and the relationship of ClO 4 − to the more well-studied atmospherically deposited anions NO 3 − and Cl − as a means to understand the prevalent processes that affect the accumulation of these species over various time scales. ClO 4 − is globally distributed in soil and groundwater in arid and semi-arid regions on Earth at concentrations ranging from 10 −1 to 10 6 μg/kg. Generally, the ClO 4 − concentration in these regions increases with aridity index, but also depends on the duration of arid conditions. In many arid and semi-arid areas, NO 3 − and ClO 4 − co-occur at molar ratios (NO 3 − /ClO 4 − ) that vary between ∼10 4 and 10 5 . We hypothesize that atmospheric deposition ratios are largely preserved in hyper-arid areas that support little or no biological activity (e.g. plants or bacteria), but can be altered in areas with more active biological processes including N 2 fixation, N mineralization, nitrification, denitrification, and microbial ClO 4 − reduction, as indicated in part by NO 3 − isotope data. In contrast, much larger ranges of Cl − /ClO 4 − and Cl − /NO 3 − ratios indicate Cl − varies independently from both ClO 4 − and NO 3 − . The general lack of correlation between Cl − and ClO 4 − or NO 3 − implies that Cl − is not a good indicator of co-deposition and should be used with care when interpreting oxyanion cycling in arid systems. The Atacama Desert appears to be unique compared to all other terrestrial locations having a NO 3 − /ClO 4 − molar ratio ∼10 3 . The relative enrichment in ClO 4 − compared to Cl − or NO 3 − and unique isotopic composition of Atacama ClO 4 − may reflect either additional in-situ production mechanism(s) or higher relative atmospheric production rates in that specific region or in the geological past. Elevated concentrations of ClO 4 − reported on the surface of Mars, and its enrichment with respect to Cl − and NO 3 − , could reveal important clues regarding the climatic, hydrologic, and potentially biologic evolution of that planet. Given the highly conserved ratio of NO 3 − /ClO 4 − in non-biologically active areas on Earth, it may be possible to use alterations of this ratio as a biomarker on Mars and for interpreting major anion cycles and processes on both Mars and Earth, particularly with respect to the less-conserved NO 3 − pool terrestrially.

Geochimica et Cosmochimica Acta↗

Stable C, O and clumped isotope systematics and 14C geochronology of carbonates from the Quaternary Chewaucan closed-basin lake system, Great Basin, USA: Implications for paleoenvironmental reconstructions using carbonates

Isotopic compositions of lacustrine carbonates are commonly used for dating and paleoenvironmental reconstructions. Here we use carbonate δ13C and δ18O, clumped (Δ47), and 14C compositions to better understand the carbonate isotope system in closed-basin lakes and trace the paleohydrologic and temperature evolution in the Chewaucan closed-basin lake system, northern Great Basin, USA, over the Last Glacial/Holocene transition. We focus on shorezone tufas to establish that they form in isotopic equilibrium with lake water and DIC, they can be dated reliably using 14C, and their clumped isotope composition can be used to reconstruct past lake temperature. Calculations of the DIC budget and reservoir age for the lake indicate residence time is short, and dominated by exchange with atmospheric CO2 at all past lake levels. Modern lake DIC and shorezone tufas yield δ13C and 14C values consistent with isotopic equilibrium with recent fossil fuel and bomb-influenced atmospheric CO2, supporting these calculations. δ13C values of fossil tufas are also consistent with isotopic equilibrium with pre-industrial atmospheric CO2 at all shoreline elevations. This indicates that the 14C reservoir effect for this material is negligible. Clumped isotope (Δ47) results indicate shorezone tufas record mean annual lake temperature. Modern (average 13 ± 2 °C) and 18 ka BP-age tufas (average 6 ± 2 °C) have significantly different temperatures consistent with mean annual temperature lowering of 7 ± 3 °C (1 SE) under full glacial conditions. For shorezone tufas and other lake carbonates, including spring mounds, mollusk shells, and ostracod tests, overall δ13C and δ18O values co-vary according to the relative contribution of spring and lacustrine end member DIC and water compositions in the drainage system, but specific isotope values depend strongly upon sample context and are not well correlated with past lake depth. This contrasts with the interpretation that carbonate isotopes in closed-basin lake systems reflect changes in DIC and water budgets connected to higher or lower lake volumes. Instead, a small overlapping range of isotope compositions characterize multiple lake levels, so that none can be identified uniquely by isotope composition alone. Relative to other lake carbonates, δ13C and δ18O values for ostracods in Ana River Canyon deposits are very strongly influenced by Ana River water, suggesting low lake level and volume characterized Summer Lake for most of the past 100,000 years. Coupled with sedimentologic observations, the Ana River deposits thus suggest dry conditions like today are close to the mean climate state in the northern Great Basin. By contrast, basin-integrating highstands such as that dating to ∼14 ka BP, during the last glacial termination, are hydrologically unique and short-lived. Overall, our results indicate carbonate isotope records must account for the specific geochemical and hydrologic characteristics of lake system in order to provide robust paleoenvironmental reconstructions.

Great Basin↗