SEARCH · USGS Science

Results for “The Earth Observer”

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 1,243 records · Page 69Linked to original sources

USGS2007

Quaternarymagmatism in the Cascades - Geologic perspectives

Foreward The Cascade magmatic arc is a belt of Quaternary volcanoes that extends 1,250 km from Lassen Peak in northern California to Meager Mountain in Canada, above the subduction zone where the Juan de Fuca Plate plunges beneath the North American Plate. This Professional Paper presents a synthesis of the entire volcanic arc, addressing all 2,300 known Quaternary volcanoes, not just the 30 or so visually prominent peaks that comprise the volcanic skyline. Study of Cascade volcanoes goes back to the geological explorers of the late 19th century and the seminal investigations of Howel Williams in the 1920s and 1930s. However, major progress and application of modern scientific methods and instrumentation began only in the 1970s with the advent of systematic geological, geophysical, and geochemical studies of the entire arc. Initial stimulus from the USGS Geothermal Research Program was enhanced by the USGS Volcano Hazards Program following the 1980 eruption of Mount St. Helens. Together, these two USGS Programs have provided more than three decades of stable funding, staffing, and analytical support. This Professional Paper summarizes the resultant USGS data sets and integrates them with the parallel contributions of other investigators. The product is based upon an all-encompassing and definitive geological database, including chemical and isotopic analyses to characterize the rocks and geochronology to provide the critical time constraints. Until now, this massive amount of data has not been summarized, and a systematic and uniform interpretation firmly grounded in geological fact has been lacking. Herein lies the primary utility of this Cascade volume. It not only will be the mandatory starting point for new workers, but also will provide essential geological context to broaden the perspectives of current investigators of specific Cascade volcanoes. Wes Hildreth's insightful understanding of volcanic processes and his uncompromising scientific integrity make him uniquely qualified to present this synthesis. During more than three decades of volcanological studies, he has carried out comprehensive investigations of Mount Adams, Mount Baker, the Three Sisters, and the Simcoe Mountains Volcanic Field. He also brings a broad experience in other volcanic arcs, having conducted integrated field and laboratory investigations at several major volcanic centers in the Andes and the Aleutian arcs. His expertise and perspective have been further enhanced by in-depth petrologic studies of caldera environments, primarily in Long Valley, California, and Yellowstone. On the basis of all these field and laboratory investigations and exhaustive literature searches, he has published three definitive petrologic syntheses addressing the passage and transformation of basaltic magmas from their mantle sources through the crust to form the many types of volcanic manifestations at the Earth's surface. A major strength of this Professional Paper is that it adheres to data first and foremost, and only then correlates these data with relevant theories. Petrological and geophysical interpretation is left to the later sections of the volume, and even there is never allowed to stray from the pertinent databases. Hildreth's interpretations are not just idle speculations, but are carefully reasoned inferences firmly based on his thorough evaluation of the observational geological data. Professional Paper 1744 should not be skimmed lightly, in the hope that the salient points will quickly rub off. Instead, every section, indeed every paragraph, presents scholarly observations and insightful interpretations that demand careful and thoughtful study. This volume will influence and guide the course of Cascade investigations for decades to come.

British Columbia, California, Oregon, Washington↗

USGS2021

The Independent Volcanic Eruption Source Parameter Archive (IVESPA, version 1.0): A new observational database to support explosive eruptive column model validation and development

