SEARCH · USGS Science

Results for “Earth Interactions”

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 469 records · Page 26Linked to original sources

USGS2014

Predictions of barrier island berm evolution in a time-varying storm climatology

Low-lying barrier islands are ubiquitous features of the world's coastlines, and the processes responsible for their formation, maintenance, and destruction are related to the evolution of smaller, superimposed features including sand dunes, beach berms, and sandbars. The barrier island and its superimposed features interact with oceanographic forces (e.g., overwash) and exchange sediment with each other and other parts of the barrier island system. These interactions are modulated by changes in storminess. An opportunity to study these interactions resulted from the placement and subsequent evolution of a 2 m high sand berm constructed along the northern Chandeleur Islands, LA. We show that observed berm length evolution is well predicted by a model that was fit to the observations by estimating two parameters describing the rate of berm length change. The model evaluates the probability and duration of berm overwash to predict episodic berm erosion. A constant berm length change rate is also predicted that persists even when there is no overwash. The analysis is extended to a 16 year time series that includes both intraannual and interannual variability of overwash events. This analysis predicts that as many as 10 or as few as 1 day of overwash conditions would be expected each year. And an increase in berm elevation from 2 m to 3.5 m above mean sea level would reduce the expected frequency of overwash events from 4 to just 0.5 event-days per year. This approach can be applied to understanding barrier island and berm evolution at other locations using past and future storm climatologies.

Louisiana↗

USGS2025

Interacting sea-level rise, sea-ice loss, storm flooding, erosion, and permafrost thaw threaten ecosystems, wildlife, and communities on the Yukon-Kuskokwim Delta

The Yukon-Kuskokwim Delta has the largest intertidal wetland in North America, is a globally critical breeding area for waterbirds, and is home to the largest regional indigenous population in the Arctic. Here, coastal tundra ecosystems, wildlife, and indigenous communities are highly vulnerable to sea-ice loss in the Bering Sea, sea-level rise, storm flooding, erosion, and collapsing ground from permafrost thaw caused by climate warming. These drivers interact in non-linear ways to increase flooding, salinization, and sedimentation, and thus, alter ecosystem trajectories and broader landscape evolution. Rapid changes in these factors over decadal time scales are highly likely to cause transformative shifts in coastal ecosystems across roughly 70% of the outer delta this century. We project saline and brackish ecotypes on the active delta floodplain with frequent sedimentation will maintain dynamic equilibrium with sea-level rise and flooding, slightly brackish ecotypes on the inactive floodplain with infrequent flooding and low sedimentation rates will be vulnerable to increased flooding and likely transition to more saline and brackish ecotypes, and fresh lacustrine and lowland ecotypes on the abandoned floodplain with permafrost plateaus will be vulnerable to thermokarst, salinization and flooding that will shift them toward brackish ecosystems. This will greatly affect bird nesting and foraging habitats, with both winners and losers. Already, some Yup'ik communities are facing relocation of their low-lying villages. The societal challenges and consequences of adapting to these changing landscapes are enormous and will require a huge societal effort.

Alaska↗

USGS2018

Opal-A in glassy pumice, acid alteration, and the 1817 phreatomagmatic eruption at Kawah Ijen (Java), Indonesia

At Kawah Ijen (Indonesia), vigorous SO 2 and HCl degassing sustains a hyperacid lake (pH ~0) and intensely alters the subsurface, producing widespread residual silica and advanced argillic alteration products. In 1817, a VEI 2 phreatomagmatic eruption evacuated the lake, depositing a widespread layer of muddy ash fall, and sending lahars down river drainages. We discovered multiple types of opaline silica in juvenile low-silica dacite pumice and in particles within co-erupted laharic sediments. Most spectacular are opal-replaced phenocrysts of plagioclase and pyroxene adjacent to pristine matrix glass and melt inclusions. Opal-bearing pumice has been found at numerous sites, including where post-eruption infiltration of acid water is unlikely. Through detailed analyses of an initial sampling of 1817 eruption products, we find evidence for multiple origins of opaline materials in pumice and laharic sediments. Evidently, magma encountered acid-altered materials in the subsurface and triggered phreatomagmatic eruptions. Syn-eruptive incorporation of opal-alunite clasts, layered opal, and fragment-filled vesicles of opal and glass, all suggest magma-rock interactions in concert with vesiculation, followed by cooling within minutes. Our experiments at magmatic temperature confirm that the opaline materials would show noticeable degradation in time periods longer than a few tens of minutes. Some glassy laharic sedimentary grains are more andesitic than the main pumice type and may represent older volcanic materials that were altered beneath the lake bottom and were forcefully ejected during the 1817 eruption. A post-eruptive origin remains likely for most of the opal-replaced phenocrysts in pumice. Experiments at 25°C and 100°C reveal that when fresh pumice is bathed in Kawah Ijen hyperacid fluid for 6 weeks, plagioclase is replaced without altering either matrix glass or melt inclusions. Moreover, lack of evidence for high-temperature annealing of the opal suggests that post-eruption alteration of pumice is more likely than pre-eruption envelopment of euhedral opal-replaced phenocrysts in dacitic melt. At Ijen and elsewhere, the ascent of magma into hydrous acid-altered mineral assemblages (e.g., opal, kaolinite, alunite) could induce rapid dehydration of hydrous minerals and amorphous materials, generating considerable steam and contributing to magmatic-hydrothermal and phreatomagmatic explosions.

