USGS Science⌕ Search

SEARCH · USGS Science

Results for “Computers and Geology”

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,063 records · Page 59Linked to original sources
USGS1973

Atlas of Mars, Amenthes quadrangle, MC-14

This mosaic is made with computer enhanced Mariner 9 pictures, rectified and scaled to a Mercator projection. Computer processing was performed by the Image Processing Laboratory (IPL) of the Jet Propulsion Laboratory (JPL). Processing included contrast enhancement and high-pass filtering which accentuates topographic features while subduing albedo variations. Photographic reproductions of the images were generated from the IPL digital magnetic tapes and mosaicked by the U.S. Geological Survey. Relative horizontal placement of pictures is controlled by spacecraft tracking data, with adjustment to fit the Mariner 9 control net of August, 1972 (Davies, M.E., 1972, Mariner 9 control net of Mars: August, 1972: Rand Corp., R-1122-JPL) where possible. Individual analytical adjustment of pictures was not performed. Triangles are shown at the correct latitude and longitude of the control points identified and numbered by Davies (ibid.). Broken crosses are used to identify the images of the control points. Standard deviation of image placement discrepancy is less than 5 km. The-projection is based on a theoretical dynamical ellipticity of 1/192. The scale is based on an equatorial radius of 3393.4 km.

Open-File Report↗
USGS1973

Atlas of Mars, Phoenicis Lacus quadrangle, MC-17

This mosaic is made with computer enhanced Mariner 9 pictures, rectified and scaled to a Mercator projection. Computer processing was performed by the Image Processing Laboratory (IPL) of the Jet Propulsion Laboratory (JPL). Processing included contrast enhancement and high-pass filtering which accentuates topographic features while subduing albedo variations. Photographic reproductions of the images were generated from the IPL digital magnetic tapes and mosaicked by the U.S. Geological Survey. Horizontal controls were derived from spacecraft tracking data, which are based on astronomical definitions of the spin axis and system of longitudes in use during the Mariner 9 mission. The projection is based on a theoretical dynamical ellipticity of 1/192. The scale is based on an equatorial radius of 3393.4 km. Standard deviation in image placement is less than 5 km at map scale.

Open-File Report↗
USGS1976

Field verification of method for distributing flow through multiple-bridge openings

Field data collected at three sites in Louisiana having multiple-bridge openings were used to check a method of distribution of flow through two or more bridges. Actual peak discharges in the various main and overflow bridges ranged from 878 to 11,055 cubic feet per second (24.9 to 313 m 2 /s). The error of the computed discharge to the actual discharge ranged between -20 and +6 percent, with the exception of one opening that had an error of +101 percent.

Louisiana↗
USGS2024

Challenges in developing vertical hazard for seismic analysis of concrete dams

The seismic response of concrete dams depends on both the horizontal (H) and vertical (V) components of ground motion (GM), and excluding the V components when conducting response history analyses (RHAs) may underestimate the seismic fragility. Although V components of GM time series could be selected to be consistent with the hazard, hazard curves for the V component (or for short, V hazard curves) are rarely available for evaluating such hazard consistency. One of the challenges in developing V hazard curves in the U.S. Geological Survey (USGS) National Seismic Hazard Model (NSHM) is the lack of ground motion models (GMMs) that directly estimate intensity measures (IMs) for the V component of GM. In this study, we compare alternative approaches of developing V hazard curves. Using a case study in California, we compare V hazard curves computed from GMMs for the V component of GM against those from applying GMMs for the V/H ratios and those from selecting multicomponent time series that are consistent with the H hazard. We discuss and invite feedback on issues such as availability of GMMs, applicability of GMMs, and choice of inputs for ground motion selection and modification (GMSM).

Conference Paper↗
USGS1996

Evaluation and modification of five techniques for estimating stormwater runoff for watersheds in west-central Florida

