Search USGS⌕ Search

SEARCH · Search USGS

Results for “Journal of Applied Statistics”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 163 records · Page 9Linked to original sources
USGS2000

Stress resistance strategy in an arid land shrub: interactions between developmental instability and fractal dimention

This paper investigates allocation of energy to mechanisms that generate and preserve architectural forms (i.e. developmental stability, complexity of branching patterns) and productivity (growth and reproduction) in response to environmental disturbances (i.e. grazing and resource availability). The statistical error in translational symmetry was used to detect random intra-individual variability during development. This can be thought of as a measure of developmental instability caused by stress. Additionally, we use changes in fractal complexity and shoot distribution of branch structures as an alternate indicator of stress. These methods were applied to Anthyllis cytisoides L., a semi-arid environment shrub, to ascertain the effect of grazing and slope exposure on developmental traits in a 2×2 factorial design. The results show that A. cytisoides maintains developmental stability at the expense of productivity. Anthyllis cytisoides was developmentally more stable when grazed and when on south-facing, as opposed to north-facing slopes. On the contrary, shoot length, leaf area, fractal dimension and reproductive-to-vegetative allocation ratio were larger in north- than in south-facing slopes. As a consequence, under extreme xeric conditions, shrub mortality increased in north-facing slopes, especially when not grazed. The removal of transpiring area and the reduction of plant competition favoured developmental stability and survival in grazed plants. Differences between grazed and ungrazed plants were most evident in more mesic (north-facing) areas.

Journal of Fish Diseases↗
USGS2012

Curie temperature isotherm analysis and tectonic implications of aeromagnetic data from Nevada

Estimates of the depth to the Curie temperature isotherm in Nevada are in accordance with other regional geologic and geophysical information and together can be explained in the context of present-day tectonism. A method to estimate the depth extent of magnetic sources from the statistical properties of magnetic anomalies was applied to a statewide compilation of aeromagnetic data from Nevada. Basal depths of magnetic sources show no apparent correlation with the so-called magnetic quiet zone, which trends northerly through the eastern part of the state, or with basin-and-range topography. However, certain correlations with published heat flow measurements are apparent and suggest that undulations in basal depth of magnetic sources are related in part to undulations in the Curie temperature isotherm. For example, an area of shallow basal depth (<10 km) near Battle Mountain corresponds to an area of exceptionally high conductive heat flow and indicates a shallow depth to the Curie temperature isotherm in this region. A narrow zone of shallow basal depth extends south from the Battle Mountain area along the 118°W meridian to at least latitude 38°N, which also is a zone of historic surface offsets and high-magnitude earthquakes. The correspondence along the 118° meridian of shallow basal depth, high heat flow, high lower crustal seismic velocities, attenuated P and S wave arrivals, historic faulting, and large earthquakes suggests that they each are related to an active north trending spreading zone in this part of the Basin and Range province.

Journal of Geophysical Research Solid Earth↗
USGS2014

A bootstrap estimation scheme for chemical compositional data with nondetects

The bootstrap method is commonly used to estimate the distribution of estimators and their associated uncertainty when explicit analytic expressions are not available or are difficult to obtain. It has been widely applied in environmental and geochemical studies, where the data generated often represent parts of whole, typically chemical concentrations. This kind of constrained data is generically called compositional data, and they require specialised statistical methods to properly account for their particular covariance structure. On the other hand, it is not unusual in practice that those data contain labels denoting nondetects, that is, concentrations falling below detection limits. Nondetects impede the implementation of the bootstrap and represent an additional source of uncertainty that must be taken into account. In this work, a bootstrap scheme is devised that handles nondetects by adding an imputation step within the resampling process and conveniently propagates their associated uncertainly. In doing so, it considers the constrained relationships between chemical concentrations originated from their compositional nature. Bootstrap estimates using a range of imputation methods, including new stochastic proposals, are compared across scenarios of increasing difficulty. They are formulated to meet compositional principles following the log-ratio approach, and an adjustment is introduced in the multivariate case to deal with nonclosed samples. Results suggest that nondetect bootstrap based on model-based imputation is generally preferable. A robust approach based on isometric log-ratio transformations appears to be particularly suited in this context. Computer routines in the R statistical programming language are provided.