Kawah Ijen↗

USGS2006

The airborne lava-seawater interaction plume at Kilauea Volcano, Hawai'i

Lava flows into the sea at Kīlauea Volcano, Hawaiʻi, and generates an airborne gas and aerosol plume. Water (H 2 O), hydrogen chloride (HCl), carbon dioxide (CO 2 ), nitrogen dioxide (NO 2 ) and sulphur dioxide (SO 2 ) gases were quantified in the plume in 2004–2005, using Open Path Fourier Transform infra-red Spectroscopy. The molar abundances of these species and thermodynamic modelling are used to discuss their generation. The range in molar HCl / H 2 O confirms that HCl is generated when seawater is boiled dry and magnesium salts are hydrolysed (as proposed by [T.M. Gerlach, J.L. Krumhansl, R.O. Fournier, J. Kjargaard, Acid rain from the heating and evaporation of seawater by molten lava: a new volcanic hazard, EOS (Trans. Am. Geophys. Un.) 70 (1989) 1421–1422]), in contrast to models of Na-metasomatism. Airborne droplets of boiled seawater brine form nucleii for subsequent H 2 O and HCl condensation, which acidifies the droplets and liberates CO 2 gas from bicarbonate and carbonate. NO 2 is derived from the thermal decomposition of nitrates in coastal seawater, which takes place as the lava heats droplets of boiled seawater brine to 350–400 °C. SO 2 is derived from the degassing of subaerial lava flows on the coastal plain. The calculated mass flux of HCl from a moderate-sized ocean entry significantly increases the total HCl emission at Kīlauea (including magmatic sources) and is comparable to industrial HCl emitters in the United States. For larger lava ocean entries, the flux of HCl will cause intense local environmental hazards, such as high localised HCl concentrations and acid rain.

Hawaii↗

USGS1998

Two-way coupling between Vesuvius eruptions and southern Apennine earthquakes, Italy, by elastic stress transfer

During the past 1000 years, eruptions of Vesuvius have often been accompanied by large earthquakes in the Apennines 50–60 km to the northeast. Statistical investigations had shown that earthquakes often preceded eruptions, typically by less than a decade, but did not provide a physical explanation for the correlation. Here, we explore elastic stress interaction between earthquakes and eruptions under the hypothesis that small stress changes can promote events when the Apennine normal faults and the Vesuvius magma body are close to failure. We show that earthquakes can promote eruptions by compressing the magma body at depth and opening suitably oriented near-surface conduits. Voiding the magma body in turns brings these same normal faults closer to Coulomb failure, promoting earthquakes. Such a coupling is strongest if the magma reservoir is a dike oriented normal to the regional extension axis, parallel to the Apennines, and the near-surface conduits and fissures are oriented normal to the Apennines. This preferred orientation suggests that the eruptions issuing from such fissures should be most closely linked in time to Apennine earthquakes. Large Apennine earthquakes since 1400 are calculated to have transferred more stress to Vesuvius than all but the largest eruptions have transferred to Apennine faults, which may explain why earthquakes more commonly lead than follow eruptions. A two-way coupling may thus link earthquakes and Vesuvius eruptions along a 100-km-long set of faults. We test the statistical significance of the earthquake-eruption correlation in the two-way coupling zone, and find a correlation significant at the 95% confidence level.

Journal of Geophysical Research B: Solid Earth↗

USGS2026

A temporal look at the influence of topographic amplification on earthquake-triggered landslides in 3D seismic simulations

