SEARCH · USGS Science

Results for “Earth System Science Data”

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 505 records · Page 28Linked to original sources

USGS2023

Plant size, latitude, and phylogeny explain within-population variability in herbivory

Interactions between plants and herbivores are central in most ecosystems, but their strength is highly variable. The amount of variability within a system is thought to influence most aspects of plant-herbivore biology, from ecological stability to plant defense evolution. Our understanding of what influences variability, however, is limited by sparse data. We collected standardized surveys of herbivory for 503 plant species at 790 sites across 116° of latitude. With these data, we show that within-population variability in herbivory increases with latitude, decreases with plant size, and is phylogenetically structured. Differences in the magnitude of variability are thus central to how plant-herbivore biology varies across macroscale gradients. We argue that increased focus on interaction variability will advance understanding of patterns of life on Earth.

Science↗

USGS2011

Processes affecting geochemistry and contaminant movement in the middle Claiborne aquifer of the Mississippi embayment aquifer system

Groundwater chemistry and tracer-based age data were used to assess contaminant movement and geochemical processes in the middle Claiborne aquifer (MCA) of the Mississippi embayment aquifer system. Water samples were collected from 30 drinking-water wells (mostly domestic and public supply) and analyzed for nutrients, major ions, pesticides, volatile organic compounds (VOCs), and transient age tracers (chlorofluorocarbons, tritium and helium-3, and sulfur hexafluoride). Redox conditions are highly variable throughout the MCA. However, mostly oxic groundwater with low dissolved solids is more vulnerable to nitrate contamination in the outcrop areas east of the Mississippi River in Mississippi and west Tennessee than in mostly anoxic groundwater in downgradient areas in western parts of the study area. Groundwater in the outcrop area was relatively young (apparent age of less than 40 years) with significantly ( p < 0.05) higher dissolved oxygen and nitrate–N concentrations and higher detections of pesticides and VOCs compared to water samples from wells in downgradient areas. Oxygen reduction and denitrification rates were low compared to other aquifers in the United States (zero order rate constants for oxygen reduction and denitrification were 4.7 and 5–10 μmol/L/year, respectively). Elevated concentrations of nitrate–N, and detections of pesticides and VOCs in some deep public supply wells (>50 m depth) indicated contaminant movement from shallow parts of the aquifer into deeper oxic zones. Given the persistence of nitrate in young oxic groundwater that was recharged several decades ago, and the lack of a confining unit, the downward movement of young contaminated water may result in higher nitrate concentrations over time in deeper parts of the aquifer containing older oxic water.

Claiborne Aquifer;Mississippi Embayment↗

USGS2024

The National Ocean Biodiversity Strategy

President Biden has been clear that the ocean is central to life on Earth. As he has proclaimed, “the ocean powers millions of jobs; feeds and sustains us; and is a rejuvenating source of inspiration, exploration, and recreation.” The Biden-Harris Administration has worked hard to fulfill the President’s goal to protect and conserve at least 30% of U.S. waters by 2030. The ocean faces increased threats from warming, overfishing, increased acidity, and loss of biodiversity. It is now more important than ever to sustain the many benefits that the ocean, coasts, and Great Lakes provide, including food, a favorable climate, recreation, physical and mental health, and for many, a sense of cultural identity. Ocean life represents an irreplaceable heritage, the foundation of a habitable planet, and a vast trove of resources. Keeping our ocean healthy requires reliable information on the changing status of these living organisms, the drivers of biodiversity change, and options for effectively addressing those drivers. Over 2 million species are estimated to live in the ocean, yet only about 240,000 species have been described by scientists. Most of the ocean’s benefits result from those diverse species interacting with one another and the environment they create. To protect and conserve the ocean, we as a nation need to make better use of existing knowledge and prioritize acquiring new biodiversity knowledge to enable better policy and management decisions. The ability to monitor ocean species and habitats has expanded dramatically over the past decade, with innovations in technology, genomics, taxonomy, big data management and sharing, artificial intelligence, and machine learning. Yet large fractions of the U.S. ocean remain almost unknown. The National Ocean Biodiversity Strategy (strategy) reflects the urgent need to leverage these advances. The goals of this strategy must be guided by the nation’s diverse voices and ways of knowing, in order to maximize effective and equitable stewardship of the ocean’s diverse life and its benefits to people. The strategy is intended as a guiding document for government to advance three overarching goals: ● Goal 1: Drive delivery of ocean biodiversity knowledge at the national scale. Objectives include developing an Implementation Plan for achieving the strategy’s three goals; establishing a coordination mechanism to manage the implementation; and documenting gaps in biodiversity knowledge and the benefits of ocean biodiversity to people and economies. ● Goal 2: Strengthen tools and institutions to deliver ocean biodiversity knowledge. Objectives include establishing a robust information pipeline to support indicators and dynamic maps of ocean biodiversity, from the coasts to the deep sea. This pipeline should include expanded observing systems and comprehensive data management; science and technology solutions to accelerate the availability of biodiversity information; and plans to leverage previous investments to rebuild and expand the nation’s human capital and infrastructure to sustain foundational taxonomy and biodiversity science. ● Goal 3: Protect, conserve, restore, and sustainably use ocean biodiversity. Objectives include expanding the collection, delivery, and use of biodiversity knowledge to inform actions that advance ocean protection, conservation, restoration, and sustainable development. Government should lead in establishing and incentivizing diverse partnerships across scales and sectors to implement those actions and should educate and involve the public to discover and value the nation’s diverse ocean life. Achieving these goals will require commitments across society: new federal and private investments, coordination across sectors to address climate and equity challenges, and engagement of Indigenous Knowledge holders and frontline communities as full partners throughout planning and implementation. The Subcommittee on Ocean Science and Technology (SOST) IWG-Biodiversity will begin developing an Implementation Plan to describe and direct specific actions to implement the strategy. Successful implementation of the strategy will harmonize and expand collection and delivery of timely knowledge on ocean life to all of society. The strategy will also enable evidence-based management and protection of the ocean. Advancing the strategy will build human and institutional capital and partnerships that support both existing mandates and new needs to rebuild and sustain biodiversity, achieve healthy ocean ecosystems, and manage living resources. Implementing the strategy will deliver knowledge for monitoring, modeling, forecasting, and assessments that support food security, public health, and cultural values, and that more effectively protect, conserve, and restore nature.

