Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Modeling Report”

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 829 records · Page 46Linked to original sources

Piloting urban ecosystem accounting for the United States

In this study, we develop urban ecosystem accounts in the U.S., using the System of Environmental-Economic Accounting Experimental Ecosystem Accounting (SEEA EEA) framework. Most ecosystem accounts focus on regional and national scales, which are appropriate for many ecosystem services. However, ecosystems provide substantial services in cities, improving quality of life and contributing to resiliency for substantial parts of the population. Our models estimate energy savings for indoor cooling resulting from heat mitigated by trees and rainfall intercepted by trees. Both models cover major cities in the contiguous U.S. and report the results through physical supply and use tables for multiple accounting periods (2011 and 2016). Using conservative assumptions, urban trees provide substantial heat mitigation (4,098 and 4,229 GWh, valued at $523 and $539 million in 2011 and 2016, respectively) and rainfall interception (2,422 and 2,627 million m 3 , valued at $434 and $425 million for 2011 and 2016, respectively). Interannual differences largely reflect variations in weather patterns. Our work shows how Earth observation data can support urban ecosystem accounting. We provide model code within a public repository to facilitate model runs elsewhere, enabling the SEEA EEA and Earth observation user communities to reuse our models and provide feedback for improvement.

Ecosystem Services↗

Trends of Abutment-Scour Prediction Equations Applied to 144 Field Sites in South Carolina

The U.S. Geological Survey conducted a study in cooperation with the Federal Highway Administration in which predicted abutment-scour depths computed with selected predictive equations were compared with field measurements of abutment-scour depth made at 144 bridges in South Carolina. The assessment used five equations published in the Fourth Edition of 'Evaluating Scour at Bridges,' (Hydraulic Engineering Circular 18), including the original Froehlich, the modified Froehlich, the Sturm, the Maryland, and the HIRE equations. An additional unpublished equation also was assessed. Comparisons between predicted and observed scour depths are intended to illustrate general trends and order-of-magnitude differences for the prediction equations. Field measurements were taken during non-flood conditions when the hydraulic conditions that caused the scour generally are unknown. The predicted scour depths are based on hydraulic conditions associated with the 100-year flow at all sites and the flood of record for 35 sites. Comparisons showed that predicted scour depths frequently overpredict observed scour and at times were excessive. The comparison also showed that underprediction occurred, but with less frequency. The performance of these equations indicates that they are poor predictors of abutment-scour depth in South Carolina, and it is probable that poor performance will occur when the equations are applied in other geographic regions. Extensive data and graphs used to compare predicted and observed scour depths in this study were compiled into spreadsheets and are included in digital format with this report. In addition to the equation-comparison data, Water-Surface Profile Model tube-velocity data, soil-boring data, and selected abutment-scour data are included in digital format with this report. The digital database was developed as a resource for future researchers and is especially valuable for evaluating the reasonableness of future equations that may be developed.

South Carolina↗

Application and evaluation of a rapid response earthquake-triggered landslide model to the 25 April 2015 Mw 7.8 Gorkha earthquake, Nepal

The 25 April 2015 M w 7.8 Gorkha earthquake produced strong ground motions across an approximately 250 km by 100 km swath in central Nepal. To assist disaster response activities, we modified an existing earthquake-triggered landslide model based on a Newmark sliding block analysis to estimate the extent and intensity of landsliding and landslide dam hazard. Landslide hazard maps were produced using Shuttle Radar Topography Mission (SRTM) digital topography, peak ground acceleration (PGA) information from the U.S. Geological Survey (USGS) ShakeMap program, and assumptions about the regional rock strength based on end-member values from previous studies. The instrumental record of seismicity in Nepal is poor, so PGA estimates were based on empirical Ground Motion Prediction Equations (GMPEs) constrained by teleseismic data and felt reports. We demonstrate a non-linear dependence of modeled landsliding on aggregate rock strength, where the number of landslides decreases exponentially with increasing rock strength. Model estimates are less sensitive to PGA at steep slopes (> 60°) compared to moderate slopes (30–60°). We compare forward model results to an inventory of landslides triggered by the Gorkha earthquake. We show that moderate rock strength inputs over estimate landsliding in regions beyond the main slip patch, which may in part be related to poorly constrained PGA estimates for this event at far distances from the source area. Directly above the main slip patch, however, the moderate strength model accurately estimates the total number of landslides within the resolution of the model (landslides ≥ 0.0162 km 2 ; observed n = 2214, modeled n = 2987), but the pattern of landsliding differs from observations. This discrepancy is likely due to the unaccounted for effects of variable material strength and local topographic amplification of strong ground motion, as well as other simplifying assumptions about source characteristics and their relationship to landsliding.

