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 361 records · Page 20Linked to original sources
USGS2024

Agricultural tile drains increase the susceptibility of streams to longer and more intense streamflow droughts

Streamflow droughts are receiving increased attention worldwide due to their impact on the environment and economy. One region of concern is the Midwestern United States, whose agricultural productivity depends on subsurface pipes known as tile drains to improve trafficability and soil conditions for crop growth. Tile drains accomplish this by rapidly transporting surplus soil moisture and shallow groundwater from fields, resulting in reduced watershed storage. However, no work has previously examined the connection between tile drainage and streamflow drought. Here, we pose the question: does the extent of watershed-level tile drainage lead to an increased susceptibly and magnitude of streamflow droughts? To answer this, we use daily streamflow data for 122 watersheds throughout the Midwestern United States to quantify streamflow drought duration, frequency, and intensity. Using spatial multiple regression models, we find that agricultural tile drainage generates statistically significant ( p < 0.05) increases in streamflow drought duration and intensity while significantly reducing drought frequency. The magnitude of the effect of tile drainage on streamflow drought characteristics is similar to that of water table depth and precipitation seasonality, both of which are known to influence streamflow droughts. Furthermore, projected changes in regional precipitation characteristics will likely drive the installation of additional tile drainage. We find that for each 10% increase in tile-drained watershed area, streamflow drought duration and intensity increase by 0.03 d and 12%, respectively, while frequency decreases by 0.10 events/year. Such increases in tile drainage may lead to more severe streamflow droughts and have a detrimental effect on the socio-environmental usage of streams throughout the Midwest.

Illinois, Indiana, Iowa, Michigan, Minnesota, Ohio↗
USGS2011

Digital hydrologic networks supporting applications related to spatially referenced regression modeling

Digital hydrologic networks depicting surface-water pathways and their associated drainage catchments provide a key component to hydrologic analysis and modeling. Collectively, they form common spatial units that can be used to frame the descriptions of aquatic and watershed processes. In addition, they provide the ability to simulate and route the movement of water and associated constituents throughout the landscape. Digital hydrologic networks have evolved from derivatives of mapping products to detailed, interconnected, spatially referenced networks of water pathways, drainage areas, and stream and watershed characteristics. These properties are important because they enhance the ability to spatially evaluate factors that affect the sources and transport of water-quality constituents at various scales. SPAtially Referenced Regressions On Watershed attributes (SPARROW), a process-based &frasl; statistical model, relies on a digital hydrologic network in order to establish relations between quantities of monitored contaminant flux, contaminant sources, and the associated physical characteristics affecting contaminant transport. Digital hydrologic networks modified from the River Reach File (RF1) and National Hydrography Dataset (NHD) geospatial datasets provided frameworks for SPARROW in six regions of the conterminous United States. In addition, characteristics of the modified RF1 were used to update estimates of mean-annual streamflow. This produced more current flow estimates for use in SPARROW modeling.

Journal of the American Water Resources Associatio↗
USGS2008

Sediment Quality and Comparison to Historical Water Quality, Little Arkansas River Basin, South-Central Kansas, 2007

