Geology Reports⌕ Search

SEARCH · Geology Reports

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 865 records · Page 48Linked to original sources

Desert tortoises ( Gopherus agassizii ) are selective herbivores that track the flowering phenology of their preferred food plants

Previous studies of desert tortoise foraging ecology in the western Mojave Desert suggest that these animals are selective herbivores, which alter their diet according to the temporal availability of preferred food plants. These studies, however, did not estimate availability of potential food plants by taking into account the spatial and temporal variability in ephemeral plant abundance that occurs within the spring season. In this study, we observed 18 free-ranging adult tortoises take 35,388 bites during the spring foraging season. We also estimated the relative abundance of potential food plants by stratifying our sampling across different phenological periods of the 3-month long spring season and by different habitats and microhabitats. This methodology allowed us to conduct statistical tests comparing tortoise diet against plant abundance. Our results show that tortoises choose food plants non-randomly throughout the foraging season, a finding that corroborates the hypothesis that desert tortoises rely on key plants during different phenological periods of spring. Moreover, tortoises only consumed plants in a succulent state until the last few weeks of spring, at which time most annuals and herbaceous perennials had dried and most tortoises had ceased foraging. Many species of food plants—including several frequently eaten species—were not detected in our plant surveys, yet tortoises located these rare plants in their home ranges. Over 50% of bites consumed were in the group of undetected species. Interestingly, tortoises focused heavily on several leguminous species, which could be nutritious foods owing to their presumably high nitrogen contents. We suggest that herbaceous perennials, which were rare on our study area but represented ~30% of tortoise diet, may be important in sustaining tortoise populations during droughts when native annuals are absent. These findings highlight the vulnerability of desert tortoises to climate change if such changes alter the availability of their preferred food plants.

California↗

Combining models of the critical streakline and the cross-sectional distribution of juvenile salmon to predict fish routing at river junctions

Because fish that enter the interior Delta have poorer survival than those emigrating via the Sacramento River, understanding the mechanisms that drive entrainment rates at side channel junctions is critically important for the management of imperiled juvenile salmon. Here, we implement a previously proposed process-based conceptual model to study entrainment rates based on three linked elements: the entrainment zone, critical streakline, and cross-sectional distribution of fish. The critical streakline is the location along a channel cross-section immediately upstream of a junction that forms the spatial divide between parcels of water that enter a side channel or remain in the main channel. The critical streakline therefore divides the main channel into entrainment zones within which fish would likely enter each channel. Combined with information about the cross-sectional distribution of fish upstream of a junction, this conceptual model provides a means to predict fish entrainment into each channel. To apply this conceptual model, we combined statistical models of the critical streakline, the cross-sectional distribution of acoustic tagged juvenile Chinook salmon, and their probability of entrainment into Georgiana Slough. We fit joint beta regression and logistic regression models to acoustic telemetry data gathered in 2011 and 2012 to estimate the cross-sectional distribution of fish upstream of the junction, and to estimate the probability of entrainment for fish on either side of the critical streakline. We show that entrainment rates can be predicted by understanding how the combination of critical streakline position and cross-sectional distribution of fish co-vary as a function of environmental covariates. By integrating over individual positions and entrainment fates to arrive at population-level entrain probability in relation to environmental covariates, our model offers managers a simple but powerful tool to evaluate how alternative actions affect migrating fish.

California↗

Magnitude and Frequency of Floods in New York

