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

USGS2022

The future of coastal monitoring through satellite remote sensing

Satellite remote sensing is transforming coastal science from a “data-poor” field into a “data-rich” field. Sandy beaches are dynamic landscapes that change in response to long-term pressures, short-term pulses, and anthropogenic interventions. Until recently, the rate and breadth of beach change have outpaced our ability to monitor those changes, due to the spatiotemporal limitations of our observational capacity. Over the past several decades, only a handful of beaches worldwide have been regularly monitored with accurate yet expensive in situ surveys. The long-term coastal-change data of these few well-monitored beaches have led to in-depth understanding of many site-specific coastal processes. However, because the best-monitored beaches are not representative of all beaches, much remains unknown about the processes and fate of the other >99% of unmonitored beaches worldwide. The fleet of Earth-observing satellites has enabled multiscale monitoring of beaches, for the very first time, by providing imagery with global coverage and up to daily frequency. The long-standing and ever-expanding archive of satellite imagery will enable coastal scientists to investigate coastal change at sites vulnerable to future sea-level rise, that is, (almost) everywhere. In the past decade, our capability to observe coastal change from space has grown substantially with computing and algorithmic power. Yet, further advances are needed in automating monitoring using machine learning, deep learning, and computer vision to fully leverage this massive treasure trove of data. Extensive monitoring and investigation of the causes and effects of coastal change at the requisite spatiotemporal scales will provide coastal managers with additional, valuable information to evaluate problems and solutions, addressing the potential for widespread beach loss due to accelerated sea-level rise, development, and reduced sediment supply. Monitoring from Earth-observing satellites is currently the only means of providing seamless data with high spatiotemporal resolution at the global scale of the impending impacts of climate change on coastal systems.

Cambridge Prisms: Coastal Futures↗

USGS2018

The Polar WRF downscaled historical and projected twenty-first century climate for the coast and foothills of Arctic Alaska

Climate change is most pronounced in the northern high latitude region. Yet, climate observations are unable to fully capture regional-scale dynamics due to the sparse weather station coverage, which limits our ability to make reliable climate-based assessments. A set of simulated data products was therefore developed for the North Slope of Alaska through a dynamical downscaling approach. The polar-optimized Weather Research and Forecast (Polar WRF) model was forced by three sources: The ERA-interim reanalysis data (for 1979–2014), the Community Earth System Model 1.0 (CESM1.0) historical simulation (for 1950–2005), and the CESM1.0 projected (for 2006–2100) simulations in two Representative Concentration Pathways (RCP4.5 and RCP8.5) scenarios. Climatic variables were produced in a 10-km grid spacing and a 3-h interval. The ERA-interim forced WRF (ERA-WRF) proves the value of dynamical downscaling, which yields more realistic topographical-induced precipitation and air temperature, as well as corrects underestimations in observed precipitation. In summary, dry and cold biases to the north of the Brooks Range are presented in ERA-WRF, while CESM forced WRF (CESM-WRF) holds wet and warm biases in its historical period. A linear scaling method allowed for an adjustment of the biases, while keeping the majority of the variability and extreme values of modeled precipitation and air temperature. CESM-WRF under RCP 4.5 scenario projects smaller increase in precipitation and air temperature than observed in the historical CESM-WRF product, while the CESM-WRF under RCP 8.5 scenario shows larger changes. The fine spatial and temporal resolution, long temporal coverage, and multi-scenario projections jointly make the dataset appropriate to address a myriad of physical and biological changes occurring on the North Slope of Alaska.

Alaska↗

USGS2026

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

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

Hawaii↗

USGS2004

Three-dimensional distribution of gas hydrate beneath southern Hydrate Ridge: Constraints from ODP Leg 204

Large uncertainties about the energy resource potential and role in global climate change of gas hydrates result from uncertainty about how much hydrate is contained in marine sediments. During Leg 204 of the Ocean Drilling Program (ODP) to the accretionary complex of the Cascadia subduction zone, we sampled the gas hydrate stability zone (GHSZ) from the seafloor to its base in contrasting geological settings defined by a 3D seismic survey. By integrating results from different methods, including several new techniques developed for Leg 204, we overcome the problem of spatial under-sampling inherent in robust methods traditionally used for estimating the hydrate content of cores and obtain a high-resolution, quantitative estimate of the total amount and spatial variability of gas hydrate in this structural system. We conclude that high gas hydrate content (30-40% of pore space or 20-26% of total volume) is restricted to the upper tens of meters below the seafloor near the summit of the structure, where vigorous fluid venting occurs. Elsewhere, the average gas hydrate content of the sediments in the gas hydrate stability zone is generally <2% of the pore space, although this estimate may increase by a factor of 2 when patchy zones of locally higher gas hydrate content are included in the calculation. These patchy zones are structurally and stratigraphically controlled, contain up to 20% hydrate in the pore space when averaged over zones ???10 m thick, and may occur in up to ???20% of the region imaged by 3D seismic data. This heterogeneous gas hydrate distribution is an important constraint on models of gas hydrate formation in marine sediments and the response of the sediments to tectonic and environmental change. ?? 2004 Published by Elsevier B.V.

