USGS Science⌕ Search

SEARCH · USGS Science

Results for “Spatial 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 1,081 records · Page 60Linked to original sources

The Everglades Depth Estimation Network (EDEN) surface-water interpolation model, version 3

The Everglades Depth Estimation Network (EDEN) is an integrated network of water-level gages, interpolation models that estimate daily water-level data at ungaged locations, and applications that generate derived hydrologic data across the freshwater part of the Greater Everglades landscape. Version 3 (V3) of the EDEN interpolation surface-water model is the most recent update, replacing the version 2 (V2) model released in 2011. The primary revision for the V3 model is the switch to the R programming language to create a more efficient and portable EDEN code relative to V2, without reliance on proprietary software. Using R, the interpolation script runs over 10 times faster and is more easily updated, for example, to accommodate changes in the gage network or to incorporate R software updates. Additional revisions made for the V3 model include updates to the interpolation model, the gage network, and groundwater-level estimations. The EDEN model domain in the Greater Everglades and Big Cypress National Preserve is divided into subdomains that are based on hydrologic boundaries. In the V3 model, the number of subdomains was increased from five to eight, which allows hydrologic boundaries, such as levees and canals, to be better represented in the interpolation scheme. Five pseudogages were added to constrain the water-level surface at subdomain boundaries. Changes made to the water-level gage network between the implementation of the V2 and V3 models are incorporated, and groundwater-level estimations are added, which are important information for hydrologic and ecological studies. Summary model performance statistics indicate similar accuracy in water-level surfaces generated by the V3 and V2 models, with a root mean square error of 4.78 centimeters for both interpolation models against independent water-level measurements. Providing stability and continuity for the EDEN user community, the V3 model closely replicates the V2 model, with a root mean square difference of 3.87 centimeters for interpolated surfaces from April 1, 2014, to March 31, 2018. The additional groundwater levels provide a realistic estimate of the saturated groundwater surface continuous with the surface-water surface for Water Conservation Areas 2A and 2B from 2000 to 2011. This continuous surface is a more accurate estimation of the spatial distribution of water in the hydrologic system than before, providing needed information for ecological studies in areas where depth to water table affects habitats. Development of the EDEN V3 model advances the tools available to scientists and resource managers for guiding large-scale field operations, describing hydrologic changes, and supporting biological and ecological assessments.

Florida↗

Temporal expansion of annual crop classification layers for the CONUS using the C5 decision tree classifier

Crop cover maps have become widely used in a range of research applications. Multiple crop cover maps have been developed to suite particular research interests. The National Agricultural Statistics Service (NASS) Cropland Data Layers (CDL) are a series of commonly used crop cover maps for the conterminous United States (CONUS) that span from 2008 to 2013. In this investigation, we sought to contribute to the availability of consistent CONUS crop cover maps by extending temporal coverage of the NASS CDL archive back eight additional years to 2000 by creating annual NASS CDL-like crop cover maps derived from a classification tree model algorithm. We used over 11 million records to train a classification tree algorithm and develop a crop classification model (CCM). The model was used to create crop cover maps for the CONUS for years 2000–2013 at 250 m spatial resolution. The CCM and the maps for years 2008–2013 were assessed for accuracy relative to resampled NASS CDLs. The CCM performed well against a withheld test data set with a model prediction accuracy of over 90%. The assessment of the crop cover maps indicated that the model performed well spatially, placing crop cover pixels within their known domains; however, the model did show a bias towards the ‘Other’ crop cover class, which caused frequent misclassifications of pixels around the periphery of large crop cover patch clusters and of pixels that form small, sparsely dispersed crop cover patches.

Remote Sensing Letters↗

Assessment of water availability in the Osage Nation using an integrated hydrologic-flow model

