SEARCH · USGS Science

Results for “Nature Based Solutions”

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

SourceResearch connectionFrom yearThrough year

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

At least 181 records · Page 10Linked to original sources

USGS2001

Interactions between fire and invasive plants in temperate grasslands of North America

A substantial number of invasive grasses, forbs and woody plants have invaded temperate grasslands in North America. Among the invading species are winter annuals, biennials, cool-season perennials, warm-season perennials, vines, shrubs, and trees. Many of these species have been deliberately introduced and widely planted; some are still used for range improvement, pastures, lawns, and as ornamentals, though many are listed as state or federal noxious weeds. Others have been greatly facilitated by widespread land disturbance. Historically, fire has been a major selective force in the evolution of temperate grasslands. Further, prescribed fire is commonly used as a method of ecological management for native grassland communities as well as in conjunction with restoration efforts. Within this context, it is important to understand how invasive species will interact with natural and prescribed fire regimes. In this paper, we consider what is known about how exotic species that invade temperate grasslands relate to fire. The primary issues addressed for each species are (1) Does fire appear to enhance colonization? (2) To what degree does fire affect the survival of plants? (3) Are plants that are burned able to regrow following fire and, if so, how rapidly can they recover? (4) How important is competition with native species to the response to fire? and (5) What effect does an invasive species have on the characteristics of the fire regime? For many species, results are preliminary, incomplete, or inconsistent among studies. For this reason, many of the conclusions drawn for individual species must be considered preliminary. Based on analyses of individual species, a conceptual framework is presented for considering how invasive plants may interact with fire when they invade an ecosystem. The major categories of influences are the native community, the fire regime, growth conditions for both invasive and native species, and influences that disturbances, human impacts, and landscape characteristics have had in the past and will have in the future. The examples considered in this paper provide support for a few, tentative generalizations. First, among our current worst invaders of temperate grasslands, adaptation to fire is quite variable. Some species are not well adapted to burning and can be easily eliminated; other species are better adapted but can still be eliminated if fire occurs during periods of particular vulnerability and/or at high frequency. There is a set of species that is extremely well adapted to fire and will not be eliminated through burning alone. Second, competitive interactions with native species play a crucial role in the success of nonnative invaders. In cases where differential burn responses between invasive and native species can be exploited, and adequate populations of native dominant species are present, fire can sometimes tip the competitive balance away from invasives. Third, there are a few invasive species that have exceptional attributes and for which there are no easy solutions. The ability of cheatgrass ( Bromus tectorum ) to enhance fire, the ability of Chinese tallow ( Triadica sebifera ) to suppress fire, and the ability of leafy spurge ( Euphorbia esula ) to resprout from repeated injury make these species exceptional threats to native diversity. Finally, the available information for many invasive species is very incomplete, partic- ularly with regard to how fire affects competitive interactions with the native community. There is much more we need to know if we are to consistently predict how invasive species will respond to fire and how burning can best be used to manage for natural diversity.

Conference Paper↗

USGS1994

Horvitz-Thompson survey sample methods for estimating large-scale animal abundance

Large-scale surveys to estimate animal abundance can be useful for monitoring population status and trends, for measuring responses to management or environmental alterations, and for testing ecological hypotheses about abundance. However, large-scale surveys may be expensive and logistically complex. To ensure resources are not wasted on unattainable targets, the goals and uses of each survey should be specified carefully and alternative methods for addressing these objectives always should be considered. During survey design, the impoflance of each survey error component (spatial design, propofiion of detected animals, precision in detection) should be considered carefully to produce a complete statistically based survey. Failure to address these three survey components may produce population estimates that are inaccurate (biased low), have unrealistic precision (too precise) and do not satisfactorily meet the survey objectives. Optimum survey design requires trade-offs in these sources of error relative to the costs of sampling plots and detecting animals on plots, considerations that are specific to the spatial logistics and survey methods. The Horvitz-Thompson estimators provide a comprehensive framework for considering all three survey components during the design and analysis of large-scale wildlife surveys. Problems of spatial and temporal (especially survey to survey) heterogeneity in detection probabilities have received little consideration, but failure to account for heterogeneity produces biased population estimates. The goal of producing unbiased population estimates is in conflict with the increased variation from heterogeneous detection in the population estimate. One solution to this conflict is to use an MSE-based approach to achieve a balance between bias reduction and increased variation. Further research is needed to develop methods that address spatial heterogeneity in detection, evaluate the effects of temporal heterogeneity on survey objectives and optimize decisions related to survey bias and variance. Finally, managers and researchers involved in the survey design process must realize that obtaining the best survey results requires an interactive and recursive process of survey design, execution, analysis and redesign. Survey refinements will be possible as further knowledge is gained on the actual abundance and distribution of the population and on the most efficient techniques for detection animals.