Earth and Planetary Science Letters↗

USGS2016

Users and uses of Landsat 8 satellite imagery—2014 survey results

Executive Summary In 2013, Landsat 8 began adding high quality, global, moderate-resolution imagery to the more than 40-year archive of Landsat imagery. To assess the potential effects of the availability of Landsat 8 imagery on users and their work, the U.S. Geological Survey (USGS) Land Remote Sensing Program (LRS) initiated a survey of Landsat users. The objectives of the survey were to 1. Characterize various Landsat user groups, such as United States (U.S.) and international users and Landsat 8 and non-Landsat 8 users; 2. Identify any differences among user groups in uses and preferences; 3. Measure the importance of and satisfaction with Landsat 8 attributes; 4. Assess the importance to users of the frequency of usable imagery; and 5. Determine any challenges in using Landsat 8. The online survey was sent to 51,617 Landsat users registered with USGS in May 2014. Almost 13,000 people responded to the survey for a response rate of 25 percent (n = 12,966). Current Landsat users (users who had used Landsat in their work in the year prior to the survey) composed 89 percent of the sample (n = 11,549) and past Landsat users composed 11 percent (n = 1,417). The results reported here apply to current Landsat users registered with the USGS Earth Resources Observation and Science (EROS) Center. Users from 161 countries responded to the survey. Of those, 19 percent were citizens or permanent residents of the United States and 81 percent resided in other countries. More than 70 percent of current users had used Landsat 8 in the year prior to the survey. The majority of Landsat 8 users (65 percent) were established users who used Landsat imagery regularly both before and after Landsat 8 imagery became available. The average current Landsat user was male, 36 years old, and highly educated, with 9 years of experience using satellite imagery or geographic information system (GIS) software. Landsat 8 users had, on average, two more years of experience than non-Landsat 8 users. Users were employed predominantly by academic institutions (65 percent), followed by private businesses (13 percent), Federal governments (10 percent), State and local governments (6 percent), and nonprofit organizations (6 percent). Of the Landsat imagery obtained in the past year by current users, on average 31 percent came from a Landsat 8 sensor. An equivalent amount came from the Landsat 7 ETM+ sensor (33 percent); slightly less came from Landsats 4 and 5 TM sensors (27 percent). Much less came from Landsats 1 through 5 MSS sensors (5 percent). Overall, more than a third of users’ work used Landsat imagery (38 percent). Of this work, on average, 37 percent of the work was operational. Landsat 8 users considered a greater proportion of their work operational than non-Landsat 8 users (39 percent compared with 29 percent). Environmental sciences and management were the most commonly selected primary applications (selected by 42 percent of users). Land use/land cover (23 percent) was the second most commonly selected primary application, followed by education (12 percent), agriculture (9 percent), and planning and development (6 percent). Landsat 8 users were asked to rank the importance of certain attributes in determining whether to use Landsat 8 imagery in their work. The archive was ranked most important, followed by cost, spatial resolution, extent of coverage, data quality, and frequency of revisit. Users were asked how satisfied they were with these same attributes as they currently apply to Landsat 8 imagery. On average, users were most satisfied with lack of cost, extent of coverage, data quality, and the archive, but they were satisfied with all attributes. Users were asked how often they needed Landsat imagery to meet various requirements for their primary application. The survey question specifically asked how often users needed usable imagery, which differs from how often they would like the Landsat satellites to acquire an image. Users were asked to identify their needed frequency of usable imagery for the following levels: 1. Threshold level—the minimum frequency of usable imagery needed to be of any value to their primary application. 2. Breakthrough level—the frequency of usable imagery that would result in a significant improvement for their primary application of the imagery. 3. Target level—the frequency of usable imagery that would only provide a limited additional increase in the expected performance for their primary application. To meet the threshold level, three-quarters of users needed usable imagery every 17 days or less frequently. At the breakthrough level, two-thirds of users (64 percent) needed a usable image every 5–16 days. The current constellation of two satellites (Landsat 7 and 8) is capable of meeting the threshold and breakthrough needs of most users at least some of the time, but a single satellite would be highly unlikely to do so. Two-fifths of users (40 percent) felt that usable imagery provided every 4 days or more frequently would meet their target level which the current Landsat constellation cannot provide. Landsat 8 users were significantly more likely than non-Landsat 8 users to need usable imagery more frequently to meet their target levels. Additionally, U.S. Landsat 8 users were significantly more likely than other Landsat users to need usable imagery more frequently in order meet both their breakthrough and target levels. To explore the effect of the availability of Landsat 8 imagery on Landsat imagery use in general, established users (those who had consistently used Landsat imagery both before and after Landsat 8 imagery became available) using Landsat 8 imagery were asked about changes in the amount of Landsat imagery they used. The majority of established users using Landsat 8 imagery (60 percent) reported an average increase of 51 percent in the number of scenes obtained after Landsat 8 imagery became available. Landsat 8 users were asked if they had encountered challenges in using Landsat 8 whereas non-Landsat 8 users were asked if such challenges had played a role in why they were not using Landsat 8 imagery. Although many users did not encounter challenges when using or trying to use Landsat 8 data, slightly less than 30 percent did encounter issues with processing the data to a usable point. The most common issue reported was not being able to create or have access to a surface reflectance corrected product. Other challenges were related to the file sizes of images being too large to download, store, or analyze. There were no statistically significant differences between Landsat 8 and non-Landsat 8 users in terms of challenges encountered when using or trying to use the imagery, which indicates that users were not unduly discouraged by the challenges they may have encountered. When asked about potential consequences of not using Landsat 8, more than half of the non-Landsat 8 users did not report detrimental effects on their work from not using the imagery. Of those who did report detrimental effects, decreased quality of work, decreased scope of work, and increased time spent on work were the most common.