Tectonophysics↗

On the use of high-resolution and deep-learning seismic catalogs for short-term earthquake forecasts: Potential benefits and current limitations

Enhanced earthquake catalogs provide detailed images of evolving seismic sequences. Currently, these data sets take some time to be released but will soon become available in real time. Here, we explore whether and how enhanced seismic catalogs feeding into established short-term earthquake forecasting protocols may result in higher predictive skill. We consider three enhanced catalogs for the 2016–2017 Central Italy sequence, featuring a bulk completeness lower by at least two magnitude units compared to the real-time catalog and an improved hypocentral resolution. We use them to inform a set of physical Coulomb Rate-and-State (CRS) and statistical Epidemic-Type Aftershock Sequence (ETAS) models to forecast the space-time occurrence of M3+ events during the first 6 months of the sequence. We track model performance using standard likelihood-based metrics and compare their skill against the best-performing CRS and ETAS models among those developed with the real-time catalog. We find that while the incorporation of the triggering contributions from new small magnitude detections of the enhanced catalogs is beneficial for both types of forecasts, these models do not significantly outperform their respective near real-time benchmarks. To explore the reasons behind this result, we perform targeted sensitivity tests that show how (a) the typical spatial discretizations of forecast experiments ( 2 km) hamper the ability of models to capture highly localized secondary triggering patterns and (b) differences in earthquake parameters (i.e., magnitude and hypocenters) reported in different catalogs can affect forecast evaluation. These findings will contribute toward improving forecast model design and evaluation strategies for next-generation seismic catalogs.

Journal of Geophysical Research--Solid Earth↗

Development of a hydrogeologic visualization model for western Sarpy County, Nebraska

Population in western Sarpy County, Nebraska, has steadily increased over the last several decades and has led to increased groundwater use for domestic purposes. To meet the increase in demand, the Papio-Missouri River Natural Resources District is seeking to use all available sources of groundwater in western Sarpy County. Additionally, elevated groundwater nitrate plus nitrite as nitrogen concentrations were detected, indicating the need to better understand the groundwater quality of the area. Although the general geology of the area is understood, the area does not have detailed information on the extent of the various aquifers, particularly the Dakota aquifer. To characterize these aquifers, the Papio-Missouri River Natural Resources District invested in airborne electromagnetic surveys of the area to better understand the subsurface geology. Although these surveys improved understanding of the groundwater systems in the area, the Papio-Missouri River Natural Resources District wanted to integrate the subsurface information with available water-quality and groundwater-level data. In response, the U.S. Geological Survey, in cooperation with the Papio-Missouri River Natural Resources District, the Nebraska Natural Resources Commission, and the Nebraska Department of Natural Resources, assembled geologic, hydrogeologic and nitrate plus nitrite as nitrogen information for the selected area into a three-dimensional visualization computer software package called GeoScene3D. The completed GeoScene3D project was assembled to provide a visualization of the groundwater systems and associated water-quality results in Sarpy County and to provide the Papio-Missouri River Natural Resources District managers with information that can be used to make more informed groundwater resource-planning decisions in the future. This report details the development of a three-dimensional model created within GeoScene3D to visualize the subsurface, particularly the Dakota Sandstone in western Sarpy County.

Nebraska↗

Simulation of an urban ground-water-flow system in the Menomonee Valley, Milwaukee, Wisconsin using analytic element modeling

