SEARCH · USGS Science

Results for “Algorithms”

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

SourceResearch connectionFrom yearThrough year

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

At least 721 records · Page 40Linked to original sources

USGS2023

Evaluating a tandem human-machine approach to labelling of wildlife in remote camera monitoring

Remote cameras (“trail cameras”) are a popular tool for non-invasive, continuous wildlife monitoring, and as they become more prevalent in wildlife research, machine learning (ML) is increasingly used to automate or accelerate the labor-intensive process of labelling (i.e., tagging) photos. Human-machine hybrid tagging approaches have been shown to greatly increase tagging efficiency (i.e., time to tag a single image). However, those potential increases hinge on the extent to which an ML model makes correct vs. incorrect predictions. We performed an experiment using a ML model that produces bounding boxes around animals, people, and vehicles in remote camera imagery (MegaDetector) to consider the impact of a ML model’s performance on its ability to accelerate human labeling. Six participants tagged trail camera images collected from 12 sites in Vermont and Maine, USA (January–September 2022) using three tagging methods (one with ML bounding box assistance and two without assistance). We used a generalized linear mixed model to examine the influence of ML model performance and tagging method on tagging efficiency. We found that ML bounding boxes offer significant improvement in tagging efficiency when labelling data compared to unassisted tagging. Additionally, the time taken to label with bounding boxes was not statistically different from an unassisted tagging approach. However, we found that gains in efficiency are contingent on the ML algorithm’s performance and that incorrect ML predictions, particularly the 4.2% false positive and 3.6% false negative predictions, can slow the tagging process compared to a non-hybrid approach. These findings indicate that although practitioners usually forgo the production of bounding boxes when selecting a data labelling process due to the increased effort, ML bounding box-assisted tagging can offer an efficient method for labeling. More broadly, ML-assisted data labelling offers an opportunity to accelerate the analysis of trail camera imagery, but an assessment of the ML model’s performance can illuminate whether the hybrid-tagging approach is ultimately a help or hinderance.

Maine, Vermont↗

USGS2013

Evapotranspiration and water balance of an anthropogenic coastal desert wetland: responses to fire, inflows and salinities

Evapotranspiration (ET) and other water balance components were estimated for Cienega de Santa Clara, an anthropogenic brackish wetland in the delta of the Colorado River in Mexico. The marsh is in the Biosphere Reserve of the Upper Gulf of California and Delta of the Colorado River, and supports a high abundance and diversity of wildlife. Over 95% of its water supply originates as agricultural drain water from the USA, sent for disposal in Mexico. This study was conducted from 2009 to 2011, before, during and after a trial run of the Yuma Desalting Plant in the USA, which will divert water from the wetland and replace it with brine from the desalting operation. The goal was to estimate the main components in the water budget to be used in creating management scenarios for this marsh. We used a remote sensing algorithm to estimate ET from meteorological data and Enhanced Vegetation Index values from the Moderate Resolution Imaging Spectrometer (MODIS) sensors on the Terra satellite. ET estimates from the MODIS method were then compared to results from a mass balance of water and salt inflows and outflows over the study period. By both methods, mean annual ET estimates ranged from 2.6 to 3.0 mm d −1 , or 50 to 60% of reference ET (ET o ). Water entered at a mean salinity of 2.6 g L −1 TDS and mean salinity in the wetland was 3.73 g L −1 TDS over the 33 month study period. Over an annual cycle, 54% of inflows supported ET while the rest exited the marsh as outflows; however, in winter when ET was low, up to 90% of the inflows exited the marsh. An analysis of ET estimates over the years 2000–2011 showed that annual ET was proportional to the volume of inflows, but was also markedly stimulated by fires. Spring fires in 2006 and 2011 burned off accumulated thatch, resulting in vigorous growth of new leaves and a 30% increase in peak summer ET compared to non-fire years. Following fires, peak summer ET estimates were equal to ET o , while in non-fire years peak ET was equal to only one-half to two-thirds of ET o . Over annual cycles, estimated ET was always lower than ET o , because T. domingensis is dormant in winter and shades the water surface, reducing direct evaporation. Thus, ET of a Typha marsh is likely to be less than an open water surface under most conditions.