Earthquakes are a primary trigger for landslides, often leading to catastrophic consequences. While numerous studies have explored the spatial distribution of earthquake-triggered landslides, understanding the interaction between seismic waves and topography remains a critical challenge. Topographic irregularities can cause seismic wave amplification, altering ground shaking, and can trigger landslides that are challenging for predictive models to anticipate. This study investigates the spatial and temporal evolution of topographically amplified landslides, focusing on coseismic landslides triggered by the Mw 7.5 mainshock of the 2018 Papua New Guinea earthquake and post-seismic landslides associated with its four aftershocks, each exceeding Mw 6.0. We employ low-frequency, three-dimensional numerical ground shaking simulations and data-driven multivariate analyses to examine how landslides evolved from the coseismic to post-seismic periods. Our findings reveal a spatial shift in landslide distribution, in which the mainshock triggered slope failures predominantly on steep hillslopes, and the aftershocks triggered landslides on gentler slopes, often near geologic boundaries. We attribute this transition partly to the earthquake legacy effect of the mainshock, where the mainshock caused weakening of these hillslopes, making them more prone to failure when aftershocks occur. Additionally, the concentration of failures along geologic contacts in the post-seismic phase suggests that site amplification, stemming from contrasts in subsurface materials, exerts a key influence on landslide occurrence. Although not explicitly captured in our current numerical simulations, this mechanism warrants further investigation for more accurate hazard modeling.

JGR Solid Earth↗

USGS2023

Competition between dissolved organic matter and freshwater plankton control methylmercury isotope fractionation during uptake and photochemical demethylation

Isotope fractionation related to photochemical reactions and planktonic uptake at the base of the food web is a major uncertainty in the biological application of mercury (Hg) stable isotopes. In freshwater systems, it is unclear how competitive interactions among methylmercury (MeHg), dissolved organic matter (DOM), and phytoplankton govern the magnitude of mass-dependent and mass-independent fractionation. This study investigated how DOM alters rates of planktonic MeHg uptake and photodegradation and corresponding Hg isotope fractionation in the presence of freshwater phytoplankton species, Raphidocelis subcapitata . Outdoor sunlight exposure experiments utilizing R. subcapitata were performed in the presence of different DOM samples using environmentally relevant ratios of MeHg-DOM thiol groups. The extent of Δ 199 Hg in phytoplankton incubations (2.99‰ St. Louis River HPOA, 1.88‰ Lake Erie HPOA) was lower compared to paired abiotic control experiments (4.29 and 2.86‰, respectively) after ∼30 h of irradiation, resulting from cell shading or other limiting factors reducing the extent of photodemethylation. Although the Δ 199 Hg/Δ 201 Hg ratio was uniform across experiments (∼1.4), Δ 199 Hg/δ 202 Hg slopes varied dramatically (from −0.96 to 15.4) across incubations with R. subcapitata and DOM. In addition, no evidence of Hg isotope fractionation was observed within R. subcapitata cells. This study provides a refined examination of Hg isotope fractionation markers for key processes occurring in the lower food web prior to bioaccumulation, critical for accurately accounting for the photochemical processing of Hg isotopes across a wide spectrum of freshwater systems.

ACS Earth and Space Chemistry↗

USGS2026

Role of volcano-tectonic interactions during early-phase magma-assisted continental rifting

We study volcanic-tectonic interactions during the early phase of rifting in the Natron Rift, a magmatic segment of the East African Rift System. We use the open-source crustal deformation modeling code, PyLith, to model deflation during the 2007–2008 explosive eruptions and inflation during 2016–2023. We compare predictions of surface displacements resulting from magma reservoir inflation from 3D models with and without topography and homogeneous and heterogeneous material properties against Global Navigation Satellite System displacements from 2016 to 2023. We prefer the model with heterogeneous material properties and topography for fitting the observed displacements. Including topography improves the fit to GNSS data and modulates stress changes from shallow magma sources. Our numerical modeling of Coulomb stress changes suggests that deflation of a magma reservoir under the active Ol Doinyo Lengai volcano could trigger oblique normal slip on the Natron Fault with Coulomb stress increases exceeding ∼0.1 MPa on the volcano-adjacent fault segment. Fault slip on the Natron Fault is less likely during inflation of the magma reservoir because normal compressive stress increases on the fault. We also demonstrate that Coulomb stress changes are sensitive to the inclusion of topography in our model, particularly because the depth of the magma reservoir is similar to the amplitude of the topographic relief.

Ol Doinyo Lenga↗

USGS2008

Eruption dynamics at Mount St. Helens imaged from broadband seismic waveforms: Interaction of the shallow magmatic and hydrothermal systems

