Investigating the relationship between extreme weather and cryptosporidiosis and giardiasis in Colorado: a multi-decade study using distributed-lag nonlinear models
====================================================================================================================================================================

* Elise N. Grover
* James L. Crooks
* Elizabeth J. Carlton
* Sara H. Paull
* William B. Allshouse
* Katherine A. James

## Abstract

**Background** Environmentally-mediated protozoan diseases like cryptosporidiosis and giardiasis are likely to be highly impacted by extreme weather, as climate-related conditions like temperature and precipitation have been linked to their survival, distribution, and overall transmission success.

**Objectives** Our aim was to investigate the relationship between extreme temperature and precipitation and cryptosporidiosis and giardiasis infection using monthly weather data and case reports from Colorado counties over a twenty-one year period.

**Methods** Data on reportable diseases and weather among Colorado counties were collected using the Colorado Electronic Disease Reporting System (CEDRS) and the Daily Surface Weather and Climatological Summaries (Daymet) Version 3 dataset, respectively. We used a conditional Poisson distributed-lag nonlinear modeling approach to estimate the lagged association (between 0 and 12-months) between relative temperature and precipitation extremes and the risk of cryptosporidiosis and giardiasis infection in Colorado counties between 1997 – 2017, relative to the risk found at average values of temperature and precipitation for a given county and month.

**Results** We found a consistent, significant increase in the relative risk of cryptosporidiosis and giardiasis infection when the maximum or minimum monthly temperature in the 4-6 months prior were low (1st percentile) for a given county and calendar month. High precipitation (90th percentile for a given county and calendar month) 2-12 months prior was associated with a significant decrease in the relative risk of cryptosporidiosis, while high precipitation (90th percentile) 12-months prior was associated with a significant increase in giardiasis infection risk.

**Discussion** Our study presents novel insights on the influence that extreme low temperatures can have on parasitic disease transmission in real-world settings. Additionally, we present preliminary evidence that the standard lag periods that are typically used in epidemiological studies to assess the impacts of extreme weather on cryptosporidiosis and giardiasis may not be capturing the entire relevant period.

## Introduction

In recent years, a diverse array of extreme weather events across the globe have been inextricably linked to human-induced climate change – including heatwaves, heavy precipitation, flooding and droughts – all of which are likely to impact the transmission of environmentally-mediated illnesses such as those transmitted by parasitic protozoa. 1–4 Protozoan diarrheal diseases such as cryptosporidiosis and giardiasis are a leading cause of the estimated 1.7 billion cases of diarrheal disease that occur every year, 8–12 with an estimated 33,900 deaths and 2.94 million disability-adjusted life-years (DALYs) lost to illness caused by enteric protozoa each year. 13 *Cryptosporidium* and *Giardia* are both zoonotic protozoan parasites that are transmitted by animals and humans via the fecal-oral route, and are most frequently acquired by consuming contaminated water or food. 14–17 As such, regions of the world that have limited access to clean water, are prone to water scarcity, or lack infrastructure to support widespread access to improved sanitation tend to be the most at risk of protozoan disease transmission. 18 Nevertheless, climate change is also expected to increase waterborne enteric disease risk in high income settings as a result of projected increases in the frequency and severity of extreme weather events. 1,4–7

Across the globe, both periods of drought and heavy rainfall can cause inefficiencies and reduced effectiveness of water treatment systems. 2,19–23 For example, studies conducted in Australia and the Netherlands have found that periods of drought were associated with decreased water quality, 25,26 and an increased risk of cryptosporidiosis. 24 Meanwhile, periods of heavy rainfall can also cause sewage system overflow, contamination of irrigation systems, wells and private water supplies, and associated food and waterborne illnesses. 2,27 Protozoa are particularly susceptible to mobilization and transmission following heavy rainfall events due to their ability to persist in soil and water for months, their capacity to infect individuals at very small doses, and their ability to be transmitted by humans, domestic animals and livestock. 12,28–30 For instance, researchers in Vancouver, Canada found that water turbidity and the risk of runoff-related cryptosporidiosis and giardiasis incidence was highest when a period of excessive dryness was followed by heavy rainfall events. 23 Exposure via runoff and waste water system failures may be particularly important for parasitic protozoan diseases, which are more resilient than most bacteria and viruses in the face of conventional water treatment methods, 32 and are less commonly the target of operational surveillance systems. 4

When it comes to the impacts of temperature extremes on diarrheal diseases, protozoan infections have been less studied than bacterial or viral infections, 33,34 though associations between maximum and minimum temperature and cryptosporidiosis and giardiasis have been suggested in a few instances. A review by Ikiroma & Pollock (2021) on the impacts of weather on cryptosporidiosis identified studies highlighting a positive association between maximum temperature and cryptosporidiosis, 35–39 and to a lesser extent, giardiasis. 38 While further mechanistic studies are needed to determine why rising temperatures are associated with an increase in protozoan disease risk in some contexts, possible explanations include increased potential for transmission via alternative vectors at higher temperatures, or the lengthening of the transmission season when temperatures are unseasonably warm during the offseason. 40–43 By contrast, very few studies have investigated the impacts of minimum temperature, and the results have been inconsistent. 39 For example, one study from Victoria, Australia found that an increase in monthly minimum temperature in the current month (no lag) was associated with an increase in cryptosporidiosis incidence in metropolitan areas, whereas the 3-month lagged monthly minimum temperature in rural areas was associated with a decrease in cryptosporidiosis incidence. 44 Further complicating our understanding of the influence of temperature on enteric protozoan disease transmission is the fact that these relationships are pathogen-specific and highly dependent on the local climate system and season. 27,38,45 For instance, a meta-analysis conducted by Jagai and colleagues (2009) demonstrated that while temperature was more strongly associated with cryptosporidiosis in temperate climates, precipitation was a stronger predictor in tropical climates, though notably, neither temperature nor precipitation was significantly associated with cryptosporidiosis in arid or semi-arid climates. 46

Although the majority of studies assessing the effects of weather and climate on cryptosporidiosis and giardiasis have assessed lag periods of between one and three months, 35–38,44 pathogen survival in the ambient environment is likely longer under optimal conditions. Survival time of *Cryptosporidium* oocysts and *Giardia* cysts will vary substantially with local environmental conditions such as ultraviolet light (UV), temperature, humidity, precipitation, soil composition, season, daylight hours, excrement moisture, water velocity, water turbidity, etc. 12,28,29,40,41,47–51 As such, determining the maximum survival time in the ambient environment can be challenging, though controlled studies have demonstrated that *Giardia* cysts submerged in river and lake water at winter temperatures maintained viability at 12-weeks post-baseline, 47 while *Cryptosporidium* oocyst infectivity can remain after being stored in water for over five months at 5-15°C, 51 and up to a year when stored in low turbidity water. 52 Given remaining uncertainty regarding the maximum survival of *Cryptosporidium* and *Giardia* and their demonstrated capacity for prolonged survival under optimal conditions, epidemiological studies that assess the potential for longer lagged effects are needed.

In light of the many influences that weather and environment have on the transmission of protozoan pathogens, as well as the remaining uncertainty surrounding the time- and location-specific conditions that may amplify or reduce risk, localized investigations of weather extremes and their lagged effects on transmission are needed. Supplemental Table S1 provides examples from the literature and a summary of some of the hypothesized mechanisms by which extreme precipitation and temperature can influence cryptosporidiosis and giardiasis transmission. Overall, despite mixed evidence on the direction of association between precipitation extremes and protozoan disease transmission, and the limited body of epidemiological evidence on the relationship between temperature extremes and cryptosporidiosis and giardiasis risk, there remain several plausible mechanisms by which weather extremes could affect the viability, distribution and overall transmission potential of *Cryptosporidium* and *Giardia*.

In this study, we investigated the lagged effects (between 0 and 12 months) that low and high temperature and precipitation values – defined in this study as the 1st, 5th, 10th, 90th, 95th and 99th percentile values of our weather variables – each have on the risk of cryptosporidiosis and giardiasis infection in Colorado counties between 1997 – 2017, relative to the risk found at average values of temperature and precipitation. In so doing, we provide valuable insight on the Colorado-specific lagged effects that precipitation and temperature extremes have on two important, but often overlooked sources of waterborne illness. As extreme weather events are likely to increase in frequency and strength in the decades to come, our findings contribute to the larger body of science on climate and health, helping to characterize the diverse and complex effects that climate change can have on human health.

