Geology Reports⌕ Search

SEARCH · Geology Reports

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 541 records · Page 30Linked to original sources
USGS2008

Infiltration from an impoundment for coal‐bed natural gas, Powder River Basin, Wyoming: Evolution of water and sediment chemistry

Development of coal‐bed natural gas (CBNG) in the Powder River Basin, Wyoming, has increased substantially in recent years. Among environmental concerns associated with this development is the fate of groundwater removed with the gas. A preferred water‐management option is storage in surface impoundments. As of January 2007, permits for more than 4000 impoundments had been issued within Wyoming. A study was conducted on changes in water and sediment chemistry as water from an impoundment infiltrated the subsurface. Sediment cores were collected prior to operation of the impoundment and after its closure and reclamation. Suction lysimeters were used to collect water samples from beneath the impoundment. Large amounts of chloride (12,300 kg) and nitrate (13,500 kg as N), most of which accumulated naturally in the sediments over thousands of years, were released into groundwater by infiltrating water. Nitrate was more readily flushed from the sediments than chloride. If sediments at other impoundment locations contain similar amounts of chloride and nitrate, impoundments already permitted could release over 48 × 10 6 kg of chloride and 52 × 10 6 kg of nitrate into groundwater in the basin. A solute plume with total dissolved solid (TDS) concentrations at times exceeding 100,000 mg/L was created in the subsurface. TDS concentrations in the plume were substantially greater than those in the CBNG water (about 2300 mg/L) and in the ambient shallow groundwater (about 8000 mg/L). Sulfate, sodium, and magnesium are the dominant ions in the plume. The elevated concentrations are attributed to cation‐exchange‐enhanced gypsum dissolution. As gypsum dissolves, calcium goes into solution and is exchanged for sodium and magnesium on clays. Removal of calcium from solution allows further gypsum dissolution.

Water Resources Research↗
USGS2011

A behavior-oriented dynamic model for sandbar migration and 2DH evolution

A nonlinear model is developed to study the time‐dependent relationship between the alongshore variability of a sandbar, a ( t ), and alongshore‐averaged sandbar position, x c ( t ). Sediment transport equations are derived from energetics‐based formulations. A link between this continuous physical representation and a parametric form describing the migration of sandbars of constant shape is established through a simple transformation of variables. The model is driven by offshore wave conditions. The parametric equations are dynamically coupled such that changes in one term (i.e., x c ) drive changes in the other (i.e., a ( t )). The model is tested on 566 days of data from Palm Beach, New South Wales, Australia. Using weighted nonlinear least squares to estimate best fit model coefficients, the model explained 49% and 41% of the variance in measured x c and a ( t ), respectively. Comparisons against a 1‐D horizontal (1DH) version of the model showed significant improvements when the 2DH terms were included (1DH and 2DH Brier skill scores were −0.12 and 0.42, respectively). Onshore bar migration was not predicted in the 1DH model, while the 2DH model correctly predicted onshore migration in the presence of 2DH morphology and allowed the bar to remain closer to shore for a given amount of breaking, providing an important hysteresis to the system. The model is consistent with observations that active bar migration occurs under breaking waves with onshore migration occurring at timescales of days to weeks and increasing 2DH morphology, while offshore migration occurs rapidly under high waves and coincides with a reduction in 2DH morphology.

Journal of Geophysical Research C: Oceans↗
USGS2011

Lunar mare deposits associated with the Orientale impact basin: New insights into mineralogy, history, mode of emplacement, and relation to Orientale Basin evolution from Moon Mineralogy Mapper (M3) data from Chandrayaan-1