A single-layer, steady-state analytic element model was constructed to simulate shallow ground-water flow in the Menomonee Valley, an old industrial center southwest of downtown Milwaukee, Wisconsin. Project objectives were to develop an understanding of the shallow ground-water flow system and identify primary receptors of recharge to the valley. The analytic element model simulates flow in a 18.3 m (60 ft) thick layer of estuarine and alluvial sediments and man-made fill that comprises the shallow aquifer across the valley. The thin, laterally extensive nature of the shallow aquifer suggests horizontal-flow predominates, thus the system can appropriately be modeled with the Dupuit-Forchheimer approximation in an analytic element model. The model was calibrated to the measured baseflow increase between two USGS gages on the Menomonee River, 90 head measurements taken in and around the valley during December 1999, and vertical gradients measured at five locations under the river and estuary in the valley. Recent construction of the Milwaukee Metropolitan Sewer District Inline Storage System (ISS) in the Silurian dolomite under the Menomonee Valley has locally lowered heads in the dolomite appreciably, below levels caused by historic pumping. The ISS is a regional hydraulic sink which removes water from the bedrock even during dry weather. The potential effect on flow directions in the shallow aquifer of dry-weather infiltration to the ISS was evaluated by adjusting the resistance of the line-sink strings representing the ISS in the model to allow infiltration from 0 to 100% of the reported 9,500 m 3 /d. The best fit to calibration targets was found between 60% (5,700 m 3 /d) and 80% (7,600 m 3 /d) of the reported dry-weather infiltration. At 60% infiltration, 65% of the recharge falling on the valley terminates at the ISS and 35% at the Menomonee River and estuary. At 80% infiltration, 73% of the recharge terminates at the ISS, and 27% at the river and estuary. Model simulations suggest that the ISS has an greater influence on the shallow ground-water flow in the eastern half of valley as compared to the western half. Preliminary three-dimensional simulations using the numerical MODFLOW code show good agreement with the single-layer simulation and supports its use in evaluating the shallow system. Copyright ASCE 2004.

Conference Paper↗

The Pothole Hydrology-Linked Systems Simulator (PHyLiSS)—Development and application of a systems model for prairie-pothole wetlands

The North American Prairie Pothole Region covers about 770,000 square kilometers of the United States and Canada (including parts of 5 States and 3 provinces: North Dakota, South Dakota, Montana, Minnesota, Iowa, Saskatchewan, Manitoba, and Alberta). The Laurentide Ice Sheet shaped the landscape of the region about 12,000 to 14,000 years ago. The retreat of the ice sheet left behind low-permeability glacial till and a landscape dotted with millions of depressions known today as prairie potholes. The wetlands that subsequently formed in these depressions, prairie-pothole wetlands, provide critical migratory-bird habitat and support dynamic aquatic communities. Extensive grasslands and productive agricultural systems surround these wetland ecosystems. In prairie-pothole wetlands, the compositions of plant, invertebrate, and vertebrate communities are highly dependent on hydrogeochemical conditions. Regional climate shifts between wet and dry periods affect the length of time that wetlands contain ponded surface water and the chemistry of that ponded water. Land-use change can exacerbate or reduce the effects of climate on wetland hydrology and water chemistry. A mechanistic understanding of the relation among climate, land use, hydrology, chemistry, and biota in prairie-pothole wetlands is needed to better understand the complex, and often interacting, effects of climate and land use on prairie-pothole wetland systems and to facilitate climate and land-use change adaptation efforts. The Pothole Hydrology-Linked Systems Simulator (PHyLiSS) model was developed to address this need. The model simulates water-surface elevation dynamics in prairie-pothole wetlands and quantifies changes in salinity. The PHyLiSS model is unique among other wetland models because it accommodates differing sizes and morphometries of wetland basins, is not dependent on a priori designations of wetland class, and allows for functional changes associated with dynamic shifts in ecohydrological states. The PHyLiSS model also has the capability to simulate wetland salinity, and potential future iterations will also simulate the effects of changing hydrology and geochemical conditions on biota. This report documents the development of the hydrological and geochemical components of the PHyLiSS model and provides example applications.

Open-File Report↗

A comparison of eDNA and visual survey methods for detection of longnose darter Percina nasuta in Missouri

The longnose darter Percina nasuta is a rare and cryptic fish that recently disappeared from much of its historic range. We developed and used an environmental DNA (eDNA) assay for longnose darter paired with visual surveys to better determine the species’ range and compare detection probability between sampling approaches in an occupancy modeling framework. We detected longnose darter eDNA further upstream in the mainstem St. Francis River than previously reported and in a tributary for the first time. Our multi-scale occupancy approach compared models where detection was constant against a model that allowed detection to vary by survey method. The constant model received the most support indicating survey method was not a strong predictor and detection was estimated at 0.70 (0.45–0.86; 95% CI) across both methods. Our study produced effective longnose darter eDNA primers and demonstrated the application of eDNA for sampling small-bodied, cryptic fish. We detected longnose darter eDNA 27 km upstream of their known range and determined that snorkel surveys are the most efficient sampling method if water clarity allows. We recommend target sample sizes to achieve various detection goals for both sample methods and our results inform future design of distributional and monitoring efforts.

