SEARCH · USGS Science

Results for “Soil Research”

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 631 records · Page 35Linked to original sources

USGS2012

Persistence of soil organic matter in eroding versus depositional landform positions

Soil organic matter (SOM) processes in dynamic landscapes are strongly influenced by soil erosion and sedimentation. We determined the contribution of physical isolation of organic matter (OM) inside aggregates, chemical interaction of OM with soil minerals, and molecular structure of SOM in controlling storage and persistence of SOM in different types of eroding and depositional landform positions. By combining density fractionation with elemental and spectroscopic analyses, we showed that SOM in depositional settings is less transformed and better preserved than SOM in eroding landform positions. However, which environmental factors exert primary control on storage and persistence of SOM depended on the nature of the landform position considered. In an annual grassland watershed, protection of SOM by physical isolation inside aggregates and chemical association of organic matter (complexation) with soil minerals, as assessed by correlation with radiocarbon concentration, were more effective in the poorly drained, lowest-lying depositional landform positions, compared to well-drained landform positions in the upper parts of the watershed. Results of this study demonstrated that processes of soil erosion and deposition are important mechanisms of long-term OM stabilization.

Journal of Geophysical Research: Biogeosciences↗

USGS2010

Flow path studies in forested watersheds of headwater tributaries of Brush Brook, Vermont

An investigation was undertaken into how headwater tributaries of Brush Brook, Vermont, could have average p H differences of almost two units (4.75 and 6.7). Sampling along four tributaries revealed that most of one tributary, below an area of seeps, had consistently higher p H, Ca 2+ , Mg 2+ , and K + , and lower Al than other sites. Bedrock mapping showed numerous fractures in vicinity of the seeps. A portion of this tributary's watershed and a portion of an acid tributary's watershed were intensively mapped for soil depth. Sampling showed the widespread existence of dense basal till in the watershed of the acid tributary but none in that of the near-neutral stream. Lateral flow, found above the dense till, was chemically similar to that of the acid tributary and to solutions sampled from soil B horizons. There were no differences in the average p H of nonseep soils sampled from either watershed. Flow paths are hypothesized to be through the B horizons in the acid tributaries and from below the soil profile in the near-neutral tributary. The acid catchment should be more sensitive to environmental change.

Vermont↗

USGS1997

Large-scale variability of wind erosion mass flux rates at Owens Lake 1. Vertical profiles of horizontal mass fluxes of wind-eroded particles with diameter greater than 50 μm

A field experiment at Owens (dry) Lake, California, tested whether and how the relative profiles of airborne horizontal mass fluxes for >50-μm wind-eroded particles changed with friction velocity. The horizontal mass flux at almost all measured heights increased proportionally to the cube of friction velocity above an apparent threshold friction velocity for all sediment tested and increased with height except at one coarse-sand site where the relative horizontal mass flux profile did not change with friction velocity. Size distributions for long-time-averaged horizontal mass flux samples showed a saltation layer from the surface to a height between 30 and 50 cm, above which suspended particles dominate. Measurements from a large dust source area on a line parallel to the wind showed that even though the saltation flux reached equilibrium &sim;650 m downwind of the starting point of erosion, weakly suspended particles were still input into the atmosphere 1567 m downwind of the starting point; thus the saltating fraction of the total mass flux decreased after 650 m. The scale length difference and ratio of 70/30 suspended mass flux to saltation mass flux at the farthest down wind sampling site confirm that suspended particles are very important for mass budgets in large source areas and that saltation mass flux can be a variable fraction of total horizontal mass flux for soils with a substantial fraction of <100-μm particles.

Journal of Geophysical Research D: Atmospheres↗

USGS2006

Using Modified Mercalli Intensities to estimate acceleration response spectra for the 1906 San Francisco earthquake

