Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Spatial Statistics”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 613 records · Page 34Linked to original sources
USGS1985

Trend analysis of weekly acid rain data, 1978-83

There are 19 stations in the National Atmospheric Deposition Program which operated over the period 1978-83 and were subsequently incorporated into the National Trends Network in 1983. The precipitation chemistry data for these stations for this period were analyzed for trend, spatial correlation, seasonality, and relationship to precipitation volume. The intent of the analysis was to provide insights on the sources of variation in precipitation chemistry and to attempt to ascertain what statistical procedures may be most useful for ongoing analysis of the National Trends Network data. The Seasonal Kendall test was used for detection of trends in raw concentrations of dissolved constituents, pH and specific conductance, and residuals of these parameters from regression analysis. Forty-one percent of the trends detected in the raw concentrations were downtrends, 4 percent were uptrends, and 55 percent showed no trends at a = 0.2. At a more restrictive significance level of a = 0.05, 24 percent of the trends detected were downtrends, 2 percent were uptrends, and 74 percent showed no trends. The two constituents of greatest interest in terms of human generated emissions and environmental effects, sulfate and nitrate, showed only downtrends, and sulfate showed the largest decreases in concentration per year of all the ions tested.

Water-Resources Investigations Report↗
USGS1985

Estimating neighborhood variability with a binary comparison matrix.

A technique which utilizes a binary comparison matrix has been developed to implement a neighborhood function for a raster format data base. The technique assigns an index value to the center pixel of 3- by 3-pixel neighborhoods. The binary comparison matrix provides additional information not found in two other neighborhood variability statistics; the function is sensitive to both the number of classes within the neighborhood and the frequency of pixel occurrence in each of the classes. Application of the function to a spatial data base from the Kenai National Wildlife Refuge, Alaska, demonstrates 1) the numerical distribution of the index values, and 2) the spatial patterns exhibited by the numerical values. -Author

Photogrammetric Engineering and Remote Sensing↗
USGS2022

A hydrologic perspective of major U.S. droughts

Drought is a recurring natural hazard that has substantial human and environmental impacts. Given continued global warming and associated climate change, there is concern that droughts could become more severe and longer lasting. To better monitor and understand drought development and persistence, it is helpful to understand the development and climatic drivers of past droughts. In this study we use monthly runoff percentiles to identify five major drought events in the conterminous United States (CONUS) from 1901 through 2020. For each drought event we examined spatial patterns of departures of mean monthly precipitation, temperature, soil moisture storage, and runoff for 2,107 hydrologic units (HUs) across the CONUS. Results indicated that precipitation deficits have been the primary driver of past major-drought events and temperature a secondary driver, even of the most recent drought event (September 1999 through September 2015) when positive temperature anomalies occurred over most of the CONUS. Additionally, negative soil moisture storage departures were more negative than runoff departures during the five drought events we examined, which emphasizes the importance of measuring both runoff and soil moisture to monitor drought conditions. We also examined the use of statistical persistence to develop short-term (i.e., 1 month) forecasts of runoff drought conditions in the CONUS by developing autoregressive integrated moving average (ARIMA) models for each HU. Results indicated that persistence can be used to predict short-term changes in the spatial pattern of drought and the areal extent of drought, but that predictions of runoff magnitude for any particular site are often poor.

International Journal of Climatology↗
USGS2020

A modeling workflow that balances automation and human intervention to inform invasive plant management decisions at multiple spatial scales

Predictions of habitat suitability for invasive plant species can guide risk assessments at regional and national scales and inform early detection and rapid-response strategies at local scales. We present a general approach to invasive species modeling and mapping that meets objectives at multiple scales. Our methodology is designed to balance trade-offs between developing highly customized models for few species versus fitting non-specific and generic models for numerous species. We developed a national library of environmental variables known to physiologically limit plant distributions and relied on human input based on natural history knowledge to further narrow the variable set for each species before developing habitat suitability models. To ensure efficiency, we used largely automated modeling approaches and human input only at key junctures. We explore and present uncertainty by using two alternative sources of background samples, including five statistical algorithms, and constructing model ensembles. We demonstrate the use and efficiency of the Software for Assisted Habitat Modeling [SAHM 2.1.2], a package in VisTrails, which performs the majority of the modeling analyses. Our workflow includes solicitation of expert feedback on model outputs such as spatial prediction results and variable response curves, and iterative improvement based on new data availability and directed field validation of initial model results. We highlight the utility of the models for decision-making at regional and local scales with case studies of two plant species that invade natural areas: fountain grass ( Pennisetum setaceum ) and goutweed ( Aegopodium podagraria ). By balancing model automation with human intervention, we can efficiently provide land managers with mapped predicted distributions for multiple invasive species to inform decisions across spatial scales.