Eruptive column models are powerful tools for investigating the transport of volcanic gas and ash, reconstructing past explosive eruptions, and simulating future hazards. However, the evaluation of these models is challenging as it requires independent estimates of the main model inputs (e.g. mass eruption rate) and outputs (e.g. column height). There exists no database of independently estimated eruption source parameters (ESPs) that is extensive, standardized, maintained, and consensus-based. This paper introduces the Independent Volcanic Eruption Source Parameter Archive (IVESPA, ivespa.co.uk), a community effort endorsed by the International Association of Volcanology and Chemistry of the Earth’s Interior (IAVCEI) Commission on Tephra Hazard Modelling. We compiled data for 134 explosive eruptive events, spanning the 1902-2016 period, with independent estimates of: i) total erupted mass of fall deposits; ii) duration; iii) eruption column height; and iv) atmospheric conditions. Crucially, we distinguish plume top versus umbrella spreading height, and the height of ash versus sulphur dioxide injection. All parameter values provided have been vetted independently by at least two experts. Uncertainties are quantified systematically, including flags to describe the degree of interpretation of the literature required for each estimate. IVESPA also includes a range of additional parameters such as total grain size distribution, eruption style, morphology of the plume (weak versus strong), and mass contribution from pyroclastic density currents, where available. We discuss the future developments and potential applications of IVESPA and make recommendations for reporting ESPs to maximize their usability across different applications. IVESPA covers an unprecedented range of ESPs and can therefore be used to evaluate and develop eruptive column models across a wide range of conditions using a standardized dataset.

Journal of Volcanology and Geothermal Research↗

USGS2016

A synthesis of the basal thermal state of the Greenland Ice Sheet

The basal thermal state of an ice sheet (frozen or thawed) is an important control upon its evolution, dynamics and response to external forcings. However, this state can only be observed directly within sparse boreholes or inferred conclusively from the presence of subglacial lakes. Here we synthesize spatially extensive inferences of the basal thermal state of the Greenland Ice Sheet to better constrain this state. Existing inferences include outputs from the eight thermomechanical ice-flow models included in the SeaRISE effort. New remote-sensing inferences of the basal thermal state are derived from Holocene radiostratigraphy, modern surface velocity and MODIS imagery. Both thermomechanical modeling and remote inferences generally agree that the Northeast Greenland Ice Stream and large portions of the southwestern ice-drainage systems are thawed at the bed, whereas the bed beneath the central ice divides, particularly their west-facing slopes, is frozen. Elsewhere, there is poor agreement regarding the basal thermal state. Both models and remote inferences rarely represent the borehole-observed basal thermal state accurately near NorthGRIP and DYE-3. This synthesis identifies a large portion of the Greenland Ice Sheet (about one third by area) where additional observations would most improve knowledge of its overall basal thermal state.

Journal of Geophysical Research F: Earth Surface↗

USGS2016

Removing sun glint from optical remote sensing images of shallow rivers

Sun glint is the specular reflection of light from the water surface, which often causes unusually bright pixel values that can dominate fluvial remote sensing imagery and obscure the water-leaving radiance signal of interest for mapping bathymetry, bottom type, or water column optical characteristics. Although sun glint is ubiquitous in fluvial remote sensing imagery, river-specific methods for removing sun glint are not yet available. We show that existing sun glint-removal methods developed for multispectral images of marine shallow water environments over-correct shallow portions of fluvial remote sensing imagery resulting in regions of unreliable data along channel margins. We build on existing marine glint-removal methods to develop a river-specific technique that removes sun glint from shallow areas of the channel without overcorrection by accounting for non-negligible water-leaving near-infrared radiance. This new sun glint-removal method can improve the accuracy of spectrally-based depth retrieval in cases where sun glint dominates the at-sensor radiance. For an example image of the gravel-bed Snake River, Wyoming, USA, observed-vs.-predicted R 2 values for depth retrieval improved from 0.66 to 0.76 following sun glint removal. The methodology presented here is straightforward to implement and could be incorporated into image processing workflows for multispectral images that include a near-infrared band.

Earth Surface Processes and Landforms↗

USGS2003

Transient rheology of the uppermost mantle beneath the Mojave Desert, California