Open-File Report↗

USGS2011

Sedimentary response to orogenic exhumation in the northern Rocky Mountain Basin and Range province, Flint Creek basin, west-central Montana

Middle Eocene through Upper Miocene sedimentary and volcanic rocks of the Flint Creek basin in western Montana accumulated during a period of significant paleoclimatic change and extension across the northern Rocky Mountain Basin and Range province. Gravity modelling, borehole data, and geologic mapping from the Flint Creek basin indicate that subsidence was focused along an extensionally reactivated Sevier thrust fault, which accommodated up to 800 m of basin fill while relaying stress between the dextral transtensional Lewis and Clark lineament to the north and the Anaconda core complex to the south. Northwesterly paleocurrent indicators, foliated metamorphic lithics, 64 Ma (40Ar/39Ar) muscovite grains, and 76 Ma (U-Pb) zircons in a ca. 27 Ma arkosic sandstone are consistent with Oligocene exhumation and erosion of the Anaconda core complex. The core complex and volcanic and magmatic rocks in its hangingwall created an important drainage divide during the Paleogene shedding detritus to the NNW and ESE. Following a major period of Early Miocene tectonism and erosion, regional drainage networks were reorganized such that paleoflow in the Flint Creek basin flowed east into an internally drained saline lake system. Renewed tectonism during Middle to Late Miocene time reestablished a west-directed drainage that is recorded by fluvial strata within a Late Miocene paleovalley. These tectonic reorganizations and associated drainage divide explain observed discrepancies in provenance studies across the province. Regional correlation of unconformities and lithofacies mapping in the Flint Creek basin suggest that localized tectonism and relative base level fluctuations controlled lithostratigraphic architecture.

Canadian Journal of Earth Sciences↗

USGS1977

Voyager imaging experiment

