SEARCH · Geology Reports

Results for “Scientific Monthly”

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 757 records · Page 42Linked to original sources

USGS2026

Analyses of meteorological and hydrological records support Tribal members’ accounts of changing climate on the Fort Apache Reservation, east–central Arizona

The Fort Apache Reservation in east–central Arizona, home to the White Mountain Apache Tribe of the Fort Apache Reservation, Arizona, contains several climate zones because of the large variation in surface elevation within the reservation. This study was carried out in cooperation with the White Mountain Apache Tribe of the Fort Apache Reservation, Arizona, to raise awareness of how the changing climate affects the Fort Apache Reservation. This report documents the evaluation of existing multidecadal meteorological and hydrological datasets for the Fort Apache Reservation, used to evaluate the effects of a changing climate on the reservation. In this evaluation, near-surface air temperature, snow depth, snow water equivalent, precipitation, and streamflow datasets were analyzed for monotonic trends indicative of changing climatic conditions during specified periods of time. The results of these trend analyses were then compared with the Tribal community's memories of the changing climate. Trend analysis of near-surface air temperatures from a U.S. Historical Climatological Network station on the Fort Apache Reservation at Whiteriver, Arizona, indicated that mean annual air temperatures have increased by an average of 2.48 degrees Fahrenheit from 1980 to 2023. Records from the same station also indicated that average monthly maximum temperatures recorded for March increased by 5.39 degrees Fahrenheit for the same time period. Annual precipitation at the five precipitation stations used in this study decreased greatly from the 1980s to 2023. The largest total decrease was 10.07 inches, or 34.7 percent. However, only one of the two precipitation stations with longer term data available prior to 1980 had a significant negative trend when data from the entire period of record, from 1901 to 2023, were analyzed. Trend analyses show a decrease in the annual maximum snow water equivalent and an earlier disappearance of the snowpack at two Natural Resources Conservation Service snow telemetry stations in the mountainous region just east of the Fort Apache Reservation from 1981 to 2023. Based on the trend analyses, the average annual maximum snow water equivalent has decreased by more than 40 percent at both stations, and the average date when the snowpack was fully melted at the stations in the spring has moved earlier in time from late April to early April or late March. However, a statistically significant trend was not determined for the early April snow water equivalent measured at a nearby Natural Resources Conservation Service snow course across its period of record, indicating that the history of mountain snowpack in this area is not fully understood. Analysis of snowfall data from a National Oceanic and Atmospheric Administration Cooperative Observer Program network station on the Fort Apache Reservation at McNary 2N, AZ (station 025412) indicated that, on average, the measured total annual snowfall at the station decreased 42.4 percent from 1935 to 2023. Streamflow data from six U.S. Geological Survey streamgages on the Fort Apache Reservation were analyzed for trends. For most streamflow gages, statistically significant trends were not determined for tested parameters when the entire streamflow period of record was used for stations with records going back to at least the 1960s. However, when the data from 1980 to 2023 was tested, most of the streamflow parameters had statistically significant negative trends. All six streamgages showed a decrease in average annual runoff of at least 50 percent from 1980 to 2023; one streamgage showed an 81.8 percent decrease. A similar statistical finding was observed in the analysis of the annual spring snowmelt peak from one of the six streamgages used in the study and located in an area receiving measurable amounts of snowmelt runoff. When data from the entire period of record (1958–2023) was used, no trend in streamflow was determined; however, a significant negative trend was determined from 1980 to 2023, indicating a decrease in average annual springtime runoff of 62.6 percent. Statistical analysis on the timing of the annual spring snowmelt peak at the same streamgage indicated the snowmelt peak is happening on average about 12 days earlier now (2023) than it did in the past. The trend results for the timing of the annual spring snowmelt peak were the same and statistically significant for both periods tested (1958–2023 and 1980–2023). Two of the streamflow records from the Fort Apache Reservation were compared to the Palmer Hydrological Drought Index computed for Arizona Climate Division 4 (East Central) by the National Centers for Environmental Information. The comparison showed that the streamflow records generally tracked the Palmer Hydrological Drought Index. In interviews, Tribal community members living on the Fort Apache Reservation described the changes in climate that they observed during their lifetimes. Common themes reported were that air temperatures have become warmer, and the weather is less predictable with changes in seasonal patterns. Drier conditions, lower snowfall, shorter winters, and lower river levels were also reported. These community member observations align with the results of this study.

Arizona↗

USGS2009

Sand resources, regional geology, and coastal processes of the Chandeleur Islands Coastal System: An evaluation of the Breton National Wildlife Refuge

Breton National Wildlife Refuge, the Chandeleur Islands chain in Louisiana, provides habitat and nesting areas for wildlife and is an initial barrier protecting New Orleans from storms. The U.S. Geological Survey (USGS) in partnership with the University of New Orleans Pontchartrain Institute for Environmental Sciences undertook an intensive study that included (1) an analysis of island change based on historical maps and remotely sensed shoreline and topographic data; (2) a series of lidar surveys at 3- to 4-month intervals after Hurricane Katrina to determine barrier island recovery potential; (3) a discussion of sea level rise and effects on the islands; (4) an analysis of sea floor evolution and sediment dynamics in the refuge over the past 150 years; (5) an assessment of the local sediment transport and sediment resource availability based on the bathymetric and subbottom data; (6) a carefully selected core collection effort to groundtruth the geophysical data and more fully characterize the sediments composing the islands and surrounds; (7) an additional survey of the St. Bernard Shoals to assess their potential as a sand resource; and (8) a modeling study to numerically simulate the potential response of the islands to the low-intensity, intermediate, and extreme events likely to affect the refuge over the next 50 years. Results indicate that the islands have become fragmented and greatly diminished in subaerial extent over time: the southern islands retreating landward as they reorganize into subaerial features, the northern islands remaining in place. Breton Island, because maintenance of the Mississippi River-Gulf Outlet (MRGO) outer bar channel requires dredging, is deprived of sand sufficient to sustain itself. Regional sediment transport trends indicate that large storms are extremely effective in transporting sand and controlling the shoreline development and barrier island geometry. Sand is transported north and south from a divergent zone near Monkey Bayou at the southern end of the Chandeleur Islands. Numerical simulation of waves and sediment transport supports the geophysical results and indicates that vast areas of the lower shoreface are affected and are undergoing erosion during storm events, that there is little or no fair weather mechanism to rework material into the littoral system, and that as a result, there is a net loss of sediment from the system. Lidar surveys revealed that the island chain immediately after Hurricane Katrina lost about 84 percent of its area and about 92 percent of its prestorm volume. Marsh platforms that supported the islands’ sand prior to the storm were reduced in width by more than one-half. Repeated lidar surveys document that in places the shoreline has retreated about 100 m under the relatively low-energy waves since Hurricanes Katrina and Rita; however, this retreat is nonuniform. Recent high-resolution geophysical surveys of the sea floor and subsurface within 5–6 km of the Chandeleur Islands during 2006 and 2007 show that, in addition to the sand that is rebuilding portions of the island chain, a large volume of sand is contained in Hewes Point, in an extensive subtidal spit platform that has formed at the northern end of the Chandeleur Islands. Hewes Point appears to be the depositional terminus of the alongshore transport system. In the southern Chandeleurs, sand is being deposited in a broad tabular deposit near Breton Island called the southern offshore sand sheet. These two depocenters account for approximately 70 percent of the estimated sediment volume located in potential borrow sites. An additional large potential source of sand for restoration lies in the St. Bernard Shoals, which are estimated to contain approximately 200 × 10 6 m 3 of sand. Successful restoration planning for the Breton National Wildlife Refuge should mimic the natural processes of early stages of barrier island evolution including lateral transport to the flanks of the island chain from a centralized sand source that will ultimately enhance the ability of the islands to naturally build backbarrier marsh, dunes, and a continuous sandy shoreline. Barrier island sediment nourishment should be executed with the understanding that gulf shoreline erosion is inevitable but that island area can be maintained and enhanced during retreat (thus significantly prolonging the life of the island chain) with strategic sand placement.

Louisiana↗

USGS2022

Long-term groundwater availability in the Waihe‘e, ‘Īao, and Waikapū aquifer systems, Maui, Hawai‘i

Groundwater levels have declined since the 1940s in the Wailuku area of central Maui, Hawai‘i, on the eastern flank of West Maui volcano, mainly in response to increased groundwater withdrawals. Available data since the 1980s also indicate a thinning of the freshwater lens and an increase in chloride concentrations of pumped water from production wells. These trends, combined with projected increases in demand for groundwater in central Maui, have led to concerns over groundwater availability and have highlighted a need to improve general understanding of the hydrologic effects of proposed groundwater withdrawals in the Waihe‘e, ‘Īao, and Waikapū areas of central Maui. A numerical groundwater model was constructed to simulate the flow and salinity of groundwater in central Maui. The model simulates the effects of changes in groundwater withdrawals and recharge on water levels, freshwater-lens thicknesses, and chloride concentrations of pumped water from production wells. The model incorporates updated water-budget estimates of groundwater recharge from infiltration and direct recharge, seepage in stream channels, and inflow from inland areas. Mean annual groundwater recharge from infiltration and direct recharge was estimated using a daily water-budget model and the most current data, including the distributions of monthly rainfall and potential evapotranspiration, for the study area for nine historical periods from 1926 through 2012: 1926–69, 1970–79, 1980–84, 1985–89, 1990–94, 1995–99, 2000–04, 2005–09, and 2010–12. The water-budget model also estimated groundwater recharge based on one hypothetical scenario that used 1980–2010 rainfall and 2017 land cover. For the nine historical periods, estimated recharge from infiltration and direct recharge within the area of the groundwater model ranged from 30.4 million gallons per day (Mgal/d) during 2010–12 to 98.7 Mgal/d during 1926–69. Variability in recharge during these periods mainly reflects changes in rainfall and irrigation over time. Between 2010 and 2014, streamflow restoration in previously diverted streams resulted in an estimated increase in recharge from seepage in stream channels of about 12.5 Mgal/d. Average groundwater inflow of about 39.6 Mgal/d from inland, dike-intruded areas to the main area of interest was estimated from an existing island-wide numerical groundwater-flow model, which is at a larger scale and incorporates a greater number of simplifying assumptions. The numerical groundwater model developed for this study was calibrated to 1926–2012 transient water levels, vertical salinity profiles, and chloride concentrations of water pumped by production wells in the study area. The model was then used to evaluate one future recharge and six selected withdrawal scenarios, developed in consultation with the Maui Department of Water Supply, in terms of long-term changes in water level and 50-percent ocean-water salinity surface. The groundwater model was also used to simulate the future salinity of water withdrawn by existing and proposed production wells. The simulations were run to steady-state conditions, providing an estimate of the long-term effects of changes in withdrawal and recharge on the groundwater resource. Results of the simulated future withdrawal scenarios indicate that, relative to 2017–18 rates, the scenarios’ long-term effect of increased withdrawals ultimately leads to lower water levels and a higher 50-percent ocean-water salinity surface indicating a thinning of the freshwater lens. Results also indicate that the increased withdrawals produce some groundwater with chloride concentration below 250 milligrams per liter and some groundwater with higher chloride concentration. The amount of drawdown near production wells and the quality of water withdrawn from production wells is dependent on the rate and spatial distribution of the withdrawals. The model was also used to evaluate how groundwater availability may be affected for a drier recharge scenario based on a published study of future climate. Model results of the future recharge scenario indicate that the rate of groundwater recharge is a controlling factor for (1) water levels, (2) the 50-percent ocean-water salinity surface, and (3) the quality of water withdrawn from production wells in the Wailuku area. Coupled with reduced groundwater recharge (with all other factors remaining equal), the modeled future withdrawals in the scenario would tend to cause lower water levels, a higher 50-percent ocean-water salinity surface, and increased salinity of water withdrawn from production wells. The three-dimensional numerical groundwater model developed for this study utilizes the latest available hydrologic and geologic information and is a useful tool for understanding the long-term hydrologic effects of additional groundwater withdrawals in central Maui. The model has several limitations, including its non-uniqueness and inability to account for local-scale heterogeneities. Short-term effects of changes in recharge and withdrawals—and optimization of pumping rates to meet increased demand for water with acceptable salinity—are possible conditions for future simulation analyses.