Missouri↗

Palaeoecology and sedimentology of the dysaerobic Bedford fauna (late Devonian), Ohio and Kentucky (USA)

Oxygen-deficient biofacies models rely on lithologic and paleontologic attributes to identify distinctive biofacies interpreted to reflect levels of oxygenation in anaerobic, dysaerobic, and aerobic parts of a stratified water column. This study of the Bedford fauna from the Bedford Shale of Ohio and Kentucky and from adjacent black-shale units reports faunal distributions different from those predicted by the accepted models. This study suggests that, although oxygenation was an important factor that determined the taxonomic makeup of the fauna, bacterially mediated nutrient recycling and substrate characteristics were more important than oxygenation in determining faunal distribution in the dysaerobic zone.

Kentucky, Ohio↗

Identifying the relative importance of water-budget information needed to quantify how land-cover change affects recharge, Hawaiian Islands

This report describes a sensitivity analysis of a water-budget model that was completed to identify the most important types of hydrologic information needed to reduce the uncertainty of model recharge estimates. The sensitivity of model recharge estimates for the Hawaiian Islands of Oʻahu and Maui was analyzed for seven model parameters potentially affected by land-cover changes within a watershed. The seven model parameters tested were canopy capacity, canopy-cover fraction, crop coefficient, fog-catch efficiency, root depth, stemflow, and trunk-storage capacity. Results of the sensitivity analysis were used to (1) quantify the relative importance of the seven model parameters to recharge assessments for three moisture zones (dry, mesic, and wet) on Oʻahu and Maui and (2) prepare a list of critical information needs for each moisture zone. The list of critical information needs was developed for three general types of land cover (forest, shrubland, and grassland) that are assumed to be affected by watershed management in the Hawaiian Islands. Identified information needs included estimates or measurements of (1) evapotranspiration processes needed to determine crop coefficients for land-cover types in all moisture zones, (2) rooting depths for land-cover types in the dry and mesic moisture zones, (3) canopy-cover fraction for forests in the wet and mesic moisture zones, (4) ratios of fog interception to rainfall for forests and shrublands in the wet moisture zone, and (5) canopy capacity for forests in the wet and mesic moisture zones. The list of information needs can guide data-collection strategies of future projects. Collection and analysis of the identified hydrologic information may help model users develop a better parameterization scheme, reduce uncertainty of values that model users assign to land-cover dependent parameters, and therefore allow future applications of the water-budget model to more accurately quantify how recharge in the Hawaiian Islands might be affected by future land-cover changes within a watershed.

Hawaii↗

Model-Based Predictions of the Effects of Harvest Mortality on Population Size and Trend of Yellow-Billed Loons

Yellow-billed loons (Gavia adamsii) breed in low densities in northern tundra habitats in Alaska, Canada, and Russia. They migrate to coastal marine habitats at mid to high latitudes where they spend their winters. Harvest may occur throughout the annual cycle, but of particular concern are recent reports of harvest from the Bering Strait region, which lies between Alaska and Russia and is an area used by yellow-billed loons during migration. Annual harvest for this region was reported to be 317, 45, and 1,077 during 2004, 2005, and 2007, respectively. I developed a population model to assess the effect of this reported harvest on population size and trend of yellow-billed loons. Because of the uncertainty regarding actual harvest and definition of the breeding population(s) affected by this harvest, I considered 25 different scenarios. Predicted trends across these 25 scenarios ranged from stability to rapid decline (24 percent per year) with halving of the population in 3 years. Through an assessment of literature and unpublished satellite tracking data, I suggest that the most likely of these 25 scenarios is one where the migrant population subjected to harvest in the Bering Strait includes individuals from breeding populations in Alaska (Arctic coastal plain and the Kotzebue region) and eastern Russia, and for which the magnitude of harvest varies among years and emulates the annual variation of reported harvest during 2004-07 (317, 45, and 1,077 yellow-billed loons). This scenario, which assumes no movement of Canadian breeders through the Bering Strait, predicts a 4.6 percent rate of annual population decline, which would halve the populations in 15 years. Although these model outputs reflect the best available information, confidence in these predictions and applicable scenarios would be greatly enhanced by more information on harvest, rates of survival and reproduction, and migratory pathways.

Open-File Report↗

