Geology Reports⌕ Search

Geology topics

Research about Wisconsin

Source-linked reports with geographic coverage including Wisconsin.

At least 37 records · Page 2Linked to original sources
USGS2025

Geomorphic habitat response units for urban stream rehabilitation, Milwaukee, Wisconsin

Urban stream rehabilitation plans can benefit from knowledge of the landscape setting and vegetative communities that were adjacent to streams prior to urbanization. Downstream to upstream connections of these characteristics can be relevant for native migratory fish species that have a range of preferred spawning habitats. Based on a need for more quantitative data on these potential connections, the U.S. Geological Survey assembled geomorphic characteristics, surficial geology, and pre-Euro-American settlement vegetation for 333 kilometers of stream segments in the Kinnickinnic River and Menomonee River subbasins of the Milwaukee River, Wisconsin. Channel slopes ranged from less than 0.3 percent to greater than 2 percent, covering at least two channel morphology and bedform types spanning low-energy irregular and pool-riffle complexes. Postglacial surficial geology ranged from coarse-grained outwash sand and gravel to lacustrine silt and clay, allowing for a range of stream substrate sizes. Presettlement riparian vegetation was mainly forest, including forested uplands, forested lowlands, and to a lesser extent, conifer-dominated wetlands in headwaters. This resulting framework of geomorphic habitat response units can be used for habitat rehabilitation projects for migratory native fish in other urban Great Lakes tributaries.

Wisconsin↗
USGS2025

Advancing climate adaptation for inland fish and fisheries

Climate change was barely a blip on the radar for inland fisheries management 20 years ago. Today, it's a central focus. A 2016 paper helped shift conversations, sparked adaptation efforts across fisheries management. The future is still uncertain, but adaptation is key to sustaining these important resources.

Wisconsin↗
USGS2025

A short-term sublethal oral exposure to microcystin-LR disrupts cecal microbiome homeostasis in mallard

Introduction: The frequency of cyanobacterial blooms seems to have increased globally in recent decades due to human induced eutrophication and climate change. Cyanobacterial blooms can produce several groups of toxins, among which microcystin-LR (MC-LR) is one of the most abundant. Effects of MC-LR on avian microbiome have not been studied and studies in laboratory murines have been limited to metabarcoding of prokaryotes. Methods: Using RNA shotgun sequencing, we compared the richness and composition of metabolically active prokaryotes, expressed virulence factors, antimicrobial resistance genes, metabolic pathways, Gene Ontology terms, enzymes, and proteins in mallards (Anas platyrhynchos) that were orally exposed to a sublethal dose of MC-LR for one week and unexposed birds. Results: Richness and composition of all compared features did not differ between exposed and control birds and none were differentially expressed between exposure groups. However, richness and/or composition of all features except virulence factors and Carbohydrate Active enzymes had multiple-fold greater dispersion in exposed birds than in controls. This effect was especially pronounced in expressed metabolic (MetaCyc) pathways. Discussion: Our results suggest that MC-LR exposure had a stochastic (rather than deterministic) effect on cecal microbiota, especially its function. Observed disturbance of the microbiota homeostasis is consistent with the Anna Karenina Principle. This principle has been documented in a wide range of eukaryotes using primarily microbial community metabarcoding. Although stochastic disturbance of microbiota function has been hypothesized, our study seems to be the first to demonstrate this in an experimental study.

Wisconsin↗
USGS2025

Assessing flood water infiltration and storage in a restored floodplain

In urban areas, floodplain restoration is gaining prominence as a strategy for restoring the natural functions of floodplain ecosystems and reducing flood risk. This has spurred research into potential interactions between floodwaters, the hyporheic zone, and the floodplain aquifer. An urban restored stream in Wisconsin, USA, was used as a case study to examine four methods to estimate floodplain infiltration and storage during overbank floods. We characterised flood-related infiltration over a 4-year period from 2018 through 2021 by simultaneously and continuously measuring groundwater levels and vertical temperature profiles with stream water levels linked to high-resolution flood inundation maps. High-resolution topographic data helped to quantify surface floodplain storage and the unsaturated soil volume relative to flood stage. Infiltration estimates from the simple methods align well with those from the more complex methods; however, the complex methods provide additional insights about the factors influencing infiltration. Results from all methods indicate that the volume of water that vertically infiltrates during floods is likely small relative to the total volume of the flood, with 0.08%–0.52% of flood water infiltrating into the floodplain, on average. Spatially variable vertical hydraulic gradients, driven by flood depth, groundwater level, and permeability, imply heterogeneous patterns of infiltration across the floodplain. Gradients favourable for infiltration typically occurred during the onset of flooding but, over the study period, were mostly (98% of the time) favourable for groundwater discharge to the channel (non-flood periods). These findings highlight the importance of considering surface-groundwater dynamics, floodplain soils, and unsaturated floodplain volume in defining the benefits of floodplain infiltration for flood attenuation.