Ecological Engineering↗

USGS2018

Estimating vegetation biomass and cover across large plots in shrub and grass dominated drylands using terrestrial lidar and machine learning

Terrestrial laser scanning (TLS) has been shown to enable an efficient, precise, and non-destructive inventory of vegetation structure at ranges up to hundreds of meters. We developed a method that leverages TLS collections with machine learning techniques to model and map canopy cover and biomass of several classes of short-stature vegetation across large plots. We collected high-definition TLS scans of 26 1-ha plots in desert grasslands and big sagebrush shrublands in southwest Idaho, USA. We used the Random Forests machine learning algorithm to develop decision tree models predicting the biomass and canopy cover of several vegetation classes from statistical descriptors of the aboveground heights of TLS points. Manual measurements of vegetation characteristics collected within each plot served as training and validation data. Models based on five or fewer TLS descriptors of vegetation heights were developed to predict the canopy cover fraction of shrubs (R 2 = 0.77, RMSE = 7%), annual grasses (R 2 = 0.70, RMSE = 21%), perennial grasses (R 2 = 0.36, RMSE = 12%), forbs (R 2 = 0.52, RMSE = 6%), bare earth or litter (R 2 = 0.49, RMSE = 19%), and the biomass of shrubs (R 2 = 0.71, RMSE = 175 g) and herbaceous vegetation (R 2 = 0.61, RMSE = 99 g) (all values reported are out-of-bag). Our models explained much of the variability between predictions and manual measurements, and yet we expect that future applications could produce even better results by reducing some of the methodological sources of error that we encountered. Our work demonstrates how TLS can be used efficiently to extend manual measurement of vegetation characteristics from small to large plots in grasslands and shrublands, with potential application to other similarly structured ecosystems. Our method shows that vegetation structural characteristics can be modeled without classifying and delineating individual plants, a challenging and time-consuming step common in previous methods applying TLS to vegetation inventory. Improving application of TLS to studies of shrub-steppe ecosystems will serve immediate management needs by enhancing vegetation inventories, environmental modeling studies, and the ability to train broader datasets collected from air and space.

Idaho↗

USGS2021

Are drought indices and climate data good indicators of ecologically relevant soil moisture dynamics in drylands?

Droughts are disproportionately impacting global dryland regions where ecosystem health and function are tightly coupled to moisture availability. Drought severity is commonly estimated using algorithms such as the standardized precipitation-evapotranspiration index (SPEI), which can estimate climatic water balance impacts at various hydrologic scales by varying computational length. However, the performance of these metrics as indicators of soil moisture dynamics at ecologically relevant scales, across soil depths, and in consideration of broader scale ecohydrological processes, requires more attention. In this study, we tested components of climatic water balance, including SPEI and SPEI computation lengths, to recreate multi-decadal and periodic soil-moisture patterns across soil profiles at 866 sites in the western United States. Modeling results show that SPEI calculated over the prior 12-months was the most predictive computation length and could recreate changes in moisture availability within the soil profile over longer periods of time and for annual recharge of deeper soil moisture stores. SPEI was slightly less successful with recreating spring surface-soil moisture availability, which is key to dryland ecosystems dominated by winter precipitation. Meteorological drought indices like SPEI are intended to be convenient and generalized indicators of meteorological water deficit. However, the inconsistent ability of SPEI to recreate ecologically relevant patterns of soil moisture at regional scales suggests that process-based models, and the larger data requirements they involve, remain an important tool for dryland ecohydrology

California, Idaho, Nevada, Oregon, Utah↗

USGS2022

Satellites quantify the spatial extent of cyanobacterial blooms across the United States at multiple scales

