SEARCH · Geology Reports

Results for “Interfaces”

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

SourceResearch connectionFrom yearThrough year

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

At least 361 records · Page 20Linked to original sources

USGS2020

Artificial nightlight alters the predator-prey dynamics of an apex carnivore

Artificial nightlight is increasingly recognized as an important environmental disturbance that influences the habitats and fitness of numerous species. However, its effects on wide‐ranging vertebrates and their interactions remain unclear. Light pollution has the potential to amplify land‐use change, and as such, answering the question of how this sensory stimulant affects behavior and habitat use of species valued for their ecological roles and economic impacts is critical for conservation and land‐use planning. Here, we combined satellite‐derived estimates of light pollution, with GPS‐data from cougars Puma concolor (n = 56), mule deer Odocoileus hemionus (n = 263) and locations of cougar‐killed deer (n = 1562 carcasses), to assess the effects of light exposure on mammal behavior and predator–prey relationships across wildland–urban gradients in the southwestern United States. Our results indicate that deer used the anthropogenic environments to access forage and were more active at night than their wildland conspecifics. Despite higher nightlight levels, cougars killed deer at the wildland–urban interface, but hunted them in the relatively darkest locations. Light had the greatest effect of all covariates on where cougars killed deer at the wildland–urban interface. Both species exhibited functional responses to light pollution at fine scales; individual cougars and deer with less light exposure increasingly avoided illuminated areas when exposed to greater radiance, whereas deer living in the wildland–urban interface selected elevated light levels. We conclude that integrating estimates of light pollution into ecological studies provides crucial insights into how the dynamic human footprint can alter animal behavior and ecosystem function across spatial scales.

Arizona, California, Colorado, Idaho, Nevada, New ↗

USGS2006

Analysis of steady-state salt-water upconing with application at Truro well field, Cape Cod, Massachusetts

Salt-water upconing describes the phenomenon where salt water is transported vertically upward under a well in response to pumpage in a fresh-water aquifer underlain by salt water. Sharp interface methods have been used successfully to describe the physics of upconing. A finite-element model is developed to simulate a sharp interface for determination of the steady-state position of the interface and maximum permissible well discharges. The model developed is compared to previous published electric-analog model results of Bennett and others (1968). Both methods are applied to a test case at Truro, Massachusetts, where maximum permissible discharges are determined by the finite-element model to range from 0.47 to 1.05 cubic feet per second for the Test Site No. 4 location.

Groundwater↗

USGS2007

Effects of sea-level rise on ground water flow in a coastal aquifer system

The effects of sea-level rise on the depth to the fresh water/salt water interface were simulated by using a density-dependent, three-dimensional numerical ground water flow model for a simplified hypothetical fresh water lens that is similar to shallow, coastal aquifers found along the Atlantic coast of the United States. Simulations of sea-level rise of 2.65 mm/year from 1929 to 2050 resulted in an increase in water levels relative to a fixed datum, yet a net decrease in water levels relative to the increased sea-level position. The net decrease in water levels was much greater near a gaining stream than farther from the stream. The difference in the change in water levels is attributed to the dampening effect of the stream on water level changes in response to sea-level rise. In response to the decreased water level altitudes relative to local sea level, the depth to the fresh water/salt water interface decreased. This reduction in the thickness of the fresh water lens varied throughout the aquifer and was greatly affected by proximity to a ground water fed stream and whether the stream was tidally influenced. Away from the stream, the thickness of the fresh water lens decreased by about 2% from 1929 to 2050, whereas the fresh water lens thickness decreased by about 22% to 31% for the same period near the stream, depending on whether the stream was tidally influenced. The difference in the change in the fresh water/salt water interface position is controlled by the difference in the net decline in water levels relative to local sea level. ?? 2007 National Ground Water Association.

Ground Water↗

USGS2011

Tidal Boundary Conditions in SEAWAT