The current eruption at Mount St. Helens is characterized by dome building and shallow, repetitive, long-period (LP) earthquakes. Waveform cross-correlation reveals remarkable similarity for a majority of the earthquakes over periods of several weeks. Stacked spectra of these events display multiple peaks between 0.5 and 2 Hz that are common to most stations. Lower-amplitude very-long-period (VLP) events commonly accompany the LP events. We model the source mechanisms of LP and VLP events in the 0.5-4 s and 8-40 s bands, respectively, using data recorded in July 2005 with a 19-station temporary broadband network. The source mechanism of the LP events includes: 1) a volumetric component modeled as resonance of a gently NNW-dipping, steam-filled crack located directly beneath the actively extruding part of the new dome and within 100 m of the crater floor and 2) a vertical single force attributed to movement of the overlying dome. The VLP source, which also includes volumetric and single-force components, is 250 m deeper and NNW of the LP source, at the SW edge of the 1980s lava dome. The volumetric component points to the compression and expansion of a shallow, magma-filled sill, which is subparallel to the hydrothermal crack imaged at the LP source, coupled with a smaller component of expansion and compression of a dike. The single-force components are due to mass advection in the magma conduit. The location, geometry and timing of the sources suggest the VLP and LP events are caused by perturbations of a common crack system.

Journal of Geophysical Research B: Solid Earth↗

USGS2022

The geochemical and textural transition between the Reef Package and its hanging wall, Stillwater Complex, Montana, USA

The highest grade Pd-Pt deposit on Earth, the J-M Reef, is hosted in coarse-grained to pegmatoidal cumulates called the Reef Package. Decades of mine development of the J-M Reef have revealed that a distinct discontinuity in rock fabric marks the top of the rock unit that hosts economic-grade sulfide mineralization. Mine geologists refer to this discontinuity as the hanging wall contact. This contact is the top of the Reef Package and is always locatable—either by the change in rock fabric or by distinctive hanging wall textures of silicate minerals—even when the reef sulfide mineralization is absent. This rather subtle textural feature is used reliably by mine geologists to follow the Reef during exploration and mine development. Although some high tenor sulfides (>1000 ppm Pd in 100% sulfide) are found sporadically in the hanging wall cumulates, these accumulations are too small to be economically viable. We present quantitative rock fabric data for four Reef Package and hanging wall intersections collected by electron back-scattered diffraction (EBSD). Plagioclase fabrics in the hanging wall are characterized by low variance in grain sizes and a strong point maximum concentration of (010) and a perpendicular girdle distribution of [100] consistent with an axial B-type fabric. These fabrics are indicative of either compaction of the crystal mush or crystal settling of nucleated crystals, the bulk magma in a chamber. Conversely, the fabrics of the Reef Package show higher variance grain in size distributions and weak to undeveloped preferred orientation of plagioclase crystals that did not undergo significant alignment or textural equilibration of plagioclase grains. The absence of foliation in the Reef Package stands in contrast both to hanging wall fabrics and to other reported EBSD datasets of plagioclase crystals orientations from the Bushveld Complex, the Skaergaard Intrusion, and the Rum Intrusion. Furthermore, plagioclase crystal size distributions for the Reef Package show flatter slopes and convex profiles with fewer crystals at small size fractions indicating the dissolution of small crystals during partial melting and textural coarsening (i.e . Ostwald ripening) and crystal growth. Crystal growth was favored over the nucleation of new crystals during prolonged interaction with a hot infiltrating melt into the resident mush resulting in the coarse-grained textures of the Reef Package cumulates. The hanging wall contact represents a boundary between partially remelted crystal mush of the Reef Package, where sulfide mineralization formed and accumulated, and an overlying essentially barren cumulate pile. The hanging wall cumulates formed following the cessation of footwall erosion and the resumption of crystal accumulation by normal magma chamber processes.

Montana↗

USGS2014

U.S. Geological Survey Karst Interest Group Proceedings, Carlsbad, New Mexico, April 29-May 2, 2014