Techniques are presented for estimating the magnitude and frequency of flood discharges on rural, unregulated streams in New York, excluding Long Island. Peak-discharge-frequency data and basin characteristics from 388 streamflow-gaging stations in New York and adjacent states were used to develop multiple linear regression equations for flood discharges with recurrence intervals ranging from 1.25 to 500 years. A generalized least-squares (GLS) procedure was used to develop the regression equations. Separate sets of equations were developed for each of six hydrologic regions of New York; standard errors of prediction range from 14 to 43 percent. Statistically significant explanatory variables in the regression equations include drainage area, main-channel slope, percent basin storage, mean annual precipitation, percent forested area, a basin lag factor, a ratio of main-channel slope to basin slope, mean annual runoff, maximum snow depth, and percentage of basin above 1,200 feet. Drainage areas for the 388 sites used in the analyses ranged from 0.41 to 4,773 square miles. Methods of computing flood discharges from the regression equations differ, depending on whether the estimate is for a gaged or ungaged basin, and whether the basin crosses hydrologic-region or state boundaries. Examples of computations are included. Discharge-frequency estimates for an additional 122 streamflow-gaging stations with significant regulation or urbanization (including Long Island) are also included as at-site estimates. Basin characteristics, log-Pearson Type III statistics, and regression and weighted estimates of the discharge-frequency relations are tabulated for the streamflow-gaging stations used in the regression analyses. Sensitivity analyses showed that mean-annual precipitation, drainage area, mean annual runoff, and maximum snow depth are the variables to which computed discharges are most sensitive in the regression equations. Included with the report is a DVD that provides computation procedures and geographic information system spatial datasets to compute basin characteristics used in the regional regression equations and flood-frequency estimates at a specified location on a stream.

Scientific Investigations Report↗

Potential use of weather radar to study movements of wintering waterfowl

To protect and restore wintering waterfowl habitat, managers require knowledge of routine wintering waterfowl movements and habitat use. During preliminary screening of Doppler weather radar data we observed biological movements consistent with routine foraging flights of wintering waterfowl known to occur near Lacassine National Wildlife Refuge (NWR), Louisiana. During the winters of 2004&ndash;2005 and 2005&ndash;2006, we conducted field surveys to identify the source of the radar echoes emanating from Lacassine NWR. We compared field data to weather radar reflectivity data. Spatial and temporal patterns consistent with foraging flight movements appeared in weather radar data on all dates of field surveys. Dabbling ducks were the dominant taxa flying within the radar beam during the foraging flight period. Using linear regression, we found a positive log-linear relationship between average radar reflectivity ( Z ) and number of birds detected over the study area ( P < 0.001, r 2 = 0.62, n = 40). Ground observations and the statistically significant relationship between radar data and field data confirm that Doppler weather radar recorded the foraging flights of dabbling ducks. Weather radars may be effective tools for wintering waterfowl management because they provide broad-scale views of both diurnal and nocturnal movements. In addition, an extensive data archive enables the study of wintering waterfowl response to habitat loss, agricultural practices, wetland restoration, and other research questions that require multiple years of data.

Louisiana↗

Historical patterns of acidification and increasing CO2 flux associated with Florida springs

Florida has one of the highest concentrations of springs in the world, with many discharging into rivers and predominantly into eastern Gulf of Mexico coast, and they likely influence the hydrochemistry of these adjacent waters; however, temporal and spatial trends have not been well studied. We present over 20 yr of hydrochemical, seasonally sampled data to identify temporal and spatial trends of pH, alkalinity, partial pressure of carbon dioxide (pCO 2 ), and CO 2 flux from five first-order-magnitude (springs that discharge greater than 2.83 m 3 s −1 ) coastal spring groups fed by the Floridan Aquifer System that ultimately discharge into the Gulf of Mexico. All spring groups had pCO 2 levels (averages 3174.3–6773.2 μatm) that were much higher than atmospheric levels of CO 2 and demonstrated statistically significant temporal decreases in pH and increases in CO 2 flux, pCO 2 , and alkalinity. Total carbon flux emissions increased from each of the spring groups by between 3.48 × 10 7 and 2.856 × 10 8 kg C yr −1 over the time period. By 2013 the Springs Groups in total emitted more than 1.1739 × 10 9 kg C yr −1 . Increases in alkalinity and pCO 2 varied from 90.9 to 347.6 μmol kg −1 and 1262.3 to 2666.7 μatm, respectively. Coastal data show higher CO 2 evasion than the open Gulf of Mexico, which suggests spring water influences nearshore waters. The results of this study have important implications for spring water quality, dissolution of the Florida carbonate platform, and identification of the effect and partitioning of carbon fluxes to and within coastal and marine ecosystems.

Florida↗

Assessing specific-capacity data and short-term aquifer testing to estimate hydraulic properties in alluvial aquifers of the Rocky Mountains, Colorado, USA