The spatial and temporal variability in streambed-sediment quality and its relation to historical water quality was assessed to provide guidance for the development of total maximum daily loads and the implementation of best-management practices in the Little Arkansas River Basin, south-central Kansas. Streambed-sediment samples were collected at 26 sites in 2007, sieved to isolate the less than 63-micron fraction (that is, the silt and clay), and analyzed for selected nutrients (total nitrogen and total phosphorus), organic and total carbon, 25 trace elements, and the radionuclides beryllium-7, cesium-137, lead-210, and radium-226. At eight sites, streambed-sediment samples also were collected and analyzed for bacteria. Particulate nitrogen, phosphorus, and organic carbon concentrations in the streambed sediment varied substantially spatially and temporally, and positive correlations among the three constituents were statistically significant. Along the main-stem Little Arkansas River, streambed-sediment concentrations of particulate nitrogen and phosphorus generally were larger at and downstream from Alta Mills, Kansas. The largest particulate nitrogen concentrations were measured in samples collected in the Emma Creek subbasin and may be related to livestock and poultry production. The largest particulate phosphorus concentrations in the basin were measured in samples collected along the main-stem Little Arkansas River downstream from Alta Mills, Kansas. Particulate nitrogen, phosphorus, and organic carbon content in the water and streambed-sediment samples typically decreased as streamflow increased. This inverse relation may be caused by an increased contribution of sediment from channel-bank sources during high flows and (or) increased particle sizes transported by the high flows. Trace element concentrations in the streambed sediment varied from site to site and typically were less than threshold-effects guidelines for possible adverse biological effects. The largest copper, lead, silver, and zinc concentrations, measured for a sample collected from Sand Creek downstream from Newton, Kansas, likely were related to urban sources of contamination. Radionuclide activities and bacterial densities in the streambed sediment varied throughout the basin. Variability in the former may be indicative of subbasin differences in the contribution of sediment from surface-soil and channel-bank sources. Streambed sediment may be useful for reconnaissance purposes to determine sources of particulate nitrogen, phosphorus, organic carbon, and other sediment-associated constituents in the basin. If flow conditions prior to streambed-sediment sampling and during water-quality sampling are considered, it may be possible to use streambed sediment as an indicator of water quality for nitrogen, phosphorus, and organic carbon. Flow conditions affect sediment-associated constituent concentrations in streambed-sediment and water samples, in part, because the sources of sediment (surface soils, channel banks) can vary with flow as can the size of the particles transported.

Scientific Investigations Report↗
USGS2011

Formulation of a correlated variables methodology for assessment of continuous gas resources with an application to the Woodford play, Arkoma Basin, eastern Oklahoma

Shale gas is a form of continuous unconventional hydrocarbon accumulation whose resource estimation is unfeasible through the inference of pore volume. Under these circumstances, the usual approach is to base the assessment on well productivity through estimated ultimate recovery (EUR). Unconventional resource assessments that consider uncertainty are typically done by applying analytical procedures based on classical statistics theory that ignores geographical location, does not take into account spatial correlation, and assumes independence of EUR from other variables that may enter into the modeling. We formulate a new, more comprehensive approach based on sequential simulation to test methodologies known to be capable of more fully utilizing the data and overcoming unrealistic simplifications. Theoretical requirements demand modeling of EUR as areal density instead of well EUR. The new experimental methodology is illustrated by evaluating a gas play in the Woodford Shale in the Arkoma Basin of Oklahoma. Differently from previous assessments, we used net thickness and vitrinite reflectance as secondary variables correlated to cell EUR. In addition to the traditional probability distribution for undiscovered resources, the new methodology provides maps of EUR density and maps with probabilities to reach any given cell EUR, which are useful to visualize geographical variations in prospectivity.

Oklahoma↗
USGS2018

Interannual snow accumulation variability on glaciers derived from repeat spatially extensive ground-penetrating radar surveys

There is significant uncertainty regarding the spatiotemporal distribution of seasonal snow on glaciers, despite being a fundamental component of glacier mass balance. To address this knowledge gap, we collected repeat, spatially extensive high-frequency ground-penetrating radar (GPR) observations on two glaciers in Alaska for five consecutive years. GPR measurements showed steep snow water equivalent (SWE) elevation gradients at both sites; continental Gulkana Glacier’s SWE gradient averaged 115 mm 100 m–1 and maritime Wolverine Glacier’s gradient averaged 440 mm 100 m–1 (over >1000 m). We extrapolated GPR point observations across the glacier surface using terrain parameters derived from digital elevation models as predictor variables in two statistical models (stepwise multivariable linear regression and regression trees). Elevation and proxies for wind redistribution had the greatest explanatory power, and exhibited relatively time-constant coefficients over the study period. Both statistical models yielded comparable estimates of glacier-wide average SWE (1 % average difference at Gulkana, 4 % average difference at Wolverine), although the spatial distributions produced by the models diverged in unsampled regions of the glacier, particularly at Wolverine. In total, six different methods for estimating the glacier-wide average agreed within ± 11 %. We assessed interannual variability in the spatial pattern of snow accumulation predicted by the statistical models using two quantitative metrics. Both glaciers exhibited a high degree of temporal stability, with ~85 % of the glacier area experiencing less than 25 % normalized absolute variability over this five-year interval. We found SWE at a sparse network (3 stakes per glacier) of long-term glaciological stake sites to be highly correlated with the GPR-derived glacier-wide average. We estimate that interannual variability in the spatial pattern of SWE is only a small component (4–10 % of glacier-wide average) of the total mass balance uncertainty and thus, our findings support the concept that sparse stake networks effectively measure interannual variability in winter balance on glaciers, rather than some spatially varying pattern of snow accumulation.