Report↗

USGS2020

Littoral sediment from rivers: Patterns, rates and processes of river mouth morphodynamics

Rivers provide important sediment inputs to many littoral cells, thereby replenishing sand and gravel of beaches around the world. However, there is limited information about the patterns and processes of littoral-grade sediment transfer from rivers into coastal systems. Here I address these information gaps by examining topographic and bathymetric data of river mouths and constructing sediment budgets to characterize time-dependent patterns of onshore, offshore, and alongshore transport. Two river deltas, which differ in their morphology, were used in this study: the Elwha River, Washington, which builds a mixed sediment Gilbert-style delta, and the Santa Clara River, California, which builds a cross-shore dispersed sand delta from hyperpycnal flows. During and after sediment discharge events, both systems exhibited a similar evolution composed of three phases: (i) submarine delta growth during offshore transport of river sediment, (ii) onshore-dominated transport from the submarine delta to a subaerial river mouth berm, and (iii) longshore-dominated transport away from the river mouth following subaerial berm development. Although stage (ii) occurred within days to weeks for the systems studied and was associated with the greatest rates of net erosion and deposition, onshore transport of sediment from submarine deposit to the beach persisted for years following the river discharge event. These morphodynamics were similar to simple equilibrium profile concepts that were modified with an onshore-dominated cross-shore transport rule. Additionally, both study sites revealed that littoral-grade sediment was initially exported to depths beyond the active littoral cell (i.e., below the depth of closure) during the stage (i). Following several years of reworking by coastal processes, bathymetric surveys suggested that 14 and 46% of the original volume of littoral-grade sediment discharged by the Santa Clara and Elwha Rivers, respectively, continued to be below the depth of closure. Combined, this suggests that integration of river sediment into a littoral cell can be a multi-year process and that the full volume of littoral-grade sediment discharged by small rivers may not be integrated into littoral cells because of sand and gravel “losses” to the continental shelf.

California, Washington↗

USGS2024

Automatic identification and quantification of volcanic hotspots in Alaska using HotLINK: The hotspot learning and identification network

An increase in volcanic thermal emissions can indicate subsurface and surface processes that precede, or coincide with, volcanic eruptions. Space-borne infrared sensors can detect hotspots—defined here as localized volcanic thermal emissions—in near-real-time. However, automatic hotspot detection systems are needed to efficiently analyze the large quantities of data produced. While hotspots have been automatically detected for over 20 years with simple thresholding algorithms, new computer vision technologies, such as convolutional neural networks (CNNs), can enable improved detection capabilities. Here we introduce HotLINK: the Hotspot Learning and Identification Network, a CNN trained to detect hotspots with a dataset of −3,800 satellite-based, Visible Infrared Imaging Radiometer Suite (VIIRS) images from Mount Veniaminof and Mount Cleveland volcanoes, Alaska. We find that our model achieves an accuracy of 96% (F1-score 0.92) when evaluated on −1,700 unseen images from the same volcanoes, and 95% (F1-score 0.67) when evaluated on −3,000 images from six additional Alaska volcanoes (Augustine Volcano, Bogoslof Island, Okmok Caldera, Pavlof Volcano, Redoubt Volcano, Shishaldin Volcano). In comparison with an existing threshold-based hotspot detection algorithm, MIROVA (Coppola et al., Geological Society, London, Special Publications, 2016, 426, 181–205), our model detects 22% more hotspots and produces 12% fewer false positives. Additional testing on −700 labeled Moderate Resolution Imaging Spectroradiometer (MODIS) images from Mount Veniaminof demonstrates that our model is applicable to this sensor’s data as well, achieving an accuracy of 98% (F1-score 0.95). We apply HotLINK to 10 years of VIIRS data and 22 years of MODIS data for the eight aforementioned Alaska volcanoes and calculate the radiative power of detected hotspots. From these time series we find that HotLINK accurately characterizes background and eruptive periods, similar to MIROVA, but also detects more subtle warming signals, potentially related to volcanic unrest. We identify three advantages to our model over its predecessors: 1) the ability to detect more subtle volcanic hotspots and produce fewer false positives, especially in daytime images; 2) probabilistic predictions provide a measure of detection confidence; and 3) its transferability, i.e., the successful application to multiple sensors and multiple volcanoes without the need for threshold tuning, suggesting the potential for global application.

