Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Nature Methods”

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 1,207 records · Page 67Linked to original sources

Using carbon dioxide in fisheries and aquatic invasive species management

To restore native fish populations, fisheries programs often depend on active removal of aquatic invasive species. Chemical removal can be an effective method of eliminating aquatic invasive species, but chemicals can induce mortality in nontarget organisms and persist in the environment. Carbon dioxide (CO 2 ) is an emerging alternative to traditional chemical control agents because it has been demonstrated to be toxic to fish, but is naturally occurring and readily neutralized. In addition, CO 2 is a commercially available gas, is highly soluble, and has high absorption efficiency. When these characteristics are paired with advances in modern, large-scale gas delivery technologies, opportunities to use CO 2 in natural or artificial (e.g., canals) waters to manage fish become increasingly feasible. Our objective is to describe the history of CO 2 use in fisheries and outline potential future applications of CO 2 to suppress and manipulate aquatic species in field and aquaculture settings.

Fisheries↗

Geologic Settings of Subsidence

This paper reviews the role of geologic processes that contribute to subsidence in order to aid those starting investigations of ground-surface subsidence. Subsidence occurs, or at least is discovered, only infrequently, and little organized information has been available. In order to assess our present state of knowledge, the author gathered fragmentary bits of information from many sources widely scattered in the literature of geology and other earth sciences. The author cites examples for each geologic process that is a potential contributor to ground-surface subsidence together with geologic evidence for diagnosing their causes. The geologic processes involved are: (1) solution of gypsum and salt and redistribution of transient fill materials through solution cavities in calcareous rocks; (2) underground erosion of uncemented or lightly cemented silt and sand through temporary underground passageways; (3) lateral plastic flow of salt, gypsum and anhydrite, shale, and clay under loading; (4) compaction of sediments by loading, drainage, vibration, and hydrocompaction; (5) tectonic movements including primary and secondary effects of earthquakes, folding, and warping; and (6) volcanic activity. Because the first four processes may be accelerated by various engineering activities, examples have been selected to illustrate subsidence both under natural conditions and under conditions modified by man's activities. Although this extensive search for existing information on the role of geologic factors in subsidence indicates that much detailed work remains to be done, the future prospects for advancing our geologic knowledge are excellent. Methods of measuring ground-surface displacements are improving rapidly. Also, broadly based investigations of known areas of major subsidence throughout the world are developing new methods of diagnosis and treatment and are yielding quantitative data that will aid our evaluation of the rates and magnitude of present-day geologic processes.

Book chapter↗

On the choice of statistical models for estimating occurrence and extinction from animal surveys

In surveys of natural animal populations the number of animals that are present and available to be detected at a sample location is often low, resulting in few or no detections. Low detection frequencies are especially common in surveys of imperiled species; however, the choice of sampling method and protocol also may influence the size of the population that is vulnerable to detection. In these circumstances, probabilities of animal occurrence and extinction will generally be estimated more accurately if the models used in data analysis account for differences in abundance among sample locations and for the dependence between site-specific abundance and detection. Simulation experiments are used to illustrate conditions wherein these types of models can be expected to outperform alternative estimators of population site occupancy and extinction. ?? 2007 by the Ecological Society of America.

Ecology↗

Robust estimation of population size when capture probabilities vary among animals

A model is given for multiple recapture studies on closed populations which allows capture probabilities to vary among individuals. The capture probability of each individual is assumed to be constant over time. Based on this model we give a nonparametric estimation procedure for population size. The estimator involves selecting one of a sequence of estimator which are each linear combinations of the capture frequencies. The individual estimators are derived from the generalized jackknife method. We also give a goodness of fit test for the model's assumption that individual capture probabilities do not change during the study. The robustness of the estimation procedure is investigated with a simulation study. By virtue of this study, and the theoretical nature of the estimator, it is judged to be robust to moderate variations in individual capture probabilities which may occur in commonly used short-term livetrapping studies.

Ecology↗

Water-resources reconnaissance of the Kwiguk (Emmonak) area, Alaska