We derive and test relations between the Modified Mercalli Intensity (MMI) and the pseudo-acceleration response spectra at 1.0 and 0.3 s - SA(1.0 s) and SA(0.3 s) - in order to map response spectral ordinates for the 1906 San Francisco earthquake. Recent analyses of intensity have shown that MMI ??? 6 correlates both with peak ground velocity and with response spectra for periods from 0.5 to 3.0 s. We use these recent results to derive a linear relation between MMI and log SA(1.0 s), and we refine this relation by comparing the SA(1.0 s) estimated from Boatwright and Bundock's (2005) MMI map for the 1906 earthquake to the SA(1.0 s) calculated from recordings of the 1989 Loma Prieta earthquake. South of San Jose, the intensity distributions for the 1906 and 1989 earthquakes are remarkably similar, despite the difference in magnitude and rupture extent between the two events. We use recent strong motion regressions to derive a relation between SA(1.0 s) and SA(0.3 s) for a M7.8 strike-slip earthquake that depends on soil type, acceleration level, and source distance. We test this relation by comparing SA(0.3 s) estimated for the 1906 earthquake to SA(0.3 s) calculated from recordings of both the 1989 Loma Prieta and 1994 Northridge earthquakes, as functions of distance from the fault. ?? 2006, Earthquake Engineering Research Institute.

Earthquake Spectra↗

USGS2008

Alaska Melilotus invasions: Distribution, origin, and susceptibility of plant communities

Melilotus alba and M. officinalis were introduced to Alaska in 1913 as potential forage crops. These species have become naturalized and are now invading large, exotic plant-free regions of Alaska. We determined distributions of M. alba and M. officinalis in Alaska from surveys conducted each summer from 2002 to 2005. Melilotus alba and M. officinalis occurred at 721 and 205 sites, respectively (39,756 total sites surveyed). The northward limit for M. alba and M. officinalis was 67.15??N and 64.87??N, respectively. Both species were strictly associated with soil disturbance. Melilotus alba extended no farther than 15 m from road edges except where M. alba on roadsides met river floodplains and dispersed downriver (Matanuska and Nenana Rivers). Melilotus has now reached the Tanana River, a tributary of the Yukon River. Populations on floodplains were most extensive on braided sections. On the Nenana River, soil characteristics did not differ between where M. alba was growing versus similar areas where it had not yet reached. The pH of river soils (7.9-8.3) was higher than highway soils (7.3). Upland taiga plant communities grow on acid soils which may protect them from invasion by Melilotus, which prefer alkaline soils; however, early succession communities on river floodplains are susceptible because soils are alkaline. ?? 2008 Regents of the University of Colorado.

Arctic, Antarctic, and Alpine Research↗

USGS2008

Selected Water- and Sediment-Quality, Aquatic Biology, and Mine-Waste Data from the Ely Copper Mine Superfund Site, Vershire, VT, 1998-2007