Geodetic data indicate that the M7.1 Hector Mine, California, earthquake was followed by a brief period (a few weeks) of rapid deformation preceding a prolonged phase of slower deformation. We find that the signal contained in continuous and campaign global positioning system data for 2.5 years after the earthquake may be explained with a transient rheology. Quantitative modeling of these data with allowance for transient (linear biviscous) rheology in the lower crust and upper mantle demonstrates that transient rheology in the upper mantle is dominant, its material properties being typified by two characteristic relaxation times ???0.07 and ???2 years. The inferred mantle rheology is a Jeffreys solid in which the transient and steady-state shear moduli are equal. Consideration of a simpler viscoelastic model with a linear univiscous rheology (2 fewer parameters than a biviscous model) shows that it consistently underpredicts the amplitude of the first ???3 months signal, and allowance for a biviscous rheology is significant at the 99.0% confidence level. Another alternative model - deep postseismic afterslip beneath the coseismic rupture - predicts a vertical velocity pattern opposite to the observed pattern at all time periods considered. Despite its plausibility, the advocated biviscous rheology model is non-unique and should be regarded as a viable alternative to the non-linear mantle rheology model for governing postseismic flow beneath the Mojave Desert. Published by Elsevier B.V.

Earth and Planetary Science Letters↗

USGS2009

Microseismicity at the North Anatolian Fault in the Sea of Marmara offshore Istanbul, NW Turkey

The North Anatolian Fault Zone (NAFZ) below the Sea of Marmara forms a “seismic gap” where a major earthquake is expected to occur in the near future. This segment of the fault lies between the 1912 Ganos and 1999 İzmit ruptures and is the only NAFZ segment that has not ruptured since 1766. To monitor the microseismic activity at the main fault branch offshore of Istanbul below the Çınarcık Basin, a permanent seismic array (PIRES) was installed on the two outermost Prince Islands, Yassiada and Sivriada, at a few kilometers distance to the fault. In addition, a temporary network of ocean bottom seismometers was deployed throughout the Çınarcık Basin. Slowness vectors are determined combining waveform cross correlation and P wave polarization. We jointly invert azimuth and traveltime observations for hypocenter determination and apply a bootstrap resampling technique to quantify the location precision. We observe seismicity rates of 20 events per month for M < 2.5 along the basin. The spatial distribution of hypocenters suggests that the two major fault branches bounding the depocenter below the Çınarcık Basin merge to one single master fault below ∼17 km depth. On the basis of a cross-correlation technique we group closely spaced earthquakes and determine composite focal mechanisms implementing recordings of surrounding permanent land stations. Fault plane solutions have a predominant right-lateral strike-slip mechanism, indicating that normal faulting along this part of the NAFZ plays a minor role. Toward the west we observe increasing components of thrust faulting. This supports the model of NW trending, dextral strike-slip motion along the northern and main branch of the NAFZ below the eastern Sea of Marmara.

Istanbul↗

USGS2012

Geodetic analysis of reservoir depletion at the Geyser steam field in northern California

Reservoir depletion at the Geysers from 1974 to 1977 is evident in measured changes in gravity, surface strain, and pore pressure drainage. The drainage area increased about 20%, the maximum gravity decrease was about −120 μGal, and the maximum elevation change was about 6 cm during this period. Since the net mass withdrawal is known, it may be combined with the gravity change to estimate a drainage volume. The maximum drainage volume is 25 km 3 . Because the depth of the caprock and extent of the pore pressure drainage is known, this volume limit implies that no significant mass withdrawal occurs below a depth of 4 km. The ratios of surface elevation changes to horizontal contraction coinciding with the drainage area imply an equant drainage geometry, assuming that fluid production produces negative dilatation. Using the same cylindrical geometry used to model the gravity, negative dilatational strain rates of 4 to 5 × 10 −5 /yr were found to produce the observed surface displacements. The likelihood of boiling in the system to produce steam combined with apparent large bulk moduli (from seismic velocities) and small pore pressure declines suggest that most of the reservoir contraction is due to cooling.

Journal of Geophysical Research Solid Earth↗

USGS2022

Diverse tsunamigenesis triggered by the Hunga Tonga-Hunga Ha’apai eruption