## Methods

### Study design

In this study, we assessed the impacts of monthly precipitation and temperatures on reported cases of cryptosporidiosis and giardiasis in Colorado counties between 1997 and 2017. Specifically, our primary objective was examining the effects of relative temperature and precipitation extremes – mean maximum monthly temperature (MAXT), mean minimum monthly temperature (MINT), and total monthly precipitation (PREC), each mean-centered by county and calendar month – on the monthly case counts of human cryptosporidiosis and giardiasis. As secondary goals, we also identified relevant lag periods over which extreme weather may impact cryptosporidiosis and giardiasis in Colorado, and evaluated whether the effects of precipitation and temperature on monthly case counts varied by season. We used a distributed lag nonlinear modeling (DLNM) approach, conditioned on total case counts within each of the twelve calendar-months for a given county. 53,54 The Colorado Multiple Institution Review Board (COMIRB) reviewed and approved of this study.

### Cryptosporidiosis and giardiasis incidence

Data on reportable diseases in Colorado were collected using the Colorado Electronic Disease Reporting System (CEDRS), which is made available upon request from the Colorado Department of Public Health and Environment (CDPHE). Our study spanned the 21-year period between January 1997 through December 2017, using monthly case reports of cryptosporidiosis and giardiasis by county as our two outcomes of interest. As both cryptosporidiosis and giardiasis are mandatory reportable diseases and CDPHE conducted lab audits to ensure proper monthly reporting, for the purposes of this study we assumed that cases of sufficient severity to warrant medical attention in Colorado counties for each month between 1997 and 2017 were captured in the CEDRS dataset.

In order to estimate incidence, we used annual county population estimates from the Colorado Department of Local Affairs. 55 Although Broomfield County was approved by voters as Colorado’s 64th county in 1998, it was not included as its own county in this analysis until 2000, as prior to this point, the Colorado Department of Local Affairs census data recorded Broomfield’s population as a part of Boulder County’s population. The natural log of each county’s annual population estimate was included as an offset term in the statistical models.

### Climatological variables

To obtain monthly estimates of temperature and precipitation for each Colorado county, we used the Daily Surface Weather and Climatological Summaries (Daymet) Version 3 dataset derived by the National Center for Atmospheric Research (NCAR), which provides daily estimates of weather parameters across a 1-kilometer grid spacing for North America. 56 For each month between January 1996 and December 2017, MAXT (°C) and the MINT (°C) were each averaged across all 1-km grid points within a given county, such that each county-month had a single estimated average value. PREC for each county-month was estimated using total daily precipitation (mm) measured at each 1-km grid point within a county’s borders, summed across all days in each month and averaged across all grid points in the county. All forms of precipitation were included (e.g., snow, sleet, rain) and converted to a water-equivalent depth in millimeters. Notably, while our outcome data begins in 1997, the data for the explanatory variables begins in January 1996 to allow for up to 12-month lagged effects for each of our weather parameters.

To make the definition of low and high most relevant to the local weather conditions found in Colorado counties, we mean-centered our weather variables by county and calendar-month, such that values could be interpreted as hot/cold or wet/dry for the specific location and time of year. All three weather variables were subsequently mean-centered by county and calendar-month, allowing MAXT, MINT and PREC to be defined relative to their location and time. For example, the average MAXT was 6.3°C for Adams County in December 2017, while the mean-centered (MC)-MAXT was +2.6°C (as compared to other Decembers between 1996 – 2017 in Adams County), indicating a warmer than average December in Adams County in 2017.

In this analysis, we were ultimately interested in the association between higher and lower than average temperature and precipitation values on cryptosporidiosis and giardiasis incidence. Thus, to define the set of relevant predictor values to be evaluated as the numerator in our relative risk calculations, we used the 1st, 5th, 10th, 90th, 95th and 99th percentile values of our mean-centered weather variables, referred to respectively as extreme low, very low, low, high, very high and extreme high values throughout this analysis, as detailed in Supplemental Table S2.

### Analysis

To assess the lagged effects of MC-MAXT, MC-MINT and MC-PREC on human cryptosporidiosis and giardiasis cases in Colorado counties between 1997-2017, we used a conditional Poisson DLNM approach. 53,54 DLNMs were used because they provide a flexible platform for exploring nonlinear outcome-predictor relationships and they permit lagged effects to be spread over time, two characteristics that are desirable when exploring the impacts of highly variable environmental conditions on health outcomes. As such, we investigated a wide range of potential lags between 0- and 12-months by varying the number of internal knots (0-6) for the spline functions used to define the lag space and comparing model performance across these formulations of the lag space. Whereas 2-month and 4-month lags were selected because they aligned with the typical lag range indicated in the literature, 35–39,44 6-month and 12-month lags was also examined in detail in order to determine whether extreme weather could have more prolonged impacts of on giardiasis or cryptosporidiosis case counts in Colorado counties between 1997 – 2017. Ultimately, this design allowed us to identify potential longer-term trends that might have otherwise been missed if only short-term lags were assessed.The impacts of MC-MAXT, MC-MINT and MC-PREC on monthly cryptosporidiosis and giardiasis county case counts were each assessed independently, making a total of six outcome-predictor combinations. The steps that were taken to conduct this analysis for each of the six outcome-predictor pairs of interest to this study are outlined in Figure 1.

![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/09/01/2023.08.31.23294911/F1.medium.gif)