Transactions of the North American Wildlife and Na↗

USGS2006

Spatial analysis of land use and shallow groundwater vulnerability in the watershed adjacent to Assateague Island National Seashore, Maryland and Virginia, USA

Spatial relations between land use and groundwater quality in the watershed adjacent to Assateague Island National Seashore, Maryland and Virginia, USA were analyzed by the use of two spatial models. One model used a logit analysis and the other was based on geostatistics. The models were developed and compared on the basis of existing concentrations of nitrate as nitrogen in samples from 529 domestic wells. The models were applied to produce spatial probability maps that show areas in the watershed where concentrations of nitrate in groundwater are likely to exceed a predetermined management threshold value. Maps of the watershed generated by logistic regression and probability kriging analysis showing where the probability of nitrate concentrations would exceed 3 mg/L (>0.50) compared favorably. Logistic regression was less dependent on the spatial distribution of sampled wells, and identified an additional high probability area within the watershed that was missed by probability kriging. The spatial probability maps could be used to determine the natural or anthropogenic factors that best explain the occurrence and distribution of elevated concentrations of nitrate (or other constituents) in shallow groundwater. This information can be used by local land-use planners, ecologists, and managers to protect water supplies and identify land-use planning solutions and monitoring programs in vulnerable areas. ?? 2006 Springer-Verlag.

Environmental Geology↗

USGS1994

Variations in the chemical and stable isotope composition of carbon and sulfur species during organic-rich sediment alteration: An experimental and theoretical study of hydrothermal activity at Guaymas Basin, Gulf of California

Organic-rich diatomaceous ooze was reacted with seawater and a Na-Ca-K-Cl fluid of seawater chlorinity at 325–400°C, 400–500 bars, and fluid/sediment mass ratios of 1.56–2.35 to constrain factors regulating the abundance and stable isotope composition of C and S species during hydrothermal alteration of sediment from Guaymas Basin, Gulf of California. Alteration of inorganic and organic sedimentary components resulted in extensive exchange reactions, the release of abundant H 2 S, CO 2 , CH 4 , and C organic , to solution, and recrystallization of the sediment to an assemblage containing albitic plagioclase, quartz, pyrrhotite, and calcite. The δ 34 S cdt "> δ 34 S cdt values of dissolved H 2 S varied from −10.9 to +4.3‰ during seawater-sediment interaction at 325 and 400°C and from −16.5 to −9.0‰ during Na-Ca-K-Cl fluid-sediment interaction at 325 and 375°C. In the absence of seawater SO 4 , H 2 S is derived from both the transformation of pyrite to pyrrhotite and S released during the degradation of organic matter. In the presence of seawater SO 4 , reduction of SO 4 contributes directly to H 2 S production. Sedimentary organic matter acts as the reducing agent during pyrite and SO 4 reduction. Requisite acidity for the reduction of SO 4 is provided by Mg fixation during early-stage sediment alteration and by albite and calcite formation in Mg-free solutions. Organically derived CH 4 was characterized by δ 13 C pdb "> δ13Cpdb values ranging between −20.8 and −23.1‰, whereas δ 13 C pdb "> δ 13 C pdb values for dissolved C organic ranged between −14.8 and −17.7%. Mass balance calculations indicate that δ 13 C values for organically derived CO 2 were ≥ − 14.8%. Residual solid sedimentary organic C showed small (≤ 0.7‰) depletions in 13 C relative to the starting sediment. The experimental results are consistent with the isotopic and chemical composition of natural hydrothermal fluids and minerals at Guaymas Basin and permit us to better constrain sources and sinks for C and S species in subseafloor hydrothermal systems at sediment-covered spreading centers. Our data show that the sulfur isotope composition of hydrothermal Sulfide minerals in Guaymas Basin can be explained by derivation of S from diagenetic sulfide and seawater sulfate. Basaltic S may also contribute to hydrothermal sulfide precipitates but is not required to explain their isotopic composition. Estimates of seawater/ sediment mass ratios based on sulfur isotopic composition of sulfide minerals and the abundance of dissolved NH 3 in vent fluids range from 3–29 during hydrothermal circulation. Sources of C in Guaymas Basin hydrothermal fluids include thermal degradation of organic matter, bacteriogenic methane production, and dissolution of diagenetic carbonate.

Baja California Sur↗

USGS1952

Geology of the Iron King Mine, Yavapai county, Arizona

