Search USGSSearch

SEARCH · Search USGS

Results for “Statistical Science”

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 451 records · Page 25Linked to original sources

Nitrate in the Mississippi River and its tributaries, 1980 to 2008: Are we making progress?

Changes in nitrate concentration and flux between 1980 and 2008 at eight sites in the Mississippi River basin were determined using a new statistical method that accommodates evolving nitrate behavior over time and produces flow-normalized estimates of nitrate concentration and flux that are independent of random variations in streamflow. The results show that little consistent progress has been made in reducing riverine nitrate since 1980, and that flow-normalized concentration and flux are increasing in some areas. Flow-normalized nitrate concentration and flux increased between 9 and 76% at four sites on the Mississippi River and a tributary site on the Missouri River, but changed very little at tributary sites on the Ohio, Iowa, and Illinois Rivers. Increases in flow-normalized concentration and flux at the Mississippi River at Clinton and Missouri River at Hermann were more than three times larger than at any other site. The increases at these two sites contributed much of the 9% increase in flow-normalized nitrate flux leaving the Mississippi River basin. At most sites, concentrations increased more at low and moderate streamflows than at high streamflows, suggesting that increasing groundwater concentrations are having an effect on river concentrations.

Illinois, Iowa, Louisiana, Missouri

Velocity structure of a bottom simulating reflector offshore Peru: Results from full waveform inversion

Much of our knowledge of the worldwide distribution of submarine gas hydrates comes from seismic observations of Bottom Simulating Reflectors (BSRs). Full waveform inversion has proven to be a reliable technique for studying the fine structure of BSRs using the compressional wave velocity. We applied a non-linear full waveform inversion technique to a BSR at a location offshore Peru. We first determined the large-scale features of seismic velocity variations using a statistical inversion technique to maximise coherent energy along travel-time curves. These velocities were used for a starting velocity model for the full waveform inversion, which yielded a detailed velocity/depth model in the vicinity of the BSR. We found that the data are best fit by a model in which the BSR consists of a thin, low-velocity layer. The compressional wave velocity drops from 2.15 km/s down to an average of 1.70 km/s in an 18m thick interval, with a minimum velocity of 1.62 km/s in a 6 m interval. The resulting compressional wave velocity was used to estimate gas content in the sediments. Our results suggest that the low velocity layer is a 6-18 m thick zone containing a few percent of free gas in the pore space. The presence of the BSR coincides with a region of vertical uplift. Therefore, we suggest that gas at this BSR is formed by a dissociation of hydrates at the base of the hydrate stability zone due to uplift and subsequently a decrease in pressure.

Earth and Planetary Science Letters

Estimates of diffuse phosphorus sources in surface waters of the United States using a spatially referenced watershed model

The statistical watershed model SPARROW (SPAtially Referenced Regression On Watershed attributes) was used to estimate the sources and transport of total phosphorus (TP) in surface waters of the United States. We calibrated the model using stream measurements of TP from 336 watersheds of mixed land use and spatial data on topography, soils, stream hydrography, and land use (agriculture, forest, shrub/grass, urban). The model explained 87% of the spatial variability in log transformed stream TP flux (kg yr -1 ). Predictions of stream yield ( kg ha -1 yr -1 ) were typically within 45% of the observed values at the monitoring sites. The model identified appreciable effects of soils, streams, and reservoirs on TP transport, The estimated aquatic rates of phosphorus removal declined with increasing stream size and rates of water flushing in reservoirs (i.e. areal hydraulic loads). A phosphorus budget for the 2.9 million km2 Mississippi River Basin provides a detailed accounting of TP delivery to streams, the removal of TP in surface waters, and the stream export of TP from major interior watersheds for sources associated with each land-use type.

Water Science and Technology

Hydrology and trout populations of cold-water rivers of Michigan and Wisconsin