On the evening of 15 January 2022, the Hunga Tonga-Hunga Ha’apai volcano 1 unleashed a violent underwater eruption, blanketing the surrounding land masses in ash and debris . The eruption generated tsunamis observed around the world. An event of this type last occurred in 1883 during the eruption of Krakatau , and thus we have the first observations of a tsunami from a large emergent volcanic eruption captured with modern instrumentation. Here we show that the explosive eruption generated waves through multiple mechanisms, including: (1) air–sea coupling with the initial and powerful shock wave radiating out from the explosion in the immediate vicinity of the eruption; (2) collapse of the water cavity created by the underwater explosion; and (3) air–sea coupling with the air-pressure pulse that circled the Earth several times, leading to a global tsunami. In the near field, tsunami impacts are strongly controlled by the water-cavity source whereas the far-field tsunami, which was unusually persistent, can be largely described by the air-pressure pulse mechanism. Catastrophic damage in some harbours in the far field was averted by just tens of centimetres, implying that a modest sea level rise combined with a future, similar event would lead to a step-function increase in impacts on infrastructure. Piecing together the complexity of this event has broad implications for coastal hazards in similar geophysical settings, suggesting a currently neglected source of global tsunamis.

Nature↗

USGS2010

Upper mantle rheology from GRACE and GPS postseismic deformation after the 2004 Sumatra-Andaman earthquake

Mantle rheology is one of the essential, yet least understood, material properties of our planet, controlling the dynamic processes inside the Earth's mantle and the Earth's response to various forces. With the advent of GRACE satellite gravity, measurements of mass displacements associated with many processes are now available. In the case of mass displacements related to postseismic deformation, these data may provide new constraints on the mantle rheology. We consider the postseismic deformation due to the M w = 9.2 Sumatra 26 December 2004 and M w = 8.7 Nias 28 March 2005 earthquakes. Applying wavelet analyses to enhance those local signals in the GRACE time varying geoids up to September 2007, we detect a clear postseismic gravity signal. We supplement these gravity variations with GPS measurements of postseismic crustal displacements to constrain postseismic relaxation processes throughout the upper mantle. The observed GPS displacements and gravity variations are well explained by a model of viscoelastic relaxation plus a small amount of afterslip at the downdip extension of the coseismically ruptured fault planes. Our model uses a 60 km thick elastic layer above a viscoelastic asthenosphere with Burgers body rheology. The mantle below depth 220 km has a Maxwell rheology. Assuming a low transient viscosity in the 60–220 km depth range, the GRACE data are best explained by a constant steady state viscosity throughout the ductile portion of the upper mantle (e.g., 60–660 km). This suggests that the localization of relatively low viscosity in the asthenosphere is chiefly in the transient viscosity rather than the steady state viscosity. We find a 8.10 18 Pa s mantle viscosity in the 220–660 km depth range. This may indicate a transient response of the upper mantle to the high amount of stress released by the earthquakes. To fit the remaining misfit to the GRACE data, larger at the smaller spatial scales, cumulative afterslip of about 75 cm at depth should be added over the period spanned by the GRACE models. It produces only small crustal displacements. Our results confirm that satellite gravity data are an essential complement to ground geodetic and geophysical networks in order to understand the seismic cycle and the Earth's inner structure.

Andaman Sea, Sumatra↗

USGS2010

Southern high latitude dune fields on Mars: Morphology, aeolian inactivity, and climate change