Previous studies indicate that cyanobacterial harmful algal bloom (cyanoHAB) frequency, extent, and magnitude have increased globally over the past few decades. However, little quantitative capability is available to assess these metrics of cyanoHABs across broad geographic scales and at regular intervals. Here, the spatial extent was quantified from a cyanobacteria algorithm applied to two European Space Agency satellite platforms—the MEdium Resolution Imaging Spectrometer (MERIS) onboard Envisat and the Ocean and Land Colour Instrument (OLCI) onboard Sentinel-3. CyanoHAB spatial extent was defined for each geographic area as the percentage of valid satellite pixels that exhibited cyanobacteria above the detection limit of the satellite sensor. This study quantified cyanoHAB spatial extent for over 2,000 large lakes and reservoirs across the contiguous United States (CONUS) during two time periods: 2008–2011 via MERIS and 2017–2020 via OLCI when cloud-, ice-, and snow-free imagery was available. Approximately 56% of resolvable lakes were glaciated, 13% were headwater, isolated, or terminal lakes, and the rest were primarily drainage lakes. Results were summarized at national-, regional-, state-, and lake-scales, where regions were defined as nine climate regions which represent climatically consistent states. As measured by satellite, changes in national cyanoHAB extent did have a strong increase of 6.9% from 2017 to 2020 (|Kendall’s tau (τ)| = 0.56; gamma (γ) = 2.87 years), but had negligible change (|τ| = 0.03) from 2008 to 2011. Two of the nine regions had moderate (0.3 ≤ |τ| < 0.5) increases in spatial extent from 2017 to 2020, and eight of nine regions had negligible (|τ| < 0.2) change from 2008 to 2011. Twelve states had a strong or moderate increase from 2017 to 2020 (|τ| ≥ 0.3), while only one state had a moderate increase and two states had a moderate decrease from 2008 to 2011. A decrease, or no change, in cyanoHAB spatial extent did not indicate a lack of issues related to cyanoHABs. Sensitivity results of randomly omitted daily CONUS scenes confirm that even with reduced data availability during a short four-year temporal assessment, the direction and strength of the changes in spatial extent remained consistent. We present the first set of national maps of lake cyanoHAB spatial extent across CONUS and demonstrate an approach for quantifying past and future changes at multiple spatial scales. Results presented here provide water quality managers information regarding current cyanoHAB spatial extent and quantify rates of change.

Ecological Indicators↗

USGS2024

Passive acoustic monitoring and convolutional neural networks facilitate high-resolution and broadscale monitoring of a threatened species

Population monitoring is an essential component of biodiversity conservation and management, but low detection probabilities for rare and/or cryptic species makes estimating abundance and occupancy challenging. Passive acoustic monitoring combined with machine learning algorithms represents a potential path forward to effectively and efficiently monitor the occurrence of rare vocalizing species across entire forest landscapes. Our objectives were to develop and implement a convolutional neural network (PNW-Cnet) to identify vocalizations of a rare and threatened forest nesting bird species – the marbled murrelet ( Brachyramphus marmoratus ) – in the Pacific Northwest, U.S.A., 2018–2021. We used PNW-Cnet predictions from broadscale passive acoustic monitoring data to examine spatiotemporal patterns in the distribution of murrelets. PNW-Cnet showed sufficiently high prediction accuracy (overall precision > 0.9) to enable broadscale population monitoring. Spatiotemporal analysis showed that annual peak murrelet call abundance occurs in ordinal weeks 28–32 (late July–Mid August) but this varied by study area. The greatest number of detections typically occurred in the Olympic Peninsula and Oregon Coast Range where late-successional forest dominates and nearer to ocean habitats. We demonstrate that passive acoustic monitoring can be used to understand intensity of use across broad scales for a rare and cryptic species in addition to the typical detection/non-detection data that are often collected. Passive acoustic monitoring combined with PNW-Cnet offers considerable promise for species distribution modeling and long-term population monitoring for rare species.

Oregon, Washington↗

USGS2006

Predicting tree species presence and basal area in Utah: A comparison of stochastic gradient boosting, generalized additive models, and tree-based methods