As part of an agreement between the Alaska Department of Natural Resources and the United States Geological Survey, Water Resources Division, a reconnaissance visit was made of the area in the vicinity of the village of Kwiguk (Emmonak), Alaska, during the period June 17-19, 1970, to evaluate the water resources and to suggest methods of developing additional sources of water.

Alaska↗

Preliminary digital map of cryptocrystalline occurrences in northern Nevada

The purpose was to identify potential cryptocrystalline material sources for tools used by indigenous people of the northern Nevada portion of the Great Basin. Cryptocrystalline occurrence data combed from the U.S. Geological Survey's Mineral Resources Data System (MRDS, 1995) were combined with sites described in Nevada rockhound guides and entered into a geographic information system (GIS). The map area encompasses northern Nevada (fig.1). This open-file report describes the methods used to convert cryptocrystalline occurrence data into a digital format, documents the file structures, and explains how to download the digital files from the U.S. Geological Survey's World Wide Web site. Uses of the spatial dataset include, but are not limited to, natural and cultural resource management, interdisciplinary activities, recreational rockhounding, and gold exploration. It is important to note that the accuracy of the spatial data varies widely, and for some purposes, field checks are advised.

Nevada↗

Integrating habitat suitability and climate constraints to predict nonnative fish invasion risk for U.S. streams and inland lakes

Introduction: Nonnative fishes in streams and inland lakes can substantially alter aquatic ecosystems by degrading habitat conditions, competing with native species, and restructuring food webs. Identifying invasion hotspots can inform proactive management and prevention of nonnative species establishment. Methods: In this study, we identified streams and lakes across the conterminous United States that are vulnerable to invasion by 18 nonnative stream and 28 nonnative lake fish species. Habitat suitability models were developed separately for stream and lake species across nine ecoregions using natural and anthropogenic predictors, and species occurrences were predicted by integrating model results with climate match scores based on comparisons between environmental conditions in species’ native and potential nonnative ranges. We classified areas with both high habitat suitability and strong climate similarity as having high invasion potential and generated hotspot maps by stacking species-level predictions. Results: Florida exhibits the highest invasion risk for both habitat types. Stream hotspots were concentrated near major metropolitan areas, whereas lake hotspots were more widely distributed across the Great Lakes basin and the northeastern U.S. Twelve species were predicted to invade both habitat types, with goldfish ( Carassius auratus ) and pirapitinga ( Piaractus brachypomus ) showing particularly high invasion potential. Discussion: These findings provide a comprehensive assessment of nonnative freshwater fish invasion risk and support targeted monitoring and management strategies.

conterminous United States↗

Contributions of Arctic PRISM to monitoring western hemispheric shorebirds

Long-term monitoring of populations is of paramount importance to understanding responses oforganisms to global environmental change and to evaluating whether conservation practices are yielding intended results through time (Wiens 2009). The population status of many shorebird species, the focus of this volume, remain poorly known. Long-distance migrant shorebirds have proven particularly difficult to monitor, in part because of their highly migratory nature and ranges that extend into highly inaccessible regions. As migrant shorebirds travel the length of the hemisphere, they congregate and disperse in ways that vary among species, locations, and years, presenting serious challenges to designing and implementing monitoring programs. Rigorous field and quantitative methods that estimate population size and monitor trends are vitally needed to direct and evaluate effective conservation measures. Many management efforts depend on unbiased population size estimates; for example, the shorebird conservation plans for both Canada and the United States seek to restore populations to levels calculated for the 1970s based on the best information available from existing surveys. Further, federal wildlife agencies within the United States and Canada have mandates to understand the state of their nations' resources under various conventions for the protection of migratory birds. Accurate estimates of population size are vital statistics for a variety of conservation activities, such as prioritizing species for conservation action and setting management targets. Areas of essential habitat, such as those designated under the Western Hemisphere Shorebird Reserve Network, the Important Bird Areas program of BirdLife International and the National Audubon Society, or Canada's National Wildlife Areas program, are all evaluated on the basis ofproportions of species' populations which they contain. The size, and trends in size, ofa species' population are considered key information for assessing its vulnerability and subsequent listing under the U.S. Endangered Species Act and the Canadian Species at Risk Act. To meet the need for information on population size and trends, shorebird biologists from Canada and the United States proposed a shared blueprint for shorebird monitoring across the Western Hemisphere in the late 1990s; this effort was undertaken in concert with the development of the Canadian and U.S. Shorebird Conservation Plans (Donaldson et al. 2000, Brown et aL 2001). Soon thereafter, partners in the monitoring effort adopted the name "Program for Regional and International Shorebird Monitoring" (PRISM). Among the primary objectives of PRISM were to estimate the population sizes and trends of breeding North American shorebirds and describe their distributions (Bart et al. 2002). PRISM members evaluated ongoing and potential monitoring approaches to address 74 taxa (including subspecies) and proposed a combination of arctic andboreal breeding surveys, temperate breeding and non-breeding surveys, and neotropical surveys.