Several traditional techniques have been used for estimating stormwater runoff from ungaged watersheds. Applying these techniques to watersheds in west-central Florida requires that some of the empirical relations be extrapolated beyond tested ranges. As a result, there is uncertainty as to the accuracy of these estimates. Sixty-six storms occurring in 15 west-central Florida watersheds were initially modeled using the Rational Method, the U.S. Geological Survey Regional Regression Equations, the Natural Resources Conservation Service TR-20 model, the U.S. Army Corps of Engineers HEC-1 model, and the U.S. Environmental Protection Agency Storm Water Management Model. The techniques were applied according to the guidelines specified in the user manuals or standard engineering textbooks as though no field data were available and the selection of input parameters was not influenced by observed data. Computed estimates were compared with observed runoff to evaluate the accuracy of the techniques. One watershed was eliminated from further evaluation when it was determined that the area contributing runoff to the stream varies with the amount and intensity of rainfall. Therefore, further evaluation and modification of the input parameters were made for only 62 storms in 14 watersheds. Runoff ranged from 1.4 to 99.3 percent of rainfall. The average runoff for all watersheds included in this study was about 36 percent of rainfall. The average runoff for the urban, natural, and mixed land use watersheds was about 41, 27, and 29 percent of rainfall, respectively. Initial estimates of peak discharge using the Rational Method produced average watershed errors that ranged from an underestimation of 50.4 percent to an overestimation of 767 percent. The coefficient of runoff ranged from 0.20 to 0.60. Calibration of the technique produced average errors that ranged from an underestimation of 3.3 percent to an over estimation of 1.5 percent. The average calibrated coefficient of runoff for each watershed ranged from 0.02 to 0.72. The average values of the coefficient of runoff necessary to calibrate the urban, natural, and mixed land use watersheds were 0.39, 0.16, and 0.08, respectively. The U.S. Geological Survey regional regression equations for determining peak discharge produced errors that ranged from an underestimation of 87.3 percent to an overestimation of 1,140 percent. The regression equations for determining runoff volume produced errors that ranged from an underestimation of 95.6 percent to an overestimation of 324 percent. Regression equations developed from data used for this study produced errors that ranged between an underestimation of 82.8 percent and an overestimation of 328 percent for peak discharge and from an underestimation of 71.2 percent to an overestimation of 241 percent for runoff volume. Use of the equations developed for west-central Florida streams produced average errors for each type of watershed that were lower than errors associated with use of the U.S. Geological Survey regional equations. Initial estimates of peak discharges and runoff volumes using the Natural Resources Conservation Service TR-20 model, produced average errors of 44.6 and 42.7 percent, respectively, for all the watersheds. Curve numbers and times of concentration were adjusted to match estimated and observed peak discharges and runoff volumes. The average change in the curve number for all the watersheds was a decrease of 2.8 percent. The average change in the time of concentration was an increase of 59.2 percent. The shape of the input dimensionless unit hydrograph also had to be adjusted to match the shape and peak time of the estimated and observed flood hydrographs. Peak rate factors for the modified input dimensionless unit hydrographs ranged from 162 to 454. The mean errors for peak discharges and runoff volumes were reduced to 18.9 and 19.5 percent, respectively, using the average calibrated input parameters for each watershed. Initial estimates of peak discharges and runoff volumes using the U.S. Army Corp of Engineers Hydrologic Engineering Center-1 model, produced average errors of 105 and 26.8 percent respectively, for all the watersheds. Curve numbers and lag times were adjusted to match estimated and observed peak discharges and runoff volumes. The average change in the curve number for all the watersheds was a decrease of 2.5 percent. The average change in the lag time was an increase of 169 percent. The mean errors for peak discharges and runoff volumes were reduced to 5.8 and 1.4 percent, respectively, using the average calibrated input parameters for each watershed. The observed and estimated peak discharges and runoff volumes could be matched by adjusting curve numbers and lag time using the U.S. Army Corp of Engineers Hydrologic Engineering Center-1 model; however, the shape of the estimated flood hydrograph and timing of the peak could not be matched. The input dimensionless unit hydrograph must also be changed to increase the accuracy of the Hydrologic Engineering Center-1 model for watersheds in west-central Florida. The source code has to be modified and recompiled to enter different dimensionless unit hydrographs into the HEC-1 program. During application of the U.S. Environmental Protection Agency Storm Water Management Model, two separate infiltration methods were evaluated. Initial estimates of peak discharges and runoff volumes produced mean errors of 46.5 and 6.8 percent, respectively, for all watersheds using the Green-Ampt infiltration method, and 48.8 and 9.5 percent, respectively, using the Horton infiltration method. The mean errors were reduced to 18 and 0.3 percent for the Green-Ampt method and 20.9 and 7.2 percent for the Horton method using the average calibrated input parameters for each watershed. Estimates of peak discharges and runoff volumes were initially made for watersheds in west-central Florida using recommended procedures, then compared to observed peak discharges and runoff volumes. Subsequently, the procedures were modified to increase accuracy for this area. The same methods used during the study could be used in other parts of the world to evaluate the accuracy of standard methods for estimating stormwater runoff.

Florida↗
USGS2019

Xenon hydrate as an analogue of methane hydrate in geologic systems out of thermodynamic equilibrium