The Osage Nation of northeastern Oklahoma, conterminous with Osage County, covers about 2,900 square miles. The area is primarily rural with 62 percent of the land being native prairie grass, and much of the area is used for cattle ranching and extraction of petroleum and natural gas. Protection of water rights are important to the Osage Nation because of its reliance on cattle ranching and the potential for impairment of water quality by petroleum extraction. Additionally, the potential for future population increases, demands for water from neighboring areas such as the Tulsa metropolitan area, and expansion of petroleum and natural-gas extraction on water resources of this area further the need for the Osage Nation to better understand its water availability. Therefore, the U.S. Geological Survey, in cooperation with the Osage Nation, completed a hydrologic investigation to assess the status and availability of surface-water and groundwater resources in the Osage Nation. A transient integrated hydrologic-flow model was constructed using the U.S. Geological Survey fully integrated hydrologic-flow model called the MODFLOW One-Water Hydrologic Model. The integrated hydrologic-flow model, called the Osage Nation Integrated Hydrologic Model (ONIHM), was constructed and uses an orthogonal grid of 276 rows and 289 columns, and each grid cell measures 1,312.34 feet (ft; 400 meters) per side, with eight variably thick vertical layers that represented the alluvial and bedrock aquifers within the study area, including the alluvial aquifer, the Vamoosa-Ada aquifer, and the minor Pennsylvanian bedrock aquifers, and the confining units. Landscape and groundwater-flow processes were simulated for two periods: (1) the 1950–2014 period from January 1950 through September 2014 and (2) the forecast period from October 2014 through December 2099. The 1950–2014 period ONIHM simulated past conditions using measured or estimated inputs, and the forecast-period ONIHM simulated three separate potential forecast conditions under constant dry, average, or wet climate conditions using calibrated input values from the 1950–2014 period ONIHM. The 1950–2014 period ONIHM was calibrated by linking the Parameter Estimation software (PEST) with the MODFLOW One-Water Hydrologic Model. PEST uses statistical parameter estimation techniques to identify the best set of parameter values to minimize the difference between measured or estimated calibration targets and their simulated equivalent values (residuals). Tikhonov regularization and singular-value decomposition-assist features of PEST were used during the calibration process. The 1950–2014 period ONIHM was calibrated to 713 measured groundwater levels at 195 wells; 95,636 estimated monthly mean groundwater levels at 124 wells; 5,307 measured streamflows at 13 streamgages; and 8,679 simulated mean monthly streamflows at 10 streamgages extracted from a surface-water model by adjusting 231 parameters. The estimated groundwater-level observations and streamflows were included as observations to improve the spatial and temporal density of observation targets during calibration. The best set of parameter values obtained during the calibration process of the 1950–2014 model was then used as the input parameter values for the forecast model simulations. A comparison of the calibration targets to their corresponding simulated values indicated that the model adequately reproduced streamflows and groundwater levels for some streamgages and wells and underestimated streamflows and groundwater levels at other locations. Measured and simulated streamflows correlated adequately with a coefficient of determination of 0.938, as did water levels with a coefficient of determination of 0.795. The 1950–2014 period ONIHM underestimated certain groundwater levels and streamflows, but generally measured or estimated calibration targets correlated well with simulated equivalents, which indicated that the model can adequately simulate the response of the hydrologic system to stresses in the 1950–2014 and forecast periods. In the 1950–2014 period ONIHM, the calibrated mean horizontal hydraulic conductivity for layer 1 alluvial aquifer was 30.7 feet per day, and the seven lower layers had a calibrated mean horizontal hydraulic conductivity of less than 3.3 feet per day. The mean calibrated groundwater-level residual was 16.6 ft, and the mean calibrated streamflow residual of the Arkansas River at Ralston, Oklahoma, streamgage (U.S. Geological Survey station 07152500) was within 6 percent (373 cubic feet per second) of mean measured streamflow for the 1950–2014 period ONIHM. The ONIHM simulated landscape fluxes of precipitation; groundwater applied by irrigation wells; evapotranspiration from precipitation, groundwater, and irrigation; runoff from precipitation; and deep percolation from precipitation. The largest loss of water from the landscape was evapotranspiration from precipitation with a calibrated mean annual outflow of 32 inches (in.): mean annual precipitation was about 36 in. Calibrated mean annual runoff and deep percolation (recharge to the water table) rates were 4.7 inches per year (in/yr) and 0.70 in/yr, respectively, for the 1950–2014 period ONIHM. The calibrated 1950–2014 period ONIHM groundwater fluxes included net farm net recharge (calculated as the difference between the inflow of recharge to the water table and the outflow of evapotranspiration from the water table such that negative values indicate that evapotranspiration from the water table was greater than deep percolation [recharge to the water table] and vice versa). Net farm net recharge was the largest flux from the groundwater system with a mean annual net outflow of 153.4 cubic feet per second. Stream leakage was the largest flux to the groundwater system with a mean annual net inflow of 152.5 cubic feet per second, indicating that, on average, the groundwater/surface-water interaction was a “losing” system where stream water leaked into the subsurface and recharged the water table. Simulated monthly trends demonstrated that net stream leakage was the largest inflow to the groundwater-flow system for 10 of the 12 months; for the other 2 months (January and March), farm net recharge (January) and net storage (March) were the largest inflow to the groundwater-flow system. A saline groundwater interface map was created for the study and compared to the water levels from the final stress period of the 1950–2014 model to identify the presence of fresh/marginal groundwater throughout the study area. Fresh/marginal groundwater was characterized as groundwater with less than 1,500 milligrams per liter of total dissolved solids. Fresh/marginal groundwater thickness ranged from 0 to 438.2 ft within the study area. The thickest regions of fresh/marginal groundwater were in the eastern part of the study area near Sand Creek, Bird Creek, and Hominy Creek and in the Arkansas River alluvial aquifer in the region downstream from the Arkansas River at Ralston, Okla. Like the 1950–2014 model, forecast model results for the landscape indicated that transpiration from precipitation was the largest flux out of the landscape for all three forecasts, constituting 77, 73, and 58 percent of precipitation for the dry, average, and wet forecasts, respectively. The dry and average forecast landscape fluxes demonstrated similar trends and magnitudes, whereas the wet forecast landscape fluxes indicated the largest changes compared to the average forecast fluxes. Most notably, runoff increased from a mean of 1.1 and 1.6 in/yr for the dry and average forecasts, respectively, to 10 in/yr for the wet forecast. Similar changes occurred for the other wet forecast landscape fluxes. The calibrated 1950–2014 period ONIHM simulated three forecasts to assess the effects of potential climatic changes on the hydrologic system from October 2014 to December 2099. The three forecasts simulated theoretical dry, average, and wet conditions using precipitation and potential evapotranspiration datasets from selected years in the calibrated 1950–2014 period ONIHM. Annual precipitation amounts were 26.89, 35.47, and 50.73 in. for the dry, average, and wet forecasts, respectively. Groundwater-flow component forecast results indicated that stream leakage is always a net inflow to the groundwater-flow system for dry, average, and wet conditions, meaning the study area stream network is always predominantly a “losing” regime where stream water infiltrates into the underlying aquifer. Storage was only a net outflow from the groundwater-flow system and indicated a replenishment to groundwater storage that resulted in an increase in groundwater levels only during the wet forecast. Further, these gains in groundwater storage for the wet forecast occurred only during February through June. Mean fresh/marginal groundwater saturated thicknesses were 125 and 126 ft for the dry and average forecast conditions, respectively, and wet forecast average thickness was 145 ft and ranged from 0 to 443 ft. The spatial extents of fresh/marginal groundwater at the end of the dry, average, and wet forecast model periods (December 2099) did not change substantially from the end of the 1950–2014 model period (September 2014).