[Figure 1.](http://medrxiv.org/content/early/2023/09/01/2023.08.31.23294911/F1)

Figure 1. Illustration of the analytical steps used in this study.
The 8-Step analytical process that was used for each outcome and predictor of interest in this study is outlined above. The two outcomes in this study were case counts of cryptosporidiosis, and case counts of giardiasis. The three predictors were the average maximum monthly temperature, the average minimum monthly temperature, and the total monthly precipitation. We therefore repeated Steps 1-8 a total of six times for this study, once for each outcome-predictor pair.

In Step 1, initial exploratory analyses were conducted to summarize between-season differences in our outcome-predictor pairs across our study period. A three-month definition of season was also used in our exploratory assessment to investigate the distribution of the data across the calendar year. Winter was defined as December – February; spring as March – May; summer as June – August; and fall as September – November.

Next, we developed several different cross-basis structures to investigate potential non-linearities in the lagged predictor-outcome relationship over time and space, using the “dlnm” package in R. 58,59 A cross-basis is a bi-dimensional space of functions that encapsulates both the shape of the relationship between the outcome and predictors (termed the “predictor space”), as well as the distributed lagged effects (i.e. the “lag space”). 53 To develop a cross-basis with suitable functions defining both the predictor space and the lag space, we investigated the fit of different third-degree spline functions (natural spine, B-spline and penalized spline) with incrementally increasing degrees of freedom (maximum of seven) and a maximum lag of 12-months for the lag space, while the predictor space was held constant using a simple linear basis function. The default settings for knot placements were used for each spline function such that internal knots were placed at equally spaced points along the lag-space while the boundary knots were set to 0- and 12-months for the natural and B-spline functions. A conditional Poisson generalized non-linear modeling (GNM) approach was used to compare model fit across the different lag space functions used in the cross-basis. All models also included an intercept, a continuous year term, an offset for county population and county-month as the stratum (i.e., conditional) variable. Each model’s Bayesian information criterion (BIC) was used to identify the three best lag space definitions for the cross-basis.

Once suitable functions were identified for the lag space, Step 3 was to develop the predictor space using a similar process for the lag space. That is, the fit of different third-degree spline functions with incrementally increasing degrees of freedom (up to a maximum of five) were compared for the predictor space, while we defined the lag space using the top three previously identified lag space formulations. In addition to assessing different variations of third-degree spline functions, threshold functions were also assessed for the predictor space. The BIC value was used again to select the three best fit models across all the different cross-basis formulations.

We looked for evidence of potential effect modification by season in Step 4. Four separate cross-basis terms (one for each three-month season) were generated, using the top performing predictor space and lag space definitions previously identified through Steps 2 and 3. In total, the procedure outlined in steps 2-4 yielded a set of 50 models for each outcome-predictor pair, which was then compared a final time using both the Akaike information criterion (AIC) and BIC values to select a final subset of six models that were reserved for further exploration. To demonstrate the analytical process that was used to build suitable cross-basis structures in Steps 2-4, Supplemental Table S3 provides an example from one outcome-predictor pair, detailing all model variants that were developed for MC-MINT and cryptosporidiosis.

The predictive capacity of each of the six AIC/BIC-selected models for a given outcome-predictor was then assessed using cross-validation in Step 5. For each of the six models, repeated training and testing was conducted using a 21-fold cross-validation. 60 Each fold contained a year of observations, such that, within each fold, 20 years of data were used to train the model while the remaining year of county-month observations was used for validation. This training-validation process was repeated a total of twenty-one times, once for each year between 1997 and 2017. The mean RMSE value summarizing the performance across all tuning and validation iterations was then used as an indicator of the overall predictive skill of the model. The model with the lowest mean RMSE was selected as the final model (Supplemental Table S4).

The final model was re-run using a conditional quasi-Poisson GNM variant to obtain adjusted confidence intervals to account for overdispersion (Step 6). As was done with for the conditional Poisson models, the quasi-conditional Poisson models included an intercept outside of the cross-basis, a continuous year term to account for long-term trends, county-month as the stratum variable, and an offset for county population. In Step 7, we used these models to generate estimates of the relative risk (RR) with 95% confidence intervals (95%CI) of cryptosporidiosis or giardiasis across a matrix of relevant lag and predictor values, relative to the risk of disease found at zero (i.e., the mean for a given county and calendar-month). The lag periods that were used in these relative risk calculations were 2-months, 4-months, 6-months, and 12-months.

Finally, to visualize our relative risk estimates and the entire predicted outcome surface across all predictor and lag values, we created a set of plots depicting the predictor-lag-outcome surface, using a reference value of zero (i.e., the mean) in all cases. We used a p-value of <0.05 to indicate statistical significance in this study. Stata 15 (Stata Statistical Software: Release 15 (2017). StataCorp LP, College Station, TX) and R Studio 4.0 (RStudio Team (2020). RStudio: Integrated Development Environment for R. RStudio, PBC, Boston, MA URL [http://www.rstudio.com/](http://www.rstudio.com/)) were used for all analyses. 59,61

## Results

Between 1997 and 2017, the total population in Colorado nearly doubled, increasing from approximately 2.9 million to 5.6 million (Table 1). For all counties, PREC was highest in the spring and summer, and lowest in the fall and winter. Cryptosporidiosis and giardiasis cases tended to be highest in the summer (38.2% of all cryptosporidiosis cases; 30.6% of all giardiasis cases) and the fall (32.1% of all cryptosporidiosis cases; 30.3% of all giardiasis cases) across all Colorado counties (Table 1). Reported cryptosporidiosis cases increased over the study period while giardiasis cases decreased (Figure 2).

![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/09/01/2023.08.31.23294911/F2.medium.gif)

[Figure 2.](http://medrxiv.org/content/early/2023/09/01/2023.08.31.23294911/F2)

Figure 2. Reported cryptosporidiosis and giardiasis cases across all Colorado counties between 1997 – 2017.
Depiction of th annual count of cryptosporidiosis and giardiasis cases across all Colorado counties between 1997 – 2017. Whereas there was a modest decrease in the total number of cases of giardiasis, there was an increase in the total number of cryptosporidiosis cases in Colorado cases between 1997 and 2017.

View this table:
[Table 1.](http://medrxiv.org/content/early/2023/09/01/2023.08.31.23294911/T1)

Table 1. Colorado county cases, population and seasonal weather for Colorado counties.

Each of the final models selected via cross-validation for modeling the impacts of weather extremes on cryptosporidiosis cases used natural or B-splines in both the lag and predictor spaces, while the models for giardiasis cases ultimately employed a natural or B-spline in the lag space and a linear function to model the predictor space. Although some of the top six AIC/BIC-selected models included interaction terms for season, none of these models were selected as the final model via cross-validation for any outcome-predictor pair. Supplemental Table S4 details the cross-basis structure and the cross-validation results for the top six models of each outcome-predictor pair.

### Cryptosporidiosis

#### Temperature

Our final models indicated that there was an increase in cryptosporidiosis case counts when relative MC-MAXT or MC-MINT in the 4-6 months prior was at an extreme low (Figure 3, Table 2). When MC-MAXT was 4.5°C lower (1st percentile) than the county and calendar-month mean 4-months prior, there was a 5.6% increase in the risk of cryptosporidiosis (RR 1.056; 95% CI 1.007 – 1.107), while at 6-months prior, there was an 8.4% increase in the risk of cryptosporidiosis relative to the mean (RR 1.084; 95% CI 1.030 – 1.141). Similarly, when MC-MINT was 3.8°C below average 6-months prior (1st percentile value), there was a 45.3% increase in reported cryptosporidiosis relative to the risk at the county and calendar-month mean (RR=1.453; 95% CI 1.347 – 1.568). At even lower temperatures (< 1st percentile), the increase in the risk of cryptosporidiosis was even more pronounced, as demonstrated in Figure 3.

![Figure 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/09/01/2023.08.31.23294911/F3.medium.gif)

[Figure 3.](http://medrxiv.org/content/early/2023/09/01/2023.08.31.23294911/F3)

Figure 3. 3D depictions of the lagged effects of temperature and precipitation on cryptosporidiosis
3D illustrations of the best performing model for the lagged relationship between cryptosporidiosis cases and MC-MAXT, MC-MINT and MC-PREC. For each predictor-outcome pair, the 3D graphic depicts the change in the relative risk of cryptosporidiosis across the range of observed temperature and precipitation values for lags of 0 – 12 months, relative to the risk found at the county and calendar-month averages for temperature and precipitation.

View this table:
[Table 2.](http://medrxiv.org/content/early/2023/09/01/2023.08.31.23294911/T2)

Table 2. The relative risk of cryptosporidiosis in Colorado counties between 1997 – 2017, comparing temperature and precipitation extremes to county averages.

By contrast, when MC-MAXT was very low or low (5th and 10th percentiles, respectively), or when MC-MINT was low (10th percentile) 2-6 months prior, there was a decrease in cryptosporidiosis risk relative to the risk at the mean (Table 2). However, when a 12-month lag was used, the relationship between cryptosporidiosis and low temperatures (10th percentile) mimicked that of extreme low temperatures, wherein the risk of cryptosporidiosis increased relative to the risk at the mean 12-months prior (MC-MAXT 10th percentile: RR 1.132; 95% CI 1.085 – 1.132; MC-MINT 10th percentile: RR 1.054; 95% CI 1.010 – 1.100).

Although the risk of cryptosporidiosis was not significantly different (p<0.05) when comparing extreme high MC-MAXT or MC-MINT (99th percentile) to the mean, when either MC-MAXT or MC-MINT was high or very high (90th and 95th percentiles) 2-6 months prior, there was a significant increase in cryptosporidiosis risk relative to the risk at the mean (Table 2).

#### Precipitation

Very high and extreme high (95th and 99th percentile values) MC-PREC in the 2-6 months prior were associated with a significant decrease in the relative risk of cryptosporidiosis (Figure 3, Table 2). For example, when MC-PREC 4-months prior was at the 99th percentile value (84 mm higher than the county and calendar-month average), the relative risk of cryptosporidiosis was 19.2% lower than the risk at the average (RR 0.808; 95% CI 0.749 – 0.872).

When MC-PREC was lower than average during the 2-6 months prior, cryptosporidiosis risk was not significantly different than the risk at the mean (Table 2). However, when a 12-month lag was used, lower than average MC-PREC (1st – 10th percentiles) was associated with a decrease in cryptosporidiosis risk as compared to the risk at the mean (e.g., MC-PREC 1st percentile: RR 0.727; 95% CI 0.659 – 0.802).

### Giardiasis

#### Temperature

Our models highlighted a negative linear relationship between giardiasis and temperature, wherein the risk of giardiasis was highest when 2-6 month lagged temperatures were lower than average, and lowest when 2-6 month lagged temperatures were higher than average (Figure 4). When MC-MAXT was low, very low or extremely low relative to the average, there was a 1-2% increase in the 2-6 month lagged risk of giardiasis (Table 3). For MC-MINT, the change in risk was more pronounced: there was a 2-5% increase in the risk of giardiasis at low to extreme low values of MC-MINT, relative to the risk at the county and calendar-month mean MC-MINT (Table 3). Similarly, when MC-MAXT and MC-MINT were high to extremely high (90th – 99th percentiles), there was a 1-2% increase in the 2-6 month lagged risk of giardiasis relative to mean MC-MAXT, and a 2-4% increase in the 2-6 month lagged risk of giardiasis relative to mean MC-MINT. The risk of giardiasis at high and low temperature extremes was not statistically significant (p < 0.05) when a 12-month lag was used.

![Figure 4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/09/01/2023.08.31.23294911/F4.medium.gif)

[Figure 4.](http://medrxiv.org/content/early/2023/09/01/2023.08.31.23294911/F4)

Figure 4. 3D depictions of the lagged effects of temperature and precipitation on giardiasis
3D illustrations of the best performing model for the lagged relationship between giardiasis cases and MC-MAXT, MC-MINT and MC-PREC. For each predictor-outcome pair, the 3D graphic depicts the change in the relative risk of giardiasis across the range of observed temperature and precipitation values for lags of 0 – 12 months, relative to the risk found at the county and calendar-month averages for temperature and precipitation.

View this table:
[Table 3.](http://medrxiv.org/content/early/2023/09/01/2023.08.31.23294911/T3)

Table 3. The relative risk of giardiasis in Colorado counties between 1997 – 2017, comparing temperature and precipitation extremes to county averages.

#### Precipitation

Neither the extreme high nor extreme low values of 2-6 month lagged MC-PREC was significantly associated with a change in giardiasis risk, relative to the county and calendar-month mean values (Table 3). However, there was consistent evidence of 12-month lagged effects for giardiasis. We found that higher than average MC-PREC 12-months prior (at the 90th, 95th and 99th percentile values) was associated with a decrease in giardiasis risk, relative to the risk at the mean (90th percentile: RR= 1.020, 95% CI 1.007 – 1.034; 95th percentile: RR= 1.030, 95% CI 1.010 – 1.051; 99th percentile: RR= 1.056, 95% CI 1.019 – 1.094). Similarly, lower than average MC-PREC 12-months prior (at the 1st, 5th and 10th percentile values) was associated with a decrease in giardiasis risk, relative to the risk at the mean (1st percentile: RR= 0.968, 95% CI 0.947 – 0.989; 5th percentile: RR= 0.977, 95% CI 0.962 – 0.992; 10th percentile: RR= 0.982, 95% CI 0.970 – 0.994).

## Discussion

In this analysis, we found that cryptosporidiosis and giardiasis responded similarly to unseasonably low temperatures in Colorado counties between 1997 and 2017. There was a significant increase in the relative risk of both cryptosporidiosis and giardiasis infection when temperatures 4-6 months prior were extremely low for a given county and calendar-month (i.e. at the 1st percentile value, corresponding with -4.5°C below average for MC-MAXT, and -3.8°C below average for MC-MINT). As far as we are aware, this is a novel finding that has not been replicated in other epidemiological studies, perhaps because both giardiasis and cryptosporidiosis incidence tend to peak in summer or late summer in the United States, 45 leading most investigations to focus primarily on warmer temperatures and transmission potential. Nevertheless, experimental studies have highlighted that both *Cryptosporidium* and *Giardia* can withstand cool temperatures, with the viability of (oo)cysts improving when stored at 4°C as compared to room temperature. 62 One study identified developmental stage parasites in mice exposed to *Cryptosporidium* oocysts and stored at -20°C, -15°C, -10°C and 5°C for between 1-168 hours before thawing, 50 highlighting oocyst persistence even when temperatures are lower that average.

By contrast, when temperatures were moderately high (90th – 95th percentiles) 2-6 months prior, we found a statistically significant decrease in giardiasis cases in Colorado counties, while there was a significant increase in cryptosporidiosis cases. For giardiasis that relationship held at the 99th percentile values of MC-MAXT and MC-MINT, while there was no evidence of an association between extreme high temperatures (99th percentile) and cryptosporidiosis cases. The fact that higher than average temperatures was associated with a decrease in giardiasis cases at all levels (90, 95th and 99th percentiles) but an increase in cryptosporidiosis at moderate levels (90th and 95th percentiles), may be related to the two organisms having different thresholds for tolerating above average temperatures. Both in vitro and in vivo experiments have suggested that *Giardia* cysts demonstrate greater temperature-dependence and higher die-off rates than *Cryptosporidium* oocysts at both high and low temperatures. 29,41 For example, one recent study focused on the potential for foodborne transmission of giardiasis and cryptosporidiosis found that *Cryptosporidium* oocysts were stable at both room temperature and refrigeration conditions (4°C), whereas *Giardia* cysts saw a 50% die-off rate at room temperature over a period of 24-hours. 62 Thus, the fact that both moderate and extreme increases in temperature relative to the county and calendar-month mean were associated with a decrease in giardiasis risk may potentially reflect the greater temperature sensitivity of *Giardia*. Meanwhile, the increase in cryptosporidiosis cases with moderately high (90th – 95th percentiles), but not extremely high temperatures (99th percentiles) could be indicative of a threshold effect, wherein moderately high temperatures lengthen the transmission season and/or allow for transmission via alternative vectors, while extremely high temperatures lead to exponential oocysts die-off and reduced transmission. 40

Another possible explanation for why *Cryptosporidium* and *Giardia* respond differently to higher-than-average temperatures is that the primary hosts associated with each pathogen are different. For example, a study conducted in Georgia, USA found that synanthropic flies living alongside wildlife and livestock could serve as mechanical vectors of *Cryptosporidium* and *Giardia*, with viable *Cryptosporidium* being isolated from 56% of sampled flies, while *Giardia* was isolated from just 8% of those sampled. 42 As warm-weather, ectothermic species, the development time of the studied fly species tends to decrease as temperatures increase, up until a threshold temperature when development will cease – for example, above 34°C for similar species of blowflies found in Asia. 63 Thus, if cattle and synanthropic flies were in fact a major part of the transmission pathway for cryptosporidiosis, but not giardiasis (as has been found in other contexts 42,64), then the association between cryptosporidiosis and moderate, but not extreme increases in maximum and minimum temperature is consistent with what would be expected. However, further study is still needed to determine whether flies are in fact an important vector in either the cryptosporidiosis or giardiasis transmission cycles in Colorado.

In our study, there was a significant decrease in the relative risk of cryptosporidiosis when total monthly precipitation in the 2-6 months prior was relatively high. However, neither relatively high nor low monthly precipitation were found to significantly impact giardiasis using a 2-, 4- or 6-month lag. One possible explanation for this difference could be that the primary underlying exposure pathways of *Cryptosporidium* and *Giardia* transmission differ between the two pathogens. For example, if surface water contamination via zoonotic livestock or wildlife species were the primary source of cryptosporidiosis transmission, but not giardiasis, then the overall concentration of cryptosporidium in surface water sources might be diluted following higher than average precipitation. As an example, this was believed to be the case in one study conducted in Canada which observed lower odds of human cryptosporidiosis when water levels were higher than average 19-20 days prior. 64 While the lag periods assessed in our study were of a considerably longer duration than those studied by Brankston and colleagues, we did find a similarly negative association between cryptosporidiosis and the 90th, 95th and 99th percentiles of total monthly precipitation at both 2-month and 4-month lags, a finding that held for the 95th and 99th percentile at the 6-month lag. Meanwhile only the 99th percentile value of precipitation was still negatively associated with cryptosporidiosis when lagged 12-months. Notably, a recent review that assessed the body of evidence in support of the concentration-dilution hypothesis – that is, that conflicting findings on the impacts of rainfall on diarrhea result from underlying differences in background rain levels – found four studies that identified dilution (specifically, rainfall following wet periods) as a potential mechanism explaining an inverse association between rainfall and diarrhea. 66 While an investigation of the impacts of rainfall following wetter than average weather was outside of the scope of this study, future investigations of the impacts of this phenomenon on giardiasis and cryptosporidiosis, as well as the underlying sources of contamination for each pathogen within Colorado counties is warranted.

Whereas other studies assessing the effects of weather and climate on cryptosporidiosis and giardiasis have most frequently assessed lag periods of between one and three months, 35–38,44 in our study, we allowed for a lag of up to 12-months. This decision was based on the body of evidence from a small set of experimental studies which have suggested that under the right conditions, *Cryptosporidium* oocysts and *Giardia* cysts may be able to survive and maintain viability for much longer than the period typically assessed in the epidemiological literature – potentially as long as six months to a year. 51,52 Whereas the relative risk of maximum monthly temperature and minimum monthly temperature extremes relative to the averages tended to only be statistically significant for giardiasis and (to a slightly lesser extent) cryptosporidiosis in the 2-6 month lag range, the 12-month lagged estimates of the relative risk of precipitation extremes (compared to the mean) were statistically significant in all but one instance (See Tables 2 & 3). Moreover, lower than average precipitation relative to the mean was only significantly associated with giardiasis or cryptosporidiosis when lagged 12-months, suggesting that the transmission of these pathogens in Colorado counties may be sensitive to longer-term periods of excess dryness or drought, as has been found with diarrheal pathogens in other regions of the world. 23,24,67 While our study was not designed to assess the impacts of rainfall after prolonged periods of drought on transmission, future studies aimed at investigating this potential in Colorado and the American West at large are greatly needed.

Although giardiasis and cryptosporidiosis cases reported in Colorado counties between 1997 – 2017 tended to follow the standard seasonal pattern that have been previously identified for temperate climates, wherein case counts tend to peak in the summer or early fall 45 (See Table 1), the final models selected via cross-validation for each of our outcome-predictor pairs did not ultimately include an interaction for season. This is notable, as season-specific effects have been demonstrated in other contexts, 31,68 which are believed to arise as a result of season-specific exposure opportunities such as local agricultural practices or recreational water uses. 68,69 While it is possible that the relationship between precipitation and temperature and cryptosporidiosis and giardiasis truly does not vary substantially by season in Colorado counties, it is also possible that our models were underpowered for detecting seasonal effects.

Thus, a limitation of this study was the potential for being underpowered when it came to the seasonal effect modification assessment, due to the panel structure of the data, our conditioning on county-month, the complex lag-predictor cross-basis structures already included in each model, and the fact that we had a relatively small number of observations for each county (252 month-years for 63 counties, and 216 for Broomfield County). Another limitation to this study was that the smallest time-unit possible given the data available was month, making it impossible for us to determine the effects of short-term weather events. Additionally, because county was the smallest geographic-unit available in this study, we were unable to investigate location-specific nuances such as variations in elevation, access to healthcare and local population factors such as popular recreational activities. While including a grouping term in our cross-basis and conditioning on county-calendar-month helps to account for unmeasured intra-county variation, we remain unable to assess inner-county differences. Another limitation to this study is that the CEDRS database that was used for obtaining monthly case counts of cryptosporidiosis and giardiasis by county does not necessarily record the location where the pathogen was acquired, only the county of residence for an infected individual.

This study relied on clinically reported cryptosporiosis and giardiasis, and, due to the variable severity of disease, many infections with these protozoa go unreported. Meanwhile, among cases of sufficient severity to warrant medical attention, changes in the use, accuracy and types of diagnostic tests that were available for the clinical detection of cryptosporidiosis and giardiasis across the study period also have the potential to bias our results. While cryptosporidiosis and giardiasis diagnostic tests remained relatively unchanged in the 90s and early 2000s, 70 the more widespread introduction of multiplex tests means that healthcare providers could more easily test for an array of diseases in patients exhibiting gastrointestinal symptoms. Being able to easily test for multiple diseases could increase detection of suspected disease, as well as increase incidental findings of cryptosporidiosis and giardiasis. In Figure 2, a plot of the total case counts by year across all Colorado counties shows that giardiasis cases appear to have a slight downward trend over time, whereas cryptosporidiosis cases have shown a steady increase across the study period. Thus, in the case of cryptosporidiosis, it is possible that increasing use of multiplex testing is influencing the upward trend in case reports. A linear term for year was included in each of our models to help account for these linear time trends.

Our analysis has shown that the risk of cryptosporidiosis and giardiasis in Colorado counties is likely to change as weather patterns in Colorado shift in the coming years. Warmer temperatures and an increase in weather extremes (i.e., blizzards, droughts, etc.) can be expected. In the last 30 years, Colorado’s average temperature has increased by 2°F (3.6 °C), with current estimates projecting that the average temperature for the state could increase another 5°F (9°C) before 2050. 72 As such, it is possible that the coming years will see a continued reduction in the transmission of giardiasis (Figure 2), as our study has highlighted that extreme high (i.e., 99th percentile) values of maximum and minimum temperatures tend to be associated with reduced giardiasis infection risk (relative to the risk at average values) in Colorado counties. However, it is possible that rising temperatures will increase cryptosporidium risk, as we found temperature extremes at the 90th and 95th percentile values were associated with a statistically significant increase in cryptosporidiosis infection risk. These discordant findings indicate that the impact of climate change may not be the same for all protozoan pathogens and highlight the need for further study of the locations-specific environmental conditions that could be simultaneously promoting the transmission potential of one protozoan pathogen while suppressing another. Further investigation into pathogen-specific environmental contamination sources in Colorado counties, as well as the potential impacts of extreme weather following prolonged periods of wet weather or prolonged periods of drought are also needed. Colorado-specific climate projections are also essential pre-requisites for quantifying the associated increased risk of infectious disease transmission over the coming decades. Ultimately, combining the findings of future climate projections and environmental risk assessments with our findings on the Colorado-specific impacts of extreme temperatures and precipitation on two important, but often overlooked parasitic diseases will help to inform the development of education and prevention messaging for high-risk populations, as well as aid in the development of tailored climate change adaptation and mitigation protocols.

## Supporting information

Supplemental Table S1 [[supplements/294911_file06.docx]](pending:yes)

Supplemental Table S2 [[supplements/294911_file07.docx]](pending:yes)

Supplemental Table S3 [[supplements/294911_file08.docx]](pending:yes)

Supplemental Table S4 [[supplements/294911_file09.docx]](pending:yes)

## Data Availability

Data on reportable diseases in Colorado were collected using the Colorado Electronic Disease Reporting System, which is made available upon request from the Colorado Department of Public Health and Environment. The Daily Surface Weather and Climatological Summaries (Daymet) Version 3 dataset that was used in this analysis is available upon request from the National Center for Atmospheric Research.

## Data Sharing

Data on reportable diseases in Colorado were collected using the Colorado Electronic Disease Reporting System, which is made available upon request from the Colorado Department of Public Health and Environment. The Daily Surface Weather and Climatological Summaries (Daymet) Version 3 dataset that was used in this analysis is available upon request from the National Center for Atmospheric Research.

## Declaration of competing financial interests

The authors declare that they have no actual or potential competing financial interests. There is no grant support to disclose.

## Acknowledgements

This work was made possible by the Colorado Department of Public Health and Environment and the National Center for Atmospheric Research, who collected and shared the data used in this analysis.

*   Received August 31, 2023.
*   Revision received August 31, 2023.
*   Accepted September 1, 2023.


*   © 2023, Posted by Cold Spring Harbor Laboratory

This pre-print is available under a Creative Commons License (Attribution 4.0 International), CC BY 4.0, as described at [http://creativecommons.org/licenses/by/4.0/](http://creativecommons.org/licenses/by/4.0/)

## REFERENCES

1.  1.1.  Masson-Delmotte, V., 
    2.   P. Zhai, 
    3.   A. Pirani, 
    4.   S.L. Connors, 
    5.   C. Péan, 
    6.   S. Berger, 
    7.   N. Caud, 
    8.   Y. Chen, 
    9.   L. Goldfarb, 
    10.  M.I. Gomis, 
    11.  M. Huang, 
    12.  K. Leitzell, 
    13.  E. Lonnoy, 
    14.  J.B.R. Matthews, 
    15.  T.K. Maycock, 
    16.  T. Waterfield, 
    17.  O. Yelekçi, 
    18.  R. Yu, and 
    19.  B. Zhou
    
    IPCC. Climate Change 2021: The Physical Science Basis. Contribution of Working Group I to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change [Masson-Delmotte, V.,  P. Zhai,  A. Pirani,  S.L. Connors,  C. Péan,  S. Berger,  N. Caud,  Y. Chen,  L. Goldfarb,  M.I. Gomis,  M. Huang,  K. Leitzell,  E. Lonnoy,  J.B.R. Matthews,  T.K. Maycock,  T. Waterfield,  O. Yelekçi,  R. Yu, and  B. Zhou (eds.)]. Cambridge University Press. In Press., 2021.
    
    

2.  2.Luber G, Lemery J. Global Climate Change and Human Health: From Science to Practice: Wiley; 2015.
    
    

3.  3.Pozio E. How globalization and climate change could affect foodborne parasites. Experimental parasitology 2020; 208: 107807.
    
    

4.  4.Fletcher SM, Stark D, Harkness J, Ellis J. Enteric protozoa in the developed world: a public health perspective. Clinical microbiology reviews 2012; 25(3): 420–49.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiY21yIjtzOjU6InJlc2lkIjtzOjg6IjI1LzMvNDIwIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjMvMDkvMDEvMjAyMy4wOC4zMS4yMzI5NDkxMS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

5.  5.Intergovernmental Panel on Climate Change. IPCC Special Report on Global Warming of 1.5 °C: Summary for Policymakers (SR15). Geneva, Switzerland, 2018.
    
    

6.  6.1.  Akhtar R
    
    Akhtar R. Introduction: Extreme Weather Events and Human Health: A Global Perspective. In: Akhtar R, ed. Extreme Weather Events and Human Health: International Case Studies. Cham: Springer International Publishing; 2020: 3-11.
    
    

7.  7.1.  Quah SR
    
    Molloy SL, Dreelin EA, Rose JB. Extreme Weather Events and Human Health. In: Quah SR, ed. International Encyclopedia of Public Health (Second Edition). Oxford: Academic Press; 2017: 70–6.
    
    

8.  8.Efstratiou A, Ongerth JE, Karanis P. Waterborne transmission of protozoan parasites: Review of worldwide outbreaks - An update 2011–2016. Water Research 2017; 114: 14–22.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.watres.2017.01.036&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28214721&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F09%2F01%2F2023.08.31.23294911.atom) 

9.  9.Lane S, Lloyd D. Current trends in research into the waterborne parasite Giardia. Crit Rev Microbiol 2002; 28(2): 123–47.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1080/1040-840291046713&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12109771&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F09%2F01%2F2023.08.31.23294911.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000177038300002&link_type=ISI) 

10. 10.Chalmers RM, Davies AP. Minireview: clinical cryptosporidiosis. Experimental parasitology 2010; 124(1): 138–46.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.exppara.2009.02.003&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19545516&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F09%2F01%2F2023.08.31.23294911.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000274134700017&link_type=ISI) 

11. 11.Liu J, Platts-Mills JA, Juma J, et al. Use of quantitative molecular diagnostic methods to identify causes of diarrhoea in children: a reanalysis of the GEMS case-control study. *Lancet (London*, England*)* 2016; 388(10051): 1291–301.
    
    

12. 12.Fayer R, Dubey JP, Lindsay DS. Zoonotic protozoa: from land to sea. Trends in parasitology 2004; 20(11): 531–6.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.pt.2004.08.008&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15471705&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F09%2F01%2F2023.08.31.23294911.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000224951400008&link_type=ISI) 

13. 13.Torgerson PR, Devleesschauwer B, Praet N, et al. World Health Organization Estimates of the Global and Regional Disease Burden of 11 Foodborne Parasitic Diseases, 2010: A Data Synthesis. PLoS medicine 2015; 12(12): e1001920.
    
    

14. 14.Robertson LJ. The potential for marine bivalve shellfish to act as transmission vehicles for outbreaks of protozoan infections in humans: A review. International Journal of Food Microbiology 2007; 120(3): 201–16.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ijfoodmicro.2007.07.058&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17928081&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F09%2F01%2F2023.08.31.23294911.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000252062800001&link_type=ISI) 

15. 15.Karanis P, Kourenti C, Smith H. Waterborne transmission of protozoan parasites: A worldwide review of outbreaks and lessons learnt. Journal of water and health 2006; 5(1): 1–38.
    
    

16. 16.Rose JB, Slifko TR. Giardia, Cryptosporidium, and Cyclospora and their impact on foods: a review. J Food Prot 1999; 62(9): 1059–70.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=10492484&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F09%2F01%2F2023.08.31.23294911.atom) 

17. 17.Dumètre A, Aubert D, Puech P-H, Hohweyer J, Azas N, Villena I. Interaction Forces Drive the Environmental Transmission of Pathogenic Protozoa. Applied and Environmental Microbiology 2012; 78(4): 905–12.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiYWVtIjtzOjU6InJlc2lkIjtzOjg6Ijc4LzQvOTA1IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjMvMDkvMDEvMjAyMy4wOC4zMS4yMzI5NDkxMS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

18. 18.Ahmed SA, Guerrero Flórez M, Karanis P. The impact of water crises and climate changes on the transmission of protozoan parasites in Africa. Pathogens and global health 2018; 112(6): 281–93.
    
    

19. 19.Nichols G, Lane C, Asgari N, Verlander NQ, Charlett A. Rainfall and outbreaks of drinking water related disease and in England and Wales. Journal of water and health 2009; 7(1): 1–8.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoicHBpd2Fqd2giO3M6NToicmVzaWQiO3M6NToiNy8xLzEiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMy8wOS8wMS8yMDIzLjA4LjMxLjIzMjk0OTExLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

20. 20.Thomas KM, Charron DF, Waltner-Toews D, Schuster C, Maarouf AR, Holt JD. A role of high impact weather events in waterborne disease outbreaks in Canada, 1975 - 2001. Int J Environ Health Res 2006; 16(3): 167–80.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1080/09603120600641326&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16611562&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F09%2F01%2F2023.08.31.23294911.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000236763800001&link_type=ISI) 

21. 21.Curriero FC, Patz JA, Rose JB, Lele S. The association between extreme precipitation and waterborne disease outbreaks in the United States, 1948-1994. Am J Public Health 2001; 91(8): 1194–9.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.2105/AJPH.91.8.1194&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11499103&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F09%2F01%2F2023.08.31.23294911.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000170345600017&link_type=ISI) 

22. 22.Schreiber C, Heinkel S-B, Zacharias N, et al. Infectious rain? Evaluation of human pathogen concentrations in stormwater in separate sewer systems. Water Science and Technology 2019; 80(6): 1022–30.
    
    

23. 23.Chhetri BK, Takaro TK, Balshaw R, et al. Associations between extreme precipitation and acute gastro-intestinal illness due to cryptosporidiosis and giardiasis in an urban Canadian drinking water system (1997-2009). Journal of water and health 2017; 15(6): 898–907.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoicHBpd2Fqd2giO3M6NToicmVzaWQiO3M6ODoiMTUvNi84OTgiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMy8wOS8wMS8yMDIzLjA4LjMxLjIzMjk0OTExLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

24. 24.Lal A, Konings P. Beyond reasonable drought: hotspots reveal a link between the ‘Big Dry’ and cryptosporidiosis in Australia’s Murray Darling Basin. Journal of water and health 2018; 16(6): 1033–7.
    
    

25. 25.Senhorst HA, Zwolsman JJ. Climate change and effects on water quality: a first impression. Water Sci Technol 2005; 51(5): 53–9.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoicHBpd2F3c3QiO3M6NToicmVzaWQiO3M6NzoiNDQvNS81MyI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIzLzA5LzAxLzIwMjMuMDguMzEuMjMyOTQ5MTEuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

26. 26.Zwolsman JJ, van Bokhoven AJ. Impact of summer droughts on water quality of the Rhine River - a preview of climate change? Water Sci Technol 2007; 56(4): 45–55.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoicHBpd2F3c3QiO3M6NToicmVzaWQiO3M6ODoiNTYvMTAvNDUiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMy8wOS8wMS8yMDIzLjA4LjMxLjIzMjk0OTExLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

27. 27.Semenza JC, Herbst S, Rechenburg A, et al. Climate Change Impact Assessment of Food- and Waterborne Diseases. Critical reviews in environmental science and technology 2012; 42(8): 857–90.
    
    

28. 28.Xiao L, Fayer R. Molecular characterisation of species and genotypes of Cryptosporidium and Giardia and assessment of zoonotic transmission. International journal for parasitology 2008; 38(11): 1239–55.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ijpara.2008.03.006&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18479685&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F09%2F01%2F2023.08.31.23294911.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000259020700004&link_type=ISI) 

29. 29.Olson ME, Goh J, Phillips M, Guselle N, McAllister TA. Giardia Cyst and Cryptosporidium Oocyst Survival in Water, Soil, and Cattle Feces. Journal of environmental quality 1999; 28: 1991–6.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.2134/jeq1999.00472425002800060040x&link_type=DOI) 