Karst aquifer systems are present throughout parts of the United States and some of its territories, and have developed in carbonate rocks (primarily limestone and dolomite) that span an interval of time encompassing more than 550 million years. The depositional environments, diagenetic processes, post-depositional tectonic events, and geochemical weathering processes that form karst aquifers are varied and complex, and involve biological, chemical, and physical changes. These factors, combined with the diverse climatic regimes under which karst development in these rocks has taken place, result in the unique dual- or triple-porosity nature of karst aquifers. These complex hydrogeologic systems typically represent challenging and unique conditions to scientists attempting to study groundwater flow and contaminant transport in these terrains. The dissolution of carbonate rocks and the subsequent development of distinct and beautiful landscapes, caverns, and springs has resulted in the most exceptional karst areas of the United States being designated as national or state parks; commercial caverns and known privately owned caves number in the tens of thousands. Both public and private properties provide access for scientists to study the flow of groundwater in situ. Likewise, the range and complexity of landforms and groundwater flow systems associated with karst terrains are enormous, perhaps more than for any other aquifer type. Karst aquifers and landscapes that form in tropical areas, such as the cockpit karst along the north coast of Puerto Rico, differ greatly from karst landforms in more arid climates, such as the Edwards Plateau in west-central Texas or the Guadalupe Mountains near Carlsbad, New Mexico, where hypogenic processes have played a major role in speleogenesis. Many of these public and private lands also contain unique flora and fauna associated with these karst hydrogeologic systems. As a result, numerous federal, state, and local agencies have a strong interest in the study of karst terrains. Many of the major springs and aquifers in the United States have developed in carbonate rocks, such as the Floridan aquifer system in Florida and parts of Alabama, Georgia, and South Carolina; the Ozark Plateaus aquifer system in parts of Arkansas, Kansas, Missouri, and Oklahoma; and the Edwards-Trinity aquifer system in west-central Texas. These aquifers, and the springs that discharge from them, serve as major water-supply sources and as unique ecological habitats. Competition for the water resources of karst aquifers is common, and urban development and the lack of attenuation of contaminants in karst areas can impact the ecosystem and water quality of these aquifers. The concept for developing a platform for interaction among scientists within the U.S. Geological Survey (USGS) working on karst-related studies evolved from the November 1999 National Ground-Water Meeting of the USGS. As a result, the Karst Interest Group (KIG) was formed in 2000. The KIG is a loose-knit, grass-roots organization of USGS and non-USGS scientists and researchers devoted to fostering better communication among scientists working on, or interested in, karst science. The primary mission of the KIG is to encourage and support interdisciplinary collaboration and technology transfer among scientists working in karst areas. Additionally, the KIG encourages collaborative studies between the different mission areas of the USGS as well as other federal and state agencies, and with researchers from academia and institutes. The KIG also encourages younger scientists by participation of students in the poster and oral sessions. To accomplish its mission, the KIG has organized a series of workshops that are held near nationally important karst areas. To date (2014) six KIG workshops, including the workshop documented in this report, have been held. The workshops typically include oral and poster sessions on selected karst-related topics and research, as well as field trips to local karst features. Proceedings of the workshops are published by the USGS and are available online at http://water.usgs.gov/ogw/karst/kig . The first KIG workshop was held in St. Petersburg, Florida, February 13–16, 2001, in the vicinity of the large springs and other karst features of the Floridan aquifer system. The second KIG workshop was held August 20–22, 2002, in Shepherdstown, West Virginia, in proximity to the carbonate aquifers of the northern Shenandoah Valley and highlighted an invited presentation on karst literature by the late Barry F. Beck of P.E. LaMoreaux and Associates. The third KIG workshop was held September 12–15, 2005, in Rapid City, South Dakota, nearby to karst features in evaporites and limestones of the Madison Group in the Black Hills of South Dakota, including Wind Cave National Park and Jewel Cave National Monument. The workshop also included a featured presentation by Thomas Casadevall, Central Region Director, USGS, on the status of earth science at the USGS and evening trips to Jewel Cave led by Mike Wiles, National Park Service (NPS) and Wind Cave led by Rod Horrocks, NPS. The fourth KIG workshop was held May 27–29, 2008, and hosted by the Hoffman Environmental Research Institute and Center for Cave and Karst Studies at Western Kentucky University in Bowling Green, Kentucky, near Mammoth Cave National Park and karst features of the Chester Upland and Pennyroyal Plateau. The workshop featured a late-night field trip into Mammoth Cave with Rickard Toomey and Rick Olsen, NPS. The fifth workshop was held April 26–29, 2011, and was a joint meeting of the USGS KIG and University of Arkansas HydroDays, hosted by the Department of Geosciences at the University of Arkansas in Fayetteville. The workshop featured an outstanding field trip to the unique karst terrain along the Buffalo National River of the southern Ozarks and a keynote presentation on paleokarst in the United States by Art and Peggy Palmer. This sixth and current 2014 KIG workshop is hosted by the National Cave and Karst Research Institute (NCKRI) in Carlsbad, New Mexico, with Director of NCKRI, George Veni, serving as co-chair of the workshop with Eve Kuniansky, USGS. The session planning committee for this sixth workshop includes Van Brahana, USGS retired and University of Arkansas Professor Emeritus; Tom Byl, USGS and Tennessee State University; Zelda Bailey, former Director of NCKRI and retired Director, National Institute of Standards and Technology, Boulder Laboratory, Colorado; Patrick Tucci, USGS retired; and Mike Bradley, Allan Clark, Geoff Delin, Daniel Doctor, James Kaufmann, Eve Kuniansky, Randy Orndorff, Larry Spangler, and Dave Weary of the USGS. The karst hydrology field trip on Thursday will be led by Lewis Land (NCKRI karst hydrologist) and the optional Friday field trip on the geology of Carlsbad Caverns National Park will be led by George Veni. The keynote speaker is Dr. Penelope Boston, Director of Cave and Karst Studies at New Mexico Tech, Socorro, and Academic Director at NCKRI, who will address the future of karst research. Additionally, there is a featured presentation “Irish karst and its management,” by Caoimhe Hickey, The Geological Survey of Ireland, preceding a panel discussion on “Collaboration During Times of Limited Resources.” The extended abstracts of USGS authors were peer reviewed and approved for publication by the U.S. Geological Survey. Articles submitted by university researchers and other federal and state agencies did not go through the formal USGS peer review and approval process, and therefore may not adhere to our editorial standards or stratigraphic nomenclature and is not research conducted or data collected by the USGS. However, all articles had at a minimum of two peer reviews, and all articles were edited for consistency of appearance in the published Proceedings. The use of trade, firm or product names in any article is for descriptive purposes only and does not imply endorsement by the U.S. Government. The USGS, Office of Groundwater, provides technical support for the Karst Interest Group website and public availability of the Proceedings from these workshops, and the USGS Groundwater Resources Program funds the publication costs. Finally, the cover illustration is the work of Ann Tihansky, USGS, used since the first KIG workshop in 2000.