Study Region : Rocky Mountains, United States Study Focus : Groundwater-flow modeling requires estimates of hydraulic properties, namely hydraulic conductivity. Hydraulic conductivity values commonly vary over orders of magnitudes however and estimation may require extensive field campaigns applying slug or pumping tests. As an alternative, specific-capacity tests can be used to estimate hydraulic properties for large areas when benchmarked with slug or pumping tests. This study combined aquifer testing with specific capacity data to estimate hydraulic properties in a large alluvial aquifer. New hydrological insights for region : In the Wet Mountain Valley, Colorado, both slug tests and pumping tests were conducted, resulting in a likely range of hydraulic-conductivity values. Aquifer-testing results were related to specific-capacity data, a more spatially distributed dataset, to expand the area of aquifer characterization beyond the distribution of wells included in aquifer testing. Specific-capacity data were used in two ways: (1) a regression was built between specific-capacity values and transmissivity derived from aquifer testing; and (2) an iterative method was utilized to estimate transmissivity from specific capacity at all sites (including sites lacking aquifer tests). Study results indicate that there is a statistically significant difference between hydraulic-conductivity values estimated using the two approaches and that the regression method yields systematically greater values. These results indicate that careful consideration of methods that use specific capacity for extrapolating aquifer properties is warranted as bias could be introduced depending on the applied methodology.

Colorado↗

Stream temperature prediction in a shifting environment: The influence of deep learning architecture

Stream temperature is a fundamental control on ecosystem health. Recent efforts incorporating process guidance into deep learning models for predicting stream temperature have been shown to outperform existing statistical and physical models. This performance is in part because deep learning architectures can actively learn spatiotemporal relationships that govern how water and energy propagate through a river network. However, exploration of how spatiotemporal awareness and process guidance influence a model's generalizability under shifting environmental conditions such as climate change is limited. Here, we use Explainable Artificial Intelligence (XAI) to interrogate how differing deep learning architectures affect a model's learned spatial and temporal dependencies, and how those learned dependencies affect a model's ability to maintain high accuracy when applied to unseen environmental conditions. Using the Delaware River Basin in the northeastern United States as a test case, we compare two spatiotemporally aware process-guided deep learning models for predicting stream temperature (a recurrent graph convolution network—RGCN, and a temporal convolution graph model—Graph WaveNet). Both models achieve equally high predictive performance when testing data are well represented in the training data (test root mean squared errors of 1.64°C and 1.65°C); however, Graph WaveNet significantly outperforms RGCN in 4 out of 5 experiments where test partitions represent different types of unseen environmental conditions. XAI results show that the architecture of Graph WaveNet leads to learned spatial relationships with greater fidelity to physical processes, and that this fidelity improves the generalizability of the model when applied to shifting and/or unseen environmental conditions.

Delaware River Basin↗

Effects of wildfire on soil hydraulic properties in the western Oregon Cascades

Wildfires can substantially impact the hydrology of forested watersheds, increasing the risk of hydrologic hazards such as flash floods and debris flows. Soil hydraulic properties related to infiltration are a key control in determining the timing and magnitude of these hydrogeomorphic events. In our study, we collected 445 soil cores from burned (216 cores) and unburned (229 cores) reference catchments and analyzed them for soil hydraulic properties 10 months after the 2022 Cedar Creek Fire in Oregon, USA. We observed significantly greater field-saturated hydraulic conductivity ( K fs ), sorptivity ( S ), and wetting front potential ( Ψ f ) in burned soils relative to unburned soils, with median ratios of 5.7, 4.4, and 5.0, respectively. Among low-, moderate-, and high burn severity groups, soil hydraulic properties were not statistically different. Reductions in median soil bulk density with increasing burn severity suggested an expansion of pore sizes, which may have been partially responsible for increasing K fs and S . Additionally, in some burned soil samples, the increase in soil hydraulic properties may have been partially related to a concurrent reduction in “natural background” water repellency that is characteristic of dry, unburned soils in the Western Cascades. We observed no evidence of spatial autocorrelation in K fs using semivariogram analysis. Principal component analysis paired with a k- means cluster analysis suggested that soil physical properties explained variations in soil hydraulic properties better than landscape attributes. Although there is a lack of regional results for comparison, our results trend in the opposite direction from drier, lower net primary productivity regions that are typically studied for post-wildfire soil hydraulic properties.