30. 30.Daniels ME, Smith WA, Schmidt WP, Clasen T, Jenkins MW. Modeling Cryptosporidium and Giardia in Ground and Surface Water Sources in Rural India: Associations with Latrines, Livestock, Damaged Wells, and Rainfall Patterns. Environmental science & technology 2016; 50(14): 7498–507.
    
    

31. 31.Muchiri JM, Ascolillo L, Mugambi M, et al. Seasonality of Cryptosporidium oocyst detection in surface waters of Meru, Kenya as determined by two isolation methods followed by PCR. Journal of water and health 2009; 7(1): 67–75.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoicHBpd2Fqd2giO3M6NToicmVzaWQiO3M6NjoiNy8xLzY3IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjMvMDkvMDEvMjAyMy4wOC4zMS4yMzI5NDkxMS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

32. 32.Omarova A, Tussupova K, Berndtsson R, Kalishev M, Sharapatova K. Protozoan Parasites in Drinking Water: A System Approach for Improved Water, Sanitation and Hygiene in Developing Countries. International journal of environmental research and public health 2018; 15(3).
    
    

33. 33.Carlton EJ, Woster AP, DeWitt P, Goldstein RS, Levy K. A systematic review and meta-analysis of ambient temperature and diarrhoeal diseases. International journal of epidemiology 2016; 45(1): 117–30.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/ije/dyv296&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26567313&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F09%2F01%2F2023.08.31.23294911.atom) 