Hawaii↗

USGS2017

Geologic map of the northeast flank of Mauna Loa volcano, Island of Hawai'i, Hawaii

Summary Mauna Loa, the largest volcano on Earth, has erupted 33 times since written descriptions became available in 1832. Some eruptions were preceded by only brief seismic unrest, while others followed several months to a year of increased seismicity. The majority of the eruptions of Mauna Loa began in the summit area (>12,000-ft elevation; Lockwood and Lipman, 1987); yet the Northeast Rift Zone (NERZ) was the source of eight flank eruptions since 1843 (table 1). This zone extends from the 13,680-ft-high summit towards Hilo (population ~60,000), the second largest city in the State of Hawaii. Although most of the source vents are farther than 30 km away, the 1880 flow from one of the vents extends into Hilo, nearly reaching Hilo Bay. The city is built entirely on flows erupted from the NERZ, most older than that erupted in 1843. Once underway, Mauna Loa's eruptions can produce lava flows that reach the sea in less than 24 hours, severing roads and utilities in their path. For example, lava flows erupted from the Southwest Rift Zone (SWRZ) in 1950 advanced at an average rate of 9.3 km per hour, and all three lobes reached the ocean within approximately 24 hours (Finch and Macdonald, 1953). The flows near the eruptive vents must have traveled even faster. In terms of eruption frequency, pre-eruption warning, and rapid flow emplacement, Mauna Loa poses an enormous volcanic-hazard threat to the Island of Hawai‘i. By documenting past activity and by alerting the public and local government officials of our findings, we can anticipate the volcanic hazards and substantially mitigate the risks associated with an eruption of this massive edifice. From the geologic record, we can deduce several generalized facts about the geologic history of the NERZ. The middle to the uppermost section of the rift zone were more active in the past 4,000 years than the lower part, perhaps due to buttressing of the lower east rift zone by Mauna Kea and Kīlauea volcanoes. The historical flows that erupted on the north flank of the rift zone, which is more vulnerable to inundation, advanced toward Hilo. Lockwood (1990) noted that the vents of historical activity are migrating to the south. The volcano appears to have a self-regulating mechanism that evenly distributes long-term activity across its flanks. The geologic record also supports this notion; the time prior to the historical period (Age Group 1, orange units, pre-A.D. 1843–1,000 yr B.P.; see map sheet 2) is dominated by activity on the south side of the NERZ. The NERZ trends N. 65° E. and is about 40 km long and 2–4 km wide, narrowing at the summit caldera. It becomes diffuse (6–7 km wide) at its down-rift terminus, at the approximately 3,400-ft elevation. Its constructional crest is marked by low spatter ramparts and by spatter cones as high as 60 m. Subparallel eruptive fissures and ground cracks cut vent deposits and flows in and near the rift crest. Lava typically flows to the north, east, or south, depending on vent location relative to the rift crest. Encompassing 1,140 km 2 of the northeast flank of Mauna Loa from the 10,880-ft elevation to sea level, the map covers the area from Hilo to Volcano on the east and includes the rift zone from Puu Ulaula quadrangle in the southwest to Hilo in the northeast. The distribution of 105 eruptive units (flows)—separated into 15 age groups ranging from more than 30,000 years B.P. to A.D. 1984—are shown, as well as the relations of volcanic and surficial sedimentary deposits. This map incorporates previously reported work published in generalized small-scale maps (Lockwood and Lipman, 1987; Buchanan-Banks, 1993; Lockwood, 1995; and Wolfe and Morris, 1996).

Hawaii↗

USGS2013

Groundwater contributions of flow, nitrate, and dissolved organic carbon to the lower San Joaquin River, California, 2006-08

The influence of groundwater on surface-water quality in the San Joaquin River, California, was examined for a 59-mile reach from the confluence with Salt Slough to Vernalis. The primary objective of this study was to quantify the rate of groundwater discharged to the lower San Joaquin River and the contribution of nitrate and dissolved organic carbon concentrations to the river. Multiple lines of evidence from four independent approaches were used to characterize groundwater contributions of nitrogen and dissolved organic carbon. Monitoring wells (in-stream and bank wells), streambed synoptic surveys (stream water and shallow groundwater), longitudinal profile surveys by boat (continuous water-quality parameters in the stream), and modeling (MODFLOW and VS2DH) provided a combination of temporal, spatial, quantitative, and qualitative evidence of groundwater contributions to the river and the associated quality. Monitoring wells in nested clusters in the streambed (in-stream wells) and on both banks (bank wells) along the river were monitored monthly from September 2006 to January 2009. Nitrate concentrations in the bank wells ranged from less than detection—that is, less than 0.01 milligrams per liter (mg/L) as nitrogen (N)—to approximately 13 mg/L as N. Nitrate was not detected at 17 of 26 monitoring wells during the study period. Dissolved organic carbon concentrations among monitoring wells were highly variable, but they generally ranged from 1 to 4 mg/L. In a previous study, 14 bank wells were sampled once in 1988 following their original installation. With few exceptions, specific conductivity and nitrate concentrations measured in this study were virtually identical to those measured 20 years ago. Streambed synoptic measurements were made by using a temporarily installed drive-point piezometer at 113 distinct transects across the stream during 4 sampling events. Nitrate concentrations exceeded the detection limit of 0.01 mg/L as N in 5 percent of groundwater samples collected from the in-stream wells as part of the synoptic surveys. Only 7 of the 113 cross-sectional transects had nitrate concentrations greater than 1 mg/L as N. In contrast, surface waters in the San Joaquin River tended to have nitrate concentrations in the 1–3 mg/L as N range. A zone of lower oxygen (less than 2 mg/L) in the streambed could limit nitrate contributions from regional groundwater flow because nitrate can be converted to nitrogen gas within this zone. Appreciable concentrations of ammonium (average concentration was 1.92 mg/L as N, and 95th percentile was 10.34 mg/L as N) in the shallow groundwater, believed to originate from anoxic mineralization of streambed sediments, could contribute nitrogen to the overlying stream as nitrate following in-stream nitrification, however. Dissolved organic carbon concentrations were highly variable in the shallow groundwater below the river (1 to 6 ft below streambed) and generally ranged between 1 and 5 mg/L, but had maximum concentrations in the 15–25 mg/L range. The longitudinal profile surveys were not particularly useful in identifying groundwater discharge areas. However, the longitudinal approach described in this report was useful as a baseline survey of measured water-quality parameters and for identifying tributary inflows that affect surface-water concentrations of nitrate. Results of the calibrated MODFLOW model indicated that the simulated groundwater discharge rate was approximately 1.0 cubic foot per second per mile (cfs/mi), and the predominant horizontal groundwater flow direction between the deep bank wells was westward beneath the river. The modeled (VS2DH) flux values (river gain versus river loss) were calculated for the irrigation and non-irrigation season, and these fluxes were an order of magnitude less than those from MODFLOW. During the irrigation season, the average river gain was 0.11 cfs/mi, and the average river loss was −0.05 cfs/mi. During the non-irrigation season, the average river gain was 0.10 cfs/mi, and the average river loss was -0.08 cfs/mi. Information on groundwater interactions and water quality collected for this study was used to estimate loads of nitrate and dissolved organic carbon from the groundwater to the San Joaquin River. Estimated loads of dissolved inorganic nitrogen and dissolved organic carbon were calculated by using concentrations measured during four streambed synoptic surveys and the estimated groundwater discharge rate to the San Joaquin River from MODFLOW of 1 cfs/mi. The estimated groundwater loads to the San Joaquin River for dissolved inorganic nitrogen and dissolved organic carbon were 300 and 350 kilograms per day, respectively. These loads represent 9 and 7 percent, respectively, of the estimated instantaneous surface-water loads for dissolved inorganic nitrogen and dissolved organic carbon at the most downstream site, Vernalis, measured during the four streambed synoptic surveys.

California↗

USGS2014

Simulation of the Lower Walker River Basin hydrologic system, west-central Nevada, using PRMS and MODFLOW models