Oregon↗

Comparing geometric differences between Landsat Collection 1 to Collection 2 level-1 products

In late 2020 the U.S. Geological Survey (USGS) began the distribution of Landsat products associated with their collection 2 reprocessing of the archive. Several changes were implemented within the Landsat Product Generation System (LPGS) and the calibration parameters applied to the Landsat imagery for the collection 2 processing. When comparing between collection 1 and collection 2 products, radiometric and geometric differences will be present. One of the most substantial changes between the two collections was an adjustment to the ground control which made the control more accurate from an absolute and relative perspective. Some of the other changes associated with collection 2 processing were to the Digital Elevation Model (DEM) used in terrain correction, the Thermal Infrared Sensor (TIRS) relative gains, TIRS absolute calibration, Operational Land Imagery (OLI) absolute gain model, OLI relative gain, and OLI bias calculation to name a few. Although these changes also have an effect on the differences between the product associated with these two collections, the change in the ground control, although typically less than one 30-meter multispectral pixel in magnitude, will have the largest effect on the differences between the products. This change in ground control is also a spatially dynamic change that although is low in spatial frequency, is nonlinear and is not a change that can be modeled on a global or even on a local Worldwide Reference System-2 (WRS-2) path and row scale. The effects of these ground control changes will be discussed and demonstrated within this paper along with examples showing their effect on specific datasets. This paper demonstrates some of these geometric differences associated with the ground control through both the registration statistics created during product generation and through an example comparison of a set of collection 1 and collection 2 products.

Conference Paper↗

Development of a reproducible method for determining quantity of water and its configuration in a marsh landscape

Coastal Louisiana is a dynamic and ever-changing landscape. From 1956 to 2010, over 3,734 km 2 of Louisiana's coastal wetlands have been lost due to a combination of natural and human-induced activities. The resulting landscape constitutes a mosaic of conditions from highly deteriorated to relatively stable with intact landmasses. Understanding how and why coastal landscapes change over time is critical to restoration and rehabilitation efforts. Historically, changes in marsh pattern (i.e., size and spatial distribution of marsh landmasses and water bodies) have been distinguished using visual identification by individual researchers. Difficulties associated with this approach include subjective interpretation, uncertain reproducibility, and laborious techniques. In order to minimize these limitations, this study aims to expand existing tools and techniques via a computer-based method, which uses geospatial technologies for determining shifts in landscape patterns. Our method is based on a raster framework and uses landscape statistics to develop conditions and thresholds for a marsh classification scheme. The classification scheme incorporates land and water classified imagery and a two-part classification system: (1) ratio of water to land, and (2) configuration and connectivity of water within wetland landscapes to evaluate changes in marsh patterns. This analysis system can also be used to trace trajectories in landscape patterns through space and time. Overall, our method provides a more automated means of quantifying landscape patterns and may serve as a reliable landscape evaluation tool for future investigations of wetland ecosystem processes in the northern Gulf of Mexico.

Louisiana↗

Building a state-space life cycle model for naturally produced Snake River fall Chinook salmon