SEAWAT, a U.S. Geological Survey groundwater flow and transport code, is increasingly used to model the effects of tidal motion on coastal aquifers. Different options are available to simulate tidal boundaries but no guidelines exist nor have comparisons been made to identify the most effective approach. We test seven methods to simulate a sloping beach and a tidal flat. The ocean is represented in one of the three ways: directly using a high hydraulic conductivity (high- K ) zone and indirect simulation via specified head boundaries using either the General Head Boundary (GHB) or the new Periodic Boundary Condition (PBC) package. All beach models simulate similar water fluxes across the upland boundary and across the sediment-water interface although the ratio of intertidal to subtidal flow is different at low tide. Simulating a seepage face results in larger intertidal fluxes and influences near-shore heads and salinity. Major differences in flow occur in the tidal flat simulations. Because SEAWAT does not simulate unsaturated flow the water table only rises via flow through the saturated zone. This results in delayed propagation of the rising tidal signal inland. Inundation of the tidal flat is delayed as is flow into the aquifer across the flat. This is severe in the high- K and PBC models but mild in the GHB models. Results indicate that any of the tidal boundary options are fine if the ocean-aquifer interface is steep. However, as the slope of that interface decreases, the high- K and PBC approaches perform poorly and the GHB boundary is preferable.

Ground Water↗

USGS1994

Crustal velocity structure of the northern Yukon-Tanana upland, central Alaska: Results from TACT refraction/wide-angle reflection data

The Fairbanks North seismic refraction/ wide-angle reflection profile, collected by the U.S. Geological Survey Trans-Alaska Crustal Transect (TACT) project in 1987, crosses the complex region between the Yukon-Tanana and Ruby terranes in interior Alaska. This region is occupied by numerous small terranes elongated in a northeast-southwest direction. These seismic data reveal a crustal velocity structure that is divided into three upper-crustal and at least two middle- to lower-crustal domains. The upper-crustal domains are delineated by two steeply dipping low-velocity anomalies that are interpreted as signatures of the Victoria Creek fault, and the Beaver Creek fault or a fault buried by the Beaver Creek fault. This tripartite upper crust extends to 8-10 km depth where a subhorizontal interface undercuts the northern and central domains. Beneath the northern domain, this interface is interpreted as the southeastwardly dipping boundary between the Tozina and Ruby terranes. The continuation of this interface beneath the central domain suggests that it may represent the detachment or basal thrust for thin-skinned tectonic amalgamation of the terranes caught between the Yukon-Tanana and Ruby terranes. The lower crust and Moho reflection exhibit differences from north to south that define at least two lower-crustal domains, interpreted as the Yukon-Tanana and Ruby terranes. Finally, the crustal thickness along the profile is nearly uniform and ranges from 31 to 34 km. Our data suggest that after initial thin-skinned amalgamation of the various terranes, this region experienced thick-skinned tectonic reorganization via strike-slip faulting. This interpretation supports a model in which at least one strand of the Tintina fault exists in this important region of Alaska.

Alaska↗

USGS2021

Subducting oceanic basement roughness impacts on upper plate tectonic structure and a backstop splay fault zone activated in the southern Kodiak aftershock region of the Mw 9.2, 1964 megathrust rupture, Alaska

In 1964, the Alaska margin ruptured in a giant Mw 9.2 megathrust earthquake, the 2nd largest during worldwide instrumental recording. The coseismic slip and aftershock region offshore Kodiak Island was surveyed in 1977 – 1981 to understand the region’s tectonics. We re-processed multichannel seismic (MCS) field data using current standard Kirchhoff depth migration and/or MCS traveltime tomography. Further surveys in 1994 added P-wave velocity structure from wide-angle seismic lines and multibeam bathymetry. Published regional gravity, backscatter, and earthquake compilations also became available at this time. Beneath the trench, rough oceanic crust is covered by ~3 to 5 km thick sediment. Sediment on the subducting plate modulates the plate interface relief. The accreted prism’s imbricate thrust faults have a complex P-wave velocity structure. Landward, an accelerated increase in P-wave velocities is marked by a backstop splay fault zone (BSFZ) that marks a transition from the prism to the higher rigidity rock beneath the middle and upper slope. Structures associated with this feature may indicate fluid flow. Further upslope, another fault extends >100 km along-strike across the middle slope. Erosion from subducting seamounts leaves embayments in the frontal prism. Plate interface roughness varies along the subduction zone. Beneath the lower and middle slope, 2.5 D plate interface images show modest relief whereas the oceanic basement image is rougher. The 1964 earthquake slip maximum coincides with the leading/landward flank of a subducting seamount and the BSFZ. The BSFZ is a potentially active structure and should be considered in tsunami hazard assessments.