Moon Mineralogy Mapper (M3) image and spectral reflectance data are combined to analyze mare basalt units in and adjacent to the Orientale multiring impact basin. Models are assessed for the relationships between basin formation and mare basalt emplacement. Mare basalt emplacement on the western nearside limb began prior to the Orientale event as evidenced by the presence of cryptomaria. The earliest post-Orientale-event mare basalt emplacement occurred in the center of the basin (Mare Orientale) and postdated the formation of the Orientale Basin by about 60-100 Ma. Over the next several hundred million years, basalt patches were emplaced first along the base of the Outer Rook ring (Lacus Veris) and then along the base of the Cordillera ring (Lacus Autumni), with some overlap in ages. The latest basalt patches are as young as some of the youngest basalt deposits on the lunar nearside. M3 data show several previously undetected mare patches on the southwestern margins of the basin interior. Regardless, the previously documented increase in mare abundance from the southwest toward the northeast is still prominent. We attribute this to crustal and lithospheric trends moving from the farside to the nearside, with correspondingly shallower density and thermal barriers to basaltic magma ascent and eruption toward the nearside. The wide range of model ages for Orientale mare deposits (3.70-1.66 Ga) mirrors the range of nearside mare ages, indicating that the small amount of mare fill in Orientale is not due to early cessation of mare emplacement but rather to limited volumes of extrusion for each phase during the entire period of nearside mare basalt volcanism. This suggests that nearside and farside source regions may be similar but that other factors, such as thermal and crustal thickness barriers to magma ascent and eruption, may be determining the abundance of surface deposits on the limbs and farside. The sequence, timing, and elevation of mare basalt deposits suggest that regional basin-related stresses exerted control on their distribution. Our analysis clearly shows that Orientale serves as an excellent example of the early stages of the filling of impact basins with mare basalt. Copyright ?? 2011 by the American Geophysical Union.

Journal of Geophysical Research E: Planets↗
USGS2012

Generation and evolution of hydrothermal fluids at Yellowstone: Insights from the Heart Lake Geyser Basin

We sampled fumaroles and hot springs from the Heart Lake Geyser Basin (HLGB), measured water and gas discharge, and estimated heat and mass flux from this geothermal area in 2009. The combined data set reveals that diverse fluids share an origin by mixing of deep solute-rich parent water with dilute heated meteoric water, accompanied by subsequent boiling. A variety of chemical and isotopic geothermometers are consistent with a parent water that equilibrates with rocks at 205°C ± 10°C and then undergoes 21% ± 2% adiabatic boiling. Measured diffuse CO 2 flux and fumarole compositions are consistent with an initial dissolved CO 2 concentration of 21 ± 7 mmol upon arrival at the caldera boundary and prior to southeast flow, boiling, and discharge along the Witch Creek drainage. The calculated advective flow from the basin is 78 ± 16 L s −1 of parent thermal water, corresponding to 68 ± 14 MW, or –1% of the estimated thermal flux from Yellowstone. Helium and carbon isotopes reveal minor addition of locally derived crustal, biogenic, and meteoric gases as this fluid boils and degasses, reducing the He isotope ratio (Rc/Ra) from 2.91 to 1.09. The HLGB is one of the few thermal areas at Yellowstone that approaches a closed system, where a series of progressively boiled waters can be sampled along with related steam and noncondensable gas. At other Yellowstone locations, steam and gas are found without associated neutral Cl waters (e.g., Hot Spring Basin) or Cl-rich waters emerge without significant associated steam and gas (Upper Geyser Basin).

Wyoming↗
USGS2018

Long-term evolution of sand transport through a river network: Relative influences of a dam versus natural changes in grain size from sand waves

Temporal and spatial nonuniformity in supplies of water and sand in a river network leads to sand transport that is in local disequilibrium with the upstream sand supply. In such river networks, sand is transported downstream as elongating waves in which coupled changes in grain size and transport occur. Depending on the magnitude of each sand‐supplying event and the interval between such events, changes in bed‐sand grain size associated with sand‐wave passage may more strongly regulate sand transport than do changes in water discharge. When sand transport is controlled more by episodic resupply of sand than by discharge, upstream dam construction may exacerbate or mitigate sand‐transport disequilibria, thus leading to complicated and difficult‐to‐predict patterns of deposition and erosion. We analyzed all historical sediment‐transport data and embarked on a 4‐year program of continuous sediment‐transport measurements to describe disequilibrium sand transport in a river network. Results indicate that sand transport in long river segments can evolve over ≥50‐year timescales following rare large sand‐supplying events. These natural changes in sand transport in distal downstream river segments can be larger than those caused by an upstream dam. Because there is no way to know a priori whether sand transport in a river has changed in response to changes in the upstream sand supply, contemporary continuous measurements of sand transport are required for accurate sand loads and budgeting. Analysis of only historical sediment‐transport measurements, as is common in the literature, may lead to incorrect conclusions with respect to current or future sediment‐transport conditions.

Journal of Geophysical Research: Earth Surface↗
USGS2018

Geomorphic evolution of a gravel‐bed river under sediment‐starved vs. sediment‐rich conditions: River response to the world's largest dam removal