The Iron King mine is about 2,000 feet west-northwest of the intersection of the 112 15 west meridian and the 34 30 north parallel in the Humboldt region in central Yavapai County , Arizona . The mine is approximately in the geographical center of the Humboldt region. Precambrian rocks form the bedrock. Late Cenozoic unconsolidated river wash and valley fill with some interbedded basalt locally mantle the Precambrian rocks, especially in the north-central part of the region. The Precambrian rocks consist of two metamorphosed volcanic formations and intrusive rocks that range in composition from quartz porphyry to gabbro. The volcanic formations originally were flows, volcanic breccias, and tuffaceous sedimentary rocks. Dynamothermal metamorphism of these rocks formed textures, structures, and mineral assemblages characteristic of low-grade ;riietamorphic rocks, but sufficient relict textures and structures remain to permit delineation of formations on the bases of their original nature. All the Precambrian rocks are foliated, except those in the interior parts of the larger intrusive masses. This foliation has two major trends, (1) north to N 20° W, and (2) N 20° E. The northeast-trending foliation is younger and locally is superimposed on the north- to northwest-trending foliation. The Precambrian rocks strike north to northwest and dip steeply, chiefly westward. Duplication of stratigraphic units and determination of tops suggest two major northeast-trending folds, which probably plunge southward. The northeast-trending foliation appears to bear an axial-plane relationship to these folds. Some masses of igneous rock probably were intruded during development of foliation. The Iron King deposit supports the only active mine in the Humboldt region. It consists of 12 steeply plunging echelon veins arranged along the footwall of a sheared and altered zone in the metamorphosed andesitic tuffaceous sedimentary rocks. Narrow zones of more intense shear probably localized the veins. Solutions first introduced quartz, pyrite, ankerite, and sericite, forming a sporadically mineralized zone (in the hanging wall of the deposit) and probably veins in the Iron King fracture system. Intra-mineralization shear strongly brecciated these early minerals and formed the structures that localized and distributed the ore minerals in the veins. After this deformation, sphalerite, galena, chalcopyrite, tennantite, arsenopyrite, pyrite, quartz, and ankerite were deposited; the last three minerals may have formed through solution and the redeposition of earlier minerals. Formation of sericite either accompanied or followed the deposition of ore-forming minerals. Silver is related closely in distribution to copper and probably is in the tennantite. Gold occurs chiefly in the pyrite. Banding, mimetic after foliation and planes developed by shearing, is pronounced in much of the vein material. Mineral zoning, generally similar in each vein, is a characteristic of the deposit. High-angle reverse faults of about 100 feet in maximum vertical separation offset the veins. These faults are nearly parallel to the veins in strike and dip.

Arizona↗

USGS2026

Determining Volcanic Risk in Auckland (DEVORA) Research Programme—A transdisciplinary approach to address the challenge of distributed volcanism in an urban environment

The Determining Volcanic Risk in Auckland (DEVORA) Research Programme was launched in 2008 to address the challenges associated with monogenetic volcanism in an urban setting and to enhance volcanic risk management in Tāmaki Makaurau Auckland in Aotearoa New Zealand. It is a multi-agency, increasingly transdisciplinary (defined here as research that transcends traditional disciplinary boundaries by integrating diverse types of knowledge, perspectives, and methods from academic and non-academic participants to create novel solutions to complex problems), and collaborative research program jointly led by Waipapa Taumata Rau University of Auckland and Earth Sciences New Zealand (ESNZ; formerly GNS Science), with core funding from Natural Hazards Commission Toka Tū Ake (NHC; formerly the Earthquake Commission, EQC) and Te Kaunihera o Tāmaki Makaurau Auckland Council (AC). The primary research focus of DEVORA is to investigate the geologic history, volcanic hazards, and risk posed by the basaltic intraplate Auckland Volcanic Field. Disruption from ash fall and gas from other Aotearoa New Zealand volcanoes is also considered. DEVORA’s work to explore exposure and vulnerability in Tāmaki Makaurau Auckland is also useful for assessing risks from other non-volcanic natural hazards, such as seismic and tsunami hazards. The greater Tāmaki Makaurau Auckland region has an ethnically and socio-economically diverse population of approximately 1.7 million, representing about one-third of the Aotearoa New Zealand population, and hosts critical infrastructure of national significance. The size and nature of the populace, consequential economic base, and important infrastructure within Tāmaki Makaurau Auckland mean that the effects of a volcanic eruption would be felt nationally, including through the disruption of air travel to Aotearoa New Zealand. The hazards from such an eruption could potentially affect hundreds of thousands of people, businesses, and lifelines (critical infrastructure). A considerable challenge for emergency and risk managers is the monogenetic nature of the volcanic field. It is not known where or when the next eruption will occur, how much warning we may get before an eruption, nor how an eruption and its effects might unfold. In this contribution, we highlight the concept and collaborative intent of the DEVORA Programme and show how it has evolved over the 16 years since its inception. We describe how DEVORA has unified more than 100 researchers (including more than 50 graduate students) and numerous stakeholders to address key issues facing Tāmaki Makaurau Auckland and describe how research findings are being implemented into policy and communicated to stakeholder agencies and the public. We also illustrate the broader influence of the DEVORA Programme and provide some learnings that might benefit others embarking on similar integrated projects, especially those focused on distributed volcanism in and near populated areas.