Wisconsin↗
USGS2025

The role of street cleaning on the water-quality performance of a stormwater treatment pond in Madison, Wisconsin

The U.S. Geological Survey, in cooperation with the Wisconsin Department of Natural Resources and the City of Madison, evaluated how street cleaning frequency influences the pollutant removal efficiency of a stormwater treatment pond in Madison, Wisconsin (2020–24). Paired influent and effluent samples were analyzed for nutrients, sediment, and chloride under a weekly and monthly street cleaning scenario. Results showed that less frequent cleaning (monthly frequency) led to higher pollutant accumulation on streets, increasing influent concentrations of nitrogen and sediment. This, in turn, allowed the pond to achieve higher overall load-reduction percentage compared to weekly cleaning, particularly for total suspended sediment, total nitrogen, and total phosphorus. Dissolved phosphorus was an exception where removal was significantly greater under weekly cleaning. One explanation could be related to internal phosphorus release from pond sediments under anoxic conditions. Nearly all events showed net export of chloride from the pond, with effluent loads exceeding influent loads for both street cleaning frequencies.

Wisconsin↗
USGS2025

A regional model comparison between MODPATH and MT3D of groundwater travel time distributions

Groundwater quality changes in wells and streams lag behind changes to land use due to groundwater travel times. Two contaminant transport methods were compared to assess differences in their simulated travel time distributions (TTDs) to streams and wells in the Wisconsin Central Sands. MODPATH simulates advective groundwater flow with particle tracking, while MT3D simulates age-mass using a finite difference solution without dispersion to allow for direct comparison of the two methods. MODPATH appropriately simulates groundwater TTDs from the water table to surface discharge but is subject to inaccuracies at weak-sink well cells due to the flow-model grid discretization and imprecise location of well discharge within well cells. MT3D better represents weak-sink well cells since it removes mass in proportion to the prescribed pumping rate, although travel time within well cells is neglected. Conversely, MT3D's treatment of surface water boundary cells is not as accurate as MODPATH because mass should be removed from the water table rather than the full cell volume. MT3D simulations of TTDs can also be confounded by the instantaneous vertical distribution of mass introduced throughout recharge cells instead of at the water table, which initiates mass along deeper flow paths. We evaluated 9 MODPATH and 13 MT3D implementations, generating differences in median travel times of up to 18 years. Both methods have strengths and weaknesses, with MT3D better representing weak-sink well cell behavior and MODPATH better representing surficial recharge and discharge. The effect of these characteristics on simulated TTDs, along with ideas for ameliorating method weaknesses, is discussed.

Wisconsin↗
USGS2025

Placing environmental DNA monitoring for new detections into perspective: Fishes in the Milwaukee River, Wisconsin

Invasive species management frameworks, such as the early detection of and rapid response to invasive species, use monitoring programs to detect new species occurrences. Resource managers use environmental DNA (eDNA) as one tool for these monitoring programs. An eDNA detection in a new location may lack perspective for resource managers and researchers because of the rarity of potential invaders and the randomness in their dispersal and detection. An example monitoring program is the eDNA-based sampling approach used by the U.S. Fish and Wildlife Service for bigheaded carps Hypophthalmichthys spp. in the upper Mississippi River and Great Lakes Basins that collects hundreds of water samples per event. The U.S. Fish and Wildlife Service detected a single positive sample for Bighead Carp Hypophthalmichthys nobilis during the spring 2021 sampling event in the Kinnickinnic River within the Milwaukee River Basin, and detected a second single positive sample for bigheaded carps during the fall 2021 sampling event in the Milwaukee River. The U.S. Fish and Wildlife Service did not detect any bigheaded carps in previous years (2015 to 2020) or in either the spring or fall 2022 sampling events. These detections lacked perspective, such as detection numbers for other species. We reanalyzed the 2021 and 2022 samples for four existing species of fish: two fairly common species (Common Carp Cyprinus carpio and Gizzard Shad Dorosoma cepedianum ) and two fairly rare species (Burbot Lota lota and Grass Carp Ctenopharyngodon idella ). We detected Common Carp during all four sampling events, Gizzard Shad during three of four sampling events, and Burbot and Grass Carp during two of four sampling events. These results demonstrated that current sampling efforts could detect other species, and bigheaded carp eDNA was not common in the Milwaukee River compared to these species. More specifically, this finding indicates bigheaded carp eDNA detections are as rare as, or rarer than, Grass Carp eDNA detections, a recent invader to the basin. Our findings also demonstrated how reanalyzing eDNA samples after positive detections for targeted species can help managers understand the context of the detections and provide perspective for the relative abundance of the targeted species. Additionally, our results highlight the importance of completing long-term eDNA-based monitoring rather than a single sampling or inventory event. These detections may have been missed in a single year or sampling event, whereas a multiyear monitoring program provides an opportunity to observe trends through time.