Scientific Investigations Report↗

USGS2007

Petrology and tectonics of Phanerozoic continent formation: From island arcs to accretion and continental arc magmatism

Mesozoic continental arcs in the North American Cordillera were examined here to establish a baseline model for Phanerozoic continent formation. We combine new trace-element data on lower crustal xenoliths from the Mesozoic Sierra Nevada Batholith with an extensive grid-based geochemical map of the Peninsular Ranges Batholith, the southern equivalent of the Sierras. Collectively, these observations give a three-dimensional view of the crust, which permits the petrogenesis and tectonics of Phanerozoic crust formation to be linked in space and time. Subduction of the Farallon plate beneath North America during the Triassic to early Cretaceous was characterized by trench retreat and slab rollback because old and cold oceanic lithosphere was being subducted. This generated an extensional subduction zone, which created fringing island arcs just off the Paleozoic continental margin. However, as the age of the Farallon plate at the time of subduction decreased, the extensional environment waned, allowing the fringing island arc to accrete onto the continental margin. With continued subduction, a continental arc was born and a progressively more compressional environment developed as the age of subducting slab continued to young. Refinement into a felsic crust occurred after accretion, that is, during the continental arc stage, wherein a thickened crustal and lithospheric column permitted a longer differentiation column. New basaltic arc magmas underplate and intrude the accreted terrane, suture, and former continental margin. Interaction of these basaltic magmas with pre-existing crust and lithospheric mantle created garnet pyroxenitic mafic cumulates by fractional crystallization at depth as well as gabbroic and garnet pyroxenitic restites at shallower levels by melting of pre-existing lower crust. The complementary felsic plutons formed by these deep-seated differentiation processes rose into the upper crust, stitching together the accreted terrane, suture and former continental margin. The mafic cumulates and restites, owing to their high densities, eventually foundered into the mantle, leaving behind a more felsic crust. Our grid-based sampling allows us to estimate an unbiased average upper crustal composition for the Peninsular Ranges Batholith. Major and trace-element compositions are very similar to global continental crust averaged over space and time, but in detail, the Peninsular Ranges are slightly lower in compatible to mildly incompatible elements, MgO, Mg#, V, Sc, Co, and Cr. The compositional similarities suggest a strong arc component in global continental crust, but the slight discrepancies suggest that additional crust formation processes are also important in continent formation as a whole. Finally, the delaminated Sierran garnet pyroxenites have some of the lowest U/Pb ratios ever measured for silicate rocks. Such material, if recycled and stored in the deep mantle, would generate a reservoir with very unradiogenic Pb, providing one solution to the global Pb isotope paradox.

Earth and Planetary Science Letters↗

USGS2021

Self-limitation of sand storage in a bedrock-canyon river arising from the interaction of flow and grain size