Auckland↗

USGS2024

Advancing at-risk species recovery planning in an era of rapid ecological change with a transparent, flexible, and expert-engaged approach

In the face of unprecedented ecological changes, the conservation community needs strategies to recover species at risk of extinction. On the Island of Maui, we collaborated with species experts and managers to assist with climate-resilient recovery planning for 36 at-risk native plant species by identifying priority areas for the management of recovery populations. To do this, we developed a tailored spatial conservation prioritization (SCP) approach distinguished by its emphasis on transparency, flexibility, and expert (TFE) engagement. Our TFE SCP approach consisted of 2 iterative steps: first, the generation of multiple candidate conservation footprints (i.e., prioritization solutions) with a flexible greedy algorithm that reflects conservation practitioners’ priorities and, second, the selection of an optimal conservation footprint based on the consideration of trade-offs in expert-agreed criteria among footprints. This process maximized buy-in by involving conservation practitioners and experts throughout, from setting goals to reviewing optimization data, defining optimization rules, and designating planning units meaningful to practitioners. We minimized the conservation footprint area necessary to meet recovery goals while incorporating species-specific measures of habitat suitability and climate resilience and retaining species-specific information for guiding recovery efforts. Our approach reduced the overall necessary conservation area by 36%, compared with selecting optimal recovery habitats for each species separately, and still identified high-quality habitat for individual species. Compared with prioritizr (an existing SCP tool), our approach identified a conservation area of equal size but with higher quality habitat. By integrating the strengths of existing techniques in a flexible and transparent design, our approach can address natural resource management constraints and provide outputs suitable for local recovery planning, consequently enhancing engagement and buy-in from conservation practitioners and experts. It demonstrates a step forward in making conservation planning more responsive to real-world complexities and helps reduce barriers to implementation for local conservation practitioners.

Hawaii↗

USGS1972

The ground-water system in southeastern Laramie County, Wyoming

Increased development of irrigation wells in southeastern Laramie County, Wyo., has caused concern about the quantity of water available. Ground water from approximately 230 large-capacity wells is used to irrigate most of the 18,165 acres under irrigation. The purpose of this study is to provide more knowledge about the character of the aquifers, quantity of water in storage, rate of withdrawal, and the effect of withdrawals on streamflow. The area studied consists of about 400 square miles in southeastern Laramie County in the extreme southeast corner of Wyoming. The White River Formation of Oligocene age and alluvium of Quaternary age are the principal aquifers. The White River Formation is made up primarily of clay, silt, and fine sand. Secondary permeability in the White River Formation accounts for it being an important aquifer. The alluvium, which Includes terrace and flood-plain deposits, consists of sand and gravel that contain some lenses of silt and clay. Existence of secondary permeability in the White River Formation has been accepted for some time although the nature of the secondary permeability has been disputed. Examination of downhole conditions with a television camera during this study revealed openings in the formation that appeared to be similar to tubes or caverns. The openings were of various sizes and shapes but only a few appeared to be associated with fracturing. Solution activity in the formation probably is an important factor in the development of secondary permeability. The study area was divided into the Pine Bluffs-Egbert area and the Carpenter area. Ground-water movement in the Pine Bluffs-Egbert area is generally eastward into Nebraska; in the Carpenter area, movement is generally southward into Colorado. Pumpage from large-capacity wells in the Pine Bluffs-Egbert area was estimated to be about 21,790 acre-feet in 1971. Water levels exhibited a declining trend annually in some areas during the period of record. Data indicate that pumpage in the Pine Bluffs-Egbert area probably is the cause of decreased base flow in Lodgepole Creek since approximately 1961. Increased pumpage, above that in 1971, will result in further reduction of discharge of Lodgepole Creek into Nebraska. In the Carpenter area, it was estimated that the terrace deposit contained about 1 million acre-feet of saturated sediments in March 1971. The amount of ground water in storage in the White River Formation in this area is unknown. Pumpage from large-capacity wells in the Carpenter area was estimated to be about 7,090 acre-feet in 1971. Recharge to the area was not estimated but there was no indication of a net decline in water levels between September 1970 and September 1971.

Open-File Report↗

USGS2013

Uniform California earthquake rupture forecast, version 3 (UCERF3): the time-independent model