Alaska↗

USGS2024

Comparing subduction ground-motion models to observations for Cascadia

We evaluate Cascadia subduction ground-motion models (GMMs), considered for the 2023 US National Seismic Hazard Model (NSHM) update, by comparing observations to model predictions. The observations comprise regional recordings from intraslab earthquakes, including contributions from 2021 and 2022 events in southern Cascadia and global records from interface earthquakes. Since the 2018 NSHM update, new GMMs for Cascadia have been published by the Next Generation Attenuation (NGA)-Subduction Project that require independent evaluation. In the regional intraslab comparisons, we highlight a characteristic frequency dependence for Cascadia data, with short periods having lower ground motions and longer periods being comparable to other subduction zones. We evaluate differences in northern and southern Cascadia and find that the NGA-Subduction GMMs developed using southern Cascadia data perform better in this region than the model that did not consider these data. We compare ground-motion variability in Cascadia with the NGA-Subduction model predictions and find differences at short periods ( T = 0.1 s) due to the use of global versus regional data in the development of these models. Moreover, the within-event component of aleatory variability from the GMMs overpredicts the standard deviation of Cascadia recordings at very short periods ( T < 0.05 s). Using global interface earthquakes as a proxy to evaluate the Cascadia GMMs, we find long-period overprediction from a simulation-based GMM and some of the empirical GMMs. When comparing recent observations, we find a similar misfit to GMMs and the 2010 and 2022 Ferndale earthquakes. Finally, we observe different basin amplification factors arising in different subsets of the data, which indicate that differences in basin factors between empirical GMMs could arise from the data selection choices by the developers. As part of evaluating the regional basin terms, we apply basin amplification factors from the magnitude 9 Cascadia earthquake simulations to the empirical GMMs for interface earthquakes. The comparisons presented in this study indicate that the NGA-Subduction GMMs for Cascadia perform well relative to observations and older subduction GMMs.

British Columbia, California, Oregon, Washington↗

USGS2006

Seismoelectric numerical modeling on a grid

Our finite-difference algorithm provides a new method for simulating how seismic waves in arbitrarily heterogeneous porous media generate electric fields through an electrokinetic mechanism called seismoelectric coupling. As the first step in our simulations, we calculate relative pore-fluid/grain-matrix displacement by using existing poroelastic theory. We then calculate the electric current resulting from the grain/fluid displacement by using seismoelectric coupling theory. This electrofiltration current acts as a source term in Poisson's equation, which then allows us to calculate the electric potential distribution. We can safely neglect induction effects in our simulations because the model area is within the electrostatic near field for the depth of investigation (tens to hundreds of meters) and the frequency ranges (10 Hz to 1 kHz) of interest for shallow seismoelectric surveys.We can independently calculate the electric-potential distribution for each time step in the poroelastic simulation without loss of accuracy because electro-osmotic feedback (fluid flow that is perturbed by generated electric fields) is at least 105 times smaller than flow that is driven by fluid-pressure gradients and matrix acceleration, and is therefore negligible. Our simulations demonstrate that, distinct from seismic reflections, the seismoelectric interface response from a thin layer (at least as thin as one-twentieth of the seismic wavelength) is considerably stronger than the response from a single interface. We find that the interface response amplitude decreases as the lateral extent of a layer decreases below the width of the first Fresnel zone. We conclude, on the basis of our modeling results and of field results published elsewhere, that downhole and/or crosswell survey geometries and time-lapse applications are particularly well suited to the seismoelectric method. ?? 2006 Society of Exploration Geophysicists.

Geophysics↗

USGS2007

Generalization of von Neumann analysis for a model of two discrete half-spaces: The acoustic case