Methane hydrate occurs naturally under pressure and temperature conditions that are not straightforward to replicate experimentally. Xenon has emerged as an attractive laboratory alternative to methane for studying hydrate formation and dissociation in multiphase systems, given that it forms hydrates under milder conditions. However, building reliable analogies between the two hydrates requires systematic comparisons, which are currently lacking. We address this gap by developing a theoretical and computational model of gas hydrates under equilibrium and nonequilibrium conditions. We first compare equilibrium phase behaviors of the Xe·H 2 O and CH 4 ·H 2 O systems by calculating their isobaric phase diagram, and then study the nonequilibrium kinetics of interfacial hydrate growth using a phase field model. Our results show that Xe·H 2 O is a good experimental analog to CH 4 ·H 2 O, but there are key differences to consider. In particular, the aqueous solubility of xenon is altered by the presence of hydrate, similar to what is observed for methane; but xenon is consistently less soluble than methane. Xenon hydrate has a wider nonstoichiometry region, which could lead to a thicker hydrate layer at the gas‐liquid interface when grown under similar kinetic forcing conditions. For both systems, our numerical calculations reveal that hydrate nonstoichiometry coupled with hydrate formation dynamics leads to a compositional gradient across the hydrate layer, where the stoichiometric ratio increases from the gas‐facing side to the liquid‐facing side. Our analysis suggests that accurate composition measurements could be used to infer the kinetic history of hydrate formation in natural settings where gas is abundant.

Geochemistry, Geophysics, Geosystems↗
USGS2004

Development of a geodatabase and conceptual model of the hydrogeologic units beneath air force plant 4 and Naval Air Station-Joint Reserve Base Carswell Field, Fort Worth, Texas

Air Force Plant 4 and adjacent Naval Air Station-Joint Reserve Base Carswell Field at Fort Worth, Texas, constitute a government-owned, contractor-operated facility that has been in operation since 1942. Contaminants from AFP4, primarily volatile organic compounds and metals, have entered the ground-water-flow system through leakage from waste-disposal sites and from manufacturing processes. The U.S. Geological Survey developed a comprehensive geodatabase of temporal and spatial environmental information associated with the hydrogeologic units (alluvial aquifer, Goodland-Walnut confining unit, and Paluxy aquifer) beneath the facility and a three-dimensional conceptual model of the hydrogeologic units integrally linked to the geodatabase. The geodatabase design uses a thematic layer approach to create layers of feature data using a geographic information system. The various features are separated into relational tables in the geodatabase on the basis of how they interact and correspond to one another. Using the geodatabase, geographic data at the site are manipulated to produce maps, allow interactive queries, and perform spatial analyses. The conceptual model for the study area comprises computer-generated, three-dimensional block diagrams of the hydrogeologic units. The conceptual model provides a platform for visualization of hydrogeologic-unit sections and surfaces and for subsurface environmental analyses. The conceptual model is based on three structural surfaces and two thickness configurations of the study area. The three structural surfaces depict the altitudes of the tops of the three hydrogeologic units. The two thickness configurations are those of the alluvial aquifer and the Goodland-Walnut confining unit. The surface of the alluvial aquifer was created using a U.S. Geological Survey 10-meter digital elevation model. The 2,130 point altitudes of the top of the Goodland-Walnut unit were compiled from lithologic logs from existing wells, available soil-boring logs, and previous studies. Data from 120 wells, primarily from existing reports, were used to create a map of the approximate altitude of the Paluxy aquifer.

Texas↗
USGS2013

Investigations of groundwater system and simulation of regional groundwater flow for North Penn Area 7 Superfund site, Montgomery County, Pennsylvania