Alaska↗
USGS2021

A machine learning approach to modeling streamflow with sparse data in ungaged watersheds on the Wyoming Range, Wyoming, 2012–17

Scant availability of streamflow data can impede the utility of streamflow as a variable in ecological models of aquatic and terrestrial species, especially when studying small streams in watersheds that lack streamgages. Streamflow data at fine resolution and broad extent were needed by collaborators for ecological research on small streams in several ungaged watersheds of southwestern Wyoming, where streamflow data are sparse. To improve the utility of sparse streamflow data to ecological research in ungaged watersheds, we developed a machine learning approach in R for modeling spatially and temporally continuous monthly streamflow from 2012 through 2017 in three semiarid montane-steppe watersheds (with drainage areas of 26–55 square miles and mean elevations of 8,031–8,455 feet) on the Wyoming Range in the upper Green River Basin. A machine learning streamflow (MLFLOW) model was calibrated and validated with 971 discrete streamflow observations and 24 static and dynamic predictor variables derived from geospatial and time series data on climatic, physiographic, and anthropogenic characteristics affecting streamflow. The predictor variables were temporally and spatially conditioned to amplify the relation of predictor variables to monthly streamflow. The MLFLOW model had satisfactory agreement between observed and predicted streamflow (coefficient of determination [ R 2 ]=0.80, Nash-Sutcliffe efficiency [NSE]=0.79, NSE with log-transformed data [logNSE]=0.82, and percent bias [PBIAS]=0.7 percent). NSE and logNSE indicated the MLFLOW model performed equally well for high and low flows, and PBIAS indicated the MLFLOW model did not overpredict or underpredict monthly streamflow. Streamflow predictions seemed to well represent the annual hydrograph within the study area during the study period. The most important variables (statistically important in the MLFLOW model) for explaining monthly streamflow were temporally and spatially conditioned dynamic climatic variables, mostly precipitation and snow water equivalent. Importance of the static and dynamic variables did not differ substantially among the three watersheds but differed considerably among the 6 years. Monthly streamflow increased with increasing precipitation, snow water equivalent, and drainage area but decreased with increasing forest cover, elevation, evapotranspiration, and temperature. The MLFLOW model was most sensitive to selection of dynamic climatic variables. Unconditioned dynamic climatic variables alone explained 54 percent of the variance ( R 2 =0.54) in monthly streamflow, whereas adding static physiographic and anthropogenic variables only explained 12 percent more of the variance ( R 2 =0.66). Also, spatial conditioning of all variables together with temporal conditioning of dynamic variables increased the variance explained in the MLFLOW model by another 14 percent ( R 2 =0.80). The MLFLOW model also had greater sensitivity to temporal than to spatial differences in the data. For the MLFLOW model trained with observations from all watersheds and years or for models trained with observations from all except one watershed or 1 year left out sequentially, performance was better in testing on observations from each watershed than from each year separately. Also, performance was better for models fitted to fewer sites than to fewer months of observations. The greatest utility of the modeling approach is the ease of use and the speed of processing input data, running the model, and interpreting the model output, whereas the greatest limitation is the need for spatially and temporally representative streamflow observations to drive the model. Although familiarity with R is necessary, only a working knowledge of hydrology (for selecting appropriate predictor variables and evaluating the quality of streamflow observations) and a rudimentary understanding of machine learning models are needed. Therefore, this modeling approach is practicable for other scientists who work with water but who are not hydrologists.

Wyoming↗
USGS2003

Strong ground-motion prediction from Stochastic-dynamic source models