In this report we present the time-independent component of the Uniform California Earthquake Rupture Forecast, Version 3 (UCERF3), which provides authoritative estimates of the magnitude, location, and time-averaged frequency of potentially damaging earthquakes in California. The primary achievements have been to relax fault segmentation assumptions and to include multifault ruptures, both limitations of the previous model (UCERF2). The rates of all earthquakes are solved for simultaneously, and from a broader range of data, using a system-level "grand inversion" that is both conceptually simple and extensible. The inverse problem is large and underdetermined, so a range of models is sampled using an efficient simulated annealing algorithm. The approach is more derivative than prescriptive (for example, magnitude-frequency distributions are no longer assumed), so new analysis tools were developed for exploring solutions. Epistemic uncertainties were also accounted for using 1,440 alternative logic tree branches, necessitating access to supercomputers. The most influential uncertainties include alternative deformation models (fault slip rates), a new smoothed seismicity algorithm, alternative values for the total rate of M≥5 events, and different scaling relationships, virtually all of which are new. As a notable first, three deformation models are based on kinematically consistent inversions of geodetic and geologic data, also providing slip-rate constraints on faults previously excluded because of lack of geologic data. The grand inversion constitutes a system-level framework for testing hypotheses and balancing the influence of different experts. For example, we demonstrate serious challenges with the Gutenberg-Richter hypothesis for individual faults. UCERF3 is still an approximation of the system, however, and the range of models is limited (for example, constrained to stay close to UCERF2). Nevertheless, UCERF3 removes the apparent UCERF2 overprediction of M6.5–7 earthquake rates and also includes types of multifault ruptures seen in nature. Although UCERF3 fits the data better than UCERF2 overall, there may be areas that warrant further site-specific investigation. Supporting products may be of general interest, and we list key assumptions and avenues for future model improvements.

California↗

USGS1988

Radarclinometry: Bootstrapping the radar reflectance function from the image pixel-signal frequency distribution and an altimetry profile

A method is derived for determining the dependence of radar backscatter on incidence angle that is applicable to the region corresponding to a particular radar image. The method is based on enforcing mathematical consistency between the frequency distribution of the image's pixel signals (histogram of DN values with suitable normalizations) and a one-dimensional frequency distribution of slope component, as might be obtained from a radar or laser altimetry profile in or near the area imaged. In order to achieve a unique solution, the auxiliary assumption is made that the two-dimensional frequency distribution of slope is isotropic. The backscatter is not derived in absolute units. The method is developed in such a way as to separate the reflectance function from the pixel-signal transfer characteristic. However, these two sources of variation are distinguishable only on the basis of a weak dependence on the azimuthal component of slope; therefore such an approach can be expected to be ill-conditioned unless the revision of the transfer characteristic is limited to the determination of an additive instrumental background level. The altimetry profile does not have to be registered in the image, and the statistical nature of the approach minimizes pixel noise effects and the effects of a disparity between the resolutions of the image and the altimetry profile, except in the wings of the distribution where low-number statistics preclude accuracy anyway. The problem of dealing with unknown slope components perpendicular to the profiling traverse, which besets the one-to-one comparison between individual slope components and pixel-signal values, disappears in the present approach. In order to test the resulting algorithm, an artificial radar image was generated from the digitized topographic map of the Lake Champlain West quadrangle in the Adirondack Mountains, U.S.A., using an arbitrarily selected reflectance function. From the same map, a one-dimensional frequency distribution of slope component was extracted. The algorithm recaptured the original reflectance function to the degree that, for the central 90% of the data, the discrepancy translates to a RMS slope error of 0.1 ???. For the central 99% of the data, the maximum error translates to 1 ???; at the absolute extremes of the data the error grows to 6 ???. ?? 1988 Kluwer Academic Publishers.

Earth, Moon and Planets↗

USGS2009

Thermal conductivity of hydrate-bearing sediments

A thorough understanding of the thermal conductivity of hydrate-bearing sediments is necessary for evaluating phase transformation processes that would accompany energy production from gas hydrate deposits and for estimating regional heat flow based on the observed depth to the base of the gas hydrate stability zone. The coexistence of multiple phases (gas hydrate, liquid and gas pore fill, and solid sediment grains) and their complex spatial arrangement hinder the a priori prediction of the thermal conductivity of hydrate-bearing sediments. Previous studies have been unable to capture the full parameter space covered by variations in grain size, specific surface, degree of saturation, nature of pore filling material, and effective stress for hydrate-bearing samples. Here we report on systematic measurements of the thermal conductivity of air dry, water- and tetrohydrofuran (THF)-saturated, and THF hydrate–saturated sand and clay samples at vertical effective stress of 0.05 to 1 MPa (corresponding to depths as great as 100 m below seafloor). Results reveal that the bulk thermal conductivity of the samples in every case reflects a complex interplay among particle size, effective stress, porosity, and fluid-versus-hydrate filled pore spaces. The thermal conductivity of THF hydrate–bearing soils increases upon hydrate formation although the thermal conductivities of THF solution and THF hydrate are almost the same. Several mechanisms can contribute to this effect including cryogenic suction during hydrate crystal growth and the ensuing porosity reduction in the surrounding sediment, increased mean effective stress due to hydrate formation under zero lateral strain conditions, and decreased interface thermal impedance as grain-liquid interfaces are transformed into grain-hydrate interfaces.