Groundwater in the vicinity of several industrial facilities in Upper Gwynedd Township and vicinity, Montgomery County, in southeast Pennsylvania has been shown to be contaminated with volatile organic compounds (VOCs), the most common of which is the solvent trichloroethylene (TCE). The 2-square-mile area was placed on the National Priorities List as the North Penn Area 7 Superfund site by the U.S. Environmental Protection Agency (USEPA) in 1989. The U.S. Geological Survey (USGS) conducted geophysical logging, aquifer testing, and water-level monitoring, and measured streamflows in and near North Penn Area 7 from fall 2000 through fall 2006 in a technical assistance study for the USEPA to develop an understanding of the hydrogeologic framework in the area as part of the USEPA Remedial Investigation. In addition, the USGS developed a groundwater-flow computer model based on the hydrogeologic framework to simulate regional groundwater flow and to estimate directions of groundwater flow and pathways of groundwater contaminants. The study area is underlain by Triassic- and Jurassic-age sandstones and shales of the Lockatong Formation and Brunswick Group in the Mesozoic Newark Basin. Regionally, these rocks strike northeast and dip to the northwest. The sequence of rocks form a fractured-sedimentary-rock aquifer that acts as a set of confined to partially confined layers of differing permeabilities. Depth to competent bedrock typically is less than 20 ft below land surface. The aquifer layers are recharged locally by precipitation and discharge locally to streams. The general configuration of the potentiometric surface in the aquifer is similar to topography, except in areas affected by pumping. The headwaters of Wissahickon Creek are nearby, and the stream flows southwest, parallel to strike, to bisect North Penn Area 7. Groundwater is pumped in the vicinity of North Penn Area 7 for industrial use, public supply, and residential supply. Results of field investigations by USGS at the site and results from other studies support, and are consistent with, a conceptual model of a layered leaky aquifer where the dip of the beds has a strong control on hydraulic connections in the groundwater system. Connections within and (or) parallel to bedding tend to be greater than across bedding. Transmissivities of aquifer intervals isolated by packers ranged over three orders of magnitude [from about 2.8 to 2,290 square feet per day (ft 2 /d) or 0.26 to 213 square meters per day (m 2 /d)], did not appear to differ much by mapped geologic unit, but showed some relation to depth being relatively smaller in the shallowest and deepest intervals (0 to 50 ft and more than 250 ft below land surface, respectively) compared to the intermediate depth intervals (50 to 250 ft below land surface) tested. Transmissivities estimated from multiple-observation well aquifer tests ranged from about 700 to 2,300 ft 2 /d (65 to 214 m 2 /d). Results of chemical analyses of water from isolated intervals or monitoring wells open to short sections of the aquifer show vertical differences in concentrations; chloride and silica concentrations generally were greater in shallow intervals than in deeper intervals. Chloride concentrations greater than 100 milligrams per liter (mg/L), combined with distinctive chloride/bromide ratios, indicate a different source of chloride in the western part of North Penn Area 7 than elsewhere in the site. Groundwater flow at a regional scale under steady-state conditions was simulated by use of a numerical model (MODFLOW-2000) for North Penn Area 7 with different layers representing saprolite/highly weathered rock near the surface and unweathered competent bedrock. The sedimentary formations that underlie the study area were modeled using dipping model layers for intermediate and deep zones of unweathered, fractured rock. Horizontal cell model size was 100 meters (m) by 100 meters (328 ft by 328 ft), and model layer thickness ranged from 6 m (19.7 ft) representing shallow weathered rock and saprolite up to 200 m (656 ft) representing deeper dipping bedrock. The model did not include detailed structure to account for local-scale differences in hydraulic properties, with the result that local-scale groundwater flow may not be well simulated. Additional detailed multi-well aquifer tests would be needed to establish the extent of interconnection between intervals at the local scale to address remediation of contamination at each source area. This regional groundwater-flow model was calibrated against measured groundwater levels (1996, 2000, and 2005) and base flow estimated from selected streamflow measurements by use of nonlinear-regression parameter-estimation algorithms to determine hydraulic conductivity and anisotropy of hydraulic conductivity, streambed hydraulic conductivity, and recharge during calibration periods. Results of the simulation using the calibrated regional model indicate that the aquifer appears to be anisotropic where hydraulic conductivity is greatest parallel to the orientation of bedding of the formations underlying the area and least in the cross-bed direction. The maximum hydraulic conductivity is aligned with the average regional strike of the formations, which is “subhorizontal” in the model because the altitudes of the beds and model cells vary in the strike, as well as dip, direction. Estimated subhorizontal hydraulic conductivities (in strike direction parallel to dipping beds) range from 0.001 to 1.67 meters per day (0.0032 to 5.5 feet per day). The ratio of minimum (dip direction) to maximum (strike direction) subhorizontal hydraulic conductivity ranges from 1/3.1 to 1/8.6, and the ratio of vertical to horizontal hydraulic conductivity ranges from 1/1 to 1/478. However, limited available field data precluded rigorous calibration of vertical anisotropy in the model. Estimated recharge rates corresponding to calibration periods in 1996, 2000, and 2005 are 150, 109, and 124 millimeters per year (5.9, 4.3, and 4.9 inches per year), respectively. The calibrated groundwater-flow model was used to simulate groundwater flow under steady-state conditions during periods of relatively high withdrawals (pumpage) (1990) and relatively low withdrawals (2000 and 2005). Groundwater-flow paths originating from recharge areas near known areas of soil contamination (sources) were simulated. Pumped industrial and production wells captured more groundwater from several of these sources during 1990 than after 1990 when pumping declined or ceased and greater amounts of contaminated groundwater moved away from North Penn Area 7 Superfund site to surrounding areas. Uncertainty in simulated groundwater-flow paths from contaminant sources and contributing areas, resulting from uncertainty in estimated hydraulic properties of the model, was illustrated through Monte Carlo simulations. The effect of uncertainty in the vertical anisotropy was not included in the Monte Carlo simulations. Contributing areas indicating the general configuration of groundwater flow towards production well MG-202 (L-22) in the study area also were simulated for the different time periods; as simulated, the flow paths do not pass through any identified contaminant source in North Penn Area 7. However, contributing areas to wells, such as MG-202, located near many pumped wells are particularly complex and, in some cases, include areas that contribute flow to streams that subsequently recharge the aquifer through stream loss. In these cases, water-quality constituents, including contaminants that are present in surface water may be drawn into the aquifer to nearby pumped wells. Results of a simulated shutdown of well MG-202 under steady-state 2005 conditions showed that the area contributing recharge for nearby production well MG-76 (L-17), when MG-202 is not pumping, shifts downstream and is similar to the area contributing recharge for MG-202 when both wells are pumping. Concentrations of constituents in groundwater samples collected in fall 2005 or spring 2006 were compared to simulated groundwater-flow paths for the year 2005 to provide a qualitative assessment of model results. The observed spatial distribution of selected constituents, including TCE, CFC-11, and CFC-113 in groundwater in 2005 and the chloride/bromide mass ratios in 2006, generally were consistent with the model results of the simulated 2005 groundwater-flow paths at North Penn Area 7, indicating the presence of several separate sources of contaminants within North Penn Area 7.