The overall objective of this experiment is exploratory reconnaissance of Jupiter, Saturn, their satellites, and Saturn's rings. Such reconnaissance, at resolutions and phase angles unobtainable from Earth, can be expected to provide much new data relevant to the atmospheric and/or surface properties of these bodies. The experiment also has the following specific objectives: Observe and characterize the global circulation of the atmospheres of Jupiter and Saturn; Determine the horizontal and vertical structure of the visible clouds and establish their relationship to the belted appearance and dynamical properties of the planetary atmospheres; Determine the vertical structure of high, optically-thin, scattering layers on Jupiter and Saturn; Determine the nature of anomalous features such as the Great Red Spot, South Equatorial Belt disturbances, etc.; Characterize the nature of the colored material in the clouds of Jupiter and Saturn, and identify the nature and sources of chromophores on Io and Titan; Perform comparative geologic studies of many satellites at less than 15-km resolution; Map and characterize the geologic structure of several satellites at high resolution (???1 km); Investigate the existence and nature of atmospheres on the satellites; Determine the mass, size, and shape of many of the satellites by direct measurement; Determine the direction of the spin axes and periods of rotation of several satellites, and establish coordinate systems for the larger satellites; Map the radial distribution of material in Saturn's rings at high resolution; Determine the optical scattering properties of the primaries, rings, and satellites at several wavelengths and phase angles; Search for novel physical phenomena, e.g., phenomena associated with the Io flux tube, meteors, aurorae, lightning, or satellite shadows. ?? 1977 D. Reidel Publishing Company.

Space Science Reviews↗

USGS2011

Self-potential investigations of a gravel bar in a restored river corridor

Self-potentials (SP) are sensitive to water fluxes and concentration gradients in both saturated and unsaturated geological media, but quantitative interpretations of SP field data may often be hindered by the superposition of different source contributions and time-varying electrode potentials. Self-potential mapping and close to two months of SP monitoring on a gravel bar were performed to investigate the origins of SP signals at a restored river section of the Thur River in northeastern Switzerland. The SP mapping and subsequent inversion of the data indicate that the SP sources are mainly located in the upper few meters in regions of soil cover rather than bare gravel. Wavelet analyses of the time-series indicate a strong, but non-linear influence of water table and water content variations, as well as rainfall intensity on the recorded SP signals. Modeling of the SP response with respect to an increase in the water table elevation and precipitation indicate that the distribution of soil properties in the vadose zone has a very strong influence. We conclude that the observed SP responses on the gravel bar are more complicated than previously proposed semi-empiric relationships between SP signals and hydraulic head or the thickness of the vadose zone. We suggest that future SP monitoring in restored river corridors should either focus on quantifying vadose zone processes by installing vertical profiles of closely spaced SP electrodes or by installing the electrodes within the river to avoid signals arising from vadose zone processes and time-varying electrochemical conditions in the vicinity of the electrodes.

Lake Constance, Neunforn, Thur catchment, Thur val↗

USGS2013

Prediction, time variance, and classification of hydraulic response to recharge in two karst aquifers

Many karst aquifers are rapidly filled and depleted and therefore are likely to be susceptible to changes in short-term climate variability. Here we explore methods that could be applied to model site-specific hydraulic responses, with the intent of simulating these responses to different climate scenarios from high-resolution climate models. We compare hydraulic responses (spring flow, groundwater level, stream base flow, and cave drip) at several sites in two karst aquifers: the Edwards aquifer (Texas, USA) and the Madison aquifer (South Dakota, USA). A lumped-parameter model simulates nonlinear soil moisture changes for estimation of recharge, and a time-variant convolution model simulates the aquifer response to this recharge. Model fit to data is 2.4% better for calibration periods than for validation periods according to the Nash–Sutcliffe coefficient of efficiency, which ranges from 0.53 to 0.94 for validation periods. We use metrics that describe the shapes of the impulse-response functions (IRFs) obtained from convolution modeling to make comparisons in the distribution of response times among sites and between aquifers. Time-variant IRFs were applied to 62% of the sites. Principal component analysis (PCA) of metrics describing the shapes of the IRFs indicates three principal components that together account for 84% of the variability in IRF shape: the first is related to IRF skewness and temporal spread and accounts for 51% of the variability; the second and third largely are related to time-variant properties and together account for 33% of the variability. Sites with IRFs that dominantly comprise exponential curves are separated geographically from those dominantly comprising lognormal curves in both aquifers as a result of spatial heterogeneity. The use of multiple IRF metrics in PCA is a novel method to characterize, compare, and classify the way in which different sites and aquifers respond to recharge. As convolution models are developed for additional aquifers, they could contribute to an IRF database and a general classification system for karst aquifers.

South Dakota, Texas↗

USGS1995

A mass proportion method for calculating melting reactions and application to melting of model upper mantle lherzolite