In 1992, Snake River basin fall Chinook salmon (Oncorhynchus tshawytscha) were listed for protection under the U.S. Endangered Species Act (NMFS 1992) and the population remained below 1000 individuals until 2000. Since then, returns from natural production has rebounded to over 20,000 spawners owing to a host of factors including reduced harvest (Peters et al. 2001), stable minimum spawning flows (Groves and Chandler 1999), summer flow augmentation (Connor et al. 2003), predator control (Beamesderfer et al. 1996), hatchery supplementation (Rosenberger et al. 2017), improved juvenile passage structures (Adams et al. 2014), summer spill operations (Perry et al. 2006; Adams et al. 2008), and periods of favorable ocean conditions and food availability (Logerwell et al. 2003; Peterson et al. 2014). Given this change in abundance coincident with numerous management actions and fluctuation in environmental drivers, quantifying which factors contributed to the observed rebound in natural production can provide critical insights into future management actions for this at-risk population. Multistage life cycle models provide a powerful analytical framework for understating how each life stage of a population contributes to population growth rate (Moussalli and Hilborn 1986; Greene and Beechie 2004). Multistage models may also be used as an analytical framework to explicitly estimate demographic parameters of a population model. This approach has an advantage over single-stage stock-recruitment models by allowing population growth rates to be partitioned among life stages rather than aggregated over an entire life cycle. Such partitioning allows for estimating 1) stage-specific density dependence, and 2) stage-specific effects of environmental factors or management actions. For example, Zabel et al. (2006) estimated parameters of a multistage model used in the context of a population viability analysis for spring/summer Chinook salmon in the Snake River, but such an approach has yet to be applied to fall Chinook salmon in the Snake River basin. Typically, data informing estimates of abundance at particular “check points” in the life cycle determines the complexity of the multistage model that can be fit to the data. For fall Chinook salmon, we are developing a two-stage model that encompasses: 1) upstream passage of spawners at Lower Granite Dam (LGR) to the subsequent downstream passage of their progeny at the dam, and 2) downstream passage of juveniles at LGR to their subsequent return from the ocean and passage at the Dam 2‒6 years later. This approach partitions the life cycle of fall Chinook salmon both spatially and temporally, which allows us to fit and compare alternative models with covariates specific to each stage. Our previous report to the ISAB (Zabel et al. 2013) detailed methods for estimating abundance of naturally produced adults and juveniles passing Lower Granite Dam, which provides the requisite data for fitting a two-stage model. The intent of this report is to describe the structure of the two-stage life cycle model, present preliminary results from fitting the model to data, and outline future directions and developments. As is clear from the diversity of models presented in this report, “life cycle models” range from very simple theoretically based population models (e.g., the Beverton-Holt stock- recruitment model) to very complex spatially explicit simulation models linked to hydrosystem hydrodynamic models (e.g., the COMPASS model for a single transition in a life cycle model, Zabel et al. 2008). We chose to develop a model of intermediate complexity that casts the two- stage life cycle model in a state-space framework (Newman et al. 2014). We chose to use a state-space framework implemented in a Bayesian framework because: • It provides both a statistical estimation framework for retrospective statistical analysis and a stochastic simulation framework for prospective analysis to evaluate alternative management actions. • Abundance estimates are uncertain. A state-space framework accounts for observation uncertainty in the abundance estimates and other data (e.g., age structure) while simultaneously estimating process uncertainty. • It allows for missing data. By drawing missing data from an appropriate probability model, uncertainty owing to missing data can be propagated without having to omit data or assume fixed values for missing data. Thus, a two-stage state-space life cycle model for fall Chinook salmon strikes an appropriate balance between model complexity, tractability, and applicability given the goals of performing both retrospective and prospective analysis to guide future management of this population.

Idaho, Oregon, Washington, Wyoming↗

Hierarchical animal movement models for population-level inference

New methods for modeling animal movement based on telemetry data are developed regularly. With advances in telemetry capabilities, animal movement models are becoming increasingly sophisticated. Despite a need for population-level inference, animal movement models are still predominantly developed for individual-level inference. Most efforts to upscale the inference to the population level are either post hoc or complicated enough that only the developer can implement the model. Hierarchical Bayesian models provide an ideal platform for the development of population-level animal movement models but can be challenging to fit due to computational limitations or extensive tuning required. We propose a two-stage procedure for fitting hierarchical animal movement models to telemetry data. The two-stage approach is statistically rigorous and allows one to fit individual-level movement models separately, then resample them using a secondary MCMC algorithm. The primary advantages of the two-stage approach are that the first stage is easily parallelizable and the second stage is completely unsupervised, allowing for an automated fitting procedure in many cases. We demonstrate the two-stage procedure with two applications of animal movement models. The first application involves a spatial point process approach to modeling telemetry data, and the second involves a more complicated continuous-time discrete-space animal movement model. We fit these models to simulated data and real telemetry data arising from a population of monitored Canada lynx in Colorado, USA.