Statistical multiple-regression analyses showed significant relationships between trout populations and hydrologic parameters. Parameters showing the higher levels of significance were temperature, hardness of water, percentage of gravel bottom, percentage of bottom vegetation, variability of streamflow, and discharge per unit drainage area. Trout populations increase with lower levels of annual maximum water temperatures, with increase in water hardness, and with increase in percentage of gravel and bottom vegetation. Trout populations also increase with decrease in variability of streamflow, and with increase in discharge per unit drainage area. Most hydrologic parameters were significant when evaluated collectively, but no parameter, by itself, showed a high degree of correlation with trout populations in regression analyses that included all the streams sampled. Regression analyses of stream segments that were restricted to certain limits of hardness, temperature, or percentage of gravel bottom showed improvements in correlation. Analyses of trout populations, in pounds per acre and pounds per mile and hydrologic parameters resulted in regression equations from which trout populations could be estimated with standard errors of 89 and 84 per cent, respectively.

Michigan, Wisconsin

Statistical analysis of the radon-222 potential of rocks in Virginia, U.S.A.

More than 3,200 indoor radon-222 (222Rn) measurements were made seasonally in an area of about 1,000 square kilometers of the Coastal Plain and Piedmont physiographic provinces in Virginia, U.S.A. Results of these measurements indicate that some geological units are associated, on the average, with twice as much indoor222Rn as other geological units, and that indoor222Rn varies seasonally. The Kruskal-Wallis test was used to test whether indoor222Rn concentrations for data gathered over the winter and summer seasons differ significantly by rock unit. The tests concluded that indoor222Rn concentrations for different rock units were not equal at the 5-percent significance level. The rocks associated with the highest median indoor222Rn concentration are specific rocks in the Mesozoic Culpeper basin, including shale and siltstone units with Jurassic diabase intrusives, and mica schists in the Piedmont physiographic province. The pre-Triassic Peters Creek Schist has the highest ranking in terms of indoor222Rn concentration. The rocks associated with the lowest indoor222Rn concentrations include coastal plain sediments, the Occoquan Granite, Falls Church Tonalite, Piney Branch Mafic and Ultramafic complex, and unnamed mafic and ultramafic inclusions, respectively. The rocks have been ranked according to observed222Rn concentration by transforming the average rank of indoor222Rn concentrations to z scores. ?? 1992 Springer-Verlag New York Inc.

Environmental Geology and Water Sciences

Disentangling the effects of multiple fires on spatially interspersed sagebrush (Artemisia spp.) communities

Questions Relative to a landscape with a mosaic of two sagebrush community types and increasing fire frequency, we asked: (a) do vegetation characteristics vary significantly with number of times burned for each sagebrush community; (b) how do vegetation responses to different fire frequencies compare between the two sagebrush communities? Location Columbia Plateau Ecoregion, Washington, USA. Methods We sampled vegetation across a landscape that burned three times over a 10‐year period in two sagebrush community types that are interspersed on unique land forms: big sagebrush ( Artemisia tridentata ) communities that occur on small “mounds” and scabland sagebrush ( A. rigida ) communities that occur on surrounding “flats.” Spatially overlapping fires permitted a balanced sampling design to assess unburned and once‐, twice‐, and thrice‐burned locations for each land form/community type. We utilized a suite of statistical analyses to determine differences among plant functional groups and biomass among unburned/burned strata by land form and compared results between land forms. Results Big sagebrush and scabland sagebrush communities responded uniquely to multiple fires, due to different fuel loadings, fire severities, succession and invasion dynamics. Big sagebrush experienced nearly complete shrub loss and conversion from exotic‐invaded shrubland to exotic annual grassland after only one fire. In contrast, scabland sagebrush retained a minor shrub component and higher relative cover of native herbaceous species, even after three fires. Both communities retained cover of native perennial grasses, including shallow‐ and deep‐rooted species, likely reflecting decreasing fire intensity with number of times burned. Conclusions Despite different community‐level responses, increasing fire frequency is transforming the entire landscape to a non‐native/native grassland mix. Quantifying unique ecosystem responses to altered wildfire regimes is critical to understanding the relative resilience of communities to disturbance and their resistance to exotic species invasion (and community type conversion). Management actions may help to maintain spatial heterogeneity of ecosystems and fire‐tolerant native species.