In the absence of sufficient data in the very near source, predictions of the intensity and variability of ground motions from future large earthquakes depend strongly on our ability to develop realistic models of the earthquake source. In this article we simulate near-fault strong ground motion using dynamic source models. We use a boundary integral method to simulate dynamic rupture of earthquakes by specifying dynamic source parameters (fracture energy and stress drop) as spatial random fields. We choose these quantities such that they are consistent with the statistical properties of slip heterogeneity found in finite-source models of past earthquakes. From these rupture models we compute theoretical strong-motion seismograms up to a frequency of 2 Hz for several realizations of a scenario strike-slip Mw 7.0 earthquake and compare empirical response spectra, spectra obtained from our dynamic models, and spectra determined from corresponding kinematic simulations. We find that spatial and temporal variations in slip, slip rise time, and rupture propagation consistent with dynamic rupture models exert a strong influence on near-source ground motion. Our results lead to a feasible approach to specify the variability in the rupture time distribution in kinematic models through a generalization of Andrews' (1976) result relating rupture speed to apparent fracture energy, stress drop, and crack length to 3D dynamic models. This suggests that a simplified representation of dynamic rupture may be obtained to approximate the effects of dynamic rupture without having to do full dynamic simulations.

Bulletin of the Seismological Society of America↗
USGS2024

First records distribution models to guide biosurveillance for non-native species

Quickly locating new populations of non-native species can reduce the ecological and economic costs of species invasions. However, the difficulty of predicting which new non-native species will establish, and where, has limited active post-border biosurveillance efforts. Because pathways of introduction underlie spatial patterns of establishment risk, an intuitive approach is to search for new non-native species in areas where many non-native species have first been detected in the past. We formalize this intuition via first records distribution models (FRDMs), which apply species distribution modeling methods to the collection of first occurrence records across species (i.e. one record per species). We define FRDMs as statistical models that quantify environmental conditions associated with species' first naturalized records to predict spatial patterns of establishment risk. We model the first records of non-native plants in the conterminous USA as a proof-of-concept. The novelty of FRDMs is that their inferences apply not just to the species that contributed data; they provide a rigorous framework for predicting hotspots of invasion for new non-native taxa that share a pathway of introduction with the modeled species. FRDMs can guide survey efforts for new non-native taxa at multiple scales and across ecosystems.

conterminous United States↗
USGS2012

Probabilistic assessment of beach and dune changes

The recent availability of spatially-dense airborne lidar data makes assessment of the vulnerability of beaches and dunes to storm impacts practical over long reaches of coast. As an initial test, elevations of the tops (D high) and bases (Dlow) of foredune ridges along a 55-km reach on the northern Outer Banks, NC were found to have considerable spatial variability suggesting that different parts of the barrier island would respond differently to storms. Comparing statistics of storm wave runup to D high and Dlow, we found that net erosion due to overwash and dune retreat should be greatest at the northern and southern ends of the study area and least in the central section. This predicted spatial pattern of storm-induced erosion is similar to the spatial pattern of long-term erosion of the shoreline which may be controlled by additional processes (such as gradients in longshore transport) as well as the cross-shore processes considered here. However, consider feedback where at erosional hot spots there is a deficit of sand (caused by gradients in longshore transport) which lead to lower dunes and enhanced erosional cross-shore processes, such as overwash. Hence, the erosional hot spots would be exacerbated, further increasing the vulnerability of the beach and dunes to net erosion.

Conference Paper↗
USGS2024

A benchmark dataset and workflow for landslide susceptibility zonation