We present a method for calculating quantitative melting reactions in systems with multiple solid solutions that accounts for changes in the mass proportions of phases between two points at different temperatures along a melting curve. This method can be applied to any data set that defines the phase proportions along a melting curve. The method yields the net change in mass proportion of all phases for the chosen melting interval, and gives an average reaction for the melting path. Instantaneous melting reactions can be approximated closely by choosing sufficiently small melting intervals. As an application of the method, reactions for melting of model upper mantle peridotite are calculated using data from the system CaO-MgO-Al 2 O 3 -SiO 2 -Na 2 O (CMASN) over the pressure interval 0.7 – 3.5 GPa. Throughout almost this entire pressure range, melting of model lherzolite involves the crystallization of one or more solid phases, and is analogous to melting at a peritectic invariant point. In addition, we show that melting reactions for small melting intervals (< 5%) along the solidus of mantle peridotite are significantly different from those calculated for large melting intervals. For large melting intervals (> 10%), reaction stoichiometries calculated in CMASN are usually in good agreement with those available for melting of natural peridotite. The coefficients of melting reactions calculated from this method can be used in equations that describe the behavior of trace elements during melting. We compare results from near-fractional melting models using (1) melting reactions and rock modes from CMASN, and (2) constant reactions representative of those used in the literature. In modeling trace element abundances in melt, significant differences arise for some elements at low degrees of melting (< 10%). In modeling element abundances in the residue, differences increase with increase in degree of melting. Reactions calculated along the model lherzolite solidus in CMASN are the only ones available at present for small degrees of melting so we recommend them for accurate trace element modeling of natural lherzolite.

Earth and Planetary Science Letters↗

USGS2001

Rhenium-osmium systematics of calcium-aluminium-rich inclusions in carbonaceous chondrites

The Re-Os isotopic systematics of calcium-aluminium-rich inclusions (CAIs) in chondrites were investigated in order to shed light on the behavior of the Re-Os system in bulk chondrites, and to constrain the timing of chemical fractionation in primitive chondrites. CAIs with relatively unfractionated rare earth element (REE) patterns (groups I, III, V, VI) define a narrow range of 187Re/188Os (0.3764-0.4443) and 187Os/188Os (0.12599-0.12717), and high but variable Re and Os abundances (3209-41,820 ppb Os). In contrast, CAIs that show depletions in highly refractory elements and strongly fractionated REE patterns (group II) also show a much larger range in 187Re/188Os (0.409-0.535) and 187Os/188Os (0.12695-0.13770), and greater than an order of magnitude lower Re and Os abundances than other groups (e.g., 75.7-680.2 ppb Os). Sixteen bulk CAIs and CAI splits plot within analytical uncertainty of a 4558 Ga reference isochron, as is expected for materials of this antiquity. Eight samples, however, plot off the isochron. Several possible reasons for these deviations are discussed. Data for multiple splits of one CAI indicate that the nonisochronous behavior for at least this CAI is the result of Re-Os reequilibration at approximately 1.6 Ga. Thus, the most likely explanation for the deviations of most of the nonisochronous CAIs is late-stage open-system behavior of Re and Os in the asteroidal environment. The 187Os/188Os-Os systematics of CAIs are consistent with previous models that indicate group II CAIs are mixtures of components that lost the bulk of their highly refractory elements in a previous condensation event and a minor second component that provided refractory elements at chondritic relative proportions. The high Re/Os of group II CAIs relative to other CAIs and chondrite bulk rocks may have been caused by variable mobilization of Re and Os during medium- to low-temperature parent body alteration ??4.5 Ga ago. This model is favored over nebular models, which pose several difficulties. The narrow range of 187Os/188Os in group I, III, V, and VI bulk CAIs, and the agreement with 187Os/188Os of whole rock carbonaceous chondrites suggest that on a bulk inclusion scale, secondary alteration only modestly fractionated Re/Os in these CAIs. The average of 187Os/188Os for group I, III, V, and VI CAIs is indistinguishable from average CI chondrites, indicating a modern solar system value for 187Os/188Os of 0.12650, corresponding to a 187Re/188Os of 0.3964. Copyright ?? 2001 Elsevier Science Ltd.

Geochimica et Cosmochimica Acta↗

USGS2013

Nature's refineries — Metals and metalloids in arc volcanoes