Environmetrics↗

Volcano deformation: Insights into magmatic systems

Volcano geodesy is the branch of geodetic science that deals with the changing shapes of volcanoes, whether large or small, deep-seated or surficial. Together with seismicity and volcanic gas flux, deformation of the ground surface can be a key indicator of subsurface conditions and processes at volcanoes—information that not only improves scientific understanding of magmatic systems but also is useful for assessing volcano hazards and mitigating their potential consequences. To take full advantage of such information requires detailed characterization of the deformation field in space and time. Currently, no single geodetic technique is capable of providing both the high spatial and temporal resolution required. However, important advances have been made recently by combining information from real-time in situ sensors such as continuous GPS, strainmeters, and tiltmeters with repeated campaign-style GPS, microgravity, interferometric synthetic aperture radar (InSAR), lidar, and photogrammetric observations. Continuous real-time data constrain the timing but not necessarily the spatial extent and pattern of deformation, whereas InSAR provides detailed spatial information but only at intervals of several days to weeks. Repeated microgravity surveys, when combined with independent measurements of surface height, are uniquely sensitive to changes in subsurface mass distribution and therefore can be used to distinguish among processes driven by magma, hydrous fluids, or gas. Simultaneous analysis of multiple geodetic datasets, especially when guided by information from global volcano databases and numerical simulations of volcanic processes and products, can provide a statistical basis for outcome prediction and serve as a guide for additional observations. In the foreseeable future, interactions among magmatic, tectonic, and hydrothermal systems will be monitored and modeled at regional scale in real time, enabling new insights into Earth’s subsurface environment.

Book chapter↗

Regional target loads of atmospheric nitrogen and sulfur deposition for the protection of stream and watershed soil resources of the Adirondack Mountains, USA

Acidic deposition contributes to a range of environmental impacts across forested landscapes, including acidification of soil and drainage water, toxic aluminum mobilization, depletion of available soil nutrient cations, and impacts to forest and aquatic species health and biodiversity. In response to decreasing levels of acidic deposition, soils and drainage waters in some regions of North America have become gradually less acidic. Thresholds of atmospheric deposition at which adverse ecological effects are manifested are called critical loads (CLs) and/or target loads (TLs). Target loads are developed based on approaches that account for spatial and temporal aspects of acidification and recovery. Exceedance represents the extent to which current or projected future levels of acidic deposition exceed the level expected to cause ecological harm. We report TLs of sulfur (S) and nitrogen (N) deposition and the potential for ecosystem recovery of watershed soils and streams in the Adirondack region of New York State, resources that have been less thoroughly investigated than lakes. Regional TLs were calculated by statistical extrapolation of hindcast and forecast simulations of 25 watersheds using the process-based model PnET-BGC coupled with empirical observations of stream hydrology and established sensitivity of sugar maple ( Acer saccharum ) to soil base saturation and brook trout ( Salvelinus fontinalis ) to stream acid neutralizing capacity (ANC). Historical impacts and the expected recovery timeline of regional soil and stream chemistry and fish community condition within the Adirondack Park were evaluated. Analysis suggests that many low-order Adirondack streams and associated watershed soils have low TLs (<40 meq/m 2 /yr of N+S deposition) to achieve specified benchmarks for recovery of soil base saturation or stream ANC. Acid-sensitive headwater and low-order streams and watershed soils in the region are expected to experience continued adverse effects from N and S deposition well into the future even under aggressive emissions reductions. Watershed soils and streams in the western Adirondack Park are particularly vulnerable to acidic deposition and currently in exceedance of TLs. The methods used for linking statistical and process-based models to consider chemical and biological response under varying flow conditions at the regional scale in this study can be applied to other areas of concern.

New York↗

Estimating effects of limiting factors with regression quantiles