Understanding river response to sediment pulses is a fundamental problem in geomorphic process studies, with myriad implications for river management. However, because large sediment pulses are rare and usually unanticipated, they are seldom studied at field scale. We examine fluvial response to a massive (~20 Mt) sediment pulse released by the largest dam removal globally, on the Elwha River, Washington, United States, in an 11‐year before‐after/control‐impact study of channel morphology and grain size. We test the hypothesis that for a given flow magnitude, greater geomorphic change occurs under sediment‐rich conditions than under sediment‐starved conditions. Channel response to flow forcing was significantly different during the sediment‐pulse peak, 1–2 years after dam removal began, than earlier or later. During peak sediment supply our hypothesis was supported; major geomorphic change occurred under low flows and unit stream power ≤60 W/m 2 . However, by 4–6 years after dam removal began, rates of geomorphic change and sensitivity to stream power had decreased substantially such that our hypothesis was no longer unequivocally supported. These findings are consistent with a two‐phase conceptual model of dam‐removal response, involving a transport‐limited state followed by a more supply‐limited state. From comparisons with other dam removals and natural sediment pulses, we infer that the longevity of sediment‐pulse signals in gravel‐bed rivers depends upon gradient, river discharge, valley morphology, and sediment grain size. Stream power associated with substantial geomorphic change varies with sediment supply, such that assigning a general threshold stream power to gravel‐bed rivers may be untenable.

Washington↗
USGS2021

Basalt geochemistry and mantle flow during early backarc basin evolution: Havre Trough and Kermadec Arc, southwest Pacific

The Havre Trough (HT) backarc basin in the southwest Pacific is in the rifting stage of development. We distinguish five types of basalt there based on their amount and kind of slab component: backarc basalts (BAB) with little or no slab component, modified BAB with slight amounts, reararc (RA) with more, remnants of the preexisting arc (Colville Ridge horsts), and arc front volcanoes within the HT. Previous subarc mantle is quickly removed and replaced by more fertile mantle with less slab component. The ambient mantle is “Pacific” isotopically, and more enriched in Nb/Yb and Nd and Hf isotope ratios north of the Central Kermadec Discontinuity at 32°S than to the south. The contrast may reflect inheritance in the south of mantle that was depleted during spreading that formed the southern South Fiji Basin and a higher degree of melting because of a wetter slab-derived flux. The slab component also differs along strike, more like a dry melt in the north and a supercritical fluid in the south. The mass fraction of slab component increases southward in the backarc as well as the arc front. RA volcanoes have the most slab component (1%–2%) and form indistinct ridges at high angles to, and <50 km behind, frontal volcanoes. Backarc basalts have less and occur throughout the basin. Slab components are distributed further into the backarc, and more irregularly, during the rifting than spreading stage of backarc basin development. The rifting stage is disorganized geochemically as well as spatially.

Geochemistry Geophysics Geosystems↗
USGS2020

The Cenozoic evolution of crustal shortening and left‐lateral shear in the central East Kunlun Shan: Implications for the uplift history of the Tibetan Plateau

The timing of crustal shortening and strike‐slip faulting along the East Kunlun Shan provides insight into the history of surface uplift and may constrain the time at which the Tibetan Plateau reached high elevations. We investigate a series of extensional basins and restraining bends along the Xidatan strand of the Kunlun strike‐slip fault, which provide an ideal setting to unravel the tectonic history of the northern plateau margin. We present new apatite (U‐Th)/He, apatite fission track, and zircon (U‐Th)/He ages and QTQt thermal modeling, 40 Ar/ 39 Ar fault gouge dating, and structural mapping from the central East Kunlun Shan. Our data suggest that the East Kunlun Shan experienced slow to negligible exhumation until late Cretaceous time, followed by an increase in rate by 65–50 Ma. Along with a ~47 Ma fault gouge age, we posit that the Paleocene–early Eocene was a time of crustal shortening along the northern plateau. Rapid exhumation along transpressional portions of the Xidatan fault initiated by 23–20 Ma, which we interpret as the local onset of strike‐slip faulting. An early Miocene transition from north‐south crustal shortening to left‐lateral shear along the East Kunlun Shan, the onset of normal and strike‐slip faulting in central and southern Tibet by 18 Ma, and lower crustal flow in eastern Tibet by 13 Ma suggest the establishment of orogen‐wide east‐west oriented extension and extrusion by the middle Miocene. The plateau‐wide shift in stress accommodation implies that high gravitational potential energy, and likely high elevation, was attained by the middle Miocene.