Washington

A novel framework to predict relative habitat selection in aquatic systems: Applying machine learning and resource selection functions to acoustic telemetry data from multiple shark species

Resource selection functions (RSFs) have been widely applied to animal tracking data to examine relative habitat selection and to help guide management and conservation strategies. While readily used in terrestrial ecology, RSFs have yet to be extensively used within marine systems. As acoustic telemetry continues to be a pervasive approach within marine environments, incorporation of RSFs can provide new insights to help prioritize habitat protection and restoration to meet conservation goals. To overcome statistical hurdles and achieve high prediction accuracy, machine learning algorithms could be paired with RSFs to predict relative habitat selection for a species within and even outside the monitoring range of acoustic receiver arrays, making this a valuable tool for marine ecologists and resource managers. Here, we apply RSFs using machine learning to an acoustic telemetry dataset of four shark species to explore and predict species-specific habitat selection within a marine protected area. In addition, we also apply this RSF-machine learning approach to investigate predator-prey relationships by comparing and averaging tiger shark relative selection values with the relative selection values derived for eight potential prey-species. We provide methodological considerations along with a framework and flexible approach to apply RSFs with machine learning algorithms to acoustic telemetry data and suggest marine ecologists and resource managers consider adopting such tools to help guide both conservation and management strategies.

Buck Island Reef National Monument, St Croix

Spatial and temporal statistical analysis of bycatch data: Patterns of sea turtle bycatch in the North Atlantic

Loggerhead (Caretta caretta) and leatherback (Dermochelys coriacea) sea turtle distributions and movements in offshore waters of the western North Atlantic are not well understood despite continued efforts to monitor, survey, and observe them. Loggerhead and leatherback sea turtles are listed as endangered by the World Conservation Union, and thus anthropogenic mortality of these species, including fishing, is of elevated interest. This study quantifies spatial and temporal patterns of sea turtle bycatch distributions to identify potential processes influencing their locations. A Ripley's K function analysis was employed on the NOAA Fisheries Atlantic Pelagic Longline Observer Program data to determine spatial, temporal, and spatio-temporal patterns of sea turtle bycatch distributions within the pattern of the pelagic fishery distribution. Results indicate that loggerhead and leatherback sea turtle catch distributions change seasonally, with patterns of spatial clustering appearing from July through October. The results from the space-time analysis indicate that sea turtle catch distributions are related on a relatively fine scale (30-200 km and 1-5 days). The use of spatial and temporal point pattern analysis, particularly K function analysis, is a novel way to examine bycatch data and can be used to inform fishing practices such that fishing could still occur while minimizing sea turtle bycatch. ?? 2008 NRC.

Canadian Journal of Fisheries and Aquatic Sciences

Modelling gully-erosion susceptibility in a semi-arid region, Iran: Investigation of applicability of certainty factor and maximum entropy models

Gully erosion susceptibility mapping is a fundamental tool for land-use planning aimed at mitigating land degradation. However, the capabilities of some state-of-the-art data-mining models for developing accurate maps of gully erosion susceptibility have not yet been fully investigated. This study assessed and compared the performance of two different types of data-mining models for accurately mapping gully erosion susceptibility at a regional scale in Chavar, Ilam, Iran. The two methods evaluated were: Certainty Factor (CF), a bivariate statistical model; and Maximum Entropy (ME), an advanced machine learning model. Several geographic and environmental factors that can contribute to gully erosion were considered as predictor variables of gully erosion susceptibility. Based on an existing differential GPS survey inventory of gully erosion, a total of 63 eroded gullies were spatially randomly split in a 70:30 ratio for use in model calibration and validation, respectively. Accuracy assessments completed with the receiver operating characteristic curve method showed that the ME-based regional gully susceptibility map has an area under the curve (AUC) value of 88.6% whereas the CF-based map has an AUC of 81.8%. According to jackknife tests that were used to investigate the relative importance of predictor variables, aspect, distance to river, lithology and land use are the most influential factors for the spatial distribution of gully erosion susceptibility in this region of Iran. The gully erosion susceptibility maps produced in this study could be useful tools for land managers and engineers tasked with road development, urbanization and other future development.