Alaska↗

USGS2025

Invited perspectives: Integrating hydrologic information into the next generation of landslide early warning systems

Although rainfall-triggered landslides are initiated by subsurface hydro-mechanical processes related to the loading, weakening, and eventual failure of slope materials, most landslide early warning systems (LEWS) have relied solely on rainfall event information. In previous decades, several studies demonstrated the value of integrating proxies for subsurface hydrologic information to improve rainfall-based forecasting of shallow landslides. More recently, broader access to commercial sensors and telemetry for real-time data transmission has invigorated new research into hydrometeorological thresholds for LEWS. Given the increasing number of studies across the globe using hydrologic monitoring, mathematical modeling, or both in combination, it is now possible to make some insights into the advantages versus limitations of this approach. The extensive progress demonstrates the value of in situ hydrologic information for reducing both failed and false alarms, through the ability to characterize infiltration during, as well as the drainage and drying processes between major storm events. There are also some areas for caution surrounding the long-term sustainability of subsurface monitoring in landslide-prone terrain, as well as unresolved questions in hillslope hydrologic modeling, which relies heavily on the assumptions of diffuse flow and vertical infiltration but often ignores preferential flow and lateral drainage. Here, we share a collective perspective based on our previous collaborative work across Europe, North America, Africa, and Asia to discuss these challenges and provide some guidelines for integrating knowledge of hydrology and climate into the next generation of LEWS. We propose that the greatest opportunity for improvement is through a measure-and-model approach to develop an understanding of landslide hydro-climatology that accounts for local controls on subsurface storage dynamics. Additionally, new efforts focused on the subsurface hydrology are complementary to existing rainfall-based methods, so leveraging these with near-term precipitation forecasts is a priority for increasing lead times.

Alaska↗

USGS2013

Topological and canonical kriging for design flood prediction in ungauged catchments: an improvement over a traditional regional regression approach?

In the United States, estimation of flood frequency quantiles at ungauged locations has been largely based on regional regression techniques that relate measurable catchment descriptors to flood quantiles. More recently, spatial interpolation techniques of point data have been shown to be effective for predicting streamflow statistics (i.e., flood flows and low-flow indices) in ungauged catchments. Literature reports successful applications of two techniques, canonical kriging, CK (or physiographical-space-based interpolation, PSBI), and topological kriging, TK (or top-kriging). CK performs the spatial interpolation of the streamflow statistic of interest in the two-dimensional space of catchment descriptors. TK predicts the streamflow statistic along river networks taking both the catchment area and nested nature of catchments into account. It is of interest to understand how these spatial interpolation methods compare with generalized least squares (GLS) regression, one of the most common approaches to estimate flood quantiles at ungauged locations. By means of a leave-one-out cross-validation procedure, the performance of CK and TK was compared to GLS regression equations developed for the prediction of 10, 50, 100 and 500 yr floods for 61 streamgauges in the southeast United States. TK substantially outperforms GLS and CK for the study area, particularly for large catchments. The performance of TK over GLS highlights an important distinction between the treatments of spatial correlation when using regression-based or spatial interpolation methods to estimate flood quantiles at ungauged locations. The analysis also shows that coupling TK with CK slightly improves the performance of TK; however, the improvement is marginal when compared to the improvement in performance over GLS.

Hydrology and Earth System Sciences↗

USGS2009

High-frequency Rayleigh-wave method

High-frequency (≥2 Hz) Rayleigh-wave data acquired with a multichannel recording system have been utilized to determine shear (S)-wave velocities in near-surface geophysics since the early 1980s. This overview article discusses the main research results of high-frequency surface-wave techniques achieved by research groups at the Kansas Geological Survey and China University of Geosciences in the last 15 years. The multichannel analysis of surface wave (MASW) method is a non-invasive acoustic approach to estimate near-surface S-wave velocity. The differences between MASW results and direct borehole measurements are approximately 15% or less and random. Studies show that simultaneous inversion with higher modes and the fundamental mode can increase model resolution and an investigation depth. The other important seismic property, quality factor ( Q ), can also be estimated with the MASW method by inverting attenuation coefficients of Rayleigh waves. An inverted model (S-wave velocity or Q ) obtained using a damped least-squares method can be assessed by an optimal damping vector in a vicinity of the inverted model determined by an objective function, which is the trace of a weighted sum of model-resolution and model-covariance matrices. Current developments include modeling high-frequency Rayleigh-waves in near-surface media, which builds a foundation for shallow seismic or Rayleigh-wave inversion in the time-offset domain; imaging dispersive energy with high resolution in the frequency-velocity domain and possibly with data in an arbitrary acquisition geometry, which opens a door for 3D surface-wave techniques; and successfully separating surface-wave modes, which provides a valuable tool to perform S-wave velocity profiling with high-horizontal resolution.