Wisconsin↗
USGS2025

Differential responses of coolwater fishes to a whole-lake reduction of a warmwater thermal guild

Climate change is transforming the ecology of lakes at a rapid pace, shifting some lakes toward warmwater-dominant habitats. As a result, warmwater fishes are increasingly becoming more prevalent in lakes where they already existed, altering the patterning and strength of species interactions. Understanding shifting species interactions (e.g., competition, predation), and the role of lake management in shaping these interactions, will be critical for lake stewardship in response to climate change. Here, we present results from an intensive 5-year experimental removal of ~285,000 warmwater fishes from a north-temperate lake. The goal of the experiment was to test whether warmwater fish reduction is effective for rewiring lake food webs to reverse undesirable conditions for coolwater species, leading to increased recruitment and abundance of coolwater fishes. Throughout the experiment, warmwater fishes were resilient to reductions, with biomass declines of 23% averaged across five species. Among coolwater fishes, the top predator walleye showed no biomass response, while yellow perch biomass increased by ~914%. Fish species biomass changes translated to food web shifts, including a yellow perch trophic position decline of 0.4, decreased zooplankton abundances, and increased zoobenthos abundances. Our results highlight differential species responses to a management action aimed at adapting to climate change. Despite similar thermal tolerances, two coolwater species responded differently to removal of warmwater fishes, highlighting the characteristics (e.g., life history strategies, adaptive capacity) that contribute to species resilience. Given the importance of biotic interactions, climate adaptation may need to go beyond a “one-size-fits-all” approach even when species have similar thermal tolerances.

Wisconsin↗
USGS2025

The breeding season and movement ecology of male white‐tailed deer in southwest Wisconsin

White-tailed deer ( Odocoileus virginianus ; hereafter, deer) have been widely studied regarding their breeding ecology and responses to hunting pressures. However, variations in defining the breeding season—its duration and timing—across studies have created uncertainty about whether regional differences in deer breeding ecology stem from ecological factors or methodological inconsistencies. This study aims to clarify the peak breeding season timing and the movement patterns of males during this period, particularly in relation to hunting seasons. Understanding how age and the timing of hunting seasons impact movement and breeding behaviors is important for wildlife managers, as these factors can affect harvest success. This study took place in southwest Wisconsin, using GPS data collected from 188 collared male deer between 15 October and 1 December from 2017 to 2020. Based on generalized linear mixed models, 2-year-old males exhibited higher hourly movement rates than other ages, and the opening weekend of the firearm hunting season had no significant effect on movement rates. In contrast, the variance in daily movement rate differed significantly between yearlings and older ages, with males 3 years and older displaying the highest variance. This suggests that older males may alternate more frequently between high-movement mate searching and lower-movement mate tending, potentially enhancing reproductive success. Similarly, 2-year-old males had larger daily ranges than both older and younger ages. Changepoint analysis of daily movement rates determined that the peak breeding season occurred between 23 October and 12 November, with little variation among ages and alternative metrics. Our findings indicate that male movement rates and ranges can reflect deer reproductive efforts and vary by age, which has important implications for reproductive success and disease transmission risk.

Wisconsin↗
USGS2025

Estimating mortality of Lake Sturgeon in the Lake Winnebago system using traditional age-based approaches and capture–recapture models