Evaluating the performance of finite-difference algorithms typically uses a technique known as von Neumann analysis. For a given algorithm, application of the technique yields both a dispersion relation valid for the discrete time-space grid and a mathematical condition for stability. In practice, a major shortcoming of conventional von Neumann analysis is that it can be applied only to an idealized numerical model - that of an infinite, homogeneous whole space. Experience has shown that numerical instabilities often arise in finite-difference simulations of wave propagation at interfaces with strong material contrasts. These interface instabilities occur even though the conventional von Neumann stability criterion may be satisfied at each point of the numerical model. To address this issue, I generalize von Neumann analysis for a model of two half-spaces. I perform the analysis for the case of acoustic wave propagation using a standard staggered-grid finite-difference numerical scheme. By deriving expressions for the discrete reflection and transmission coefficients, I study under what conditions the discrete reflection and transmission coefficients become unbounded. I find that instabilities encountered in numerical modeling near interfaces with strong material contrasts are linked to these cases and develop a modified stability criterion that takes into account the resulting instabilities. I test and verify the stability criterion by executing a finite-difference algorithm under conditions predicted to be stable and unstable. ?? 2007 Society of Exploration Geophysicists.

Geophysics↗

USGS2003

Empirical ground-motion relations for subduction-zone earthquakes and their application to Cascadia and other regions

Ground-motion relations for earthquakes that occur in subduction zones are an important input to seismic-hazard analyses in many parts of the world. In the Cascadia region (Washington, Oregon, northern California, and British Columbia), for example, there is a significant hazard from megathrust earthquakes along the subduction interface and from large events within the subducting slab. These hazards are in addition to the hazard from shallow earthquakes in the overlying crust. We have compiled a response spectra database from thousands of strong-motion recordings from events of moment magnitude ( M ) 5-8.3 occurring in subduction zones around the world, including both interface and in-slab events. The 2001 M 6.8 Nisqually and 1999 M 5.9 Satsop earthquakes are included in the database, as are many records from subduction zones in Japan (Kyoshin-Net data), Mexico (Guerrero data), and Central America. The size of the database is four times larger than that available for previous empirical regressions to determine ground-motion relations for subduction-zone earthquakes. The large dataset enables improved determination of attenuation parameters and magnitude scaling, for both interface and in-slab events. Soil response parameters are also better determined by the data.

British Columbia, California, Oregon, Washington↗

USGS2015

Pathogen exposure varies widely among sympatric populations of wild and domestic felids across the United States

Understanding how landscape, host, and pathogen traits contribute to disease exposure requires systematic evaluations of pathogens within and among host species and geographic regions. The relative importance of these attributes is critical for management of wildlife and mitigating domestic animal and human disease, particularly given rapid ecological changes, such as urbanization. We screened >1,000 samples from sympatric populations of puma (Puma concolor), bobcat (Lynx rufus) and domestic cat (Felis catus) across urban gradients in six sites, representing three regions, in North America for exposure to a representative suite of bacterial, protozoal and viral pathogens (Bartonella sp., Toxoplasma gondii, feline herpesvirus-1, feline panleukopenea virus, feline calicivirus, feline immunodeficiency virus). We evaluated prevalence within each species, and examined host trait and land cover determinants of exposure-providing an unprecedented analysis of factors relating to potential for infections in domesticated and wild felids. Prevalence differed among host species (highest for puma and lowest for domestic cat) and was greater for indirectly transmitted pathogens. Sex was inconsistently predictive of exposure to directly transmitted pathogens only, and age infrequently predictive of both direct and indirectly transmitted pathogens. Determinants of pathogen exposure were widely divergent between the wild felid species. For puma, suburban landuse predicted increased exposure to Bartonella sp. in southern California, and FHV-1 exposure increased near urban edges in Florida. This may suggest inter-specific transmission with domestic cats via flea vectors (California) and direct contact (Florida) around urban boundaries. Bobcats captured near urban areas had increased exposure to T. gondii in Florida, suggesting an urban source of prey. Bobcats captured near urban areas in Colorado and Florida had higher FIV exposure, possibly suggesting increased intra-specific interactions through pile-up of home ranges. Beyond these regional and pathogen specific relationships, proximity to the wildland urban interface did not generally increase the probability of disease exposure in wild or domestic felids, emphasizing the importance of local ecological determinants. Indeed, pathogen exposure was often negatively associated with the wildland urban interface for all felids. Our analyses suggest cross-species pathogen transmission events around this interface may be infrequent, but followed by self-sustaining propagation within the new host species.

California, Colorado,Florida↗

USGS1970