Many efforts are underway to produce broad-scale forest attribute maps by modelling forest class and structure variables collected in forest inventories as functions of satellite-based and biophysical information. Typically, variants of classification and regression trees implemented in Rulequest's?? See5 and Cubist (for binary and continuous responses, respectively) are the tools of choice in many of these applications. These tools are widely used in large remote sensing applications, but are not easily interpretable, do not have ties with survey estimation methods, and use proprietary unpublished algorithms. Consequently, three alternative modelling techniques were compared for mapping presence and basal area of 13 species located in the mountain ranges of Utah, USA. The modelling techniques compared included the widely used See5/Cubist, generalized additive models (GAMs), and stochastic gradient boosting (SGB). Model performance was evaluated using independent test data sets. Evaluation criteria for mapping species presence included specificity, sensitivity, Kappa, and area under the curve (AUC). Evaluation criteria for the continuous basal area variables included correlation and relative mean squared error. For predicting species presence (setting thresholds to maximize Kappa), SGB had higher values for the majority of the species for specificity and Kappa, while GAMs had higher values for the majority of the species for sensitivity. In evaluating resultant AUC values, GAM and/or SGB models had significantly better results than the See5 models where significant differences could be detected between models. For nine out of 13 species, basal area prediction results for all modelling techniques were poor (correlations less than 0.5 and relative mean squared errors greater than 0.8), but SGB provided the most stable predictions in these instances. SGB and Cubist performed equally well for modelling basal area for three species with moderate prediction success, while all three modelling tools produced comparably good predictions (correlation of 0.68 and relative mean squared error of 0.56) for one species. ?? 2006 Elsevier B.V. All rights reserved.

Ecological Modelling↗

USGS2008

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↗

USGS2011

Resolving structural uncertainty in natural resources management using POMDP approaches

In recent years there has been a growing focus on the uncertainties of natural resources management, and the importance of accounting for uncertainty in assessing management effectiveness. This paper focuses on uncertainty in resource management in terms of discrete-state Markov decision processes (MDP) under structural uncertainty and partial observability. It describes the treatment of structural uncertainty with approaches developed for partially observable resource systems. In particular, I show how value iteration for partially observable MDPs (POMDP) can be extended to structurally uncertain MDPs. A key difference between these process classes is that structurally uncertain MDPs require the tracking of system state as well as a probability structure for the structure uncertainty, whereas with POMDPs require only a probability structure for the observation uncertainty. The added complexity of the optimization problem under structural uncertainty is compensated by reduced dimensionality in the search for optimal strategy. A solution algorithm for structurally uncertain processes is outlined for a simple example in conservation biology. By building on the conceptual framework developed for POMDPs, natural resource analysts and decision makers who confront structural uncertainties in natural resources can take advantage of the rapid growth in POMDP methods and approaches, and thereby produce better conservation strategies over a larger class of resource problems.

Ecological Modelling↗

USGS2016

Effects of model complexity and priors on estimation using sequential importance sampling/resampling for species conservation

We examined the effects of complexity and priors on the accuracy of models used to estimate ecological and observational processes, and to make predictions regarding population size and structure. State-space models are useful for estimating complex, unobservable population processes and making predictions about future populations based on limited data. To better understand the utility of state space models in evaluating population dynamics, we used them in a Bayesian framework and compared the accuracy of models with differing complexity, with and without informative priors using sequential importance sampling/resampling (SISR). Count data were simulated for 25 years using known parameters and observation process for each model. We used kernel smoothing to reduce the effect of particle depletion, which is common when estimating both states and parameters with SISR. Models using informative priors estimated parameter values and population size with greater accuracy than their non-informative counterparts. While the estimates of population size and trend did not suffer greatly in models using non-informative priors, the algorithm was unable to accurately estimate demographic parameters. This model framework provides reasonable estimates of population size when little to no information is available; however, when information on some vital rates is available, SISR can be used to obtain more precise estimates of population size and process. Incorporating model complexity such as that required by structured populations with stage-specific vital rates affects precision and accuracy when estimating latent population variables and predicting population dynamics. These results are important to consider when designing monitoring programs and conservation efforts requiring management of specific population segments.

Ecological Modelling↗

USGS2022

A three-dimensional Lagrangian particle tracking model for predicting transport of eggs of rheophilic-spawning carps in turbulent rivers