Objective The Lake Winnebago system in Wisconsin supports a popular winter spear fishery for Lake Sturgeon Acipenser fulvescens . Setting harvest caps for this fishery relies on estimating instantaneous natural mortality rate ( M ), which can be done using age-based approaches or capture–recapture models that incorporate recoveries of fish with passive integrated transponder (PIT) tags or detections of fish with acoustic transmitters. Our objectives were to determine (1) if recent estimates of exploitation ( u ) have exceeded the 5% harvest cap, (2) if M and total mortality rates are similar among estimation methods that rely on age estimates or capture–recapture methods, and (3) if potential differences in mortality estimates would affect harvest caps. Methods Harvest of PIT-tagged fish was used to evaluate u from 2010 to 2019. Catch curves incorporating corrected fin ray ages were used to estimate total mortality and M for fish collected from 2010 to 2019. Capture–recapture models were used to estimate annual survival and M from detections of fish with acoustic transmitters from 2007 to 2019 and recoveries of PIT-tagged fish from 1999 to 2020. Mortality estimates were used to calculate and compare sex-specific harvest caps among estimation methods. Results Observed u did not exceed 5% for either sex between 2010 and 2019. Estimates of M varied among methods (males: M = 0.001–0.134; females: M = 0.001–0.131), with PIT-based models consistently providing the lowest and telemetry-based models providing the highest estimates. Simulations indicated that female u has limited potential to exceed 5% if M from fin ray ages or telemetry is used to set harvest caps, while PIT-based simulations showed no indication of cap exceedance. Conclusions Harvest management practices in the Lake Winnebago system appear to have kept Lake Sturgeon exploitation below the 5% harvest cap from 2010 to 2019. Capture–recapture models relying on PIT tags appear to provide the most precise approach for setting harvest caps for this fishery.

Wisconsin↗
USGS2025

Comparative toxicity of aircraft deicers, runway deicers, and road salt in winter airport runoff

Aircraft deicer and anti-icer formulations and airfield pavement deicer formulations are commonly used at airports that experience freezing precipitation. In addition, road salt is often applied in urban areas surrounding airports. This mixture of deicers in receiving streams can be toxic to aquatic organisms. A 17-year study of deicers and toxicity in runoff samples from sites in and near Milwaukee Mitchell International Airport included bioassay testing of Ceriodaphnia dubia , Pimephales promelas , Raphidocelis subcapitata , and Aliivibrio fischeri concurrently with quantification of freezing-point depressants that represent presence of the various deicers. Propylene glycol was monitored to represent aircraft deicers and anti-icers, acetate and formate were monitored to represent airfield pavement deicers, and chloride was monitored to represent road salt. Toxicity endpoints (e.g., lethal concentrations and inhibitive concentrations) were compiled for 19 airport deicer products and sodium chloride. Comparison of freezing-point depressant concentrations with toxicity endpoints indicated that propylene glycol, acetate, and chloride exceeded endpoint concentrations in multiple samples: more than 70% of 492 runoff samples exceeded at least one endpoint concentration with exceedances of propylene glycol and acetate-based deicer endpoints as well as road salt endpoints. Validation of these results were observed with measured toxic effects on C. dubia , P. promelas , and R. subcapitata in bioassays conducted on undiluted samples: samples with multiple deicers exceeding endpoint thresholds had greater toxicity than samples with individual deicer endpoint exceedances, and samples with no endpoint exceedances typically exhibited less toxicity than samples with only one endpoint exceedance. A comparison of additivity of toxicity quotients with bioassay results fit a log-logistic regression curve, indicating that the condition of tested organisms degraded with increasing deicer presence. These findings suggest that concurrent reduction of runoff from aircraft, airfield pavement, and road salt deicers will be needed to eliminate aquatic toxicity in adjacent watersheds.

Wisconsin↗
USGS2025

Seasonal stratification drives bioaccumulation of pelagic mercury sources in eutrophic lakes