Earth Science↗

USGS2018

Bias correction of simulated historical daily streamflow at ungauged locations by using independently estimated flow duration curves

In many simulations of historical daily streamflow distributional bias arising from the distributional properties of residuals has been noted. This bias often presents itself as an underestimation of high streamflow and an overestimation of low streamflow. Here, 1168 streamgages across the conterminous USA, having at least 14 complete water years of daily data between 1 October 1980 and 30 September 2013, are used to explore a method for rescaling simulated streamflow to correct the distributional bias. Based on an existing approach that separates the simulated streamflow into components of temporal structure and magnitude, the temporal structure is converted to simulated nonexceedance probabilities and the magnitudes are rescaled using an independently estimated flow duration curve (FDC) derived from regional regression. In this study, this method is applied to a pooled ordinary kriging simulation of daily streamflow coupled with FDCs estimated by regional regression on basin characteristics. The improvement in the representation of high and low streamflows is correlated with the accuracy and unbiasedness of the estimated FDC. The method is verified by using an idealized case; however, with the introduction of regionally regressed FDCs developed for this study, the method is only useful overall for the upper tails, which are more accurately and unbiasedly estimated than the lower tails. It remains for future work to determine how accurate the estimated FDCs need to be to be useful for bias correction without unduly reducing accuracy. In addition to its potential efficacy for distributional bias correction, this particular instance of the methodology also represents a generalization of nonlinear spatial interpolation of daily streamflow using FDCs. Rather than relying on single index stations, as is commonly done to reflect streamflow timing, this approach to simulation leverages geostatistical tools to allow a region of neighbors to reflect streamflow timing.

Hydrology and Earth System Sciences↗

USGS2015

What lies beneath: geophysical mapping of a concealed Precambrian intrusive complex along the Iowa–Minnesota border

Large-amplitude gravity and magnetic highs over northeast Iowa are interpreted to reflect a buried intrusive complex composed of mafic–ultramafic rocks, the northeast Iowa intrusive complex (NEIIC), intruding Yavapai province (1.8–1.72 Ga) rocks. The age of the complex is unproven, although it has been considered to be Keweenawan (∼1.1 Ga). Because only four boreholes reach the complex, which is covered by 200–700 m of Paleozoic sedimentary rocks, geophysical methods are critical to developing a better understanding of the nature and mineral resource potential of the NEIIC. Lithologic and cross-cutting relations interpreted from high-resolution aeromagnetic and airborne gravity gradient data are presented in the form of a preliminary geologic map of the basement Precambrian rocks. Numerous magnetic anomalies are coincident with airborne gravity gradient (AGG) highs, indicating widespread strongly magnetized and dense rocks of likely mafic–ultramafic composition. A Yavapai-age metagabbro unit is interpreted to be part of a layered intrusion with subvertical dip. Another presumed Yavapai unit has low density and weak magnetization, observations consistent with felsic plutons. Northeast-trending, linear magnetic lows are interpreted to reflect reversely magnetized diabase dikes and have properties consistent with Keweenawan rocks. The interpreted dikes are cut in places by normally magnetized mafic–ultramafic rocks, suggesting that the latter represent younger Keweenawan rocks. Distinctive horseshoe-shaped magnetic and AGG highs correspond with a known gabbro, and surround rocks with weaker magnetization and lower density. Here, informally called the Decorah complex, the source body has notable geophysical similarities to Keweenawan alkaline ring complexes, such as the Coldwell and Killala Lake complexes, and Mesoproterozoic anorogenic complexes, such as the Kiglapait, Hettasch, and Voisey’s Bay intrusions in Labrador. Results presented here suggest that much of the NEIIC is composed of such complexes, and broadly speaking, may be a discontinuous group of several intrusive bodies. Most units are cut by suspected northwest-trending faults imaged as magnetic lineaments, and one produces apparent sinistral fault separation of a dike in the eastern part of the survey area. The location, trend, and apparent sinistral sense of motion are consistent with the suspected faults being part of the Belle Plaine fault zone, a complex transform fault zone within the Midcontinent rift system that is here proposed to correspond with a major structural discontinuity.

Iowa, Minnesota↗

USGS2016

Holocene climate changes in eastern Beringia (NW North America) – A systematic review of multi-proxy evidence