Similarities, differences, and some genetic problems of the Wyoming and Colorado plateau types of uranium deposits in sandstone

Uranium deposits of the Wyoming roll type and the Colorado Plateau peneconcordant type are the principal domestic resources of uranium. Both types occur in lenticular sandstone beds of continental origin, have a similar suite of elements and minerals, are associated with mildly altered rock, and are thought to have formed by reduction precipitation from ground waters before significant regional deformation. They differ in several respects. The Wyoming deposits are in unlithified sandstone that is highly arkosic, whereas the Colorado Plateau deposits are in lithified sandstone that is only slightly to moderately arkosic. Uranium is the only ore metal in the Wyoming deposits, whereas vanadium or copper is more abundant than uranium in some of the Plateau deposits. The Wyoming deposits are elongate crescent-shaped bodies that extend vertically through, or partly through, a sandstone unit and which are scattered, like widely spaced beads on a string, along miles-long interfaces between oxidized (altered) and unoxidized sandstone, whereas the Plateau deposits are thin tabular layers that are nearly concordant to bedding and which occur as discrete bodies, like raisins in raisin bread, enveloped in rock altered by reduction. The Wyoming ore rolls and interfaces were dynamic, having been pushed downdip by downward-moving oxygen-bearing water that passed through the interfaces and deposited the ore minerals on the reducing side, whereas the Plateau deposits seemingly formed as static bodies, localized by intensive reducing "patches" in a mildly reducing environment. These differences focus attention on genetic problems relating to the Eh of the ore-bearing and altering solutions, the shape and localization of deposits, and the source of the uranium.

Arizona, Colorado, New Mexico, South Dakota, Utah,↗

USGS2008

Modeling soil moisture processes and recharge under a melting snowpack

Recharge into granitic bedrock under a melting snowpack is being investigated as part of a study designed to understand hydrologic processes involving snow at Yosemite National Park in the Sierra Nevada Mountains of California. Snowpack measurements, accompanied by water content and matric potential measurements of the soil under the snowpack, allowed for estimates of infiltration into the soil during snowmelt and percolation into the bedrock. During portions of the snowmelt period, infiltration rates into the soil exceeded the permeability of the bedrock and caused ponding to be sustained at the soil-bedrock interface. During a 5-d period with little measured snowmelt, drainage of the ponded water into the underlying fractured granitic bedrock was estimated to be 1.6 cm d?1, which is used as an estimate of bedrock permeability. The numerical simulator TOUGH2 was used to reproduce the field data and evaluate the potential for vertical flow into the fractured bedrock or lateral flow at the bedrock-soil interface. During most of the snowmelt season, the snowmelt rates were near or below the bedrock permeability. The field data and model results support the notion that snowmelt on the shallow soil overlying low permeability bedrock becomes direct infiltration unless the snowmelt rate greatly exceeds the bedrock permeability. Late in the season, melt rates are double that of the bedrock permeability (although only for a few days) and may tend to move laterally at the soil-bedrock interface downgradient and contribute directly to streamflow. ?? Soil Science Society of America.

Conference Paper↗

USGS2007

Surface-source downhole seismic analysis in R

This report discusses a method for interpreting a layered slowness or velocity model from surface-source downhole seismic data originally presented by Boore (2003). I have implemented this method in the statistical computing language R (R Development Core Team, 2007), so that it is freely and easily available to researchers and practitioners that may find it useful. I originally applied an early version of these routines to seismic cone penetration test data (SCPT) to analyze the horizontal variability of shear-wave velocity within the sediments in the San Francisco Bay area (Thompson et al., 2006). A more recent version of these codes was used to analyze the influence of interface-selection and model assumptions on velocity/slowness estimates and the resulting differences in site amplification (Boore and Thompson, 2007). The R environment has many benefits for scientific and statistical computation; I have chosen R to disseminate these routines because it is versatile enough to program specialized routines, is highly interactive which aids in the analysis of data, and is freely and conveniently available to install on a wide variety of computer platforms. These scripts are useful for the interpretation of layered velocity models from surface-source downhole seismic data such as deep boreholes and SCPT data. The inputs are the travel-time data and the offset of the source at the surface. The travel-time arrivals for the P- and S-waves must already be picked from the original data. An option in the inversion is to include estimates of the standard deviation of the travel-time picks for a weighted inversion of the velocity profile. The standard deviation of each travel-time pick is defined relative to the standard deviation of the best pick in a profile and is based on the accuracy with which the travel-time measurement could be determined from the seismogram. The analysis of the travel-time data consists of two parts: the identification of layer-interfaces, and the inversion for the velocity of each layer. The analyst usually picks layer-interfaces by visual inspection of the travel-time data. I have also developed an algorithm that automatically finds boundaries which can save a significant amount of the time when analyzing a large number of sites. The results of the automatic routines should be reviewed to check that they are reasonable. The interactivity of these scripts allows the user to add and to remove layers quickly, thus allowing rapid feedback on how the residuals are affected by each additional parameter in the inversion. In addition, the script allows many models to be compared at the same time.

