USGS Science⌕ Search

SEARCH · USGS Science

Results for “Machine Learning with Applications”

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 91 records · Page 5Linked to original sources
USGS2019

Context matters: Using reinforcement learning to develop human-readable, state-dependent outbreak response policies

The number of all possible epidemics of a given infectious disease that could occur on a given landscape is large for systems of real-world complexity. Furthermore, there is no guarantee that the control actions that are optimal, on average, over all possible epidemics are also best for each possible epidemic. Reinforcement learning (RL) has been used to develop machine-readable context-dependent solutions for complex problems with many possible realisations ranging from video-games to the game of Go. RL could be a valuable tool to generate context-dependent policies for outbreak response, though translating the resulting policies into simple rules that can be read and interpreted by human decision-makers remains a challenge. Here we illustrate the application of RL to the development of context-dependent outbreak response policies to minimise outbreaks of foot-and-mouth disease. We show that control based on the resulting context-dependent policies, which adapt interventions to the specific outbreak, result in smaller outbreaks than static policies. We further illustrate two approaches for translating the complex machine-readable policies into simple heuristics that can be evaluated by human decision-makers.

Philosophical Transactions of the Royal Society B:↗
USGS2021

Global cropland-extent product at 30-m resolution (GCEP30) derived from Landsat satellite time-series data for the year 2015 using multiple machine-learning algorithms on Google Earth Engine cloud