Book chapter↗

Research priorities for migratory birds under climate change—A qualitative value of information assessment

The mission of the U.S. Geological Survey National Climate Adaptation Science Center is to provide actionable, management-relevant research on climate change effects on ecosystems and wildlife to U.S. Department of the Interior bureaus. Providing this kind of useful scientific information requires understanding how natural-resource managers make decisions and identifying research priorities that support those decision-making processes. Migratory bird management and conservation of migratory bird habitat are central components of the U.S. Department of the Interior’s mission. In particular, the U.S. Fish and Wildlife Service has an intensive, complex decision-making process for identifying high-priority parcels of land that will contribute to migratory bird conservation through permanent acquisition or easement. Climate change introduces several uncertainties into this decision-making process, and additional climate change research should help to support more informed decision making regarding habitat acquisition. Not all climate change related uncertainties, however, will have a meaningful effect on acquisition decisions; therefore, understanding which uncertainties have the most potential to alter decision making is crucial. This document summarizes a multiyear effort to clarify the major sources of climate change uncertainty that affect migratory bird management and to articulate related research priorities. We worked with U.S. Fish and Wildlife Service staff to assess the primary ways in which climate change is likely to affect migratory birds and their habitats; to clarify uncertainties surrounding these effects; and to assess how uncertainties may affect habitat acquisition decisions. Using a modified structured decision-making approach, we assessed a set of hypotheses about how climate change will affect migratory birds and their habitats. Then, we used a qualitative value of information assessment to rank the most important topics for future research. The ranking process was built on an assessment of three primary characteristics: the magnitude of uncertainty, the topic’s relevance to habitat acquisition decision making, and the feasibility of reducing the uncertainty. Based on the results of this process, high-priority topics for future research include the following: The effects of rising temperatures on spatial distributions of migratory birds during the breeding and nonbreeding seasons; Climate-driven changes to avian community composition through homogenization and loss of specialists; The effects of decreased precipitation on abundance in the breeding season; and The effects of rising temperatures on abundance in the nonbreeding season. In addition to describing high-priority research needs, this document provides a summary of the methodology used to identify, assess, and rank uncertainties. This method was developed for a climate change related topic where a full quantitative value of information approach may not be feasible. The results and methodology described here may be useful for U.S. Geological Survey and other science-funding agencies interested in improving the applicability of their research to natural-resource management decision making.

Circular↗

Method for estimating potential wetland extent by utilizing streamflow statistics and flood-inundation mapping techniques: Pilot study for land along the Wabash River near Terre Haute, Indiana