Journal of Chemometrics↗
USGS2023

Identification of seismo-volcanic regimes at Whakaari/White Island (New Zealand) via systematic tuning of an unsupervised classifier

We present an algorithm based on Self-Organizing Maps (SOM) and k-means clustering to recognize patterns in a continuous 12.5-year tremor time series recorded at Whakaari/White Island volcano, New Zealand (hereafter referred to as Whakaari). The approach is extendable to a variety of volcanic settings through systematic tuning of the classifier. Hyperparameters are evaluated by statistical means, yielding a combination of “ideal” SOM parameters for the given data set. Extending from this, we applied a Kernel Density Estimation approach to automatically detect changes within the observed seismicity. We categorize the Whakaari seismic time series into regimes representing distinct volcano-seismic states during recent unrest episodes at Whakaari (2012/2013, 2016, and 2019). There is a clear separation in classification results between background regimes and those representing elevated levels of unrest. Onset of unrest is detected by the classifier 6 weeks before the August 2012 eruption, and ca. 3.5 months before the December 2019 eruption, respectively. Regime changes are corroborated by changes in commonly monitored tremor proxies as well as with reported volcanic activity. The regimes are hypothesized to represent diverse mechanisms including: system pressurization and depressurization, degassing, and elevated surface activity. Labeling these regimes improves visualization of the 2012/2013 and 2019 unrest and eruptive episodes. The pre-eruptive 2016 unrest showed a contrasting shape and nature of seismic regimes, suggesting differing onset and driving processes. The 2016 episode is proposed to result from rapid destabilization of the shallow hydrothermal system, while rising magmatic gases from new injections of magma better explain the 2012/2013 and 2019 episodes.

Whakaari, White Island↗
USGS2013

Potential climate change impacts on temperate forest ecosystem processes

Large changes in atmospheric CO2, temperature and precipitation are predicted by 2100, yet the long-term consequences for carbon, water, and nitrogen cycling in forests are poorly understood. We applied the PnET-CN ecosystem model to compare the long-term effects of changing climate and atmospheric CO2 on productivity, evapotranspiration, runoff, and net nitrogen mineralization in current Great Lakes forest types. We used two statistically downscaled climate projections, PCM B1 (warmer and wetter) and GFDL A1FI (hotter and drier), to represent two potential future climate and atmospheric CO2 scenarios. To separate the effects of climate and CO2, we ran PnET-CN including and excluding the CO2 routine. Our results suggest that, with rising CO2 and without changes in forest type, average regional productivity could increase from 67% to 142%, changes in evapotranspiration could range from –3% to +6%, runoff could increase from 2% to 22%, and net N mineralization could increase 10% to 12%. Ecosystem responses varied geographically and by forest type. Increased productivity was almost entirely driven by CO2 fertilization effects, rather than by temperature or precipitation (model runs holding CO2 constant showed stable or declining productivity). The relative importance of edaphic and climatic spatial drivers of productivity varied over time, suggesting that productivity in Great Lakes forests may switch from being temperature to water limited by the end of the century.

Michigan;Minnesota;Wisconsin↗
USGS2009

Survival and breeding of polar bears in the southern Beaufort Sea in relation to sea ice

1. Observed and predicted declines in Arctic sea ice have raised concerns about marine mammals. In May 2008, the US Fish and Wildlife Service listed polar bears (Ursus maritimus) - one of the most ice-dependent marine mammals - as threatened under the US Endangered Species Act. 2. We evaluated the effects of sea ice conditions on vital rates (survival and breeding probabilities) for polar bears in the southern Beaufort Sea. Although sea ice declines in this and other regions of the polar basin have been among the greatest in the Arctic, to date population-level effects of sea ice loss on polar bears have only been identified in western Hudson Bay, near the southern limit of the species' range. 3. We estimated vital rates using multistate capture-recapture models that classified individuals by sex, age and reproductive category. We used multimodel inference to evaluate a range of statistical models, all of which were structurally based on the polar bear life cycle. We estimated parameters by model averaging, and developed a parametric bootstrap procedure to quantify parameter uncertainty. 4. In the most supported models, polar bear survival declined with an increasing number of days per year that waters over the continental shelf were ice free. In 2001-2003, the ice-free period was relatively short (mean 101 days) and adult female survival was high (0 ∙ 96-0 ∙ 99, depending on reproductive state). In 2004 and 2005, the ice-free period was longer (mean 135 days) and adult female survival was low (0 ∙ 73-0 ∙ 79, depending on reproductive state). Breeding rates and cub litter survival also declined with increasing duration of the ice-free period. Confidence intervals on vital rate estimates were wide. 5. The effects of sea ice loss on polar bears in the southern Beaufort Sea may apply to polar bear populations in other portions of the polar basin that have similar sea ice dynamics and have experienced similar, or more severe, sea ice declines. Our findings therefore are relevant to the extinction risk facing approximately one-third of the world's polar bears.