Kansas, Oklahoma↗

Water-quality data analysis of the upper Gunnison River watershed, Colorado, 1989-99

Water-quality data from October 1969 to December 1999 for both surface water and ground water in the upper Gunnison River watershed were retrieved and compiled from the U.S. Geological Survey National Water Information System and the U.S. Environmental Protection Agency Storage and Retrieval databases. Analyses focused primarily on a subset of these data from October 1989 to December 1999. The upper Gunnison River watershed is located west of the Continental Divide in the Southern Rocky Mountains physiographic province. Surface-water-quality data were compiled for 482 sites in the upper Gunnison River watershed. Most values of surface-water temperature, dissolved oxygen, and pH were within Colorado Department of Public Health and Environment (CDPHE) in-stream standards. Calcium bicarbonate type water was the most spatially dominant water type in the basin. Nutrients were most commonly sampled along the Slate River and East River near Crested Butte and along the Gunnison River from the confluence of the East and Taylor Rivers to the western edge of the watershed. Median ammonia concentrations were low, with many concentrations less than laboratory reporting levels. All nitrate concentrations met the CDPHE in-stream standard of 10 milligrams per liter. More than 30 percent of stream sites with total phosphorus data (23 of 61 sites) had concentrations greater than the U.S. Environmental Protection Agency (USEPA) recommendation for controlling eutrophication. Ammonia concentrations at a site on the Slate River near Crested Butte had a statistically significant upward trend for the 1995–99 period. The Slate River near Crested Butte site is located immediately downstream from the towns of Crested Butte and Mount Crested Butte and may reflect recent population growth or other land-use changes. However, the rate of change of the trend is small (0.017 milligram per liter per year). Although a multiple comparison test showed nitrate concentrations were statistically different between agriculture and forest sites and between agriculture and urban land-use classified sites, median concentrations were low among all land-use settings. Median concentrations of total phosphorus were greatest in rangeland areas and least in urban areas. No significant differences were identified for median concentrations of total phosphorus in agriculture and forest land-use areas. Median concentrations of arsenic, lead, mercury, selenium, and silver were low or below reporting levels throughout the watershed. Aluminum, cadmium, copper, lead, manganese, and zinc concentrations were elevated near the town of Crested Butte and on Henson Creek upstream from Lake City, which may be explained by upstream areas of historical mining. Samples for six trace elements exceeded standards: cadmium, copper, lead, manganese, silver, and zinc. A downward trend (3 micrograms per liter per year) was identified for the dissolved iron concentration at a site on the Gunnison River at County Road 32 downstream from the city of Gunnison. Streambed-sediment samples from areas affected by historical mining also had elevated concentrations of some trace elements. Chlorophyll- a concentrations in samples from Blue Mesa Reservoir and streams in the Crested Butte and Gunnison areas were typical of unenriched to moderately enriched conditions. Median concentrations of 5-day biochemical oxygen demand concentrations for sites between Crested Butte and Blue Mesa Reservoir were less than 2 milligrams per liter. Occasional high (greater than 200 counts per 100 milliliters) concentrations for fecal coliform were determined at selected sites within the study area. However, median concentrations were less than 100 counts per 100 milliliters except for the Squaw Creek and Cimarron River areas in the western part of the watershed. Ground-water-quality data have been collected by the U.S. Geological Survey from 99 wells. Many wells were completed in aquifers composed of Holocene-age valley fill and alluvium. Most field properties were within the USEPA Secondary Drinking Water Regulations (SDWR) range for treated drinking water, except for 2 (of 40) pH samples. Calcium bicarbonate was the predominant water type in nearly all aquifers except for the aquifers composed of volcanic rock, which had more sodium and sulfate mixed water types. Wells with sulfate concentrations exceeding the SDWR of 250 milligrams per liter were completed in aquifers composed of volcanic rock near Lake City. Dissolution and oxidation of sulfide minerals in these aquifers may explain the elevated sulfate concentrations in ground water at these locations. Nutrient concentrations in ground water were generally low, and median concentrations for ammonia, nitrite, and dissolved phosphorus were below reporting levels. All nitrate concentrations in the samples were below the USEPA drinking-water maximum contaminant level of 10 mg/L. No statistical difference was found in nitrate concentrations among the four land-use classifications (agriculture, forest, rangeland, and urban). Trace elements in ground water were generally below the USEPA SDWR. Three iron samples exceeded the USEPA SDWR of 300 micrograms per liter at two wells located near the city of Gunnison and at a well south of the town of Powderhorn near the Cebolla River. Nine of 39 manganese samples exceeded the USEPA SDWR of 50 micrograms per liter and were collected from aquifers composed of Holocene-age valley fill and alluvium near Gunnison and Crested Butte and in one well near the Cebolla River. Radon gas is a natural radioactive decay product of uranium. All 39 radon samples collected from ground water in the watershed exceeded the proposed USEPA drinking-water maximum contaminant level of 300 picocuries per liter and ranged from 426 to 3,830 picocuries per liter.

Colorado↗

Spatial and temporal changes in microbial community structure associated with recharge-influenced chemical gradients in a contaminated aquifer