Reconstructing climates of the past relies on a variety of evidence from a large number of sites to capture the varied features of climate and the spatial heterogeneity of climate change. This review summarizes available information from diverse Holocene paleoenvironmental records across eastern Beringia (Alaska, westernmost Canada and adjacent seas), and it quantifies the primary trends of temperature- and moisture-sensitive records based in part on midges, pollen, and biogeochemical indicators (compiled in the recently published Arctic Holocene database, and updated here to v2.1). The composite time series from these proxy records are compared with new summaries of mountain-glacier and lake-level fluctuations, terrestrial water-isotope records, sea-ice and sea-surface-temperature analyses, and peatland and thaw-lake initiation frequencies to clarify multi-centennial- to millennial-scale trends in Holocene climate change. To focus the synthesis, the paleo data are used to frame specific questions that can be addressed with simulations by Earth system models to investigate the causes and dynamics of past and future climate change. This systematic review shows that, during the early Holocene (11.7–8.2 ka; 1 ka = 1000 cal yr BP), rather than a prominent thermal maximum as suggested previously, temperatures were highly variable, at times both higher and lower than present (approximate mid-20th-century average), with no clear spatial pattern. Composited pollen, midge and other proxy records average out the variability and show the overall lowest summer and mean-annual temperatures across the study region during the earliest Holocene, followed by warming over the early Holocene. The sparse data available on early Holocene glaciation show that glaciers in southern Alaska were as extensive then as they were during the late Holocene. Early Holocene lake levels were low in interior Alaska, but moisture indicators show pronounced differences across the region. The highest frequency of both peatland and thaw-lake initiation ages also occurred during the early Holocene. During the middle Holocene (8.2–4.2 ka), glaciers retreated as the regional average temperature increased to a maximum between 7 and 5 ka, as reflected in most proxy types. Following the middle Holocene thermal maximum, temperatures decreased starting between 4 and 3 ka, signaling the onset of Neoglacial cooling. Glaciers in the Brooks and Alaska Ranges advanced to their maximum Holocene extent as lakes generally rose to modern levels. Temperature differences for averaged 500-year time steps typically ranged by 1–2 °C for individual records in the Arctic Holocene database, with a transition to a cooler late Holocene that was neither abrupt nor spatially coherent. The longest and highest-resolution terrestrial water isotope records previously interpreted to represent changes in the Aleutian low-pressure system around this time are here shown to be largely contradictory. Furthermore, there are too few records with sufficient resolution to identify sub-centennial-scale climate anomalies, such as the 8.2 ka event. The review concludes by suggesting some priorities for future paleoclimate research in the region.

Quaternary Science Reviews↗

USGS2024

Influence of inherited structure on flexural extension in foreland basin systems: Evidence from the northern Arkoma basin and southern Ozark dome, USA

Extensional faults are key components of foreland basin systems. They form within the upper crust in response to flexure of the lithosphere and accommodate subsidence within the foredeep and forebulge depozones. Such faults are excellent proxies for orogenic system evolution and control the distribution of natural resources and hazards. However, the spatiotemporal evolution of flexural extension has not been documented previously at a regional scale, thereby limiting our understanding of underlying geodynamic controls. Here, we resolve late Paleozoic flexural extension in the northern Arkoma basin and southern Ozark dome, USA. We synthesize a large database of previous mapping, existing research, subsurface data, and geophysical data into 3D geologic and 2D kinematic models. Mesh surfaces representing several key horizons from the Carboniferous Period (ca. 335-306 Ma) were constructed. These surfaces were built from oil and gas well tops (n = ∼10,000) and surface geologic map contacts using an advanced kriging method. The mesh surfaces are offset by a complex 3D fault network, allowing detailed analysis of along-strike and down-dip variations in fault displacement. Analysis of the 3D model reveals a regular and repeated fault segmentation pattern wherein E -W striking, left- and foreland-stepping en échelon normal faults are segmented by inherited NE striking basement faults. Maximum vertical separation along the E -W normal faults is generally focused between the inherited NE-trending faults. This suggests that the inherited basement faults delocalized extensional strain during late Paleozoic normal faulting. Maximum vertical separation and fault localization may correlate to areas with high-amplitude positive magnetic anomalies interpreted as Mesoproterozoic granitic rocks. Speculative covariance of magnetic anomalies and fault displacements implies that the relatively strong basement granite concentrated stress, leading to localized faulting within the relatively thin sedimentary cover. Lastly, we show that flexural extension migrated southeast to northwest from the Chesterian-Morrowan (ca. 335-319 Ma) to the Desmoinesian (ca. 306 Ma). The migratory flexural extension may be explained by diachronous loading during Pangean assembly, or by synchronous loading but variable load compensation due to inherent factors.

Arkansas, Oklahoma↗

USGS2000

Stratigraphic and structural evolution of the Selenga Delta Accommodation Zone, Lake Baikal Rift, Siberia

