SEARCH · Search USGS

Results for “Analytical Letters”

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 37 records · Page 2Linked to original sources

USGS2004

Assessing the resolution-dependent utility of tomograms for geostatistics

Geophysical tomograms are used increasingly as auxiliary data for geostatistical modeling of aquifer and reservoir properties. The correlation between tomographic estimates and hydrogeologic properties is commonly based on laboratory measurements, co-located measurements at boreholes, or petrophysical models. The inferred correlation is assumed uniform throughout the interwell region; however, tomographic resolution varies spatially due to acquisition geometry, regularization, data error, and the physics underlying the geophysical measurements. Blurring and inversion artifacts are expected in regions traversed by few or only low-angle raypaths. In the context of radar traveltime tomography, we derive analytical models for (1) the variance of tomographic estimates, (2) the spatially variable correlation with a hydrologic parameter of interest, and (3) the spatial covariance of tomographic estimates. Synthetic examples demonstrate that tomograms of qualitative value may have limited utility for geostatistics; moreover, the imprint of regularization may preclude inference of meaningful spatial statistics from tomograms.

Geophysical Research Letters↗

USGS1986

Maps and tables showing data and analyses of semiquantitative emmission spectrometry and atomic-absorption spectrophotometry of rock samples, Ugashik, Bristol Bay, and part of Karluk quadrangles, Alaska

The accompanying maps and tables show analytical data and data analyses from rock samples collected in conjunction with geologic mapping in the Ugashik, Bristol Bay and western Karluck quadrangles from 1979 through 1981. This work was conducted under the auspices of the Alaska Mineral Resource Assessment Program (AMRAP). A total of 337 samples were collected for analysis, primarily in areas of surficial alteration. The sample locations are shown on sheet 1: they are concentrated along the Pacific Ocean side of the area because the Bristol Bay lowlands part of the map is predominantly unconsolidated Quaternary deposits. Sample collection was by the following people, with their respective two letter identifying code shown in parentheses: W.H. Allaway (AY), J.E. Case (CE), D.P. Cox (CX), R.L. Detterman, (DT), T.G. Theodore (MK), F.H. Wilson (WS), and M.E. Yount (YB).

Alaska↗

USGS2007

Moment inference from tomograms

Time-lapse geophysical tomography can provide valuable qualitative insights into hydrologic transport phenomena associated with aquifer dynamics, tracer experiments, and engineered remediation. Increasingly, tomograms are used to infer the spatial and/or temporal moments of solute plumes; these moments provide quantitative information about transport processes (e.g., advection, dispersion, and rate-limited mass transfer) and controlling parameters (e.g., permeability, dispersivity, and rate coefficients). The reliability of moments calculated from tomograms is, however, poorly understood because classic approaches to image appraisal (e.g., the model resolution matrix) are not directly applicable to moment inference. Here, we present a semi-analytical approach to construct a moment resolution matrix based on (1) the classic model resolution matrix and (2) image reconstruction from orthogonal moments. Numerical results for radar and electrical-resistivity imaging of solute plumes demonstrate that moment values calculated from tomograms depend strongly on plume location within the tomogram, survey geometry, regularization criteria, and measurement error.

Geophysical Research Letters↗

USGS2005

National parks and protected areas: Appoaches for balancing social, economic, and ecological values

The balance of nature in any strict sense has been upset long ago…The only option we have is to create a new balance objectively determined for each area in accordance with the intended use of that area.” --Aldo Leopold, 1927, in a letter to the Superintendent of Glacier National Park The planning and management staff of state/national parks and protected areas face a complex set of management problems and an uncertain future. Today, internal and external forces are combining with changes in management philosophy (from a species to an ecosystem focus) to create a new direction in natural resources management. National Parks and Protected Areas: Approaches for Balancing Social, Economic and Ecological Values is peerless in its unified treatment of the issues surrounding this subject. From decision-making for planning and management to the principles of ecology and economics, this text examines the analytical methods, information technologies, and planning and management problems associated with protected area planning and management. Protected area managers and students in undergraduate and graduate courses in natural resource management will appreciate this highly readable book.