Journal of Geophysical Research B: Solid Earth↗

USGS2007

Ensuring confidence in radionuclide-based sediment chronologies and bioturbation rates

Sedimentary records of naturally occurring and fallout-derived radionuclides are widely used as tools for estimating both the ages of recent sediments and rates of sedimentation and bioturbation. Developing these records to the point of data interpretation requires careful sample collection, processing, analysis and data modeling. In this work, we document a number of potential pitfalls that can impact sediment core records and their interpretation. This paper is not intended as an exhaustive treatment of these potential problems. Rather, the emphasis is on potential problems that are not well documented in the literature, as follows: (1) the mere sampling of sediment cores at a resolution that is too coarse can result in an apparent diffusive mixing of the sedimentary record at rates comparable to diffusive bioturbation rates observed in many locations; (2) 210 Pb profiles in slowly accumulating sediments can easily be misinterpreted to be driven by sedimentation, when in fact bioturbation is the dominant control. Multiple isotopes of different half lives and/or origin may help to distinguish between these two possible interpretations; (3) apparent mixing can occur due simply to numerical artifacts inherent in the finite difference approximations of the advection diffusion equation used to model sedimentation and bioturbation. Model users need to be aware of this potential problem. Solutions to each of these potential pitfalls are offered to ensure the best possible sediment age estimates and/or sedimentation and bioturbation rates can be obtained.

Estuarine, Coastal and Shelf Science↗

USGS2021

SiteOpt: An open-source R-package for site selection and portfolio optimization

Conservation planning involves identifying and selecting actions to best achieve objectives for managing natural, social and cultural resources. Conservation problems are often high dimensional when specified as combinatorial or portfolio problems and when multiple competing objectives are considered at varying spatial and temporal scales. Although analytical techniques such as modern portfolio theory (MPT) have been developed to address these complex problems, open source computational platforms for executing these approaches are not readily available. We present a user-friendly R-package called SiteOpt for optimization of binary decisions while explicitly considering environmental or economic uncertainty and the risk tolerance of decision makers. We illustrate the package with spatially-explicit site selection problems (i.e. spatial conservation planning), including an option for divestment (i.e. selling assets), when accounting for future uncertainties in designing conservation areas. The tool is applicable to both spatial and non-spatial problems, such as budget allocation or species selection. Constraints for spatial design and spatial dependencies (e.g. connectivity among sites) can also be specified in SiteOpt. Users can optimize site selection based on two competing objectives by solving for the Nash bargaining solution. Importantly, by quantifying uncertainty and asset spatial correlation, a measure of risk can be included as one such objective to be traded off against portfolio benefits. Thus, SiteOpt can be used to explicitly manage risk in portfolio-based spatial optimization. This tool facilitates decisions in a variety of problem settings, including reserve selection, invasive species management, allocation of law enforcement activities for conservation, budget allocation and asset selection under uncertainty and risk.

Ecography↗

USGS2024

Grassland Effectiveness Monitoring (GEM): A tiered approach for habitat treatment assessment across private lands incentive programs