The data contained in this report are a compilation of selected water- and sediment-quality, aquatic biology, and mine-waste data collected at the Ely Copper Mine Superfund site in Vershire, VT, from August 1998 through May 2007. The Ely Copper Mine Superfund site is in eastern, central Vermont (fig. 1) within the Vermont Copper Belt (Hammarstrom and others, 2001). The Ely Copper Mine site was placed on the U.S. Environmental Protection Agency (USEPA) National Priorities List in 2001. Previous investigations conducted at the site documented that the mine is contributing metals and highly acidic waters to local streams (Hammarstrom and others, 2001; Holmes and others, 2002; Piatak and others, 2003, 2004, and 2006). The U.S. Geological Survey (USGS), in cooperation with the USEPA, compiled selected data from previous investigations into uniform datasets that will be used to help characterize the extent of contamination at the mine. The data may be used to determine the magnitude of biological impacts from the contamination and in the development of remediation activities. This report contains analytical data for samples collected from 98 stream locations, 6 pond locations, 21 surface-water seeps, and 29 mine-waste locations. The 98 stream locations are within 3 streams and their tributaries. Ely Brook flows directly through the Ely Copper Mine then into Schoolhouse Brook (fig. 2), which joins the Ompompanoosuc River (fig. 1). The six pond locations are along Ely Brook Tributary 2 (fig. 2). The surface-water seeps and mine-waste locations are near the headwaters of Ely Brook (fig. 2 and fig. 3). The datasets 'Site_Directory' and 'Coordinates' contain specific information about each of the sample locations including stream name, number of meters from the mouth of stream, geographic coordinates, types of samples collected (matrix of sample), and the figure on which the sample location is depicted. Data have been collected at the Ely Copper Mine Superfund site by the USEPA, the Vermont Department of Environmental Conservation (VTDEC), and the USGS. Data also have been collected on behalf of USEPA by the following agencies: Arthur D. Little Incorporated (ADL), U.S. Army Cold Region Research and Engineering Laboratory (CRREL), URS Corporation (URS), USEPA, and USGS. These data provide information about the aquatic communities and their habitats, including chemical analyses of surface water, pore water, sediments, and fish tissue; assessments of macroinvertebrate and fish assemblages; physical characteristics of sediments; and chemical analyses of soil and soil leachate collected in and around the piles of mine waste.

Data Series↗

USGS1999

Explaining spatial variability in mean annual runoff in the conterminous United States

The hydrologic concepts needed in a water-balance model to estimate the spatial variation in mean annual runoff for the 344 climate divisions in the conterminous United States (U.S.) were determined. The concepts that were evaluated were the climatic supply of water (precipitation), climatic demand for water (potential evapotranspiration), seasonality in supply and demand, and soil-moisture-storage capacity. Most (91%) of the spatial variability in mean annual runoff for the climate divisions in the conterminous U.S. was explained by the spatial variability of mean annual precipitation minus mean annual potential evapotranspiration. When soil-moisture-storage capacity and seasonality in supply and demand were added to the water balance, the explained variance in mean annual runoff increased slightly, and the error in estimated mean annual runoff decreased significantly. Adding soil-moisture-storage capacity and seasonality in supply and demand provided the most improvement in areas where seasonal supply and demand are out of phase.

conterminous United States↗

USGS2005

A hierarchical perspective of plant diversity

Predictive models of plant diversity have typically focused on either a landscapea??s capacity for richness (equilibrium models), or on the processes that regulate competitive exclusion, and thus allow species to coexist (nonequilibrium models). Here, we review the concepts and purposes of a hierarchical, multiscale model of the controls of plant diversity that incorporates the equilibrium model of climatic favorability at macroscales, nonequilibrium models of competition at microscales, and a mixed model emphasizing environmental heterogeneity at mesoscales. We evaluate the conceptual model using published data from three spatially nested datasets: (1) a macroscale analysis of ecoregions in the continental and western U.S.; (2) a mesoscale study in California; and (3) a microscale study in the Siskiyou Mountains of Oregon and California. At the macroscale (areas from 3889 km2 to 638,300 km2), climate (actual evaporation) was a strong predictor of tree diversity (R2 = 0.80), as predicted by the conceptual model, but area was a better predictor for vascular plant diversity overall (R2 = 0.38), which suggests different types of plants differ in their sensitivity to climatic controls. At mesoscales (areas from 1111 km2 to 15,833 km2 ), climate was still an important predictor of richness (R2 = 0.52), but, as expected, topographic heterogeneity explained an important share of the variance (R2 = 0.19), showed positive correlations with diversity of trees, shrubs, and annual and perennial herbs, and was the primary predictor of shrub and annual plant species richness. At microscales (0.1 ha plots), spatial patterns of diversity showed a clear unimodal pattern along a climatea??driven productivity gradient and a negative relationship with soil fertility. The strong decline in understory and total diversity at the most productive sites suggests that competitive controls, as predicted, can override climatic controls at this scale. We conclude that this hierarchical, multiscale model provides a sound basis to understand and analyze plant species diversity. Specifically, future research should employ the principles in this paper to explore climatic controls on species richness of different life forms, better quantify environmental heterogeneity in landscapes, and analyze how these largea??scale factors interact with local nonequilibrium dynamics to maintain plant diversity.