34. 34.Levy K, Woster AP, Goldstein RS, Carlton EJ. Untangling the Impacts of Climate Change on Waterborne Diseases: a Systematic Review of Relationships between Diarrheal Diseases and Temperature, Rainfall, Flooding, and Drought. Environmental science & technology 2016; 50(10): 4905–22.
    
    

35. 35.Hu W, Mengersen K, Tong S. Risk factor analysis and spatiotemporal CART model of cryptosporidiosis in Queensland, Australia. BMC infectious diseases 2010; 10(1): 311.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21029426&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F09%2F01%2F2023.08.31.23294911.atom) 

36. 36.Hu W, Mengersen K, Fu S-Y, Tong S. The use of ZIP and CART to model cryptosporidiosis in relation to climatic variables. International Journal of Biometeorology 2010; 54(4): 433–40.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s00484-009-0294-4&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20084521&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F09%2F01%2F2023.08.31.23294911.atom) 

37. 37.Hu W, Tong S, Mengersen K, Connell D. Weather Variability and the Incidence of Cryptosporidiosis: Comparison of Time Series Poisson Regression and SARIMA Models. Annals of epidemiology 2007; 17(9): 679–88.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.annepidem.2007.03.020&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17604645&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F09%2F01%2F2023.08.31.23294911.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000249293100005&link_type=ISI) 