Chemical data for fumaroles and for atmospheric gas and ash plumes from active arc volcanoes provide glimpses of the rates of release of metal and metalloids, such as Tl and Cd, from shallow and mid-crust magmas. Data from copper deposits formed in ancient volcanoes at depths of up to about 1500 m in the fractures below paleo-fumaroles, and at around 2000–4000 m in association with sub-volcanic intrusions (porphyry copper deposits) provide evidence of sub-surface deposition of Cu–Au–Ag–Mo and a range of other minor elements including Te, Se, As and Sb. These deposits, or ‘sinks’, of metals consistently record sustained histories of magmatic gas streaming through volcanic systems interspersed by continuing intrusive and eruptive activity. Here we integrate data from ancient and modern volcanic systems and show that the fluxes of metals and metalloids are controlled by a) the maintenance of fracture permeability in the stressed crust below volcanoes and b) the chemical processes that are triggered as magmatic gas, initially undersaturated with metals and metalloids, expands from lithostatic to very low pressure conditions through fracture arrays. The recognition of gas streaming may also account for the phenomenon of ‘excess degassing’, and defines an integral, but generally understated, component of active volcanic systems – a volcanic gas core – that is likely to be integral to the progression of eruptions to Plinean state. Destabilization of solvated molecular metal and metalloid species in magmatic gas mixtures and changes in their redox state are triggered, as it expands to the surface by abrupt pressure drops, or throttles' in the fracture array that guides expansion to the surface. The electronically harder, low electronegativity metals, such as copper and iron, deposit rapidly in response to expansion followed more slowly by arsenic with antimony as sulfosalts. Heavy, large radius, softer elements such as bismuth, lead, and thallium along with cadmium are strongly fractionated along the way, eventually venting their excess along with SO 2 , CO 2 , and other components of the carrier gas, into the atmosphere. These elements, many of which are toxic, may also be dispersed by mixing with groundwater in the permeable crust below volcanoes and generate potential health risks due to Hg, As, and Se contamination of drinking water resources.

Earth-Science Reviews↗

USGS2021

Detrital zircon record of magmatism and sediment dispersal across the North American Cordilleran arc system (28-48°N)

As zircon U-Pb geochronology has become a leading method in sediment provenance studies and basin analysis over the past 20 years, the volume of detrital zircon data made available in published literature has enabled researchers to go beyond source-to-sink provenance studies to explore increasingly complex geologic problems. In this review, we utilize the growing body of detrital zircon data acquired from Jurassic-Paleocene forearc and foreland basin strata of the North American Cordillera to investigate the Mesozoic to earliest Cenozoic evolution of the arc and its associated basins between 28°N and 48°N. Our compilation includes 830 detrital zircon samples (101,898 individual ages) from 70 studies published between 2000 and 2020. For comparative purposes, we also compile 1307 igneous zircon U-Pb ages that characterize the magmatic history of the arc. We place primary emphasis on detrital zircon ages between 251 and 56 Ma that we infer to be uniquely derived from magmatic sources in the arc. Informed by existing knowledge of magmatic, structural, and sedimentological processes that acted on the orogen, we investigate spatial and temporal trends in these “arc-derived zircon” to establish a detrital record of arc magmatism, investigate source-to-sink relationships between the arc and adjacent basins, and discuss controls on sediment dispersal across the orogen. Our review shows that compilations of detrital zircon data from the Cordilleran forearc and foreland basin systems are excellent proxies for arc magmatism because the basins are enriched in arc-derived zircon and compilations provide space- and time-integrated records of crystallization ages. The compiled detrital zircon data support a history of continuous arc magmatism throughout Mesozoic and earliest Cenozoic time, characterized by low-volume magmatism from Triassic-Early Jurassic time (~251–174 Ma) and episodic higher-volume magmatism from Middle Jurassic-Late Cretaceous time (~174–66 Ma). These trends elucidate the initiation and timing of magmatic events at the orogen-scale and corroborate our understanding of cyclic arc behavior. Detrital zircon distributions are spatially and temporally variable both within and across basins, which we discuss relative to topographic development of the orogen and attendant responses of sediment dispersal systems. Detrital zircon distributions in the forearc signal rapid transfer of sediment from the arc to basins dominantly via fluvial processes. In contrast, detrital zircon distributions across the foreland reflect the presence of topographic barriers in the hinterland region of the arc that effectively isolated parts of the foreland. The presence of hinterland topography in turn highlights the important role of ash-fall events in delivering arc-derived zircon to the foreland, underscoring the need to consider ash-fall processes in paleodrainage reconstructions. These broad regional trends, and in general the close linkage between orogenic process and sediment dispersal, emerge from our compilation because it averages out much of the local variability observed in studies of more limited geographic or temporal extent.

North American Cordilleran arc system↗

USGS2011

EAARL Coastal Topography--Cape Canaveral, Florida, 2009: First Surface