The Quarterly Review of Biology↗

USGS1984

LANDSAT M. S. S. IMAGE MOSAIC OF TUNISIA.

The Landsat mosaic of Tunisia funded by USAID for the Remote Sensing Laboratory, Soils Division, Ministry of Agriculture, Tunisia, was completed by the USGS in September 1983. It is a mixed mosaic associating digital corrections and enhancements to manual mosaicking and corresponding to the Tunisian request for high resolution and the limited available funds. The scenes were processed by the Environmental Research Institute of Michigan, resampling the data geodesically corrected to fit the Universal Transverse Mercator projection using control points from topographic maps at 1:50,000 and 1:100,000 scales available in the U. S. The mosaicking was done in the Eastern Mapping Center under the supervision of the Graphic Arts System Section. The three black and white mosaics were made at the 1:1,000,000 scale and various products generated. They included color film positives at 1:2,000,000 and 1:4,000,000 scales reproducible in the Remote Sensing Laboratory in Tunis and corresponding color prints as well as tricolor prints at various scales from 1:500,000 to 1:2,000,000.

Conference Paper↗

USGS2021

A U.S.-China EcoPartnership study of disturbed wetland vegetation in West Dongting Lake, China

West Dongting Lake in China is important for human livelihoods and habitat of migratory waterfowl and other wildlife. The waterway re-engineering and agriculture intensification have contributed to changes in hydrology, sediment, and vegetation on the floodplain. This paper describes an EcoPartnership program conducted by the U.S. Geological Survey, Wetland and Aquatic Research Center, and Beijing Forestry University. It focused on the development of a wetland ecosystem network in West Dongting Lake with technical support from the U.S. partner using a number of related studies to examine wetland vegetation dynamics from upstream to downstream along the tributaries. The results of U.S. studies showed that the regeneration potential of species might be altered by changes in climate and local environment, and seed bank depletion by germination may be a major conservation threat in a future with recurring droughts in swamps of the southeastern United States. In the monsoonal wetlands of West Dongting Lake, the soil seed bank could be used as a seed source for revegetation after hydrologic restoration with the introduction of certain foundational species and the removal of poplar plantations. Also, West Dongting Lake is at high ecological risk of mercury pollution. Wetland ecosystem monitoring may allow managers to use the information to predict effects of climate change, water level and flow changes on sedimentation, and to manage for desired vegetation to support waterfowl and ecosystem services. The cooperation of two countries through the EcoPartnership program is now well established and poised for extensive research projects in the future.

West Dongting Lake↗

USGS2011

An analysis of modern pollen rain from the Maya lowlands of northern Belize

In the lowland Maya area, pollen records provide important insights into the impact of past human populations and climate change on tropical ecosystems. Despite a long history of regional paleoecological research, few studies have characterized the palynological signatures of lowland ecosystems, a fact which lowers confidence in ecological inferences made from palynological data. We sought to verify whether we could use pollen spectra to reliably distinguish modern ecosystem types in the Maya lowlands of Central America. We collected 23 soil and sediment samples from eight ecosystem types, including upland, riparian, secondary, and swamp ( bajo ) forests; pine savanna; and three distinct wetland communities. We analyzed pollen spectra with non-metric multidimensional scaling (NMDS), and found significant compositional differences in ecosystem types' pollen spectra. Forested sites had spectra dominated by Moraceae/Urticaceae pollen, while non-forested sites had significant portions of Poaceae, Asteraceae, and Amaranthaceae pollen. Upland, bajo , and riparian forest differed in representation of Cyperaceae, Bactris -type, and Combretaceae/Melastomataceae pollen. High percentages of pine ( Pinus ), oak ( Quercus ), and the presence of Byrsonima characterized pine savanna. Despite its limited sample size, this study provides one of the first statistical analyses of modern pollen rain in the Maya lowlands. Our results show that pollen assemblages can accurately reflect differences between ecosystem types, which may help refine interpretations of pollen records from the Maya area.