Pennsylvania↗
USGS1992

Sampling design for spatially distributed hydrogeologic and environmental processes

A methodology for the design of sampling networks over space is proposed. The methodology is based on spatial random field representations of nonhomogeneous natural processes, and on optimal spatial estimation techniques. One of the most important results of random field theory for physical sciences is its rationalization of correlations in spatial variability of natural processes. This correlation is extremely important both for interpreting spatially distributed observations and for predictive performance. The extent of site sampling and the types of data to be collected will depend on the relationship of subsurface variability to predictive uncertainty. While hypothesis formulation and initial identification of spatial variability characteristics are based on scientific understanding (such as knowledge of the physics of the underlying phenomena, geological interpretations, intuition and experience), the support offered by field data is statistically modelled. This model is not limited by the geometric nature of sampling and covers a wide range in subsurface uncertainties. A factorization scheme of the sampling error variance is derived, which possesses certain atttactive properties allowing significant savings in computations. By means of this scheme, a practical sampling design procedure providing suitable indices of the sampling error variance is established. These indices can be used by way of multiobjective decision criteria to obtain the best sampling strategy. Neither the actual implementation of the in-situ sampling nor the solution of the large spatial estimation systems of equations are necessary. The required values of the accuracy parameters involved in the network design are derived using reference charts (readily available for various combinations of data configurations and spatial variability parameters) and certain simple yet accurate analytical formulas. Insight is gained by applying the proposed sampling procedure to realistic examples related to sampling problems in two dimensions.

Advances in Water Resources↗
USGS2010

An efficient deterministic-probabilistic approach to modeling regional groundwater flow: 2. Application to Owens Valley, California

The applicability of a deterministic-probabilistic model for predicting water tables in southern Owens Valley, California, is evaluated. The model is based on a two-layer deterministic model that is cascaded with a two-point probability model. To reduce the potentially large number of uncertain variables in the deterministic model, lumping of uncertain variables was evaluated by sensitivity analysis to reduce the total number of uncertain variables to three variables: hydraulic conductivity, storage coefficient or specific yield, and source-sink function. Results demonstrate that lumping of uncertain parameters reduces computational effort while providing sufficient precision for the case studied. Simulated spatial coefficients of variation for water table temporal position in most of the basin is small, which suggests that deterministic models can predict water tables in these areas with good precision. However, in several important areas where pumping occurs or the geology is complex, the simulated spatial coefficients of variation are over estimated by the two-point probability method.

California↗
USGS2026

As above, so below? A framework for integrating long-term water quantity trends reveals divergent patterns in groundwater and low streamflow across the United States