In a contaminated water‐table aquifer, we related microbial community structure on aquifer sediments to gradients in 24 geochemical and contaminant variables at five depths, under three recharge conditions. Community amplified ribsosomal DNA restriction analysis (ARDRA) using universal 16S rDNA primers and denaturing gradient gel electrophoresis (DGGE) using bacterial 16S rDNA primers indicated: (i) communities in the anoxic, contaminated central zone were similar regardless of recharge; (ii) after recharge, communities at greatest depth were similar to those in uncontaminated zones; and (iii) after extended lack of recharge, communities at upper and lower aquifer margins differed from communities at the same depths on other dates. General aquifer geochemistry was as important as contaminant or terminal electron accepting process (TEAP) chemistry in discriminant analysis of community groups. The Shannon index of diversity ( H ) and the evenness index ( E ), based on DGGE operational taxonomic units (OTUs), were statistically different across community groups and aquifer depths. Archaea or sulphate‐reducing bacteria 16S rRNA abundance was not clearly correlated with TEAP chemistry indicative of methanogenesis or sulphate reduction. Eukarya rRNA abundance varied by depth and date from 0 to 13% of the microbial community. This contaminated aquifer is a dynamic ecosystem, with complex interactions between physical, chemical and biotic components, which should be considered in the interpretation of aquifer geochemistry and in the development of conceptual or predictive models for natural attenuation or remediation.

Environmental Microbiology↗

Resolving model parameter values from carbon and nitrogen stock measurements in a wide range of tropical mature forests using nonlinear inversion and regression trees

Objectively assessing the performance of a model and deriving model parameter values from observations are critical and challenging in landscape to regional modeling. In this paper, we applied a nonlinear inversion technique to calibrate the ecosystem model CENTURY against carbon (C) and nitrogen (N) stock measurements collected from 39 mature tropical forest sites in seven life zones in Costa Rica. Net primary productivity from the Moderate-Resolution Imaging Spectroradiometer (MODIS), C and N stocks in aboveground live biomass, litter, coarse woody debris (CWD), and in soils were used to calibrate the model. To investigate the resolution of available observations on the number of adjustable parameters, inversion was performed using nine setups of adjustable parameters. Statistics including observation sensitivity, parameter correlation coefficient, parameter sensitivity, and parameter confidence limits were used to evaluate the information content of observations, resolution of model parameters, and overall model performance. Results indicated that soil organic carbon content, soil nitrogen content, and total aboveground biomass carbon had the highest information contents, while measurements of carbon in litter and nitrogen in CWD contributed little to the parameter estimation processes. The available information could resolve the values of 2-4 parameters. Adjusting just one parameter resulted in under-fitting and unacceptable model performance, while adjusting five parameters simultaneously led to over-fitting. Results further indicated that the MODIS NPP values were compressed as compared with the spatial variability of net primary production (NPP) values inferred from inverse modeling. Using inverse modeling to infer NPP and other sensitive model parameters from C and N stock observations provides an opportunity to utilize data collected by national to regional forest inventory systems to reduce the uncertainties in the carbon cycle and generate valuable databases to validate and improve MODIS NPP algorithms.

Ecological Modelling↗

On the use of high-resolution and deep-learning seismic catalogs for short-term earthquake forecasts: Potential benefits and current limitations

Enhanced earthquake catalogs provide detailed images of evolving seismic sequences. Currently, these data sets take some time to be released but will soon become available in real time. Here, we explore whether and how enhanced seismic catalogs feeding into established short-term earthquake forecasting protocols may result in higher predictive skill. We consider three enhanced catalogs for the 2016–2017 Central Italy sequence, featuring a bulk completeness lower by at least two magnitude units compared to the real-time catalog and an improved hypocentral resolution. We use them to inform a set of physical Coulomb Rate-and-State (CRS) and statistical Epidemic-Type Aftershock Sequence (ETAS) models to forecast the space-time occurrence of M3+ events during the first 6 months of the sequence. We track model performance using standard likelihood-based metrics and compare their skill against the best-performing CRS and ETAS models among those developed with the real-time catalog. We find that while the incorporation of the triggering contributions from new small magnitude detections of the enhanced catalogs is beneficial for both types of forecasts, these models do not significantly outperform their respective near real-time benchmarks. To explore the reasons behind this result, we perform targeted sensitivity tests that show how (a) the typical spatial discretizations of forecast experiments ( 2 km) hamper the ability of models to capture highly localized secondary triggering patterns and (b) differences in earthquake parameters (i.e., magnitude and hypocenters) reported in different catalogs can affect forecast evaluation. These findings will contribute toward improving forecast model design and evaluation strategies for next-generation seismic catalogs.

Journal of Geophysical Research--Solid Earth↗

Water-quality characteristics, including sodium-adsorption ratios, for four sites in the Powder River drainage basin, Wyoming and Montana, water years 2001-2004