These remotely sensed, geographically referenced elevation measurements of lidar-derived first-surface (FS) topography datasets were produced collaboratively by the U.S. Geological Survey (USGS), St. Petersburg Coastal and Marine Science Center, St. Petersburg, FL, and the National Aeronautics and Space Administration (NASA), Kennedy Space Center, FL. This project provides highly detailed and accurate datasets of a portion of the eastern Florida coastline beachface, acquired on May 28, 2009. 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 Advanced Airborne 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 red-green-blue (RGB) digital camera, a high-resolution multispectral 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 aircraft, but the instrument was deployed on a Pilatus PC-6. 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, and the resulting data were then 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 presurvey 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 used routinely 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.

Data Series↗

USGS1985

Evaluating earthquake hazards in the Los Angeles region— An earth-science perspective

Potentially destructive earthquakes are inevitable in the Los Angeles region of California, but hazards prediction can provide a basis for reducing damage and loss. This volume identifies the principal geologically controlled earthquake hazards of the region (surface faulting, strong shaking, ground failure, and tsunamis), summarizes methods for characterizing their extent and severity, and suggests opportunities for their reduction. Two systems of active faults generate earthquakes in the Los Angeles region: northwest-trending, chiefly horizontal-slip faults, such as the San Andreas, and west-trending, chiefly vertical-slip faults, such as those of the Transverse Ranges. Faults in these two systems have produced more than 40 damaging earthquakes since 1800. Ninety-five faults have slipped in late Quaternary time (approximately the past 750,000 yr) and are judged capable of generating future moderate to large earthquakes and displacing the ground surface. Average rates of late Quaternary slip or separation along these faults provide an index of their relative activity. The San Andreas and San Jacinto faults have slip rates measured in tens of millimeters per year, but most other faults have rates of about 1 mm/yr or less. Intermediate rates of as much as 6 mm/yr characterize a belt of Transverse Ranges faults that extends from near Santa Barbara to near San Bernardino. The dimensions of late Quaternary faults provide a basis for estimating the maximum sizes of likely future earthquakes in the Los Angeles region: moment magnitude .(M) 8 for the San Andreas, M 7 for the other northwest-trending elements of that fault system, and M 7.5 for the Transverse Ranges faults. Geologic and seismologic evidence along these faults, however, suggests that, for planning and designing noncritical facilities, appropriate sizes would be M 8 for the San Andreas, M 7 for the San Jacinto, M 6.5 for other northwest-trending faults, and M 6.5 to 7 for the Transverse Ranges faults. The geologic and seismologic record indicates that parts of the San Andreas and San Jacinto faults have generated major earthquakes having recurrence intervals of several tens to a few hundred years. In contrast, the geologic evidence at points along other active faults suggests recurrence intervals measured in many hundreds to several thousands of years. The distribution and character of late Quaternary surface faulting permit estimation of the likely location, style, and amount of future surface displacements. An extensive body of geologic and geotechnical information is used to evaluate areal differences in future levels of shaking. Bedrock and alluvial deposits are differentiated according to the physical properties that control shaking response; maps of these properties are prepared by analyzing existing geologic and soils maps, the geomorphology of surficial units, and. geotechnical data obtained from boreholes. The shear-wave velocities of near-surface geologic units must be estimated for some methods of evaluating shaking potential. Regional-scale maps of highly generalized shearwave velocity groups, based on the age and texture of exposed geologic units and on a simple two-dimensional model of Quaternary sediment distribution, provide a first approximation of the areal variability in shaking response. More accurate depictions of near-surface shear-wave velocity useful for predicting ground-motion parameters take into account the thickness of the Quaternary deposits, vertical variations in sediment .type, and the correlation of shear-wave velocity with standard penetration resistance of different sediments. A map of the upper Santa Ana River basin showing shear-wave velocities to depths equal to one-quarter wavelength of a 1-s shear wave demonstrates the three-dimensional mapping procedure. Four methods for predicting the distribution and strength of shaking from future earthquakes are presented. These techniques use different measures of strong-motion

California↗

USGS2009

EAARL Coastal Topography-Pearl River Delta 2008: First Surface