Climate, land-use, and disturbance drive long-term global trends in groundwater levels and streamflow. At large scales, these trends are typically considered separately, despite the well-established concept that groundwater and surface water comprise a single resource. Joint trend assessment at national scales is challenging because it requires pairing and aggregating data from spatially disparate streamflow and groundwater monitoring sites for which no established framework exists. Here, we evaluate alternative approaches for integrating groundwater and streamflow data to enable joint trend analysis—a critical step toward understanding how water-budget components respond concurrently and interactively to environmental drivers. Mann–Kendall trends were computed for individual groundwater (annual mean depth) and streamflow (annual low of 7 d averages) sites across the U.S over 21- (2000–2020), 31- (1990–2020), and 41-year (1980–2020) periods. Regional Kendall trends were calculated using five spatially contiguous and noncontiguous regional classifications for aggregation based on subsurface (e.g. aquifer, geology) and surface (e.g. watershed, landscape) characteristics. Site-level results revealed contrasting trends, with tendencies toward increasing low flows (wetting) and increasing groundwater depths (drying). Agreement between streamflow and groundwater trends increased with regional aggregation and longer timeframes, though persistent skew toward streamflow wetting and groundwater drying remained. Results varied by region and trend period, with notable consistencies: unified drying in the West/Southwest and wetting in the Upper Midwest. Directional mismatches in long-term trends were prominent in the High Plains and Mississippi Alluvial Plain, whereas near-term mismatches were most evident in the Northwest. Aggregation by hydrologic landscape regions (HLR) yielded the greatest agreement between groundwater and streamflow trends. These findings indicate that coupled responses may represent combined influences of climate, relief, and geology, as captured by HLR, more strongly than geography or geology alone. Integrated water availability assessments may benefit from a multi-characteristic classification framework to treat groundwater and surface water as a unified resource.

Environmental Research: Water↗
USGS2012

Flood inundation map library, Fort Kent, Maine

Severe flooding occurred in northern Maine from April 28 to May 1, 2008, and damage was extensive in the town of Fort Kent (Lombard, 2010). Aroostook County was declared a Federal disaster area on May 9, 2008. The extent of flooding on both the Fish and St. John Rivers during this event showed that the current Federal Emergency Management Agency (FEMA) Flood Insurance Study (FIS) and Flood Insurance Rate Map (FIRM) (Federal Emergency Management Agency, 1979) were out of date. The U.S. Geological Survey (USGS) conducted a study to develop a flood inundation map library showing the areas and depths for a range of flood stages from bankfull to the flood of record for Fort Kent to complement an updated FIS (Federal Emergency Management Agency, in press). Hydrologic analyses that support the maps include computer models with and without the levee and with various depths of backwater on the Fish River. This fact sheet describes the methods used to develop the maps and describes how the maps can be accessed.

Maine↗
USGS2017

Hydrologic and hydraulic analyses of Great Meadow wetland, Acadia National Park, Maine

The U.S. Geological Survey completed hydrologic and hydraulic analyses of Cromwell Brook and the Sieur de Monts tributary in Acadia National Park, Maine, to better understand causes of flooding in complex hydrologic and hydraulic environments, like those in the Great Meadow wetland and Sieur de Monts Spring area. Regional regression equations were used to compute peak flows with from 2 to 100-year recurrence intervals at seven locations. Light detection and ranging data were adjusted for bias caused by dense vegetation in the Great Meadow wetland; and then combined with local ground surveys used to define the underwater topography and hydraulic structures in the study area. Hydraulic modeling was used to evaluate flood response in the study area to a variety of hydrologic and hydraulic scenarios. Hydraulic modeling indicates that enlarging the culvert at Park Loop Road could help mitigate flooding near the Sieur de Monts Nature Center that is caused by streamflows with large recurrence intervals; however, hydraulic modeling also indicates that the Park Loop Road culvert does not aggravate flooding near the Nature Center caused by the more frequent high intensity rainstorms. That flooding is likely associated with overland flow resulting from (1) quick runoff from the steep Dorr Mountain hitting the lower gradient Great Meadow wetland area and (2) poor drainage aggravated by beaver dams holding water in the wetland. Rapid geomorphic assessment data collected in June 2015 and again in April 2016 indicate that Cromwell Brook has evidence of aggradation, degradation, and channel widening throughout the drainage basin. Two of five reference cross sections developed for this report also indicate channel aggradation.

Maine↗
USGS2020

Bathymetry of Morris Lake (Newton Reservoir), New Jersey, 2018