Tectonics↗
USGS2021

An eddy-resolving numerical model to study turbulent flow, sediment and bed evolution using detached eddy simulation in a lateral separation zone at the field-scale

Turbulence-resolving simulations elucidate key elements of fluid dynamics and sediment transport in fluvial environments. This research presents a feasible strategy for applying state-of-the-art computational fluid mechanics to the study of sediment transport and morphodynamic processes in lateral separation zones, which are common features in canyon rivers where massive lateral flow separation causes large-scale turbulence that controls sediment erosion and deposition. An eddy-resolving model was developed and tested at the field-scale, coupling a viscous flow and sediment transport solver using Detached Eddy Simulation techniques. A morphodynamic model was applied to the viscous flow/sediment solver to calculate erosion and deposition. A simulation of turbulence was performed at the grid resolution for a straight channel to determine the relative contributions of modeled and resolved diffusivity. The time-dependent, energetically important, correlative, non-stationary signals of the simulated quantities were captured at the lateral separation zone. Strong periodic signals featured by high amplitude were found at the separation zone, while low frequency pulsations were observed at the reattachment zone of the lateral separation zone. Interactions between the eddies and the loose bed boundaries resulted in erosion of sediment at the main channel followed by deposition at the primary eddy and eddy bars. tions elucidate key elements of fluid dynamics and sediment transport in fluvial environments. This research presents a feasible strategy for applying state-of-the-art computational fluid mechanics to the study of sediment transport and morphodynamic processes in lateral separation zones, which are common features in canyon rivers where massive lateral flow separation causes large-scale turbulence that controls sediment erosion and deposition. An eddy-resolving model was developed and tested at the field scale, coupling viscous flow and sediment transport solver using Detached Eddy Simulation (DES) techniques. A morphodynamic model was applied to the viscous flow/sediment solver to calculate erosion and deposition. A simulation of turbulence was performed at the grid resolution for a straight channel to determine the relative contributions of modeled and resolved diffusivity. The time-dependent, energetically important, correlative non-stationary signals of the simulated quantities were captured at the lateral separation zone. Strong periodic signals featured by high amplitude were found at the separation zone, while low frequency pulsations were observed at the reattachment zone of the lateral separation zone. Interactions between the eddies and the loose bed boundaries resulted in massive erosion of sediment at the main channel followed by deposition at the primary eddy and eddy bars.

Arizona↗
USGS2022

Multi-decadal simulation of marsh topography evolution under sea level rise and episodic sediment loads

Coastal marsh within Mediterranean climate zones is exposed to episodic watershed runoff and sediment loads that occur during storm events. Simulating future marsh accretion under sea level rise calls for attention to: (a) physical processes acting over the time scale of storm events and (b) biophysical processes acting over time scales longer than storm events. Using the upper Newport Bay in Southern California as a case study, we examine the influence of event-scale processes on simulated change in marsh topography by comparing: (a) a biophysical model that integrates with an annual time step and neglects event-scale processes (BP-Annual), (b) a physical model that resolves event-scale processes but neglects biophysical interactions (P-Event), and (c) a biophysical model that resolves event-scale physical processes and biophysical processes at annual and longer time scales (BP-Event). A calibrated BP-Event model shows that large (>20-year return period) episodic storm events are major drivers of marsh accretion, depositing up to 30 cm of sediment in one event. Greater deposition is predicted near fluvial sources and tidal channels and less on marshes further from fluvial sources and tidal channels. In contrast, the BP-Annual model poorly resolves spatial structure in marsh accretion as a consequence of neglecting event-scale processes. Furthermore, the P-Event model significantly overestimates marsh accretion as a consequence of neglecting marsh surface compaction driven by annual scale biophysical processes. Differences between BP-Event and BP-Annual models translate up to 20 cm per century in marsh surface elevation.

California↗
USGS2023

Future marsh evolution due to tidal changes induced by human adaptation to sea level rise