Introduction The decline of North American grasslands is a topic of increasing interest as agencies and organizations work to address subsequent declines in wildlife species, including grassland birds (Rosenberg et al. 2019), pronghorn (Antilocapra americana) (Gedir et al. 2015), and other grassland-dependent taxa. In response to grassland habitat loss, conservation programs and policies have been developed to provide biologists and landowners mechanisms to restore grassland habitat on private lands. These range from federal programs such as the U.S. Fish and Wildlife Service’s (USFWS) Partners for Fish and Wildlife program (PFW) and the Natural Resources Conservation Service’s (NRCS) Environmental Quality Incentive Program, state programs such as Texas Parks and Wildlife Department’s Pastures for Upland Birds, and partnership-based programs such as the Oaks and Prairies Joint Venture’s Grassland Restoration Incentive Program (GRIP). When managing and restoring grassland habitat, these programs typically utilize the same set of practices: prescribed grazing, prescribed fire, herbicide, brush management, and range planting. To restore grasslands at scale, large sums of funding have been invested into these programs. For example, approximately $1 million is spent by PFW in Texas annually (D. Wilhelm, USFWS, personal communication). However, it is rare for these conservation programs to have associated monitoring efforts that are consistently used to assess the effectiveness of the aforementioned practices in achieving program and project objectives. As such, it is difficult to communicate to funders and the public whether these conservation programs are successfully addressing grassland decline. Similar objectives and practices across programs, however, provide a unique opportunity for developing an innovative monitoring methodology that focuses on providing a collaborative solution to this dilemma. Our objective is to develop a grassland monitoring program that can be used across agencies, organizations, and conservation programs to assess the effectiveness of practices and programs in restoring healthy grasslands. In addition, we seek to develop a program that maintains consistency to allow for regional to national-scale reporting while also allowing flexibility for achieving local and partner monitoring objectives. Our vision is for this program to be used across North American grasslands to improve our ability to address the questions: “Is it actually working?” and “How can we improve habitat delivery?”

Conference Paper↗

USGS2007

Stability and uncertainty of finite-fault slip inversions: Application to the 2004 Parkfield, California, earthquake

The 2004 Parkfield, California, earthquake is used to investigate stability and uncertainty aspects of the finite-fault slip inversion problem with different a priori model assumptions. We utilize records from 54 strong ground motion stations and 13 continuous, 1-Hz sampled, geodetic instruments. Two inversion procedures are compared: a linear least-squares subfault-based methodology and a nonlinear global search algorithm. These two methods encompass a wide range of the different approaches that have been used to solve the finite-fault slip inversion problem. For the Parkfield earthquake and the inversion of velocity or displacement waveforms, near-surface related site response (top 100 m, frequencies above 1 Hz) is shown to not significantly affect the solution. Results are also insensitive to selection of slip rate functions with similar duration and to subfault size if proper stabilizing constraints are used. The linear and nonlinear formulations yield consistent results when the same limitations in model parameters are in place and the same inversion norm is used. However, the solution is sensitive to the choice of inversion norm, the bounds on model parameters, such as rake and rupture velocity, and the size of the model fault plane. The geodetic data set for Parkfield gives a slip distribution different from that of the strong-motion data, which may be due to the spatial limitation of the geodetic stations and the bandlimited nature of the strong-motion data. Cross validation and the bootstrap method are used to set limits on the upper bound for rupture velocity and to derive mean slip models and standard deviations in model parameters. This analysis shows that slip on the northwestern half of the Parkfield rupture plane from the inversion of strong-motion data is model dependent and has a greater uncertainty than slip near the hypocenter.

California↗

USGS2007

Revisiting the cape cod bacteria injection experiment using a stochastic modeling approach

Bromide and resting-cell bacteria tracer tests conducted in a sandy aquifer at the U.S. Geological Survey Cape Cod site in 1987 were reinterpreted using a three-dimensional stochastic approach. Bacteria transport was coupled to colloid filtration theory through functional dependence of local-scale colloid transport parameters upon hydraulic conductivity and seepage velocity in a stochastic advection - dispersion/attachment - detachment model. Geostatistical information on the hydraulic conductivity (K) field that was unavailable at the time of the original test was utilized as input. Using geostatistical parameters, a groundwater flow and particle-tracking model of conservative solute transport was calibrated to the bromide-tracer breakthrough data. An optimization routine was employed over 100 realizations to adjust the mean and variance ofthe natural-logarithm of hydraulic conductivity (InK) field to achieve best fit of a simulated, average bromide breakthrough curve. A stochastic particle-tracking model for the bacteria was run without adjustments to the local-scale colloid transport parameters. Good predictions of mean bacteria breakthrough were achieved using several approaches for modeling components of the system. Simulations incorporating the recent Tufenkji and Elimelech (Environ. Sci. Technol. 2004, 38, 529-536) correlation equation for estimating single collector efficiency were compared to those using the older Rajagopalan and Tien (AIChE J. 1976, 22, 523-533) model. Both appeared to work equally well at predicting mean bacteria breakthrough using a constant mean bacteria diameter for this set of field conditions. Simulations using a distribution of bacterial cell diameters available from original field notes yielded a slight improvement in the model and data agreement compared to simulations using an average bacterial diameter. The stochastic approach based on estimates of local-scale parameters for the bacteria-transport process reasonably captured the mean bacteria transport behavior and calculated an envelope of uncertainty that bracketed the observations in most simulation cases.

Massachusetts↗

USGS2009

Atmospheric correction at AERONET locations: A new science and validation data set