Book↗

USGS2003

How should environmental stress affect the population dynamics of disease?

We modelled how stress affects the population dynamics of infectious disease. We were specifically concerned with stress that increased susceptibility of uninfected hosts when exposed to infection. If such stresses also reduced resources, fecundity and/or survivorship, there was a reduction in the host carrying capacity. This lowered the contact between infected and uninfected hosts, thereby decreasing transmission. In addition, stress that increased parasite mortality decreased disease. The opposing effects of stress on disease dynamics made it difficult to predict the response of disease to environmental stress. We found analytical solutions with negative, positive, convex and concave associations between disease and stress. Numerical simulations with randomly generated parameter values suggested that the impact of host-specific diseases generally declined with stress while the impact of non-specific (or open) diseases increased with stress. These results help clarify predictions about the interaction between environmental stress and disease in natural populations.

Ecology Letters↗

USGS2025

Estimates of global surface water dynamics harnessing near real-time land cover observations and open science geospatial capabilities

Spatio-temporal changes to our world’s surface water resources are escalating. Translating how these changes impact communities and ecosystems requires time-varying data of Global Surface Water Extents (GSWE). Traditionally, GSWE mapping has been limited to static estimates, with recent efforts focusing on annual averages, frequency and occurrence of long-term variations. Building upon these foundational capabilities, we harnessed remotely sensed Sentinel-2 based near real-time Dynamic World (DW) land cover products to produce the first-of-its-kind 10 m resolution GSWE dataset representing 2015–2023. Our dataset estimated 2.5 million km 2 of permanent waters and 8 million km 2 of seasonal waters worldwide. Comparing our Sentinel-2 based data to contemporary Landsat-based GSWE, we observed that our data mapped less water within the >50% probability of occurrence range, suggesting a lower presence of open permanent water especially in high latitudes, deviating from what we previously learnt from Landsat data. Statistical analysis compared to well-established observational products and widely used GSWE datasets across some of the world’s most ecologically significant regions, including Pantanal in South America and Haor in South Asia, supports the overall physical realism of our data in predicting global open surface water dynamics. Our key contribution is a prototype Open Science operational framework that extracts routinely available DW products, runs geospatial analytics, and creates actionable water information for educators, researchers, and stakeholders at any scale of practical interest. We present examples of this operational capability through instant mapping of flood in Spain and drought in Lake Urmia, Central Asia, frequent monitoring of river extent changes at the Ganges–Brahmaputra confluence, and above all, interoperability with other existing GSWE applications.

Environmental Research Letters↗

USGS2023

New insights into the relationship between mass eruption rate and volcanic column height based on the IVESPA dataset

Rapid and simple estimation of the mass eruption rate (MER) from column height is essential for real-time volcanic hazard management and reconstruction of past explosive eruptions. Using 134 eruptive events from the new Independent Volcanic Eruption Source Parameter Archive (IVESPA, v1.0), we explore empirical MER-height relationships for four measures of column height: spreading level, sulfur dioxide height, and top height from direct observations and as reconstructed from deposits. These relationships show significant differences and highlight limitations of empirical models currently used in operational and research applications. The roles of atmospheric stratification, wind, and humidity remain challenging to detect across the wide range of eruptive conditions spanned in IVESPA, ultimately resulting in empirical relationships outperforming analytical models that account for atmospheric conditions. This finding highlights challenges in constraining the MER-height relation using heterogeneous observations and empirical models, which reinforces the need for improved eruption source parameter data sets and physics-based models.

Geophysical Research Letters↗

USGS2014

Metrics for comparing dynamic earthquake rupture simulations