Seismic reflection profiles from the Lake Baikal Rift reveal extensive details about the sediment thickness, structural geometry and history of extensional deformation and syn-rift sedimentation in this classic continental rift. The Selenga River is the largest single source of terrigenous input into Lake Baikal, and its large delta sits astride the major accommodation zone between the Central and South basins of the lake. Incorporating one of the world's largest lacustrine deltas, this depositional system is a classic example of the influence of rift basin structural segmentation on a major continental drainage. More than 3700 km of deep basin-scale multi-channel seismic reflection (MCS) data were acquired during the 1989 Russian and the 1992 Russian–American field programs. The seismic data image most of the sedimentary section, including pre-rift basement in several localities. The MCS data reveal that the broad bathymetric saddle between these two major half-graben basins is underlain by a complex of severely deformed basement blocks, and is not simply a consequence of long-term deltaic deposition. Maximum sediment thickness is estimated to be more than 9 km in some areas around the Selenga Delta. Detailed stratigraphic analyses of the Selenga area MCS data suggest that modes of deposition have shifted markedly during the history of the delta. The present mode of gravity- and mass-flow sedimentation that dominates the northern and southern parts of the modern delta, as well as the pronounced bathymetric relief in the area, are relatively recent developments in the history of the Lake Baikal Rift. Several episodes of major delta progradation, each extending far across the modern rift, can be documented in the MCS data. The stratigraphic framework defined by these prograding deltaic sequences can be used to constrain the structural as well as depositional evolution of this part of the Baikal Rift. An age model has been established for this stratigraphy, by tying the delta sequences to the site of the Baikal Drilling Project 1993 Drill Hole. Although the drill hole is only 100 m deep, and the base of the cores is only ∼670 ka in age, ages were extrapolated to deeper stratigraphic intervals using the Reflection-Seismic-Radiocarbon method of Cohen et al. (1993). The deep prograding delta sequences now observed in the MCS data probably formed in response to major fluctuations in sediment supply, rather than in response to shifts in lake level. This stratigraphic framework and age model suggest that the deep delta packages developed at intervals of approximately 400 ka and may have formed as a consequence of climate changes affiliated with the northern hemisphere glaciations. The stratigraphic analysis also suggests that the Selenga Basin and Syncline developed as a distinct depocentre only during the past ∼2–3 Ma.

Lake Baikal Rift↗

USGS2013

Assessing winter cover crop nutrient uptake efficiency using a water quality simulation model

Winter cover crops are an effective conservation management practice with potential to improve water quality. Throughout the Chesapeake Bay Watershed (CBW), which is located in the Mid-Atlantic US, winter cover crop use has been emphasized and federal and state cost-share programs are available to farmers to subsidize the cost of winter cover crop establishment. The objective of this study was to assess the long-term effect of planting winter cover crops at the watershed scale and to identify critical source areas of high nitrate export. A physically-based watershed simulation model, Soil and Water Assessment Tool (SWAT), was calibrated and validated using water quality monitoring data and satellite-based estimates of winter cover crop species performance to simulate hydrological processes and nutrient cycling over the period of 1991–2000. Multiple scenarios were developed to obtain baseline information on nitrate loading without winter cover crops planted and to investigate how nitrate loading could change with different winter cover crop planting scenarios, including different species, planting times, and implementation areas. The results indicate that winter cover crops had a negligible impact on water budget, but significantly reduced nitrate leaching to groundwater and delivery to the waterways. Without winter cover crops, annual nitrate loading was approximately 14 kg ha −1 , but it decreased to 4.6–10.1 kg ha −1 with winter cover crops resulting in a reduction rate of 27–67% at the watershed scale. Rye was most effective, with a potential to reduce nitrate leaching by up to 93% with early planting at the field scale. Early planting of winter cover crops (~30 days of additional growing days) was crucial, as it lowered nitrate export by an additional ~2 kg ha −1 when compared to late planting scenarios. The effectiveness of cover cropping increased with increasing extent of winter cover crop implementation. Agricultural fields with well-drained soils and those that were more frequently used to grow corn had a higher potential for nitrate leaching and export to the waterways. This study supports the effective implement of winter cover crop programs, in part by helping to target critical pollution source areas for winter cover crop implementation.

Chesapeake Bay Watershed↗

USGS2008

(U-Th)/He dating of kimberlites-A case study from north-eastern Kansas