Maya lowlands↗

USGS2018

Analysis of vegetation recovery surrounding a restored wetland using the normalized difference infrared index (NDII) and normalized difference vegetation index (NDVI)

Watershed restoration efforts seek to rejuvenate vegetation, biological diversity, and land productivity at Cienega San Bernardino, an important wetland in southeastern Arizona and northern Sonora, Mexico. Rock detention and earthen berm structures were built on the Cienega San Bernardino over the course of four decades, beginning in 1984 and continuing to the present. Previous research findings show that restoration supports and even increases vegetation health despite ongoing drought conditions in this arid watershed. However, the extent of restoration impacts is still unknown despite qualitative observations of improvement in surrounding vegetation amount and vigor. We analyzed spatial and temporal trends in vegetation greenness and soil moisture by applying the normalized difference vegetation index (NDVI) and normalized difference infrared index (NDII) to one dry summer season Landsat path/row from 1984 to 2016. The study area was divided into zones and spectral data for each zone was analyzed and compared with precipitation record using statistical measures including linear regression, Mann– Kendall test, and linear correlation. NDVI and NDII performed differently due to the presence of continued grazing and the effects of grazing on canopy cover; NDVI was better able to track changes in vegetation in areas without grazing while NDII was better at tracking changes in areas with continued grazing. Restoration impacts display higher greenness and vegetation water content levels, greater increases in greenness and water content through time, and a decoupling of vegetation greenness and water content from spring precipitation when compared to control sites in nearby tributary and upland areas. Our results confirm the potential of erosion control structures to affect areas up to 5 km downstream of restoration sites over time and to affect 1 km upstream of the sites.

Cuenca Los Ojos, San Bernardino National Wildlife ↗

USGS2023

Investigating permafrost carbon dynamics in Alaska with artificial intelligence

Positive feedbacks between permafrost degradation and the release of soil carbon into the atmosphere impact land–atmosphere interactions, disrupt the global carbon cycle, and accelerate climate change. The widespread distribution of thawing permafrost is causing a cascade of geophysical and biochemical disturbances with global impacts. Currently, few earth system models account for permafrost carbon feedback (PCF) mechanisms. This research study integrates artificial intelligence (AI) tools and information derived from field-scale surveys across the tundra and boreal landscapes in Alaska. We identify and interpret the permafrost carbon cycling links and feedback sensitivities with GeoCryoAI, a hybridized multimodal deep learning (DL) architecture of stacked convolutionally layered, memory-encoded recurrent neural networks (NN). This framework integrates in-situ measurements and flux tower observations for teacher forcing and model training. Preliminary experiments to quantify, validate, and forecast permafrost degradation and carbon efflux across Alaska demonstrate the fidelity of this data-driven architecture. More specifically, GeoCryoAI logs the ecological memory and effectively learns covariate dynamics while demonstrating an aptitude to simulate and forecast PCF dynamics—active layer thickness (ALT), carbon dioxide flux (CO 2 ), and methane flux (CH 4 )—with high precision and minimal loss (i.e. ALT RMSE : 1.327 cm [1969–2022]; CO 2 RMSE : 0.697 µ molCO 2 m −2 s −1 [2003–2021]; CH 4 RMSE : 0.715 nmolCH 4 m −2 s −1 [2011–2022]). ALT variability is a sensitive harbinger of change, a unique signal characterizing the PCF, and our model is the first characterization of these dynamics across space and time.

Alaska↗

USGS2022

Acetylenotrophic and diazotrophic Bradyrhizobium sp. strain I71 from TCE-contaminated soils