With sea level rise threatening coastal development, decision-makers are beginning to act by modifying shorelines. Previous research has shown that hardening or softening shorelines may change the tidal range under future sea level rise. Tidal range can also be changed by natural factors. Coastal marshes, which humans increasingly depend on for shoreline protection, are ecologically sensitive to tidal range. Therefore, it is critical to examine how changes in tidal range could influence marsh processes. A marsh accretion model was used to investigate the ecological response of a San Francisco Bay, California, USA marsh to multiple tidal range scenarios and sea level rise from 2010 to 2100. The scenarios include a baseline scenario with no shoreline modifications in the estuary, a shoreline hardening scenario that amplifies the tidal range, and 14 tidal range scenarios as a sensitivity analysis that span tidal amplification and reduction of the baseline scenario. The modeling results expose key tradeoffs to consider when planning for sea level rise. Compared to the baseline, the hardening scenario shows minor differences. However, further tidal amplification prolongs marsh survival but decreases Sarcocornia pacifica cover, an important species for certain threatened wildlife and an effective attenuator of wave energy. Conversely, tidal reduction precipitates marsh drowning but shows gains in Sarcocornia pacifica cover. These mixed impacts of tidal amplification and reduction shown by the model indicate potential tradeoffs in relation to marsh survival, habitat characteristics, and shoreline protection. This study suggests the need for a cross-sectoral, regional approach to sea level rise adaptation.

California↗
USGS2024

Evolution of a lake margin recorded in the Sutton Island member of the Murray formation, Gale crater, Mars

This study uses data from the Mars Science Laboratory Curiosity rover to document the facies of the Sutton Island member of the Murray formation, interpret paleoenvironments, and establish key stratigraphic transitions at Gale crater. Two facies associations were identified: Facies Association 1 (FA1) and Facies Association 2 (FA2). Individual facies in FA1 include planar-laminated mudstone with minor intervals of planar sandstone, ripple cross-laminated sandstone, cross-stratified sandstone, and alternating laminated sandstone and mudstone. Meter-thick packages of planar-laminated mudstone in FA1 are interpreted to represent deposition in low-energy ponded environments along the lake margin. Straight- and curve-crested ripple cross-laminated facies are interpreted to represent current-influenced deposition. Cross-stratified sandstone facies consist of dm-thick sets that represent deposition in distal channels. Intercalated mm-scale mudstone and sandstone laminae represent waning flow conditions and possible channel abandonment. Facies in FA1 collectively represent deposition in a distal delta plain. FA2 is comprised of planar-laminated mudstone with minor sandstone and is interpreted to represent deposition in a lacustrine-basin setting by suspension settling linked to density flows. FA1 transitions upward into FA2, defining a rapid transgression substantial enough to facilitate the deposition of distal lake facies above delta plain facies. The abrupt transition from FA2 back to FA1 deltaic deposits is suggestive of forced regression. Facies in FA1 and FA2 are consistent with the prevalence of aqueous environments recorded in other Murray formation members and extend our understanding of the dynamic sedimentary processes that characterized ancient lacustrine systems at Gale crater.

Journal of Geophysical Research - Planets↗
USGS2024

Formation and evolution of the Pacific-North American (San Andreas) plate boundary: Constraints from the crustal architecture of northern California

The northward migration of the Mendocino triple junction (MTJ) drives a fundamental plate boundary transformation from convergence to translation; producing a series of strike-slip faults, that become the San Andreas plate boundary. We find that the 3-D structure of the Pacific plate lithosphere in the vicinity of the MTJ controls the location of San Andreas plate boundary formation. At the time of initiation of the Pacific-North America plate boundary (∼30 Ma), the sequential interaction with the western margin of North America of the Pioneer Fracture Zone, soon followed by the Mendocino Fracture Zone, led to the capture of a small segment of partially subducted Farallon lithosphere by the Pacific plate, termed the Pioneer Fragment (PF). Since that time, the PF has translated with the Pacific Plate along the western margin of North America. Recently developed, high-resolution seismic-tomographic imagery of northern California indicates that (a) the PF is extant, occupying the western half of the slab window, immediately south of the MTJ; (b) the eastern edge of the PF lies beneath the newly forming Maacama fault system, which develops to become the locus for the primary plate boundary structure after approximately 6–10 Ma; and (c) the location of the translating PF adjacent to the asthenosphere of the slab window generates a shear zone within and below the crust that develops into the plate boundary faults. As a result, the San Andreas plate boundary forms interior to the western margin of North America, rather than at its western edge.

California↗
USGS2024

Zircon constraints on the eruptive sequence and magma evolution of rhyolites at South Sister volcano, Oregon