Morris Lake, also known as Newton Reservoir, has been the source of drinking water for the Town of Newton, New Jersey, since the early 1900s. Although Morris Lake has been used as a source of drinking water for many years, its capacity was previously uncertain. In April 2018, the U.S. Geological Survey and the New Jersey Department of Environmental Protection conducted a bathymetric survey of Morris Lake using a multibeam echosounder to map the reservoir. The points measured with the multibeam echosounder were combined with light detection and ranging data above the water surface and processed to create a 3.3-foot (1 meter) raster grid of the bathymetric surface, bathymetric contours at 2-foot intervals of depth and elevation, and an elevation-area-capacity table. The results of the bathymetric survey show that Morris Lake has a maximum depth of just over 119 feet with an average depth of 42 feet. Like the surrounding topography, parts of the reservoir are extremely steep. The capacity of the reservoir at full spillway level is 1,980 million gallons, with a corresponding surface area of 145 acres. The accuracy of the mapped multibeam echosounder bathymetric data was evaluated using a quality assurance dataset collected with a single-beam echosounder; 9,386 quality assurance points were spatially joined with the mapped raster surface to compute measurement errors. The calculated median point error for Morris Lake was 0.23 foot, the median absolute error was 0.35 foot, and the 95-percent accuracy was 2.68 feet. The largest errors occurred in the steepest areas of the reservoir and in unmeasured areas. Geospatial files of the bathymetry data, including the mapped bathymetric surface, contours, and capacity tables, quality assurance points, and associated metadata are available for download as part of an accompanying U.S. Geological Survey data release.

New Jersey↗
USGS2016

Suspended-sediment concentrations, bedload, particle sizes, surrogate measurements, and annual sediment loads for selected sites in the lower Minnesota River Basin, water years 2011 through 2016

Accurate measurements of fluvial sediment are important for assessing stream ecological health, calculating flood levels, computing sediment budgets, and managing and protecting water resources. Sediment-enriched rivers in Minnesota are a concern among Federal, State, and local governments because turbidity and sediment-laden waters are the leading impairments and affect more than 6,000 miles of rivers in Minnesota. The suspended sediment in the lower Minnesota River is deleterious, contributing about 75 to 90 percent of the suspended sediment being deposited into Lake Pepin. The Saint Paul District of the U.S. Army Corps of Engineers and the Lower Minnesota River Watershed District collaborate to maintain a navigation channel on the lower 14.7 miles of the Minnesota River through scheduled dredging operations. The Minnesota Pollution Control Agency has adopted a sediment-reduction strategy to reduce sediment in the Minnesota River by 90 percent by 2040. The U.S. Geological Survey, in cooperation with the U.S. Army Corps of Engineers, the Minnesota Pollution Control Agency, and the Lower Minnesota River Watershed District, collected suspended-sediment, bedload, and particle-size samples at five sites in the lower Minnesota River Basin during water years 2011 through 2014 and surrogate measurements of acoustic backscatter at one of these sites on the lower Minnesota River during water years 2012 through 2016 to quantify sediment loads and improve understanding of sediment-transport relations. Annual sediment loads were computed for calendar years 2011 through 2014. Data collected from water years 2011 through 2014 indicated that two tributaries, Le Sueur River and High Island Creek, had the highest sediment yield and concentrations of suspended sediment. These tributaries also had greater stream gradients than the sites on the Minnesota River. Suspended fines were greater than suspended sand at all sites in the study area. The range of median particle sizes matched the range for stream gradients from greatest to smallest. Bedload ranged from 3 to 20 percent of the total load at the Le Sueur River, Minnesota River at Mankato, and High Island Creek and was less than 1 percent of the total load at the Minnesota River near Jordan and at Fort Snelling State Park. The reach of the Minnesota River between Mankato and Jordan is a major source of sediment, with the sediment yield at Jordan being two and a half times greater than at Mankato. Between Jordan and Fort Snelling, the sediment yield decreases substantially, which indicates that the Minnesota River in this reach is a sink for sediment. Surrogate measurements (acoustic backscatter) collected with suspended-sediment concentration data from water years 2012 through 2016 from the Minnesota River at Fort Snelling State Park indicated strong relations between the acoustic backscatter and suspended-sediment concentrations. These results point to the dynamic nature of sediment aggradation, degradation, and transport in the Minnesota River Basin. The analyses described in this report will improve the understanding of sediment-transport relations and sediment budgets in the Minnesota River Basin.

Minnesota↗
USGS2013

Long-range hazard assessment of volcanic ash dispersal for a Plinian eruptive scenario at Popocatépetl volcano (Mexico): implications for civil aviation safety