Potential wetland extents were estimated for a 14-mile reach of the Wabash River near Terre Haute, Indiana. This pilot study was completed by the U.S. Geological Survey in cooperation with the U.S. Department of Agriculture, Natural Resources Conservation Service (NRCS). The study showed that potential wetland extents can be estimated by analyzing streamflow statistics with the available streamgage data, calculating the approximate water-surface elevation along the river, and generating maps by use of flood-inundation mapping techniques. Planning successful restorations for Wetland Reserve Program (WRP) easements requires a determination of areas that show evidence of being in a zone prone to sustained or frequent flooding. Zone determinations of this type are used by WRP planners to define the actively inundated area and make decisions on restoration-practice installation. According to WRP planning guidelines, a site needs to show evidence of being in an "inundation zone" that is prone to sustained or frequent flooding for a period of 7 consecutive days at least once every 2 years on average in order to meet the planning criteria for determining a wetland for a restoration in agricultural land. By calculating the annual highest 7-consecutive-day mean discharge with a 2-year recurrence interval (7MQ2) at a streamgage on the basis of available streamflow data, one can determine the water-surface elevation corresponding to the calculated flow that defines the estimated inundation zone along the river. By using the estimated water-surface elevation ("inundation elevation") along the river, an approximate extent of potential wetland for a restoration in agricultural land can be mapped. As part of the pilot study, a set of maps representing the estimated potential wetland extents was generated in a geographic information system (GIS) application by combining (1) a digital water-surface plane representing the surface of inundation elevation that sloped in the downstream direction of flow and (2) land-surface elevation data. These map products from the pilot study will aid the NRCS and its partners with the onsite inundation-zone verification in agricultural land for a potential restoration and will assist in determining at what elevation to plant hardwood trees for increased survivability on ground above frequently flooded terraces.

Indiana↗

Flood-frequency estimation for very low annual exceedance probabilities using historical, paleoflood, and regional information with consideration of nonstationarity

Streamflow estimates for floods with an annual exceedance probability of 0.001 or lower are needed to accurately portray risks to critical infrastructure, such as nuclear powerplants and large dams. However, extrapolating flood-frequency curves developed from at-site systematic streamflow records to very low annual exceedance probabilities (less than 0.001) results in large uncertainties in the streamflow estimates. Traditionally, methods for statistically estimating flood frequency have relied on the systematic streamflow record, which provides a time series of annual maximum flood peaks, often including some historical peaks. However, most peak-flow records are less than 100 years, and uncertainties are large when trying to extrapolate magnitudes of very low annual exceedance probability events. Other data may be available that extend the record beyond the systematic dataset. Historical data are defined as data from outside the period of systematic records but within the period of human records. Examples of historical information include flood estimates from other agencies and newspaper accounts that can be translated to flood magnitude point estimates, interval estimates, or perception thresholds (such as a statement that an 1880 flood was the largest since 1869). Paleoflood data, which may also extend the dataset, include a broad range of information about flood occurrence or magnitude from sources like sediment deposits or tree rings. Several assumptions are made in flood-frequency analysis, and an understanding of whether the data conform to these assumptions is desired. A particularly difficult assumption to evaluate for flood-frequency analysis is the underlying assumption that the flood series is stationary—the assumption that a time series of peak flow varies around a constant mean within a particular range of values (constant variance). As the hydrologic community’s understanding of natural systems and anthropogenic effects on streamflows has evolved, the community has come to understand that many surface-water systems exhibit one or more forms of nonstationarity, and thus the stationarity assumption is often violated to some degree. However, there is currently (2020) no consensus among hydrologists regarding the most appropriate flood-frequency-analysis methods for nonstationary systems, and this topic remains an active area of research. A literature review was completed to summarize the state of the science of flood frequency. The literature review highlights tools available to detect nonstationarities and identifies approaches that include external information to inform flood-frequency analysis. To demonstrate methods for initial data analysis and for incorporating historical and paleoflood information in flood-frequency analysis, five sites were selected: the Red River of the North at James Avenue Pumping Station, Winnipeg, Manitoba, Canada; lower reach, Rapid Creek, South Dakota; Spring Creek, South Dakota; Cherry Creek near Melvin, Colorado; and Escalante River near Escalante, Utah. The sites were chosen for the availability of published historical and paleoflood data and for their geographic diversity and unique characteristics, which highlighted issues such as autocorrelation, change points, trends, outlier peaks, or short periods of record. An initial data analysis that involved examining records for autocorrelation, change points, and trends was completed for all sites. The flood-frequency analysis completed for this study used version 7.2 of the U.S. Geological Survey PeakFQ program. Multiple analyses were done on each site documenting the change in the flood-frequency curve when additional historical or paleoflood data were added. When other flood-frequency studies were available, their results were compared to the results here. The comparisons in some cases simply show the effect of additional years of data, whereas other comparisons show results from probability distributions or fitting methods other than those used in PeakFQ. For the Red River of the North, flood-frequency analysis shows that paleoflood data appear necessary to reasonably estimate very low annual exceedance probabilities. For the analysis of the lower reach of Rapid Creek and Spring Creek, paleoflood information helped put a high outlier from the systematic period in context; however, very low annual exceedance probabilities at these sites still had extraordinarily large confidence bounds. These sites also showed that paleoflood information might be transferred from one site to another, with the caveat that this is a case where we had existing paleoflood data to test the transfer of paleoflood information—this is not the case at many sites, and transferring paleoflood information requires assumptions about the comparability of floods at the sites. The Cherry Creek analysis affirmed the result of an earlier study that showed that the generalized Pareto distribution was not a good distribution for estimating very low annual exceedance probabilities. The Escalante River analysis showed that adding paleoflood information might increase uncertainty for very low annual exceedance probabilities, compared to analysis with the systematic period of record information only, when the paleoflood peaks are of much larger magnitudes than the systematic record.