The U.S. Geological Survey, in cooperation with the Wyoming Department of Environmental Quality, monitors streams throughout the Powder River structural basin in Wyoming and parts of Montana for potential effects of coalbed natural gas development. Specific conductance and sodium-adsorption ratios may be larger in coalbed waters than in stream waters that may receive the discharge waters. Therefore, continuous water-quality instruments for specific conductance were installed and discrete water-quality samples were collected to characterize water quality during water years 2001-2004 at four sites in the Powder River drainage basin: Powder River at Sussex, Wyoming; Crazy Woman Creek near Arvada, Wyoming; Clear Creek near Arvada, Wyoming; and Powder River at Moorhead, Montana. During water years 2001-2004, the median specific conductance of 2,270 microsiemens per centimeter at 25 degrees Celsius (?S/cm) in discrete samples from the Powder River at Sussex, Wyoming, was larger than the median specific conductance of 1,930 ?S/cm in discrete samples collected downstream from the Powder River at Moorhead, Montana. The median specific conductance was smallest in discrete samples from Clear Creek (1,180 ?S/cm), which has a dilution effect on the specific conductance for the Powder River at Moorhead, Montana. The daily mean specific conductance from continuous water-quality instruments during the irrigation season showed the same spatial pattern as specific conductance values for the discrete samples. Dissolved sodium, sodium-adsorption ratios, and dissolved solids generally showed the same spatial pattern as specific conductance. The largest median sodium concentration (274 milligrams per liter) and the largest range of sodium-adsorption ratios (3.7 to 21) were measured in discrete samples from the Powder River at Sussex, Wyoming. Median concentrations of sodium and sodium-adsorption ratios were substantially smaller in Crazy Woman Creek and Clear Creek, which tend to decrease sodium concentrations and sodium-adsorption ratios at the Powder River at Moorhead, Montana. Dissolved-solids concentrations in discrete samples were closely correlated with specific conductance values; Pearson's correlation coefficients were 0.98 or greater for all four sites. Regression equations for discrete values of specific conductance and sodium-adsorption ratios were statistically significant (p-values <0.001) at all four sites. The strongest relation (R2=0.92) was at the Powder River at Sussex, Wyoming. Relations on Crazy Woman Creek (R2=0.91) and Clear Creek (R2=0.83) also were strong. The relation between specific conductance and sodium-adsorption ratios was weakest (R2=0.65) at the Powder River at Moorhead, Montana; however, the relation was still significant. These data indicate that values of specific conductance are useful for estimating sodium-adsorption ratios. A regression model called LOADEST was used to estimate dissolved-solids loads for the four sites. The average daily mean dissolved-solids loads varied among the sites during water year 2004. The largest average daily mean dissolved-solids load was calculated for the Powder River at Moorhead, Montana. Although the smallest concentrations of dissolved solids were in samples from Clear Creek, the smallest average daily mean dissolved-solids load was calculated for Crazy Woman Creek. The largest loads occurred during spring runoff, and the smallest loads occurred in late summer, when streamflows typically were smallest. Dissolved-solids loads may be smaller than average during water years 2001-2004 because of smaller than average streamflow as a result of drought conditions.

Scientific Investigations Report↗

A water-leach procedure for estimating bioaccessibility of elements in soils from transects across the United States and Canada

An objective of the North American Soil Geochemical Landscapes Project is to provide relevant data concerning bioaccessible concentrations of elements in soil to government and other institutions undertaking environmental studies. A protocol was developed that employs a 1-g soil sample agitated overnight with 40 mL of reverse-osmosis de-ionized water for 20 h, and determination of 63 elements following three steps of centrifugation by inductively coupled plasma&ndash;atomic emission spectrometry and inductively coupled plasma&ndash;mass spectrometry the following day. Statistical summaries are presented for those 48 elements (Ag, Al, As, B, Ba, Be, Br, Ca, Cd, Ce, Co, Cr, Cs, Cu, Dy, Er, Eu, Fe, Ga, Gd, Ge, Hf, Ho, I, K, La, Li, Lu, Mg, Mn, Mo, Na, Nb, Nd, Ni, P, Pb, Pr, Rb, Re, S, Sb, Si, Sm, Sn, Sr, Tb, Ti, Tl, Tm, U, V, W, Y, Yb, Zn, Zr, and pH) for which <20% of their data were reported as below the detection limit. The resulting data set contains analyses for 161 A-horizon soils collected along two transects, one along the 38th parallel across the USA and the other from northern Manitoba to the USA&ndash;Mexico border. The spatial distribution of three selected elements (Ca, Cu, and Pb) along the two transects is discussed in this paper both as absolute amounts liberated by the leach and expressed as a percentage of the total, or near-total, amounts determined for the elements. The Ca data reflect broad trends in soil parent materials, their weathering, and subsequent soil development. Calcium concentrations are generally found to be lower in the older soils of the eastern USA. The Cu data are higher in the eastern half of the USA, correlating with soil organic C, with which it is sequestered. The Pb data exhibit little regional variability due to natural sources, but are influenced by anthropogenic sources. Based on the Pb results, the percentage water-extractable data demonstrate promise as a tool for identifying anthropogenic components. The soil&ndash;water partition (distribution) coefficients, K d s (L/kg), were determined and their relevance to estimating bioaccessible amounts of elements to soil fauna and flora is discussed. Finally, a possible link between W concentrations in human urine and water-extractable W levels in Nevada soils is discussed.

Applied Geochemistry↗

Patterns and drivers of cliff erosion in Big Sur, California, USA using repeat photogrammetry, 2017–2023

Seacliff erosion in steep terrain poses major risks to transportation and critical infrastructure. In Big Sur, California, USA, seacliff erosion threatens the sustainability of the central coast stretch of California State Route 1, a transportation corridor that is critical to the region's economy. Published cliff retreat rates for the region range from 1 to 40 cm yr −1 , highlighting that high-resolution, process-based studies could enhance understanding of the causes of spatial and temporal variability. We quantified cliff erosion and investigated its drivers along ∼13 km of the Big Sur coastline at week–month timescales during the late fall to early spring wet seasons between January 2017 and June 2023 by analyzing 3D point clouds developed from aerial imagery using four-dimensional structure-from-motion (4D SfM) photogrammetry techniques. We calculated cliff face retreat rates of 2.23 ± 3.06 cm yr −1 (mean ±1 σ ), an order of magnitude lower than long-term estimated rates for the region (which included large deep-seated landslides), but in line with short-term rates reported across California. Change detection imagery comparison, cliff profiles through time, and statistical analysis reveal a cyclical cliff evolution process in which erosion by wave action at the cliff base destabilizes the cliff and primes it for subsequent failure during precipitation events. Although more erosion by volume could be attributed to precipitation-induced increases in soil moisture (784 m 3 km −1 yr −1 ) compared with erosion attributed to wave power (282 m 3 km −1 yr −1 ), our observations underscore the coupled nature of these processes in driving cliff evolution, consistent with established theory and observations.