PLoS ONE↗
USGS2019

Global phylodynamic analysis of avian paramyxovirus-1 provides evidence of inter-host transmission and intercontinental spatial diffusion

Background Avian avulavirus (commonly known as avian paramyxovirus-1 or APMV-1) can cause disease of varying severity in both domestic and wild birds. Understanding how viruses move among hosts and geography would be useful for informing prevention and control efforts. A Bayesian statistical framework was employed to estimate the evolutionary history of 1602 complete fusion gene APMV-1 sequences collected from 1970 to 2016 in order to infer viral transmission between avian host orders and diffusion among geographic regions. Ancestral states were estimated with a non-reversible continuous-time Markov chain model, allowing transition rates between discrete states to be calculated. The evolutionary analyses were stratified by APMV-1 classes I ( n = 198) and II ( n = 1404), and only those sequences collected between 2006 and 2016 were allowed to contribute host and location information to the viral migration networks. Results While the current data was unable to assess impact of host domestication status on APMV-1 diffusion, these analyses supported the sharing of APMV-1 among divergent host taxa. The highest supported transition rate for both classes existed from domestic chickens to Anseriformes (class I:6.18 transitions/year, 95% highest posterior density (HPD) 0.31–20.02, Bayes factor (BF) = 367.2; class II:2.88 transitions/year, 95%HPD 1.9–4.06, BF = 34,582.9). Further, among class II viruses, domestic chickens also acted as a source for Columbiformes (BF = 34,582.9), other Galliformes (BF = 34,582.9), and Psittaciformes (BF = 34,582.9). Columbiformes was also a highly supported source to Anseriformes (BF = 322.0) and domestic chickens (BF = 402.6). Additionally, our results provide support for the diffusion of viruses among continents and regions, but no interhemispheric viral exchange between 2006 and 2016. Among class II viruses, the highest transition rates were estimated from South Asia to the Middle East (1.21 transitions/year; 95%HPD 0.36–2.45; BF = 67,107.8), from Europe to East Asia (1.17 transitions/year; 95%HPD 0.12–2.61; BF = 436.2) and from Europe to Africa (1.06 transitions/year, 95%HPD 0.07–2.51; BF = 169.3). Conclusions While migration appears to occur infrequently, geographic movement may be important in determining viral diversification and population structure. In contrast, inter-order transmission of APMV-1 may occur readily, but most events are transient with few lineages persisting in novel hosts.

BMC Evolutionary Biology↗
USGS2007

Spatial design and strength of spatial signal: Effects on covariance estimation

In a spatial regression context, scientists are often interested in a physical interpretation of components of the parametric covariance function. For example, spatial covariance parameter estimates in ecological settings have been interpreted to describe spatial heterogeneity or “patchiness” in a landscape that cannot be explained by measured covariates. In this article, we investigate the influence of the strength of spatial dependence on maximum likelihood (ML) and restricted maximum likelihood (REML) estimates of covariance parameters in an exponential-with-nugget model, and we also examine these influences under different sampling designs—specifically, lattice designs and more realistic random and cluster designs—at differing intensities of sampling ( n =144 and 361). We find that neither ML nor REML estimates perform well when the range parameter and/or the nugget-to-sill ratio is large—ML tends to underestimate the autocorrelation function and REML produces highly variable estimates of the autocorrelation function. The best estimates of both the covariance parameters and the autocorrelation function come under the cluster sampling design and large sample sizes. As a motivating example, we consider a spatial model for stream sulfate concentration.

Journal of Agricultural, Biological, and Environme↗
USGS2007

Wave-driven spatial and temporal variability in sea-floor sediment mobility in the Monterey Bay, Cordell Bank, and Gulf of the Farallones National Marine Sanctuaries