Journal of Animal Ecology↗
USGS2015

Effects of freezing on white perch Morone americana (Gmelin, 1789): Implications for multivariate morphometrics

This study tested the hypothesis that duration of freezing differentially affects whole-body morphometrics of a derived teleost. Whole-body morphometrics are frequently analyzed to test hypotheses of different species, or stocks within a species, of fishes. Specimens used for morphometric analyses are typically fixed or preserved prior to analysis, yet little research has been done on how fixation or preservation methods or duration of preservation of specimens might affect outcomes of multivariate statistical analyses of differences in shape. To determine whether whole-body morphometrics changed as a result of freezing, 23 whole-body morphometrics of age-1 white perch ( Morone americana ) from western Lake Erie (n = 211) were analyzed immediately after capture, after being held on ice overnight, and after freezing for 100 or 200 days. Discriminant function analysis revealed that all four groups differed significantly from one another (P < 0.0001). The first canonical axis reflected long-axis morphometrics, where there was a clear pattern of positive translation along this axis with duration of preservation. Re-classification analysis demonstrated fish were typically assigned to their original preservation class except for fish frozen 100 days, which assigned mostly to frozen 200 days. Morphometric comparisons using frozen fish must be done on fish frozen for identical periods of time to avoid biases related to the length of time they were frozen. Similar experiments should be conducted on other species and also using formalin- and alcohol-preserved specimens.

Western Lake Erie↗
USGS2012

Tools for quantifying isotopic niche space and dietary variation at the individual and population level.

Ecologists are increasingly using stable isotope analysis to inform questions about variation in resource and habitat use from the individual to community level. In this study we investigate data sets from 2 California sea otter ( Enhydra lutris nereis ) populations to illustrate the advantages and potential pitfalls of applying various statistical and quantitative approaches to isotopic data. We have subdivided these tools, or metrics, into 3 categories: IsoSpace metrics, stable isotope mixing models, and DietSpace metrics. IsoSpace metrics are used to quantify the spatial attributes of isotopic data that are typically presented in bivariate (e.g., δ 13 C versus δ 15 N) 2-dimensional space. We review IsoSpace metrics currently in use and present a technique by which uncertainty can be included to calculate the convex hull area of consumers or prey, or both. We then apply a Bayesian-based mixing model to quantify the proportion of potential dietary sources to the diet of each sea otter population and compare this to observational foraging data. Finally, we assess individual dietary specialization by comparing a previously published technique, variance components analysis, to 2 novel DietSpace metrics that are based on mixing model output. As the use of stable isotope analysis in ecology continues to grow, the field will need a set of quantitative tools for assessing isotopic variance at the individual to community level. Along with recent advances in Bayesian-based mixing models, we hope that the IsoSpace and DietSpace metrics described here will provide another set of interpretive tools for ecologists.

Journal of Mammalogy↗
USGS2001

Processes controlling the episodic streamwater transport of atrazine and other agrichemicals in an agricultural watershed