Continuous hydrologic simulation and flood-frequency, hydraulic, and flood-hazard analysis of the Blackberry Creek watershed, Kane County, Illinois

Results of hydrologic model, flood-frequency, hydraulic model, and flood-hazard analysis of the Blackberry Creek watershed in Kane County, Illinois, indicate that the 100-year and 500-year flood plains range from approximately 25 acres in the tributary F watershed (a headwater subbasin at the northeastern corner of the watershed) to almost 1,800 acres in Blackberry Creek main stem. Based on 1996 land-cover data, most of the land in the 100-year and 500-year flood plains was cropland, forested and wooded land, and grassland. A relatively small percentage of urban land was in the flood plains. The Blackberry Creek watershed has undergone rapid urbanization in recent decades. The population and urbanized lands in the watershed are projected to double from the 1990 condition by 2020. Recently, flood-induced damage has occurred more frequently in urbanized areas of the watershed. There are concerns about the effect of urbanization on flood peaks and volumes, future flood-mitigation plans, and potential effects on the water quality and stream habitats. This report describes the procedures used in developing the hydrologic models, estimating the flood-peak discharge magnitudes and recurrence intervals for flood-hazard analysis, developing the hydraulic model, and the results of the analysis in graphical and tabular form. The hydrologic model, Hydrological Simulation Program-FORTRAN (HSPF), was used to perform the simulation of continuous water movements through various patterns of land uses in the watershed. Flood-frequency analysis was applied to an annual maximum series to determine flood quantiles in subbasins for flood-hazard analysis. The Hydrologic Engineering Center-River Analysis System (HEC-RAS) hydraulic model was used to determine the 100-year and 500-year flood elevations, and to determine the 100-year floodway. The hydraulic model was calibrated and verified using high water marks and observed inundation maps for the July 17-18, 1996, flood event. Digital maps of the 100-year and 500-year flood plains and the 100-year floodway for each tributary and the main stem of Blackberry Creek were compiled.

Scientific Investigations Report↗

Regional groundwater-flow model of the Lake Michigan Basin in support of Great Lakes Basin water availability and use studies

A regional groundwater-flow model of the Lake Michigan Basin and surrounding areas has been developed in support of the Great Lakes Basin Pilot project under the U.S. Geological Survey's National Water Availability and Use Program. The transient 2-million-cell model incorporates multiple aquifers and pumping centers that create water-level drawdown that extends into deep saline waters. The 20-layer model simulates the exchange between a dense surface-water network and heterogeneous glacial deposits overlying stratified bedrock of the Wisconsin/Kankakee Arches and Michigan Basin in the Lower and Upper Peninsulas of Michigan; eastern Wisconsin; northern Indiana; and northeastern Illinois. The model is used to quantify changes in the groundwater system in response to pumping and variations in recharge from 1864 to 2005. Model results quantify the sources of water to major pumping centers, illustrate the dynamics of the groundwater system, and yield measures of water availability useful for water-resources management in the region. This report is a complete description of the methods and datasets used to develop the regional model, the underlying conceptual model, and model inputs, including specified values of material properties and the assignment of external and internal boundary conditions. The report also documents the application of the SEAWAT-2000 program for variable-density flow; it details the approach, advanced methods, and results associated with calibration through nonlinear regression using the PEST program; presents the water-level, drawdown, and groundwater flows for various geographic subregions and aquifer systems; and provides analyses of the effects of pumping from shallow and deep wells on sources of water to wells, the migration of groundwater divides, and direct and indirect groundwater discharge to Lake Michigan. The report considers the role of unconfined conditions at the regional scale as well as the influence of salinity on groundwater flow. Lastly, it describes several categories of limitations and discusses ways of extending the regional model to address issues at the local scale. Results of the simulations portray a regional groundwater-flow system that, over time, has largely maintained its natural predevelopment configuration but that locally has been strongly affected by well withdrawals. The quantity of rainfall in the Lake Michigan Basin and adjacent areas supports a dense surface-water network and recharge rates consistent with generally shallow water tables and predominantly shallow groundwater flow. At the regional scale, pumping has not caused major modifications of the shallow flow system, but it has resulted in decreases in base flow to streams and in direct discharge to Lake Michigan (about 2 percent of the groundwater discharged and about 0.5 cubic foot per second per mile of shoreline). On the other hand, well withdrawals have caused major reversals in regional flow patterns around pumping centers in deep, confined aquifers - most noticeably in the Cambrian-Ordovician aquifer system on the west side of Lake Michigan near the cities of Green Bay and Milwaukee in eastern Wisconsin, and around Chicago in northeastern Illinois, as well as in some shallow bedrock aquifers (for example, in the Marshall aquifer near Lansing, Mich.). The reversals in flow have been accompanied by large drawdowns with consequent local decrease in storage. On the west side of Lake Michigan, groundwater withdrawals have caused appreciable migration of the deep groundwater divides. Before the advent of pumping, the deep Lake Michigan groundwater-basin boundaries extended west of the Lake Michigan surface-water basin boundary, in some places by tens of miles. Over time, the pumping centers have replaced Lake Michigan as the regional sink for the deep flow system. The regional model is intended to support the framework pilot study of water availability and use for the Great Lakes Basin (Reeves, in press).