Wind and wave patterns affect many aspects of continental shelves and shorelines geomorphic evolution. Although our understanding of the processes controlling sediment suspension on continental shelves has improved over the past decade, our ability to predict sediment mobility over large spatial and temporal scales remains limited. The deployment of robust operational buoys along the U.S. West Coast in the early 1980s provides large quantities of high-resolution oceanographic and meteorologic data. By 2006, these data sets were long enough to clearly identify long-term trends and compute statistically significant probability estimates of wave and wind behavior during annual and interannual climatic cycles (that is, El Niño and La Niña). Wave-induced sediment mobility on the shelf and upper slope off central California was modeled using synthesized oceanographic and meteorologic data as boundary input for the Delft SWAN model, sea-floor grain-size data provided by the usSEABED database, and regional bathymetry. Differences in waves (heights, periods, and directions) and winds (speeds and directions) between El Niño and La Niña months cause temporal and spatial variations in peak wave-induced bed shear stresses. These variations, in conjunction with spatially heterogeneous unconsolidated sea-floor sedimentary cover, result in predicted sediment mobility widely varying in both time and space. These findings indicate that these factors have significant consequences for both geological and biological processes.

California↗
USGS2009

Assessing rates of forest change and fragmentation in Alabama, USA, using the vegetation change tracker model

Forest change is of great concern for land use decision makers and conservation communities. Quantitative and spatial forest change information is critical for addressing many pressing issues, including global climate change, carbon budgets, and sustainability. In this study, our analysis focuses on the differences in geospatial patterns and their changes between federal forests and nonfederal forests in Alabama over the time period 1987–2005, by interpreting 163 Landsat Thematic Mapper (TM) scenes using a vegetation change tracker (VCT) model. Our analysis revealed that for the most part of 1990 s and between 2000 and 2005, Alabama lost about 2% of its forest on an annual basis due to disturbances, but much of the losses were balanced by forest regeneration from previous disturbances. The disturbance maps revealed that federal forests were reasonably well protected, with the fragmentation remaining relatively stable over time. In contrast, nonfederal forests, which are predominant in area share (about 95%), were heavily disturbed, clearly demonstrating decreasing levels of fragmentation during the time period 1987–1993 giving way to a subsequent accelerating fragmentation during the time period 1994–2005. Additionally, the identification of the statistical relationships between forest fragmentation status and forest loss rate and forest net change rate in relation to land ownership implied the distinct differences in forest cutting rate and cutting patterns between federal forests and nonfederal forests. The forest spatial change information derived from the model has provided valuable insights regarding regional forest management practices and disturbance regimes, which are closely associated with regional economics and environmental concerns.

Forest Ecology and Management↗
USGS2020

The predictive skills of elastic Coulomb rate-and-state aftershock forecasts during the 2019 Ridgecrest, California, earthquake sequence

Operational earthquake forecasting protocols commonly use statistical models for their recognized ease of implementation and robustness in describing the short-term spatiotemporal patterns of triggered seismicity. However, recent advances on physics-based aftershock forecasting reveal comparable performance to the standard statistical counterparts with significantly improved predictive skills when fault and stress field heterogeneities are considered. Here, we perform a pseudo-prospective forecasting experiment during the first month of the 2019 Ridgecrest (California) earthquake sequence. We develop seven Coulomb rate-and-state models that couple static stress change estimates with continuum mechanics expressed by the rate-and-state friction laws. Our model parametrization supports a gradually increasing complexity; we start from a preliminary model implementation with simplified slip distributions and spatially homogeneous receiver faults to reach an enhanced one featuring optimized fault constitutive parameters, finite-fault slip models, secondary triggering effects, and spatially heterogenous planes informed by pre-existing ruptures. The data-rich environment of Southern California allows us to test whether incorporating data collected in near real-time during an unfolding earthquake sequence boosts our predictive power. We assess the absolute and relative performance of the forecasts by means of statistical tests used within the Collaboratory for the Study of Earthquake Predictability (CSEP) and compare their skills against a standard benchmark ETAS model for the short (24 hours after the two Ridgecrest mainshocks) and intermediate-term (one month). Stress-based forecasts expect heightened rates along the whole near-fault region and increased expected seismicity rates in Central Garlock Fault. Our comparative model evaluation supports that faulting heterogeneities coupled with secondary triggering effects are the most critical success components behind physics-based forecasts, but also underlines the importance of model updates incorporating near real-time available aftershock data reaching better performance than ETAS models. We explore the physical basis behind our results by investigating the localized shut down of pre-existing normal faults in the Ridgecrest near-source area.