Scientific Investigations Report↗

Quality and quantity of terrestrial landscape connectivity in the US Great Lakes Coastal Zone

Context The Great Lakes Coastal Zone (GLCZ), the world’s longest continuous freshwater coastline, consists of interspersed natural and developed landcover and is inhabited by millions of people. Anthropogenic change fragments the GLCZ, decreasing landscape connectivity. Objectives We evaluated functional landscape connectivity of the U.S. GLCZ, assessing regional connectivity, protected areas’ effects on connectivity, and identified pinchpoints where connectivity is most vulnerable. Methods We modeled landscape connectivity to describe corridors and redundant pathways for two scenarios: (1) continuous terrestrial habitat patches and (2) protected areas. We compared these two scenarios by total corridor area and quality at two spatial scales (coastal zone and ecoregion). We used random forest analyses to examine how ecological factors influenced the connectivity rank and pinchpoint values. A principal components analysis identified how land cover influenced connectivity among ecoregions. Results The GLCZ was 69% connected. Although 91% of the GLCZ’s high-quality habitat was within corridors, less than 50% of the corridor area was high-quality (low resistance to movement). While the northernmost region retained high-quality connectivity, connectivity loss and degradation increased with development further south. Corridors between protected areas provided less connectivity, of worse quality, in areas of development and agriculture. Lower connectivity was associated with higher impervious surface cover. Conclusions Protected area placement affected both the quantity and quality of connectivity. Improving and maintaining functional connectivity could be accomplished by strategically restoring and protecting habitats most vulnerable to connectivity loss, predominantly in urban and agricultural regions. Our connectivity improvement index is a useful tool to identify areas where connectivity may be improved.

Illinois, Indiana, Michigan, Minnesota, New York, ↗

Elevation-based probabilistic mapping of irregularly flooded wetlands along the northern Gulf of Mexico coast