Illinois, Indiana, Michigan, Wisconsin↗

Development of an online reporting format to facilitate the inclusion of ecosystem services into Conservation Reserve Enhancement Program reports

The Conservation Reserve Enhancement Program is a program administered by the U.S. Department of Agriculture’s Farm Service Agency. The Secretary of Agriculture is required to submit an annual report to Congress on Conservation Reserve Enhancement Program agreements that, among other things, reports on the progress made towards fulfilling commitments outlined in the agreements. The U.S. Geological Survey developed an online reporting form designed to ensure that consistent information is submitted to the Farm Service Agency from Conservation Reserve Enhancement Program State partners. Combined with the automated importation of text from partner-provided forms to word-processing documents, individual State reports and annual reports to Congress can now be produced efficiently and in a standardized format. Use of a standardized reporting format will also assist the Farm Service Agency in collecting information needed to support ecosystem service quantifications that go beyond the quantifications required from partners to document progress towards meeting the specific purposes and objectives identified in each agreement. Addition of these overarching conservation effect quantifications builds upon past ecosystem services modeling efforts based on the Integrated Valuation of Ecosystem Services and Tradeoffs suite of open-source software models; these offer a spatially explicit means to quantify additional ecosystem services across diverse partners in a consistent manner. Data sources are currently available to provide much of the information needed to run these models and complete simulations that would facilitate the quantification and reporting of the societal values of conservation actions taken under the Conservation Reserve Enhancement Program. It is the aim of this report to provide the information needed to move towards widescale monitoring of the Nation’s ecosystem services in a natural accounting framework, similar to the framework used to value financial and human capital.

Open-File Report↗

Vegetation projections for Wind Cave National Park with three future climate scenarios: Final report in completion of Task Agreement J8W07100052