California↗
USGS2012

Ecosystem engineering varies spatially: a test of the vegetation modification paradigm for prairie dogs

Colonial, burrowing herbivores can be engineers of grassland and shrubland ecosystems worldwide. Spatial variation in landscapes suggests caution when extrapolating single-place studies of single species, but lack of data and the need to generalize often leads to ‘model system’ thinking and application of results beyond appropriate statistical inference. Generalizations about the engineering effects of prairie dogs ( Cynomys sp.) developed largely from intensive study at a single complex of black-tailed prairie dogs C. ludovicianus in northern mixed prairie, but have been extrapolated to other ecoregions and prairie dog species in North America, and other colonial, burrowing herbivores. We tested the paradigm that prairie dogs decrease vegetation volume and the cover of grasses and tall shrubs, and increase bare ground and forb cover. We sampled vegetation on and off 279 colonies at 13 complexes of 3 prairie dog species widely distributed across 5 ecoregions in North America. The paradigm was generally supported at 7 black-tailed prairie dog complexes in northern mixed prairie, where vegetation volume, grass cover, and tall shrub cover were lower, and bare ground and forb cover were higher, on colonies than at paired off-colony sites. Outside the northern mixed prairie, all 3 prairie dog species consistently reduced vegetation volume, but their effects on cover of plant functional groups varied with prairie dog species and the grazing tolerance of dominant perennial grasses. White-tailed prairie dogs C. leucurus in sagebrush steppe did not reduce shrub cover, whereas black-tailed prairie dogs suppressed shrub cover at all complexes with tall shrubs in the surrounding habitat matrix. Black-tailed prairie dogs in shortgrass steppe and Gunnison's prairie dogs C. gunnisoni in Colorado Plateau grassland both had relatively minor effects on grass cover, which may reflect the dominance of grazing-tolerant shortgrasses at both complexes. Variation in modification of vegetation structure may be understood in terms of the responses of different dominant perennial grasses to intense defoliation and differences in foraging behavior among prairie dog species. Spatial variation in the engineering role of prairie dogs suggests spatial variation in their keystone role, and spatial variation in the roles of other ecosystem engineers. Thus, ecosystem engineering can have a spatial component not evident from single-place studies.

Ecography: Pattern and Diversity in Ecology↗
USGS2009

A landscape indicator approach to the identification and articulation of the consequences of land-cover change in the Mid-Atlantic Region, 1973-2001

Landscape indicators, derived from land-use and land-cover data, hydrology, nitrate deposition, and elevation data, were used by Jones and others (2001a) to calculate the ecological consequences of land-cover change. Nitrate loading and physical bird habitat were modeled from 1973 and 1992 land-cover and other spatial data for the Mid-Atlantic region. Utilizing the same methods, this study extends the analysis another decade with the use of the 2001 National Land Cover Dataset. Land-cover statistics and trends are calculated for three time periods: 1973-1992, 1992-2001 and 1973-2001. In addition, high-resolution aerial photographs (1 meter or better ground-sample distance) were acquired and analyzed for thirteen pairs of adjacent USGS 7.5 minute quadrangle maps in areas where distinct positive or negative changes to nitrogen loading and bird habitat were previously calculated. During the entire 30 year period, the data show that there was extensive loss of agriculture and forest area and a major increase in urban land-cover classes. However, the majority of the conversion of other classes to urban occurred during the 1992-2001 period. During the 1973-1992 period, there was only moderate increase in urban area, while there was an inverse relationship between agricultural change and forest change. In general, forest gain and agricultural loss was found in areas of improving landscape indicators, and forest loss and agricultural gain was found to occur in areas of declining indicators related to habitat and nitrogen loadings, which was generally confirmed by the aerial photographic analysis. In terms of the specific model results, bird habitat, which is mainly related to the extent of forest cover, declined overall with forest extent, but was also affected more in the decline of habitat quality. Nitrate loading, which is mainly related to agricultural land cover actually improved from 1992-2001, and in the overall study, mainly due to the conversion of agriculture to forests and urban. The high-resolution imagery analysis was significant in that it confirmed, at a very local level, the specific land-cover changes that were driving the landscape metrics and model results that were calculated from moderate resolution land-cover data and models. These were generally subtle changes in patch size of agriculture, forest, and urban areas, but had substantial effects on bird habitat and nitrogen loadings. This analysis of high-resolution imagery demonstrates and confirms the important ability of moderate-resolution land-cover data to capture significant landscape-level activity that is directly related to specific metrics of ecological significance. It also demonstrates consistent landscape-scale relationships between data derived from high-resolution, moderate-resolution and landscape-model sources. Finally, many of the areas of improvement and decline in bird habitat and nitrogen loadings appear to be potentially regional in nature and likely reflect some local trend in landscape activity. Although the use of ecoregions as sampling units has been criticized in recent years, these results show that basic changes in Level 1 land-cover categories, such as forest and agriculture, may still reflect ecoregional patterns and considerations at some scale of mapping and analysis. This is a potentially important area for future landscape-indicator research. This and other follow-on research opportunities are discussed.