38. 38.Naumova EN, Jagai JS, Matyas B, DeMaria A, Jr., MacNeill IB, Griffiths JK. Seasonality in six enterically transmitted diseases and ambient temperature. Epidemiology and infection 2007; 135(2): 281–92.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1017/S0950268806006698&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17291363&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F09%2F01%2F2023.08.31.23294911.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000244652800013&link_type=ISI) 

39. 39.Ikiroma IA, Pollock KG. Influence of weather and climate on cryptosporidiosis-A review. Zoonoses Public Health 2021; 68(4): 285–98.
    
    

40. 40.Lal A, Baker MG, Hales S, French NP. Potential effects of global environmental changes on cryptosporidiosis and giardiasis transmission. Trends in parasitology 2013; 29(2): 83–90.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.pt.2012.10.005&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23219188&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F09%2F01%2F2023.08.31.23294911.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000315180200005&link_type=ISI) 

41. 41.Alum A, Absar IM, Asaad H, Rubino JR, Ijaz MK. Impact of environmental conditions on the survival of cryptosporidium and giardia on environmental surfaces. Interdiscip Perspect Infect Dis 2014; 2014: 210385.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1155/2014/210385&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25045350&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F09%2F01%2F2023.08.31.23294911.atom) 

42. 42.Conn DB, Weaver J, Tamang L, Graczyk TK. Synanthropic flies as vectors of Cryptosporidium and Giardia among livestock and wildlife in a multispecies agricultural complex. *Vector borne and zoonotic diseases (Larchmont*, NY*)* 2007; 7(4): 643–51.
    
    