Open-File Report↗

USGS2014

Continuous resistivity profiling and seismic-reflection data collected in April 2010 from Indian River Bay, Delaware

A geophysical survey to delineate the fresh-saline groundwater interface and associated sub-bottom sedimentary structures beneath Indian River Bay, Delaware, was carried out in April 2010. This included surveying at higher spatial resolution in the vicinity of a study site at Holts Landing, where intensive onshore and offshore studies were subsequently completed. The total length of continuous resistivity profiling (CRP) survey lines was 145 kilometers (km), with 36 km of chirp seismic lines surveyed around the perimeter of the bay. Medium-resolution CRP surveying was performed using a 50-meter streamer in a baywide grid. Results of the surveying and data inversion showed the presence of many buried paleochannels beneath Indian River Bay that generally extended perpendicular from the shoreline in areas of modern tributaries, tidal creeks, and marshes. An especially wide and deep paleochannel system was imaged in the southeastern part of the bay near White Creek. Many paleochannels also had high-resistivity anomalies corresponding to low-salinity groundwater plumes associated with them, likely due to the presence of fine-grained estuarine mud and peats in the channel fills that act as submarine confining units. Where present, these units allow plumes of low-salinity groundwater that was recharged onshore to move beyond the shoreline, creating a complex fresh-saline groundwater interface in the subsurface. The properties of this interface are important considerations in construction of accurate coastal groundwater flow models. These models are required to help predict how nutrient-rich groundwater, recharged in agricultural watersheds such as this one, makes its way into coastal bays and impacts surface-water quality and estuarine ecosystems.

Delaware↗

USGS2011

Physical and chemical characteristics including total and geochemical forms of phosphorus in sediment from the top 30 centimeters of cores collected in October 2006 at 26 sites in Upper Klamath Lake, Oregon

μThis study of phosphorus (P) cycling in eutrophic Upper Klamath Lake (UKL), Oregon, was conducted by the U.S. Geological Survey in cooperation with the U.S. Bureau of Reclamation. Lakebed sediments from the upper 30 centimeters (cm) of cores collected from 26 sites were characterized. Cores were sampled at 0.5, 1.5, 2.5, 3.5, 4.5, 10, 15, 20, 25, and 30 cm. Prior to freezing, water content and sediment pH were determined. After being freeze-dried, all samples were separated into greater than 63-micron (μm) particle-size (coarse) and less than 63-μm particle-size (fine) fractions. In the surface samples (0.5 to 4.5 cm below the sediment water interface), approximately three-fourths of the particles were larger than 63-μm. The ratios of the coarse particle-size fraction (>63 μm) and the fine particle-size fraction (<63 μm) were approximately equal in samples at depths greater than 10 cm below the sediment water interface. Chemical analyses included both size fractions of freeze-dried samples. Chemical analyses included determination of total concentrations of aluminum (Al), calcium (Ca), carbon (C), iron (Fe), poorly crystalline Fe, nitrogen (N), P, and titanium (Ti). Total Fe concentrations were the largest in sediment from the northern portion of UKL, Howard Bay, and the southern portion of the lake. Concentrations of total Al, Ca, and Ti were largest in sediment from the northern, central, and southernmost portions of the lake and in sediment from Howard Bay. Concentrations of total C and N were largest in sediment from the embayments and in sediment from the northern arm and southern portion of the lake in the general region of Buck Island. Concentrations of total C were larger in the greater than 63-μm particle-size fraction than in the less than 63-μm particle-size fraction. Sediments were sequentially extracted to determine concentrations of inorganic forms of P, including loosely sorbed P, P associated with poorly crystalline Fe oxides, and P associated with mineral phases. The difference between the concentration of total P and sum of the concentrations of inorganic forms of P is referred to as residual P. Residual P was the largest fraction of P in all of the sediment samples. In UKL, the correlation between concentrations of total P and total Fe in sediment is poor (R2<0.1). The correlation between the concentrations of total P and P associated with poorly crystalline Fe oxides is good (R2=0.43) in surface sediment (0.5-4.5 cm below the sediment water interface) but poor (R2<0.1) in sediments at depths between 10 cm and 30 cm. Phosphorus associated with poorly crystalline Fe oxides is considered bioavailable because it is released when sediment conditions change from oxidizing to reducing, which causes dissolution of Fe oxides.