California↗

Distribution and trends in reference evapotranspiration in the North China plain

The distribution and trends in reference evapotranspiration (ET(o)) are extremely important to water resources planning for agriculture, and it is widely believed that rates of ET(o) will increase with global warming. This is a big concern in China, where water deficits are common in the North China Plain (NCP). In this study, Penman-Monteith reference evapotranspiration at 26 meteorological stations during 1961-2006 in and around the NCP was calculated. The temporal variations and spatial distribution of ET(o) were analyzed and the causes for the variations were discussed. The results showed that: (1) the NCP was divided into two climatic regions based on aridity values: a semiarid region that accounts for 69% of the area and subhumid regions that made of the remaining area; (2) over the entire NCP, the highest annual ET(o) occurred in the central and western areas and the lowest total ET(o) was observed in the east. Comparing the mean monthly ET(o) and annual ET(o) distributions, the high ET(o) values from May through July mainly determined the annual ET(o) distribution; (3) for the whole NCP, annual ET(o) showed a statistically significant decrease of 11.92 mm/decade over the 46 years of data collection in the NCP or approximately a 5% total decrease compared to the ET(o) values in 1961; (4) to determine which variable has the greatest effect on the decrease in ET(o), decadal changes were observed for daily values of maximum air temperature (+0.16 degrees C), minimum air temperature (+0.35 degrees C), net radiation (-0.13 MJ m(-2)), and mean wind speed (-0.09 m s(-1)). These results indicate that the decreasing net radiation and wind speed had a bigger impact on ET(o) rates than the increases observed by the maximum and minimum temperatures.

North China Plain↗

Estimated probabilities, volumes, and inundation area depths of potential postwildfire debris flows from Carbonate, Slate, Raspberry, and Milton Creeks, near Marble, Gunnison County, Colorado

During 2009, the U.S. Geological Survey, in cooperation with Gunnison County, initiated a study to estimate the potential for postwildfire debris flows to occur in the drainage basins occupied by Carbonate, Slate, Raspberry, and Milton Creeks near Marble, Colorado. Currently (2010), these drainage basins are unburned but could be burned by a future wildfire. Empirical models derived from statistical evaluation of data collected from recently burned basins throughout the intermountain western United States were used to estimate the probability of postwildfire debris-flow occurrence and debris-flow volumes for drainage basins occupied by Carbonate, Slate, Raspberry, and Milton Creeks near Marble. Data for the postwildfire debris-flow models included drainage basin area; area burned and burn severity; percentage of burned area; soil properties; rainfall total and intensity for the 5- and 25-year-recurrence, 1-hour-duration-rainfall; and topographic and soil property characteristics of the drainage basins occupied by the four creeks. A quasi-two-dimensional floodplain computer model (FLO-2D) was used to estimate the spatial distribution and the maximum instantaneous depth of the postwildfire debris-flow material during debris flow on the existing debris-flow fans that issue from the outlets of the four major drainage basins. The postwildfire debris-flow probabilities at the outlet of each drainage basin range from 1 to 19 percent for the 5-year-recurrence, 1-hour-duration rainfall, and from 3 to 35 percent for 25-year-recurrence, 1-hour-duration rainfall. The largest probabilities for postwildfire debris flow are estimated for Raspberry Creek (19 and 35 percent), whereas estimated debris-flow probabilities for the three other creeks range from 1 to 6 percent. The estimated postwildfire debris-flow volumes at the outlet of each creek range from 7,500 to 101,000 cubic meters for the 5-year-recurrence, 1-hour-duration rainfall, and from 9,400 to 126,000 cubic meters for the 25-year-recurrence, 1-hour-duration rainfall. The largest postwildfire debris-flow volumes were estimated for Carbonate Creek and Milton Creek drainage basins, for both the 5- and 25-year-recurrence, 1-hour-duration rainfalls. Results from FLO-2D modeling of the 5-year and 25-year recurrence, 1-hour rainfalls indicate that the debris flows from the four drainage basins would reach or nearly reach the Crystal River. The model estimates maximum instantaneous depths of debris-flow material during postwildfire debris flows that exceeded 5 meters in some areas, but the differences in model results between the 5-year and 25-year recurrence, 1-hour rainfalls are small. Existing stream channels or topographic flow paths likely control the distribution of debris-flow material, and the difference in estimated debris-flow volume (about 25 percent more volume for the 25-year-recurrence, 1-hour-duration rainfall compared to the 5-year-recurrence, 1-hour-duration rainfall) does not seem to substantially affect the estimated spatial distribution of debris-flow material. Historically, the Marble area has experienced periodic debris flows in the absence of wildfire. This report estimates the probability and volume of debris flow and maximum instantaneous inundation area depths after hypothetical wildfire and rainfall. This postwildfire debris-flow report does not address the current (2010) prewildfire debris-flow hazards that exist near Marble.

Colorado↗

Decadal-scale export of nitrogen, phosphorus, and sediment from the Susquehanna River basin, USA: Analysis and synthesis of temporal and spatial patterns