Dating kimberlite intrusive rocks by radiogenic isotope geochronology often is a difficult task, complicated by both the lack of dateable minerals within kimberlite as well as significant sample alteration that can degrade samples and alter parent-daughter ratios. This study presents a new geochronologic tool for timing the emplacement of kimberlites using the (U-Th)/He system to date the cooling of common kimberlite phenocrystic and xenocrystic minerals. To demonstrate the use of this technique, new apatite, titanite, zircon, magnetite and garnet (U-Th)/He ages constrain the timing of emplacement for the Stockdale, Tuttle, Baldwin Creek, Bala, and Leonardville kimberlite pipes, located in Riley County, Kansas. Zircon from the Tuttle pipe and titanite from the Stockdale pipe yield (U-Th)/He ages of 108.6 ?? 9.6??Ma and 106.4 ?? 3.1??Ma, respectively. These data are consistent with new Tuttle kimberlite Rb-Sr analyses of phlogopite megacrysts that give a five point isochron age of 106.6 ?? 1.0??Ma. Similarly, an apatite (U-Th)/He age of 85.3 ?? 2.3??Ma from the Baldwin Creek kimberlite is in agreement with a Rb-Sr phlogopite age of 88.4 ?? 2.7??Ma. These dates demonstrate that (U-Th)/He thermochronometry provides reliable timing constraints on the cooling of common kimberlite xenocrystic phases, thereby timing kimberlite emplacement. In addition to the use of more commonly used apatite and zircon (U-Th)/He thermochronometers, we produced reliable emplacement ages of 103.0 ?? 7.5??Ma for the Bala kimberlite using (U-Th)/He dating of phenocrystic magnetite and an age of 98.8 ?? 8.9??Ma for the Tuttle kimberlite using (U-Th)/He dating of megacrystic garnet. In contrast, kimberlitic apatite (U-Th)/He ages from the Stockdale, Bala, Tuttle, and Leonardville kimberlites yield ages ranging from 67.3 ?? 4.4??Ma to 64.3 ?? 5.6??Ma, suggesting a local, possibly hydrothermal reheating event resulting in resetting of the apatite (U-Th)/He clock in latest Cretaceous to earliest Tertiary time. Additional (U-Th)/He analyses of apatite from nearby sandstone and basement rocks suggest regional cooling below ~ 70????C at ~ 165??Ma. These (U-Th)/He and Rb-Sr age data imply that the kimberlites were emplaced over a period of time from ~ 85-110??Ma with several pipes subjected to local reheating at ~ 65??Ma. ?? 2008 Elsevier B.V. All rights reserved.

Earth and Planetary Science Letters↗

USGS2014

Three-dimensional imaging, change detection, and stability assessment during the centerline trench levee seepage experiment using terrestrial light detection and ranging technology, Twitchell Island, California, 2012

A full scale field seepage test was conducted on a north-south trending levee segment of a now bypassed old meander belt on Twitchell Island, California, to understand the effects of live and decaying root systems on levee seepage and slope stability. The field test in May 2012 was centered on a north-south trench with two segments: a shorter control segment and a longer seepage test segment. The complete length of the trench area measured 40.4 meters (m) near the levee centerline with mature trees located on the waterside and landside of the levee flanks. The levee was instrumented with piezometers and tensiometers to measure positive and negative porewater pressures across the levee after the trench was flooded with water and held at a constant hydraulic head during the seepage test—the results from this component of the experiment are not discussed in this report. We collected more than one billion three-dimensional light detection and ranging (lidar) data points before, during, and after the centerline seepage test to assess centimeter-scale stability of the two trees and the levee crown. During the seepage test, the waterside tree toppled (rotated 20.7 degrees) into the water. The landside tree rotated away from the levee by 5 centimeters (cm) at a height of 2 m on the tree. The paved surface of the levee crown had three regions that showed subsidence on the waterside of the trench—discussed as the northern, central, and southern features. The northern feature is an elongate region that subsided 2.1 cm over an area with an average width of 1.35 m that extends 15.8 m parallel to the trench from the northern end of the trench to just north of the trench midpoint, and is associated with a crack 1 cm in height that formed during the seepage test on the trench wall. The central subsidence feature is a semicircular region on the waterside of the trench that subsided by as much as 6.2 cm over an area 3.4 m wide and 11.2 m long. The southern feature is an elongate region that has a maximum subsidence of 3.5 cm over an area 0.75 m wide and 8.1 m long and is associated with a number of small fractures in the pavement that are predominately north-south-trending and parallel to the trench. We determined that there was no significant motion of the levee flank during the last week of the seepage test. We also determined biomorphic parameters for the landside tree, such as the 3D positioning on the levee, tree height, levee parallel/perpendicular cross sectional area, and canopy centroid. These biomorphic parameters were requested to support a University of California Berkeley team studying seepage and stability on the levee. A gridded, 2-cm bare-earth digital elevation model of the levee crown and the landside levee flank from the final terrestrial lidar (T-Lidar) survey provided detailed topographic data for future assessment. Because the T-Lidar was not integrated into the project design, other than an initial courtesy dataset to help characterize the levee surface, our ability to contribute to the overall science goals of the seepage test was limited. Therefore, our analysis focused on developing data collection and processing methodology necessary to align ultra high-resolution T-Lidar data (with an average spot spacing 2–3 millimeters on the levee crown) from several instrument setup locations to detect, measure, and characterize dynamic centimeter-scale deformation and surface changes during the seepage test.

California↗

USGS2008

EAARL coastal topography - Northern Gulf of Mexico