Executive Summary Global food and water security analysis and management require precise and accurate global cropland-extent maps. Existing maps have limitations, in that they are (1) mapped using coarse-resolution remote-sensing data, resulting in the lack of precise mapping location of croplands and their accuracies; (2) derived by collecting and collating national statistical data that are often subjective, leading to substantial uncertainties in cropland-area estimates, as well as their locations; and (3) extracted from one or more classes of a land use–land cover product in which cropland classes are not the focus of mapping, leading to their mixing with other classes and creating significant errors of omission and commission. These limitations can be overcome by producing high-resolution cropland-extent maps using satellite-sensor data, such as Landsat 30-m resolution or higher. The most fundamental cropland product is the high-resolution cropland-extent map because all higher level cropland products, such as crop-watering method (that is, whether crops are irrigated or rainfed), crop types, cropping intensities, cropland fallows, crop productivity, and crop-water productivity, are dependent on a precise and accurate cropland-extent product. Given these realities, the overarching goal of this study was to produce a Landsat satellite-derived global cropland-extent product at 30-m resolution. The work, which involved a paradigm shift in how global cropland-extent maps are produced, involved the following five key steps: (1) petabyte-scale computing that involved multiyear, 8- to 16-day, time-series Landsat 30-m resolution data for the global land surface; (2) composition of analysis-ready data (ARD) cubes; (3) creation of a large global-reference data hub for machine learning; (4) use of multiple machine-learning algorithms (MLAs) by writing software and computing in the cloud; and (5) Google Earth Engine (GEE) cloud computing. The five key steps involved nine distinct phases. First, the world was segmented into 74 agroecological zones (AEZs). Second, Landsat 8- to 16-day data were used to time-composite 10-band (blue, green, red, near-infrared, short-wave infrared band 1, short-wave infrared band 2, thermal infrared, enhanced vegetation index, normalized difference water index, and normalized difference vegetation index) Landsat 30-m resolution data cubes for every 2- to 4-month time period during 3- to 4-year periods (stated as nominal-year 2015 or, simply, 2015), along with two additional 30-m resolution bands (Shuttle Radar Topography Mission elevation, and slope) in each of the 74 AEZs. Third, more than 100,000 reference-training data samples were collected using ground data (some of which were collected using a mobile application), as well as submeter- to 5-m-resolution, very high-resolution imagery sourced from other reliable sources. Fourth, reference-training data were used to create a knowledge base for separating cropland from noncropland. Fifth, MLAs such as the pixel-based supervised random forest and support-vector machines were written on the GEE using Python and JavaScript. Sixth, object-based recursive hierarchical segmentation algorithm was used, in addition to MLAs, to overcome uncertainties. Seventh, MLAs used the knowledge base to classify and separate cropland from noncropland. Eighth, accuracy assessment was conducted by generating error matrices for each of the 74 AEZs using 19,171 independent validation-data samples. Ninth, cropland areas were computed for all countries of the world and compared with United Nation’s (UN’s) Food and Agricultural Organization (FAO) and other national statistics. The outcome was a Landsat-derived global cropland-extent product at 30-m resolution (GCEP30), which has an overall accuracy of 91.7 percent. For the cropland class, producer’s accuracy was 83.4 percent, and user’s accuracy was 78.3 percent. GCEP30 calculated (using direct pixel count) the global net-cropland area (GNCA) for the year 2015 as 1.873 billion hectares (~12.6 percent of the Earth’s terrestrial area). The continental cropland distribution as a percentage of GNCA was Asia, 33 percent; Europe, 25.5 percent; Africa, 16.7 percent; North America, 14.4 percent; South America, 8.1 percent; and Australia and Oceania, 2.4 percent. The worldwide cropland areas in GCEP30 for 2015 were higher by 236 to 299 million hectares (Mha) compared to national statistics reported elsewhere for the same year (for example, in Food and Agriculture Organization’s corporate statistical database [FAOSTAT] and in the monthly irrigated and rainfed crop areas [MIRCA] database). The global cropland area reported for 2015 increased by 344 Mha (22.5 percent), compared to the year 2000. During the same period (2000–2015), the world’s population increased by 20 percent. Whereas some of these areal increases are real increases in cropland areas, others are due to the types of data, methods, and approaches used. Using the highest known resolution (compared to previous coarse-resolution global products) enabled this study to capture fragmented croplands. Coarse-resolution data compute areas on the basis of subpixels, which, for a large proportion of certain land use–land cover classes, will show only a certain percentage of the total pixel area as actual area. Subpixel areas can lead to substantial uncertainties in area computation, as determining the exact fraction of cropland areas within a coarse-resolution pixel is resource intensive and subject to errors. Other innovations in GCEP30 include reference-data hubs, machine learning, and cloud computing. Cropland areas in 214 countries, territories, departments, and regions were calculated for the year 2015 using GCEP30, on the basis of UN’s global administrative unit layers (GAUL) boundaries. The 10 leading countries in terms of cropland area (as a percentage of the GNCA) were India (9.6 percent), United States (8.95 percent), China (8.82 percent), Russia (8.32 percent), Brazil (3.42 percent), Ukraine (2.32 percent), Canada (2.29 percent), Argentina (2.05 percent), Indonesia (2 percent), and Nigeria (1.91 percent). Together, these 10 countries occupy 50 percent of the global cropland, and they have 52 percent of the global population. Their combined cropland area increased by 2 percent between 2000 and 2015, compared to the substantial increase in population of 517 million (15.5 percent). Together, India, United States, China, and Russia encompass 36 percent of the total area. In the United States and Canada, from 2000 to 2015, cropland decreased by about 2 percent, whereas their populations increased by 14 and 13 percent, respectively. The additional food requirements in these 10 countries, which are caused by increased populations, as well as increasing nutritional demands, are met by production increases in existing cropland or through virtual food trade, or both. More than 18 countries, territories, departments, or regions had 60 percent or more of their geographic area as cropland: Republic of Moldova, San Marino, and Hungary had more than 80 percent of the country’s area as cropland; Denmark, Ukraine, Ireland, and Bangladesh, 70 to 80 percent; and Uruguay, Netherlands, United Kingdom, Spain, Lithuania, Poland, Gaza Strip, Czechia, Italy, India, and Azerbaijan, 60 to 70 percent. Europe and South Asia can be considered agricultural capitals of the world, on the basis of their percentages of geographic area as cropland. United States, China, and Russia, which all have high cropland areas, are ranked second, third, and fourth in the world; India is ranked first. However, the amount of cropland as a percentage of the country’s geographic area is relatively very low for United States (18.3 percent), China (17.7 percent), and Russia (9.5 percent), whereas it is 60.5 percent for India. Most African and South American countries, territories, departments, or regions have less than 15 percent of their geographic area as cropland. China and India together house 36 percent of the world’s population; however, between 2000 and 2015, the amount of China’s cropland area fell by 18.9 percent, owing to urban expansion and the abandonment of farmlands caused by demographic changes (that is, the movement of population from villages to cities). In contrast, China’s population grew by 10 percent. The amount of India’s cropland increased by 8.5 percent, whereas its population grew by 20 percent. This study showed that, out of the 10 leading cropland countries, Ukraine, Nigeria, Russia, and Indonesia showed an 18 to 31 percent increase in cropland areas, on the basis of GCEP30 by the year 2015, compared to 2000. Nigeria’s cropland area increased by 25 percent, and its population increased by 31 percent in the same period. In these countries, food security is maintained by cropland expansion, productivity increases, and virtual food trade. Nevertheless, this trend of increasing net-cropland area and productivity will likely become difficult to maintain, owing to diminishing arable lands and plateauing of 50 years of continual yield increases, requiring policymakers to explore novel and data-supported approaches to solving future food security issues. The GCEP30 product, which can be browsed at full resolution at www.croplands.org , has been released for public download and use through U.S. Geological Survey (USGS)–National Aeronautics and Space Administration (NASA) Land Processes Distributed Active Archive Center (see https://lpdaac.usgs.gov/news/release-of-gfsad-30-meter-cropland-extent-products/ ).

Professional Paper↗
USGS2018

A database of natural monthly streamflow estimates from 1950 to 2015 for the conterminous United States

Quantifying and understanding the natural streamflow regime, defined as expected streamflow that would occur in the absence of anthropogenic modification to the hydrologic system, is critically important for the development of management strategies aimed at protecting aquatic ecosystems. Water balance models have been applied frequently to estimate natural flows, but are limited in the number of predictor variables that can be included. Here, a statistical machine learning technique — random forest modeling — was applied to estimate natural flows at a monthly time‐step from 1950 to 2015 for >2.5 million stream reaches in the conterminous United States (U.S.) using 200 potential predictor variables. We describe the development and documentation of this dataset and assess model performance. Model fit statistics (mean Nash–Sutcliffe efficiency = 0.85; observed/expected ratio = 0.94) indicate good correspondence between predicted and observed flows at nearly 2,000 streamgages. As an example application of the dataset, the observed streamflow record at a site prior to and after the construction of an upstream reservoir was compared with estimated natural flows to demonstrate the magnitude of seasonal depletions in streamflow due to the reservoir. This dataset can be applied to quantify natural and anthropogenic processes contributing to streamflow depletion or augmentation, and assess associated ecological effects.

Journal of the American Water Resources Associatio↗
USGS2019

Modelling gully-erosion susceptibility in a semi-arid region, Iran: Investigation of applicability of certainty factor and maximum entropy models

Gully erosion susceptibility mapping is a fundamental tool for land-use planning aimed at mitigating land degradation. However, the capabilities of some state-of-the-art data-mining models for developing accurate maps of gully erosion susceptibility have not yet been fully investigated. This study assessed and compared the performance of two different types of data-mining models for accurately mapping gully erosion susceptibility at a regional scale in Chavar, Ilam, Iran. The two methods evaluated were: Certainty Factor (CF), a bivariate statistical model; and Maximum Entropy (ME), an advanced machine learning model. Several geographic and environmental factors that can contribute to gully erosion were considered as predictor variables of gully erosion susceptibility. Based on an existing differential GPS survey inventory of gully erosion, a total of 63 eroded gullies were spatially randomly split in a 70:30 ratio for use in model calibration and validation, respectively. Accuracy assessments completed with the receiver operating characteristic curve method showed that the ME-based regional gully susceptibility map has an area under the curve (AUC) value of 88.6% whereas the CF-based map has an AUC of 81.8%. According to jackknife tests that were used to investigate the relative importance of predictor variables, aspect, distance to river, lithology and land use are the most influential factors for the spatial distribution of gully erosion susceptibility in this region of Iran. The gully erosion susceptibility maps produced in this study could be useful tools for land managers and engineers tasked with road development, urbanization and other future development.

Science of the Total Environment↗
USGS2025

Analyzing multi-year nitrate concentration evolution in Alabama aquatic systems using a machine learning model

Rising nitrate contamination in water systems poses significant risks to public health and ecosystem stability, necessitating advanced modeling to understand nitrate dynamics more accurately. This study applies the long short-term memory (LSTM) modeling to investigate the hydrologic and environmental factors influencing nitrate concentration dynamics in rivers and aquifers across the state of Alabama in the southeast of the United States. By integrating dynamic data such as streamflow and groundwater levels with static catchment attributes, the machine learning model identifies primary drivers of nitrate fluctuations, offering detailed insights into the complex interactions affecting multi-year nitrate concentrations in natural aquatic systems. In addition, a novel LSTM-based approach utilizes synthetic surface water nitrate data to predict groundwater nitrate levels, helping to address monitoring gaps in aquifers connected to these rivers. This method reveals potential correlations between surface water and groundwater nitrate dynamics, which is particularly meaningful given the lack of water quality observations in many aquifers. Field applications further show that, while the LSTM model effectively captures seasonal trends, limitations in representing extreme nitrate events suggest areas for further refinement. These findings contribute to data-driven water quality management, enhancing understanding of nitrate behavior in interconnected water systems.

Alabama↗
USGS2026

Classifying and mapping wetland vegetation assemblages in coastal Louisiana with Landsat imagery, 1985-2025

The tidal wetlands of Louisiana comprise about 25% of those found throughout the conterminous United States yet estimates of wetland loss rates in the region between 1932 and 2016 have exceeded 60 km 2 yr -1 . To mitigate further degradation and wetland loss in the region, a globally unprecedented $50B, 50-year plan for coastal Louisiana is driving restoration efforts, and demand exists from multiple stakeholders for regularly updated, regional-scale, accurate land cover information. We used machine learning (random forests; RF) and cloud computing to develop a new Landsat-based, marsh vegetation community geospatial dataset. The dataset depicts wetland vegetation community types defined in a previous study at annual (1985–2025) time steps at 30-m resolution. An RF algorithm was used to integrate training samples with feature variables derived from Landsat imagery, and the resulting geospatial data product achieved an overall correct classification rate of 78%. The approach for development of the land cover dataset presented here has potential for application in other coastal wetland habitats throughout the world.

BioRxiv↗
USGS2019

Digital mapping of ecological land units using a nationally scalable modeling framework

Ecological site descriptions (ESDs) and associated state-and-transition models (STMs) provide a nationally consistent classification and information system for defining ecological land units for management applications in the United States. Current spatial representations of ESDs, however, occur via soil mapping and are therefore confined to the spatial resolution used to map soils within a survey area. Land management decisions occur across a range of spatial scales and therefore require ecological information that spans similar scales. Digital mapping provides an approach for optimizing the spatial scale of modeling products to best serve decision makers and have the greatest impact in addressing land management concerns. Here, we present a spatial modeling framework for mapping ecological sites using machine learning algorithms, soil survey field observations, soil survey geographic databases, ecological site data, and a suite of remote sensing-based spatial covariates (e.g., hyper-temporal remote sensing, terrain attributes, climate data, land-cover, lithology). Based on the theoretical association between ecological sites and landscape biophysical properties, we hypothesized that the spatial distribution of ecological sites could be predicted using readily available geospatial data. This modeling approach was tested at two study areas within the western United States, representing 6.1 million ha on the Colorado Plateau and 7.5 million ha within the Chihuahuan Desert. Results show our approach was effective in mapping grouped ecological site classes (ESGs), with 10-fold cross-validation accuracies of 70% in the Colorado Plateau based on 1405 point observations across eight expertly-defined ESG classes and 79% in the Chihuahuan Desert based on 2589 point observations across nine expertly-defined ESG classes. Model accuracies were also evaluated using external-validation datasets; resulting in 56 and 44% correct classification for the Colorado Plateau and Chihuahuan Desert, respectively. National coverage of the training and covariate data used in this study provides opportunities for a consistent national-scale mapping effort of ecological sites.

Chihuahuan Desert, Colorado Plateau↗
USGS2025

Total uncertainty quantification in inverse solutions with deep learning surrogate models

We propose an approximate Bayesian method for quantifying the total uncertainty in inverse partial differential equation (PDE) solutions obtained with machine learning surrogate models, including operator learning models. The proposed method accounts for uncertainty in the observations, PDE, and surrogate models. First, we use the surrogate model to formulate a minimization problem in the reduced space for the maximum a posteriori (MAP) inverse solution. Then, we randomize the MAP objective function and obtain samples of the posterior distribution by minimizing different realizations of the objective function. We test the proposed framework by comparing it with the iterative ensemble smoother and deep ensembling methods for a nonlinear diffusion equation with an unknown space-dependent diffusion coefficient. Among other applications, this equation describes the flow of groundwater in an unconfined aquifer. Depending on the training dataset and ensemble sizes, the proposed method provides similar or more descriptive posteriors of the parameters and states than the iterative ensemble smoother method. Deep ensembling underestimates uncertainty and provides less-informative posteriors than the other two methods. Our results show that, despite inherent uncertainty, surrogate models can be used for parameter and state estimation as an alternative to the inverse methods relying on (more accurate) numerical PDE solvers.

Journal of Computational Physics↗
USGS2024

Predictive understanding of stream salinization in a developed watershed using machine learning

Stream salinization is a global issue, yet few models can provide reliable salinity estimates for unmonitored locations at the time scales required for ecological exposure assessments. Machine learning approaches are presented that use spatially limited high-frequency monitoring and spatially distributed discrete samples to estimate the daily stream-specific conductance across a watershed. We compare the predictive performance of space- and time-unaware Random Forest models and space- and time-aware Recurrent Graph Convolution Neural Network models (KGE: 0.67 and 0.64, respectively) and use explainable artificial intelligence methods to interpret model predictions and understand salinization drivers. These models are applied to the Delaware River Basin, a developed watershed with diverse land uses that experiences anthropogenic salinization from winter deicer applications. These models capture seasonality for the winter first flush of deicers, and the streams with elevated predictions correspond well with indicators of deicer application. This result suggests that these models can be used to identify potential salinity-impaired streams for winter best management practices. Daily salinity predictions are driven primarily by land cover (urbanization) trends that may represent anthropogenic salinization processes and weather at time scales up to three months. Such modeling approaches are likely transferable to other watersheds and can be applied to further understand salinization risks and drivers.

Environmental Science and Technology↗
USGS2025

Using imaging spectroscopy and elevation in machine learning to estimate soil salinity in intermittently tidal wetlands

Coastal soil salinization patterns are changing due to drought, sea level rise (SLR), and changing freshwater inflow. These changes are expected to impact coastal wetland plant health and ecosystem function, such as changes to biomass and productivity. These impacts have led to greater interest in how we monitor soil salinization across spatial and temporal scales. Remote sensing is a promising tool for estimating soil salinity at the spatial scales required for decision making by land managers. However, the development of a remote sensing estimation approach for wetland soil salinity must account for two factors: (1) the high spatial and temporal heterogeneity of coastal wetlands and (2) the fact that soil salinity is the result of multiple historical land use, hydrological, and geomorphic processes. In spring 2022, a combined airborne-field campaign, known as SHIFT, collected a weekly time series of airborne visible to shortwave infrared (VSWIR) image spectroscopy data. This dataset provides a unique opportunity to assess the application of fine spatial (5 m) and temporal (weekly) resolution VSWIR data to estimate root zone soil salinity; when combined with environmental variables such as elevation, these data can account for some of these factors. In this study, we utilized VSWIR and elevation datasets in a random forest regression to predict and map soil salinity in an intermittently tidal estuary, Devereux Slough, located in Santa Barbara County, California. The final model combined spectral indices with elevation to better capture soil salinity dynamics despite lower correlation ( r = 0.85) than solely using elevation ( r = 0.92). This research demonstrates the utility of remote sensing datasets, namely, elevation and the modified Anthocyanin Reflectance Index (mARI), for predicting root zone soil salinity in intermittently tidal coastal wetlands. These findings are an important step in advancing coastal remote sensing by creating a gridded salinity dataset that can be used for salinity monitoring and other coastal applications, such as modeling change in vegetation communities or ecosystems facing the impacts of climatic variability and change.

California↗
USGS2023

Application of a lipopolysaccharide (LPS)-stimulated mitogenesis assay in smallmouth bass (Micropterus dolomieu) to augment wild fish health studies

The utility of a functional immune assay for smallmouth bass ( Micropterus dolomieu ) lymphocyte mitogenesis was evaluated. Wild populations in the Potomac River have faced disease and mortality with immunosuppression from exposure to chemical contaminants a suspected component. However, a validated set of immune parameters to screen for immunosuppression in wild fish populations is not available. Prior to use in ecotoxicology studies, ancillary factors influencing the mitogenic response need to be understood. The assay was field-tested with fish collected from three sites in West Virginia as part of health assessments occurring in spring (pre-spawn; April–May) and fall (recrudescence; October–November). Anterior kidney leukocytes were exposed to lipopolysaccharide (LPS) from E.coli O111:B4 or mitogen-free media and proliferation was measured using imaging flow cytometry with advanced machine learning to distinguish lymphocytes. An anti-smallmouth bass IgM monoclonal antibody was used to identify IgM+ lymphocytes. Lymphocyte mitogenesis, or proliferative responses, varied by site and season and positively and negatively correlated with factors such as sex, age, tissue parasites, and macrophage aggregates. Background proliferation of IgM− lymphocytes was negatively correlated to LPS-induced proliferation in both seasons at all sites, but only in spring for IgM+ lymphocytes. The results demonstrate that many factors, in addition to chemical contaminants, may influence lymphocyte proliferation.

Virginia, West Virginia↗
USGS2023

Approaching the upper boundary of driver-response relationships: Identifying factors using a novel framework integrating quantile regression with interpretable machine learning

The identification of factors that may be forcing ecological observations to approach the upper boundary provides insight into potential mechanisms affecting driver-response relationships, and can help inform ecosystem management, but has rarely been explored. In this study, we propose a novel framework integrating quantile regression with interpretable machine learning. In the first stage of the framework, we estimate the upper boundary of a driver-response relationship using quantile regression. Next, we calculate “potentials” of the response variable depending on the driver, which are defined as vertical distances from the estimated upper boundary of the relationship to observations in the driver-response variable scatter plot. Finally, we identify key factors impacting the potential using a machine learning model. We illustrate the necessary steps to implement the framework using the total phosphorus (TP)-Chlorophyll a (CHL) relationship in lakes across the continental US. We found that the nitrogen to phosphorus ratio (N:P), annual average precipitation, total nitrogen (TN), and summer average air temperature were key factors impacting the potential of CHL depending on TP. We further revealed important implications of our findings for lake eutrophication management. The important role of N:P and TN on the potential highlights the co-limitation of phosphorus and nitrogen and indicates the need for dual nutrient criteria. Future wetter and/or warmer climate scenarios can decrease the potential which may reduce the efficacy of lake eutrophication management. The novel framework advances the application of quantile regression to identify factors driving observations to approach the upper boundary of driver-response relationships.

Frontiers of Environmental Science & Engineering↗
USGS2024

Chlorophyll a in lakes and streams of the United States (2005–2022)

The concentration of chlorophyll a in phytoplankton and periphyton represents the amount of algal biomass. We compiled an 18-year record (2005–2022) of pigment data from water bodies across the United States (US) to support efforts to develop process-based, machine learning, and remote sensing models for prediction of harmful algal blooms (HABs). To our knowledge, this dataset of nearly 84,000 sites and over 1,374,000 pigment measurements is the largest compilation of harmonized discrete, laboratory-extracted chlorophyll data for the US. These data were compiled from the Water Quality Portal (WQP) and previously unpublished U.S. Geological Survey’s National Water Quality Laboratory (NWQL) data. Data were harmonized for reporting units, pigment type, duplicate values, collection depth, site name, negative values, and some extreme values. Across the country, data show great variation by state in sampling frequency, distribution, and methods. Uses for such data include the calibration of models, calibration of field sensors, examination of relationship to nutrients and other drivers, evaluation of temporal trends, and other applications addressing local to national scale concerns.

Nature Scientific Data↗
USGS2022

Using landscape genomics to delineate future adaptive potential for climate change in the Yosemite toad (Anaxyrus canorus)

An essential goal in conservation biology is delineating population units that maximize the probability of species persisting into the future and adapting to future environmental change. However, future-facing conservation concerns are often addressed using retrospective patterns that could be irrelevant. We recommend a novel landscape genomics framework for delineating future “Geminate Evolutionary Units” (GEUs) in a focal species: (1) identify loci under environmental selection, (2) model and map adaptive conservation units that may spawn future lineages, (3) forecast relative selection pressures on each future lineage, and (4) estimate their fitness and likelihood of persistence using geo-genomic simulations. Using this process, we delineated conservation units for the Yosemite toad ( Anaxyrus canorus ), a U.S. federally threatened species that is highly vulnerable to climate change. We used a genome-wide dataset, redundancy analysis, and Bayesian association methods to identify 24 candidate loci responding to climatic selection ( R 2 ranging from 0.09 to 0.52), after controlling for demographic structure. Candidate loci included genes such as MAP3K5, involved in cellular response to environmental change. We then forecasted future genomic response to climate change using the multivariate machine learning algorithm Gradient Forests. Based on all available evidence, we found three GEUs in Yosemite National Park, reflecting contrasting adaptive optima: YF-North (high winter snowpack with moderate summer rainfall), YF-East (low to moderate snowpack with high summer rainfall), and YF-Low-Elevation (low snowpack and rainfall). Simulations under the RCP 8.5 climate change scenario suggest that the species will decline by 29% over 90 years, but the highly diverse YF-East lineage will be least impacted for two reasons: (1) geographically it will be sheltered from the largest climatic selection pressures, and (2) its standing genetic diversity will promote a faster adaptive response. Our approach provides a comprehensive strategy for protecting imperiled non-model species with genomic data alone and has wide applicability to other declining species.

California↗
USGS2010

New software methods in radar ornithology using WSR-88D weather data and potential application to monitoring effects of climate change on bird migration

Radar ornithology has provided tools for studying the movement of birds, especially related to migration. Researchers have presented qualitative evidence suggesting that birds, or at least migration events, can be identified using large broad scale radars such as the WSR-88D used in the NEXRAD weather surveillance system. This is potentially a boon for ornithologists because such data cover a large portion of the United States, are constantly being produced, are freely available, and have been archived since the early 1990s. A major obstacle to this research, however, has been that identifying birds in NEXRAD data has required a trained technician to manually inspect a graphically rendered radar sweep. A single site completes one volume scan every five to ten minutes, producing over 52,000 volume scans in one year. This is an immense amount of data, and manual classification is infeasible. We have developed a system that identifies biological echoes using machine learning techniques. This approach begins with training data using scans that have been classified by experts, or uses bird data collected in the field. The data are preprocessed to ensure quality and to emphasize relevant features. A classifier is then trained using this data and cross validation is used to measure performance. We compared neural networks, naive Bayes, and k-nearest neighbor classifiers. Empirical evidence is provided showing that this system can achieve classification accuracies in the 80th to 90th percentile. We propose to apply these methods to studying bird migration phenology and how it is affected by climate variability and change over multiple temporal scales.

Conference Paper↗
USGS2022

Integration of vegetation classification with land cover mapping: Lessons from regional mapping efforts in the Americas

Aims : Natural resource management and biodiversity conservation rely on inventories of vegetation that span multiple management or political jurisdictions. However, while remote sensing data and analytical tools have enabled production of maps at increasing spatial resolution and reliability, there are limited examples where national or continental-scaled maps are produced to represent vegetation at high thematic detail. We illustrate two examples that have bridged the gap between traditional land cover mapping and modern vegetation classification. Study area : Our two case studies include national ( USA ) and continental (North and South America) vegetation and land cover mapping. These studies span conditions from subpolar to tropical latitudes of the Americas. Methods : Both case studies used a supervised modeling approach with the International Vegetation Classification ( IVC ) to produce maps that provide for greater thematic detail. Georeferenced locations for these vegetation types are used by machine learning algorithms to train a predictive model and generate a distribution map. Results : The USA LANDFIRE (Landscape Fire and Resource Management Planning Tools Project) case study illustrates how a history of vegetation-based classification and availability of key inputs can come together to generate standard map products covering more than 9.8 million km 2 that are unsurpassed anywhere in the world in terms of spatial and thematic resolution. That being said, it also remains clear that mapping at the thematic resolution of the IVC Group and finer resolution require very large and spatially balanced inputs of georeferenced samples. Even with extensive prior data collection efforts, these remain a key limitation. The NatureServe effort for the Americas - encompassing 22% of the global land surface - demonstrates methods and outputs suitable for worldwide application at continental scales. Conclusions : Continued collection of input data used in the case studies could enable mapping at these spatial and thematic resolutions around the globe.

Vegetation Classification and Survey↗
USGS2022

Impact of spectral resolution on quantifying cyanobacteria in lakes and reservoirs: A machine-learning assessment

Cyanobacterial harmful algal blooms are an increasing threat to coastal and inland waters. These blooms can be detected using optical radiometers due to the presence of phycocyanin (PC) pigments. The spectral resolution of best-available multispectral sensors limits their ability to diagnostically detect PC in the presence of other photosynthetic pigments. To assess the role of spectral resolution in the determination of PC, a large (N = 905) database of colocated in situ radiometric spectra and PC are employed. We first examine the performance of selected widely used machine-learning (ML) models against that of benchmark algorithms for hyperspectral remote sensing reflectance ( R r s ) spectra resampled to the spectral configuration of the Hyperspectral Imager for the Coastal Ocean (HICO) with a full-width at half-maximum (FWHM) of < 6 nm. Results show that the multilayer perceptron (MLP) neural network applied to HICO spectral configurations (median errors < 65%) outperforms other ML models. This model is subsequently applied to R r s spectra resampled to the band configuration of existing satellite instruments and of the one proposed for the next Landsat sensor. These results confirm that employing MLP models to estimate PC from hyperspectral data delivers tangible improvements compared with retrievals from multispectral data and benchmark algorithms (with median errors between ~73% and 126%) and shows promise for developing a globally applicable cyanobacteria measurement approach.

IEEE Transactions in Geoscience and Remote Sensing↗
USGS2022

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

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

Journal of Applied Ecology↗