In a study area spanning the martian surface poleward of 50?? S., 1190 dune fields have been identified, mapped, and categorized based on dune field morphology. Dune fields in the study area span ??? 116400km2, leading to a global dune field coverage estimate of ???904000km2, far less than that found on Earth. Based on distinct morphological features, the dune fields were grouped into six different classes that vary in interpreted aeolian activity level from potentially active to relatively inactive and eroding. The six dune field classes occur in specific latitude zones, with a sequence of reduced activity and degradation progressing poleward. In particular, the first signs of stabilization appear at ???60?? S., which broadly corresponds to the edge of high concentrations of water-equivalent hydrogen content (observed by the Neutron Spectrometer) that have been interpreted as ground ice. This near-surface ground ice likely acts to reduce sand availability in the present climate state on Mars, stabilizing high latitude dunes and allowing erosional processes to change their morphology. As a result, climatic changes in the content of near-surface ground ice are likely to influence the level of dune activity. Spatial variation of dune field classes with longitude is significant, suggesting that local conditions play a major role in determining dune field activity level. Dune fields on the south polar layered terrain, for example, appear either potentially active or inactive, indicating that at least two generations of dune building have occurred on this surface. Many dune fields show signs of degradation mixed with crisp-brinked dunes, also suggesting that more than one generation of dune building has occurred since they originally formed. Dune fields superposed on early and late Amazonian surfaces provide potential upper age limits of ???100My on the south polar layered deposits and ???3Ga elsewhere at high latitudes. No craters are present on any identifiable dune fields, which can provide a lower age limit through crater counting: assuming all relatively stabilized dune fields represent a single noncontiguous surface of uniform age, their estimated crater retention age is <???10000years. An average-sized uncratered dune field (94km2) has a crater retention age <???8My. This apparent youth suggests that present-day climate conditions are responsible for the observed degradation and reduced level of aeolian activity. A lack of observed transport pathways and the absence of large dune fields in the largest basins (Hellas and Argyre Planitiae) are consistent with the previously proposed idea that dune sands are not typically transported far from their source regions on Mars. ?? 2009 Elsevier B.V.

Geomorphology↗

USGS2007

U/Th series radionuclides as coastal groundwater tracers

The study of coastal groundwater has recently surfaced as an active interdisciplinary area of research, driven foremost by its importance as a poorly quantified pathway for subsurface material transport into coastal ecosystems. Key issue in coastal groundwater research include a complete geochemical characterization of the groundwater(s); quantification of the kinetics of subsurface transport, including rock-water interactions; determination of groundwater ages; tracing of groundwater discharge into coastal waters using radiochemical fingerprints; and an assessment of the potential ecological impact of such subsurface flow to a reviving water body. For such applications, the isotopic systemics of select naturally occurring radionucludes in the U/Th series has proven to be particularly useful. These radionuclides (e.g., U, Th, Ram and Rn) are ubiquitous in all groundwaters ad are represented by several isotopes with widely different half-lives and chemistries (Figure 1). As a result, varied biogeochemical processes occurring over a broad range of time scales can be studied. In source rock, most U/Th series isotopes in secular equilibrium; that is, the rate of decay of a daughter isotope is equal to that of it radiogenic parent, and so will have equal activities (in this context, the specific activity is simply a measure of the amount of radioactivity per unit amount). In contrast, these nuclides exhibit strong fractionations within the surrounding groundwaters because of their respective physiochemical differences. Disequilibria in U/Th series radionuclides can thus be used to identify distinct water masses, quantify release rates from source rocks, assess groundwater migration rates, and assess groundwater discharge rates in coastal waters., Large isotopic variations also have the potential for providing precise fingerprints for groundwaters from specific aquifers and have been explored as a means for calculating groundwater ages and estuarine water mass transit times. The highly fractionated nature of U/Th series nuclides in groundwater is illustrated by the range in some measured activities. highest activities are typically observed for 222 Rn, reflecting the inert nature of this noble gas. Groundwater 222 Rn (t 1/2 =3.8) activities are thus controlled only by rapid in situ decay (Table 1) and production within host rocks, without the added complications of reversible removal via absorption or precipitation. Uranium, which is soluble as U(VI) in oxidizing waters, is present in intermediate activities in groundwaters that are moderated by redox-initiated removal onto aquifer rocks. The alkaline earth Ra and, to a greater extent, the less soluble actinide Th are readily removed from groundwater by water -- rock interactions and so are strongly depleted. Both of these elements have very short-lived as well as longer-lived isotopes, and so isotopes compositions reflect processes over a range of time scales. Many studies have evaluated and behavior of select radionuclides in groundwater and surface water systems. Recent advances in high-=precision mass spectrometry have opted new possibilities for more subtle interpretations in select long-lived U/Th series isotopes, such as U, Ra, Pa, and Th. However, these techniques have yet to be fully developed, ahns as a consequence, such data remain largely scarce and underutilized. Although many different approaches have been developed to study radionucluide behavior in groundwater, all are based on principles of radioactive production and decay and knowledge of source terms from weathering and recoil processes, as well as removal terms from the interaction with aquifer host rock surface by sorption and precipitation. This review is structured to present first a brief description of the background, driving forces, scales, and ecological significance of submarine groundwater discharge. Following this, a description of the geochemistry and behavior of select radionuclides in groundwater will be presented, and their application to tracing submarine groundwater discharge will be discussed.