Open-File Report↗
USGS2015

Machine learning for predicting soil classes in three semi-arid landscapes

Mapping the spatial distribution of soil taxonomic classes is important for informing soil use and management decisions. Digital soil mapping (DSM) can quantitatively predict the spatial distribution of soil taxonomic classes. Key components of DSM are the method and the set of environmental covariates used to predict soil classes. Machine learning is a general term for a broad set of statistical modeling techniques. Many different machine learning models have been applied in the literature and there are different approaches for selecting covariates for DSM. However, there is little guidance as to which, if any, machine learning model and covariate set might be optimal for predicting soil classes across different landscapes. Our objective was to compare multiple machine learning models and covariate sets for predicting soil taxonomic classes at three geographically distinct areas in the semi-arid western United States of America (southern New Mexico, southwestern Utah, and northeastern Wyoming). All three areas were the focus of digital soil mapping studies. Sampling sites at each study area were selected using conditioned Latin hypercube sampling (cLHS). We compared models that had been used in other DSM studies, including clustering algorithms, discriminant analysis, multinomial logistic regression, neural networks, tree based methods, and support vector machine classifiers. Tested machine learning models were divided into three groups based on model complexity: simple, moderate, and complex. We also compared environmental covariates derived from digital elevation models and Landsat imagery that were divided into three different sets: 1) covariates selected a priori by soil scientists familiar with each area and used as input into cLHS, 2) the covariates in set 1 plus 113 additional covariates, and 3) covariates selected using recursive feature elimination. Overall, complex models were consistently more accurate than simple or moderately complex models. Random forests (RF) using covariates selected via recursive feature elimination was consistently the most accurate, or was among the most accurate, classifiers between study areas and between covariate sets within each study area. We recommend that for soil taxonomic class prediction, complex models and covariates selected by recursive feature elimination be used. Overall classification accuracy in each study area was largely dependent upon the number of soil taxonomic classes and the frequency distribution of pedon observations between taxonomic classes. Individual subgroup class accuracy was generally dependent upon the number of soil pedon observations in each taxonomic class. The number of soil classes is related to the inherent variability of a given area. The imbalance of soil pedon observations between classes is likely related to cLHS. Imbalanced frequency distributions of soil pedon observations between classes must be addressed to improve model accuracy. Solutions include increasing the number of soil pedon observations in classes with few observations or decreasing the number of classes. Spatial predictions using the most accurate models generally agree with expected soil–landscape relationships. Spatial prediction uncertainty was lowest in areas of relatively low relief for each study area.

New Mexico, Utah, Wyoming↗
USGS2019

Climate, streamflow, and lake-level trends in the Great Lakes Basin of the United States and Canada, water years 1960–2015