Science of the Total Environment

Specific Conductance and Dissolved-Solids Characteristics for the Green River and Muddy Creek, Wyoming, Water Years 1999-2008

Southwestern Wyoming is an area of diverse scenery, wildlife, and natural resources that is actively undergoing energy development. The U.S. Department of the Interior's Wyoming Landscape Conservation Initiative is a long-term science-based effort to assess and enhance aquatic and terrestrial habitats at a landscape scale, while facilitating responsible energy development through local collaboration and partnerships. Water-quality monitoring has been conducted by the U.S. Geological Survey on the Green River near Green River, Wyoming, and Muddy Creek near Baggs, Wyoming. This monitoring, which is being conducted in cooperation with State and other Federal agencies and as part of the Wyoming Landscape Conservation Initiative, is in response to concerns about potentially increased dissolved solids in the Colorado River Basin as a result of energy development. Because of the need to provide real-time dissolved-solids concentrations for the Green River and Muddy Creek on the World Wide Web, the U.S. Geological Survey developed regression equations to estimate dissolved-solids concentrations on the basis of continuous specific conductance using relations between measured specific conductance and dissolved-solids concentrations. Specific conductance and dissolved-solids concentrations were less varied and generally lower for the Green River than for Muddy Creek. The median dissolved-solids concentration for the site on the Green River was 318 milligrams per liter, and the median concentration for the site on Muddy Creek was 943 milligrams per liter. Dissolved-solids concentrations ranged from 187 to 594 milligrams per liter in samples collected from the Green River during water years 1999-2008. Dissolved-solids concentrations ranged from 293 to 2,485 milligrams per liter in samples collected from Muddy Creek during water years 2006-08. The differences in dissolved-solids concentrations in samples collected from the Green River compared to samples collected from Muddy Creek reflect the different basin characteristics. Relations between specific conductance and dissolved-solids concentrations were statistically significant for the Green River (p-value less than 0.001) and Muddy Creek (p-value less than 0.001); therefore, specific conductance can be used to estimate dissolved-solids concentrations. Using continuous specific conductance values to estimate dissolved solids in real-time on the World Wide Web increases the amount and improves the timeliness of data available to water managers for assessing dissolved-solids concentrations in the Colorado River Basin.

Scientific Investigations Report

Predictive analysis using chemical-gene interaction networks consistent with observed endocrine activity and mutagenicity of U.S. streams

In a recent U.S. Geological Survey/U.S. Environmental Protection Agency study assessing >700 organic compounds in 38 streams, in vitro assays indicated generally low estrogen, androgen, and glucocorticoid receptor activities, but identified 13 surface waters with 17β estradiol equivalent (E2Eq) activities greater than the 1 ng/L level of concern for feminization of male fish. Among the 36 samples assayed for mutagenicity in the Salmonella bioassay (reported here), 25% were considered mutagenic (statistically significant slope and at least a two-fold increase in revertants/plate). Endocrine and mutagenic activities of the water samples were well correlated with each other and with the total number and cumulative concentrations of detected chemical contaminants. To test the predictive utility of knowledgebase-leveraging approaches, site-specific predicted chemical-gene (pCGA) and predicted analogous pathway-linked (pPLA) association networks identified in the Comparative Toxicogenomics Database were compared with observed endocrine/mutagenic bioactivities. We evaluated pCGA/pPLA patterns among sites by cluster analysis and principal component analysis and grouped the pPLA into broad mode-of-action classes. Measured E2Eq and mutagenic activities correlated well with predicted pathways. The pPLA analysis also revealed correlations with signaling, metabolic, and regulatory groups, suggesting that other effects pathways may be associated with chemical contaminants in these waters and indicating the need for broader bioassay coverage to assess potential adverse impacts.

Environmental Science & Technology

Atmospheric nitrogen in the Mississippi River Basin: Amissions, deposition and transport