These remotely sensed, geographically referenced elevation measurements of Lidar-derived first surface (FS) topography were produced as a collaborative effort between the U.S. Geological Survey (USGS), Florida Integrated Science Center (FISC), St. Petersburg, FL; the University of New Orleans (UNO), Pontchartrain Institute for Environmental Sciences (PIES), New Orleans, LA; and the National Aeronautics and Space Administration (NASA), Wallops Flight Facility, VA. This project provides highly detailed and accurate datasets of a portion of the Pearl River Delta in Louisiana and Mississippi, acquired March 9-11, 2008. 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 Advanced Airborne 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 red-green-blue (RGB) 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 submeter 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, and the resulting data were then 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 presurvey 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 used routinely to create maps that represent submerged or first surface topography. Specialized filtering algorithms have been implemented to determine the 'bare earth' under vegetation from a point cloud of last return elevations.

Data Series↗

USGS2018

Characterization and evaluation of controls on post-fire streamflow response across western U.S. watersheds

This research investigates the impact of wildfires on watershed flow regimes, specifically focusing on evaluation of fire events within specified hydroclimatic regions in the western United States, and evaluating the impact of climate and geophysical variables on response. Eighty-two watersheds were identified with at least 10 years of continuous pre-fire daily streamflow records and 5 years of continuous post-fire daily flow records. Percent change in annual runoff ratio, low flows, high flows, peak flows, number of zero flow days, baseflow index, and Richards–Baker flashiness index were calculated for each watershed using pre- and post-fire periods. Independent variables were identified for each watershed and fire event, including topographic, vegetation, climate, burn severity, percent area burned, and soils data. Results show that low flows, high flows, and peak flows increase in the first 2 years following a wildfire and decrease over time. Relative response was used to scale response variables with the respective percent area of watershed burned in order to compare regional differences in watershed response. To account for variability in precipitation events, runoff ratio was used to compare runoff directly to PRISM precipitation estimates. To account for regional differences in climate patterns, watersheds were divided into nine regions, or clusters, through k-means clustering using climate data, and regression models were produced for watersheds grouped by total area burned. Watersheds in Cluster 9 (eastern California, western Nevada, Oregon) demonstrate a small negative response to observed flow regimes after fire. Cluster 8 watersheds (coastal California) display the greatest flow responses, typically within the first year following wildfire. Most other watersheds show a positive mean relative response. In addition, simple regression models show low correlation between percent watershed burned and streamflow response, implying that other watershed factors strongly influence response. Spearman correlation identified NDVI, aridity index, percent of a watershed's precipitation that falls as rain, and slope as being positively correlated with post-fire streamflow response. This metric also suggested a negative correlation between response and the soil erodibility factor, watershed area, and percent low burn severity. Regression models identified only moderate burn severity and watershed area as being consistently positively/negatively correlated, respectively, with response. The random forest model identified only slope and percent area burned as significant watershed parameters controlling response. Results will help inform post-fire runoff management decisions by helping to identify expected changes to flow regimes, as well as facilitate parameterization for model application in burned watersheds.

Hydrology and Earth System Sciences↗

USGS2016

Joint analysis of geodetic and earthquake fault-plane solution data to constrain magmatic sources: A case study from Kīlauea Volcano

A joint analysis of geodetic and seismic datasets from Kīlauea Volcano during a period of magmatic unrest in 2006 demonstrates the effectiveness of this combination for testing and constraining models of magma dynamics for a complex, multi-source system. At the end of 2003, Kīlauea's summit began a four-year-long period of inflation due to a surge in magma supply to the volcano. In 2006, for the first time since 1982, Kīlauea's Southwest Rift Zone (SWRZ) also experienced inflation. To investigate the characteristics of active magma sources and the nature of their interactions with faults in the SWRZ during 2006, we integrate, through Coulomb stress modeling, contemporary geodetic data from InSAR and GPS with a new catalogue of double-couple fault-plane solutions for volcano-tectonic earthquakes. We define two periods of inflation during 2006 based on the rate of deformation measured in daily GPS data, spanning February to 15 March 2006 (Period 1) and 16 March to 30 September 2006 (Period 2). InSAR data for these two periods are inverted to determine the position, change in size, and shape of inflation sources in each period. Our new models are consistent with microseismic activity from each period. They suggest that, during Period 1, deformation in the SWRZ can be explained by pressurization of magma in a spherical reservoir beneath the south caldera, and that, during Period 2, magma was also aseismically intruded farther to the southwest into the SWRZ along a sub-horizontal plane. Our Coulomb stress analysis shows that the microseismicity recorded in the SWRZ is induced by overpressurization of the south caldera reservoir, and not by magma intrusion into the SWRZ. This study highlights the importance of a joint analysis of independent geophysical datasets to fully constrain the nature of magma accumulation.

Hawaii↗