Episodic streamwater transport of atrazine (a common agricultural herbicide) and nutrients has been observed throughout agricultural watersheds in the United States and poses a serious threat to the quality of its water resources. Catchment-scale atrazine and nutrient transport processes after agricultural application are still poorly understood, and predicting episodic streamwater composition remains an elusive goal. We instrumented a 1.2-km 2 agricultural catchment near Harrisonburg, Virginia, and examined streamwater, overland flow, soil water, groundwater, and rainfall during the summer of 1998. Storm chemographs demonstrated different patterns for constituents derived primarily from weathering (silica and calcium), compared to constituents derived primarily from early spring land applications (nitrate, atrazine, DOC, potassium, chloride, and sulfate). During storms, the concentrations of silica and calcium decreased, the atrazine response was variable, and the concentrations of nitrate, DOC, potassium, chloride, and sulfate increased; the elevated nitrate signal lagged several hours behind the other elevated constituents. Graphical and statistical analyses indicated a relatively stable spring-fed baseflow was modified by a mixture of overland flow and soil water. A rapid, short-duration overland-flow pulse dominated the streamflow early in the event and contributed most of the potassium, DOC, chloride, suspended sediment, and atrazine. A longer-duration soil–water pulse dominated the streamflow later in the event and contributed the nitrate as well as additional potassium, DOC, sulfate, and atrazine. The contributions to the episodic streamflow were quantified using a flushing model in which overland-flow and soil–water concentrations decreased exponentially with time during an episode. Flushing time constants for the overland-flow and soil–water reservoirs were calculated on a storm-by-storm basis using separate tracers for each time-variable reservoir. Initial component concentrations were estimated through regression analyses. Mass-balance calculations were used for flow separations and to predict the observed streamwater composition. Model forecasts indicated that reduced fertilizer and pesticide application (rather than elimination of overland-flow or soil–water contributions) was necessary to improve the episodic streamwater composition. This study provides important additional understanding of the catchment-scale processes by which land-applied pesticides and nutrients can move through agricultural systems.

Virginia↗
USGS2015

Implications of the methodological choices for hydrologic portrayals of climate change over the contiguous United States: Statistically downscaled forcing data and hydrologic models

Continental-domain assessments of climate change impacts on water resources typically rely on statistically downscaled climate model outputs to force hydrologic models at a finer spatial resolution. This study examines the effects of four statistical downscaling methods [bias-corrected constructed analog (BCCA), bias-corrected spatial disaggregation applied at daily (BCSDd) and monthly scales (BCSDm), and asynchronous regression (AR)] on retrospective hydrologic simulations using three hydrologic models with their default parameters (the Community Land Model, version 4.0; the Variable Infiltration Capacity model, version 4.1.2; and the Precipitation–Runoff Modeling System, version 3.0.4) over the contiguous United States (CONUS). Biases of hydrologic simulations forced by statistically downscaled climate data relative to the simulation with observation-based gridded data are presented. Each statistical downscaling method produces different meteorological portrayals including precipitation amount, wet-day frequency, and the energy input (i.e., shortwave radiation), and their interplay affects estimations of precipitation partitioning between evapotranspiration and runoff, extreme runoff, and hydrologic states (i.e., snow and soil moisture). The analyses show that BCCA underestimates annual precipitation by as much as −250 mm, leading to unreasonable hydrologic portrayals over the CONUS for all models. Although the other three statistical downscaling methods produce a comparable precipitation bias ranging from −10 to 8 mm across the CONUS, BCSDd severely overestimates the wet-day fraction by up to 0.25, leading to different precipitation partitioning compared to the simulations with other downscaled data. Overall, the choice of downscaling method contributes to less spread in runoff estimates (by a factor of 1.5–3) than the choice of hydrologic model with use of the default parameters if BCCA is excluded.

Journal of Hydrometeorology↗
USGS2010

On the resolution of shallow mantle viscosity structure using post-earthquake relaxation data: Application to the 1999 Hector Mine, California, earthquake