We present 230 Th- 238 U crystallization ages and trace element compositions for zircons spanning the late Pleistocene to Holocene rhyolite eruptive record at South Sister volcano in the central Oregon Cascade Range. Most zircon ages are between 100 and 20 ka, with very few in secular equilibrium (>350 ka). The weighted mean of zircon ages for the two oldest South Sister rhyolites, 31.5 ± 2.1 and 39.1 ± 2.4 ka, are significantly younger than the associated 40 Ar/ 39 Ar ages, 47.4 ± 9.7 and 51.4 ± 9.7 ka. We propose that these 40 Ar/ 39 Ar dates, performed on plagioclase separates, are compromised by a subtle amount of excess Ar and therefore the younger weighted mean zircon ages yield more reliable eruption ages. These results imply that the interval of rhyolite eruption at South Sister during the late Pleistocene was both shorter and more productive than previously thought and that eruption at South Sister initiated after Middle Sister. Compositionally, zircons from the Pleistocene rhyolites are broadly similar and show down-temperature zircon and plagioclase crystallization trends. However, we argue that destabilized amphibole and titanite in a common mush also exert leverage on the Pleistocene zircon trace element compositions. Divergence in the Eu/Eu* ratio between the Pleistocene and Holocene lavas implies chemically distinct magma reservoirs originating from the Pleistocene rhyolite eruptive sequence and the Holocene eruptive sequence. This work suggests a higher flux of rhyolite volcanism than previously thought and characterizes magmatic storage distinctions between the Pleistocene and Holocene rhyolites, aiding in the assessment of future eruptive hazards at South Sister volcano.

Oregon↗
USGS2026

Structural evolution and slip rate variations through time of the Puente-Hills blind-thrust fault beneath Los Angeles: Implications for seismic hazard and folding kinematics