The export of nitrogen (N), phosphorus (P), and suspended sediment (SS) is a long-standing management concern for the Chesapeake Bay watershed, USA. Here we present a comprehensive evaluation of nutrient and sediment loads over the last three decades at multiple locations in the Susquehanna River basin (SRB), Chesapeake's largest tributary watershed. Sediment and nutrient riverine loadings, including both dissolved and particulate fractions, have generally declined at all sites upstream of Conowingo Dam (non-tidal SRB outlet). Period-of-record declines in riverine yield are generally smaller than those in source input, suggesting the possibility of legacy contributions. Consistent with other watershed studies, these results reinforce the importance of considering lag time between the implementation of management actions and achievement of river quality improvement. Whereas flow-normalized loadings for particulate species have increased recently below Conowingo Reservoir, those for upstream sites have declined, thus substantiating conclusions from prior studies about decreased reservoir trapping efficiency. In regard to streamflow effects, statistically significant log-linear relationships between annual streamflow and annual constituent load suggest the dominance of hydrological control on the inter-annual variability of constituent export. Concentration-discharge relationships revealed general chemostasis and mobilization effects for dissolved and particulate species, respectively, both suggesting transport-limitation conditions. In addition to affecting annual export rates, streamflow has also modulated the relative importance of dissolved and particulate fractions, as reflected by its negative correlations with dissolved P/total P, dissolved N/total N, particulate P/SS, and total N/total P ratios. For land-use effects, period-of-record median annual yields of N, P, and SS all correlate positively with the area fraction of non-forested land but negatively with that of forested land under all hydrological conditions. Overall, this work has informed understanding with respect to four major factors affecting constituent export ( i.e. , source input, reservoir modulation, streamflow, and land use) and demonstrated the value of long-term river monitoring.

Maryland, New York, Pennsylvania↗

Depth-duration frequency of precipitation for Texas

The U.S. Geological Survey, in cooperation with the Texas Department of Transportation, conducted a study of the depth-duration frequency of precipitation for Texas. Depth-duration frequency is an estimate of the depth of precipitation for a specified duration and frequency or recurrence interval. For this report, precipitation durations of 15, 30, and 60 minutes; 1, 2, 3, 6, 12, and 24 hours; and 1, 2, 3, 5, and 7 days were investigated. The recurrence intervals for the frequencies range from 2 to 500 years. The time series of precipitation annual maxima for 173 fifteen-minute, 274 hourly, and 865 daily National Weather Service precipitation stations with at least 10 years of record in Texas provide the basis of depth-duration frequency for each identified duration. In total, about 3,030; 10,160; and 38,120 cumulative years of record are available for the 15-minute, hourly, and daily stations, respectively. L-moment statistics of the precipitation annual maxima were calculated for each duration and for each station using unbiased L-moment estimators. The statistics calculated were the mean, L-scale, L-coefficient of variation, L-skew, and L-kurtosis. The mean for each station and duration was corrected for the bias associated with fixed-interval recording of precipitation. The generalized logistic distribution was determined, using L-moment ratio diagrams, as an appropriate probability distribution for modeling the frequency of annual maxima for durations of 15 minutes to 24 hours; whereas, the generalized extreme-value distribution was determined as appropriate for durations of 1 to 7 days. The location, scale, and shape parameters of the distributions for each duration and each station were calculated from the L-moments. These parameters were contoured using spatial interpolation, based on the geostatistical method of kriging, to produce 37 maps that depict the spatial variation and magnitude of each parameter. Contour maps of the shape parameter for the generalized extreme-value distribution for durations of 1 to 7 days are not presented; the root mean square errors of preliminary maps for 1- to 7-day shape parameters were not appreciably smaller than the statewide standard deviation. Therefore, a single statewide mean shape parameter was used for 1- to 7-day durations. The depth-duration frequency for any location in Texas can be estimated using the contour maps and the equation of the corresponding distribution.

Texas↗

Trade-offs of predation and foraging explain sexual segregation in African buffalo

1. Many studies have investigated why males and females segregate spatially in sexually dimorphic species. These studies have focused primarily on temperate zone ungulates in areas lacking intact predator communities, and few have directly assessed predation rates in different social environments. 2. Data on the movement, social affiliation, mortality and foraging of radio-collared African buffalo (Syncerus caffer) were collected from 2001-06 in the Kruger National Park, South Africa. 3. The vast majority of mortality events were due to lion (Panthera leo) predation, and the mortality hazard associated with being an adult male buffalo in a male-only 'bachelor' group was almost four times higher than for adult females in mixed herds. The mortality rates of adult males and females within mixed herds were not statistically different. Mortality sites of male and female buffalo were in areas of low visibility similar to those used by bachelor groups, while mixed herds tended to use more open habitats. 4. Males in bachelor groups ate similar or higher quality food (as indexed by percentage faecal nitrogen), and moved almost a third less distance per day compared with mixed herds. As a result, males in bachelor groups gained more body condition than did males in breeding herds. 5. Recent comparative analyses suggest the activity-budget hypothesis as a common underlying cause of social segregation. However, our intensive study, in an area with an intact predator community showed that male and female buffalo segregated by habitat and supported the predation-risk hypothesis. Male African buffalo appear to trade increased predation risk for additional energy gains in bachelor groups, which presumably leads to increased reproductive success. ?? 2008 The Authors.

Journal of Animal Ecology↗

Evaluation of mercury in rainbow trout collected from Duck Valley Indian Reservation reservoirs, southwestern Idaho and northern Nevada, 2007, 2009, and 2013