Earthquakes are complex events that involve a myriad of interactions among multiple geologic features and processes. One of the tools that is available to assist with their study is computer simulation, particularly dynamic rupture simulation. A dynamic rupture simulation is a numerical model of the physical processes that occur during an earthquake. Starting with the fault geometry, friction constitutive law, initial stress conditions, and assumptions about the condition and response of the near‐fault rocks, a dynamic earthquake rupture simulation calculates the evolution of fault slip and stress over time as part of the elastodynamic numerical solution (Ⓔ see the simulation description in the electronic supplement to this article). The complexity of the computations in a dynamic rupture simulation make it challenging to verify that the computer code is operating as intended, because there are no exact analytic solutions against which these codes’ results can be directly compared. One approach for checking if dynamic rupture computer codes are working satisfactorily is to compare each code’s results with the results of other dynamic rupture codes running the same earthquake simulation benchmark. To perform such a comparison consistently, it is necessary to have quantitative metrics. In this paper, we present a new method for quantitatively comparing the results of dynamic earthquake rupture computer simulation codes.

Seismological Research Letters↗

USGS2011

Contribution of global groundwater depletion since 1900 to sea-level rise

Removal of water from terrestrial subsurface storage is a natural consequence of groundwater withdrawals, but global depletion is not well characterized. Cumulative groundwater depletion represents a transfer of mass from land to the oceans that contributes to sea-level rise. Depletion is directly calculated using calibrated groundwater models, analytical approaches, or volumetric budget analyses for multiple aquifer systems. Estimated global groundwater depletion during 1900–2008 totals ~4,500 km 3 , equivalent to a sea-level rise of 12.6 mm (>6% of the total). Furthermore, the rate of groundwater depletion has increased markedly since about 1950, with maximum rates occurring during the most recent period (2000–2008), when it averaged ~145 km 3 /yr (equivalent to 0.40 mm/yr of sea-level rise, or 13% of the reported rate of 3.1 mm/yr during this recent period).

Geophysical Research Letters↗

USGS2009

Landsliding in partially saturated materials

Rainfall‐induced landslides are pervasive in hillslope environments around the world and among the most costly and deadly natural hazards. However, capturing their occurrence with scientific instrumentation in a natural setting is extremely rare. The prevailing thinking on landslide initiation, particularly for those landslides that occur under intense precipitation, is that the failure surface is saturated and has positive pore‐water pressures acting on it. Most analytic methods used for landslide hazard assessment are based on the above perception and assume that the failure surface is located beneath a water table. By monitoring the pore water and soil suction response to rainfall, we observed shallow landslide occurrence under partially saturated conditions for the first time in a natural setting. We show that the partially saturated shallow landslide at this site is predictable using measured soil suction and water content and a novel unified effective stress concept for partially saturated earth materials.

Geophysical Research Letters↗

USGS2020

Goals and development of the Alaska Volcano Observatory Seismic Network and application to forecasting and detecting volcanic eruptions

The Alaska Volcano Observatory (AVO) seismic network has been in operation since 1988 and during this time has grown from 29 to 217 seismic stations providing real-time monitoring of 32 active volcanoes in Alaska, as well as useful data for regional earthquake monitoring. Since 1988, AVO has detected 59 volcanic eruptions at Aleutian arc volcanoes, and 31 of these have been captured by local seismic instrumentation. As part of this monitoring effort, AVO has cataloged more than 120,000 earthquake hypocenters and magnitudes associated with volcanic processes throughout the arc. This high rate of volcanic activity provides an excellent opportunity to study seismicity associated with magmatic and eruptive processes and develop and refine analytical techniques to track volcanic seismicity and warn of hazardous eruptions. The network is currently undergoing an extensive upgrade, replacing aging short-period analog seismometers with digital broadband instruments. These are expected to improve AVO’s seismic capability and further facilitate other geophysical instrumentation such as continuous Global Positioning System receivers, infrasound sensors, and web cams.

Seismological Research Letters↗

USGS2012

Predicted variation of stress orientation with depth near an active fault: Application to the Cajon Pass Scientific Drillhole, southern California