These remotely sensed, geographically referenced elevation measurements of Lidar-derived coastal topography were produced as a collaborative effort between the U.S. Geological Survey (USGS), Florida Integrated Science Center (FISC), St. Petersburg, FL and the National Aeronautics and Space Administration (NASA), Wallops Flight Facility, VA. One objective of this research is to create techniques to survey areas for the purposes of geomorphic change studies following major storm events. The USGS Coastal and Marine Geology Program's National Assessment of Coastal Change Hazards project is a multi-year undertaking to identify and quantify the vulnerability of U.S. shorelines to coastal change hazards such as effects of severe storms, sea-level rise, and shoreline erosion and retreat. Airborne Lidar surveys conducted during periods of calm weather are compared to surveys collected following extreme storms in order to quantify the resulting coastal change. Other applications of high-resolution topography include habitat mapping, ecological monitoring, volumetric change detection, and event assessment. The purpose of this project is to provide highly detailed and accurate datasets of the northern Gulf of Mexico coastal areas, acquired on September 19, 2004, immediately following Hurricane Ivan. The datasets are made available for use as a management tool to research scientists and natural resource managers. An innovative airborne Lidar instrument originally developed at the NASA Wallops Flight Facility, and known as the Experimental Airborne Advanced Research Lidar (EAARL), was used during data acquisition. The EAARL system is a raster-scanning, waveform-resolving, green-wavelength (532 nanometer) Lidar designed to map near-shore bathymetry, topography, and vegetation structure simultaneously. The EAARL sensor suite includes the raster-scanning, water-penetrating full-waveform adaptive Lidar, a down-looking RGB (red-green-blue) digital camera, a high-resolution multi-spectral color infrared (CIR) camera, two precision dual-frequency kinematic carrier-phase GPS receivers and an integrated miniature digital inertial measurement unit which provide for sub-meter georeferencing of each laser sample. The nominal EAARL platform is a twin-engine Cessna 310 aircraft, but the instrument may be deployed on a range of light aircraft. A single pilot, a Lidar operator, and a data analyst constitute the crew for most survey operations. This sensor has the potential to make significant contributions in measuring sub-aerial and submarine coastal topography within cross-environmental surveys. Elevation measurements were collected over the survey area using the EAARL system on September 19, 2004. The survey resulted in the acquisition of 3.2 gigabytes of data. The data were processed using the Airborne Lidar Processing System (ALPS), a custom-built processing system developed in a NASA-USGS collaboration. ALPS supports the exploration and processing of Lidar data in an interactive or batch mode. Modules for pre-survey flight line definition, flight path plotting, Lidar raster and waveform investigation, and digital camera image playback have been developed. Processing algorithms have been developed to extract the range to the first and last significant return within each waveform. ALPS is routinely used to create maps that represent submerged or sub-aerial topography. Specialized filtering algorithms have been implemented to determine the 'bare earth' under vegetation from a point cloud of 'last return' elevations.

Alabama, Florida, Mississippi↗

USGS2022

Forecasting explosions at Sinabung Volcano, Indonesia, based on SO2 emission rates

Dome-building volcanic eruptions are often associated with frequent Vulcanian explosions, which constitute a substantial threat to proximal communities. One proposed mechanism driving such explosions is the sealing of the shallow volcanic system followed by pressurization due to gas accumulation beneath the seal. We investigate this hypothesis at Sinabung Volcano (Sumatra, Indonesia), which has been in a state of eruption since August 2010. In 2013, the volcano began erupting a lava dome and lava flow, and frequent explosions produced eruptive columns that rose many kilometers into the atmosphere and at times sent pyroclastic density currents down the southeast flanks. A network of scanning Differential Optical Absorption Spectrometers (DOAS) was installed on the volcano’s eastern flank in 2016 to continuously monitor SO 2 emission rates during daytime hours. Analysis of the DOAS data from October 2016 to September 2017 revealed that passive SO 2 emissions were generally lower in the 5 days leading up to explosive events (∼100 t/d) than was common in 5-day periods leading up to days on which no explosions occurred (∼200 t/d). The variability of passive SO 2 emissions, expressed as the standard deviation, also took on a slightly wider range of values before days with explosions (0–103 t/d at 1-sigma) than before days without explosions (43–117 t/d). These observations are consistent with the aforementioned seal-failure model, where the sealing of the volcanic conduit blocks gas emissions and leads to pressurization and potential Vulcanian explosions. We develop a forecasting methodology that allows calculation of a relative daily explosion probability based solely on measurements of the SO 2 emission rate in the preceding days. We then calculate forecast explosion probabilities for the remaining SO 2 emissions dataset (October 2017—September 2021). While the absolute accuracy of forecast explosion probabilities is variable, the method can inform the probability of an explosion occurring relative to that on other days in each test period. This information can be used operationally by volcano observatories to assess relative risk. The SO 2 emissions-based forecasting method is likely applicable to other open vent volcanoes experiencing dome-forming eruptions.

Frontiers in Earth Science↗