Using seismic reflection profiles, historical well logging data, and luminescence and radiocarbon ages, we determine a Pleistocene-Holocene slip history for the central, Santa Fe Springs segment of the Puente Hills blind-thrust fault (PHT), a major seismogenic fault situated beneath the urbanized Los Angeles metropolitan region. We analyze the geometry of correlative stratigraphic units in the forelimb and backlimb of the overlying growth-fold, the Santa Fe Springs anticline, and determine the uplift of seven age-correlative markers. Uplift measurements are converted to thrust displacements on the underlying PHT using a structural method laid out by Don et al. (2022, https://doi.org/10.1785/0120220048 ) that accounts for the geometry of the fault. These data indicate that deep thrust displacement on the PHT is partially consumed updip in the creation of a hanging-wall fault-bend fold, with forelimb growth strata recording <80% of the slip documented within the backlimb. Chronological data from growth strata yield age constraints for folding and faulting on the underlying PHT, providing a detailed incremental slip history derived from both the forelimb and backlimb folding for seven discrete growth horizons, spanning the past 1.4 million years. The resulting six incremental slip rates demonstrate that fault slip has varied through time from the middle Pleistocene to the Holocene. Moreover, these results reveal synchronous acceleration of both the central, Santa Fe Springs and western, Los Angeles segments of the fault system since late Pleistocene time (after 200 ka) and slip rates of greater than 2 mm/yr on the downdip, backlimb, fault-ramp below the anticline.

California↗
USGS2026

Impacts of transient critical zone evolution on landslide susceptibility and hazard Across the Utuado batholith, Puerto Rico

Landscapes in humid climates often experience rapid chemical and mechanical weathering that can change groundwater fluxes, drainage pathways, and frictional strength via alterations in subsurface mineral composition and matrix porosity. These changes to hillslope critical zone structure and strength should control landslide susceptibility, size, and failure style. Here we test the impact of differential weathering on landslide failure processes by exploiting a natural experiment along the dioritic Utuado batholith in Puerto Rico's Cordillera Central, where a sharp weathering gradient follows an upstream-propagating knickzone into an uplifted shoreline platform that was affected by Hurricane Maria in 2017. Through a combination of landslide mapping, field geophysical and geotechnical surveys, topographic analysis, calibrated stream profile inversion, 3D back-analysis of hillslope strength, and debris-flow inundation assessments, we show that differential weathering induces changes to critical zone structure and consequently landslide susceptibility and hazard. Our mapping reveals that the hurricane-induced landslides preferentially cluster along less weathered hillslopes despite a marked increase in frictional strength. Jointly inverted seismic refraction and electrical resistivity profiles show that the frictionally weak, heavily weathered hillslopes are highly porous and dry, which we posit acts to stabilize slopes through limiting deeper hydrologic connectivity that can produce shallow water tables with destabilizing lateral seepage forces during storms. Debris-flow inundation mapping and analysis reveals that major contributors to debris-flow growth, channel steepness and landslide coalescence, scale linearly with hillslope friction angle. This leads to higher-mobility debris flows along knickzone reaches, which disproportionately affects communities at the base of the plateau.

Puerto Rico↗
USGS2012

Evolution of the caldera‐forming eruption at Crater Lake, Oregon, indicated by component analysis of lithic fragments

Crater Lake caldera (8 × 10 km), formed 6845 years B. P. ( 14 C age) during the climactic eruption of the volcanic edifice known as Mount Mazama, is intermediate in size between small calderas associated with central vent eruptions and large calderas that have ring fracture vent systems. Our quantitative study of lithic fragments in the ejecta confirms the existing model of changes in vent configuration during the climactic eruption of Mount Mazama. Initial activity was from a single vent that produced a rhyodacite pumice fall from a Plinian column. Altered preexisting volcanic rocks are the predominant lithic type in the Plinian deposit, and their extensive hydrothermal alteration is considered as evidence of their relatively deep origin. The Wineglass Welded Tuff lies atop the Plinian deposit and contains a higher proportion of fresh volcanic rocks, suggesting enlargement of the single vent by slumping of its walls. This same vent enlargement caused the Plinian eruption column to collapse and feed valley‐hugging pyroclastic flows that deposited the Wineglass Welded Tuff. When enough material was erupted from the shallow magma chamber that its roof was no longer adequately supported, Mount Mazama collapsed to form the caldera, while highly energetic pyroclastic flows produced the climactic ignimbrite. A lag breccia that represents the proximal facies of the compositionally zoned climactic ignimbrite lies atop the Wineglass Welded Tuff and contains predominantly altered volcanic rocks of deeper origin, accompanied by minor granitoids from the magma chamber walls. Azimuthal differences in lithic component proportions in the lag breccia correlate well with the geology of the caldera walls, indicating that the climactic ignimbrite was ejected by multiple vents along a ring fracture system. Systematic lithic component changes within the lag breccia suggest different quarrying levels that reflect waxing and waning of the discharge rate during the caldera collapse phase of the climactic eruption. Our lithic component analysis demonstrates that calderas that may be too small to experience structural resurgence, such as Crater Lake, nevertheless may form by syneruptive subsidence along ring fractures.

Oregon↗
USGS2010

Chemical evolution of groundwater near a sinkhole lake, northern Florida: 1. Flow patterns, age of groundwater, and influence of lakewater leakage

Leakage from sinkhole lakes significantly influences recharge to the Upper Floridan aquifer in poorly confined sediments in northern Florida. Environmental isotopes (oxygen 18, deuterium, and tritium), chlorofluorocarbons (CFCs: CFC-11, CCl 3 F; CFC-12, CCl 2 F 2 ; and CFC-113, C 2 Cl 3 F 3 ), and solute tracers were used to investigate groundwater flow patterns near Lake Barco, a seepage lake in a mantled karst setting in northern Florida. Stable isotope data indicated that the groundwater downgradient from the lake contained 11–67% lake water leakage, with a limit of detection of lake water in groundwater of 4.3%. The mixing fractions of lake water leakage, which passed through organic-rich sediments in the lake bottom, were directly proportional to the observed methane concentrations and increased with depth in the groundwater flow system. In aerobic groundwater upgradient from Lake Barco, CFC-modeled recharge dates ranged from 1987 near the water table to the mid 1970s for water collected at a depth of 30 m below the water table. CFC-modeled recharge dates (based on CFC-12) for anaerobic groundwater downgradient from the lake ranged from the late 1950s to the mid 1970s and were consistent with tritium data. CFC-modeled recharge dates based on CFC-11 indicated preferential microbial degradation in anoxic waters. Vertical hydraulic conductivities, calculated using CFC-12 modeled recharge dates and Darcy's law, were 0.17, 0.033, and 0.019 m/d for the surficial aquifer, intermediate confining unit, and lake sediments, respectively. These conductivities agreed closely with those used in the calibration of a three-dimensional groundwater flow model for transient and steady state flow conditions.

Florida↗