Most models of lower crust/mantle viscosity inferred from postearthquake relaxation assume one or two uniform-viscosity layers. A few existing models possess apparently significant radially variable viscosity structure in the shallow mantle (e.g., the upper 200 km), but the resolution of such variations is not clear. We use a geophysical inverse procedure to address the resolving power of inferred shallow mantle viscosity structure using postearthquake relaxation data. We apply this methodology to 9 years of GPS-constrained crustal motions after the 16 October 1999 M = 7.1 Hector Mine earthquake. After application of a differencing method to isolate the postearthquake signal from the “background” crustal velocity field, we find that surface velocities diminish from ∼20 mm/yr in the first few months to ≲2 mm/yr after 2 years. Viscoelastic relaxation of the mantle, with a time-dependent effective viscosity prescribed by a Burgers body, provides a good explanation for the postseismic crustal deformation, capturing both the spatial and temporal pattern. In the context of the Burgers body model (which involves a transient viscosity and steady state viscosity), a resolution analysis based on the singular value decomposition reveals that at most, two constraints on depth-dependent steady state mantle viscosity are provided by the present data set. Uppermost mantle viscosity (depth ≲ 60 km) is moderately resolved, but deeper viscosity structure is poorly resolved. The simplest model that explains the data better than that of uniform steady state mantle viscosity involves a linear gradient in logarithmic viscosity with depth, with a small increase from the Moho to 220 km depth. However, the viscosity increase is not statistically significant. This suggests that the depth-dependent steady state viscosity is not resolvably different from uniformity in the uppermost mantle.

California↗
USGS2021

Racial, ethnic, and social patterns in the recreation specialization of birdwatchers: An analysis of United States eBird registrants

Although birdwatchers comprise a large and growing proportion of the American public, there is a lack of racial and ethnic diversity in the birdwatching community. Previous research suggests that this homogeneity is self-perpetuating, as ethno-racial minorities are less likely to pursue activities in which no one they know participates. However, it is unclear whether this trend in birdwatching participation also applies to degree of subsequent participant involvement. Using a national online survey of US birdwatchers, we measured the degree of recreation specialization among birdwatchers along affective, cognitive, and behavioral dimensions. We also determined whether respondents had social connections (acquaintances, close friends, or relatives) who birdwatch. We then used logistic regression to determine which ethno-racial groups were more likely to have birdwatcher social connections, and multiple linear regression to investigate how our measures of recreation specialization varied by ethno-racial group. As expected, the ethno-racial composition of the birdwatchers we studied was significantly less diverse than that of the American public. Of the 29,380 respondents who reported their ethno-racial group, 5.2% were Black, Indigenous, and/or people of color (including Native American, Black, Asian, Pacific Islander, Hispanic/Latino, or multiracial), while 94.8% were non-Hispanic White. However, we observed no statistically significant ethno-racial patterns in overall degree of recreation specialization, even when controlling for social connection and demographic characteristics. Considering the three dimensions of specialization individually, we found that some ethno-racial predictors were statistically significant, but coefficients were too small to be practically significant. We conclude that while some ethno-racial groups are underrepresented among birdwatchers, there is insufficient evidence that they are also under-specialized.

Journal of Outdoor Recreation and Tourism↗
USGS2022

NABat ML: Utilizing deep learning to enable crowdsourced development of automated, scalable solutions for documenting North American bat populations

Bats play crucial ecological roles and provide valuable ecosystem services, yet many populations face serious threats from various ecological disturbances. The North American Bat Monitoring Program (NABat) aims to use its technology infrastructure to assess status and trends of bat populations, while developing innovative and community-driven conservation solutions. Here, we present NABat ML , an automated machine-learning algorithm that improves the scalability and scientific transparency of NABat acoustic monitoring. This model combines signal processing techniques and convolutional neural networks (CNNs) to detect and classify recorded bat echolocation calls. We developed our CNN model with internet-based computing resources (‘cloud environment’), and trained it on >600,000 spectrogram images. We also incorporated species range maps to improve the robustness and accuracy of the model for future ‘unseen’ data. We evaluated model performance using a comprehensive, independent, holdout dataset. NABat ML successfully distinguished 31 classes (30 species and a noise class) with overall weighted-average accuracy and precision rates of 92%, and ≥90% classification accuracy for 19 of the bat species. Using a single cloud-environment computing instance, the entire model training process took <16 h. Synthesis and applications . Our convolutional neural network (CNN)-based model, NABat ML , classifies 30 North American bat species using their recorded echolocation calls with an overall accuracy of 92%. In addition to providing highly accurate species-level classification, NABat ML and its outputs are compatible with Bayesian and other statistical techniques for measuring uncertainty in classification. Our model is open-source and reproducible, enabling future implementations as software on end-user devices and cloud-based web applications. These qualities make NABat ML highly suitable for applications ranging from grassroots community science initiatives to big-data methods developed and implemented by researchers and professional practitioners. We believe the transparency and accessibility of NABat ML will encourage broad-scale participation in bat monitoring, and enable development of innovative solutions needed to conserve North American bat species.