Abstract Acetylene (C 2 H 2 ) is a molecule rarely found in nature, with very few known natural sources, but acetylenotrophic microorganisms can use acetylene as their primary carbon and energy source. As of 2018 there were 15 known strains of aerobic and anaerobic acetylenotrophs; however, we hypothesize there may yet be unrecognized diversity of acetylenotrophs in nature. This study expands the known diversity of acetylenotrophs by isolating the aerobic acetylenotroph, Bradyrhizobium sp. strain I71, from trichloroethylene (TCE)-contaminated soils. Strain I71 is a member of the class Alphaproteobacteria and exhibits acetylenotrophic and diazotrophic activities, the only two enzymatic reactions known to transform acetylene. This unique capability in the isolated strain may increase the genus’ economic impact beyond agriculture as acetylenotrophy is closely linked to bioremediation of chlorinated contaminants. Computational analyses indicate that the Bradyrhizobium sp. strain I71 genome contains 522 unique genes compared to close relatives. Moreover, applying a novel hidden Markov model of known acetylene hydratase (AH) enzymes identified a putative AH enzyme. Protein annotation with I-TASSER software predicted the AH from the microbe Syntrophotalea acetylenica as the closest structural and functional analog. Furthermore, the putative AH was flanked by horizontal gene transfer (HGT) elements, like that of AH in anaerobic acetylenotrophs, suggesting an unknown source of acetylene or acetylenic substrate in the environment that is selecting for the presence of AH. Importance The isolation of Bradyrhizobium strain I71 expands the distribution of acetylene-consuming microbes to include a group of economically important microorganisms. Members of Bradyrhizobium are well studied for their abilities to improve plant health and increase crop yields by providing bioavailable nitrogen. Additionally, acetylene-consuming microbes have been shown to work in tandem with other microbes to degrade soil contaminants. Based on genome, cultivation, and protein prediction analysis, the ability to consume acetylene is likely not widespread within the genus Bradyrhizobium . These findings suggest that the suite of phenotypic capabilities of strain I71 may be unique and make it a good candidate for further study in several research avenues.

California↗

USGS2004

The new GFDL global atmosphere and land model AM2-LM2: Evaluation with prescribed SST simulations

The configuration and performance of a new global atmosphere and land model for climate research developed at the Geophysical Fluid Dynamics Laboratory (GFDL) are presented. The atmosphere model, known as AM2, includes a new gridpoint dynamical core, a prognostic cloud scheme, and a multispecies aerosol climatology, as well as components from previous models used at GFDL. The land model, known as LM2, includes soil sensible and latent heat storage, groundwater storage, and stomatal resistance. The performance of the coupled model AM2-LM2 is evaluated with a series of prescribed sea surface temperature (SST) simulations. Particular focus is given to the model's climatology and the characteristics of interannual variability related to El Nin??o-Southern Oscillation (ENSO). One AM2-LM2 integration was perfor med according to the prescriptions of the second Atmospheric Model Intercomparison Project (AMIP II) and data were submitted to the Program for Climate Model Diagnosis and Intercomparison (PCMDI). Particular strengths of AM2-LM2, as judged by comparison to other models participating in AMIP II, include its circulation and distributions of precipitation. Prominent problems of AM2-LM2 include a cold bias to surface and tropospheric temperatures, weak tropical cyclone activity, and weak tropical intraseasonal activity associated with the Madden-Julian oscillation. An ensemble of 10 AM2-LM 2 integrations with observed SSTs for the second half of the twentieth century permits a statistically reliable assessment of the model's response to ENSO. In general, AM2-LM2 produces a realistic simulation of the anomalies in tropical precipitation and extratropical circulation that are associated with ENSO. ?? 2004 American Meteorological Society.

Journal of Climate↗

USGS1982

Earthquake engineering research program in Chile