Water levels in the Great Lakes fluctuate substantially because of complex interactions among inputs (precipitation and streamflow), outputs (evaporation and outflow), and other factors. This report by the U.S. Geological Survey in cooperation with the Great Lakes Restoration Initiative was completed to describe trends in climate, streamflow, lake levels, and major water-budget components within the Great Lakes Basin for water years (WYs) 1960–2015 (study period). Resulting trends are applicable only to the study period and should not be considered indicative of longer-term trends. Analyses of climate trends used monthly data from the Parameter-elevation Regressions on Independent Slopes Model, which are available only for the United States. Trend tests were completed for annual and seasonal time series of monthly means for total precipitation, daily minimum air temperature ( T min ), and daily maximum air temperature ( T max ). Statistical significance for all time-trend tests (climate, streamflow, and lake levels) was determined using the Mann‑Kendall test for probability values less than or equal to 0.10. Trend analyses were completed without adjustments for serial correlation; however, a modified Mann-Kendall test was subsequently used to examine potential effects of short-term persistence in time-series data. Effects of short-term persistence were considered inconsequential for climate data and minor for streamflow data; however, the presence of short-term persistence in water-budget components had more substantial effects on trend analyses. Spatial distributions of trends in climatic data for WYs 1960–2015 for the U.S. part of the Great Lakes Basin (land only) indicate (1) generally ubiquitous upward trends in T min and (2) a sharp transition from neutral or downward trends in precipitation northwest of Lake Michigan to generally upward trends east of Lake Michigan. Trends in T max were not statistically significant. Analyses of annual climatic data aggregated for the U.S. land part of the Great Lakes Basin indicated statistically significant upward trends for precipitation and T min , and similar statistically significant trends existed for all the individual lake subbasins except Lake Superior. Of 103 U.S. Geological Survey streamgages analyzed for streamflow trends, 71 had significant annual trends (54 upward and 17 downward). Downward trends in annual streamflow are concentrated northwest of Lake Michigan (16 streamgages), and upward trends are concentrated east of Lake Michigan (53 streamgages). Of the 71 streamgages with significant annual trends, 70 had at least one season with a significant trend that matched the annual trend direction. Of 35 Environment and Climate Change Canada streamgages analyzed, 22 had significant upward trends in annual streamflow, and all but 1 of these 22 had at least one season with a significant upward trend. None of the Environment and Climate Change Canada streamgages had significant downward annual trends, and only one had a significant downward seasonal trend. Trends in lake levels and several major water-budget components affecting lake levels were analyzed for the study period. Significant downward trends in lake level and outflow for Lake Superior are driven primarily by low lake levels and outflows during WYs 1998–2014. A significant downward trend in runoff from the contributing drainage area also is indicated, which is consistent with numerous streamgages northwest of Lake Michigan with significant downward trends in annual streamflow. A significant upward trend in annual overlake evaporation also is indicated, which is consistent with the spatially distributed upward trends in annual T min . The sum of overlake precipitation and runoff from the contributing drainage area for each of the Great Lakes, less overlake evaporation, composes a variable called net basin supply (NBS). A significant downward trend in NBS is indicated for Lake Superior, which is consistent with significant trends for individual components of runoff (downward) and evaporation (upward) that contributed to a significant downward trend for lake outflow. Statistically significant upward trends in NBS for Lake Saint Clair and Lake Ontario offset the downward trend for Lake Superior and combine with nonsignificant upward trends in NBS for Lakes Michigan and Huron and Lake Erie to produce a neutral trend in NBS for the basin. A predictable pattern in monthly mean lake levels is noted for Lake Superior, with the minimum for each year usually during or near March and the maximum commonly during or near September or October. When an October lake level is in a period of substantial decline, potential for an ensuing short-term period of below-mean lake levels is enhanced. Downstream from Lake Superior, monthly lake levels have sawtooth patterns that somewhat resemble those for Lake Superior but with decreased predictability in timing. Similar to Lake Superior, Lakes Michigan and Huron, Lake Saint Clair, and Lake Erie all have a prolonged period of low lake levels around WYs 1998–2014; however, a significant downward trend is indicated only for Lakes Michigan and Huron. All these lakes also have a period of low lake levels before about WY 1968, when minimum lake levels were lower than during WYs 1998–2014. The significant downward trend of outflow from Lake Superior is carried downstream into Lakes Michigan and Huron; however, trends in outflow from the next three lakes downstream (Lakes Saint Clair, Erie, and Ontario) are offset by increased precipitation and runoff and are not significant.

Scientific Investigations Report↗
USGS2007

Landscape factors influencing the spatial distribution and abundance of mosquito vector Culex quinquefasciatus (Diptera: Culicidae) in a mixed residential-agricultural community in Hawai'i