Oregon↗

USGS2012

Preliminary investigation of the effects of sea-level rise on groundwater levels in New Haven, Connecticut

Global sea level rose about 0.56 feet (ft) (170 millimeters (mm)) during the 20th century. Since the 1960s, sea level has risen at Bridgeport, Connecticut, about 0.38 ft (115 mm), at a rate of 0.008 ft (2.56 mm + or - 0.58 mm) per year. With regional subsidence, and with predicted global climate change, sea level is expected to continue to rise along the northeast coast of the United States through the 21st century. Increasing sea levels will cause groundwater levels in coastal areas to rise in order to adjust to the new conditions. Some regional climate models predict wetter climate in the northeastern United States under some scenarios. Scenarios for the resulting higher groundwater levels have the potential to inundate underground infrastructure in lowlying coastal cities. New Haven is a coastal city in Connecticut surrounded and bisected by tidally affected waters. Monitoring of water levels in wells in New Haven from August 2009 to July 2010 indicates the complex effects of urban influence on groundwater levels. The response of groundwater levels to recharge and season varied considerably from well to well. Groundwater temperatures varied seasonally, but were warmer than what was typical for Connecticut, and they seem to reflect the influence of the urban setting, including the effects of conduits for underground utilities. Specific conductance was elevated in many of the wells, indicating the influence of urban activities or seawater in Long Island Sound. A preliminary steady-state model of groundwater flow for part of New Haven was constructed using MODFLOW to simulate current groundwater levels (2009-2010) and future groundwater levels based on scenarios with a rise of 3 ft (0.91 meters (m)) in sea level, which is predicted for the end of the 21st century. An additional simulation was run assuming a 3-ft rise in sea level combined with a 12-percent increase in groundwater recharge. The model was constructed from existing hydrogeologic information for the New Haven area and from new information on groundwater levels collected during October 2009-June 2010. For the scenario with a 3-ft rise in sea level and no increase in recharge, simulated groundwater levels near the coast rose 3 ft; this increased water level tapered off toward a discharge area at the only nontidal stream in the study area. Simulated stream discharge increased at the nontidal stream because of the increased gradient. Although groundwater levels rose, the simulated difference between the groundwater levels in the aquifer and the increased sea level declined, indicating that the depth to the interface between freshwater and saltwater may possibly decline. Simulated water levels were affected by rise in sea level even in areas where the water table was at 17-24 ft (5.2-7.3 m) above current (2011) sea level. For the scenario with increased recharge, simulated groundwater levels were as much as an additional foot higher at some locations in the study area. The results of this preliminary investigation indicate that groundwater levels in coastal areas can be expected to rise and may rise higher if groundwater recharge also increases. This finding has implications for the disposal of stormwater through infiltration, a low-impact development practice designed to improve water quality and reduce overland peak discharge. Other implications include increased risk of basement flooding and increased groundwater seepage into underground sewer pipes and utility corridors in some areas. These implications will present engineering challenges to New Haven and Yale University. The preliminary model developed for this study can be the starting point for further simulation of future alternative scenarios for sea-level rise and recharge. Further simulations could identify those areas of New Haven where infrastructure may be at greatest risk from rising levels of groundwater. The simulations described in this report have limitations due to the preliminary scope of the work. Approaches to improve simulations include but are not limited to incorporating: * The variable density of seawater into the model in order to understand the current and future location of the interface between freshwater and saltwater; * Collection of additional data in order to better resolve temporal and spatial patterns in water levels in the aquifer; * Improved estimates of recharge through direct and indirect measurements of freshwater discharge from the study area; and * Transient simulations for greater understanding of the amount of time required for water levels and the position of the interface between freshwater and saltwater to adjust to changes in sea level and recharge.