Chemical Reviews↗

USGS2023

Investigating permafrost carbon dynamics in Alaska with artificial intelligence

Positive feedbacks between permafrost degradation and the release of soil carbon into the atmosphere impact land–atmosphere interactions, disrupt the global carbon cycle, and accelerate climate change. The widespread distribution of thawing permafrost is causing a cascade of geophysical and biochemical disturbances with global impacts. Currently, few earth system models account for permafrost carbon feedback (PCF) mechanisms. This research study integrates artificial intelligence (AI) tools and information derived from field-scale surveys across the tundra and boreal landscapes in Alaska. We identify and interpret the permafrost carbon cycling links and feedback sensitivities with GeoCryoAI, a hybridized multimodal deep learning (DL) architecture of stacked convolutionally layered, memory-encoded recurrent neural networks (NN). This framework integrates in-situ measurements and flux tower observations for teacher forcing and model training. Preliminary experiments to quantify, validate, and forecast permafrost degradation and carbon efflux across Alaska demonstrate the fidelity of this data-driven architecture. More specifically, GeoCryoAI logs the ecological memory and effectively learns covariate dynamics while demonstrating an aptitude to simulate and forecast PCF dynamics—active layer thickness (ALT), carbon dioxide flux (CO 2 ), and methane flux (CH 4 )—with high precision and minimal loss (i.e. ALT RMSE : 1.327 cm [1969–2022]; CO 2 RMSE : 0.697 µ molCO 2 m −2 s −1 [2003–2021]; CH 4 RMSE : 0.715 nmolCH 4 m −2 s −1 [2011–2022]). ALT variability is a sensitive harbinger of change, a unique signal characterizing the PCF, and our model is the first characterization of these dynamics across space and time.

Alaska↗

USGS2023

First observation of the ground-state electron-capture of 40K

Potassium-40 is a widespread, naturally occurring isotope whose radioactivity impacts estimated geological ages spanning billions of years, nuclear structure theory, and subatomic rare-event searches—including those for dark matter and neutrinoless double-beta decay. The decays of this long-lived isotope must be precisely known for its use as a geochronometer, and to account for its presence in low-background experiments. There are several known decay modes for potassium-40, but a predicted electron-capture decay directly to the ground state of argon-40 has never been observed. The existence of this decay mode impacts several fields, while theoretical predictions span an order of magnitude. Here we report on the first, successful observation of this rare decay mode, obtained by the KDK (potassium decay) Collaboration using a novel combination of a low-threshold x-ray detector surrounded by a tonne-scale, high-efficiency γ -ray tagger at Oak Ridge National Laboratory. A blinded analysis reveals a distinctly nonzero ratio of intensities of ground-state electron-captures ( I EC 0 ) over excited-state ones ( I EC ∗ ) of I EC 0 / I EC ∗ = 0.0095 stat ± 0.0022 sys ± 0.0010 (68% CL), with the null hypothesis rejected at 4 σ [Stukel et al. , Phys. Rev. Lett. 131 , 052503 (2023) ]. In terms of branching ratio, this unambiguous signal yields I EC 0 = 0.098 % stat ± 0.023 % sys ± 0.010 % , roughly half of the commonly used prediction. This first observation of a third-forbidden unique electron capture improves our understanding of low-energy backgrounds in dark-matter searches and has implications for nuclear-structure calculations. For example, a shell-model based theoretical estimate for the neutrinoless double-beta decay half-life of calcium-48 is increased by a factor of 7 + 3 − 2 . Our nonzero measurement shifts geochronological ages by up to a percent; implications are illustrated for Earth and solar system chronologies.