Grass carp, bighead carp, and silver carp spawn in flowing water. Their eggs, and then larvae, develop while drifting. Hydraulic conditions and water temperature control spawning locations, egg survival, and the downstream distance traveled before the hatched larvae can swim for low velocity nursery habitats. Existing egg drift models simulate the fluvial transport of carp eggs but have limitations in capturing the effect of localized turbulence on egg transport due to inadequate dimensions of hydrodynamics and/or empirical parameterization of river dispersion. We present a three-dimensional Lagrangian particle tracking model that uses fully resolved river hydrodynamics and a continuous random walk algorithm driven by local turbulent kinetic energy and its dissipation rate. We incorporate a new set of equations to compute evolving egg characteristics with fully resolved 3-D hydrodynamics. To demonstrate the performance of the model, we conducted a case study in an eight-kilometer reach of the Missouri River at the discharge of approximately 25% daily flow exceedance. Three-dimensional river hydrodynamics was modeled, calibrated, and evaluated with measurement data. Egg drift was modeled and compared using fully three-dimensional, depth-averaged two-dimensional, and zone-averaged one-dimensional hydrodynamics. The comparison shows a generally good agreement among models of downstream egg transport due to advection but a different dispersion pattern of eggs in the river, as a result of turbulent diffusion and shear induced dispersion.

Ecological Modelling↗

USGS2018

Machine learning for ecosystem services

Recent developments in machine learning have expanded data-driven modelling (DDM) capabilities, allowing artificial intelligence to infer the behaviour of a system by computing and exploiting correlations between observed variables within it. Machine learning algorithms may enable the use of increasingly available ‘big data’ and assist applying ecosystem service models across scales, analysing and predicting the flows of these services to disaggregated beneficiaries. We use the Weka and ARIES software to produce two examples of DDM: firewood use in South Africa and biodiversity value in Sicily, respectively. Our South African example demonstrates that DDM (64–91% accuracy) can identify the areas where firewood use is within the top quartile with comparable accuracy as conventional modelling techniques (54–77% accuracy). The Sicilian example highlights how DDM can be made more accessible to decision makers, who show both capacity and willingness to engage with uncertainty information. Uncertainty estimates, produced as part of the DDM process, allow decision makers to determine what level of uncertainty is acceptable to them and to use their own expertise for potentially contentious decisions. We conclude that DDM has a clear role to play when modelling ecosystem services, helping produce interdisciplinary models and holistic solutions to complex socio-ecological issues.

Ecosystem Services↗

USGS2006

Seafloor habitat mapping of the New York Bight incorporating sidescan sonar data

The efficacy of using sidescan sonar imagery, image classification algorithms and geographic information system (GIS) techniques to characterize the seafloor bottom of the New York Bight were assessed. The resulting seafloor bottom type map was compared with fish trawl survey data to determine whether there were any discernable habitat associations. An unsupervised classification with 20 spectral classes was produced using the sidescan sonar imagery, bathymetry and secondarily derived spatial heterogeneity to characterize homogenous regions within the study area. The spectral classes, geologic interpretations of the study region, bathymetry and a bottom landform index were used to produce a seafloor bottom type map of 9 different bottom types. Examination of sediment sample data by bottom type indicated that each bottom type class had a distinct composition of sediments. Analysis of adult summer flounder, Paralichthys dentatus, and adult silver hake, Merluccius bilinearis, presence/absence data from trawl surveys did not show evidence of strong associations between the species distributions and seafloor bottom type. However, the absence of strong habitat associations may be more attributable to the coarse scale and geographic uncertainty of the trawl sampling data than conclusive evidence that no habitat associations exist for these two species. ?? 2006 Elsevier Ltd. All rights reserved.

Estuarine, Coastal and Shelf Science↗

USGS2015

Assessing the impacts of climate and land use and land cover change on the freshwater availability in the Brahmaputra River basin

Study Region: Brahmaputra River basin in South Asia. Study Focus: The Soil and Water Assessment Tool was used to evaluate sensitivities and patterns in freshwater availability due to projected climate and land use changes in the Brahmaputra basin. The daily observed discharge at Bahadurabad station in Bangladesh was used to calibrate and validate the model and analyze uncertainties with a sequential uncertainty fitting algorithm. The sensitivities and impacts of projected climate and land use changes on basin hydrological components were simulated for the A1B and A2 scenarios and analyzed relative to a baseline scenario of 1988–2004. New hydrological insights for the region: Basin average annual ET was found to be sensitive to changes in CO 2 concentration and temperature, while total water yield, streamflow, and groundwater recharge were sensitive to changes in precipitation. The basin hydrological components were predicted to increase with seasonal variability in response to climate and land use change scenarios. Strong increasing trends were predicted for total water yield, streamflow, and groundwater recharge, indicating exacerbation of flooding potential during August–October, but strong decreasing trends were predicted, indicating exacerbation of drought potential during May–July of the 21st century. The model has potential to facilitate strategic decision making through scenario generation integrating climate change adaptation and hazard mitigation policies to ensure optimized allocation of water resources under a variable and changing climate.