In a recent Concepts paper in Ecology, Thomson et al. emphasized that assumptions of conventional correlation and regression analyses fundamentally conflict with the ecological concept of limiting factors, and they called for new statistical procedures to address this problem. The analytical issue is that unmeasured factors may be the active limiting constraint and may induce a pattern of unequal variation in the biological response variable through an interaction with the measured factors. Consequently, changes near the maxima, rather than at the center of response distributions, are better estimates of the effects expected when the observed factor is the active limiting constraint. Regression quantiles provide estimates for linear models fit to any part of a response distribution, including near the upper bounds, and require minimal assumptions about the form of the error distribution. Regression quantiles extend the concept of one-sample quantiles to the linear model by solving an optimization problem of minimizing an asymmetric function of absolute errors. Rank-score tests for regression quantiles provide tests of hypotheses and confidence intervals for parameters in linear models with heteroscedastic errors, conditions likely to occur in models of limiting ecological relations. We used selected regression quantiles (e.g., 5th, 10th, ..., 95th) and confidence intervals to test hypotheses that parameters equal zero for estimated changes in average annual acorn biomass due to forest canopy cover of oak (Quercus spp.) and oak species diversity. Regression quantiles also were used to estimate changes in glacier lily (Erythronium grandiflorum) seedling numbers as a function of lily flower numbers, rockiness, and pocket gopher (Thomomys talpoides fossor) activity, data that motivated the query by Thomson et al. for new statistical procedures. Both example applications showed that effects of limiting factors estimated by changes in some upper regression quantile (e.g., 90-95th) were greater than if effects were estimated by changes in the means from standard linear model procedures. Estimating a range of regression quantiles (e.g., 5-95th) provides a comprehensive description of biological response patterns for exploratory and inferential analyses in observational studies of limiting factors, especially when sampling large spatial and temporal scales.

Ecology↗

Status and trends of total nitrogen and total phosphorus concentrations, loads, and yields in streams of Mississippi, water years 2008–18

To assess the status and trends of conditions of surface waters throughout Mississippi, the U.S. Geological Survey, in cooperation with the Mississippi Department of Environmental Quality (MDEQ), summarized concentrations and estimated loads, yields, trends, and spatial and temporal patterns of total nitrogen (TN) and total phosphorus (TP) at 20 stream sites in MDEQ’s ambient water-quality monitoring network and 2 stream sites in the U.S. Geological Survey’s National Water-Quality Assessment Project’s monitoring network. Comparison of streamflow at the time of water-quality sample collection to flow-duration curves for each site showed that samples were relatively evenly spread over a wide range of flows, indicating that load estimations were representative of a wide range of flows. Relation of streamflow to concentrations of TN and TP varied among sites and land use. Sites with high agriculture land use in the drainage basin tended to have a positive correlation between streamflow and concentration, suggesting influence of event-driven nonpoint-source runoff. Sites near urban (developed) areas tended to have a negative correlation between streamflow and concentration, suggesting chronic point-source influences during low-flow conditions. Sites with high forest land use and lower agriculture and urban (developed) land use showed little to no association between streamflow and concentration. Seasonal distributions of concentrations of TN and TP also corresponded closely with variations in land use. Sites near urban (developed) land had the highest concentrations in late summer and fall, sites with a high percentage of agricultural land had the highest concentrations in the spring, and sites that were primarily forested or with little developed land did not exhibit substantial changes in concentration across seasons. Eight sites had statistical likelihoods for upward trends of TN loads, and seven sites had statistical likelihoods for downward trends. Trends in TN loads at six sites were considered “about as likely as not,” meaning that a site has an equal chance of having an upward or downward trend. Trend results of mean annual flow-normalized loads of TP for the period of analysis (2008–18) showed that 16 sites had upward trends, 3 sites had downward trends, and 2 sites were considered “about as likely as not.” Results from our study were compared to results from existing regional models to assess accuracy of predictions at a local scale. Comparisons of yields predicted from 2012 regional-scale SPAtially Referenced Regressions on Watershed attributes (SPARROW) to results from this study showed the 2012 SPARROW-predicted estimates varied in consistency with results from this study. The 2012 SPARROW-prediction model underestimated TN yields, more often and by a slightly larger degree, more than it overestimated TN yields. The 2012 SPARROW-predicted model tended to underestimate yields at study sites with higher yields. All four sites in the predominantly agricultural area of northwest Mississippi, locally known as the Mississippi Delta, were underestimated by 2012 SPARROW. For TP, yield comparisons at sites with lower yields were consistent, yields at sites with midrange yields tended to be overestimated by SPARROW, and yields at sites with high yields tended to be underestimated by SPARROW. TP yields at four sites in the Mississippi Delta were underestimated by the 2012 SPARROW-predicted model. Results of select sites from our study were also compared to other published load estimates from an earlier time period to evaluate possible trends. Comparison of TN yields at four sites and TP yields at three sites from the study-derived estimates to estimates made from data spanning 1993–2004 showed decreasing TN yields at all four sites and decreasing TP yields at two of three sites, with increasing yields of TP at the Yazoo River lower site. Also, a third comparison of the TN and TP yields of the Yazoo River lower site of this study to estimates made from data spanning 1996–97 showed decreasing TN yields but similar TP yields. This suggests that TN yields may have decreased over the last 20–30 years, but TP yields remain constant or are increasing.