Preliminary measurements of the stress orientation at a depth of 2 km in the Cajon Pass Scientific Drillhole, about 3.6 km from the San Andreas fault in southern California, indicate that the local direction of maximum compression is nearly normal to the fault zone. The measurements have been interpreted to indicate that the regional orientation of the maximum compression is normal to the fault, and taken as evidence for a very weak fault. The orientation of the regional maximum compression expected from plate tectonic arguments is about 66° NE from the strike of the fault. Geodetic data indicate that the orientation of maximum compressive strain rate is about 43° NE from the strike of the fault, and show nearly pure right-lateral shear acting parallel to the fault. These apparent conflicts in the inferred orientation of the axis of maximum compression may be explained in part by a model in which the fault zone is locked over a depth interval in the range of 2-5 to 15 km, but is very weak above and below that interval. An analytic solution for this model predicts orientations of the principal stresses that are nearly parallel or normal to the fault, for observations made close to the fault and above and below the locked zone, but an orientation consistent with plate tectonic estimates at depths corresponding to that of the locked interval. This solution does require, however, a few mm/yr of creep at the surface on the San Andreas or nearby sub-parallel faults (such as the San Jacinto), which has not yet been observed, or a shallow zone near the faults of distributed deformation.

Geophysical Research Letters↗

USGS2017

Inflation leading to a Slow Slip Event and volcanic unrest at Mt. Etna in 2016: Insights from CGPS data

Global Positioning System (CGPS) data from Mount Etna between May 2015 and September 2016 show intense inflation and a concurrent Slow Slip Event (SSE) from 11 December 2015 to 17 May 2016. In May 2016, an eruptive phase started from the summit craters, temporarily stopping the ongoing inflation. The CGPS data presented here give us the opportunity to determine (1) the source of the inflating body, (2) the strain rate parameters highlighting shear strain rate accumulating along NE Rift and S Rift, (3) the magnitude of the SSE, and (4) possible interaction between modeled sources and other flank structures through stress calculations. By analytical inversion, we find an inflating source 5.5 km under the summit (4.4 km below sea level) and flank slip in a fragmented shallow structure accommodating displacements equivalent to a magnitude M w 6.1 earthquake. These large displacements reflect a complex mechanism of rotations indicated by the inversion of CGPS data for strain rate parameters. At the scale of the volcano, these processes can be considered precursors of seismic activity in the eastern flank of the volcano but concentrated mainly on the northern boundary of the mobile eastern flank along the Pernicana Fault and in the area of the Timpe Fault System.

Mt. Etna↗

USGS2017

Imaging shear strength along subduction faults

Subduction faults accumulate stress during long periods of time and release this stress suddenly, during earthquakes, when it reaches a threshold. This threshold, the shear strength, controls the occurrence and magnitude of earthquakes. We consider a 3-D model to derive an analytical expression for how the shear strength depends on the fault geometry, the convergence obliquity, frictional properties, and the stress field orientation. We then use estimates of these different parameters in Japan to infer the distribution of shear strength along a subduction fault. We show that the 2011 M w 9.0 Tohoku earthquake ruptured a fault portion characterized by unusually small variations in static shear strength. This observation is consistent with the hypothesis that large earthquakes preferentially rupture regions with relatively homogeneous shear strength. With increasing constraints on the different parameters at play, our approach could, in the future, help identify favorable locations for large earthquakes.

Geophysical Research Letters↗

USGS2011

Rotational motions for teleseismic surface waves

[1] We report the findings for the first teleseismic six degree‐of‐freedom (6‐DOF) measurements including three components of rotational motions recorded by a sensitive rotation‐rate sensor (model R‐1, made by eentec) and three components of translational motions recorded by a traditional seismometer (STS‐2) at the NACB station in Taiwan. The consistent observations in waveforms of rotational motions and translational motions in sections of Rayleigh and Love waves are presented in reference to the analytical solution for these waves in a half space of Poisson solid. We show that additional information (e.g., Rayleigh wave phase velocity, shear wave velocity of the surface layer) might be exploited from six degree‐of‐freedom recordings of teleseismic events at only one station. We also find significant errors in the translational records of these teleseismic surface waves due to the sensitivity of inertial translation sensors (seismometers) to rotational motions. The result suggests that the effects of such errors need to be counted in surface wave inversions commonly used to derive earthquake source parameters and Earth structure.

Geophysical Research Letters↗

USGS2024

Earth’s free surface complicates inference of absolute stress from earthquake-Induced stress rotations