Atmospheric deposition of nitrogen has been cited as a major factor in the nitrogen saturation of forests in the north-eastern United States and as a contributor to the eutrophication of coastal waters, including the Gulf of Mexico near the mouth of the Mississippi River. Sources of nitrogen emissions and the resulting spatial patterns of nitrogen deposition within the Mississippi River Basin, however, have not been fully documented. An assessment of atmospheric nitrogen in the Mississippi River Basin was therefore conducted in 1998-1999 to: (1) evaluate the forms in which nitrogen is deposited from the atmosphere; (2) quantify the spatial distribution of atmospheric nitrogen deposition throughout the basin; and (3) relate locations of emission sources to spatial deposition patterns to evaluate atmospheric transport. Deposition data collected through the NADP/NTN (National Atmospheric Deposition Program/National Trends Network) and CASTNet (Clean Air Status and Trends Network) were used for this analysis. NO(x) Tier 1 emission data by county was obtained for 1992 from the US Environmental Protection Agency (Emissions Trends Viewer CD, 1985-1995, version 1.0, September 1996) and NH3 emissions data was derived from the 1992 Census of Agriculture (US Department of Commerce. Census of Agriculture, US Summary and County Level Data, US Department of Commerce, Bureau of the Census. Geographic Area series, 1995:1b) or the National Agricultural Statistics Service (US Department of Agriculture. National Agricultural Statistics Service Historical Data. Accessed 7/98 at URL, 1998. http://www.usda.gov/nass/pubs/hisdata.htm). The highest rates of wet deposition of NO3- were in the north-eastern part of the basin, downwind of electric utility plants and urban areas, whereas the highest rates of wet deposition of NH4+ were in Iowa, near the center of intensive agricultural activities in the Midwest. The lowest rates of atmospheric nitrogen deposition were on the western (windward) side of the basin, which suggests that most of the nitrogen deposited within the basin is derived from internal sources. Atmospheric transport eastward across the basin boundary is greater for NO3- than NH4+, but a significant amount of NH4+ is likely to be transported out of the basin through the formation of (NH4)2SO4 and NH4NO3 particles - a process that greatly increases the atmospheric residence time of NH4+. This process is also a likely factor in the atmospheric transport of nitrogen from the Midwest to upland forest regions in the North-East, such as the western Adirondack region of New York, where NH4+ constitutes 38% of the total wet deposition of N.

Science of Total Environment

Hazard potential of compound flooding from rainfall, storm surge, and groundwater in coastal New York and Connecticut

Compound flood events, the co-occurrence of multiple flood drivers, can result in flood hazard potential exceeding that of any single driver alone. To evaluate compound flooding in a semi-urbanized coastal area, historical records dating back to 1970 are used to study the co-occurrences of high precipitation, storm surge, and shallow groundwater conditions along the coastlines of New York and Connecticut. Joint return periods for coincident precipitation-surge events were computed using statistical dependence models and compared to the assumption of independence as a ratio, referred to here as a return period adjustment. Results indicate distinct seasonality where compound events in the area disproportionately occur in the cold season between October and April. Return period adjustments range from a factor of 1 to almost 9, demonstrating the range in precipitation-storm surge dependence across the study area. Across all 24 station triad locations, groundwater levels were elevated during times of precipitation- surge co-occurrence, reflecting the tendency for coastal storms and shallow groundwater conditions to co-occur seasonally. The result is a pseudo-trivariate compound flood hazard score and corresponding hazard map that integrates dependence between daily precipitation-surge events and overall monthly groundwater levels (as a precondition) into a relative compound hazard score. The location with the highest compound flood hazard score is on the south shore of Long Island, as well as locations across coastal Connecticut where groundwater levels compound the co-occurrence of heavy precipitation and storm surge.

Connecticut, New York

Characterization of salinity and selenium loading and land-use change in Montrose Arroyo, western Colorado, from 1992 to 2010

Salinity and selenium are naturally occurring and perva-sive in the lower Gunnison River Basin of Colorado, includ-ing the watershed of Montrose Arroyo. Although some of the salinity and selenium loading in the Montrose Arroyo study area is from natural sources, additional loading has resulted from the introduction of intensive irrigation in the water-shed. With increasing land-use change and the conversion from irrigated agricultural to urban land, land managers and stakeholders need information about the long-term effects of land-use change on salinity and selenium loading. In response to the need to advance salinity and selenium science, the U.S. Geological Survey, in cooperation with the Bureau of Reclamation, Colorado River Basin Salinity Control Forum, and Colorado River Water Conservation District, developed a study to characterize salinity and selenium loading and how salinity and selenium sources may relate to land-use change in Montrose Arroyo. This report characterizes changes in salinity and selenium loading to Montrose Arroyo from March 1992 to February 2010 and the magnitude of land-use change between unirrigated desert, irrigated agricultural, and urban land-use/land-cover types, and discusses how the respective loads may relate to land-use change. Montrose Arroyo is an approximately 8-square-mile watershed in Montrose County in western Colorado. Salinity and selenium were studied in Montrose Arroyo in a 2001 study as part of a salinity- and selenium-control lateral project. The robust nature of the historical dataset indicated that Montrose Arroyo was a prime watershed for a follow-up study. Two sites from the 2001 study were used to monitor salinity and selenium loads in Montrose Arroyo in the follow-up study. Over the period of 2 water years and respective irrigation seasons (2008-2010), 27 water-quality samples were collected and streamflow measurements were made at the historical sites MA2 and MA4. Salinity and selenium concen-trations, loads, and streamflow were compared between the pre-lateral-project and post-growth periods and between the post-lateral-project and post-growth periods. No significant differences in streamflow, salinity (concen-tration and load), or selenium (concentration and load) were found at MA4 between the pre-lateral project and post-growth periods or between the post-lateral-project and post-growth periods. The statistical analysis indicated no significant dif-ferences in streamflow or salinity (both concentration and load) between the pre-lateral-project and post-growth periods or between the post-lateral-project and post-growth periods at MA2; however, selenium concentrations and loads were significantly greater between the pre-lateral-project and post-growth periods and between the post-lateral-project and post-growth periods at MA2. Land-use change between MA4 and MA2 may have contributed to the determined differences in selenium values, but the specific mechanisms causing the increases between periods are unknown. The size of the urbanized area in Montrose Arroyo was quantified for 1993, 2002, and 2009 by using a geographic information system (GIS) with imagery from the specified years. The greatest change in land use from 1993 to 2009 was the increase of urban land due to conversion from irrigated agricultural land. The conversion of previously unirrigated desert to urban land or irrigated agriculture could become more common if urbanization and development expands into the eastern part of the watershed because a majority of the un-urbanized land in eastern Montrose Arroyo is unirrigated desert. By applying GIS to the City of Montrose 2008 com-prehensive growth plan, it was estimated that approximately 786 acres of previously irrigated agricultural land will be converted to urban land and 689 acres of unirrigated desert will be converted to urban land under the plan scenario. New development on previously unirrigated land in shale areas would likely increase the potential for mobilization of sele-nium and salinity from new sources to Montrose Arroyo and the Lower Gunnison River Basin.

Scientific Investigations Report

Latency and geofence testing of wireless emergency alerts intended for the ShakeAlert® earthquake early warning system for the West Coast of the United States of America

ShakeAlert, the earthquake early warning (EEW) system for the West Coast of the United States, attempts to provides crucial warnings before strong shaking occurs. However, because the alerts are triggered only when an earthquake is already in progress, and the alert latencies and delivery times are platform dependent, the time between these warnings and the arrival of shaking is variable. The ShakeAlert system uses, among other public alerting platforms like a mobile phone operating system, smartphone apps, and the Federal Emergency Management Agency Integrated Public Alert & Warning System (IPAWS). IPAWS sends Wireless Emergency Alerts (WEAs) informing people via their smartphones and other mobile devices about various events, such as natural hazards, child abductions, or public health information about COVID-19. However, little is known about the IPAWS delivery latencies. Given that people may have only a few seconds of notice after they receive an alert to take a protective action before they feel earthquake shaking, quantifying latencies is critical to understanding whether the IPAWS system is useful for EEW. In this study, we developed new methods to test the IPAWS distribution system's performance, both with devices in a controlled environment and as well as with a 2019 community-based feedback form, in Oakland and San Diego County, California, respectively. The controlled environment test used mobile phones (including smart and non-smart phones) and associated devices to determine alert receipt times; the community research form had participants self-report their receipt times. By triangulating the data between the controlled test environment and the community research, we determined the latency statistics as well as whether the geofence (the geographic area where the alert was intended to be sent) held broadly. We found that the latencies were similar between the two tests despite the large differences in population sizes. WEA messages were received within a median time frame of 6–12 s, and the geofence held with only a few exceptions. We use this latency to assess how the system would have performed in two large earthquakes, the 1989 M6.9 Loma Prieta and 2019 M7.1 Ridgecrest earthquakes, which both occurred near our WEA test locations. Our analysis revealed that had IPAWS been available during those earthquakes, particularly Loma Prieta, it would have provided crucial seconds of notice that damaging shaking was imminent in some locations relatively far from the epicenter. Further, we find affordable non-smart phones can receive WEAs as fast as smartphones. Finally, our new method can be used for latency and geospatial testing going forward for IPAWS and other similar alerting systems.

California

Landslide initiation thresholds in data-sparse regions: Application to landslide early warning criteria in Sitka, Alaska, USA

Probabilistic models to inform landslide early warning systems often rely on rainfall totals observed during past events with landslides. However, these models are generally developed for broad regions using large catalogs, with dozens, hundreds, or even thousands of landslide occurrences. This study evaluates strategies for training landslide forecasting models with a scanty record of landslide-triggering events, which is a typical limitation in remote, sparsely populated regions. We evaluate 136 statistical models trained on a precipitation dataset with five landslide-triggering precipitation events recorded near Sitka, Alaska, USA, as well as > 6000 d of non-triggering rainfall (2002–2020). We also conduct extensive statistical evaluation for three primary purposes: (1) to select the best-fitting models, (2) to evaluate performance of the preferred models, and (3) to select and evaluate warning thresholds. We use Akaike, Bayesian, and leave-one-out information criteria to compare the 136 models, which are trained on different cumulative precipitation variables at time intervals ranging from 1 h to 2 weeks, using both frequentist and Bayesian methods to estimate the daily probability and intensity of potential landslide occurrence (logistic regression and Poisson regression). We evaluate the best-fit models using leave-one-out validation as well as by testing a subset of the data. Despite this sparse landslide inventory, we find that probabilistic models can effectively distinguish days with landslides from days without slide activity. Our statistical analyses show that 3 h precipitation totals are the best predictor of elevated landslide hazard, and adding antecedent precipitation (days to weeks) did not improve model performance. This relatively short timescale of precipitation combined with the limited role of antecedent conditions likely reflects the rapid draining of porous colluvial soils on the very steep hillslopes around Sitka. Although frequentist and Bayesian inferences produce similar estimates of landslide hazard, they do have different implications for use and interpretation: frequentist models are familiar and easy to implement, but Bayesian models capture the rare-events problem more explicitly and allow for better understanding of parameter uncertainty given the available data. We use the resulting estimates of daily landslide probability to establish two decision boundaries that define three levels of warning. With these decision boundaries, the frequentist logistic regression model incorporates National Weather Service quantitative precipitation forecasts into a real-time landslide early warning “dashboard” system ( https://sitkalandslide.org/ , last access: 9 October 2023). This dashboard provides accessible and data-driven situational awareness for community members and emergency managers.

Alaska

Lithogeochemistry of Carlin-type gold mineralization in the Gold Bar district, Battle Mountain-Eureka trend, Nevada

The Gold Bar district contains five Carlin-type gold deposits and four resources for a combined gold endowment of 1.6 M oz [50 t]. The gold deposits are hosted in Devonian carbonate rocks below parautochthonous and allochthonous Paleozoic siliciclastic rocks emplaced during the Early Mississippian Antler orogeny. The district is in the Battle Mountain-Eureka trend, a long-lived structural feature that localized intrusions and ore deposits of different types and ages. The whole-rock geochemistry of four different mineralized and unmineralized Devonian carbonate rock units (two favorable and two unfavorable) were determined and interpreted in the context of the regional geology. A combination of basic statistics, R-mode factor analysis, isocon plots, and alteration diagrams were utilized to (1) identify favorable geochemical attributes of the host rocks, (2) characterize alteration and associated element enrichments and depletions, and (3) identify the mechanism of gold precipitation. This approach also led to the recognition of other types of alteration and mineralization in host rocks previously thought to be solely affected by Carlin-type mineralization. Unit 2 of the Upper Member of the Denay Formation, with the highest Al2O3, Fe2O3 and SiO2 contents and the lowest CaO content, is the most favorable host rock. Based on the high regression coefficients of data arrays on X-Y plots that project toward the origin, Al2O3 and TiO2 were immobile and K2O and Fe2O3 were relatively immobile during alteration and mineralization. Specific element associations identified by factor analysis are also prominent on isocon diagrams that compare the composition of fresh and altered equivalents of the same rock units. The most prominent associations are: Au, As, Sb, SiO2, TI, -CaO and -LOI, the main gold mineralizing event and related silicification and decalcification; Cd, Zn, Ag, P, Ni and Tl, an early base metal event; and MgO, early dolomitization. Alteration diagrams, consisting of X-Y plots of SiO2/Al2O3, K2O/Al2O3, CO2/Al2O3, that documented in classic Carlin-type gold deposits in the region, but the size of the deposits and the intensity of alteration and mineralization are less. The presence of other types of mineralization in the Gold Bar district is also common to most of the other Carlin-type districts located in major mineral belts. The approach used in this study is well suited to the interpretation of multi-element geochemical data from other study areas with superimposed alteration and mineralization. ?? 2002 Elsevier Science B.V. All rights reserved.

Ore Geology Reviews

Factors contributing to pesticide contamination in riverine systems: The role of wastewater and landscape sources

Wastewater treatment plant (WWTP) discharges can be a source of organic contaminants, including pesticides, to rivers. An integrated model was developed for the Potomac River watershed (PRW) to determine the amount of accumulated wastewater percentage of streamflow (ACCWW) and calculate predicted environmental concentrations (PECs) for 14 pesticides in non-tidal National Hydrography Dataset Plus Version 2.1 stream segments. Predicted environmental concentrations were compared to measured environmental concentrations (MECs) from 32 stream sites that represented a range of ACCWW and land use to evaluate model performance and to assess possible non-WWTP loading sources. Statistical agreement between PECs and MECs was strongest for insecticides, followed by fungicides and herbicides. Principal component analysis utilizing optical fluorescence and ancillary water quality data identified wastewater and urban runoff sources. Pesticides that indicated relatively larger sources from WWTPs included dinotefuran, fipronil, carbendazim, thiabendazole, and prometon whereas imidacloprid, azoxystrobin, propiconazole, tebuconazole, and diuron were more related to urban runoff. In addition, PECs generally comprised a low proportion of MECs, which indicates possible dominant loading sources beyond WWTP discharges. Cumulative potential toxicity was higher for sites with greater ACCWW and/or located in developed areas. Imidacloprid, fipronil, and carbendazim accounted for the largest portion of predicted potential toxicity across sites. The chronic aquatic life toxicity benchmarks for freshwater invertebrates were exceeded for 82 % of the imidacloprid detections ( n = 28) and 47 % of the fipronil detections ( n = 19). These results highlight the ecological implications of pesticide contamination from WWTP discharges and also the potential legacy effects from accumulated soil and groundwater sources. Pesticide management strategies that mitigate both current and historical impacts may improve the health of aquatic ecosystems.

Potomac River watershed