Walker Lake is a terminal lake in west-central Nevada with almost all outflow occurring through evaporation. Diversions from Walker River since the early 1900s have contributed to a substantial reduction in flow entering Walker Lake. As a result, the lake is receding, and salt concentrations have increased to a level in which Oncorhynchus clarkii henshawi (Lahontan Cutthroat trout) are no longer present, and the lake ecosystem is threatened. Consequently, there is a concerted effort to restore the Walker Lake ecosystem and fishery to a level that is more sustainable. However, Walker Lake is interlinked with the lower Walker River and adjacent groundwater system which makes it difficult to understand the full effect of upstream water-management actions on the overall hydrologic system including the lake level, volume, and dissolved-solids concentrations of Walker Lake. To understand the effects of water-management actions on the lower Walker River Basin hydrologic system, a watershed model and groundwater flow model have been developed by the U.S. Geological Survey in cooperation with the Bureau of Reclamation and the National Fish and Wildlife Foundation. The watershed model was developed using the precipitation runoff modeling system (PRMS) and the groundwater flow model was constructed using the MODular groundwater FLOW model (MODFLOW) and both were calibrated for the lower Walker River Basin. These models can be incorporated in an integrated Groundwater and Surface-water FLOW (GSFLOW) model of the lower Walker River Basin. Additionally, the MODFLOW model developed for this study is useful for efficiently simulating long-term and large-scale effects of water-management actions on groundwater hydrology, streamflow, and Walker Lake level, volume, and dissolved-solids concentrations. The lower Walker River Basin PRMS model (LWR_PRMS) was constructed using a subbasin approach to aid in development and calibration, and simulates a 30-year period from 1978 to 2007 using daily time steps. The LWR_PRMS was used to estimate the distribution of groundwater recharge specified in the MODFLOW model. The highest rates of groundwater recharge occur in the Wassuk Range beneath perennial and ephemeral stream channels, whereas lower rates of recharge occur beneath alluvial fans along mountain fronts. The total groundwater recharge estimated using PRMS was about 25,000 acre-feet per year. The lower Walker River Basin MODFLOW (LWR_MF) model simulates an 89-year period using monthly time steps. The LWR_MF was constructed with an initial steady-state simulation to represent dynamic equilibrium conditions from 1908 to 1918 and then a transient simulation representing the period 1919–2007. The model was calibrated using a combination of manual and automated methods of adjusting model parameters to minimize errors between model simulated results and weighted observations of groundwater levels, streamflows, and lake level. Hydrologic conditions simulated with the LWR_MF include the movement and change in storage of groundwater, and the water budgets for Walker River, Walker Lake, and the groundwater system. The LWR_MF computed dissolved-solids concentrations for Walker Lake using simulated lake volume and an assumed constant internal salt mass of 37.2 million tons. Effects of potential changes in water management on future conditions (scenarios) of the lower Walker River Basin hydrologic system and Walker Lake from 2011 to 2070 were evaluated. Several water-management scenarios were considered, including a baseline scenario that represents no changes in system management, improved irrigation efficiencies for the Walker River Indian Irrigation Project (WRIIP), a range of increased streamflows entering the lower Walker River Basin, and, the fallowing of fields on the WRIIP. For the baseline scenario, it was assumed that streamflow conditions from 1981 to 2010 will be repeated in the future. Results indicate that Walker Lake level and volume continue to decline but at a slower rate as the surface area of the lake becomes smaller and lake evaporation decreases. Dissolved-solids concentrations in Walker Lake continue to increase and increase much more rapidly during periods when minimal flows reach the lake due to a diminished lake volume. Alternatively, in years with high runoff, lake level increases are greater and dissolved-solids decreases are greater, compared with equivalent runoffs experienced during 1981–2010. The simulated effects of improving WRIIP efficiencies on Walker River streamflows, Walker Lake inflow, level, and dissolved-solids concentrations, and crop consumptive use, are compared with the baseline reference scenario for a range of irrigation efficiency improvements from 0 to 25 percent over 60 years. Results indicate that water is conserved through a reduction in irrigation-induced groundwater recharge and subsequent groundwater discharge through evapotranspiration. The conserved water mostly goes to increased streamflow to Walker Lake, followed by increased crop consumptive use, then increased evaporation from Weber Reservoir. The simulated effects of increased streamflows at Walker River at Wabuska streamgage (10301500) on Walker Lake inflow, level, and dissolved-solids concentrations, and crop consumptive use, are compared with the baseline scenario after 60 years under two different management methods for Weber Reservoir. Results indicate Walker Lake level and dissolved-solids concentrations stabilized with increased irrigation-season streamflow of about 40,000 acre-feet per year at the Walker River at Wabuska streamgage. Walker Lake level increased, and dissolved-solids concentration decreased, with increased flows of 50,000 acre-feet per year or more. After 60 years with additional irrigation-season streamflows of 50,000 acre-feet per year, Walker Lake level increased by about 48 feet, and lake dissolved-solids concentrations decreased by about 3,000 milligrams per liter (mg/L). With 75,000 acre-feet per year of additional streamflow, Walker Lake level increased by 70 feet, and dissolved-solids concentration decreased by 7,600 milligrams per liter. The effects of fallowing of Walker River Indian Irrigation Project fields from 2007 to 2010 on Walker Lake inflow, level, and dissolved solids were evaluated. Fallowing resulted in a near doubling of Walker River inflow to Walker Lake during this period, an increase in Walker Lake level of about 1.4 feet, and a decrease in dissolved-solids concentration of about 540 mg/L.

Nevada↗

USGS2014

Water-quality variability and constituent transport and processes in streams of Johnson County, Kansas, using continuous monitoring and regression models, 2003-11

The population of Johnson County, Kansas increased by about 24 percent between 2000 and 2012, making it one of the most rapidly developing areas of Kansas. The U.S. Geological Survey, in cooperation with the Johnson County Stormwater Management Program, began a comprehensive study of Johnson County streams in 2002 to evaluate and monitor changes in stream quality. The purpose of this report is to describe water-quality variability and constituent transport for streams representing the five largest watersheds in Johnson County, Kansas during 2003 through 2011. The watersheds ranged in urban development from 98.3 percent urban (Indian Creek) to 16.7 percent urban (Kill Creek). Water-quality conditions are quantified among the watersheds of similar size (50.1 square miles to 65.7 square miles) using continuous, in-stream measurements, and using regression models developed from continuous and discrete data. These data are used to quantify variability in concentrations and loads during changing streamflow and seasonal conditions, describe differences among sites, and assess water quality relative to water-quality standards and stream management goals. Water quality varied relative to streamflow conditions, urbanization in the upstream watershed, and contributions from wastewater treatment facilities and storm runoff. Generally, as percent impervious surface (a measure of urbanization) increased, streamflow yield increased. Water temperature of Indian Creek, the most urban site which is also downstream from wastewater facility discharges, was higher than the other sites about 50 percent of the time, particularly during winter months. Dissolved oxygen concentrations were less than the Kansas Department of Health and Environment minimum criterion of 5 milligrams per liter about 15 percent of the time at the Indian Creek site. Dissolved oxygen concentrations were less than the criterion about 10 percent of the time at the rural Blue River and Kill Creek sites, and less than 5 percent of the time at the other sites. Low dissolved oxygen at all sites generally coincided with lowest streamflow and warmer water temperatures. Hourly dissolved oxygen concentrations less than 5 milligrams per liter were measured at all sites every year, indicating that even under normal climate conditions in non-urban watersheds such as Kill Creek, dissolved oxygen concentrations may not meet State aquatic-life criterion. Specific conductance was nearly always highest in Indian and Mill Creeks, which were the most urban streams with the largest upstream discharges from wastewater treatment facilities. The largest chloride concentrations and variability were recorded at urban sites and during winter. Each winter during the study period, chloride concentrations in the most urban site, Indian Creek, exceeded the U.S. Environmental Protection Agency-recommended criterion of 230 milligrams per liter for at least 10 consecutive days. The U.S. Environmental Protection Agency-recommended ecoregion criterion for turbidity was exceeded 30 (Indian Creek) to 50 (Blue River) percent of the time. The highest average annual streamflow-weighted suspendedsediment concentration during the study period was in Mill Creek, which has undergone rapid development that likely contributed to higher sediment concentrations. One of the largest suspended-sediment load events in Indian Creek was recorded in early May 2007 when about 25 percent of the total annual sediment load was transported during a period of about 2.25 days. A simultaneous load event was recorded in Kill Creek, when about 75 percent of the total annual sediment load was transported. Sediment yields generally increased as percent impervious surface increased. Computed hourly total nitrogen and total phosphorus concentrations and yields and streamflow-weighted concentrations generally were largest in Indian and Mill Creeks. Annual percent contribution of total nitrogen in the Blue River from wastewater treatment facility discharges ranged from 19 percent in 2010 to 60 percent in 2006. Annual percent contribution of total nitrogen in Indian Creek from wastewater treatment facility discharges ranged from 35 percent in 2010 to 93 percent in 2006. The largest percent nutrient contributions from wastewater discharges coincided with the smallest annual precipitation and streamflow volume, resulting in less contribution originating from runoff. Fecal indicator bacteria Escherichia coli density at the urban Indian Creek site was usually the largest of the five monitoring sites, with an annual median density that consistently exceeded the State primary contact criterion value but was less than the secondary contact criterion. Less than 1 percent of the total annual bacteria load in the Blue River and Indian Creek originated from wastewater discharges, except during 2006 when about 6 percent of the Indian Creek load originated from wastewater. Continuous water-quality monitoring provides a foundation for comprehensive evaluation and understanding of variability and loading characteristics in streams in Johnson County. Because several directly measured parameters are strongly correlated with particular constituents of interest, regression models provide a valuable tool for evaluating variability and loading on the basis of computed continuous data. Continuous data are particularly useful for characterizing nonpoint-source contributions from stormwater runoff. Transmission of continuous data in real-time makes it possible to rapidly detect and respond to potential environmental concerns. As monitoring technologies continue to improve, so does the ability to monitor additional constituents of interest, with smaller measurement error, and at lower operational cost. Continuous water-quality data including model information and computed concentrations and loads during the study period are available at http://nrtwq.usgs.gov/ks/ .

Kansas↗

USGS2014

Simulation of the effects of rainfall and groundwater use on historical lake water levels, groundwater levels, and spring flows in central Florida