Mosquito-borne avian diseases, principally avian malaria (Plasmodium relictum Grassi and Feletti) and avian pox (Avipoxvirus sp.) have been implicated as the key limiting factor associated with recent declines of endemic avifauna in the Hawaiian Island archipelago. We present data on the relative abundance, infection status, and spatial distribution of the primary mosquito vector Culex quinquefasciatus Say (Diptera: Culicidae) across a mixed, residential-agricultural community adjacent to Hawai'i Volcanoes National Park on Hawai'i Island. We modeled the effect of agriculture and forest fragmentation in determining relative abundance of adult Cx. quinquefasciatus in Volcano Village, and we implement our statistical model in a geographic information system to generate a probability of mosquito capture prediction surface for the study area. Our model was based on biweekly captures of adult mosquitoes from 20 locations within Volcano Village from October 2001 to April 2003. We used mixed effects logistic regression to model the probability of capturing a mosquito, and we developed a set of 17 competing models a priori to specifically evaluate the effect of agriculture and fragmentation (i.e., residential landscapes) at two spatial scales. In total, 2,126 mosquitoes were captured in CO 2-baited traps with an average probability of 0.27 (SE = 0.10) of capturing one or more mosquitoes per trap night. Twelve percent of mosquitoes captured were infected with P. relictum. Our data indicate that agricultural lands and forest fragmentation significantly increase the probability of mosquito capture. The prediction surface identified areas along the Hawai'i Volcanoes National Park boundary that may have high relative abundance of the vector. Our data document the potential of avian malaria transmission in residential-agricultural landscapes and support the need for vector management that extends beyond reserve boundaries and considers a reserve's spatial position in a highly heterogeneous landscape.

Journal of Medical Entomology↗
USGS2017

Modeling waterfowl habitat selection in the Central Valley of California to better understand the spatial relationship between commercial poultry and waterfowl

Wildlife researchers frequently study resource and habitat selection of wildlife to understand their potential habitat requirements and to conserve their populations. Understanding wildlife spatial-temporal distributions related to habitat have other applications such as to model interfaces between wildlife and domestic food animals in order to mitigate disease transmission to food animals. The highly pathogenic avian influenza (HPAI) virus represents a significant risk to the poultry industry. The Central Valley of California offers a unique geographical confluence of commercial poultry and wild waterfowl, which are thought to be a key reservoir of avian influenza (AI). Therefore, understanding spatio-temporal distributions of waterfowl could improve our understanding of potential risk of HPAI exposure from a commercial poultry perspective. Using existing radio-telemetry data on waterfowl (U.S. Geological Survey) in combination with habitat and vegetation data based on Geographic Information Systems (GIS), we are developing GIS-based statistical models that predict the probability of waterfowl presence (Habitat Suitability Mapping). Near-real-time application can be developed using recent habitat data derived from Landsat imagery (acquired by satellites and publicly available through the U.S. Geological Survey) to predict temporally- and spatially-varying distributions of waterfowl in the Central Valley. These results could be used to provide decision support for the poultry industry in addressing potential risk of HPAI exposure related to waterfowl proximity.

California↗
USGS2014

Methods for estimating magnitude and frequency of floods in Arizona, developed with unregulated and rural peak-flow data through water year 2010