43. 43.King BJ, Monis PT. Critical processes affecting Cryptosporidium oocyst survival in the environment. Parasitology 2007; 134(Pt 3): 309–23.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1017/S0031182006001491&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17096874&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F09%2F01%2F2023.08.31.23294911.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000244657800001&link_type=ISI) 

44. 44.Kent L, McPherson M, Higgins N. A positive association between cryptosporidiosis notifications and ambient temperature, Victoria, Australia, 2001-2009. Journal of water and health 2015; 13(4): 1039–47.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoicHBpd2Fqd2giO3M6NToicmVzaWQiO3M6OToiMTMvNC8xMDM5IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjMvMDkvMDEvMjAyMy4wOC4zMS4yMzI5NDkxMS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

45. 45.Lal A, Hales S, French N, Baker MG. Seasonality in human zoonotic enteric diseases: a systematic review. PloS one 2012; 7(4): e31883.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0031883&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22485127&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F09%2F01%2F2023.08.31.23294911.atom) 

46. 46.Jagai JS, Castronovo DA, Monchak J, Naumova EN. Seasonality of cryptosporidiosis: A meta-analysis approach. Environmental research 2009; 109(4): 465–78.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19328462&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F09%2F01%2F2023.08.31.23294911.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000265897700012&link_type=ISI) 

47. 47.deRegnier DP, Cole L, Schupp DG, Erlandsen SL. Viability of Giardia cysts suspended in lake, river, and tap water. Appl Environ Microbiol 1989; 55(5): 1223–9.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiYWVtIjtzOjU6InJlc2lkIjtzOjk6IjU1LzUvMTIyMyI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIzLzA5LzAxLzIwMjMuMDguMzEuMjMyOTQ5MTEuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

48. 48.Petersen HH, Enemark HL, Olsen A, Amin MG, Dalsgaard A. Transport of Cryptosporidium parvum oocysts in soil columns following applications of raw and separated liquid slurries. Appl Environ Microbiol 2012; 78(17): 5994–6000.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiYWVtIjtzOjU6InJlc2lkIjtzOjEwOiI3OC8xNy81OTk0IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjMvMDkvMDEvMjAyMy4wOC4zMS4yMzI5NDkxMS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

49. 49.Robertson LJ, Campbell AT, Smith HV. Survival of Cryptosporidium parvum oocysts under various environmental pressures. Appl Environ Microbiol 1992; 58(11): 3494–500.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiYWVtIjtzOjU6InJlc2lkIjtzOjEwOiI1OC8xMS8zNDk0IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjMvMDkvMDEvMjAyMy4wOC4zMS4yMzI5NDkxMS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