Connecticut↗

USGS2012

Groundwater flow and water budget in the surficial and Floridan aquifer systems in east-central Florida

A numerical transient model of the surficial and Floridan aquifer systems in east-central Florida was developed to (1) increase the understanding of water exchanges between the surficial and the Floridan aquifer systems, (2) assess the recharge rates to the surficial aquifer system from infiltration through the unsaturated zone and (3) obtain a simulation tool that could be used by water-resource managers to assess the impact of changes in groundwater withdrawals on spring flows and on the potentiometric surfaces of the hydrogeologic units composing the Floridan aquifer system. The hydrogeology of east-central Florida was evaluated and used to develop and calibrate the groundwater flow model, which simulates the regional fresh groundwater flow system. The U.S. Geological Survey three-dimensional groundwater flow model, MODFLOW-2005, was used to simulate transient groundwater flow in the surficial, intermediate, and Floridan aquifer systems from 1995 to 2006. The east-central Florida transient model encompasses an actively simulated area of about 9,000 square miles. Although the model includes surficial processes-rainfall, irrigation, evapotranspiration, runoff, infiltration, lake water levels, and stream water levels and flows-its primary purpose is to characterize and refine the understanding of groundwater flow in the Floridan aquifer system. Model-independent estimates of the partitioning of rainfall into evapotranspiration, streamflow, and aquifer recharge are provided from a water-budget analysis of the surficial aquifer system. The interaction of the groundwater flow system with the surface environment was simulated using the Green-Ampt infiltration method and the MODFLOW-2005 Unsaturated-Zone Flow, Lake, and Streamflow-Routing Packages. The model is intended to simulate the part of the groundwater system that contains freshwater. The bottom and lateral boundaries of the model were established at the estimated depths where the chloride concentration is 5,000 milligrams per liter in the Floridan aquifer system. Potential flow across the interface represented by this chloride concentration is simulated by the General Head Boundary Package. During 1995 through 2006, there were no major groundwater withdrawals near the freshwater and saline-water interface, making the general head boundary a suitable feature to estimate flow through the interface. The east-central Florida transient model was calibrated using the inverse parameter estimation code, PEST. Steady-state models for 1999 and 2003 were developed to estimate hydraulic conductivity (K) using average annual heads and spring flows as observations. The spatial variation of K was represented using zones of constant values in some layers, and pilot points in other layers. Estimated K values were within one order of magnitude of aquifer performance test data. A simulation of the final two years (2005-2006) of the 12-year model, with the K estimates from the steady-state calibration, was used to guide the estimation of specific yield and specific storage values. The final model yielded head and spring-flow residuals that met the calibration criteria for the 12-year transient simulation. The overall mean residual for heads, defining residual as simulated minus measured value, was -0.04 foot. The overall root-mean square residual for heads was less than 3.6 feet for each year in the 1995 to 2006 simulation period. The overall mean residual for spring flows was -0.3 cubic foot per second. The spatial distribution of head residuals was generally random, with some minor indications of bias. Simulated average evapotranspiration (ET) over the 1995 to 2006 period was 34.5 inches per year, compared to the calculated average ET rate of 36.6 inches per year from the model-independent water-budget analysis. Simulated average net recharge to the surficial aquifer system was 3.6 inches per year, compared with the calculated average of 3.2 inches per year from the model-independent waterbudget analysis. Groundwater withdrawals from the Floridan aquifer system averaged about 800 million gallons per day, which is equivalent to about 2 inches per year over the model area and slightly more than half of the simulated average net recharge to the surficial aquifer system over the same period. Annual net simulated recharge rates to the surficial aquifer system were less than the total groundwater withdrawals from the Floridan aquifer system only during the below-average rainfall years of 2000 and 2006.

Florida↗