Physical Review C↗

USGS1996

Application of ground-penetrating radar to investigation of near-surface fault properties in the San Francisco Bay region

In many geologic environments, ground-penetrating radar (GPR) provides high-resolution images of near-surface Earth structure. GPR data collection is nondestructive and very economical. The scale of features detected by GPR lies between those imaged by high-resolution seismic reflection surveys and those exposed in trenches and is therefore potentially complementary to traditional techniques for fault location and mapping. Sixty-two GPR profiles were collected at 12 sites in the San Francisco Bay region. Results show that GPR data correlate with large-scale features in existing trench observations, can be used to locate faults where they are buried or where their positions are not well known, and can identify previously unknown fault segments. The best data acquired were on a profile across the San Andreas fault, traversing Pleistocene terrace deposits south of Olema in Marin County; this profile shows a complicated multi-branched fault system from the ground surface down to about 40 m, the maximum depth for which data were recorded.

California↗

USGS2022

Medium spatial resolution mapping of global land cover and land cover change across multiple decades from Landsat

Land cover maps are essential for characterizing the biophysical properties of the Earth’s land areas. Because land cover information synthesizes a rich array of information related to both the ecological condition of land areas and their exploitation by humans, they are widely used for basic and applied research that requires information related to land surface properties (e.g., terrestrial carbon models, water balance models, weather, and climate models) and are core inputs to models and analyses used by natural resource scientists and land managers. As the Earth’s global population has grown over the last several decades rates of land cover change have increased dramatically, with enormous impacts on ecosystem services (e.g., biodiversity, water supply, carbon sequestration, etc.). Hence, accurate information related to land cover is essential for both managing natural resources and for understanding society’s ecological, biophysical, and resource management footprint. To address the need for high-quality land cover information we are using the global record of Landsat observations to compile annual maps of global land cover from 2001 to 2020 at 30 m spatial resolution. To create these maps we use features derived from time series of Landsat imagery in combination with ancillary geospatial data and a large database of training sites to classify land cover at annual time step. The algorithm that we apply uses temporal segmentation to identify periods with stable land cover that are separated by breakpoints in the time series. Here we provide an overview of the methods and data sets we are using to create global maps of land cover. We describe the algorithms used to create these maps and the core land cover data sets that we are creating through this effort, and we summarize our approach to accuracy assessment. We also present a synthesis of early results and discuss the strengths and weaknesses of our early map products and the challenges that we have encountered in creating global land cover data sets from Landsat. Initial accuracy assessment for North America shows good overall accuracy (77.0 ± 2.0% correctly classified) and 79.8% agreement with the European Space Agency (ESA) WorldCover product. The land cover mapping results we report provide the foundation for robust, repeatable, and accurate mapping of global land cover and land cover change across multiple decades at 30 m spatial resolution from Landsat.

Frontiers in Remote Sensing↗

USGS2017

Reviews and syntheses: Field data to benchmark the carbon cycle models for tropical forests

For more accurate projections of both the global carbon (C) cycle and the changing climate, a critical current need is to improve the representation of tropical forests in Earth system models. Tropical forests exchange more C, energy, and water with the atmosphere than any other class of land ecosystems. Further, tropical-forest C cycling is likely responding to the rapid global warming, intensifying water stress, and increasing atmospheric CO 2 levels. Projections of the future C balance of the tropics vary widely among global models. A current effort of the modeling community, the ILAMB (International Land Model Benchmarking) project, is to compile robust observations that can be used to improve the accuracy and realism of the land models for all major biomes. Our goal with this paper is to identify field observations of tropical-forest ecosystem C stocks and fluxes, and of their long-term trends and climatic and CO 2 sensitivities, that can serve this effort. We propose criteria for reference-level field data from this biome and present a set of documented examples from old-growth lowland tropical forests. We offer these as a starting point towards the goal of a regularly updated consensus set of benchmark field observations of C cycling in tropical forests.