Popocatépetl is one of Mexico’s most active volcanoes threatening a densely populated area that includes Mexico City with more than 20 million inhabitants. The destructive potential of this volcano is demonstrated by its Late Pleistocene–Holocene eruptive activity, which has been characterized by recurrent Plinian eruptions of large magnitude, the last two of which destroyed human settlements in pre-Hispanic times. Popocatépetl’s reawakening in 1994 produced a crisis that culminated with the evacuation of two villages on the northeastern flank of the volcano. Shortly after, a monitoring system and a civil protection contingency plan based on a hazard zone map were implemented. The current volcanic hazards map considers the potential occurrence of different volcanic phenomena, including pyroclastic density currents and lahars. However, no quantitative assessment of the tephra hazard, especially related to atmospheric dispersal, has been performed. The presence of airborne volcanic ash at low and jet-cruise atmospheric levels compromises the safety of aircraft operations and forces re-routing of aircraft to prevent encounters with volcanic ash clouds. Given the high number of important airports in the surroundings of Popocatépetl volcano and considering the potential threat posed to civil aviation in Mexico and adjacent regions in case of a Plinian eruption, a hazard assessment for tephra dispersal is required. In this work, we present the first probabilistic tephra dispersal hazard assessment for Popocatépetl volcano. We compute probabilistic hazard maps for critical thresholds of airborne ash concentrations at different flight levels, corresponding to the situation defined in Europe during 2010, and still under discussion. Tephra dispersal mode is performed using the FALL3D numerical model. Probabilistic hazard maps are built for a Plinian eruptive scenario defined on the basis of geological field data for the “Ochre Pumice” Plinian eruption (4965 14 C yr BP). FALL3D model input eruptive parameters are constrained through an inversion method carried out with the semi-analytical HAZMAP model and are varied by sampling them using probability density functions. We analyze the influence of seasonal variations on ash dispersal and estimate the average persistence of critical ash concentrations at relevant locations and airports. This study assesses the impact that a Plinian eruption similar to the Ochre Pumice eruption would have on the main airports of Mexico and adjacent areas. The hazard maps presented here can support long-term planning that would help minimize the impacts of such an eruption on civil aviation.

Popocatépetl Volcano↗
USGS2021

Changes in liquefaction severity in the San Francisco Bay Area with sea-level rise

This paper studies the impacts of sea-level rise on liquefaction triggering and severity around the San Francisco Bay Area, California, for the M 7.0 “HayWired” earthquake scenario along the Hayward fault. This work emerged from stakeholder engagement for the US Geological Survey releases of the HayWired earthquake scenario and the Coastal Storm Modeling System projects, in which local planners and engineers asked where, why, and by how much liquefaction hazards may change due to sea-level rise in the future. We assess the impacts of sea-level rise on liquefaction by computing changes in liquefaction potential index (LPI) for over 400 cone penetration test (CPT) soundings around the San Francisco Bay for groundwater table models developed for current and increased sea levels of up to 5 m. For the M 7.0 HayWired earthquake scenario, we find that while the majority of sites are insensitive to sea-level changes of less than 1 m, some sites are highly sensitive to small changes in water levels. We then repeat these analyses for a uniform shaking scenario to isolate the effects of sea-level rise and we find similar patterns of change. For both earthquake scenarios, modest changes in overall LPI are expected for increases in sea level, but individual sites may see significant increases in liquefaction likelihood and severity.

Conference Paper↗
USGS1998

The probability of lava inundation at the proposed and existing Kulani prison sites

The State of Hawai`i has proposed building a 2,300-bed medium-security prison about 10 km downslope from the existing Kulani medium-security correctional facility. The proposed and existing facilities lie on the northeast rift zone of Mauna Loa, which last erupted in 1984 in this same general area. We use the best available geologic mapping and dating with GIS software to estimate the average recurrence interval between lava flows that inundate these sites. Three different methods are used to adjust the number of flows exposed at the surface for those flows that are buried to allow a better representation of the recurrence interval. Probabilities are then computed, based on these recurrence intervals, assuming that the data match a Poisson distribution. The probability of lava inundation for the existing prison site is estimated to be 11- 12% in the next 50 years. The probability of lava inundation for the proposed sites B and C are 2- 3% and 1-2%, respectively, in the same period. The probabilities are based on estimated recurrence intervals for lava flows, which are approximately proportional to the area considered. The probability of having to evacuate the prison is certainly higher than the probability of lava entering the site. Maximum warning times between eruption and lava inundation of a site are estimated to be 24 hours for the existing prison site and 72 hours for proposed sites B and C. Evacuation plans should take these times into consideration.

Open-File Report↗