Landslide susceptibility shows the spatial likelihood of landslide occurrence in a specific geographical area and is a relevant tool for mitigating the impact of landslides worldwide. As such, it is the subject of countless scientific studies. Many methods exist for generating a susceptibility map, mostly falling under the definition of statistical or machine learning. These models try to solve a classification problem: given a collection of spatial variables, and their combination associated with landslide presence or absence, a model should be trained, tested to reproduce the target outcome, and eventually applied to unseen data. Contrary to many fields of science that use machine learning for specific tasks, no reference data exist to assess the performance of a given method for landslide susceptibility. Here, we propose a benchmark dataset consisting of 7360 slope units encompassing an area of about 4,100 km 2 "> 4,100 km 2 in Central Italy. Using the dataset, we tried to answer two open questions in landslide research: (1) what effect does the human variability have in creating susceptibility models; (2) how can we develop a reproducible workflow for allowing meaningful model comparisons within the landslide susceptibility research community. With these questions in mind, we released a preliminary version of the dataset, along with a “call for collaboration,” aimed at collecting different calculations using the proposed data, and leaving the freedom of implementation to the respondents. Contributions were different in many respects, including classification methods, use of predictors, implementation of training/validation, and performance assessment. That feedback suggested refining the initial dataset, and constraining the implementation workflow. This resulted in a final benchmark dataset and landslide susceptibility maps obtained with many classification methods. Values of area under the receiver operating characteristic curve obtained with the final benchmark dataset were rather similar, as an effect of constraints on training, cross–validation, and use of data. Brier score results show larger variability, instead, ascribed to different model predictive abilities. Correlation plots show similarities between results of different methods applied by the same group, ascribed to a residual implementation dependence. We stress that the experiment did not intend to select the “best” method but only to establish a first benchmark dataset and workflow, that may be useful as a standard reference for calculations by other scholars. The experiment, to our knowledge, is the first of its kind for landslide susceptibility modeling. The data and workflow presented here comparatively assess the performance of independent methods for landslide susceptibility and we suggest the benchmark approach as a best practice for quantitative research in geosciences.

Earth-Science Reviews↗
USGS1997

Spatial analysis of temperature (BHT/DST) data and consequences for heat-flow determination in sedimentary basins

Large numbers of bottom-hole temperatures (BHTs) and temperatures measured during drill-stem tests (DSTs) are available in areas explored for hydrocarbons, but their usefulness for estimating geothermal gradients and heat-flow density is limited. We investigated a large data set of BHT and DST measurements taken in boreholes in the American Midcontinent, a geologically uniform stable cratonic area, and propose an empirical correction for BHTs based on relationships between BHTs, DSTs, and thermal logs. This empirical correction is compared with similar approaches determined for other areas. The data were analyzed by multivariate statistics prior to the BHT correction to identify anomalous measurements and quantify external influences. Spatial patterns in temperature measurements for major stratigraphic units outline relations to regional structure. Comparision of temperature and structure trend-surface residuals reveals a relationship between temperature highs and local structure highs. The anticlines, developed by continuous but intermittent movement of basement fault blocks in the Late Paleozoic, are subtle features having closures of 10-30 m and contain relatively small hydrocarbon reservoirs. The temperature anomalies of the order of 5-7 ??C may reflect fluids moving upward along fractures and faults, rather than changes in thermal conductivity resulting from different pore fluids. ?? Springer-Verlag 1997.

Conference Paper↗
USGS2011

Effects of uncertain topographic input data on two-dimensional flow modeling in a gravel-bed river

Many applications in river research and management rely upon two-dimensional (2D) numerical models to characterize flow fields, assess habitat conditions, and evaluate channel stability. Predictions from such models are potentially highly uncertain due to the uncertainty associated with the topographic data provided as input. This study used a spatial stochastic simulation strategy to examine the effects of topographic uncertainty on flow modeling. Many, equally likely bed elevation realizations for a simple meander bend were generated and propagated through a typical 2D model to produce distributions of water-surface elevation, depth, velocity, and boundary shear stress at each node of the model's computational grid. Ensemble summary statistics were used to characterize the uncertainty associated with these predictions and to examine the spatial structure of this uncertainty in relation to channel morphology. Simulations conditioned to different data configurations indicated that model predictions became increasingly uncertain as the spacing between surveyed cross sections increased. Model sensitivity to topographic uncertainty was greater for base flow conditions than for a higher, subbankfull flow (75% of bankfull discharge). The degree of sensitivity also varied spatially throughout the bend, with the greatest uncertainty occurring over the point bar where the flow field was influenced by topographic steering effects. Uncertain topography can therefore introduce significant uncertainty to analyses of habitat suitability and bed mobility based on flow model output. In the presence of such uncertainty, the results of these studies are most appropriately represented in probabilistic terms using distributions of model predictions derived from a series of topographic realizations. Copyright 2011 by the American Geophysical Union.

Water Resources Research↗
USGS2011

Developing effective sampling designs for monitoring natural resources in Alaskan national parks: an example using simulations and vegetation data