The urbanization of central Florida has progressed substantially in recent decades, and the total population in Lake, Orange, Osceola, Polk, and Seminole Counties more than quadrupled from 1960 to 2010. The Floridan aquifer system is the primary source of water for potable, industrial, and agricultural purposes in central Florida. Despite increases in groundwater withdrawals to meet the demand of population growth, recharge derived by infiltration of rainfall in the well-drained karst terrain of central Florida is the largest component of the long-term water balance of the Floridan aquifer system. To complement existing physics-based groundwater flow models, artificial neural networks and other data-mining techniques were used to simulate historical lake water level, groundwater level, and spring flow at sites throughout the area. Historical data were examined using descriptive statistics, cluster analysis, and other exploratory analysis techniques to assess their suitability for more intensive data-mining analysis. Linear trend analyses of meteorological data collected by the National Oceanic and Atmospheric Administration at 21 sites indicate 67 percent of sites exhibited upward trends in air temperature over at least a 45-year period of record, whereas 76 percent exhibited downward trends in rainfall over at least a 95-year period of record. Likewise, linear trend analyses of hydrologic response data, which have varied periods of record ranging in length from 10 to 79 years, indicate that water levels in lakes (307 sites) were about evenly split between upward and downward trends, whereas water levels in 69 percent of wells (out of 455 sites) and flows in 68 percent of springs (out of 19 sites) exhibited downward trends. Total groundwater use in the study area increased from about 250 million gallons per day (Mgal/d) in 1958 to about 590 Mgal/d in 1980 and remained relatively stable from 1981 to 2008, with a minimum of 559 Mgal/d in 1994 and a maximum of 773 Mgal/d in 2000. The change in groundwater-use trend in the early 1980s and the following period of relatively slight trend is attributable to the concomitant effects of increasing public-supply withdrawals and decreasing use of water by the phosphate industry and agriculture. On the basis of available historical data and exploratory analyses, empirical lake water-level, groundwater-level, and spring-flow models were developed for 22 lakes, 23 wells, and 6 springs. Input time series consisting of various frequencies and frequency-band components of daily rainfall (1942 to 2008) and monthly total groundwater use (1957 to 2008) resulted in hybrid signal-decomposition artificial neural network models. The final models explained much of the variability in observed hydrologic data, with 43 of the 51 sites having coefficients of determination exceeding 0.6, and the models matched the magnitude of the observed data reasonably well, such that models for 32 of the 51 sites had root-mean-square errors less than 10 percent of the measured range of the data. The Central Florida Artificial Neural Network Decision Support System was developed to integrate historical databases and the 102 site-specific artificial neural network models, model controls, and model output into a spreadsheet application with a graphical user interface that allows the user to simulate scenarios of interest. Overall, the data-mining analyses indicate that the Floridan aquifer system in central Florida is a highly conductive, dynamic, open system that is strongly influenced by external forcing. The most important external forcing appears to be rainfall, which explains much of the multiyear cyclic variability and long-term downward trends observed in lake water levels, groundwater levels, and spring flows. For most sites, groundwater use explains less of the observed variability in water levels and flows than rainfall. Relative groundwater-use impacts are greater during droughts, however, and long-term trends in water levels and flows were identified that are consistent with historical groundwater-use patterns. The sensitivity of the hydrologic system to rainfall is expected, owing to the well-drained karst terrain and relatively thin confinement of the Floridan aquifer system in much of central Florida. These characteristics facilitate the relatively rapid transmission of infiltrating water from rainfall to the water table and contribute to downward leakage of water to the Floridan aquifer system. The areally distributed nature of rainfall, as opposed to the site-specific nature of groundwater use, and the generally high transmissivity and low storativity properties of the semiconfined Floridan aquifer system contribute to the prevalence of water-level and flow patterns that mimic rainfall patterns. In general, the data-mining analyses demonstrate that the hydrologic system in central Florida is affected by groundwater use differently during wet periods, when little or no system storage is available (high water levels), compared to dry periods, when there is excess system storage (low water levels). Thus, by driving the overall behavior of the system, rainfall indirectly influences the degree to which groundwater use will effect persistent trends in water levels and flows, with groundwater-use impacts more prevalent during periods of low water levels and spring flows caused by low rainfall and less prevalent during periods of high water levels and spring flows caused by high rainfall. Differences in the magnitudes of rainfall and groundwater use during wet and dry periods also are important determinants of hydrologic response. An important implication of the data-mining analyses is that rainfall variability at subannual to multidecadal timescales must be considered in combination with groundwater use to provide robust system-response predictions that enhance sustainable resource management in an open karst aquifer system. The data-driven approach was limited, however, by the confounding effects of correlation between rainfall and groundwater use, the quality and completeness of the historical databases, and the spatial variations in groundwater use. The data-mining analyses indicate that available historical data when used alone do not contain sufficient information to definitively quantify the related individual effects of rainfall and groundwater use on hydrologic response. The knowledge gained from data-driven modeling and the results from physics-based modeling, when compared and used in combination, can yield a more comprehensive assessment and a more robust understanding of the hydrologic system than either of the approaches used separately.

Florida↗

USGS2021

Characterization of factors affecting groundwater levels in and near the former Lake Traverse Indian Reservation, South Dakota, water years 1956–2017

The U.S. Geological Survey (USGS), in cooperation with the Sisseton Wahpeton Oyate, completed a study to characterize water-level fluctuations in observation wells relative to driving factors that affect water levels in and near the historical 1867 boundary of the Lake Traverse Indian Reservation. The study investigated concerns regarding potential effects of groundwater withdrawals and climate conditions on groundwater levels within an area that includes the historical boundary of the reservation and a surrounding area that extends 10 miles in all directions within South Dakota. Characterization of water-level fluctuations in observation wells and relative driving factors was accomplished by statistical trend analysis. Monthly data from the Parameter-elevation Regressions on Independent Slopes Model (PRISM) were aggregated to obtain annual and seasonal datasets for total precipitation, minimum air temperature ( T min ), and maximum air temperature ( T max ) for the study area and a surrounding buffer area. Trend tests for gridded data for total precipitation, T min , and T max were completed for annual and seasonal time series for water years 1956–2017, which is about 2 years before the earliest available water-level measurements. A 2-year offset was arbitrarily selected because scrutiny of water-level and precipitation data indicated that responses of groundwater levels for many of the observation wells lagged major changes in precipitation patterns by about 2 years. Statistically significant upward trends were detected for annual precipitation and annual T min for most of the study area and the surrounding buffer area. Statistically significant downward trends in T max were detected for only a few 2.5 arc-minute grid cells; however, the sparsity of the spatial coverage reduces confidence that these are true trends, in contrast to the near completeness of the spatial coverage in upward trends for T min . Spatial distributions of statistically significant trends in seasonal climate data were generally similar to the annual trends, but with substantial differences in the spatial density of the trends. Potential interactions among water levels in observation wells and streamflow were examined through correlation analyses of the annual median water level for each of 76 observation wells versus the annual mean streamflow for each of four area streamgages. Potential interactions among water levels in observation wells and lake levels were examined through correlation analyses involving 25 area lakes. Resulting correlation coefficients were used as part of an approach for selecting a lake to be plotted in conjunction with water-level and precipitation data for each observation well. Groundwater trends for 76 observation wells were analyzed for three separate water-level parameters (minimum, median, and maximum) because wells are measured sporadically, and data are biased towards more frequent measurements during periods of heaviest irrigation demand. Trends in the time series of annual precipitation (from PRISM) starting 2 years earlier than the associated water-level trend also were analyzed for the location of each individual observation well. Sen’s slope and Mann-Kendall p -values were computed for the three water-level parameters and for the annual precipitation time series. Graphs showing results of trend analyses for each observation well also showed changes with time in the sum of licensed groundwater withdrawals within six specified radii (0.5, 1.0, 2.0, 3.0, 4.0, and 5.0 miles) of each well as a qualitative indicator of proximal groundwater demand. Trends in groundwater levels in observation wells in the study area are predominantly upward, with 43 of 76 wells having significant upward trends for at least one of the three water-level parameters and only 8 wells having significant downward trends for at least one water-level parameter. The upward groundwater trends are driven by predominantly upward precipitation trends, with 43 wells (not all the same wells) also having significant upward trends and no wells having significant downward trends. Significant upward precipitation trends were detected for only two of the eight wells with significant downward groundwater trends. Groundwater levels in some observation wells likely are also substantially affected by interactions with surface water, especially with lakes. Water levels in many area lakes increased in response to wet conditions of the early 1990s and have maintained high water levels ever since. It is recognized that in many cases lakes that were selected for plotting with groundwater hydrographs likely are not hydraulically connected with a groundwater system or aquifer associated with an individual well; however, interactions also are plausible for numerous other lakes for which water-level records are not available.

South Dakota↗

USGS2008

Drier forest composition associated with hydrologic change in the Apalachicola River floodplain, Florida

Forests of the Apalachicola River floodplain had shorter flood durations, were drier in composition, and had 17 percent fewer trees in 2004 than in 1976. The change to drier forest composition is expected to continue for at least 80 more years. Floodplain drying was caused by large declines in river levels resulting from erosion of the river channel after 1954 and from decreased flows in spring and summer months since the 1970s. Water-level declines have been greatest at low and medium flows, which are the most common flows (occurring about 80 percent of the time). Water levels have remained relatively unchanged during large floods which continue to occur about three times per decade. A study conducted by the U.S. Geological Survey compared temporal changes in hydrologic conditions, forest composition, forest characteristics, and individual species of trees, as well as estimated the potential for change in composition of floodplain forests in the nontidal reach of the Apalachicola River. The study was conducted with the cooperation of the Florida Department of Environmental Protection and the Northwest Florida Water Management District. Forest composition and field observations from studies conducted in 1976-1984 (termed “1976 data”) were used as baseline data for comparison with data from plots sampled in 2004-2006 (“2004 data”). Flood durations were shorter in all periods subsequent to 1923-1976. The periods of record used to calculate flood durations for forest data were subsets of the complete record available (1923-2004). At sampled plots in all forest types and reaches combined, flood durations changed an average of more than 70 percent toward the baseline flood duration of the next drier forest type. For all forest types, changes in flood durations toward the next drier type were greatest in the upper reach (95.9 percent) and least in the lower reach (42.0 percent). All forests are expected to be 38.2 percent drier in species composition by 2085, the year when the median age of surviving 2004 subcanopy trees will reach the median age (99 years) of the 2004 large canopy trees. The change will be greatest for forests in the upper reach (45.0 percent). Forest composition changes from pre-1954 to 2085 were calculated using Floodplain Indices from 1976 and 2004 tree-size classes and replicate plots. Species composition in high bottomland hardwood forests is expected to continue to change, and some low bottomland hardwood forests are expected to become high bottomland hardwood forests. Organisms associated with floodplain forests will be affected by the changes in tree species, which will alter the timing of leaf-out, fruiting, and leaf-drop, the types of fruit and debris produced, and soil chemistry. Swamps will contain more bottomland hardwood species, but will also have an overall loss of tree density. The density of trees in swamps significantly decreased by 37 percent from 1976 to 2004. Of the estimated 4.3 million (17 percent) fewer trees that existed in the nontidal floodplain in 2004 than in 1976, 3.3 million trees belonged to four swamp species: popash, Ogeechee tupelo, water tupelo, and bald cypress. Water tupelo, the most important tree in the nontidal floodplain in terms of basal area and density, has declined in number of trees by nearly 20 percent since 1976. Ogeechee tupelo, the species valuable to the tupelo honey industry, has declined in number of trees by at least 44 percent. Greater hydrologic variability in recent years may be the reason swamps have had a large decrease in tree density. Drier conditions are detrimental for the growth of swamp species, and periodic large floods kill invading bottomland hardwood trees. The loss of canopy density in swamps may result in the swamp floor being exposed to more light with an increase in the amount of ground cover present, which in turn, would reduce tree replacement. The microclimate of the swamp floor would become warmer due to the decrease in shade and inundation. Soils would become dehydrated more quickly in dry periods and debris would decompose more quickly. A loss of tree density in swamps would lead to a decrease in tree and leaf litter biomass, which would have additional effects on swamp organisms. The loss of litter would result in a loss of substrate for benthic organisms in the floodplain and, ultimately, in the downstream waters of the river and estuary.

Florida↗

USGS2015

Updated numerical model with uncertainty assessment of 1950-56 drought conditions on brackish-water movement within the Edwards aquifer, San Antonio, Texas