Journal of Applied Ecology↗
USGS2015

Effects of climate and land cover on hydrology in the southeastern U.S.: Potential impacts on watershed planning

The hydrologic response to statistically downscaled general circulation model simulations of daily surface climate and land cover through 2099 was assessed for the Apalachicola-Chattahoochee-Flint River Basin located in the southeastern United States. Projections of climate, urbanization, vegetation, and surface-depression storage capacity were used as inputs to the Precipitation-Runoff Modeling System to simulate projected impacts on hydrologic response. Surface runoff substantially increased when land cover change was applied. However, once the surface depression storage was added to mitigate the land cover change and increases of surface runoff (due to urbanization), the groundwater flow component then increased. For hydrologic studies that include projections of land cover change (urbanization in particular), any analysis of runoff beyond the change in total runoff should include effects of stormwater management practices as these features affect flow timing and magnitude and may be useful in mitigating land cover change impacts on streamflow. Potential changes in water availability and how biota may respond to changes in flow regime in response to climate and land cover change may prove challenging for managers attempting to balance the needs of future development and the environment. However, these models are still useful for assessing the relative impacts of climate and land cover change and for evaluating tradeoffs when managing to mitigate different stressors.

Alabama, Florida, Georgia↗
USGS2015

Microscopical characterization of carbon materials derived from coal and petroleum and their interaction phenomena in making steel electrodes, anodes and cathode blocks for the Microscopy of Carbon Materials Working Group of the ICCP

This paper describes the evaluation of petrographic textures representing the structural organization of the organic matter derived from coal and petroleum and their interaction phenomena in the making of steel electrodes, anodes and cathode blocks. This work represents the results of the Microscopy of Carbon Materials Working Group in Commission III of the International Committee for Coal and Organic Petrology between the years 2009 and 2013. The round robin exercises were run on photomicrograph samples. For textural characterization of carbon materials the existing ASTM classification system for metallurgical coke was applied. These round robin exercises involved 15 active participants from 12 laboratories who were asked to assess the coal and petroleum based carbons and to identify the morphological differences, as optical texture (isotropic/anisotropic), optical type (punctiform, mosaic, fibre, ribbon, domain), and size. Four sets of digital black and white microphotographs comprising 151 photos containing 372 fields of different types of organic matter were examined. Based on the unique ability of carbon to form a wide range of textures, the results showed an increased number of carbon occurrences which have crucial role in the chosen industrial applications. The statistical method used to evaluate the results was based on the “raw agreement indices”. It gave a new and original view on the analysts' opinion by not only counting the correct answers, but also all of the knowledge and experience of the participants. Comparative analyses of the average values of the level of overall agreement performed by each analyst in the exercises during 2009–2013 showed a great homogeneity in the results, the mean value being 90.36%, with a minimum value of 83% and a maximum value of 95%.

International Journal of Coal Geology↗
USGS2009

Disease, predation and demography: Assessing the impacts of bovine tuberculosis on African buffalo by monitoring at individual and population levels