Monitoring natural resources in Alaskan national parks is challenging because of their remoteness, limited accessibility, and high sampling costs. We describe an iterative, three-phased process for developing sampling designs based on our efforts to establish a vegetation monitoring program in southwest Alaska. In the first phase, we defined a sampling frame based on land ownership and specific vegetated habitats within the park boundaries and used Path Distance analysis tools to create a GIS layer that delineated portions of each park that could be feasibly accessed for ground sampling. In the second phase, we used simulations based on landcover maps to identify size and configuration of the ground sampling units (single plots or grids of plots) and to refine areas to be potentially sampled. In the third phase, we used a second set of simulations to estimate sample size and sampling frequency required to have a reasonable chance of detecting a minimum trend in vegetation cover for a specified time period and level of statistical confidence. Results of the first set of simulations indicated that a spatially balanced random sample of single plots from the most common landcover types yielded the most efficient sampling scheme. Results of the second set of simulations were compared with field data and indicated that we should be able to detect at least a 25% change in vegetation attributes over 31. years by sampling 8 or more plots per year every five years in focal landcover types. This approach would be especially useful in situations where ground sampling is restricted by access.

Alaska↗
USGS2022

Geomorphic controls on floodplain connectivity, ecosystem services, and sensitivity to climate change: An example from the lower Missouri River

Floodplains of large rivers are exploited for agricultural production, industrial and municipal development, and transportation infrastructure. Recently, increased frequency of costly floods has prompted consideration of whether offsetting benefits might accrue from management of floodplains for ecosystem services. We employed a simple inundation model for 800 km of the Lower Missouri River, USA, to evaluate spatial and temporal distributions of ecological floodplain inundation metrics and how those distributions might vary with levee removal and climatic change. The model evaluates inundation at 30 × 30 m resolution on a daily basis over 82 years of record. We quantified provisioning of waterfowl habitat and potential denitrification. Spatial variability is affected by ongoing geomorphic adjustments that affect floodplain connectivity. Statistical models indicate that available floodplain area and recent aggradation are predictive of most inundation metrics. Connectivity is sensitive to climate-change scenarios that predict increased floodplain inundation during spring waterfowl migrations; the greatest sensitivity to future climate exists where channel-floodplain geomorphology presently enhances floodplain connectivity. Evaluation of floodplain denitrification indicates that on average, the nonleveed part of the floodplain could denitrify 0.05%–1.7% of the mean annual nitrogen load of the river. Levee removal could increase this rate to only 3.6% of the nitrogen load. The capacity of floodplain connectivity to influence certain ecosystem services is highly variable in space along the Lower Missouri River and may be appreciably influenced by climate change. Hence, decisions to optimize management of large-river floodplains are likely to be highly location dependent.

Illinois, Iowa, Kansas, Missouri, Nebraska, South ↗
USGS2009

Using geometrical, textural, and contextual information of land parcels for classification of detailed urban land use

Detailed urban land use data are important to government officials, researchers, and businesspeople for a variety of purposes. This article presents an approach to classifying detailed urban land use based on geometrical, textural, and contextual information of land parcels. An area of 6 by 14 km in Austin, Texas, with land parcel boundaries delineated by the Travis Central Appraisal District of Travis County, Texas, is tested for the approach. We derive fifty parcel attributes from relevant geographic information system (GIS) and remote sensing data and use them to discriminate among nine urban land uses: single family, multifamily, commercial, office, industrial, civic, open space, transportation, and undeveloped. Half of the 33,025 parcels in the study area are used as training data for land use classification and the other half are used as testing data for accuracy assessment. The best result with a decision tree classification algorithm has an overall accuracy of 96 percent and a kappa coefficient of 0.78, and two naive, baseline models based on the majority rule and the spatial autocorrelation rule have overall accuracy of 89 percent and 79 percent, respectively. The algorithm is relatively good at classifying single-family, multifamily, commercial, open space, and undeveloped land uses and relatively poor at classifying office, industrial, civic, and transportation land uses. The most important attributes for land use classification are the geometrical attributes, particularly those related to building areas. Next are the contextual attributes, particularly those relevant to the spatial relationship between buildings, then the textural attributes, particularly the semivariance texture statistic from 0.61-m resolution images.

Annals of the Association of American Geographers↗
USGS2012

Foraging optimally for home ranges