This paper describes an Aerosol Robotic Network (AERONET)-based Surface Reflectance Validation Network (ASRVN) and its data set of spectral surface bidirectional reflectance and albedo based on Moderate Resolution Imaging Spectroradiometer (MODIS) TERRA and AQUA data. The ASRVN is an operational data collection and processing system. It receives 50 ?? 50 km 2 ; subsets of MODIS level 1B (L1B) data from MODIS adaptive processing system and AERONET aerosol and water-vapor information. Then, it performs an atmospheric correction (AC) for about 100 AERONET sites based on accurate radiative-transfer theory with complex quality control of the input data. The ASRVN processing software consists of an L1B data gridding algorithm, a new cloud-mask (CM) algorithm based on a time-series analysis, and an AC algorithm using ancillary AERONET aerosol and water-vapor data. The AC is achieved by fitting the MODIS top-of-atmosphere measurements, accumulated for a 16-day interval, with theoretical reflectance parameterized in terms of the coefficients of the Li SparseRoss Thick (LSRT) model of the bidirectional reflectance factor (BRF). The ASRVN takes several steps to ensure high quality of results: 1) the filtering of opaque clouds by a CM algorithm; 2) the development of an aerosol filter to filter residual semitransparent and subpixel clouds, as well as cases with high inhomogeneity of aerosols in the processing area; 3) imposing the requirement of the consistency of the new solution with previously retrieved BRF and albedo; 4) rapid adjustment of the 16-day retrieval to the surface changes using the last day of measurements; and 5) development of a seasonal backup spectral BRF database to increase data coverage. The ASRVN provides a gapless or near-gapless coverage for the processing area. The gaps, caused by clouds, are filled most naturally with the latest solution for a given pixel. The ASRVN products include three parameters of the LSRT model (kL, kG, and kV), surface albedo, normalized BRF (computed for a standard viewing geometry, VZA = 0, SZA = 45??), and instantaneous BRF (or one-angle BRF value derived from the last day of MODIS measurement for specific viewing geometry) for the MODIS 500-m bands 17. The results are produced daily at a resolution of 1 km in gridded format. We also provide a cloud mask, a quality flag, and a browse bitmap image. The ASRVN data set, including 6 years of MODIS TERRA and 1.5 years of MODIS AQUA data, is available now as a standard MODIS product (MODASRVN) which can be accessed through the Level 1 and Atmosphere Archive and Distribution System website ( http://ladsweb.nascom.nasa.gov/data/search.html). It can be used for a wide range of applications including validation analysis and science research. ?? 2006 IEEE.

IEEE Transactions on Geoscience and Remote Sensing↗

USGS2026

Developing a low-cost drone-based method for deploying and retrieving autonomous recording units in inaccessible areas

Autonomous Recording Units (ARUs) are increasingly used in research to support passive acoustic monitoring, but deployment in remote or inaccessible locations remains logistically challenging. Traditional manual placement is labor-intensive, time-consuming, potentially hazardous, and often limited to habitat edges, creating sampling biases and restricting spatial coverage. To address placement limitations in wetland systems, we developed and field-tested a low-cost, lightweight floating platform paired with a drone-based deployment and retrieval system. The drone-deployable ARU platform, constructed from inexpensive, off-the-shelf materials (~US $21 per unit), weighed ~560 g and was mounted with an ARU (AudioMoth) at ~1.2 m above water. The platform was designed for stability, portability, and compliance with Federal Aviation Administration Part 107 regulations. Field trials were conducted across 50 operations in various types of historical rice-field impoundments along the South Carolina coast, which provide critical habitat for secretive marsh birds such as rails (e.g ., Rallus, Laterallus , and Porzana ) and bitterns (e.g ., Botaurus and Ixobrychus ). All 50 deployments and retrievals were successful, demonstrating the robustness of the method under diverse hydrologic and vegetative conditions. Deployment times were significantly shorter than retrieval times (median 6 vs. 9 min), with retrieval requiring greater precision for hook engagement. Operation times scaled predictably with distance: deployment increased by 0.93 min/100 m (R 2 = 0.94) and retrieval by 1.17 min/100 m (R 2 = 0.95). No platform damage or displacement was observed, as emergent vegetation was likely providing natural anchoring. Our approach offers an affordable and effective solution for expanding ARU coverage in large, inaccessible wetlands, reducing sampling bias and enhancing biodiversity monitoring. With continued advances in drone payload capacity, battery endurance, and beyond-visual-line-of-sight flight regulations, our workflow could be adapted for deploying multi-sensor platforms (e.g., ARUs, eDNA samplers, water and air quality sensors, plant sample collectors, insect traps, and trail cameras) to support integrated biodiversity monitoring in challenging landscapes worldwide.

South Carollina↗