Introduction The effects of climate change on the natural resources protected by Parks will likely be substantial, but geographically variable, due to local variation in climate trajectories and differences among ecosystems in their vulnerability to climate change. The projections of general circulation models (GCMs) indicate the possible magnitude and direction of future climate change for a region, but the utility of these projections for more local scales, those of individual National Park Service (NPS) units, are more uncertain because the coarse-scale GCMs lack much of the topographic detail that alters local climates. In addition, complex, interacting effects of temperature, precipitation, atmospheric CO 2 concentrations, fire, and herbivores on the vegetation that is the foundational natural resource of many NPS units present challenges in assessing the effects of projected future climates on plant and animal assemblages managed by the NPS. In spring 2009, Wind Cave National Park (WICA) served as a case study in a workshop assessing the use of scenario planning as a tool for park management planning in the face of rapidly changing climate. One outcome of the workshop was the recognized need for quantitative models to better understand the range of possible vegetation changes under different future climates and management decisions. This report addresses this need; it describes our adaptation of a dynamic global vegetation model (DGVM) to WICA vegetation and the resulting projections of future vegetation under three future climate scenarios and 11 management scenarios determined by Park natural resource managers. Wind Cave National Park lies along a narrow transition zone between the ponderosa pine (Pinus ponderosa) forests of the Black Hills and the mixed grass prairie that once extended with few interruptions over the lower, gentler terrain, subject to warmer, drier climate to the east and south of the Park. The location and character of this transition is strongly influenced by fire frequency and intensity (Brown and Sieg 1999). Furthermore, the mixed grass prairie occupies a broader transition zone between eastern tallgrass prairie and the shortgrass prairie of the western Great Plains. The dominance of species characteristic of these two prairie types varies with soil moisture availability, evaporative demand, and recent grazing history (Cogan et al. 1999). In addition, Wind Cave lies near the midpoint of a long gradient of C 3 (cool season) grass dominance to the north and C 4 (warm season) grass dominance to the south. The ecotonal position of WICA may make it particularly sensitive to climate change. For example, small changes in fire frequency and/or intensity and the vigor of trees vs. grass could dramatically shift the proportions of these two lifeforms. The Park hydrology is also sensitive to changes in the balance between infiltration of precipitation and evapotranspiration, as on average, only a small fraction of annual precipitation reaches the deeper soil layers that feed permanent streamflow. The resources at risk at Wind Cave NP include the Cave itself, as well as small backcountry caves, a genetically important bison herd, and other prairie species including the black-tailed prairie dog and endangered black-footed ferrets. All of these resources will be directly affected by climate change impacts on vegetation and hydrology. Natural resource management challenges at WICA are substantial, diverse, and intertwined. Aboveground, the park has been recognized as exemplary for its high quality vegetation (Marriot et al. 1999), though the park is relatively small for the diversity of vegetation types and species that it supports. Even without a changing climate, maintaining the integrity of the plant communities is complicated by the park’s legislated responsibility to maintain viable populations of bison, elk and pronghorn. In addition, the federally endangered black-footed ferret was recently re-introduced to the park. This species requires large extents of prairie dog towns for prey and habitat. Prairie dogs impact vegetation by constant clipping, grazing and soil disturbance, thereby affecting plant composition and productivity. Moreover, naturally high interannual climate variability and the strong influence of precipitation on grass productivity in this region combine to yield high interannual variability in the amount of forage available for the wildlife that the park is tasked to maintain. Finally, fire, which is now primarily controlled by WICA and NPS Northern Great Plains fire management programs, is intertwined with all other natural resource issues at WICA, as it can impact prairie dog colony and forest expansion, ungulate foraging behavior, invasive plant species, and hydrological processes. Although not capable of capturing all of these complexities, dynamic vegetation models do provide a means for quantitatively projecting vegetation futures in future climates under plausible fire and grazing regimes. Our work uses the DGVM MC1 to simulate the effects of future climate projections and management practices on the vegetation of WICA. MC1 is designed to project potential vegetation as influenced by natural processes and hence is appropriate for national parks, where conservation of native biota and ecosystems is of great importance. Since the initial application of MC1 to a small portion of WICA (Bachelet et al. 2000), the model has been altered to improve model performance with the inclusion of dynamic fire. Applying this improved version to WICA required substantial recalibration, during which we have made a number of improvements to MC1 that will be incorporated as permanent changes. In this report we document these changes and our calibration procedure following a brief overview of the model. We compare the projections of current vegetation to the current state of the park and present projections of vegetation dynamics under future climates downscaled from three GCMs selected to represent the existing range in available GCM projections. In doing so, we examine the consequences of different management options regarding fire and grazing, major aspects of biotic management at Wind Cave.

South Dakota↗

Functional response models to estimate feeding rates of wading birds

Forager (predator) abundance may mediate feeding rates in wading birds. Yet, when modeled, feeding rates are typically derived from the purely prey-dependent Holling Type II (HoII) functional response model. Estimates of feeding rates are necessary to evaluate wading bird foraging strategies and their role in food webs; thus, models that incorporate predator dependence warrant consideration. Here, data collected in a mangrove swamp in Puerto Rico in 1994 were reanalyzed, reporting feeding rates for mixed-species flocks after comparing fits of the HoII model, as used in the original work, to the Beddington-DeAngelis (BD) and Crowley-Martin (CM) predator-dependent models. Model CM received most support (AIC c wi = 0.44), but models BD and HoII were plausible alternatives (AIC c ??? 2). Results suggested that feeding rates were constrained by predator abundance. Reductions in rates were attributed to interference, which was consistent with the independently observed increase in aggression as flock size increased (P < 0.05). Substantial discrepancies between the CM and HoII models were possible depending on flock sizes used to model feeding rates. However, inferences derived from the HoII model, as used in the original work, were sound. While Holling's Type II and other purely prey-dependent models have fostered advances in wading bird foraging ecology, evaluating models that incorporate predator dependence could lead to a more adequate description of data and processes of interest. The mechanistic bases used to derive models used here lead to biologically interpretable results and advance understanding of wading bird foraging ecology.

Waterbirds↗