Economic models predict behavior of animals based on the presumption that natural selection has shaped behaviors important to an animal's fitness to maximize benefits over costs. Economic analyses have shown that territories of animals are structured by trade-offs between benefits gained from resources and costs of defending them. Intuitively, home ranges should be similarly structured, but trade-offs are difficult to assess because there are no costs of defense, thus economic models of home-range behavior are rare. We present economic models that predict how home ranges can be efficient with respect to spatially distributed resources, discounted for travel costs, under 2 strategies of optimization, resource maximization and area minimization. We show how constraints such as competitors can influence structure of homes ranges through resource depression, ultimately structuring density of animals within a population and their distribution on a landscape. We present simulations based on these models to show how they can be generally predictive of home-range behavior and the mechanisms that structure the spatial distribution of animals. We also show how contiguous home ranges estimated statistically from location data can be misleading for animals that optimize home ranges on landscapes with patchily distributed resources. We conclude with a summary of how we applied our models to nonterritorial black bears (Ursus americanus) living in the mountains of North Carolina, where we found their home ranges were best predicted by an area-minimization strategy constrained by intraspecific competition within a social hierarchy. Economic models can provide strong inference about home-range behavior and the resources that structure home ranges by offering falsifiable, a priori hypotheses that can be tested with field observations.

Journal of Mammalogy↗
USGS2020

Informing amphibian conservation efforts with abundance-based metapopulation models

Science-based management strategies are needed to halt or reverse the global decline of amphibians. In many cases, sound management requires reliable models built using monitoring data. Historically, monitoring and statistical modeling efforts have focused on estimating occupancy using detection–nondetection data. Spatial occupancy models are useful for studying colonization–extinction dynamics, but richer insights can be gained from estimating abundance and density-dependent demographic rates. We developed an integrated abundance-based metapopulation model of the processes contributing to spatiotemporal variation in patch population density. We fit our model to a combination of detection–nondetection and count data from a 14-yr study of a reintroduced metapopulation of federally threatened Chiricahua Leopard Frogs (Lithobates chiricahuensis). Pond-specific population growth rate was influenced by pond hydroperiod and frog density, such that permanent and semipermanent ponds with low densities of adult frogs experienced the highest annual population growth rates. Immigration rate declined as the distance among ponds increased. After reintroduction in 2003, metapopulation-level abundance increased and appeared to stabilize around 1300 adult frogs (95% CI = 1192–1471) by year 2015. Further, changes in metapopulation abundance were driven mostly by changes in abundance at a few ponds. These high-density populations, which would not have been identifiable with traditional occupancy-based metapopulation models, are likely especially important for species recovery in the area. Abundance-based metapopulation models can be widely applied to inform conservation efforts, by providing higher quality information needed to prioritize habitat patches for management and can be used to make more accurate predictions of metapopulation extinction risk.

Herpetologica↗
USGS2014

HydroClimATe: hydrologic and climatic analysis toolkit

The potential consequences of climate variability and climate change have been identified as major issues for the sustainability and availability of the worldwide water resources. Unlike global climate change, climate variability represents deviations from the long-term state of the climate over periods of a few years to several decades. Currently, rich hydrologic time-series data are available, but the combination of data preparation and statistical methods developed by the U.S. Geological Survey as part of the Groundwater Resources Program is relatively unavailable to hydrologists and engineers who could benefit from estimates of climate variability and its effects on periodic recharge and water-resource availability. This report documents HydroClimATe, a computer program for assessing the relations between variable climatic and hydrologic time-series data. HydroClimATe was developed for a Windows operating system. The software includes statistical tools for (1) time-series preprocessing, (2) spectral analysis, (3) spatial and temporal analysis, (4) correlation analysis, and (5) projections. The time-series preprocessing tools include spline fitting, standardization using a normal or gamma distribution, and transformation by a cumulative departure. The spectral analysis tools include discrete Fourier transform, maximum entropy method, and singular spectrum analysis. The spatial and temporal analysis tool is empirical orthogonal function analysis. The correlation analysis tools are linear regression and lag correlation. The projection tools include autoregressive time-series modeling and generation of many realizations. These tools are demonstrated in four examples that use stream-flow discharge data, groundwater-level records, gridded time series of precipitation data, and the Multivariate ENSO Index.

Techniques and Methods↗