Biogeosciences↗

USGS2014

Time-averaged discharge rate of subaerial lava at Kīlauea Volcano, Hawai‘i, measured from TanDEM-X interferometry: Implications for magma supply and storage during 2011-2013

Differencing digital elevation models (DEMs) derived from TerraSAR add-on for Digital Elevation Measurements (TanDEM-X) synthetic aperture radar imagery provides a measurement of elevation change over time. On the East Rift Zone (EZR) of Kīlauea Volcano, Hawai‘i, the effusion of lava causes changes in topography. When these elevation changes are summed over the area of an active lava flow, it is possible to quantify the volume of lava emplaced at the surface during the time spanned by the TanDEM-X data—a parameter that can be difficult to measure across the entirety of an ~100 km 2 lava flow field using ground-based techniques or optical remote sensing data. Based on the differences between multiple TanDEM-X-derived DEMs collected days to weeks apart, the mean dense-rock equivalent time-averaged discharge rate of lava at Kīlauea between mid-2011 and mid-2013 was approximately 2 m 3 /s, which is about half the long-term average rate over the course of Kīlauea's 1983–present ERZ eruption. This result implies that there was an increase in the proportion of lava stored versus erupted, a decrease in the rate of magma supply to the volcano, or some combination of both during this time period. In addition to constraining the time-averaged discharge rate of lava and the rates of magma supply and storage, topographic change maps derived from space-based TanDEM-X data provide insights into the four-dimensional evolution of Kīlauea's ERZ lava flow field. TanDEM-X data are a valuable complement to other space-, air-, and ground-based observations of eruptive activity at Kīlauea and offer great promise at locations around the world for aiding with monitoring not just volcanic eruptions but any hazardous activity that results in surface change, including landslides, floods, earthquakes, and other natural and anthropogenic processes.

Hawai'i↗

USGS2011

Bedform response to flow variability

Laboratory observations and computational results for the response of bedform fields to rapid variations in discharge are compared and discussed. The simple case considered here begins with a relatively low discharge over a flat bed on which bedforms are initiated, followed by a short high-flow period with double the original discharge, during which the morphology of the bedforms adjusts, followed in turn by a relatively long period of the original low discharge. For the grain size and hydraulic conditions selected, the Froude number remains subcritical during the experiment, and sediment moves predominantly as bedload. Observations show rapid development of quasi-two-dimensional bedforms during the initial period of low flow with increasing wavelength and height over the initial low-flow period. When the flow increases, the bedforms rapidly increase in wavelength and height, as expected from other empirical results. When the flow decreases back to the original discharge, the height of the bedforms quickly decreases in response, but the wavelength decreases much more slowly. Computational results using an unsteady two-dimensional flow model coupled to a disequilibrium bedload transport model for the same conditions simulate the formation and initial growth of the bedforms fairly accurately and also predict an increase in dimensions during the high-flow period. However, the computational model predicts a much slower rate of wavelength increase, and also performs less accurately during the final low-flow period, where the wavelength remains essentially constant, rather than decreasing. In addition, the numerical results show less variability in bedform wavelength and height than the measured values; the bedform shape is also somewhat different. Based on observations, these discrepancies may result from the simplified model for sediment particle step lengths used in the computational approach. Experiments show that the particle step length varies spatially and temporally over the bedforms during the evolution process. Assuming a constant value for the step length neglects the role of flow alterations in the bedload sediment-transport process, which appears to result in predicted bedform wavelength changes smaller than those observed. However, observations also suggest that three-dimensional effects play at least some role in the decrease of bedform wavelength, so incorporating better models for particle hop lengths alone may not be sufficient to improve model predictions.

Earth Surface Processes and Landforms↗