Irregularly flooded wetlands are found above the mean high water tidal datum and are exposed to tides and saltwater less frequently than daily. These wetlands provide important ecosystem services, such as providing habitat for fish and wildlife, enhancing water quality, ameliorating flooding impacts, supporting coastal food webs, and protecting upslope areas from erosion. Mapping irregularly flooded wetlands is challenging given their expansive coverage and dynamic nature. Furthermore, coastal wetlands are expected to change over the coming century due to sea-level rise and changes in the frequency and intensity of extreme storms. Consequently, coastal managers need baseline information on the spatial distribution of wetlands along with efficient and repeatable methods for observing changes. In this study, we used coastal wetlands from existing land use land cover data, best available lidar-derived digital elevation models, and Monte Carlo simulations to incorporate elevation uncertainty to create a probabilistic map of irregularly flooded wetlands along the northern Gulf of Mexico coast (USA). Our approach integrated findings from a review of coastal wetland elevation error in lidar datasets and an analysis of spatial autocorrelations of wetland elevation. We found a positive correlation ( r = 0.563, p < 0.0001) when comparing the probability estimated from a digital elevation model and in situ elevation observations. The differences in probability had a mean bias error of −0.04 (i.e., digital elevation model-based probability tends to be slightly lower), a mean absolute error of 0.20, and a root mean square error of 0.26. Beyond this overall validation, we explored error metrics for land cover classes and lidar collection details. To quantify areal coverage of the probabilistic output, we classified the probability values into equal bins using an interval of 0.33. The areal coverage of the lowest probability bin (“unlikely”; probability ≤0.33) was separated into the upper and lower portions of the irregularly flooded wetland zone. Of the coastal wetlands along the northern Gulf of Mexico coast about 38% were classified as unlikely and low with the greatest coverage in south Louisiana and the Everglades and around 33% were classified as unlikely and high with the greatest coverage in the Everglades and Texas. The relative coverage within the highest probability bin (“likely”; probability >0.66) covered around 13%, with the greatest coverage in south Florida, south Louisiana, and Texas. The framework developed in this study can be transferred to other coastal wetland areas and updated to observe changes with sea-level rise.

Gulf of Mexico coastline↗

Field and laboratory evaluations of visible light as a cue for guiding downstream-migrating juvenile Sea Lamprey

Objective We evaluated white light as a potential guidance cue for juvenile Sea Lamprey Petromyzon marinus in a natural setting as well as the effect of water velocity (0.25-, 0.50-, 0.75-, and 1.0-m/s test velocities) on light guidance behavior in a controlled laboratory flume, and characterized emigration timing and movement rates in a small stream (∼10 m wide and 0.7 m deep). Methods Behaviors and rates of downstream movement were monitored using PIT telemetry in both studies. Results In the field study, downstream movement by juveniles released during October 30–November 27 appeared to be cued by precipitation-induced flow events when water temperatures ranged between 4°C and 8°C. Juveniles expressed lateral attraction to a short, bank-mounted linear light array, but the guidance effect was not strong or consistent between bank light locations. Downstream movement rates decreased slightly when juveniles were exposed to the light cue. In the laboratory flume experiment, at water velocities of 0.25 and 0.75 m/s, lamprey were 2.8 and 3.3 times more likely to be detected at antennas along a wall with a linear light array compared with other antennas across the width of the flume. Significant changes in distribution were detected farther upstream in the flume during 0.25- and 0.50-m/s water velocity trials compared with 0.75-m/s trials. Further, the rate of downstream movement through the length of the flume decreased under artificial lighting compared with dark controls under the 0.25- and 0.75-m/s velocity conditions. Conclusions The results suggest lamprey exhibit a behavioral response to the light cue in both lab and field, but water velocity influences how effectively juveniles can respond to light cues.

Massachusetts↗

Assessing the risk persistent drought using climate model simulations and paleoclimate data

Projected changes in global rainfall patterns will likely alter water supplies and ecosystems in semiarid regions during the coming century. Instrumental and paleoclimate data indicate that natural hydroclimate fluctuations tend to be more energetic at low (multidecadal to multicentury) than at high (interannual) frequencies. State-of-the-art global climate models do not capture this characteristic of hydroclimate variability, suggesting that the models underestimate the risk of future persistent droughts. Methods are developed here for assessing the risk of such events in the coming century using climate model projections as well as observational (paleoclimate) information. Where instrumental and paleoclimate data are reliable, these methods may provide a more complete view of prolonged drought risk. In the U.S. Southwest, for instance, state-of-the-art climate model projections suggest the risk of a decade-scale megadrought in the coming century is less than 50%; the analysis herein suggests that the risk is at least 80%, and may be higher than 90% in certain areas. The likelihood of longer-lived events (>35 yr) is between 20% and 50%, and the risk of an unprecedented 50-yr megadrought is nonnegligible under the most severe warming scenario (5%&ndash;10%). These findings are important to consider as adaptation and mitigation strategies are developed to cope with regional impacts of climate change, where population growth is high and multidecadal megadrought&mdash;worse than anything seen during the last 2000 years&mdash;would pose unprecedented challenges to water resources in the region.