Brahmaputra River basin↗

USGS2024

Improving crop-specific groundwater use estimation in the Mississippi Alluvial Plain: Implications for integrated remote sensing and machine learning approaches in data-scarce regions

Study region The Mississippi Alluvial Plain (MAP) in the United States (US). Study focus Understanding local-scale groundwater use, a critical component of the water budget, is necessary for implementing sustainable water management practices. The MAP is one of the most productive agricultural regions in the US and extracts more than 11 km 3 /year for irrigation activities. Consequently, groundwater-level declines in the MAP region pose a substantial challenge to water sustainability, and hence, we need reliable groundwater pumping monitoring solutions to manage this resource appropriately. New hydrological insights for the region We incorporate remote sensing datasets and machine learning to improve an existing lookup table-based model of groundwater use previously developed by the U.S. Geological Survey (USGS). Here, we employ Distributed Random Forests, an ensemble machine learning algorithm to predict annual and monthly groundwater use (2014–2020) throughout this region at 1-km resolution, using pumping data from existing flowmeters in the Mississippi Delta. Our model compares favorably with the existing USGS model, with higher R 2 (0.51 compared to 0.42 in the previous model), and lower root mean square error (RMSE) and mean absolute error (MAE)— 0.14 m and 0.09 m, respectively in our model, compared to 0.15 m and 0.1 m in the previous model. Therefore, this work advances our ability to predict groundwater use in regions with scarce or limited in-situ groundwater withdrawal data availability.

Journal of Hydrology Regional Studies↗

USGS2021

Establishment of baseline cytology metrics in nestling American kestrels (Falco sparverius): Immunomodulatory effects of the flame retardant isopropylated triarylphosphate isomers

Avian populations must mount effective immune responses upon exposure to environmental stressors such as avian influenza and xenobiotics. Although multiple immune assays have been tested and applied to various avian species, antibody-mediated immune responses in non-model avian species are not commonly reported due to the lack of commercially available species-specific antibodies. The objectives of the present study were to advance methods for studying wild bird immune responses and to apply these to the evaluation of cytological responses after exposure of American kestrels, Falco sparverius, to a commercial flame retardant mixture containing isopropylated triarylphosphate isomers (ITP). Hatchlings were gavaged daily with safflower oil or 1.5 ug/g bw/day of ITP suspended in safflower oil, then bled on days 9, 17, and 21. The ITP treatment group ( n = 18) and a subset of controls (Poly I:C treatment group; n = 10) were injected on days 9 and 15 with a synthetic analog of viral double-stranded RNA, polyinosinic:polycytidylic acid (Poly I:C), a toll-like receptor ligand and synthetic viral mimic, and responses compared to a sham injected control group (n = 8). The hypotheses tested whether kestrels showed immunological differences among treatment groups, genetic sex, and/or white blood cell (WBC) subpopulation type over time. A flow cytometry (FCM) gating strategy categorized heterophils (H), lymphocytes (L), and monocytes (M) and their proportions, and measured relative fluorescence in response to anti-chicken CD4 binding. Fluorescent cell surfaces and some granular/vacuolar inclusions were visualized by epifluorescence microscopy. A fourth subpopulation with higher levels of granularity than M but less than H became increasingly apparent with time and was gated along with the H subpopulation; its frequency of occurrence was lowest in the ITP group ( P = 0.0023). The percentages of cells differed among treatment groups, days, and sexes ( P = 0.0001). For both sexes, percentages of H and L were higher than M in control and Poly I:C. In the ITP group, L percentages were higher than H and M ( P = 0.0457), and H and L were higher than M on days 9 and 21 ( P = 0.0001). The ratios of H:L and H:WBC, indicators of robust immunity, were also higher on days 9 and 21 than on 17 ( P = 0.0079). For each sex, the highest levels of activity measured by FCM geometric means (GEO) of fluorescence (indicative of antibody binding) were observed on day 9 ( P = 0.0001 female, and P = 0.0011 male) in H over both L and M ( P < 0.0001 for each). In males, GEO of the Poly I:C group was higher than that of the ITP group ( P = 0.0374), with no difference observed among females over all days. By using a FCM algorithm for population comparisons of fluorescence to investigate binding within H, the T(x) scores indicated higher fluorescence in control and Poly I:C groups over ITP ( P = 0.0001). Unlike chickens, Gallus gallus , which express CD4 primarily on L, kestrels bound the commercial antibody primarily within the gated H subpopulation, suggesting an immunophenotypic difference between taxa, despite a ~60% identity of Falco CD4 amino acid sequences with chicken CD4. The emergent cell subset within the gated H presented dendritic-like cell (DLC) morphological and functional properties, apparently serving as an effector cell. This study adds interpretive context to ecological investigations of infection and of potential immunomodulation by emerging compounds, whereby the early innate responses are mediated by the various cell subsets serving as useful quantitative markers of immunological condition. Data showed that dietary exposure to ITP was immunosuppressive for male and female kestrels over the course of the experiment, reducing DLC frequency compared to the Poly I:C controls. Heterophils and DLC were important in facilitating innate immunological responses.