Increased lake eutrophication, influenced by changing climate and land use, alters aquatic cycling and bioaccumulation of mercury (Hg). Additionally, seasonally dynamic lake circulation and plankton community composition can confound our ability to predict changes in biological Hg concentrations and sources. To assess temporal variation, we examined seasonal total Hg (THg) and methylmercury (MeHg) concentrations and stable isotope values in seston, waters, sediments, and fish from two adjacent urban eutrophic lakes in Madison, Wisconsin. In Lake Monona, surface sediment THg concentrations were elevated due to comparably higher urbanization and historical industrial inputs, whereas Lake Mendota sediments had lower concentrations corresponding with largely agricultural and suburban surrounding watershed. Surface water THg and MeHg were similar between lakes and seasonally dynamic, but water profiles exhibited elevated concentrations in the meta- and hypolimnion, highlighting water column MeHg production. Seston MeHg concentrations were often highest at shoulder seasons, possibly owing to metalimnetic MeHg delivery, but also differences in biomass and water clarity. The ∆199Hg and δ202Hg values in seston were similar between lakes, despite differing sediment THg concentrations and values, suggesting a shared bioaccumulated source of MeHg. Measurement of MeHg stable isotopes further elucidated that seston and fish predominantly bioaccumulated pelagic-sourced MeHg.

Wisconsin↗
USGS2025

Sediment nutrient dynamics in selected Milwaukee metropolitan area streams, Wisconsin, 2022

The U.S. Geological Survey and Milwaukee Metropolitan Sewerage District in Wisconsin have an ongoing partnership to monitor water quality in streams in the Milwaukee metropolitan area and to assess the effects of stream restoration on habitat and water quality. Because sediment nutrient dynamics can improve or further impair water quality, we measured sediment nitrogen and phosphorus concentrations, potential nitrogen removal, and potential phosphorus retention at 32 streams sites in the Milwaukee metropolitan area in summer 2022. Four of the sites were in rehabilitated stream reaches. Based on the results from this project, we provided a rating of good, fair, or poor for the sediment nutrient dynamics at each site. Sediment nitrogen removal and phosphorus retention in stream reaches increased as the proportion of fine particles increased. Bioavailable nitrogen and phosphorus concentrations increased in stream reaches as particle size decreased, especially at locations with more silt and clay. Particle size typically decreased from upstream to downstream, and several of the sites with finer particles were in downstream parts of the study area, especially in the Milwaukee Estuary sites. The sites that had more fine sediment and higher bioavailable nutrients also had elevated rates of nitrification and denitrification enzyme activity, which is a measurement of the potential nitrogen removal when nitrogen and carbon do not limit the denitrifying bacteria. Additionally, sites with fine sediment typically had the highest potential to retain phosphorus as fine sediment provides many binding sites for phosphorus. The binding sites can become saturated with phosphorus, however, increasing the potential for the phosphorus to be released to stream water. Five of the 32 sampled sites were potentially saturated with phosphorus. Sites that contained more agricultural land in their drainage areas were at higher risk of having sediment that was saturated or near-saturation; however, the sites that had more agricultural land in their drainage areas also had higher nitrification rates. Results from this study indicate that stream rehabilitation projects that promote sediment deposition and accumulation of organic matter in the stream channel can increase nitrogen removal and phosphorus retention.

Wisconsin↗
USGS2025

Sulfate reduction drives elevated methylmercury formation in water column of eutrophic freshwater lake

Mercury (Hg) contamination of aquatic food webs is controlled in part by the formation and accumulation of toxic and bioaccumulative methylmercury (MeHg). MeHg production is mediated by metabolically diverse microorganisms carrying the hgcAB gene pair, while the demethylation reaction is mediated by several biotic and abiotic processes. However, the relative importance of these two processes on MeHg accumulation and the environmental factors that influence them are poorly characterized, especially in eutrophic environments. In this study, both Hg methylation and MeHg demethylation in a eutrophic freshwater lake were linked to ambient MeHg concentrations and hgcA abundance and expression. High methylation rate potentials indicated in situ MeHg formation was a key source of MeHg to the water column, driven by high hgcA abundance and transcription. Molybdate treatment decreased methylation rate potentials, highlighting the importance of sulfate reduction in driving MeHg formation. Sulfate-reducing bacteria accounted for over 50% of the hgcA gene transcription, despite representing less than 10% of the hgcA -carrying microbial community. An arsR -like transcriptional regulator preceded many hgcA sequences; these were transcriptionally active and linked to lower hgcA expression. Overall, this study elucidates the microbial and biogeochemical processes that influence the in situ formation of MeHg in understudied eutrophic freshwater environments.

Wisconsin↗
USGS2025

Responses of larval fish and zooplankton to whole-lake 2,4-D herbicide treatments used to control Eurasian watermilfoil in northern Wisconsin lakes