Earthquake engineering research in Chile has been carried out for more than 30 years. Systematic research is done at the university of Chile in Santiago. Other universities such as the Catholic University, university of Concepcion, and the Federico Santa Maria Technical University have begun to teach and conduct research in earthquake engineering in recent years. The fields of research at the University of Chile include the following: Theoretical and empirical characterization of strong round motion. Installation of strong-motion accelerographs. Capacity of strong ground motion to produce structural damage. Seismic risk. Seismic behavior of reinforced masonry houses and buildings. Seismic design of high-rise buildings. Soil dynamics with emphasis on soil liquefaction and soil stability.

Earthquake Information Bulletin (USGS)↗

USGS2017

Integrating landslide and liquefaction hazard and loss estimates with existing USGS real-time earthquake information products

The U.S. Geological Survey (USGS) has made significant progress toward the rapid estimation of shaking and shakingrelated losses through their Did You Feel It? (DYFI), ShakeMap, ShakeCast, and PAGER products. However, quantitative estimates of the extent and severity of secondary hazards (e.g., landsliding, liquefaction) are not currently included in scenarios and real-time post-earthquake products despite their significant contributions to hazard and losses for many events worldwide. We are currently running parallel global statistical models for landslides and liquefaction developed with our collaborators in testing mode, but much work remains in order to operationalize these systems. We are expanding our efforts in this area by not only improving the existing statistical models, but also by (1) exploring more sophisticated, physics-based models where feasible; (2) incorporating uncertainties; and (3) identifying and undertaking research and product development to provide useful landslide and liquefaction estimates and their uncertainties. Although our existing models use standard predictor variables that are accessible globally or regionally, including peak ground motions, topographic slope, and distance to water bodies, we continue to explore readily available proxies for rock and soil strength as well as other susceptibility terms. This work is based on the foundation of an expanding, openly available, case-history database we are compiling along with historical ShakeMaps for each event. The expected outcome of our efforts is a robust set of real-time secondary hazards products that meet the needs of a wide variety of earthquake information users. We describe the available datasets and models, developments currently underway, and anticipated products.

Conference Paper↗

USGS2014

Linking bovine tuberculosis on cattle farms to white-tailed deer and environmental variables using Bayesian hierarchical analysis

Bovine tuberculosis is a bacterial disease caused by Mycobacterium bovis in livestock and wildlife with hosts that include Eurasian badgers ( Meles meles ), brushtail possum ( Trichosurus vulpecula ), and white-tailed deer ( Odocoileus virginianus ). Risk-assessment efforts in Michigan have been initiated on farms to minimize interactions of cattle with wildlife hosts but research on M. bovis on cattle farms has not investigated the spatial context of disease epidemiology. To incorporate spatially explicit data, initial likelihood of infection probabilities for cattle farms tested for M. bovis , prevalence of M. bovis in white-tailed deer, deer density, and environmental variables for each farm were modeled in a Bayesian hierarchical framework. We used geo-referenced locations of 762 cattle farms that have been tested for M. bovis , white-tailed deer prevalence, and several environmental variables that may lead to long-term survival and viability of M. bovis on farms and surrounding habitats (i.e., soil type, habitat type). Bayesian hierarchical analyses identified deer prevalence and proportion of sandy soil within our sampling grid as the most supported model. Analysis of cattle farms tested for M. bovis identified that for every 1% increase in sandy soil resulted in an increase in odds of infection by 4%. Our analysis revealed that the influence of prevalence of M. bovis in white-tailed deer was still a concern even after considerable efforts to prevent cattle interactions with white-tailed deer through on-farm mitigation and reduction in the deer population. Cattle farms test positive for M. bovis annually in our study area suggesting that the potential for an environmental source either on farms or in the surrounding landscape may contributing to new or re-infections with M. bovis . Our research provides an initial assessment of potential environmental factors that could be incorporated into additional modeling efforts as more knowledge of deer herd factors and cattle farm prevalence is documented.

Michigan↗