Bedrock-canyon rivers tend to be supply limited because they are efficient transporters of sediment and not because the upstream supply of sediment is small. A byproduct of this supply limitation is that the finer alluvium stored in these rivers has shorter residence times and smaller volumes than in alluvial rivers. To improve our understanding of disequilibrium sediment transport and its effect on sand storage in bedrock-canyon rivers, we undertook a 20-year study, synthesized herein, of the Colorado River in Grand Canyon. Despite the large loads for which it was renowned, this river exhibited evidence of natural sand-supply limitation and became the perfect natural laboratory for studying sand transport in a bedrock canyon after upstream dam construction exacerbated this supply limitation. During our study, we made and analyzed an unprecedented ∼2.5 million measurements of the suspended and bed sediment. Results indicate that sand storage in this bedrock-canyon river is self-limiting owing to the physical controls of flow and grain size causing negative feedbacks that likely also operate in other bedrock-canyon rivers. Following episodic tributary floods that supply finer sand, sand migrates quickly downstream in the form of a wave in which large systematic changes in bed-sand grain size occur. These grain-size changes cause discharge-independent systematic changes in suspended-sand concentration in excess of a factor of 20. Although the tributary supply of sand increases the amount of sand storage, it also greatly increases the downstream sand transport by causing bed-sand fining, thus limiting the residence time and volume of sand storage.

Arizona↗

USGS2008

Categorizing natural disaster damage assessment using satellite-based geospatial techniques

Remote sensing of a natural disaster's damage offers an exciting backup and/or alternative to traditional means of on-site damage assessment. Although necessary for complete assessment of damage areas, ground-based damage surveys conducted in the aftermath of natural hazard passage can sometimes be potentially complicated due to on-site difficulties (e.g., interaction with various authorities and emergency services) and hazards (e.g., downed power lines, gas lines, etc.), the need for rapid mobilization (particularly for remote locations), and the increasing cost of rapid physical transportation of manpower and equipment. Satellite image analysis, because of its global ubiquity, its ability for repeated independent analysis, and, as we demonstrate here, its ability to verify on-site damage assessment provides an interesting new perspective and investigative aide to researchers. Using one of the strongest tornado events in US history, the 3 May 1999 Oklahoma City Tornado, as a case example, we digitized the tornado damage path and co-registered the damage path using pre- and post-Landsat Thematic Mapper image data to perform a damage assessment. We employed several geospatial approaches, specifically the Getis index, Geary's C , and two lacunarity approaches to categorize damage characteristics according to the original Fujita tornado damage scale (F-scale). Our results indicate strong relationships between spatial indices computed within a local window and tornado F-scale damage categories identified through the ground survey. Consequently, linear regression models, even incorporating just a single band, appear effective in identifying F-scale damage categories using satellite imagery. This study demonstrates that satellite-based geospatial techniques can effectively add spatial perspectives to natural disaster damages, and in particular for this case study, tornado damages.

Natural Hazards and Earth System Sciences↗

USGS2023

Time-lapse seafloor surveys reveal how turbidity currents and internal tides in Monterey Canyon interact with the seabed at centimeter-scale

Here we show how ultra-high resolution seabed mapping using new technology can help to understand processes that sculpt submarine canyons. Time-lapse seafloor surveys were conducted in the axis of Monterey Canyon, ∼50 km from the canyon head (∼1,840 m water depth) over an 18-month period. These surveys comprised 5-cm resolution multibeam bathymetry, 1-cm resolution lidar bathymetry, and 2-mm resolution stereophotographic imagery. Bathymetry data reveal centimeter-scale textures that would be undetectable by more traditional survey methods. Upward-looking Acoustic Doppler Current Profilers at the site recorded the flow character of internal tides and the passage of three turbidity currents, while sediment cores collected from the site record flow deposits. Combined with flow and core data, the bathymetry shows how turbidity currents and internal tides modify the seabed. The turbidity currents drape sediment across the site, infilling bedform troughs and smoothing erosional features carved by the internal tides (e.g., rippled scours). Turbidity currents with speeds of 0.9–3.3 m/s failed to cause notable bedform movement, which is surprising given that flows with similar speeds produced rapid bedform migration elsewhere, including the upper Monterey Canyon. The lack of migration may be related to the character of the underlying substrate or indicate that turbidity currents at the site lack dense, near-bed layers. The scale of scours produced by the internal tides (≤0.7 m/s) approaches the scale of features recorded in the ancient rock record. Thus, these results illustrate how the scale gap between seabed mapping technology and the rock record may eventually be bridged.

California↗

USGS2018

Technical note: False low turbidity readings from optical probes during high suspended-sediment concentrations