Flooding is among the worst natural disasters responsible for loss of life and property in Arizona, underscoring the importance of accurate estimation of flood magnitude for proper structural design and floodplain mapping. Twenty-four years of additional peak-flow data have been recorded since the last comprehensive regional flood frequency analysis conducted in Arizona. Periodically, flood frequency estimates and regional regression equations must be revised to maintain the accurate estimation of flood frequency and magnitude. Annual peak-flow data collected through water year 2010 were compiled from 448 unregulated streamflow-gaging stations, hereafter referred to as streamgages, in Arizona having a minimum of 10 years of record. Flood frequency estimates were first computed with station (or at-site) skew using the Expected Moments Algorithm with a multiple Grubbs-Beck test to identify multiple potentially influential low flows to fit a Pearson Type III distribution. Next, a multiple step Bayesian least-squares-regression approach was used to determine a new statewide regional skew of −0.09. No basin characteristics analyzed were statistically significant in explaining the variation in skew and as a result, the constant model was chosen as the best regional skew model for the Arizona study area. The mean square error used in Bulletin 17B (B17B) of the Interagency Advisory Committee on Water Data is used to describe the precision of the regional skew. The constant model had a mean square error equal to 0.08, which corresponds to an effective record length of 85 years. This is a marked improvement over a previous Arizona regional skew analysis, with a reported mean square error of 0.31, for a corresponding effective record length of around 17 years. Thus the new regional model had almost five times the information content (as measured by effective record length) of that calculated in USGS Water Supply Paper 2433, published in 1997, or the value of 0.302 reported in the B17B generalized skew map. The flood frequency estimates were recalculated using a weighted skew of the station and regional skew. Station flood frequency estimates for each streamgage are presented for the 50-, 20-, 10-, 4-, 2-, 1-, 0.5-, and 0.2-percent annual exceedance probabilities. Geographical information systems were used to compute basin characteristic information for each streamgage for the purpose of developing regional equations to estimate flood statistics at ungaged basins. Five hydrologic flood regions in Arizona were defined in a multivariate regionalization process based on mean basin elevation, mean annual precipitation, and soil permeability. A regional generalized least-squares-regression analysis was used to develop five sets of equations from 344 nonredundant streamgages, corresponding to five regions, for estimating the 50-, 20-, 10-, 4-, 2-, 1-, 0.5-, and 0.2-percent annual exceedance probabilities at ungaged basins in Arizona. The regression equations developed for these five regions were based on one or more of the statistically significant explanatory variables: drainage area, mean basin elevation, and mean annual precipitation. Average standard errors of prediction for the regression regions for the five regions ranged from 27 to 122 percent and the pseudo-coefficients of determination (pseudo-R 2 ), a measure of the proportion of peak-flow variation that is explained by the basin characteristics, ranged from 68 to 98 percent. Regression equations for Central Highlands (region 4) had the lowest model error and the greatest pseudo-R 2 metrics. The equations for Colorado Plateau (region 2) regression equations generally had greater model error and lower pseudo-R 2 metrics. The improvement of regional regression equation model error and pseudo-R 2 metrics was related to higher numbers of streamgages, longer period of record, and even spatial coverage within a region. The regional regression equations were integrated into the U.S. Geological Survey’s StreamStats program. The StreamStats program is a national map-based web application that allows the public to easily access published flood frequency and basin characteristic statistics. The interactive web application allows a user to select a point within a watershed (gaged or ungaged) and retrieve flood-frequency estimates derived from the current regional regression equations and geographic information system data within the selected basin. StreamStats provides users with an efficient and accurate means for retrieving the most up to date flood frequency and basin characteristic data. StreamStats is intended to provide consistent statistics, minimize user error, and reduce the need for large datasets and costly geographic information system software.

Arizona↗
USGS2011

Influences of the El Niño Southern Oscillation and the Pacific Decadal Oscillation on the timing of the North American spring

Detrended, modelled first leaf dates for 856 sites across North America for the period 1900–2008 are used to examine how the El Niño Southern Oscillation (ENSO) and the Pacific Decadal Oscillation (PDO) separately and together might influence the timing of spring. Although spring (mean March through April) ENSO and PDO signals are apparent in first leaf dates, the signals are not statistically significant (at a 95% confidence level (p < 0.05)) for most sites. The most significant ENSO/PDO signal in first leaf dates occurs for El Niño and positive PDO conditions. An analysis of the spatial distributions of first leaf dates for separate and combined ENSO/PDO conditions features a northwest–southeast dipole that is significantly (at p < 0.05) different than the distributions for neutral conditions. The nature of the teleconnection between Pacific SST's and first leaf dates is evident in comparable composites for detrended sea level pressure (SLP) in the spring months. During positive ENSO/PDO, there is an anomalous flow of warm air from the southwestern US into the northwestern US and an anomalous northeasterly flow of cold air from polar regions into the eastern and southeastern US. These flow patterns are reversed during negative ENSO/PDO. Although the magnitudes of first leaf date departures are not necessarily significantly related to ENSO and PDO, the spatial patterns of departures are significantly related to ENSO and PDO. These significant relations and the long-lived persistence of SSTs provide a potential tool for forecasting the tendencies for first leaf dates to be early or late.

International Journal of Climatology↗