In 2010, the U.S. Geological Survey, in cooperation with the San Antonio Water System, began a study to assess the brackish-water movement within the Edwards aquifer (more specifically the potential for brackish-water encroachment into wells near the interface between the freshwater and brackish-water transition zones, referred to in this report as the transition-zone interface) and effects on spring discharge at Comal and San Marcos Springs under drought conditions using a numerical model. The quantitative targets of this study are to predict the effects of higher-than-average groundwater withdrawals from wells and drought-of-record rainfall conditions of 1950–56 on (1) dissolved-solids concentration changes at production wells near the transition-zone interface, (2) total spring discharge at Comal and San Marcos Springs, and (3) the groundwater head (head) at Bexar County index well J-17. The predictions of interest, and the parameters implemented into the model, were evaluated to quantify their uncertainty so the results of the predictions could be presented in terms of a 95-percent credible interval. The model area covers the San Antonio and Barton Springs segments of the Edwards aquifer; the history-matching effort was focused on the San Antonio segment. A previously developed diffuse-flow model of the Edwards aquifer, which forms the basis for the model in this assessment, is primarily based on a conceptualization in which flow in the aquifer is predominately through a network of numerous small fractures and openings. Primary updates to this model include an extension of the active area downdip, a conversion to an 8-layer SEAWAT variable-density flow and transport model to simulate dissolved-solids concentration effects on water density, history matching to 1999–2009 conditions, and parameter estimation in a highly parameterized context using automated methods in PEST (a model-independent Parameter ESTimation code). In addition to the best-fit parameter values derived from history matching, the uncertainty of model parameters was also estimated by using linear uncertainty analysis. Comparison of “prior” (before history matching) and “posterior” (after history matching) variances of parameters indicate that the information within the observation dataset used for history matching informs many parameters. The concentration threshold parameters were well-informed by the observation dataset as their posterior distributions were much narrower than their prior distributions. The transition-zone scaling parameters of hydraulic conductivity, effective porosity, and specific storage were all informed by the observation dataset, as evidenced by the difference between the prior and posterior variances. Saline-zone scaling parameters, alternatively, were not informed by the observation dataset for effective porosity and specific storage. Resulting posterior drier-month, wetter-month, and annual recharge multiplier parameter variances are important to understanding how well recharge is estimated and implemented within the model. The shifts of the posterior distributions left and right indicate that there were zones where less or more water was needed in the model. The widths of the distributions were not decreased substantially, indicating that many of the best-fit recharge parameters are not statistically different from the initial values specified in the history-matching effort. Recharge from rainfall is the driving force behind groundwater flow and heads in the aquifer; therefore, an increase in understanding of this process would benefit model development by potentially decreasing the uncertainty of this parameter. The history-matching effort was most helpful in informing the parameters in the model that control discharge at springs, namely, the spring orifice (drain) altitude and drain conductance parameters for each spring. The uncertainty assessment of the predictive model (a hypothetical recurrence of 1950–56 drought conditions and higher-than-average groundwater withdrawals from wells) provided insights into the potential effects of these conditions on dissolved-solids concentration changes at production wells near the transition-zone interface, discharges at Comal and San Marcos Springs, and heads at Bexar County index well J-17. Results at the 25 production wells near the transition-zone interface indicate that the uncertainty of model input parameters based on expert knowledge yielded an upper bound of the 95-percent credible interval of dissolved-solids concentrations that exceeds the secondary drinking water standards of 1,000 milligrams per liter (mg/L) of the Texas Commission on Environmental Quality (TCEQ) for many wells. However, the history-matching process provided key information to inform prediction-sensitive model parameters and therefore, contributed to a substantial decrease of the upper bound of the 95-percent credible interval to below the secondary drinking water standards. Reductions in dissolved-solids concentration changes were on the order of 400 mg/L to 1,300 mg/L. The reduction in uncertainty in regards to this prediction implies that this prediction of dissolved-solids concentration change can be made with some certainty using this current model and that those parameters that control this prediction are informed by the observation dataset. Even though predictive uncertainty was reduced for this prediction, dissolved-solids concentration changes were still greater than zero, indicating a minimal increase in concentration at these 25 production wells during the 7-year simulation period is likely. However, this minimal concentration increase indicates a small potential for movement of the brackish-water transition zone near these wells during the 7-year simulation period of drought-ofrecord (1950–56) rainfall conditions with higher-than-average groundwater withdrawals by wells. Predictive results of total spring discharge during the 7-year period, as well as head predictions at Bexar County index well J-17, were much different than the dissolved-solids concentration change results at the production wells. These upper bounds are an order of magnitude larger than the actual prediction which implies that (1) the predictions of total spring discharge at Comal and San Marcos Springs and head at Bexar County index well J-17 made with this model are not reliable, and (2) parameters that control these predictions are not informed well by the observation dataset during historymatching, even though the history-matching process yielded parameters to reproduce spring discharges and heads at these locations during the history-matching period. Furthermore, because spring discharges at these two springs and heads at Bexar County index well J-17 represent more of a cumulative effect of upstream conditions over a larger distance (and longer time), many more parameters (with their own uncertainties) are potentially controlling these predictions than the prediction of dissolved-solids concentration change at the prediction wells, and therefore contributing to a large posterior uncertainty.

Texas↗

USGS2015

California State Waters Map Series: Offshore of Refugio Beach, California

In 2007, the California Ocean Protection Council initiated the California Seafloor Mapping Program (CSMP), designed to create a comprehensive seafloor map of high-resolution bathymetry, marine benthic habitats, and geology within California’s State Waters. The CSMP approach is to create highly detailed seafloor maps through collection, integration, interpretation, and visualization of swath sonar data, acoustic backscatter, seafloor video, seafloor photography, high-resolution seismic-reflection profiles, and bottom-sediment sampling data. The map products display seafloor morphology and character, identify potential marine benthic habitats, and illustrate both the surficial seafloor geology and shallow (to about 100 m) subsurface geology. The Offshore of Refugio Beach map lies within the western Santa Barbara Channel region of the Southern California Bight. This geologically complex region forms a major biogeographic transition zone, separating the cold-temperate Oregonian province north of Point Conception from the warm-temperate California province to the south. The map area is in the southern part of the Western Transverse Ranges geologic province, which is north of the California Continental Borderland. Significant clockwise rotation—at least 90°—since the early Miocene has been proposed for the Western Transverse Ranges province, and geodetic studies indicate that the region is presently undergoing north-south shortening. Uplift rates (as much as 0.5 mm/yr) that are based on studies of onland marine terraces provide further evidence of significant shortening. The coastal zone of this map area lies at the steep flank of the Santa Ynez Mountains. The crest of the range, which lies about 8 km from the shoreline (north of the map area), has a maximum elevation of about 780 m. This area is largely open space, partly used for livestock grazing, with no significant towns or population centers. Highway 101 crosses the map area, adjacent to and within a few hundred meters of the shoreline. The most significant developments are the recreational state beaches at El Capitan Beach and Refugio Beach. The beaches are subject to erosion each winter during storm-wave attack, and then they undergo gradual recovery or accretion during the more gentle wave climate of the late spring, summer, and fall months. The Offshore of Refugio Beach map area lies in the west-central part of the Santa Barbara littoral cell, which is characterized by west-to-east transport of sediment from Point Arguello on the northwest to Hueneme and Mugu Canyons on the southeast. Longshore drift rates have been reported to range from about 160,000 to 800,000 tons/yr, averaging 400,000 tons/yr. Sediment supply to the western and central part of the littoral cell, including the map area, is mainly from relatively small coastal watersheds. Within the map area, these coastal watersheds include (from east to west) Cañada del Capitan, Tajiquas Creek, Arroyo Hondo, Cañada del Molino, and several unnamed canyons and creeks. The Santa Ynez and Santa Maria Rivers, the mouths of which are 80 to 120 km northwest of the map area, are not considered significant sediment sources because Point Conception and Point Arguello provide obstacles to downcoast sediment transport and also because much of their sediment load is trapped in dams. The Ventura and Santa Clara Rivers, the mouths of which are about 70 km to the southeast of the map area, are not sediment sources for the map area. The offshore part of the map area consists of relatively flat and shallow continental shelf, which dips gently seaward (about 0.8° to 1.0°) so that water depths at the shelf break, roughly coincident with the California’s State Waters limit, are about 80 to 100 m. This part of the Santa Barbara Channel is relatively well protected from large Pacific swells from the north and northwest by Point Conception and from the south and southwest by offshore islands and banks. The shelf is underlain by variable amounts of upper Quaternary marine and fluvial sediments deposited as sea level fluctuated in the late Pleistocene. In the map area, the shelf break is at depths of about 90 m and lies about 5.6 to 6.4 km offshore. Beyond the shelf break, the slope is steep (as much as about 7°) and unstable. Several submarine landslides, including the large (130 km 2 ) Goleta landslide complex, have been documented offshore of Goleta, a few kilometers east of the map area. This compound slump complex may have been initiated more than 200,000 year ago, but it also includes three recent failures that may have been generated 8,000 to 10,000 years ago. A local, 10-m-high tsunami may have been generated from these failure events. Small folds related to local faulting are superimposed on the homocline that makes up the south flank of the Santa Ynez Mountains. One of these superimposed anticlines hosts the Molino gas field, which was discovered in 1962 and subsequently developed through directional drilling from onshore wells. The Oligocene Sespe and Vaqueros Formations are the reservoirs in the Molino gas field, and the map area includes numerous seafloor hydrocarbon seeps. Seafloor habitats in the broad Santa Barbara Channel region consist of significant amounts of soft, unconsolidated sediment interspersed with isolated areas of rocky habitat that support kelp-forest communities nearshore and rocky-reef communities in deep water. The potential marine benthic habitat types mapped in the Offshore of Refugio Beach map area are directly related to its Quaternary geologic history, geomorphology, and active sedimentary processes. These potential habitats, which lie primarily within the Shelf (continental shelf) but also partly within the Flank (basin flank or continental slope) megahabitats, primarily are composed of soft sediment interrupted by a few carbonate mounds. This homogeneous seafloor of sediment and low-relief bedrock provides promising habitat for groundfish, crabs, shrimp, and other marine benthic organisms.

California↗

USGS2017

Water quality and quantity and simulated surface-water and groundwater flow in the Laurel Hill Creek Basin, southwestern Pennsylvania, 1991–2007