Turbidity, a measure of water clarity, is monitored for a variety of purposes including (1) to help determine whether water is safe to drink, (2) to establish background conditions of lakes and rivers and detect pollution caused by construction projects and stormwater discharge, (3) to study sediment transport in rivers and erosion in catchments, (4) to manage siltation of water reservoirs, and (5) to establish connections with aquatic biological properties, such as primary production and predator–prey interactions. Turbidity is typically measured with an optical probe that detects light scattered from particles in the water. Probes have defined upper limits of the range of turbidity that they can measure. The general assumption is that when turbidity exceeds this upper limit, the values of turbidity will be constant, i.e., the probe is pegged ; however, this assumption is not necessarily valid. In rivers with limited variation in the physical properties of the suspended sediment, at lower suspended-sediment concentrations, an increase in suspended-sediment concentration will cause a linear increase in turbidity. When the suspended-sediment concentration in these rivers is high, turbidity levels can exceed the upper measurement limit of an optical probe and record a constant pegged value. However, at extremely high suspended-sediment concentrations, optical turbidity probes do not necessarily stay pegged at a constant value. Data from the Colorado River in Grand Canyon, Arizona, USA, and a laboratory experiment both demonstrate that when turbidity exceeds instrument-pegged conditions, increasing suspended-sediment concentration (and thus increasing turbidity) may cause optical probes to record decreasing false turbidity values that appear to be within the valid measurement range of the probe. Therefore, under high-turbidity conditions, other surrogate measurements of turbidity (e.g., acoustic-attenuation measurements or suspended-sediment samples) are necessary to correct these low false turbidity measurements and accurately measure turbidity.

Hydrology and Earth System Sciences↗

USGS2008

Probabilistic seismic hazard in the San Francisco Bay area based on a simplified viscoelastic cycle model of fault interactions

We construct a viscoelastic cycle model of plate boundary deformation that includes the effect of time-dependent interseismic strain accumulation, coseismic strain release, and viscoelastic relaxation of the substrate beneath the seismogenic crust. For a given fault system, time-averaged stress changes at any point (not on a fault) are constrained to zero; that is, kinematic consistency is enforced for the fault system. The dates of last rupture, mean recurrence times, and the slip distributions of the (assumed) repeating ruptures are key inputs into the viscoelastic cycle model. This simple formulation allows construction of stress evolution at all points in the plate boundary zone for purposes of probabilistic seismic hazard analysis (PSHA). Stress evolution is combined with a Coulomb failure stress threshold at representative points on the fault segments to estimate the times of their respective future ruptures. In our PSHA we consider uncertainties in a four-dimensional parameter space: the rupture peridocities, slip distributions, time of last earthquake (for prehistoric ruptures) and Coulomb failure stress thresholds. We apply this methodology to the San Francisco Bay region using a recently determined fault chronology of area faults. Assuming single-segment rupture scenarios, we find that fature rupture probabilities of area faults in the coming decades are the highest for the southern Hayward, Rodgers Creek, and northern Calaveras faults. This conclusion is qualitatively similar to that of Working Group on California Earthquake Probabilities, but the probabilities derived here are significantly higher. Given that fault rupture probabilities are highly model-dependent, no single model should be used to assess to time-dependent rupture probabilities. We suggest that several models, including the present one, be used in a comprehensive PSHA methodology, as was done by Working Group on California Earthquake Probabilities.

Journal of Geophysical Research B: Solid Earth↗

USGS2015

Coevolution of bed surface patchiness and channel morphology: 1. Mechanisms of forced patch formation

Riverbeds frequently display a spatial structure where the sediment mixture composing the channel bed has been sorted into discrete patches of similar grain size. Even though patches are a fundamental feature in gravel bed rivers, we have little understanding of how patches form, evolve, and interact. Here we present a two-dimensional morphodynamic model that is used to examine in greater detail the mechanisms responsible for the development of forced bed surface patches and the coevolution of bed morphology and bed surface patchiness. The model computes the depth-averaged channel hydrodynamics, mixed-grain-size sediment transport, and bed evolution by coupling the river morphodynamic model Flow and Sediment Transport with Morphological Evolution of Channels (FaSTMECH) with a transport relation for gravel mixtures and the mixed-grain-size Exner equation using the active layer assumption. To test the model, we use it to simulate a flume experiment in which the bed developed a sequence of alternate bars and temporally and spatially persistent forced patches with a general pattern of coarse bar tops and fine pools. Cross-stream sediment flux causes sediment to be exported off of bars and imported into pools at a rate that balances downstream gradients in the streamwise sediment transport rate, allowing quasi-steady bar-pool topography to persist. The relative importance of lateral gravitational forces on the cross-stream component of sediment transport is a primary control on the amplitude of the bars. Because boundary shear stress declines as flow shoals over the bars, the lateral sediment transport is increasingly size selective and leads to the development of coarse bar tops and fine pools.

Journal of Geophysical Research F: Earth Surface↗