PERSISTENTLY HIGHEST RISK AREAS FOR HANTAVIRUS PULMONARY SYNDROME: POTENTIAL SITES FOR REFUGIA
|
|
- Kelly Doyle
- 5 years ago
- Views:
Transcription
1 Ecological Applications, 17(1), 2007, pp Ó 2007 by the Ecological Society of America PERSISTENTLY HIGHEST RISK AREAS FOR HANTAVIRUS PULMONARY SYNDROME: POTENTIAL SITES FOR REFUGIA GREGORY E. GLASS, 1,4 TIMOTHY SHIELDS, 1 BIN CAI, 1 TERRY L. YATES, 2 AND ROBERT PARMENTER 2,3 1 The W. Harry Feinstone Department of Molecular Microbiology and Immunology, The Johns Hopkins Bloomberg School of Public Health, 615 North Wolfe Street, Baltimore, Maryland USA 2 The Department of Biology and the Museum of Southwestern Biology, The University of New Mexico, Albuquerque, New Mexico USA 3 The Valles Caldera National Trust, State Highway 4, P.O. Box 359, Jemez Springs, New Mexico USA Abstract. Interannual variation in the number of cases of human disease caused by hantaviruses in North America has been hypothesized to reflect environmental changes that influence rodent reservoir populations. This hypothesis postulates that when cases are rare reservoir populations are geographically restricted in patches of suitable habitat. Identifying these sites, which is needed to test the hypothesis, has proven to be a challenge. Satellite imagery of the U.S. Southwest has shown associations among the likelihood of human hantaviral disease and increases in the rodent populations, as well as increased prevalence of Sin Nombre virus (SNV) in rodent populations. In this study we characterize local areas that had environmental signatures that persisted as predicted highest risk sites for human disease through much of the 1990s. These areas represent a small percentage (0.3%) of the region. Exploratory analyses indicate that these areas were not randomly distributed, but were associated with certain landscape characteristics. Characteristics of elevation, slope, aspect, and land cover were associated with persistent high risk. Using multivariate Poisson regression to control for confounding effects, sites with deciduous- or mixed-forest land cover on moderate to steep slopes (.58) above 2130 m elevation were associated with increasing numbers of years at highest risk. These are candidate locations for refugia. Sites associated with cleared ground or shrubland were less often associated with high risk compared to reference conditions. The seasonal patterns of vegetation growth in persistently high-risk areas were compared to matched locations using MODIS (moderate resolution imaging spectroradiometer) NDVI (normalized difference vegetation index) during a time of a severe drought in the region from 2002 to Despite the drought and regardless of land cover, the NDVI in persistently highest risk areas had an early onset, with significantly higher levels of green vegetation that lasted longer than at comparable sites. These observations identify locations that can be monitored for the abundance of P. maniculatus and presence of SNV. If these sites are refugia, we predict they will be occupied by infected deer mice when other monitored sites are unoccupied. Key words: deer mouse; GIS; hantavirus; hantavirus pulmonary syndrome; Peromyscus maniculatus; remote sensing; Sin Nombre virus; trophic cascade hypothesis. INTRODUCTION Hantavirus Pulmonary Syndrome (HPS) is a human disease caused by infection with various hantaviruses. These viruses are zoonotic agents that naturally infect and are maintained in murine rodents (Yates et al. 2002). Extensive surveys, coupled with phylogenetic analyses of rodents and viruses, implicate single species of rodents as hosts of single hantaviruses with occasional spillover to other species, including humans (Childs et al Mills et al. 1997, Yates et al. 2002). HPS was recognized following an outbreak of pulmonary disease among humans in the U.S. Southwest in 1993, and infection with Sin Nombre virus (SNV) was traced to its Manuscript received 25 April 2005; revised 7 March 2006; accepted 3 May Corresponding Editor: N. T. Hobbs. 4 gglass@jhsph.edu 129 reservoir, Peromyscus maniculatus, the deer mouse (Nichol et al. 1993, Childs et al. 1994). The number of cases of HPS varies both geographically and among years (e.g. Hjelle and Glass 2000). The causes of these fluctuations are not clear but are presumed to be due, generally, to changes in contacts between infectious mice and humans (Childs et al. 1995, Hjelle and Glass 2000). Parmenter and colleagues (1993) postulated that one of the main drivers for the interannual variation in HPS cases was changes in environmental conditions that were thought to be due to fluctuations in the patterns of precipitation and temperature. These fluctuations affected the abundance of deer mice and subsequently the human risk of HPS (Parmenter et al. 1993, Yates et al. 2002). They hypothesized that as precipitation increased and temperatures were ameliorated net primary productivity (NPP) increased with attendant increases in ground-
2 130 GREGORY E. GLASS ET AL. Ecological Applications Vol. 17, No. 1 dwelling arthropod populations. Increased NPP acted both directly and indirectly on local deer mouse populations, resulting in larger populations. Increased rodent densities were associated with increased viral transmission among rodent population members and larger numbers of dispersing animals that, in turn, led to increased transmission to humans who came in contact with these animals (Yates et al. 2002). This ultimately became known as the Trophic Cascade Hypothesis (TCH). The TCH is being evaluated in several studies, though much work remains (Mills et al. 1997, Yates et al. 2002). According to the TCH, as environmental conditions become more severe, many local rodent populations decline and become extinct. However, in some locations the environment allows survival of sufficiently large rodent populations that the virus persists. It is from these sites that SNV later reappears and spreads. A key step in evaluating the TCH is to identify the conditions needed to maintain large reservoir populations that carry SNV and where they occur. Large reservoir populations are needed because hantaviruses are horizontally transmitted agents (agents that are spread among members of the population after birth). In the absence of vertical transmission (from mother to offspring, typically in utero), sufficient numbers of susceptible individuals must be recruited and acquire virus from infectious individuals so that the virus does not go extinct locally. Parmenter et al. (1993) proposed that locations where these environmental conditions existed served as refugia for the virus. Studies of rodent reservoir populations using traditional approaches have had difficulty testing the TCH because of challenges in locating refugia. Many of the local rodent populations monitored as part of the longitudinal studies have been subject to repeated extinction and only sporadic SNV infection (e.g., Calisher et al. 2001, Yates et al and references therein), suggesting that these sites do not fit the conditions for refugia. A primary challenge with testing the TCH is that the environmental conditions needed to maintain P. maniculatus and SNV are not defined by the theory, but only by their existence. In earlier work (Glass et al. 2000), the risk of HPS was associated with environmental conditions characterized by satellite imagery nearly 10 months prior to disease outbreak. The best linear classifier, using Landsat Thematic Mapper (available online) 5 (TM) imagery, combined TM bands 1, 5, and 7 and elevation (.2094 m). Subsequent studies of local rodent populations during an El Nin o/southern Oscillation (ENSO) event in examined areas where the environmental signature indicated, for humans, high risk, low risk, or areas that fluctuated between low and high risk during a two-year period. The prevalence of SNV infection 5 hhttp://eros.usgs.gov/products/satellite/tm.htmli started from a common baseline in all cases, and the characteristics of local rodent populations in sites where human risk fluctuated between years were similar to low-risk sites. Sites that were predicted to be high risk for two consecutive years were characterized by portions of the mouse population known to be at highest risk for infection. The prevalence of SNV infection in deer mice from these sites exceeded 30% by the second year of study, even though they had started with a prevalence of infection comparable to low-risk sites (Glass et al. 2002). The apparent need for sites to persist as high risk for at least two years before infection prevalence increased was consistent with the interpretation that virus had been introduced into these local populations as suitable conditions persisted during the study. Here, we characterize physical conditions at sites where high risk of HPS persists for extended periods of time (persistently highest risk areas). The goals were to identify locations that could meet the criteria for refugia and to characterize their physical features with the intent of subsequently monitoring local rodent populations at these sites. We use the onset of a severe, long-term, regional drought to examine the seasonal dynamics of vegetative patterns, as measured by normalized difference vegetation index (NDVI), in persistently highest risk areas, compared to the remainder of the region. The rationale for characterizing the environment based on human risk is based on studies by Childs and colleagues (1995) who showed that the only demonstrable difference between households where HPS occurred and those where it didn t were the number of infected P. maniculatus captured. Thus, sites that are at highest risk for HPS over many years also should be sites with environmental signatures consistently associated with large numbers of infected P. maniculatus. If the TCH is true, then environmental patterns in these areas could correspond to conditions that favor the survival of deer mice populations. METHODS Persistently highest risk (PHR) sites were identified in a region of the U.S. Southwest by extending previously described methods (Glass et al. 2000). Human risk was modeled from the case-control data set with logisticregression analyses, in which linear combinations of Landsat TM 5 bands and elevation.2094 meters were used. Annual models of risk were generated for five years, symmetrically spanning most of the 1990s (1992, 1993, 1995, 1997, and 1998). The annual HPS risk maps were coregistered and overlaid. The highest risk areas within a single year were located by examining the geographic extent for predicted values from the logistic-regression analysis. Thresholds for the predicted value corresponding to ;2.5% of the study area were selected. This choice was arbitrary but was used because it provided good concordance between the receiver operator characteristic (ROC) of NDVI and the ROC of predicted human risk in 1992
3 January 2007 HANTAVIRUS REFUGIA 131 (Glass et al. 2000). We also evaluated alternative thresholds of 5%, 10%, and 20%, which represented more extensive areas, but these did not alter the qualitative results. Exploratory analyses Each year s highest risk areas were dichotomized above and below the 2.5% threshold, and persistently highest risk areas were identified by overlaying the annual risk maps. Sites that were in 2.5% highest risk categories throughout the study period were categorized as PHR. The spatial stability of PHR sites was examined by recording the number of pixels that were in the highest risk areas 0, 1, 2, 3, 4, or 5 years, and these were compared with the observed numbers expected under a binomial distribution with P ¼ If highest risk sites persisted, then an excess of sites in the highest risk categories with four or five years recurrence and a paucity of sites with no or one year in the highest risk categories was expected. Failure to reject the hypothesis would indicate that pixels in highest risk categories varied geographically among years. As part of the exploratory analyses, the spatial relationships of PHR areas with physical and landcover characteristics of the region were examined to determine if potential refugia corresponded to recognizable landscape features. PHR sites were overlaid with digital raster graphics (DRGs) generated from U.S. Geological Survey 1: scale topographic maps. Images were visually examined to identify aspects of the landscape associated with PHR areas. Digital elevation models (DEMs), obtained from the U.S. Geological Survey (USGS) National Elevation Database in 7.5-min data, were used to estimate the elevation, slope, and aspect of sites. The DEM files were coregistered to the PHR areas. Elevation and aspect of PHR areas were compared with the remaining study area by Komolgorov-Smirnov two-sample test and Watson s U 2 test, respectively. The null hypotheses were that PHR areas were random samples of the underlying patterns of elevation and aspect within the study region. Land-cover characteristics were derived from U.S. EPA National Land Cover Data (NLCD) from the Multi-Resolution Land Characteristics (MRLC) Consortium. NLCD level II 1992 classifications were extracted for all pixels in PHR areas, as well as the entire study region. The proportions of the study area and PHR areas were compared in a 2 3 C table (two rows by the number of columns necessary for the comparison). Land-cover categories that individually represented,0.5% of the study region were aggregated into an other class for the exploratory analysis. Multivariate analyses It was anticipated that many environmental variables, such as elevation and land cover, would be correlated among themselves, making it difficult to control for confounding effects. Multivariate Poisson regression was performed to reduce confounding. Outcomes were the numbers of years (0 5) that an individual pixel was scored as highest risk. The predictor variables were derived from the environmental variables used in the exploratory analyses. Indicator variables were generated based on these variables and effects defined relative to reference classes. The large number of pixels in the study area (;98.5 million pixels), made it infeasible to perform a single analysis. Instead, a Monte Carlo analysis was conducted to obtain point and interval estimates of the regression coefficients. Approximately pixels were randomly selected without replacement in the study area and their environmental variables extracted. Poisson regression was performed and the coefficient for each variable was generated while including the effects of all other variables. The process was repeated 200 times, with replacement between samples. The median value of each coefficient derived from the 200 repetitions was chosen as the point estimate for the environmental parameters. The empirical interval estimate excluded the extreme 2.5% of values for each parameter. Parameters with estimates that overlapped 0 were not considered to differ significantly from the reference category but were retained in the model to control for their impact on the remaining predictor variables. Environmental dynamics in PHR areas Seasonal dynamics in vegetation greenness of PHR areas were characterized by the normalized difference vegetation index (NDVI). NDVI values were generated as a data product (MOD13) from the Moderate Resolution Imaging Spectroradiometer (MODIS) sensor for 2002, 2003, and the first quarter Spatial resolution of the MODIS NDVI was 250 m. Temporal resolution was 16 days. NDVI, a commonly used indicator of vegetation growth/productivity (Lillesand et al. 2004), was previously compared with the HPS risk algorithm (Glass et al. 2000). It compared favorably with the risk algorithm in highest and lowest risk areas, though it was not as sensitive in the intermediate range of risk. However, intermediate risk areas were not included in the PHR sites. The PHR image that was overlaid on the NLCD file was used to stratify PHR areas and the remaining portion of the study site by land-cover type. We selected a series of locations from PHR sites and randomly selected lower risk sites stratified into the four most common land-cover categories: evergreen forest, shrubland, grasslands/herbaceous, and pasture/hay accounted for more than 95% of the region. Site selection was restricted by requiring that the location be at least 10 pixels from boundaries between PHR areas and the remaining image in each of the land-cover classifications. This was done to reduce the potential misclassification of the selected locations. The changes in NDVI for PHR and remaining sites were compared, within land-cover categories, using time series analyses to determine if the onset, intensity, or
4 132 GREGORY E. GLASS ET AL. Ecological Applications Vol. 17, No. 1 sites. We also expected that SNV should be introduced into sites closest to PHR areas than sites farther away. PLATE 1. Persistently high risk area associated with Hantavirus Pulmonary Syndrome. Photo credit: Douglas E. Norris. duration of vegetation greenness differed. Seasonal Autoregressive Integrated Moving Average (SARIMA) models were fit to the data. Differences between PHR and lower risk sites within strata were identified by the need to postulate alternative models to characterize the comparable time series, such as different model structures or significant differences in the coefficients and lag structures for individual models. Yearly changes in abundance of local deer mouse (P. maniculatus) populations near PHR sites were examined using data from the 38 locations which had been previously sampled for small rodent populations and tested for SNV during 1998 and Approximately half the sites were in high-risk areas during two years of the study, while the remaining locations were either in low-risk areas or sites that varied in risk during the two year study. Counts of deer mice trapped using standard protocols and the presence of mice with antibodies to SNV, as tested by ELISA, were used as measures of population abundance and the presence of virus. Methods of sampling and testing are described in detail elsewhere (Glass et al. 2002). We expected that if PHR sites were refugia, that populations nearer PHR sites were more likely to retain resident populations of deer mice and were more likely to be recolonized than sites in lower quality habitat that were further away from PHR RESULTS Persistently high risk (PHR) sites represented a small fraction of the study area, accounting for 0.3% (311 km 2 out of km 2 ) of the region. Despite their limited geographic extent, PHR areas had a coherent structure related to physical features, elevation, land cover, and seasonal vegetative growth patterns (see Plate 1). Sites that had previously been sampled for deer mice near PHR locations were more likely, over the two year period, to acquire and retain deer mouse populations that were infected with SNV. PHR sites were spatially stable across years. Sites that were at highest risk during three, four, or five years of the study were approximately 138, 1300, and times more likely, respectively, than expected than if the highest risk pixels were randomly distributed (v 2 ¼ , df ¼ 5, P, 0.001) among years (Fig. 1). These sites frequently corresponded to well-defined physical features on the digital raster graphics (DRGs). Often, they occurred along arroyos or canyons associated with mountains and mesas (Fig. 2). However, even within the same canyon, only portions were PHR areas. These sites occurred adjacent to streams, rivers, or other water bodies. Near these features, they occurred at slightly higher elevations along the drainages rather than immediately adjacent to the watercourses themselves. PHR sites were predominantly restricted to a narrow range of elevations compared to the study region as a whole (Komolgorov-Smirnov two-sample test ¼ 0.76; P, ). PHR sites did not occur below 1600 m elevation, even though more than one third of the study region was lower (Fig. 3). Approximately 80% of the PHR sites were within a narrow 400 m elevation band between 2600 m and 3000 m, although only 6% of the study region occurred within this range. PHR sites also were not found above 3200 m, even though parts of the study region reached 3600 m. In exploratory analyses, PHR sites were more frequent in areas with east FIG. 1. Temporal persistence of highest risk sites. The histogram shows observed number of pixels rated as the highest risk category/expected number of highest risk pixels if they were under random distribution. Sites classified in the highest risk category tended to persist across years.
5 January 2007 HANTAVIRUS REFUGIA 133 FIG. 2. PHR sites (red circles) were frequently associated with arroyos and canyons on mesas or mountains (panel A; DRG NJ 12-11, Series V502 quadrangle) or along streams and rivers (panel B; DRG NJ 12-12, Series V502). Locations are indicated by Digital Raster Graphics (DRG) identifiers assigned by the U.S. Geological Survey. northeast east southeast aspects than the study region as a whole (16.4% vs. 12.2%) (Watson s U 2 ¼ 5.78, P, 0.001) and they were less common on flat areas than expected (27.8% vs. 40.9%). This, in part, reflected their absence in lower elevation areas, where much of the flat ground occurred. PHR areas were disproportionately associated with evergreen forest classification in the NLCD and negatively associated with shrubland; these two classes accounted for ;80% of the study area (v 2 ¼ 82.45, df ¼ 1, P ). Most of the remaining PHR areas were grassland or pasture, though in proportions consistent with their occurrence in the region (Table 1). Although more than two-thirds of PHR areas were associated with evergreen forest, this still represented only 0.65% of all the land with this classification (see Plate 1). In contrast, some land-cover classes that occupied vanishingly small portions of the region, such as herbaceous wetlands
6 134 GREGORY E. GLASS ET AL. Ecological Applications Vol. 17, No. 1 FIG. 3. Comparison of elevations at PHR sites (persistently highest risk) and those within the study region (background elevation). PHR sites occurred within a narrow range of elevations. Approximately 80% of the PHR areas occurred between 2600 and 3000 m, while only 6% of the study area was at these elevations. (accounting for 0.004% of the region), were frequently (18.3%) classified as PHR. In multivariate analyses, the number of years a site was in the highest risk level was related to a combination of land cover, elevation, and slope (Table 2). Referent sites were agricultural or pasture land, below 1801 m elevation on flat ground. Bare ground and shrubland were associated with significantly fewer numbers of years at highest risk, compared to referent sites, while mixed and deciduous forest classified areas were associated with significantly more years as highest risk. Evergreen forest was no more likely to contribute additional years at highest risk than the referent landcover class. Increasing elevation was associated with increasing years at highest risk, with the effect markedly increasing above 2130 m (Table 2). Similarly, compared to flat ground, moderate and steep slopes (.48) were associated with more years at highest risk. After controlling for land cover, elevation, and slope, aspect did not contribute significantly to the numbers of years at highest risk (Table 2). Eight to twelve sites, in each of the four major landcover classes (evergreen forest, shrubland, grassland, or pasture) had NDVI values extracted from MODIS NDVI for 848 days from January 2002 through April 2004 (Fig. 4A D). Regardless of land cover, PHR sites showed a marked seasonal change in NDVI values with a rapid increase beginning about day 64 of each year (Fig. 4, times 4, 29, and 53). NDVI values in PHR sites peaked over a broad time period, between day 160 and 320 of each year, and then declined. In contrast, lower risk sites showed a slower greening through the growing season and only occasionally reached NDVI values approaching those in PHR areas (Fig. 4; Appendix A: Table A1). Autocorrelation and partial autocorrelation plots were used to identify appropriate models for NDVI each category of sites (Appendix B: Table B1). The time series were differenced to achieve stationarity. The most common class of models was seasonal-mixed-autoregressive moving-average models in which the autoregressive and moving-average components were first order and the seasonal-moving-average component also was first order, with a seasonal component between 352 and 384 days. The best models fitting the seasonal changes in NDVI at PHR sites in evergeen forests, shrublands, and grasslands differed substantially from corresponding sites that were not highest risk (Appendix B: Table B1). Human-modified pastureland sites, whether PHR or not, (Fig. 4D; Appendix B: Table B1), were most similar to one another with similar AR(1) structure, though they differed in their seasonal MA coefficients. Sites sampled for small mammals and tested for SNV in 1998 and 1999 were examined relative to PHR areas. None of the previously trapped sites were within any PHR areas, although 13 of the trapped locations were within a kilometer, and three were within 200 m (Fig. 5). On average, trapped high-risk sites were nearer (2.0 km km; mean 6 SD) than trapped low-risk sites (13.5 km km) to PHR areas. No deer mice were trapped at seven sites in Four sites were low-risk and three were high-risk areas (Fig. 5). During resampling in 1999, only a single deer mouse was captured at one of these four low-risk sites and it was uninfected with SNV, indicating colonization by mice and introduction of SNV was rare at these locations. In addition, seven of the 11 low-risk sites that had deer mice in 1998 lost their mice by 1999 (Fig. 5), indicative of high local extinction rates in these habitats. By 1999, 10 out of 18 low-risk sites were unoccupied. Infected deer mice were found at only one of the low-risk sites. Each of the three unoccupied high-risk sites in 1998 yielded deer mice in 1999, and at least one mouse at each Comparison of persistently highest risk sites and the remaining portion of the study area by major land-cover classes. TABLE 1. Land cover Persistently highest risk (%) Not highest risk (%) Evergreen forest Shrubland Grassland Pasture Other Notes: Values represent the percentage of the study area that fell into a risk category, stratified by land cover. All classes with.0.5% coverage were categorized separately. Land-cover data were based on the National Land Cover Data (NLCD) 1992 database. Other: Open water, low-density residential, commercial, bare rock, quarries, deciduous forest, mixed forest, orchards, row crops, small grains, parks, woody wetlands, herbaceous wetlands.
7 January 2007 HANTAVIRUS REFUGIA 135 Effect of selected environmental variables on the number of years pixels were identified as persistently highest risk. TABLE 2. Variable Referent class Regression coefficient 95% interval estimate Residential/industrial agriculture/pasture , 0.93 Bare ground/quarries agriculture/pasture , 0.09 Deciduous/mixed forest agriculture/pasture , 0.68 Evergreen agriculture/pasture , 0.12 Shrubland agriculture/pasture , 0.30 Wetlands agriculture/pasture , 1.09 Elevation ( m) elevation 1800 m , 0.10 Elevation ( m) elevation 1800 m , 0.44 Elevation ( m) elevation 1800 m , 3.65 Elevation (.2800 m) elevation 1800 m , 4.43 Slope (18 48) flat ground , 0.37 Slope (58 258) flat ground , 0.68 Slope (.258) flat ground , 0.84 Aspect ( ) flat ground , 0.24 Aspect (18 448) flat ground , 0.34 Aspect ( ) flat ground , 0.38 Aspect ( ) flat ground , 0.38 Aspect ( ) flat ground , 0.25 Aspect ( ) flat ground , 0.25 Aspect ( ) flat ground ; , 0.25 Aspect ( ) flat ground , 0.24 Notes: Multivariate Poisson regression was used to control for confounding effects by estimating coefficients from a randomly selected pixels in the region. The models were repeated 200 times, and empirical point and interval estimates were derived as the median value for each variable and the 95% interval. Variables that differ significantly from zero (as derived by Monte Carlo estimates; see Methods: Multivariate analysis) are indicated by daggers ( ). All changes are in comparison to a selected category (referent class). of these sites had antibodies to SNV (Fig. 5). All of the high-risk sites had deer mice by 1999 and 18 out of 20 had at least one deer mouse infected with SNV, suggesting successful population persistence and that viral introduction was common in suitable habitats near PHR sites. Sites with infected deer mice were, with one exception, within 8 km of PHR areas. Sites at farther distances either no longer had deer mice by the second year, or those that were present were not infected. DISCUSSION Sites with persistent environmental signatures similar to where people were at highest risk for HPS occurred in very limited locations within a geographic region that is viewed as a major epicenter for disease. If PHR sites represent refugia, their limited extent could explain previous workers difficulty locating them. PHR sites had well-defined physical characteristics and were focused in a narrow elevational range with moderate to steep slopes. Land-cover categories were not good descriptors of increasing risk, though this may reflect the habitat generalist nature of the rodent reservoir or the coarse categorization of land classes in the data set. Deer mice occur in diverse habitats throughout much of North America, and in the U.S. Southwest, they range from desert to alpine habitats (Mills et al. 1997). Despite this wide range of tolerances, Mills and colleagues found that the abundance of P. maniculatus varied with altitude and vegetation community, as did the prevalence of Sin Nombre Virus (SNV) in the local rodent populations. They found deer mice less commonly at the lowest and highest elevations and prevalence of SNV decreased at the elevation extremes. Nearly half (20/41 sites) of their locations where deer mice occurred yielded no mice with antibodies to SNV. Occasionally this was due to small sample sizes, but many of their sites had either a disproportionate overabundance or paucity of mice with SNV. The focal nature of infection was unrelated to population density and the bases for the patterns remained unexplained. However, their results mirror the altitudinal and land cover distributions of PHR sites, with a disproportionate fraction of sites above 2100 m associated with a variety of land-cover categories. PHR and adjacent areas also may account for the focal distribution of SNV infection in mouse populations, as areas adjacent to PHR sites were more likely to acquire and retain infected deer mice than sites at greater distances (Fig. 5). Similarly, Calisher and colleagues (2001) have hypothesized that local edaphic and vegetation conditions drive reservoir population dynamics to favor either an excess or paucity of infected population members by affecting rates of reproduction and altering age structure so that in some sites long-lived individuals serve as the sources for SNV infection and maintain transmission. The primary advantage of the current analysis is that it provides the detailed characterization of these sites, as well as their geographic locations. A key factor distinguishing PHR sites from other sites was the seasonal pattern of vegetative growth (Fig. 4).
8 136 GREGORY E. GLASS ET AL. Ecological Applications Vol. 17, No. 1 FIG. 4. Seasonal changes in normalized difference vegetation index (NDVI) derived from the Moderate Resolution Imaging Spectroradiometer (MODIS) between 2002 and the first quarter of 2004 at selected PHR sites (dotted lines) and sites that were not PHR (solid lines). Colors of lines were assigned to individual locations arbitrarily. Sites are stratified by land cover: (A) evergreen forest, (B) shrubland, (C) grassland, and (D) pasture. These four land cover categories accounted for.98% of the study area. PHR sites tended to show an earlier, more rapid and more extensive vegetative growth than did comparable lower risk sites. Even during an extensive drought, PHR sites had an early season onset of vegetation growth (early March), with an extensive growing season (until late October to early November), before a rapid decline with the onset of winter conditions, typically in mid- to late November. The ability of vegetation in PHR sites to show such vigorous and long-lasting growing seasons often corresponded to locations along arroyos, drainage systems where moisture is likely to accumulate and persist. The association with increased moisture also is consistent with the increase in years at risk with deciduous and mixed deciduous land cover (Table 2) which, in this region, is predominantly associated with oaks (Quercus spp). Thus, although the region is generally viewed as high risk for HPS, the real area of consistent, long-term, high risk is quite limited. The specific characteristics of these sites may explain why during most years there are few human cases, while during some years larger numbers of HPS are observed. The increased number of years at risk associated with sites.2130 m (Table 2) occurs near the maximum elevation occupied by most people in the region. The altitudinal distribution of the human control sites observed by Glass et al. (2000) was 2370 m (1962 6
9 January 2007 HANTAVIRUS REFUGIA 137 FIG. 4. Continued. 198 m; mean 6 SD). PHR sites, in particular, were predominantly within the m band (Fig. 3), above the elevation of human residence. Because few people reside at elevations with PHR sites, relatively few are likely to become infected during most years. If PHR areas are refugia, then during most years mice would need to disperse from these locations to lower elevations before reaching human habitations. Distance from PHR sites might be a surrogate of human risk and the timing of pathogen acquisition. The current analysis does not establish that refugia for SNV exist or that the PHR areas have identified them. However, it shows that areas with environmental signatures similar to those that distinguish high-risk locations for HPS from lower risk sites may persist for many years. These sites may represent refugia because the physical and ecological conditions are consistent with environmental features that favor survival and reproduction in deer mouse populations (Yates et al. 2002). These sites also have characteristics similar to those favoring SNV persistence (Childs et al. 1995, Mills et al. 1997). The observation that many of the highest risk sites remained as such for four or five of the years studied (Fig. 1) suggests that they could serve as sites where SNV reaches high levels in local reservoir populations. A single year as a high-risk area is not sufficient for SNV to reach high levels in reservoir populations, but within
10 138 GREGORY E. GLASS ET AL. Ecological Applications Vol. 17, No. 1 FIG. 5. Number of deer mice captured in 1998 (dark bars) and 1999 (light bars) with increasing distance (km) from PHR sites. Classified high-risk sites (indicated by daggers above bars) tended to be closer to PHR sites than low-risk sites. Low-risk sites were more likely to lose deer mice between 1998 and two years as a high-risk site, the average prevalence of infection in local populations can reach more than 30% (Glass et al. 2002), especially near PHR areas. This is consistent with the original epidemiologic investigations of HPS cases (Childs et al. 1995) where, in case households, the prevalence of SNV infection among mice was more than 30%. By implication, the Landsat image classifier identifies environmental conditions that are associated with large deer mouse populations that have high prevalence of SNV. If PHR sites are refugia, then we expect that sampling in locations with these environmental signatures will yield deer mice at times when sampling in other locations finds few if any of this species. In addition, local populations of deer mice sampled in these sites should be more likely infected with SNV than samples of mice collected in sites that are not persistently highest risk locations. The loss of deer mice from low-risk sites in 1999 (Fig. 5) and the colonization of high-risk sites near PHR areas, with the introduction of SNV, is consistent with this hypothesis. If the trophic cascade hypothesis (TCH) is true, then sampling with increasing distance from PHR should be able to document the movement of individual mice from PHR to lower elevations as environmental conditions become more favorable and net primary productivity (NPP) increases. According to the TCH, the movement of mice from refugia should conform to metapopulation theory (Hanski 1999), and we expect that locations away from refugia should be colonized as a decreasing function of distance from the refugia. Local extinction among these colonized sites should be inversely related to local habitat quality and patch size and should occur as environmental conditions become more severe. Further detailed studies are needed to characterize the ecological dynamics of PHR areas, compared to similar areas that are not PHR, to better understand the dynamics that influence deer mouse populations. Additional studies also are needed to evaluate the quality of the environmental data layers used in these analyses. A primary advantage of using these environmental data is the well-characterized metadata. However, groundtruthing is needed to clarify some issues. For example, the PHR areas in NLCD-classified (National Land Cover Data) pasture and grassland, though representing small geographic regions, were unexpected because these are not usually considered suitable habitats for deer mice. This could indicate misclassification in the landcover database. However, Mills and colleagues (1997) reported deer mice from some of their sites in these habitats and Parmenter (R. Parmenter, unpublished data) has observed relatively large local populations during the autumn in these habitats in northern New Mexico. Future ground-based studies could target these apparently anomalous locations to determine the reasons for the observations and evaluate if these locations truly are suitable for reservoir species. Regardless, combining remotely sensed environmental data with ground studies provides a rational strategy for targeted field studies of local populations of deer mice that are responsible for transmission of SNV. By identifying the physical structure and spatial context of PHR locations, we can use the geographic extent and contiguity of local populations to create testable
11 January 2007 HANTAVIRUS REFUGIA 139 hypotheses concerning the persistence and movement of individuals responsible for viral transmission. ACKNOWLEDGMENTS We thank the many individuals who assisted with the field aspects of the study, especially those individuals who conducted the sampling efforts in 1998 and Financial support for this study was provided by funding from the NSF and NIH program in Ecology of Infectious Diseases DEB , NASA NCC5-305, and NOAA NA96GP0419. Comments from M. L. Farnsworth and an anonymous reviewer were especially helpful and improved the quality of the final manuscript. LITERATURE CITED Calisher, C. H., J. N. Mills, W. P. Sweeney, J. R. Choate, D. E. Sharp, K. M. Canestorp, and B. J. Beaty Do unusual site-specific population dynamics of rodent reservoirs provide clues to the natural history of hantaviruses? Journal of Wildlife Diseases 37: Childs, J. E., et al Serological and genetic identification of Peromyscus maniculatus as the primary rodent reservoir for a new hantavirus in the southwestern United States. Journal of Infectious Diseases 169: Childs, J. E., et al A household-based, case-control study of environmental factors associated with hantavirus pulmonary syndrome in the southwestern United States. American Journal of Tropical Medicine and Hygiene 52: Glass, G. E., et al Anticipating risk areas for hantavirus pulmonary syndrome with remotely sensed data: re-examination of the 1993 outbreak. Emerging Infectious Diseases 6: Glass, G. E., et al Satellite imagery characterizes local animal reservoir populations of Sin Nombre virus in Southwestern United States. Proceedings of the National Academy of Sciences (USA) 99: Hanski, I Metapopulation ecology. Oxford University Press, Oxford, UK. Hjelle, B., and G. E. Glass Outbreak of hantavirus infection in the Four Corners Region of the US in the wake of the El Niño Southern Oscillation. Journal of Infectious Diseases 181: Lillesand, T. M., R. W. Kieffer, and J. W. Chipman Remote sensing and image interpretation. Fifth edition. John Wiley and Sons, Hoboken, New Jersey, USA. Mills, J. N., et al Patterns of association with host and habitat: antibody reactive with Sin Nombre virus in small mammals in the major biotic communities of the southwestern United States. American Journal of Tropical Medicine and Hygiene 56: Nichol, S. T., C. F. Spiropoulou, S. Morzunov, P. E. Rollin, T. G. Ksiazek, H. Feldmann, A. Sanchez, S. Zaki, J. Childs, and C. J. Peters Genetic identification of a novel hantavirus associated with an outbreak of acute respiratory illness in the southwestern United States. Science 262: Parmenter, R. R., J. W. Brunt, D. I. Moore, and M. S. Ernest The Hantavirus epidemic in the Southwest: rodent population dynamics and the implications for transmission of hantavirus-associated adult respiratory distress syndrome (HARDS) in the Four Corners Region. Pages 1 45 in Sevilleta Publication No. 41. Sevilleta Long Term Ecological Research (LTER) Project, University of New Mexico, Albuquerque, New Mexico, USA Yates, T. L., et al The ecology and evolutionary history of an emergent disease: Hantavirus Pulmonary Syndrome. BioScience 52: APPENDIX A Descriptive statistics for NDVI in persistently highest risk (refugia) and other (not refugia) areas classified as evergreen forest, shrubland, grassland, and pasture (Ecological Archives A A1). APPENDIX B Coefficients from SARIMA models for average NDVI values in four major land-cover classes (evergreen forest, grassland, shrubland, pasture) in highest risk areas (refugia) and lower risk areas (not refugia) (Ecological Archives A A2).
APPENDIX. Normalized Difference Vegetation Index (NDVI) from MODIS data
APPENDIX Land-use/land-cover composition of Apulia region Overall, more than 82% of Apulia contains agro-ecosystems (Figure ). The northern and somewhat the central part of the region include arable lands
More informationSpatiotemporal patterns in the Hantavirus infection
Spatiotemporal patterns in the Hantavirus infection G. Abramson 1,2, * and V. M. Kenkre 1, 1 Center for Advanced Studies and Department of Physics and Astronomy, University of New Mexico, Albuquerque,
More informationTraveling Waves of Infection in the Hantavirus Epidemics
Bulletin of Mathematical Biology (2003) 65, 519 534 doi:10.1016/s0092-8240(03)00013-2 Traveling Waves of Infection in the Hantavirus Epidemics G. ABRAMSON AND V. M. KENKRE Center for Advanced Studies and
More informationData Fusion and Multi-Resolution Data
Data Fusion and Multi-Resolution Data Nature.com www.museevirtuel-virtualmuseum.ca www.srs.fs.usda.gov Meredith Gartner 3/7/14 Data fusion and multi-resolution data Dark and Bram MAUP and raster data Hilker
More informationCHAPTER 1 THE UNITED STATES 2001 NATIONAL LAND COVER DATABASE
CHAPTER 1 THE UNITED STATES 2001 NATIONAL LAND COVER DATABASE Collin Homer*, Jon Dewitz, Joyce Fry, and Nazmul Hossain *U.S. Geological Survey (USGS) Center for Earth Resources Observation and Science
More informationSensitivity of AERSURFACE Results to Study Area and Location. Paper 2009-A-127-AWMA
Sensitivity of AERSURFACE Results to Study Area and Location Paper 2009-A-127-AWMA Prepared by: Anthony J. Schroeder, CCM Senior Consultant George J. Schewe, CCM, QEP Principal Consultant Trinity Consultants
More information2011 National Seasonal Assessment Workshop for the Eastern, Southern, & Southwest Geographic Areas
2011 National Seasonal Assessment Workshop for the Eastern, Southern, & Southwest Geographic Areas On January 11-13, 2011, wildland fire, weather, and climate met virtually for the ninth annual National
More informationREMOTE SENSING ACTIVITIES. Caiti Steele
REMOTE SENSING ACTIVITIES Caiti Steele REMOTE SENSING ACTIVITIES Remote sensing of biomass: Field Validation of Biomass Retrieved from Landsat for Rangeland Assessment and Monitoring (Browning et al.,
More informationPrincipals and Elements of Image Interpretation
Principals and Elements of Image Interpretation 1 Fundamentals of Photographic Interpretation Observation and inference depend on interpreter s training, experience, bias, natural visual and analytical
More informationSpatial Process VS. Non-spatial Process. Landscape Process
Spatial Process VS. Non-spatial Process A process is non-spatial if it is NOT a function of spatial pattern = A process is spatial if it is a function of spatial pattern Landscape Process If there is no
More informationClassification of Erosion Susceptibility
GEO327G: GIS & GPS Applications in Earth Sciences Classification of Erosion Susceptibility Denali National Park, Alaska Zehao Xue 12 3 2015 2 TABLE OF CONTENTS 1 Abstract... 3 2 Introduction... 3 2.1 Universal
More informationThe Relationship between Vegetation Changes and Cut-offs in the Lower Yellow River Based on Satellite and Ground Data
Journal of Natural Disaster Science, Volume 27, Number 1, 2005, pp1-7 The Relationship between Vegetation Changes and Cut-offs in the Lower Yellow River Based on Satellite and Ground Data Xiufeng WANG
More informationDROUGHT ASSESSMENT USING SATELLITE DERIVED METEOROLOGICAL PARAMETERS AND NDVI IN POTOHAR REGION
DROUGHT ASSESSMENT USING SATELLITE DERIVED METEOROLOGICAL PARAMETERS AND NDVI IN POTOHAR REGION Researcher: Saad-ul-Haque Supervisor: Dr. Badar Ghauri Department of RS & GISc Institute of Space Technology
More informationQuality Assessment of Shuttle Radar Topography Mission Digital Elevation Data. Thanks to. SRTM Data Collection SRTM. SRTM Galapagos.
Quality Assessment of Shuttle Radar Topography Mission Digital Elevation Data Third International Conference on Geographic Information Science College Park, Maryland, October 20-23 Ashton Shortridge Dept.
More informationEpiMAN-TB, a decision support system using spatial information for the management of tuberculosis in cattle and deer in New Zealand
EpiMAN-TB, a decision support system using spatial information for the management of tuberculosis in cattle and deer in New Zealand J.S. McKenzie 1, R.S. Morris 1, C.J. Tutty 2, D.U. Pfeiffer 1 Dept of
More informationThe Application of Remote Sensing and GIS Tools in the Study of Lyme Disease Risk Prediction Kathryn Berger
The Application of Remote Sensing and GIS Tools in the Study of Lyme Disease Risk Prediction Kathryn Berger Lyme disease has become one of the most prevalent vector borne diseases in the United States
More informationNR402 GIS Applications in Natural Resources. Lesson 9: Scale and Accuracy
NR402 GIS Applications in Natural Resources Lesson 9: Scale and Accuracy 1 Map scale Map scale specifies the amount of reduction between the real world and the map The map scale specifies how much the
More informationDROUGHT RISK EVALUATION USING REMOTE SENSING AND GIS : A CASE STUDY IN LOP BURI PROVINCE
DROUGHT RISK EVALUATION USING REMOTE SENSING AND GIS : A CASE STUDY IN LOP BURI PROVINCE K. Prathumchai, Kiyoshi Honda, Kaew Nualchawee Asian Centre for Research on Remote Sensing STAR Program, Asian Institute
More informationPROGRESS REPORT for COOPERATIVE BOBCAT RESEARCH PROJECT. Period Covered: 1 January 31 May Prepared by
PROGRESS REPORT for COOPERATIVE BOBCAT RESEARCH PROJECT Period Covered: 1 January 31 May 2011 Prepared by John A. Litvaitis, Derek Broman, and Marian K. Litvaitis Department of Natural Resources University
More informationPutative Canada Lynx (Lynx canadensis) Movements across I-70 in Colorado
Putative Canada Lynx (Lynx canadensis) Movements across I-70 in Colorado INTRODUCTION March 8, 2012 Jake Ivan, Mammals Researcher Colorado Parks and Wildlife 317 W. Prospect Fort Collins, CO 80526 970-472-4310
More informationRemote Sensing and Geospatial Application for Wetlands Mapping, Assessment, and Mitigation
Remote Sensing and Geospatial Application for Wetlands Mapping, Assessment, and Mitigation Hydrology Soils MSU Seminar Series Remote Sensing and Geospatial Applications September 4, 2002 Vegetation NEPA
More informationThe New Normal or Was It?
The New Normal or Was It? by Chuck Coffey The recent drought has caused many to reflect upon the past and wonder what is in store for the future. Just a couple of years ago, few agricultural producers
More informationAn Internet-based Agricultural Land Use Trends Visualization System (AgLuT)
An Internet-based Agricultural Land Use Trends Visualization System (AgLuT) Second half yearly report 01-01-2001-06-30-2001 Prepared for Missouri Department of Natural Resources Missouri Department of
More informationOBJECT-BASED CLASSIFICATION USING HIGH RESOLUTION SATELLITE DATA AS A TOOL FOR MANAGING TRADITIONAL JAPANESE RURAL LANDSCAPES
OBJECT-BASED CLASSIFICATION USING HIGH RESOLUTION SATELLITE DATA AS A TOOL FOR MANAGING TRADITIONAL JAPANESE RURAL LANDSCAPES K. Takahashi a, *, N. Kamagata a, K. Hara b a Graduate School of Informatics,
More informationUSING GIS CARTOGRAPHIC MODELING TO ANALYSIS SPATIAL DISTRIBUTION OF LANDSLIDE SENSITIVE AREAS IN YANGMINGSHAN NATIONAL PARK, TAIWAN
CO-145 USING GIS CARTOGRAPHIC MODELING TO ANALYSIS SPATIAL DISTRIBUTION OF LANDSLIDE SENSITIVE AREAS IN YANGMINGSHAN NATIONAL PARK, TAIWAN DING Y.C. Chinese Culture University., TAIPEI, TAIWAN, PROVINCE
More information7.1 INTRODUCTION 7.2 OBJECTIVE
7 LAND USE AND LAND COVER 7.1 INTRODUCTION The knowledge of land use and land cover is important for many planning and management activities as it is considered as an essential element for modeling and
More informationRating of soil heterogeneity using by satellite images
Rating of soil heterogeneity using by satellite images JAROSLAV NOVAK, VOJTECH LUKAS, JAN KREN Department of Agrosystems and Bioclimatology Mendel University in Brno Zemedelska 1, 613 00 Brno CZECH REPUBLIC
More informationIMPROVING REMOTE SENSING-DERIVED LAND USE/LAND COVER CLASSIFICATION WITH THE AID OF SPATIAL INFORMATION
IMPROVING REMOTE SENSING-DERIVED LAND USE/LAND COVER CLASSIFICATION WITH THE AID OF SPATIAL INFORMATION Yingchun Zhou1, Sunil Narumalani1, Dennis E. Jelinski2 Department of Geography, University of Nebraska,
More informationGRAPEVINE LAKE MODELING & WATERSHED CHARACTERISTICS
GRAPEVINE LAKE MODELING & WATERSHED CHARACTERISTICS Photo Credit: Lake Grapevine Boat Ramps Nash Mock GIS in Water Resources Fall 2016 Table of Contents Figures and Tables... 2 Introduction... 3 Objectives...
More informationLecture Topics. 1. Vegetation Indices 2. Global NDVI data sets 3. Analysis of temporal NDVI trends
Lecture Topics 1. Vegetation Indices 2. Global NDVI data sets 3. Analysis of temporal NDVI trends Why use NDVI? Normalize external effects of sun angle, viewing angle, and atmospheric effects Normalize
More informationAppendix J Vegetation Change Analysis Methodology
Appendix J Vegetation Change Analysis Methodology Regional Groundwater Storage and Recovery Project Draft EIR Appendix-J April 2013 APPENDIX J- LAKE MERCED VEGETATION CHANGE ANALYSIS METHODOLOGY Building
More informationChitra Sood, R.M. Bhagat and Vaibhav Kalia Centre for Geo-informatics Research and Training, CSK HPKV, Palampur , HP, India
APPLICATION OF SPACE TECHNOLOGY AND GIS FOR INVENTORYING, MONITORING & CONSERVATION OF MOUNTAIN BIODIVERSITY WITH SPECIAL REFERENCE TO MEDICINAL PLANTS Chitra Sood, R.M. Bhagat and Vaibhav Kalia Centre
More informationVCS MODULE VMD0018 METHODS TO DETERMINE STRATIFICATION
VMD0018: Version 1.0 VCS MODULE VMD0018 METHODS TO DETERMINE STRATIFICATION Version 1.0 16 November 2012 Document Prepared by: The Earth Partners LLC. Table of Contents 1 SOURCES... 2 2 SUMMARY DESCRIPTION
More informationWeather and Climate Summary and Forecast February 2018 Report
Weather and Climate Summary and Forecast February 2018 Report Gregory V. Jones Linfield College February 5, 2018 Summary: For the majority of the month of January the persistent ridge of high pressure
More informationPRMS WHITE PAPER 2014 NORTH ATLANTIC HURRICANE SEASON OUTLOOK. June RMS Event Response
PRMS WHITE PAPER 2014 NORTH ATLANTIC HURRICANE SEASON OUTLOOK June 2014 - RMS Event Response 2014 SEASON OUTLOOK The 2013 North Atlantic hurricane season saw the fewest hurricanes in the Atlantic Basin
More informationCase Study: Ecological Integrity of Grasslands in the Apache Highlands Ecoregion
Standard 9: Screen all target/biodiversity element occurrences for viability or ecological integrity. Case Study: Ecological Integrity of Grasslands in the Apache Highlands Ecoregion Summarized from: Marshall,
More informationIntroduction to Geographic Information Systems (GIS): Environmental Science Focus
Introduction to Geographic Information Systems (GIS): Environmental Science Focus September 9, 2013 We will begin at 9:10 AM. Login info: Username:!cnrguest Password: gocal_bears Instructor: Domain: CAMPUS
More informationWill it be a Good Ski Season? Correlation between El Niño and U.S. Weather
Will it be a Good Ski Season? Correlation between El Niño and U.S. Weather Lillian Petersen 1 Abstract This study investigates the correlation between El Niño and the weather across the U.S. El Niño is
More informationPreliminary Runoff Outlook February 2018
Preliminary Runoff Outlook February 2018 Prepared by: Flow Forecasting & Operations Planning Water Security Agency General Overview The Water Security Agency (WSA) is preparing for 2018 spring runoff including
More informationThe Documentation of Extreme Hydrometeorlogical Events: Two Case Studies in Utah, Water Year 2005
The Documentation of Extreme Hydrometeorlogical Events: Two Case Studies in Utah, Water Year 2005 Tim Bardsley1*, Mark Losleben2, Randy Julander1 1. USDA, NRCS, Snow Survey Program, Salt Lake City, Utah.
More information1 Introduction: 2 Data Processing:
Darren Janzen University of Northern British Columbia Student Number 230001222 Major: Forestry Minor: GIS/Remote Sensing Produced for: Geography 413 (Advanced GIS) Fall Semester Creation Date: November
More informationResolving habitat classification and structure using aerial photography. Michael Wilson Center for Conservation Biology College of William and Mary
Resolving habitat classification and structure using aerial photography Michael Wilson Center for Conservation Biology College of William and Mary Aerial Photo-interpretation Digitizing features of aerial
More informationChapter 6. Fundamentals of GIS-Based Data Analysis for Decision Support. Table 6.1. Spatial Data Transformations by Geospatial Data Types
Chapter 6 Fundamentals of GIS-Based Data Analysis for Decision Support FROM: Points Lines Polygons Fields Table 6.1. Spatial Data Transformations by Geospatial Data Types TO: Points Lines Polygons Fields
More informationNational Wildland Significant Fire Potential Outlook
National Wildland Significant Fire Potential Outlook National Interagency Fire Center Predictive Services Issued: September, 2007 Wildland Fire Outlook September through December 2007 Significant fire
More informationRelationship between weather factors and survival of mule deer fawns in the Peace Region of British Columbia
P E A C E R E G I O N T E C H N I C A L R E P O R T Relationship between weather factors and survival of mule deer fawns in the Peace Region of British Columbia by: Nick Baccante and Robert B. Woods Fish
More informationTrait-based Australian mammal distribution patterns and extinction risks
Trait-based Australian mammal distribution patterns and extinction risks Jessica Berryman Salt Lake Community College Society of Conservation GIS Conference Pacific Grove, CA July 27, 2015 Extinction in
More informationUNIT 5: ECOLOGY Chapter 15: The Biosphere
CORNELL NOTES Directions: You must create a minimum of 5 questions in this column per page (average). Use these to study your notes and prepare for tests and quizzes. Notes will be stamped after each assigned
More informationRESEARCH METHODOLOGY
III. RESEARCH METHODOLOGY 3.1. Time and Research Area The field work was taken place in primary forest around Toro village in Lore Lindu National Park, Indonesia. The study area located in 120 o 2 53 120
More informationPriority areas for grizzly bear conservation in western North America: an analysis of habitat and population viability INTRODUCTION METHODS
Priority areas for grizzly bear conservation in western North America: an analysis of habitat and population viability. Carroll, C. 2005. Klamath Center for Conservation Research, Orleans, CA. Revised
More informationEstimating Probability of Success Rate
Estimating Probability of Success Rate Don Ferguson, PhD A member of Appalachian Search and Rescue Conference Objective Provide SAR planners with a tool to help optimize effort allocation based on a number
More informationIntroduction to risk assessment
Introduction to risk assessment Inception workshop of the project Strengthening the regional preparedness, prevention and response against lumpy skin disease in Belarus, Moldova and Ukraine (TCP/RER/3605)
More informationPRINCIPAL COMPONENT ANALYSIS OF MONTH-TO-MONTH PRECIPITATION VARIABILITY FOR NCDC CALIFORNIA CLIMATIC DIVISIONS, ( THROUGH SEASONS)
JP1.6 PRINCIPAL COMPONENT ANALYSIS OF MONTH-TO-MONTH PRECIPITATION VARIABILITY FOR NCDC CALIFORNIA CLIMATIC DIVISIONS, (1895-6 THROUGH 2000-1 SEASONS) Charles J. Fisk * NAWCWPNS-Point Mugu, CA 1. INTRODUCTION
More informationWeather and Climate Summary and Forecast March 2019 Report
Weather and Climate Summary and Forecast March 2019 Report Gregory V. Jones Linfield College March 2, 2019 Summary: Dramatic flip from a mild winter to a top five coldest February on record in many locations
More informationA Basic Introduction to Geographic Information Systems (GIS) ~~~~~~~~~~
A Basic Introduction to Geographic Information Systems (GIS) ~~~~~~~~~~ Rev. Ronald J. Wasowski, C.S.C. Associate Professor of Environmental Science University of Portland Portland, Oregon 3 September
More informationThe Wide Dynamic Range Vegetation Index and its Potential Utility for Gap Analysis
Summary StatMod provides an easy-to-use and inexpensive tool for spatially applying the classification rules generated from the CT algorithm in S-PLUS. While the focus of this article was to use StatMod
More informationQuality and Coverage of Data Sources
Quality and Coverage of Data Sources Objectives Selecting an appropriate source for each item of information to be stored in the GIS database is very important for GIS Data Capture. Selection of quality
More informationSeasonal and interannual relations between precipitation, soil moisture and vegetation in the North American monsoon region
Seasonal and interannual relations between precipitation, soil moisture and vegetation in the North American monsoon region Luis A. Mendez-Barroso 1, Enrique R. Vivoni 1, Christopher J. Watts 2 and Julio
More informationComparison of Imputation Procedures for Replacing Denied-access Plots
Comparison of Imputation Procedures for Replacing Denied-access Plots Susan L. King 1 Abstract. In forest inventories, missing plots are caused by hazardous terrain, inaccessible locations, or denied access.
More informationSpatial Effects on Current and Future Climate of Ipomopsis aggregata Populations in Colorado Patterns of Precipitation and Maximum Temperature
A. Kenney GIS Project Spring 2010 Amanda Kenney GEO 386 Spring 2010 Spatial Effects on Current and Future Climate of Ipomopsis aggregata Populations in Colorado Patterns of Precipitation and Maximum Temperature
More informationDelineation of high landslide risk areas as a result of land cover, slope, and geology in San Mateo County, California
Delineation of high landslide risk areas as a result of land cover, slope, and geology in San Mateo County, California Introduction Problem Overview This project attempts to delineate the high-risk areas
More informationEagle Creek Post Fire Erosion Hazard Analysis Using the WEPP Model. John Rogers & Lauren McKinney
Eagle Creek Post Fire Erosion Hazard Analysis Using the WEPP Model John Rogers & Lauren McKinney Columbia River Gorge at Risk: Using LiDAR and GIS-based predictive modeling for regional-scale erosion susceptibility
More informationCLIMATE AND LAND USE DRIVERS OF MALARIA RISK IN THE PERUVIAN AMAZON,
CLIMATE AND LAND USE DRIVERS OF MALARIA RISK IN THE PERUVIAN AMAZON, 2001-2012 Ben Zaitchik Johns Hopkins University Beth Feingold, Denis Valle, Alex Sandoval, Carlos Alvarez Antonio, Rosa Patricia Zegarra
More informationAn ENSO-Neutral Winter
An ENSO-Neutral Winter This issue of the Blue Water Outlook newsletter is devoted towards my thoughts on the long range outlook for winter. You will see that I take a comprehensive approach to this outlook
More informationEPIDEMIOLOGY FOR URBAN MALARIA MAPPING
TELE-EPIDEMIOLOGY EPIDEMIOLOGY FOR URBAN MALARIA MAPPING @IRD/M Dukhan Vanessa Machault Observatoire Midi-Pyrénées, Laboratoire d Aérologie Pleiades days 17/01/2012 The concept of Tele-epidemiology The
More informationKimberly J. Mueller Risk Management Solutions, Newark, CA. Dr. Auguste Boissonade Risk Management Solutions, Newark, CA
1.3 The Utility of Surface Roughness Datasets in the Modeling of United States Hurricane Property Losses Kimberly J. Mueller Risk Management Solutions, Newark, CA Dr. Auguste Boissonade Risk Management
More informationENVE203 Environmental Engineering Ecology (Nov 05, 2012)
ENVE203 Environmental Engineering Ecology (Nov 05, 2012) Elif Soyer Ecosystems and Living Organisms Population Density How Do Populations Change in Size? Maximum Population Growth Environmental Resistance
More informationPipeline Routing Using Geospatial Information System Analysis
Pipeline Routing Using Geospatial Information System Analysis Mahmoud Reza 1 Delavar and Fereydoon 2 Naghibi 1-Assistance Professor, Dept. of Surveying and Geomatic Eng., Eng. Faculty, University of Tehran,
More informationUse of Thematic Mapper Satellite Imagery, Hemispherical Canopy Photography, and Digital Stream Lines to Predict Stream Shading
Use of Thematic Mapper Satellite Imagery, Hemispherical Canopy Photography, and Digital Stream Lines to Predict Stream Shading David Nagel GIS Analyst Charlie Luce Research Hydrologist Bárbara Gutiérrez
More informationThe Current SLE/WN Epidemic Assesment
FMEL Arboviral Epidemic Risk Assessment: First Update for 2014 Week 18 (May 1, 2014) The Current SLE/WN Epidemic Assesment Funding for the Florida Medical Entomology Laboratory Epidemic Risk Model ended
More informationQuick Response Report #126 Hurricane Floyd Flood Mapping Integrating Landsat 7 TM Satellite Imagery and DEM Data
Quick Response Report #126 Hurricane Floyd Flood Mapping Integrating Landsat 7 TM Satellite Imagery and DEM Data Jeffrey D. Colby Yong Wang Karen Mulcahy Department of Geography East Carolina University
More informationReview Using the Geographical Information System and Remote Sensing Techniques for Soil Erosion Assessment
Polish J. of Environ. Stud. Vol. 19, No. 5 (2010), 881-886 Review Using the Geographical Information System and Remote Sensing Techniques for Soil Erosion Assessment Nuket Benzer* Landscape Architecture
More informationSouthwest LRT Habitat Analysis. May 2016 Southwest LRT Project Technical Report
Southwest LRT Habitat Analysis Southwest LRT Project Technical Report This page intentionally blank. Executive Summary This technical report describes the habitat analysis that was performed to support
More informationSampling The World. presented by: Tim Haithcoat University of Missouri Columbia
Sampling The World presented by: Tim Haithcoat University of Missouri Columbia Compiled with materials from: Charles Parson, Bemidji State University and Timothy Nyerges, University of Washington Introduction
More informationA Small Migrating Herd. Mapping Wildlife Distribution 1. Mapping Wildlife Distribution 2. Conservation & Reserve Management
A Basic Introduction to Wildlife Mapping & Modeling ~~~~~~~~~~ Rev. Ronald J. Wasowski, C.S.C. Associate Professor of Environmental Science University of Portland Portland, Oregon 8 December 2015 Introduction
More informationISO MODIS NDVI Weekly Composites for Canada South of 60 N Data Product Specification
ISO 19131 MODIS NDVI Weekly Composites for South of 60 N Data Product Specification Revision: A Data specification: MODIS NDVI Composites for South of 60 N - Table of Contents - 1. OVERVIEW... 3 1.1. Informal
More informationHarrison 1. Identifying Wetlands by GIS Software Submitted July 30, ,470 words By Catherine Harrison University of Virginia
Harrison 1 Identifying Wetlands by GIS Software Submitted July 30, 2015 4,470 words By Catherine Harrison University of Virginia cch2fy@virginia.edu Harrison 2 ABSTRACT The Virginia Department of Transportation
More informationP7: Limiting Factors in Ecosystems
P7: Limiting Factors in Ecosystems Purpose To understand that physical factors temperature and precipitation limit the growth of vegetative ecosystems Overview Students correlate graphs of vegetation vigor
More informationAssessing Drought in Agricultural Area of central U.S. with the MODIS sensor
Assessing Drought in Agricultural Area of central U.S. with the MODIS sensor Di Wu George Mason University Oct 17 th, 2012 Introduction: Drought is one of the major natural hazards which has devastating
More informationA case study for self-organized criticality and complexity in forest landscape ecology
Chapter 1 A case study for self-organized criticality and complexity in forest landscape ecology Janine Bolliger Swiss Federal Research Institute (WSL) Zürcherstrasse 111; CH-8903 Birmendsdorf, Switzerland
More informationBryan F.J. Manly and Andrew Merrill Western EcoSystems Technology Inc. Laramie and Cheyenne, Wyoming. Contents. 1. Introduction...
Comments on Statistical Aspects of the U.S. Fish and Wildlife Service's Modeling Framework for the Proposed Revision of Critical Habitat for the Northern Spotted Owl. Bryan F.J. Manly and Andrew Merrill
More informationANALYSIS AND VALIDATION OF A METHODOLOGY TO EVALUATE LAND COVER CHANGE IN THE MEDITERRANEAN BASIN USING MULTITEMPORAL MODIS DATA
PRESENT ENVIRONMENT AND SUSTAINABLE DEVELOPMENT, NR. 4, 2010 ANALYSIS AND VALIDATION OF A METHODOLOGY TO EVALUATE LAND COVER CHANGE IN THE MEDITERRANEAN BASIN USING MULTITEMPORAL MODIS DATA Mara Pilloni
More informationP6: Global Patterns in Green-Up and Green-Down
P6: Global Patterns in Green-Up and Green-Down Purpose To investigate the annual cycle of plant growth and decline using visualizations and graphs Overview Students will analyze visualizations and graphs
More informationOverview of The application of GIS and. Remote Sensing in Vector-Borne Disease Ecology
Steven Engborg NRS 509 December 14, 2017 Final Project Overview of The application of GIS and Remote Sensing in Vector-Borne Disease Ecology Anyone who learns how to create, manage, or use a Geographic
More informationProgress Report Year 2, NAG5-6003: The Dynamics of a Semi-Arid Region in Response to Climate and Water-Use Policy
Progress Report Year 2, NAG5-6003: The Dynamics of a Semi-Arid Region in Response to Climate and Water-Use Policy Principal Investigator: Dr. John F. Mustard Department of Geological Sciences Brown University
More informationModelling spatio-temporal patterns of disease
Modelling spatio-temporal patterns of disease Peter J Diggle CHICAS combining health information, computation and statistics References AEGISS Brix, A. and Diggle, P.J. (2001). Spatio-temporal prediction
More informationWeather and Climate Summary and Forecast January 2018 Report
Weather and Climate Summary and Forecast January 2018 Report Gregory V. Jones Linfield College January 5, 2018 Summary: A persistent ridge of high pressure over the west in December produced strong inversions
More informationLecture 9: Reference Maps & Aerial Photography
Lecture 9: Reference Maps & Aerial Photography I. Overview of Reference and Topographic Maps There are two basic types of maps? Reference Maps - General purpose maps & Thematic Maps - maps made for a specific
More informationMozambique. General Climate. UNDP Climate Change Country Profiles. C. McSweeney 1, M. New 1,2 and G. Lizcano 1
UNDP Climate Change Country Profiles Mozambique C. McSweeney 1, M. New 1,2 and G. Lizcano 1 1. School of Geography and Environment, University of Oxford. 2.Tyndall Centre for Climate Change Research http://country-profiles.geog.ox.ac.uk
More informationAn Introduction to Geographic Information System
An Introduction to Geographic Information System PROF. Dr. Yuji MURAYAMA Khun Kyaw Aung Hein 1 July 21,2010 GIS: A Formal Definition A system for capturing, storing, checking, Integrating, manipulating,
More informationUtility of National Spatial Data for Conservation Design Projects
Utility of National Spatial Data for Conservation Design Projects Steve Williams Biodiversity and Spatial Information Center North Carolina State University PIF CDW St. Louis, MO April 11, 2006 Types of
More informationClimate Outlook through 2100 South Florida Ecological Services Office Vero Beach, FL September 9, 2014
Climate Outlook through 2100 South Florida Ecological Services Office Vero Beach, FL September 9, 2014 Short Term Drought Map: Short-term (
More informationGlossary. The ISI glossary of statistical terms provides definitions in a number of different languages:
Glossary The ISI glossary of statistical terms provides definitions in a number of different languages: http://isi.cbs.nl/glossary/index.htm Adjusted r 2 Adjusted R squared measures the proportion of the
More informationMichigan State University, East Lansing, MI USA. Lansing, MI USA.
On-line Supporting Information for: Using Cost-Effective Targeting to Enhance the Efficiency of Conservation Investments in Payments for Ecosystem Services Xiaodong Chen1,*, Frank Lupi2, Andrés Viña1,
More informationSummary Description Municipality of Anchorage. Anchorage Coastal Resource Atlas Project
Summary Description Municipality of Anchorage Anchorage Coastal Resource Atlas Project By: Thede Tobish, MOA Planner; and Charlie Barnwell, MOA GIS Manager Introduction Local governments often struggle
More informationClimate Outlook through 2100 South Florida Ecological Services Office Vero Beach, FL January 13, 2015
Climate Outlook through 2100 South Florida Ecological Services Office Vero Beach, FL January 13, 2015 Short Term Drought Map: Short-term (
More informationModeling Disease Transmission in Long-tailed Macaques on Bali
Modeling Disease Transmission in Long-tailed Macaques on Bali Kelly Lane Gerhard Niederwieser Ryan Kennedy University of Notre Dame Macaque Background Coexisted in temples across the island for at least
More informationDevelopment of statewide 30 meter winter sage grouse habitat models for Utah
Development of statewide 30 meter winter sage grouse habitat models for Utah Ben Crabb, Remote Sensing and Geographic Information System Laboratory, Department of Wildland Resources, Utah State University
More informationTHE WILDLIFE SOCIETY CONFERENCE 2015
THE WILDLIFE SOCIETY CONFERENCE 2015 Evaluation of Desert Tortoise Habitat Suitability Model Predictions in the Soda Mountain Area, San Bernardino County Susanne Heim Introduction Overview of Species Habitat
More informationReport. Developing a course component on disaster management
Report Developing a course component on disaster management By Chira Prangkio Tawee Chaipimonplin Department of Geography, Faculty of Social Sciences, Chiang Mai University Thailand Presented at Indian
More informationGrant Opportunity Monitoring Bi-State Sage-grouse Populations in Nevada
Grant Opportunity Monitoring Bi-State Sage-grouse Populations in Nevada Proposals are due no later than November 13, 2015. Grant proposal and any questions should be directed to: Shawn Espinosa @ sepsinosa@ndow.org.
More information