Laurel Hill Creek is considered one of the most pristine waterways in southwestern Pennsylvania and has high recreational value as a high-quality cold-water fishery; however, the upper parts of the basin have documented water-quality impairments. Groundwater and surface water are withdrawn for public water supply and the basin has been identified as a Critical Water Planning Area (CWPA) under the State Water Plan. The U.S. Geological Survey, in cooperation with the Somerset County Conservation District, collected data and developed modeling tools to support the assessment of water-quality and water-quantity issues for a basin designated as a CWPA. Streams, springs, and groundwater wells were sampled for water quality in 2007. Streamflows were measured concurrent with water-quality sampling at main-stem sites on Laurel Hill Creek and tributaries in 2007. Stream temperatures were monitored continuously at five main-stem sites from 2007 to 2010. Water usage in the basin was summarized for 2003 and 2009 and a Water-Analysis Screening Tool (WAST) developed for the Pennsylvania State Water Plan was implemented to determine whether the water use in the basin exceeded the “safe yield” or “ the amount of water that can be withdrawn from a water resource over a period of time without impairing the long-term utility of a water resource .” A groundwater and surface-water flow (GSFLOW) model was developed for Laurel Hill Creek and calibrated to the measured daily streamflow from 1991 to 2007 for the streamflow-gaging station near the outlet of the basin at Ursina, Pa. The CWPA designation requires an assessment of current and future water use. The calibrated GSFLOW model can be used to assess the hydrologic effects of future changes in water use and land use in the basin. Analyses of samples collected for surface-water quality during base-flow conditions indicate that the highest nutrient concentrations in the main stem of Laurel Hill Creek were at sites in the northeastern part of the basin where agricultural activity is prominent. All of the total nitrogen (N) and a majority of the total phosphorus (P) concentrations in the main stem exceeded regional nutrient criteria levels of 0.31 and 0.01 milligrams per liter (mg/L), respectively. The highest total N and total P concentrations in the main stem were 1.42 and 0.06 mg/L, respectively. Tributary sites with the highest nutrient concentrations are in subbasins where treated wastewater is discharged, such as Kooser Run and Lost Creek. The highest total N and total P concentrations in subbasins were 3.45 and 0.11 mg/L, respectively. Dissolved chloride and sodium concentrations were highest in the upper part of the basin downstream from Interstate 76 because of road deicing salts. The mean base-flow concentrations of dissolved chloride and sodium were 117 and 77 mg/L, respectively, in samples from the main stem just below Interstate 76, and the mean concentrations in Clear Run were 210 and 118 mg/L, compared to concentrations less than 15 mg/L in tributaries that were not affected by highway runoff. Water quality in forested tributary subbasins underlain by the Allegheny and Pottsville Formations was influenced by acidic precipitation and, to a lesser extent, the underlying geology as indicated by pH values less than 5.0 and corresponding specific conductance ranging from 26 to 288 microsiemens per centimeter at 25 degrees Celsius for some samples; in contrast, pH values for main stem sites ranged from 6.6 to 8.5. Manganese (Mn) was the only dissolved constituent in the surface-water samples that exceeded the secondary maximum contaminant level (SMCL). More than one-half the samples from the main stem had Mn concentrations exceeding the SMCL level of 50 micrograms per liter (μg/L), whereas only 19 percent of samples from tributaries exceeded the SMCL for Mn. Stream temperatures along the main stem of Laurel Hill Creek became higher moving downstream. During the summer months of June through August, the daily mean temperatures at the five sites exceeded the limit of 18.9 degrees Celsius (°C) for a cold-water fishery. The maximum instantaneous values for each site ranged from 27.2 to 32.8 °C. Water-quality samples collected at groundwater sites (wells and springs) indicate that wells developed within the Mauch Chunk Formation had the best water quality, whereas wells developed within the Allegheny and Pottsville Formations yielded the poorest water quality. Waters from the Mauch Chunk Formation had the highest median pH (7.6) and alkalinity (80 mg/L calcium carbonate) values. The lowest pH and alkalinity median values were in waters from the Allegheny and Pottsville Formations. Groundwater samples collected from wells in the Allegheny and Pottsville Formations also had the highest concentrations of dissolved iron (Fe) and dissolved Mn. Seventy-eight percent of the groundwater samples collected from the Allegheny Formation exceeded the SMCL of 300 μg/L for Fe and 50 μg/L for Mn. Forty-three and 62 percent of the groundwater samples collected from the Pottsville Formation exceeded the SMCL for iron and Mn, respectively. The highest Fe and Mn concentrations for surface waters were measured for tributaries draining the Pottsville Formation. The highest median Fe concentration for tributaries was in samples from streams draining the Allegheny Formation. During base-flow conditions, the streamflow per unit area along the main stem of Laurel Hill Creek was lowest in the upper parts of the basin [farthest upstream site 0.07 cubic foot per second per square mile (ft 3 /s/mi 2 )] and highest (two sites averaging about 0.20 (ft 3 /s/mi 2 ) immediately downstream from Laurel Hill Lake in the center of the basin. Tributaries with the highest streamflow per unit area were those subbasins that drain the western ridge of the Laurel Hill Creek Basin. The mean streamflow per unit area for tributaries draining areas that extend into the western ridge and draining eastern or central sections was 0.24 and 0.05 ft 3 /s/mi 2 , respectively. In general, as the drainage area increased for tributary basins, the streamflow per unit area increased. Criteria established by the Pennsylvania Department of Environmental Protection indicate that the safe yield of water withdrawals from the Laurel Hill Creek Basin is 1.43 million gallons per day (Mgal/d). Water-use data for 2009 indicate that net (water withdrawals subtracted by water discharges) water withdrawals from groundwater and surface-water sources in the basin were approximately 1.93 Mgal/d. Water withdrawals were concentrated in the upper part of the basin with approximately 80 percent of the withdrawals occurring in the upper 36 mi 2 of the basin. Three subbasins—Allen Creek, Kooser Run, and Shafer Run— in the upper part were affected the most by water withdrawals such that safe yields were exceeded by more than 1,000 percent in the first two and more than 500 percent in the other. In the subbasin of Shafer Run, intermittent streamflow characterizes sections that historically have been perennial. The GSFLOW model of the Laurel Hill Creek Basin is a simple one-layer representation of the groundwater flow system. The GSFLOW model was primarily calibrated to reduce the error term associated with base-flow periods. The total amount of observed streamflow at the Laurel Hill Creek at Ursina, Pa. streamflow-gaging station and the simulated streamflow were within 0.1 percent over the entire modeled period; however, annual differences between simulated and observed streamflow showed a range of -27 to 24 percent from 1992 to 2007 with nine of the years having less than a 10-percent difference. The primary source of simulated streamflow in the GSFLOW model was the subsurface (interflow; 62 percent), followed by groundwater (25 percent) and surface runoff (13 percent). Most of the simulated subsurface flow that reached the stream was in the form of slow flow as opposed to preferential (fast) interflow.

Pennsylvania↗

USGS2019

Geologic map of the central-southeast flank of Mauna Loa Volcano, Island of Hawaii, Hawaii

Mauna Loa, the largest volcano on Earth, has erupted 33 times since written descriptions became available in 1832. Some eruptions began with only brief seismic unrest, while others followed several months to a year of increased seismicity. Once underway, its eruptions can produce lava flows that may reach the sea in less than 24 hours, severing roads and utilities. For example, lava flows erupted from the Southwest Rift Zone (SWRZ) in 1950 advanced at an average rate of 9.3 km per hour, and all three lobes reached the ocean within approximately 24 hours (Finch and Macdonald, 1953). Near the eruptive vents, the flows must have traveled even faster. In terms of eruption frequency, pre-eruption warning, and rapid flow emplacement, Mauna Loa poses an enormous volcanic-hazard threat to the Island of Hawai‘i. Volcanic hazards on Mauna Loa may be anticipated, and risk substantially mitigated, by documenting the past activity to refine our knowledge of the hazards and by alerting the public and local government officials of our findings and their implications for hazards assessments and risk. From the geologic record, we may deduce several generalized facts about the geologic history of the Northeast Rift Zone (NERZ). The middle to uppermost segments of the rift zone were more active in the past 4,000 years than the lower portion of the rift zone. This may be due to buttressing of the lower east rift zone by Mauna Kea and Kīlauea volcanoes. The historical flows that erupted on the north side of the rift zone advanced toward Hilo. This flank of the volcano may be more vulnerable to inundation. Lockwood (1990) noted that the vents of historical activity are migrating to the south. The volcano appears to have a self-regulating mechanism that evenly distributes long-term activity across its flanks. The geologic record also supports this notion; the time prior to the historical period (Age Group 1, pre-A.D. 1832 to 1,000 yrs B.P.; orange units) is dominated by activity on the south side of the NERZ. Although most Mauna Loa eruptions begin in the summit area at the 12,000-ft elevation (Lockwood and Lipman, 1987), the central-southeast flank has not been the source of any activity. All flows originated from the summit or the upper reaches of the Northeast Rift Zone (NERZ) or the Southwest Rift Zone (SWRZ). The NERZ was the source of eight flank eruptions since 1843. The NERZ extends from the 13,680-ft-high summit towards Hilo (population ~60,000; second-largest city in State of Hawaii). The northern portion of the map area is built entirely on flows erupted from the NERZ. The SWRZ extends from the summit towards Kalae (South Point) at sea level. The southern portion of the map area is built entirely on flows erupted from the SWRZ. The map area extends from the 10,350 -ft elevation on Mauna Loa’s east flank toward the Hawaii Volcanoes National Park and the town of Volcano (population approx. 2,000) in the northeast. At the south boundary of the map area is the town of Pāhala (population approx. 900). This map includes areas adjacent to and downslope of the NERZ and regions east of and directly downslope of Moku‘āweoweo, Mauna Loa’s summit caldera. The map encompasses 506 km 2 of the southeast flank (fig. 1) of Mauna Loa from 10,350 -ft elevation to sea level. The map of the central-southeast flank of Mauna Loa shows the distribution and relations of volcanic and surficial sedimentary deposits separated into 15 age groups ranging from a period greater than 50,000 yr B.P. to A.D. 1984. It incorporates previously reported work published in generalized small-scale maps (Lockwood and Lipman, 1987; Lockwood, 1995; Wolfe and Morris, 1996). This map is the second in a series of five maps that will cover Mauna Loa volcano. See SIM 2932-A at https://doi.org/10.3133/sim2932A . NOTE: Map sheet 1 contains lines and type with overprint. This feature may be turned on or off in the Adobe Acrobat page display preferences.

Hawaii↗

USGS2004

Hydrology and cycling of nitrogen and phosphorus in Little Bean Marsh: A remnant riparian wetland along the Missouri River in Platte County, Missouri, 1996–97

The lack of concurrent water-quality and hydrologic data on riparian wetlands in the Midwestern United States has resulted in a lack of knowledge about the water-quality functions that these wetlands provide. Therefore, Little Bean Marsh, a remnant riparian wetland along the Missouri River, was investigated in 1996 and 1997 primarily to determine the magnitude and character of selected water-quality benefits that can be produced in such a wetland and to identify critical processes that can be managed in remnant or restored riparian wetlands for amelioration of water quality. Little Bean Marsh averages 69 hectares in size, has a maximum depth of about 1 meter, and the majority of the marsh is covered by macrophytes. In 1997, 41 percent of the water received by Little Bean Marsh was from direct precipitation, 14 percent was from ground-water seepage, 30 percent from watershed runoff, and 15 percent was backflow from Bean Lake. Although, Little Bean Marsh was both a ground-water recharge and discharge area, discharge to the marsh was three times the recharge to ground water. Ground-water levels closely tracked marsh water levels indicating a strong hydraulic connection between ground water and the marsh. Reduced surface runoff and ground-water availability are stabilizing influences on marsh hydrology and probably contribute to the persistence of emergent vegetation. The rapid hydraulic connection between Little Bean Marsh and ground water indicates that the hydrologic regime of most wetlands along the lower Missouri River is largely a function of the altitude of the marsh bottom relative to the altitude of the water table. More water was lost from the marsh through evapotranspiration (59 percent) than all other pathways combined. This is partially because the transpiration process of abundant macrophytes can greatly contribute to the evapotranspiration above that lost from open water surfaces. Surface outflow accounted for 36 percent and ground-water seepage accounted for only 5 percent of the losses. Large residence times allows the marsh to greatly affect water quality before water escapes as ground-water recharge or surface outflow. The shallowness of Little Bean Marsh and ion exclusion during ice formation caused the highest specific conductances of 1,100 to 1,300 microsiemens per centimeter at 25 degrees Celsius to occur during the winter. This concentration of dissolved solutes under ice can make wetlands more vulnerable to toxic contaminants than deeper surface-water bodies. Dissolved oxygen was less than 5 mg/L (milligrams per liter) for 3 to 4 months and near 0 mg/L for about 1 month in summer. Despite depths of less than 1 meter, temperature stratification persisted more than 3 months during the summers of 1996 and 1997, preventing mixing and contributing to periods of anoxia. Shallow depths and extended periods of anoxia in the marsh limit the ability of some organisms to escape high-temperature stress. Turbidity in Little Bean Marsh usually was low for several reasons: sediment loadings from the largely flood-plain drainage were low, emergent vegetation shade out algae and shield the water from wind, and high concentrations of bivalent cations increase flocculation rates of inorganic suspended material. The high concentrations of bivalent cations was largely because of a substantial amount of ground-water seepage into the marsh. Dissolved organic nitrogen was the dominant nitrogen species in Little Bean Marsh. Denitrification and biotic uptake kept more than 62 percent of nitrate (NO3) and 43 percent of ammonium (NH4) concentrations in marsh samples less than a detection limit of 0.005 mg/L. This contrasts with the Missouri River where inorganic NO3 dominates. Consequently, artificial flood-plain drainage that bypasses riparian wetlands likely deliver substantially more biotically available inorganic nitrogen to receiving waters than surface water that has been routed through a remnant wetland. Average total nitrogen concentrations in Little Bean Marsh were substantially less than those at other Missouri River wetlands, roughly one-half the mean concentrations in the Missouri River, but roughly twice the average nitrogen values in reservoirs of the glaciated plains of Missouri. The largest concentrations of nearly all species of nitrogen and phosphorus and the most intense period of hypereutrophy coincided with a phytoplankton bloom and senescence of River Bulrush (Scirpus fluviatilis) and common cattail (Typha latifolia) in September 1997. The rapid leaching of nitrogen that occurs soon after macrophyte senescence combined with a recent destratification of the marsh probably provided nitrogen to the nitrogen-limited open-water areas and triggered a phytoplankton bloom. Despite the rarity of runoff events, surface runoff from the watershed, combined with atmospheric deposition, contributed more than seven times the 530 kg (kilograms) of nitrogen that escaped Little Bean Marsh in surface outflow during 1997. Atmospheric deposition alone was more than 530 kg. Seepage to ground water contained less than 1.5 percent of the nitrogen leaving the marsh in surface outflow. The slow decay rate of Scirpus fluviatilis and reducing conditions in bottom sediments make burial of organic nitrogen a substantial sink of nitrogen. Denitrification experiments indicate that denitrification rates were limited by NO3 in the water column. Consequently, decomposition and nitrification of NH4 and organic nitrogen are the rate limiting steps of nitrogen removal in Little Bean Marsh. The NO3-limited rates of denitrification also indicate that Little Bean Marsh has a large unused capacity for nitrogen removal. These data indicate that the vast extent of riparian marshes along the Missouri and Mississippi Rivers may have had a substantial role in limiting NO3 loads to the Gulf of Mexico before agricultural development of flood plains. Drainage and removal of riparian marshes may be a major cause of the increased NO3 loads to the Gulf of Mexico. Periods of anoxia had much larger effects on phosphorus release than the other variables. The largest concentrations of phosphorus occurred in late summer and corresponded with senescing macrophytes, periods of anoxia, and a large algal bloom in Little Bean Marsh. Low water levels prevented the escape of phosphorus in surface outflow during these periods of highest phosphorus concentrations. Dry weather in late summer is typical and probably makes the correspondence of low water levels, anoxia, and consequent low phosphorus release a common occurrence in marshes along the Missouri River. Little Bean Marsh retained more than 95 percent of the phosphorus it received. The amount of phosphorus in surface inflows to the marsh were more than one order of magnitude greater than that escaping in surface outflows. The long hydraulic residence time of the marsh and large contributions of iron from ground water (that provide many sorption sites for phosphorus) make the marsh an effective sediment and phosphorus trap.

Missouri↗

USGS2025

Assessment of channel morphology, hydraulics, and bedload transport along the Siletz River, western Oregon

Significant Findings Chinook salmon ( Oncorhynchus tshawytscha ) and Pacific lamprey ( Entosphenus tridentatus ) are native, anadromous fish species in the Siletz River Basin, western Oregon, that face many threats to their survival in freshwater and the ocean. The Confederated Tribes of Siletz Indians of Oregon seek to mitigate freshwater threats to Chinook salmon and Pacific lamprey, where possible, with habitat conservation and restoration efforts. This study was conducted to assist the Confederated Tribes of Siletz Indians of Oregon in documenting and understanding the hydrogeomorphic processes shaping present-day habitat conditions and assessing future habitat implications for Chinook salmon and Pacific lamprey along the main-stem Siletz River. As such, this study focused on understanding geomorphic processes and patterns of channel change, including lateral and vertical adjustments in channel position and changes in bed-material sediment (sands, gravels, and cobbles that mantle the channel bed), which collectively determine overall patterns of channel morphology and fluvial habitats. Objective One was to evaluate lateral changes in channel position, vertical changes in bed elevation, and longitudinal patterns in bed-material particle size along the Siletz River using detailed channel maps developed from aerial photographs collected from 1939 to 2016, long-term records of stage and discharge collected by the U.S. Geological Survey (USGS) near the City of Siletz, and sediment particle size data. Objective Two was to assess hydraulic conditions using one- and two-dimensional hydraulic models and transport capacity of bed-material sediment using bedload transport models and sediment particle size data for a range of discharge conditions. Objective Three was to identify potential burrowing habitat for lamprey larvae (PBH) along the Siletz River network and provide insights in local factors influencing PBH along the main-stem Siletz River. The overall findings are synthesized to describe habitat implications for Chinook salmon and Pacific lamprey under present-day and future conditions. Results of Objective One, an evaluation of changes in channel position and bed elevations and longitudinal patterns in bed-material particle size along the Siletz River, include the following From 1939 to 2016, channel planform and the mapped area of gravel bars did not change considerably along the 97.2-kilometers (km) of Siletz River between Elk Creek and Millport Slough, except for in short sections generally less than 1 km long. This inherent lateral channel stability results from the resistant bedrock and terraces that bound most of the channel and limit lateral changes in channel position. Intermittent sections along the study reach where the active channel widened at channel bends displayed noticeable planform changes and increases in mapped bar area. From water year (WY; a 12-month period from October 1 through September 30 and named for the year in which it ends) 1906 to 2021, changes in the stage-discharge relation interpreted as rising and lowering channel bed levels were observed at the USGS streamgage 14305500 (Siletz River at Siletz, OR) in response to floods (such as high flows in February 1996 and January 2002 that exceeded 0.1 and 0.667 annual exceedance probabilities [AEP] events, respectively). However, the rating curve representing the stage-discharge relation did not change in response to high-magnitude floods between 2007 and 2021. Along the approximately 54-km of the Siletz River between Moonshine Park and the Bulls Bag area, surficial particle distributions varied considerably between sampling sites in response to changes in channel width and gradient, sediment inputs from tributaries, and basin geology. Despite this variability, median particle sizes tended to decrease in size in the downstream direction over the 54 km. Results of Objective Two, an evaluation of hydraulic and bedload transport conditions along the Siletz River, include the following The most substantial increases in maximum and mean water velocity and bedload sediment transport capacity occurred at events between the 0.995 and 0.50 AEPs. Events of these magnitudes occur approximately every 1–2 years. Smaller events (0.995 AEP) are generally contained by the banks of the main channel, whereas larger events (0.50 AEP) generally spill over the banks and inundate high-elevation bar and low-elevation floodplain surfaces. Multiple smaller floods within a WY that exceed the 0.995 and 0.50 AEPs (such as in WY 1996) can transport as much or more bed-material sediment than a single, higher magnitude event (such as the maximum event recorded in WY 2000 with an estimated AEP of around 0.002). Bedload transport capacity generally exceeds sediment supply (greater than 2 millimeters [mm]) for most of the study area from Wildcat Creek to the City of Siletz as evidenced by substantial in-channel bedrock. Despite overall conditions of limited sediment supply, transport capacity still varies considerably within and between years with discharge magnitude and spatially in relation to local hydraulics imposed by bedrock, channel morphology, and human infrastructure (such as bridges). Results of Objective Three, an analysis of PBH for lamprey larvae, include the following About 28 percent of the Siletz River network meets the mean annual suspended sediment loads and channel slope criteria for PBH. Along the main-stem Siletz River, in-channel bedrock outcrops and high transport capacity are expected to further constrain PBH. Together, these results suggest that most of the Siletz River between Wildcat Creek and the City of Siletz has had only modest vertical and lateral change between the 1930s and 2010s because of the bedrock in and along the main channel and the river’s relatively high transport capacity relative to bed-material sediment supply. However, localized sections of the Siletz River where the active channel widens, particularly at channel bends, exhibited some change in channel planform and the locations and area of gravel bars. In the future, moderate increases in autumn-winter discharge may not result in substantial changes in coarse gravel bars along the Siletz River but may result in selective transport of finer bed-material sediment (gravel, sands, and silts) that provide spawning habitats for Chinook salmon and Pacific lamprey and burrowing habitats for lamprey larvae. Assuming no substantial changes in bed-material sediment supply, increased bedload transport capacity may cause frequent entrainment of lamprey larvae that are burrowed in coarse sand deposits, suspension and downstream transport of salmon eggs incubating in gravels, and reductions in the areas of spawning gravels for Chinook salmon and Pacific lamprey. Exact implications of current and future discharge conditions for these species along the Siletz River depends on many factors, including sediment supply, local hydraulics, and the timing of flood events relative to fish life stages.

Oregon↗

USGS2020

Spatial and temporal patterns in streamflow, water chemistry, and aquatic macroinvertebrates of selected streams in Fairfax County, Virginia, 2007–18

Urbanization substantially alters the landscape in ways that can impact stream hydrology, water chemistry, and the health of aquatic communities. Stormwater best management practices (BMPs) are the primary tools used to mitigate the effects of urban stressors such as increased runoff, decreased baseflow, and increased nutrient and sediment transport. To date, Fairfax County Virginia’s stormwater management program has made substantial investments into the implementation of both structural and nonstructural BMPs aimed at restoring and protecting watersheds. The U.S. Geological Survey (USGS), in cooperation with Fairfax County, Virginia, established a long-term water-resources monitoring program to evaluate the watershed-scale effects of these investments. Monitoring began at 14 stations in 2007 and was expanded to 20 stations in 2013. This report utilized the first 10 years of data collection to (1) assess water quantity and quality, as well as ecological condition; (2) compute annual nutrient and sediment loads; and (3) evaluate trends in streamflow, water quality, and ecological condition. Efforts are underway to link the biotic and abiotic patterns described herein to watershed management practices as well as factors such as land use change, public works infrastructure, and climate. Hydrologic, chemical, and benthic macroinvertebrate community conditions in the streams monitored were similar to those observed in other studies of urban streams. Multidecadal trends in baseflow indices and runoff ratios at long-term Chesapeake Bay Non-tidal Network streamgages (CB-NTN) indicate a decrease in groundwater recharge and increase in storm runoff as a result of urbanization. Streamflow yields varied spatially with land cover, geology, and soil characteristics, whereas flashiness was positively related to impervious area. Dissolved oxygen typically was lowest in the Coastal Plain and across all Triassic Lowlands streams, and highest in the Piedmont. Dissolved oxygen concentrations generally were above Virginia’s minimum criterion of 4.0 milligrams per liter (mg/L), most violations occurred at Paul Spring Branch in the Coastal Plain during the warmest months of the year owing to increased chemical and biological oxygen demand. Typical pH values of the monitored streams centered on neutrality (pH = 7); however, diurnal fluctuations were most prevalent in the continuous pH data at Flatlick Branch (FLAT; a Triassic Lowlands station), as a result of increased photosynthesis catalyzed by phosphorus-rich geology. Specific conductance (SC) varied spatially owing to geology (highest at Triassic Lowlands stations) and anthropogenic disturbance (watersheds with high impervious land cover). Specific conductance typically was inversely related to streamflow except in winter months following deicing road salt applications, when values increased by several orders of magnitude. A significant increase in SC of about 2 percent per year was observed from the combined trend result of all monitoring stations over the 10-year period. Significant SC increases occurred at nearly all monitoring stations. Increasing trends were observed during winter and nonwinter months, which suggests that salts applied to deice roadways and other impervious surfaces are stored in the environment and released year-round. Suspended-sediment (SS) concentrations in monthly samples did not vary significantly between most stations, but typically were highest in the spring and lowest in the fall as a result of seasonal differences in streamflow and climate. Suspended-sediment yields ranged from 62 to 1,428 tons per square mile (ton/mi 2 ), with a median of 302 ton/mi 2 . Annual loads were greatest during the wettest water years (October 1-September 30; 2008, 2011, and 2014), with the greatest interannual variability occurring at Difficult Run above Fox Lake (DIFF) and South Fork Little Difficult Run (SFLIL). Suspended sediment was primarily composed of silts and clays; however, the proportion of sand in suspended sediment was related positively to streamflow. Cross-correlation analyses suggested the dominant sources of SS were streambank erosion and resuspension of in-channel material at DIFF and FLAT; whereas, upland sources and erosion of upper streambanks were more common at Dead Run (DEAD), Long Branch (LONG), and SFLIL. Median total phosphorus (TP) concentrations ranged from 0.016 to 0.077 mg/L, with a networkwide median of 0.022 mg/L, were highest in the warm season (April-September), and were composed primarily of dissolved phosphorous. Although TP concentrations were relatively low across the network, the highest concentrations were consistently at stations located in the Triassic Lowlands, owing to phosphorous-rich geology, and in the Coastal Plain, owing to the low-phosphorous sorptive capacity of those soils. A significant increase in TP concentration occurred in a few stations, but the combined trend results from all stations demonstrated a significant increase of about 4 percent per year. Networkwide increases were also observed in total dissolved phosphorus, orthophosphate, and total particulate phosphorus. The composition of TP shifted from dissolved to particulate as streamflow increased and for this reason loads primarily were composed of particulate phosphorous. Median annual TP loads were highest at FLAT and DEAD and ranged from 247 to 642 pounds per square mile (lbs/mi 2 ) networkwide. Interannual variability in phosphorous yields was apparent at most stations; the highest loading years were also the wettest years during the study period and coincident with the highest peak annual flows. Total nitrogen (TN) concentrations typically were low throughout the network with exceptions occurring at stations located in watersheds with a high density of septic infrastructure. Elevated TN concentrations also were observed in some watersheds without a high density of septic systems and may be attributable to geologic and soil properties that limit denitrification as well as other unknown anthropogenic inputs. Total nitrogen typically was dominated by nitrate during baseflows; however, the proportion of particulate nitrogen increased during stormflows. Total nitrogen yields were similar across stations, with medians ranging from about 3,600 to 6,300 lbs/mi 2 and were related to annual streamflow volume. Total nitrogen concentrations and flow-normalized concentrations decreased over the 10-year period at 7 stations, with median reductions of about 2.5 percent. Increasing trends were observed at the two stations with the highest median TN concentration (Captain Hickory Run and SFLIL, 3–5 mg/L), both watersheds contain a high density of septic infrastructure. The combined trend results from all stations revealed no trend in TN and a declining trend in nitrate of about 2 percent per year. Overall, benthic community metrics indicated that streams throughout Fairfax County were initially of poor health; however, many metrics show an improving trend (from poor to fair based on the Fairfax County Index of Biological Integrity [IBI]). Significant increasing trends in IBI occurred at the network-scale and at 4 individual stations; additionally, scores improved by at least 1 qualitative category (for example, poor to fair, fair to good) at 11 of the 14 stations between 2009 (the first year all 14 stations were sampled) and 2017. Changes in all metrics suggest that the biodiversity, function, and condition of streams in Fairfax County are improving, but some of these improvements are driven by increased diversity and percent composition of organisms that are tolerant of the urban environment.

Virginia↗

USGS2016

Simulated effects of groundwater withdrawals from aquifers in Ocean County and vicinity, New Jersey

Rapid population growth since the 1930s in Ocean County and vicinity, New Jersey, has placed increasing demands upon the area’s freshwater resources. To examine effects of groundwater withdrawals, a three-dimensional groundwater-flow model was developed to simulate the groundwater-flow systems of five area aquifers: the unconfined Kirkwood-Cohansey aquifer system and Vincentown aquifer, and three confined aquifers— the Rio Grande water-bearing zone, the Atlantic City 800-foot sand, and the Piney Point aquifer. The influence of withdrawals is evaluated by using transient groundwater-flow model simulations that incorporate three withdrawal schemes. These are (1) no-withdrawal conditions; (2) 2000–03 withdrawal conditions, using reported monthly withdrawals at all production wells from January 2000 through December 2003; and (3) maximum-allocation withdrawal conditions using the maximum withdrawal allowed by New Jersey Department of Environmental Protection permits at each well. Particle tracking analysis, using results from model simulations, delineated particle flow paths from production wells to the point of recharge, and estimated particle travel times. Compared with no-withdrawal conditions, 2000–03 withdrawal conditions reduced the amount of groundwater flow out of the Kirkwood-Cohansey aquifer system into streams, increased the net flow of water into other layers, reduced net flow into or out of storage, and reduced flow from the Kirkwood-Cohansey aquifer system to constant head cells. Freshwater discharging to the Barnegat Bay-Little Egg Harbor estuary from streams and groundwater is essential to maintaining the ecology of the bay. Examination of selected stress periods indicates that simulated base flow in streams flowing into the Barnegat Bay-Little Egg Harbor estuary is reduced by as much as 49 cubic feet per second for 2000 to 2003 withdrawal conditions when compared with no-withdrawal conditions. In the three confined aquifers, water levels during periods of low recharge and high withdrawals, and high recharge and low withdrawals, were examined to determine seasonal effects on the confined flow systems. The simulated potentiometric surface of the Rio Grande water-bearing zone and the Atlantic City 800-foot sand during selected stress periods indicates substantial declines from no-withdrawal conditions to 2000–03 conditions as a result of groundwater withdrawals. Cones of depression in Toms River Township, Seaside Heights and Seaside Park Boroughs, and Barnegat Light Borough developed in the potentiometric surface of the Piney Point aquifer in response to withdrawals. Maximum-allocation withdrawals decreased flow out of the Kirkwood-Cohansey aquifer system to constant head cells, increased flow out of the aquifer system to adjacent and lower layers, and reduced groundwater discharge to streams when compared with 2000–03 withdrawal conditions. Increases in withdrawals from the Rio Grande water-bearing zone, the Atlantic City 800-foot sand, and the Piney Point aquifer result in an increase in simulated net groundwater flow into these aquifers. Base-flow reduction from 2000–03 conditions to maximum-allocation conditions of 25 to 29 cubic feet per second in all streams draining to the Barnegat Bay-Little Egg Harbor also is indicated. Potentiometric surfaces of the Rio Grande water-bearing zone, Atlantic City 800-foot sand, and the Piney Point aquifer during two stress periods of simulated maximum-allocation withdrawal conditions indicated the expansion of several cones of depression developed during 2000–03 withdrawals. Simulation of average 2000–03 withdrawal conditions indicated the extent to which the groundwater-flow system is susceptible to potential saltwater intrusion into near-shore wells. Travel time from recharge to discharge location ranged from 11 to approximately 50,700 years in near-shore Kirkwood-Cohansey aquifer system wells. Those in Seaside Heights Borough, in Island Beach State Park (Berkeley Township), and in Ship Bottom Borough have particle travel times from 140 to 12,000 years and flow paths that originated under Barnegat Bay or the Atlantic Ocean from the simulation of average maximum-allocation withdrawal conditions. Travel time along flow paths to wells screened in the Rio Grande water-bearing zone and the Atlantic City 800-foot sand from recharge to discharge point ranged from nearly 530 years to greater than 3.73 million years from the simulation of average 2000–03 withdrawal conditions. Particle tracking indicated that most wells screened in these aquifers derived a large part of their recharge from the Oswego River Basin, with a small portion of flow originating either beneath Barnegat Bay or to the east beneath the Atlantic Ocean. Travel time along flow paths that start beneath either Barnegat Bay or the Atlantic Ocean ranged from 2,300 to approximately 134,000 years from the simulation of average maximum-allocation withdrawal conditions."

New Jersey↗