In North America, Eurasian watermilfoil Myriophyllum spicatum is one of the most widespread non-native aquatic plant species in freshwater ecosystems. Applications of 2,4-dichlorophenoxyacetic acid (2,4-D) herbicides are often utilized to control Eurasian watermilfoil. Herbicide applications may have unintended effects on non-target organisms like zooplankton and fish, but these effects are poorly understood. Our objectives were to determine if whole-lake 2,4-D herbicide treatments (epilimnetic target rate = 0.3 ppm) used to control Eurasian watermilfoil affected (1) density, diversity, and size of zooplankton and (2) density, diversity, and growth of larval fishes in three northern Wisconsin lakes when compared to trends observed in three untreated reference lakes. Average peak concentrations of 2,4-D following applications ranged from 0.152 to 0.257 ppm. We found no statistically significant evidence that the applications influenced the zooplankton and larval fish metrics we examined, but subtle responses to the herbicides may not have been detected given the inherent variation in our response metrics. Additional replication and research are needed to assess responses to higher concentrations and repeated applications over time and to document responses across a broader range of lake types and fish communities.

Wisconsin↗
USGS2025

Characterization of the long-distance dispersal kernel of white-tailed deer and evaluating its impact on chronic wasting disease spread in Wisconsin

Chronic wasting disease (CWD) is a fatal neurodegenerative disease infecting cervids. It is highly contagious and caused by misfolded prions that propagate via templated conformational conversion of the cervid’s normal prion protein. Prevalence of CWD in free-ranging deer in North America is mostly low, but in some regions local prevalence has reached 80%. CWD prions can be transmitted via direct contact with infected individuals or indirectly through the environment. Infected individuals shed prions through feces, urine, saliva or carcasses, and prions have long environmental persistence. Long-distance dispersal of infected deer poses a significant risk for CWD spread. We propose an integrodifference equation (IDE) model to capture CWD dynamics and the consequences of long-distance dispersal behavior in white-tailed deer (WTD, Odocoileus virginianus ). A diffusion-settling model characterizes long-distance dispersal kernels, accommodating hypothetical dispersal behaviors through time-dependent settling rate functions. Three new closed-form dispersal kernels are approximated using Laplace’s method and parameterized with GPS location data collected from WTD in Wisconsin, USA. Settling rates reflecting ongoing sensitivity to stimuli which prompt deer to disperse from their natal home range give the most supported long-distance dispersal kernel. Impact of long-distance dispersal on CWD spread is quantified using the IDE model. At high population densities, long-distance dispersal can magnify CWD spread by a factor of four. At lower population densities single infected individuals cannot initiate an outbreak, but CWD may still spread due to the accumulation of environmental hazard from prions behind the wave of invasion, possibly presenting substantial management challenges.

Wisconsin↗
USGS2025

Spatiotemporal causal inference with mechanistic ecological models: Evaluating targeted culling on chronic wasting disease dynamics in cervids

Spatiotemporal causal inference methods are needed to detect the effect of interventions on indirectly measured epidemiological outcomes that go beyond studying spatiotemporal correlations. Chronic wasting disease (CWD) causes neurological degeneration and eventual death to white-tailed deer ( Odocoileus virginianus ) in Wisconsin. Targeted culling involves removing deer after traditional hunting seasons in areas with high CWD prevalence. The evaluation of the causal effects of targeted culling in the spread and growth of CWD is an important unresolved research and CWD management question that can guide surveillance efforts. Reaction–diffusion partial differential equations (PDEs) can be used to mechanistically model the underlying spatiotemporal dynamics of wildlife diseases, like CWD, allowing researchers to make inference about unobserved epidemiological quantities. These models indirectly regress spatiotemporal covariates on diffusion and growth rates parameterizing such PDEs, obtaining associational conclusions. In this work we develop an innovative method to obtain causal estimators for the effect of targeted culling interventions on CWD epidemiological processes using an inverse-probability-of-treatment-weighted technique by means of marginal structural models embedded in the PDE fitting process. Additionally we establish a novel scheme for sensitivity analysis under unmeasured confounder for testing the hypothesis of a significant causal effect in the indirectly measured epidemiological outcomes. Our methods can be broadly used to study the impact of spatiotemporal interventions and treatment exposures in the epidemiological evolution of infectious diseases that can help to inform future efforts to mitigate public health implications and wildlife disease burden.

Wisconsin↗