The U.S. Geological Survey, in cooperation with the Shoshone-Paiute Tribes of the Duck Valley Indian Reservation, analyzed mercury (Hg) concentration in rainbow trout ( Oncorhynchus mykiss ) collected from three reservoirs on the reservation (Mountain View, Lake Billy Shaw, and Sheep Creek) during sampling events in 2007, 2009, and 2013, to determine the risk of Hg exposure to Tribal members and the general public. Mercury concentration in predatory fish tends to increase with fish length, and this tendency was true for rainbow trout in the reservoirs on the reservation (r 2 = 0.44&ndash;0.70). Mean (average) and median Hg concentrations in fish tissue were determined for each reservoir for each sample year. All Hg concentrations were less than the U.S. Environmental Protection Agency&rsquo;s water-quality criterion of 0.30 milligram per kilogram (mg/kg wet weight [ww]) and the Idaho Department of Environmental Quality&rsquo;s reasonable potential to exceed threshold of 0.24 mg/kg. Idaho Department of Health and Welfare toxicologists determined that the Hg concentrations in rainbow trout in this study would not warrant a fish-consumption advisory for this species. Throughout this report, statistical findings with a p -value of less than 0.05 are referred to as &ldquo;significant.&rdquo; Mean Hg concentrations in fish-tissue samples collected from Mountain View Reservoir were higher in 2007 (0.12 mg/kg ww) than in 2009 and 2013 (0.07 and 0.06 mg/kg ww, respectively), indicating a significant mean decrease. Mean Hg concentrations in fish-tissue samples collected from Lake Billy Shaw showed no significant differences among sample years (2007, 0.12 mg/kg ww; 2009, 0.07 mg/kg ww; 2013, 0.09 mg/kg ww). Mean Hg concentrations in fish-tissue samples collected from Sheep Creek Reservoir significantly increased in 2013 (0.10 mg/kg ww) from concentrations in 2007 and 2009 (0.06 and 0.05 mg/kg ww, respectively). These temporal and spatial variations are not unexpected, as each body of water may differ in the factors and conditions affecting the rate of methylation and demethylation. Coupled with the dynamic put-and-take fishery, the outcomes reflect the system complexities among reservoirs despite their fairly close proximity to one another. The influence of these other factors is evident when the analysis of atmospheric Hg deposition at Mercury Deposition Network site NV02 in northern Nevada showed no significant linear trend in wet Hg deposition rates for 2003&ndash;2013 (average 3.02 micrograms per square meter).

Idaho, Nevada↗

Mapping of compositional properties of coal using isometric log-ratio transformation and sequential Gaussian simulation – A comparative study for spatial ultimate analyses data

Chemical properties of coal largely determine coal handling, processing, beneficiation methods, and design of coal-fired power plants. Furthermore, these properties impact coal strength, coal blending during mining, as well as coal's gas content, which is important for mining safety. In order for these processes and quantitative predictions to be successful, safer, and economically feasible, it is important to determine and map chemical properties of coals accurately in order to infer these properties prior to mining. Ultimate analysis quantifies principal chemical elements in coal. These elements are C, H, N, S, O, and, depending on the basis, ash, and/or moisture. The basis for the data is determined by the condition of the sample at the time of analysis, with an “as-received” basis being the closest to sampling conditions and thus to the in-situ conditions of the coal. The parts determined or calculated as the result of ultimate analyses are compositions, reported in weight percent, and pose the challenges of statistical analyses of compositional data. The treatment of parts using proper compositional methods may be even more important in mapping them, as most mapping methods carry uncertainty due to partial sampling as well. In this work, we map the ultimate analyses parts of the Springfield coal from an Indiana section of the Illinois basin, USA, using sequential Gaussian simulation of isometric log-ratio transformed compositions. We compare the results with those of direct simulations of compositional parts. We also compare the implications of these approaches in calculating other properties using correlations to identify the differences and consequences. Although the study here is for coal, the methods described in the paper are applicable to any situation involving compositional data and its mapping.

Indiana↗

Assessing the impact of site-specific BMPs using a spatially explicit, field-scale SWAT model with edge-of-field and tile hydrology and water-quality data in the Eagle Creek watershed, Ohio

The Eagle Creek watershed, a small subbasin (125 km 2 ) within the Maumee River Basin, Ohio, was selected as a part of the Great Lakes Restoration Initiative (GLRI) “Priority Watersheds” program to evaluate the effectiveness of agricultural Best Management Practices (BMPs) funded through GLRI at the field and watershed scales. The location and quantity of BMPs were obtained from the U.S. Department of Agriculture-Natural Resources Conservation Service National Conservation Planning (NCP) database. A Soil and Water Assessment Tool (SWAT) model was built and calibrated for this predominantly agricultural Eagle Creek watershed, incorporating NCP BMPs and monitoring data at the watershed outlet, an edge-of-field (EOF), and tile monitoring sites. Input air temperature modifications were required to induce simulated tile flow to match monitoring data. Calibration heavily incorporated tile monitoring data to correctly proportion surface and subsurface flow, but calibration statistics were unsatisfactory at the EOF and tile monitoring sites. At the watershed outlet, satisfactory to very good calibration statistics were achieved over a 2-year calibration period, and satisfactory statistics were found in the 2-year validation period. SWAT fixes parameters controlling nutrients primarily at the watershed level; a refinement of these parameters at a smaller-scale could improve field-level calibration. Field-scale modeling results indicate that filter strips (FS) are the most effective single BMPs at reducing dissolved reactive phosphorus, and FS typically decreased sediment and nutrient yields when added to any other BMP or BMP combination. Cover crops were the most effective single, in-field practice by reducing nutrient loads over winter months. Watershed-scale results indicate BMPs can reduce sediment and nutrients, but reductions due to NCP BMPs in the Eagle Creek watershed for all water-quality constituents were less than 10%. Hypothetical scenarios simulated with increased BMP acreages indicate larger investments of the appropriate BMP or BMP combination can decrease watershed level loads.

Ohio↗