Journal of Climate↗

Sublake geologic structure from high-resolution seismic-reflection data from four sinkhole lakes in the Lake Wales Ridge, central Florida

Seismic-reflection profiles from Lake Wales, Blue Lake, Lake Letta, and Lake Apthorp located along the Lake Wales Ridge in central Florida provide local detail within the regional hydrogeologic framework as described by litho- and hydrostratigraphic cross sections. Lakes located with the mantled karst region have long been considered to be sinkhole lakes, originating from subsidence activity. High-resolution seismic- reflection data confirm this origin for these four lakes. The geologic framework of the Lake Wales Ridge has proven to be a suitable geologic setting for continuous high-resolution seismic-reflection profiling in lakes; however, the nature of the lake-bottom sediments largely controls the quality of the seismic data. In lakes with significant organic-rich bottom deposits, interpretable record was limited to areas where organic deposits were minimal. In lakes with clean, sandy bottoms, the seismic-reflection methods were highly successful in obtaining data that can be correlated with sublake subsidence features. These techniques are useful in examining sublake geology and providing a better understanding of how confining units are affected by subsidence in a region where their continuity is of significant importance to local lake hydrology. Although local geologic control around each lake generally corresponds to the regional geologic framework, local deviations from regional geologic trends occur in sublake areas affected by subsidence activity. Each of the four lakes examined represents a unique set of geologic controls and provides some degree of structural evidence of subsidence activity. Sublake geologic structures identified include: (1) marginal lake sediments dipping into bathymetric lows, (2) lateral discontinuity of confining units including sags and breaches, (3) the disruption and reworking of overlying unconsolidated siliciclastic sediments as they subside into the underlying irregular limestone surface, and (4) sublake regions where confining units appear to remain intact and unaffected by nearby subsidence activity. Each lake likely is underlain by several piping features rather than one large subsidence feature.

Florida↗

Documentation and verification of VST2D; a model for simulating transient, Variably Saturated, coupled water-heat-solute Transport in heterogeneous, anisotropic 2-Dimensional, ground-water systems with variable fluid density

This report describes a model for simulating transient, Variably Saturated, coupled water-heatsolute Transport in heterogeneous, anisotropic, 2-Dimensional, ground-water systems with variable fluid density (VST2D). VST2D was developed to help understand the effects of natural and anthropogenic factors on quantity and quality of variably saturated ground-water systems. The model solves simultaneously for one or more dependent variables (pressure, temperature, and concentration) at nodes in a horizontal or vertical mesh using a quasi-linearized general minimum residual method. This approach enhances computational speed beyond the speed of a sequential approach. Heterogeneous and anisotropic conditions are implemented locally using individual element property descriptions. This implementation allows local principal directions to differ among elements and from the global solution domain coordinates. Boundary conditions can include time-varying pressure head (or moisture content), heat, and/or concentration; fluxes distributed along domain boundaries and/or at internal node points; and/or convective moisture, heat, and solute fluxes along the domain boundaries; and/or unit hydraulic gradient along domain boundaries. Other model features include temperature and concentration dependent density (liquid and vapor) and viscosity, sorption and/or decay of a solute, and capability to determine moisture content beyond residual to zero. These features are described in the documentation together with development of the governing equations, application of the finite-element formulation (using the Galerkin approach), solution procedure, mass and energy balance considerations, input requirements, and output options. The VST2D model was verified, and results included solutions for problems of water transport under isohaline and isothermal conditions, heat transport under isobaric and isohaline conditions, solute transport under isobaric and isothermal conditions, and coupled water-heat-solute transport. The first three problems considered in model verification were compared to either analytical or numerical solutions, whereas the coupled problem was compared to measured laboratory results for which no known analytic solutions or numerical models are available. The test results indicate the model is accurate and applicable for a wide range of conditions, including when water (liquid and vapor), heat (sensible and latent), and solute are coupled in ground-water systems. The cumulative residual errors for the coupled problem tested was less than 10–8 cubic centimeter per cubic centimeter, 10-5 moles per kilogram, and 102 calories per cubic meter for liquid water content, solute concentration and heat content, respectively. This model should be useful to hydrologists, engineers, and researchers interested in studying coupled processes associated with variably saturated transport in ground-water systems.