Environment International↗

USGS2005

Parallel Fortran-MPI software for numerical inversion of the Laplace transform and its application to oscillatory water levels in groundwater environments

A parallel Fortran-MPI (Message Passing Interface) software for numerical inversion of the Laplace transform based on a Fourier series method is developed to meet the need of solving intensive computational problems involving oscillatory water level's response to hydraulic tests in a groundwater environment. The software is a parallel version of ACM (The Association for Computing Machinery) Transactions on Mathematical Software (TOMS) Algorithm 796. Running 38 test examples indicated that implementation of MPI techniques with distributed memory architecture speedups the processing and improves the efficiency. Applications to oscillatory water levels in a well during aquifer tests are presented to illustrate how this package can be applied to solve complicated environmental problems involved in differential and integral equations. The package is free and is easy to use for people with little or no previous experience in using MPI but who wish to get off to a quick start in parallel computing. ?? 2004 Elsevier Ltd. All rights reserved.

Environmental Modelling and Software↗

USGS2014

A computer program for uncertainty analysis integrating regression and Bayesian methods

This work develops a new functionality in UCODE_2014 to evaluate Bayesian credible intervals using the Markov Chain Monte Carlo (MCMC) method. The MCMC capability in UCODE_2014 is based on the FORTRAN version of the differential evolution adaptive Metropolis (DREAM) algorithm of Vrugt et al. (2009), which estimates the posterior probability density function of model parameters in high-dimensional and multimodal sampling problems. The UCODE MCMC capability provides eleven prior probability distributions and three ways to initialize the sampling process. It evaluates parametric and predictive uncertainties and it has parallel computing capability based on multiple chains to accelerate the sampling process. This paper tests and demonstrates the MCMC capability using a 10-dimensional multimodal mathematical function, a 100-dimensional Gaussian function, and a groundwater reactive transport model. The use of the MCMC capability is made straightforward and flexible by adopting the JUPITER API protocol. With the new MCMC capability, UCODE_2014 can be used to calculate three types of uncertainty intervals, which all can account for prior information: (1) linear confidence intervals which require linearity and Gaussian error assumptions and typically 10s–100s of highly parallelizable model runs after optimization, (2) nonlinear confidence intervals which require a smooth objective function surface and Gaussian observation error assumptions and typically 100s–1,000s of partially parallelizable model runs after optimization, and (3) MCMC Bayesian credible intervals which require few assumptions and commonly 10,000s–100,000s or more partially parallelizable model runs. Ready access allows users to select methods best suited to their work, and to compare methods in many circumstances.

Environmental Modelling and Software↗