Summary Understanding the effects of disease is critical to determining appropriate management responses, but estimating those effects in wildlife species is challenging. We used bovine tuberculosis (BTB) in the African buffalo Syncerus caffer population of Kruger National Park, South Africa, as a case study to highlight the issues associated with estimating chronic disease effects in a long‐lived host. We used known and radiocollared buffalo, aerial census data, and a natural gradient in pathogen prevalence to investigate if: (i) at the individual level, BTB infection reduces reproduction; (ii) BTB infection increases vulnerability to predation; and (iii) at the population level, increased BTB prevalence causes reduced population growth. There was only a marginal reduction in calving success associated with BTB infection, as indexed by the probability of sighting a known adult female with or without a calf ( P = 0·065). Since 1991, BTB prevalence increased from 27 to 45% in the southern region and from 4 to 28% in the central region of Kruger National Park. The prevalence in the northern regions was only 1·5% in 1998. Buffalo population growth rates, however, were neither statistically different among regions nor declining over time. Lions Panthera leo did not appear to preferentially kill test‐positive buffalo. The best (Akaike's Information Criterion corrected for small sample size) AIC c model with BTB as a covariate [exp(β) = 0·49; 95% CI = (0·24–1·02)] suggested that the mortality hazard for positive individuals was no greater than for test‐negative individuals. Synthesis and applications . Test accuracy, time‐varying disease status, and movement among populations are some of the issues that make the detection of chronic disease impacts challenging. For these reasons, the demographic impacts of bovine tuberculosis in the Kruger National Park remain undetectable despite 6 years of study on known individuals and 40 years of population counts. However, the rainfall and forage conditions during this study were relatively good and the impacts of many chronic diseases may be a non‐linear function of environmental conditions such that they are only detectable in stressful periods.

Kruger National Park↗
USGS2005

Applying petrophysical models to radar travel time and electrical resistivity tomograms: Resolution-dependent limitations

[1] Geophysical imaging has traditionally provided qualitative information about geologic structure; however, there is increasing interest in using petrophysical models to convert tomograms to quantitative estimates of hydrogeologic, mechanical, or geochemical parameters of interest (e.g., permeability, porosity, water content, and salinity). Unfortunately, petrophysical estimation based on tomograms is complicated by limited and variable image resolution, which depends on (1) measurement physics (e.g., electrical conduction or electromagnetic wave propagation), (2) parameterization and regularization, (3) measurement error, and (4) spatial variability. We present a framework to predict how core‐scale relations between geophysical properties and hydrologic parameters are altered by the inversion, which produces smoothly varying pixel‐scale estimates. We refer to this loss of information as “correlation loss.” Our approach upscales the core‐scale relation to the pixel scale using the model resolution matrix from the inversion, random field averaging, and spatial statistics of the geophysical property. Synthetic examples evaluate the utility of radar travel time tomography (RTT) and electrical‐resistivity tomography (ERT) for estimating water content. This work provides (1) a framework to assess tomograms for geologic parameter estimation and (2) insights into the different patterns of correlation loss for ERT and RTT. Whereas ERT generally performs better near boreholes, RTT performs better in the interwell region. Application of petrophysical models to the tomograms in our examples would yield misleading estimates of water content. Although the examples presented illustrate the problem of correlation loss in the context of near‐surface geophysical imaging, our results have clear implications for quantitative analysis of tomograms for diverse geoscience applications.

Journal of Geophysical Research B: Solid Earth↗
USGS2025

Concentration-discharge relations and transient metal loads reveal spatiotemporal variability in solute-generation mechanisms in a mine-affected watershed

Concentration-discharge (CQ) relations are commonly used to understand geochemical and hydrologic controls on the generation of solutes in watersheds. Despite the widespread application of CQ relations, this technique has been infrequently applied to acid mine drainage (AMD) sites, but the CQ framework may allow mechanistic understanding of remedial outcomes such as impoundment of water within underground mines. Results of CQ analyses and changes in metal loads in an AMD affected watershed in Colorado, USA indicate that dissolved loads increased at many individual locations following water impoundment within mine workings. Although increased loads were observed at most individual locations, these increases were offset by a large decrease in loading from the largest mine. A loading analysis that included data from an instream monitoring location showed a statistically significant decrease in Fe and Zn after bulkhead emplacement, indicating a net positive effect of bulkheads. Streams generally displayed dilution CQ patterns whereas mines and springs showed either flushing or chemostatic patterns prior to bulkheading, which transitioned to chemostatic patterns following bulkheading, indicating a transition from dynamic to equilibrium geochemical processes. Saturation indices for sulfide and secondary minerals indicated that mines and springs were near equilibrium for phases including schwertmannite, fluorite, and gypsum. Saturation indices vary through time for mines suggesting progressive leaching of sulfide minerals as the mass of available minerals in the mine workings decreases. Together, these diverse analyses provide an integrated understanding of the variability in solute generating processes in this watershed and may inform remediation plans for similarly affected sites by indicating the nature of mineralogic controls on water quality.

Colorado↗