50. 50.Fayer R, Nerad T. Effects of low temperatures on viability of Cryptosporidium parvum oocysts. Appl Environ Microbiol 1996; 62(4): 1431–3.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiYWVtIjtzOjU6InJlc2lkIjtzOjk6IjYyLzQvMTQzMSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIzLzA5LzAxLzIwMjMuMDguMzEuMjMyOTQ5MTEuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

51. 51.Fayer R, Trout JM, Jenkins MC. Infectivity of Cryptosporidium parvum Oocysts Stored in Water at Environmental Temperatures. The Journal of parasitology 1998; 84(6): 1165–9.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.2307/3284666&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=9920307&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F09%2F01%2F2023.08.31.23294911.atom) 

52. 52.Badenoch J. Cryptosporidium in water supplies. Report of the group of experts: HMSO; 1990.
    
    

53. 53.Gasparrini A, Armstrong B, Kenward MG. Distributed lag non-linear models. Statistics in medicine 2010; 29(21): 2224–34.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/sim.3940&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20812303&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F09%2F01%2F2023.08.31.23294911.atom) 

54. 54.Armstrong BG, Gasparrini A, Tobias A. Conditional Poisson models: a flexible alternative to conditional logistic case cross-over analysis. BMC medical research methodology 2014; 14(1): 122.
    
    

55. 55.Colorado Department of Local Affairs. County and Municipal Population Timeseries. 2019. [https://demography.dola.colorado.gov/population/data/county-muni-timeseries/](https://demography.dola.colorado.gov/population/data/county-muni-timeseries/) (accessed March 1, 2019.
    
    

56. 56.Thornton PE, Thornton MM, Mayer BW, et al. Daymet: Daily Surface Weather Data on a 1-km Grid for North America, Version 3. ORNL Distributed Active Archive Center; 2017.
    
    

57. 57.Colorado Rural Health Center. Colorado: County Designations, 2016. 2016. [https://www.colorado.gov/pacific/sites/default/files/PCO\_CHSC\_CountyDesignations\_2016.pdf](https://www.colorado.gov/pacific/sites/default/files/PCO_CHSC_CountyDesignations_2016.pdf) (accessed March, 29 2019).
    
    

58. 58.Gasparrini A. Distributed Lag Linear and Non-Linear Models in R: The Package dlnm. Journal of Statistical Software 2011; 43(8): 1–20.
    
    

59. 59.RStudio Team. RStudio: Integrated Development for R. RStudio, PBC, Boston, MA, 2020. 2020.
    
    

60. 60.Kuhn M. Building Predictive Models in R Using the caret Package. Journal of Statistical Software 2008; 28(5): 1–26.
    
    

61. 61.StataCorp. Stata Statistical Software: Release 15. College Station, TX: StataCorp LP; 2015.
    
    

62. 62.Utaaker KS, Skjerve E, Robertson LJ. Keeping it cool: Survival of Giardia cysts and Cryptosporidium oocysts on lettuce leaves. Int J Food Microbiol 2017; 255: 51–7.
    
    

63. 63.Zhang Y, Wang Y, Sun J, et al. Temperature-dependent development of the blow fly Chrysomya pinguis and its significance in estimating postmortem interval. Royal Society Open Science 2019; 6(9): 190003.
    
    

64. 64.Brankston G, Boughen C, Ng V, Fisman DN, Sargeant JM, Greer AL. Assessing the impact of environmental exposures and Cryptosporidium infection in cattle on human incidence of cryptosporidiosis in Southwestern Ontario, Canada. PloS one 2018; 13(4): e0196573–e.
    
    

65. 65.Tolouei S, Dewey R, Snodgrass WJ, et al. Assessing microbial risk through event-based pathogen loading and hydrodynamic modelling. The Science of the total environment 2019; 693: 133567.
    
    

66. 66.Kraay ANM, Man O, Levy MC, Levy K, Ionides E, Eisenberg JNS. Understanding the Impact of Rainfall on Diarrhea: Testing the Concentration-Dilution Hypothesis Using a Systematic Review and Meta-Analysis. Environmental health perspectives 2020; 128(12): 126001.
    
    

67. 67.Bush KF, O’Neill MS, Li S, et al. Associations between extreme precipitation and gastrointestinal-related hospital admissions in Chennai, India. Environmental health perspectives 2014; 122(3): 249–54.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24345350&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F09%2F01%2F2023.08.31.23294911.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000333670500014&link_type=ISI) 

68. 68.Odoi A, Martin SW, Michel P, Holt J, Middleton D, Wilson J. Geographical and temporal distribution of human giardiasis in Ontario, Canada. International journal of health geographics 2003; 2(1): 5.
    
    

69. 69.Learmonth J, Ionas G, Pita A, Cowie R. Seasonal shift in Cryptosporidium parvum transmission cycles in New Zealand. Journal of Eukaryotic Microbiology 2002; (Supplement): 34S–5S.
    
    

70. 70.Ricciardi A, Ndao M. Diagnosis of parasitic infections: what’s going on? Journal of biomolecular screening 2015; 20(1): 6–21.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1177/1087057114548065&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25170017&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F09%2F01%2F2023.08.31.23294911.atom) 

71. 71.Plowright RK, Sokolow SH, Gorman ME, Daszak P, Foley JE. Causal inference in disease ecology: investigating ecological drivers of disease emergence. Frontiers in Ecology and the Environment 2008; 6(8): 420–9.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1890/070086&link_type=DOI) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000260158300018&link_type=ISI) 

72. 72.Colorado Health Institute. Colorado’s Climate and Colorado’s Health: Examining the Connection: CHI, 2017.
    
    

73. 73.1.  Motarjemi Y
    
    Dixon BR. Protozoa: Cryptosporidium spp. In: Motarjemi Y, ed. Encyclopedia of Food Safety. Waltham: Academic Press; 2014: 18-22.
    
    

74. 74.Chaidez C, Soto M, Gortares P, Mena K. Occurrence of Cryptosporidium and Giardia in irrigation water and its impact on the fresh produce industry. Int J Environ Health Res 2005; 15(5): 339–45.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1080/09603120500289010&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16416751&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F09%2F01%2F2023.08.31.23294911.atom) 

75. 75.Davies-Colley R, Nagels J, Lydiard E. Stormflow-dominated loads of faecal pollution from an intensively dairy-farmed catchment. Water Sci Technol 2008; 57(10): 1519–23.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoicHBpd2F3c3QiO3M6NToicmVzaWQiO3M6MTA6IjU3LzEwLzE1MTkiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMy8wOS8wMS8yMDIzLjA4LjMxLjIzMjk0OTExLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

76. 76.Ajeagah GA, Njine T, Bilong Bilong C, et al. Seasonal distribution of enteric opportunistic cryptosporidium spp. oocysts and giardia spp. cysts in a tropical water basin cameroon. Water 2010; 2: 44–57.
    
    

77. 77.Spanakos G, Biba A, Mavridou A, Karanis P. Occurrence of Cryptosporidium and Giardia in recycled waters used for irrigation and first description of Cryptosporidium parvum and C. muris in Greece. Parasitology research 2015; 114(5): 1803–10.
    
    

78. 78.Boithias L, Choisy M, Souliyaseng N, et al. Hydrological Regime and Water Shortage as Drivers of the Seasonal Incidence of Diarrheal Diseases in a Tropical Montane Environment. PLOS Neglected Tropical Diseases 2016; 10(12): e0005195.
    
    

79. 79.Boxall AB, Hardy A, Beulke S, et al. Impacts of climate change on indirect human exposure to pathogens and chemicals from agriculture. Environmental health perspectives 2009; 117(4): 508–14.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1289/ehp.0800084&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19440487&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F09%2F01%2F2023.08.31.23294911.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000264704500022&link_type=ISI) 

80. 80.Miller WA, Lewis DJ, Pereira MDG, et al. Farm Factors Associated with Reducing Cryptosporidium Loading in Storm Runoff from Dairies. Journal of environmental quality 2008; 37(5): 1875–82.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.2134/jeq2007.0413&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18689749&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F09%2F01%2F2023.08.31.23294911.atom)