Mississippi↗

A guide to Bayesian model checking for ecologists

Checking that models adequately represent data is an essential component of applied statistical inference. Ecologists increasingly use hierarchical Bayesian statistical models in their research. The appeal of this modeling paradigm is undeniable, as researchers can build and fit models that embody complex ecological processes while simultaneously accounting for observation error. However, ecologists tend to be less focused on checking model assumptions and assessing potential lack of fit when applying Bayesian methods than when applying more traditional modes of inference such as maximum likelihood. There are also multiple ways of assessing the fit of Bayesian models, each of which has strengths and weaknesses. For instance, Bayesian P values are relatively easy to compute, but are well known to be conservative, producing P values biased toward 0.5. Alternatively, lesser known approaches to model checking, such as prior predictive checks, cross‐validation probability integral transforms, and pivot discrepancy measures may produce more accurate characterizations of goodness‐of‐fit but are not as well known to ecologists. In addition, a suite of visual and targeted diagnostics can be used to examine violations of different model assumptions and lack of fit at different levels of the modeling hierarchy, and to check for residual temporal or spatial autocorrelation. In this review, we synthesize existing literature to guide ecologists through the many available options for Bayesian model checking. We illustrate methods and procedures with several ecological case studies including (1) analysis of simulated spatiotemporal count data, (2) N‐mixture models for estimating abundance of sea otters from an aircraft, and (3) hidden Markov modeling to describe attendance patterns of California sea lion mothers on a rookery. We find that commonly used procedures based on posterior predictive P values detect extreme model inadequacy, but often do not detect more subtle cases of lack of fit. Tests based on cross‐validation and pivot discrepancy measures (including the “sampled predictive P value”) appear to be better suited to model checking and to have better overall statistical performance. We conclude that model checking is necessary to ensure that scientific inference is well founded. As an essential component of scientific discovery, it should accompany most Bayesian analyses presented in the literature.

Ecological Monographs↗

A high-resolution, daily hindcast (1990-2021) of Alaskan river discharge and temperature from coupled and optimized physical models

Water quality and freshwater ecosystems are affected by river discharge and temperature. Models are frequently used to estimate river temperature on large spatial and temporal scales due to limited observations of discharge and temperature. In this study, we use physically based river routing and temperature models to simulate daily discharge and river temperature for rivers in 138 basins in Alaska, including the entire Yukon River basin, from 1990–2021. The river temperature model was optimized for ice free months using a surrogate-based model optimization method, improving model performance at uncalibrated river gages. A common statistical model relating local air and water temperature was used as a benchmark. The physically based river temperature model exhibited superior performance compared to the benchmark statistical model after optimization, suggesting river temperature model optimization could become more routine. The river temperature model demonstrated high sensitivity to air temperature and model parameterization, and lower sensitivity to discharge. Validation of the models showed a Kling-Gupta Efficiency of 0.46 for daily river discharge and a root mean square error of 2.04°C for daily river temperature, improving on the non-optimized physical model and the benchmark statistical model, which had root mean square errors of 3.24 and 2.97°C, respectively. The simulation shows that rivers in northern Alaska have higher maximum summer temperatures and more variability than rivers in the Central and Southern regions. Furthermore, this framework can be readily adapted for use across models and regions.

Alaska↗