Water-Resources Investigations Report↗

Methods for estimating low-flow statistics for Massachusetts streams

Methods and computer software are described in this report for determining flow duration, low-flow frequency statistics, and August median flows. These low-flow statistics can be estimated for unregulated streams in Massachusetts using different methods depending on whether the location of interest is at a streamgaging station, a low-flow partial-record station, or an ungaged site where no data are available. Low-flow statistics for streamgaging stations can be estimated using standard U.S. Geological Survey methods described in the report. The MOVE.1 mathematical method and a graphical correlation method can be used to estimate low-flow statistics for low-flow partial-record stations. The MOVE.1 method is recommended when the relation between measured flows at a partial-record station and daily mean flows at a nearby, hydrologically similar streamgaging station is linear, and the graphical method is recommended when the relation is curved. Equations are presented for computing the variance and equivalent years of record for estimates of low-flow statistics for low-flow partial-record stations when either a single or multiple index stations are used to determine the estimates. The drainage-area ratio method or regression equations can be used to estimate low-flow statistics for ungaged sites where no data are available. The drainage-area ratio method is generally as accurate as or more accurate than regression estimates when the drainage-area ratio for an ungaged site is between 0.3 and 1.5 times the drainage area of the index data-collection site. Regression equations were developed to estimate the natural, long-term 99-, 98-, 95-, 90-, 85-, 80-, 75-, 70-, 60-, and 50-percent duration flows; the 7-day, 2-year and the 7-day, 10-year low flows; and the August median flow for ungaged sites in Massachusetts. Streamflow statistics and basin characteristics for 87 to 133 streamgaging stations and low-flow partial-record stations were used to develop the equations. The streamgaging stations had from 2 to 81 years of record, with a mean record length of 37 years. The low-flow partial-record stations had from 8 to 36 streamflow measurements, with a median of 14 measurements. All basin characteristics were determined from digital map data. The basin characteristics that were statistically significant in most of the final regression equations were drainage area, the area of stratified-drift deposits per unit of stream length plus 0.1, mean basin slope, and an indicator variable that was 0 in the eastern region and 1 in the western region of Massachusetts. The equations were developed by use of weighted-least-squares regression analyses, with weights assigned proportional to the years of record and inversely proportional to the variances of the streamflow statistics for the stations. Standard errors of prediction ranged from 70.7 to 17.5 percent for the equations to predict the 7-day, 10-year low flow and 50-percent duration flow, respectively. The equations are not applicable for use in the Southeast Coastal region of the State, or where basin characteristics for the selected ungaged site are outside the ranges of those for the stations used in the regression analyses. A World Wide Web application was developed that provides streamflow statistics for data collection stations from a data base and for ungaged sites by measuring the necessary basin characteristics for the site and solving the regression equations. Output provided by the Web application for ungaged sites includes a map of the drainage-basin boundary determined for the site, the measured basin characteristics, the estimated streamflow statistics, and 90-percent prediction intervals for the estimates. An equation is provided for combining regression and correlation estimates to obtain improved estimates of the streamflow statistics for low-flow partial-record stations. An equation is also provided for combining regression and drainage-area ratio estimates to obtain improved e

Massachusetts↗