The stress redistribution from an earthquake can produce localized measurable rotations of the principal stress axes if the absolute level of differential stress in the crust in on the order of the earthquake stress drop. Two simple analytic solutions have been developed to estimate the differential stress from an observed stress rotation. However, each has assumptions that may not be accurate near Earth’s free surface. I model synthetic earthquakes in an elastic half-space, and show that the assumptions of the methods are accurate for strike-slip earthquakes, and for deep dip-slip earthquakes. However, they are incorrect for shallow dip-slip earthquakes. I introduce a free surface correction for one of the methods for dip-slip earthquakes. I revise an analysis of stress rotations due to great subduction zone earthquakes, including this correction. The results support the original conclusion of near complete stress drop for many shallow subduction zone earthquakes.

Geophysical Research Letters↗

USGS2017

Identifying species conservation strategies to reduce disease-associated declines

Emerging infectious diseases (EIDs) are a salient threat to many animal taxa, causing local and global extinctions, altering communities and ecosystem function. The EID chytridiomycosis is a prominent driver of amphibian declines, which is caused by the fungal pathogen Batrachochytrium dendrobatidis (Bd). To guide conservation policy, we developed a predictive decision-analytic model that combines empirical knowledge of host-pathogen metapopulation dynamics with expert judgment regarding effects of management actions, to select from potential conservation strategies. We apply our approach to a boreal toad ( Anaxyrus boreas boreas ) and Bd system, identifying optimal strategies that balance tradeoffs in maximizing toad population persistence and landscape-level distribution, while considering costs. The most robust strategy is expected to reduce the decline of toad breeding sites from 53% to 21% over 50 years. Our findings are incorporated into management policy to guide conservation planning. Our online modeling application provides a template for managers of other systems challenged by EIDs.

Conservation Letters↗

USGS2016

Radiometric dating of marine-influenced coal using Re–Os geochronology

Coal deposits are integral to understanding the structural evolution and thermal history of sedimentary basins and correlating contemporeous estuarine and fluvial delatic strata with marine sections. While marine shales may readily lend themselves to Re–Os dating due to the dominance of hydrogenous Re and Os, the lack of a chronometer for near-shore sedimentary environments hampers basinwide correlations in absolute time. Here, we employ the Re–Os geochronometer, along with total organic carbon (TOC) and Rock–Eval data, to determine the timing and conditions of a marine incursion at the top of the Matewan coal bed, Kanawha Formation, Pottsville Group, West Virginia, USA. The observed range for hydrogen index (HI: 267–290 mg hydrocarbon/gram total organic carbon) for these coal samples suggests dominance of aliphatic hydrocarbons with low carbon (<C 19 ) chain length. Average Re ( 107.6±16.4 ng/g ) and Os ( 0.52±0.09 ng/g ) concentrations of the marine-influenced Matewan coal are higher by few orders of magnitude than published data for terrestrial coal. A Re–Os isochron for the Matewan coal provides an age of 325±14 Ma (Model 3; MSWD = 12; n=19 ; 2 σ ). This is the first Re–Os age derived from coal samples; the age overlaps a new composite Re–Os age of 317±2 Ma for the immediately overlying Betsie Shale Member. External precision for replicate Os analyses carried out for several Matewan coal samples shows a positive correlation with their HI. The HI, which is low in terrestrial organic matter, reflects the degree of marine influence. Thus, samples with the most profound marine influence also have the best analytical reproducibility. Equilibration of Os isotopes with seawater under marine conditions overwhelms variability inherited from terrestrial plant debris, decreasing scatter on the isochron. The 187 Re/ 188 Os ratios of the Matewan coal (&sim;3300–5135) are higher than most of those previously published for Phanerozoic black shale (mostly <2000). Mass balance calculations based on Re/TOC and Os/TOC ratios for the Matewan coal indicate that both Re and Os are primarily marine in origin, and their high 187 Re/ 188 Os ratios confirm efficient removal of both elements from a sulfidic water column into the coal. We show that Re–Os geochronology of marine-influenced coal can be a viable tool for constraining depositional ages.

West Virginia↗