Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Inland Waters”

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

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

628 records · Page 35Linked to original sources

Use of a Numerical Model to Simulate the Hydrologic System and Transport of Contaminants Near Joint Base Cape Cod, Western Cape Cod, Massachusetts

Historical training and operational activities at Joint Base Cape Cod (JBCC) on western Cape Cod, Massachusetts, have resulted in the release of contaminants into an underlying glacial aquifer that is the sole source of water to the surrounding communities. Remedial systems have been installed to contain and remove contamination from the aquifer. Groundwater withdrawals for public supply are expected to increase as the region continues to urbanize. Increases in water-supply withdrawals and wastewater return flow likely will affect the hydrologic system around JBCC and could affect the transport of any contamination that may remain in the aquifer following remediation of contamination from the JBCC. The U.S. Geological Survey, in cooperation with the Air Force Civil Engineer Center, developed a numerical, steady-state regional model of the Sagamore flow lens on western Cape Cod and evaluated the potential effects of future (2030) groundwater withdrawals on water levels, streamflows, hydraulic gradients, and advective transport near the JBCC. The aquifer consists generally of sandy sediments underlain by impermeable bedrock and is bounded laterally by a freshwater/saltwater interface. Data on the altitude of the bedrock surface, position of the freshwater/saltwater interface, lithology of the aquifer, spatial distribution of recharge, and hydrologic boundaries were incorporated into the three-dimensional, finite-difference groundwater flow model. Some inputs into the numerical model—aquifer properties, leakances, and recharge—are represented as parameters to facilitate estimation of optimal parameter values in an inverse calibration. A hybrid parameterization scheme, with both zones of piecewise constancy and pilot points, is used to represent hydraulic conductivity; other adjustable parameters include recharge, boundary leakance, and porosity. Data on water levels, the distribution of subsurface contamination, and groundwater ages were compiled, evaluated, and used to develop observations of long-term average hydraulic gradients and advective-transport patterns. These observations of steady-state hydrologic conditions were combined with the parameterized groundwater model in an inverse calibration to estimate model parameters that best fit the observations. Current (2010) and future (2030) conditions were simulated in the calibrated model to characterize the groundwater flow system and to determine potential effects of increased groundwater withdrawals on advective-transport patterns at the JBCC. Groundwater flow and advective transport are radially outward from a water-table divide in the northern part of the JBCC; flow diverges from the divide toward all points of the compass. Most groundwater flow and contaminant transport occur in shallow parts of the aquifer. On average, about one-half of the groundwater flux occurs in the shallowest 20 percent of the saturated thickness; shallow flow is even more predominant near streams and lakes. Projected (2030) increases in groundwater withdrawals decrease water levels by a maximum of about 1.2 feet in the northern part of the JBCC; drawdowns exceeding 1 foot generally are limited to areas near the largest increases in withdrawals, such as in the northern part of the JBCC, near Long Pond in Falmouth, and in eastern Barnstable. Streamflow decreases average about 6 percent; the largest decreases are in areas with the largest drawdowns. Changes in hydraulic-gradient directions at the water table exceed 1 degree in about 13 percent of the aquifer, generally near groundwater divides where gradient magnitudes are small and near large groundwater withdrawals. Predictions of advective transport from randomly selected locations at the water table are similar for current (2010) and future (2030) groundwater withdrawals. The results indicate that projected increases in groundwater withdrawals affect water levels and streamflows, but effects on hydraulic gradients and advective transport at the JBCC likely are small. Several underlying assumptions inherent in the model, including observations and weights used in the calibration, representation of local-scale heterogeneity, and simulation of the freshwater/saltwater interface, could affect model calibration and predictions; these assumptions were evaluated with alternative models and alternative inverse calibrations. Eight alternative calibrations were performed in which different, but reasonable, observations and weights were used. The preferred calibrated model had the best overall fit to the observations. Fine-grained silty sediments occur in many parts of the aquifer, and silt lenses can locally affect hydraulic gradients. A set of alternative models in which silts were represented with different correlation distances and hydraulic conductivities indicated that explicitly representing silt lenses could affect model calibration but that the implicit representation of local-scale heterogeneity may be sufficient at the regional scale to represent regional-scale hydraulic gradients. For the coastal boundary, two alternative models representing silty and sandy seabeds and their associated interface positions were developed to test the importance of the assumed coastal-boundary condition. The two alternative models resulted in different predictions of streamflow—streamflows increase with smaller (silty) seabed leakances. However, predictions of advective transport, particularly near the JBCC, generally were similar between the alternative and preferred calibrated models, indicating that the seabed leakance and associated interface position at the coastal boundary does not affect simulations of advective transport in inland parts of the aquifer.

Massachusetts↗

Characterization of peak streamflows and flood inundation of selected areas in southeastern Texas and southwestern Louisiana from the August and September 2017 flood resulting from Hurricane Harvey

Hurricane Harvey made landfall near Rockport, Texas, on August 25, 2017, as a Category 4 hurricane with wind gusts exceeding 150 miles per hour. As Harvey moved inland, the forward motion of the storm slowed down and produced tremendous rainfall amounts over southeastern Texas, with 8-day rainfall amounts exceeding 60 inches in some locations, which is about 15 inches more than average annual amounts of rainfall for eastern Texas and the Texas coast. Historic flooding occurred in Texas as a result of the widespread, heavy rainfall; wind and flood damages were estimated to be $125 billion, and the storm resulted in at least 68 direct fatalities. In the immediate aftermath of the Harvey-related flood event, the U.S. Geological Survey (USGS) and the Federal Emergency Management Agency initiated a cooperative study to evaluate the magnitude of the flood, determine the probability of occurrence, and map the extent of the flood in Texas. Seventy-four USGS streamflow-gaging stations in Texas with at least 15 years of record and no large data gaps in the period of record had a 2017 annual peak streamflow related to Harvey ranking in the top five of all annual peaks for each given station. New peaks of record streamflow were recorded at 40 of the 74 USGS streamflow-gaging stations. The number of years of peak streamflow record for the 74 analyzed streamflow-gaging stations ranged from 18 to 105, with a mean number of 55 years. The annual exceedance probability estimates for the analyzed streamflow-gaging stations ranged from less than 0.2 to 14.0 percent. USGS field crews surveyed 2,123 high-water marks to obtain water-surface elevations, in feet above the North American Vertical Datum of 1988. In some locations, several water-surface elevations were averaged to obtain 1 water-surface elevation, resulting in 1,258 water-surface elevations. Some of these high-water marks were used, along with peak-stage data from USGS streamflow-gaging stations, to create 19 inundation maps to document the areal extent of the maximum depth of the flooding. Digital datasets of the inundation area, modeling boundary, water-depth rasters, and final map products are available from the USGS data release associated with this report ( https://doi.org/10.5066/F7VH5N3N ).

Arkansas, Louisiana, Texas↗

Environmental quality and preservation; reefs, corals, and carbonate sands; guides to reef-ecosystem health and environment

Introduction In recent years, the health of the entire coral reef ecosystem that lines the outer shelf off the Florida Keys has declined markedly. In particular, loss of those coral species that are the building blocks of solid reef framework has significant negative implications for economic vitality of the region. What are the reasons for this decline? Is it due to natural change, or are human activities (recreational diving, ship groundings, farmland runoff, nutrient influx, air-borne contaminants, groundwater pollutants) a contributing factor and if so, to what extent? At risk of loss are biologic resources of the reefs, including habitats for endangered species in shoreline mangroves, productive marine and wetland nurseries, and economic fisheries. A healthy reef ecosystem builds a protective offshore barrier to catastrophic wave action and storm surges generated by tropical storms and hurricanes. In turn, a healthy reef protects the homes, marinas, and infrastructure on the Florida Keys that have been designed to capture a lucrative tourism industry. A healthy reef ecosystem also protects inland agricultural and livestock areas of South Florida whose produce and meat feed much of the United States and other parts of the world. In cooperation with the National Oceanic and Atmospheric Administration's (NOAA) National Marine Sanctuary Program, the U.S. Geological Survey (USGS) continues longterm investigations of factors that may affect Florida's reefs. One of the first steps in distinguishing between natural change and the effects of human activities, however, is to determine how coral reefs have responded to past environmental change, before the advent of man. By so doing, accurate scientific information becomes available for Marine Sanctuary management to understand natural change and thus to assess and regulate potential human impact better. The USGS studies described here evaluate the distribution (location) and historic vitality (thickness) of Holocene reefs in South Florida, relative to type of underlying bedrock morphology, and their varied natural response to rising sea level. These studies also assess movement and accumulation of sands, relative to direction of prevailing energy, and origin of the component sand grains. Geophysical data collected with highresolution sound-wave instruments that provide pictures of the sediment and bedrock are used to interpret sediment thickness. Reef thickness is determined by collecting limestone rock cores by drilling. Drill cores through reefs are used to identify the coral species that built them and to determine how reefs reacted to rising sea level. These data are supplemented by using isotope-dating techniques to derive the carbon-14 (C14) age of the corals and mangrove peat in the cores. Mangrove peat forms in very shallow water and at the shoreline but is found today buried beneath offshore reefs.

Open-File Report↗

FluOil: A novel tool for modeling the transport of oil-particle aggregates in inland waterways

Spilled oil in inland waterways can aggregate with mineral and organic particles to form oil-particle aggregates (OPAs). OPAs can be transported in suspension or deposited to the bed. Modeling the fate and transport of OPAs can provide useful information for making mitigation decisions. A novel open-source tool, FluOil, is developed to predict where OPAs may deposit and when they arrive in affected river/lake reaches by implementing the random walk particle tracking algorithm to represent the advection, diffusion, deposition, and resuspension of OPAs. The usability of FluOil is demonstrated with the 2010 Kalamazoo River oil spill case study. An unsteady hydrodynamic model simulates the river hydraulics and provides hydraulic data for use in FluOil. Settling velocity and critical shear stress for resuspension are the most important OPA properties concerning the transport and deposition of OPAs. Settling velocity determines the vertical distribution of OPAs and, thus, the travel speed, whereas critical shear stress determines where and when OPAs are deposited and resuspended.

Michigan↗

Using heat as a tracer to determine groundwater seepage in the Indian River Lagoon, Florida, April–November, 2017

The U.S. Geological Survey, in cooperation with the St. Johns River Water Management District, conducted a study to examine water fluxes in two small study areas in the Indian River Lagoon. Vertical arrays of temperature sensors were placed at multiple locations in the lagoon bed to measure temperature time series in the vertical profile. These data at one of the study areas, Eau Gallie, were used in two numerical models, 1DTempPro and VFLUX, to estimate seepage flux rates into the lagoon. 1DTempPro uses an inverse-modeling approach to calibrate groundwater flux to the measured temperature time series. VFLUX isolates the fundamental frequency signal in the temperature data and utilizes the resulting amplitude and phase differences between sensor locations to determine vertical water flux. Field measurements were made during two time periods, March 23 to April 28, 2017, and June 1 to November 3, 2017. Simulating the first, drier period at one location with 1DTempPro helped determine reasonable seepage fluctuations and provided guidelines for choosing which temperature sensor pairs used in the VFLUX simulations would produce the best results. VFLUX simulations at eight locations indicated daily average seepage flux rates of less than 20 centimeters per day (cm/d) and substantial seepage flux out to a distance of at least 110 meters from shore. The spatial variation in average seepage flux rates within 40 meters of shore seemed large, ranging from about 3 to 20 cm/d. In the VFLUX application using the June 1–November 3, 2017 data, the seepage flux has a higher magnitude and fluctuation than the first simulation period, making the isolation of the fundamental temperature frequency signal in the temperature data difficult. However, useful partial or full simulations were achieved at 6 of the 10 locations. The storm surge of Hurricane Irma on September 10, 2017, changed the depths of the sensors relative to the lagoon bed and disrupted the ability of VFLUX to compute seepage flux for the posthurricane period. The June 1 to November 3, 2017, computed seepage flux rates were higher than those for the March 24 to April 28, 2017, period and were sometimes as great as 40 cm/d, and more than 60 cm/d at one location. The seepage time-series data collected during Hurricane Irma indicates a downward seepage flux as a result of the storm surge, followed by upwelling from precipitation recharge inland. The average seepage flux rates are higher than those during the March–April period and are over 25 cm/d near the coast and about 20 cm/d 130 meters offshore.

Florida↗

Lithostratigraphic, geophysical, and hydrogeologic observations from a boring drilled to bedrock in glacial sediments near Nantucket Sound in East Falmouth, Massachusetts

In spring 2016, a 310-foot-deep boring (named MA–FSW 750) was drilled by the U.S. Geological Survey near Nantucket Sound in East Falmouth, Massachusetts, to investigate the hydrogeology of the southern coast of western Cape Cod. Few borings that are drilled to bedrock exist in the area, and the study area was selected to fill a gap between comprehensive geologic datasets inland to the north and marine geophysical data from beneath Nantucket Sound to the south. A permanent monitoring well (MA–FSW 750–0100) was installed in the boring upon the completion of the drilling and core collection. Observations from sediment cores and surface and borehole geophysical measurements were used to delineate three zones relevant to understanding groundwater flow at the study location. Shallow sands and gravels (0–107 feet [ft] below land surface [bls]) underlain by silt-rich fine and very fine sand (107–175 ft bls) form a zone of high permeability underlain by a zone of relatively lower permeability, referred to as the “shallow high-permeability” and “low-permeability” zones, respectively. A sharp lithological contact separating the shallow high-permeability and low-permeability zones may affect vertical flow of groundwater. Fine to coarse sand with intervals of clay and silt from 175 to 300 ft bls represent a deep zone of relatively high permeability, referred to as the “deep high-permeability” zone. A compacted, nonsorted unit (identified as basal till) and the bedrock surface were encountered at 300 and 305 ft bls, respectively. Hydraulic conductivity estimates from nuclear magnetic resonance logs and sediment grain-size distribution analyses indicated that the shallow high-permeability zone contributes substantially to the capacity of the aquifer to transmit groundwater at the study location. Results from geophysical surveys indicate a gradual transition from fresh to saline groundwater in the interval from 105 to 160 ft bls. Freshwater at the study site is present in the saturated unconsolidated sediments only in the 75 ft between 30 ft (the water table) and 105 ft bls in the shallow high-permeability zone. Sediments shallower than 175 ft bls closely resemble the downward fining post-Wisconsinan age deltaic and lacustrine deposits present in many parts of western Cape Cod; sediments deeper than 175 ft appear to be the product of earlier depositional processes more local to the southern coast of western Cape Cod. This study highlights how high-resolution observations of cored material coupled with a multitool geophysical approach can characterize a single boring to help better understand regional glacial history and hydrogeology.

Massachusetts↗

Inverse modeling for seawater intrusion in coastal aquifers: Insights about parameter sensitivities, variances, correlations and estimation procedures derived from the Henry problem

Inverse modeling studies employing data collected from the classic Henry seawater intrusion problem give insight into several important aspects of inverse modeling of seawater intrusion problems and effective measurement strategies for estimation of parameters for seawater intrusion. Despite the simplicity of the Henry problem, it embodies the behavior of a typical seawater intrusion situation in a single aquifer. Data collected from the numerical problem solution are employed without added noise in order to focus on the aspects of inverse modeling strategies dictated by the physics of variable-density flow and solute transport during seawater intrusion. Covariances of model parameters that can be estimated are strongly dependent on the physics. The insights gained from this type of analysis may be directly applied to field problems in the presence of data errors, using standard inverse modeling approaches to deal with uncertainty in data. Covariance analysis of the Henry problem indicates that in order to generally reduce variance of parameter estimates, the ideal places to measure pressure are as far away from the coast as possible, at any depth, and the ideal places to measure concentration are near the bottom of the aquifer between the center of the transition zone and its inland fringe. These observations are located in and near high-sensitivity regions of system parameters, which may be identified in a sensitivity analysis with respect to several parameters. However, both the form of error distribution in the observations and the observation weights impact the spatial sensitivity distributions, and different choices for error distributions or weights can result in significantly different regions of high sensitivity. Thus, in order to design effective sampling networks, the error form and weights must be carefully considered. For the Henry problem, permeability and freshwater inflow can be estimated with low estimation variance from only pressure or only concentration observations. Permeability, freshwater inflow, solute molecular diffusivity, and porosity can be estimated with roughly equivalent confidence using observations of only the logarithm of concentration. Furthermore, covariance analysis allows a logical reduction of the number of estimated parameters for ill-posed inverse seawater intrusion problems. Ill-posed problems may exhibit poor estimation convergence, have a non-unique solution, have multiple minima, or require excessive computational effort, and the condition often occurs when estimating too many or co-dependent parameters. For the Henry problem, such analysis allows selection of the two parameters that control system physics from among all possible system parameters.

Advances in Water Resources↗

Bright spots for inland fish and fisheries to guide future hydropower development

Hydropower production is one of the greatest threats to fluvial ecosystems and freshwater biodiversity. Now that we have entered the Anthropocene , there is an opportunity to reflect on what might constitute a ‘sustainable’ Anthropocene in the context of hydropower and riverine fish populations. Considering elements of existing practices that promote favorable social-ecological outcomes (i.e., ‘bright spots’) is timely given that there are plans to expand hydropower capacity in previously undammed rivers, intensify dam development in some of the world's largest river systems, and re-license existing facilities. We approach this from a pragmatic perspective: for the foreseeable future, hydropower will likely remain an important source of renewable electricity. To offer support for moving toward a more ‘sustainable’ Anthropocene, we provide syntheses of best practices during the siting, design, construction, operation, and compensation phases of hydropower development to minimize impacts on inland fish. For each phase, we offer positive examples (or what might be considered ‘bright spots’) pertaining to some of the approaches described within our syntheses, acknowledging that these projects may not be viewed as without ecological and (or) societal detriment by all stakeholders. Our findings underscore the importance of protecting critical habitat and free-flowing river reaches through careful site selection and basin-scale planning, infrastructure designs that minimize reservoir effects and facilitate safe passage of fish, construction of hydropower plants using best practices that minimize long-term damage, operating guidelines that mimic natural flow conditions, and compensation that is lasting, effective, inclusive, and locally relevant. Learning from these ‘bright spots’ may require engagement of diverse stakeholders, professionals, and governments at scales that extend well beyond a given site, river, or even basin. Indeed, environmental planning that integrates hydropower development into broader discussions of conserving regional biodiversity and ecosystem services will be of utmost importance.

Water Biology and Security↗

Differing modes of biotic connectivity within freshwater ecosystem mosaics

We describe a collection of aquatic and wetland habitats in an inland landscape, and their occurrence within a terrestrial matrix, as a “freshwater ecosystem mosaic” (FEM). Aquatic and wetland habitats in any FEM can vary widely, from permanently ponded lakes, to ephemerally ponded wetlands, to groundwater‐fed springs, to flowing rivers and streams. The terrestrial matrix can also vary, including in its influence on flows of energy, materials, and organisms among ecosystems. Biota occurring in a specific region are adapted to the unique opportunities and challenges presented by spatial and temporal patterns of habitat types inherent to each FEM. To persist in any given landscape, most species move to recolonize habitats and maintain mixtures of genetic materials. Species also connect habitats through time if they possess needed morphological, physiological, or behavioral traits to persist in a habitat through periods of unfavorable environmental conditions. By examining key spatial and temporal patterns underlying FEMs, and species‐specific adaptations to these patterns, a better understanding of the structural and functional connectivity of a landscape can be obtained. Fully including aquatic, wetland, and terrestrial habitats in FEMs facilitates adoption of the next generation of individual‐based models that integrate the principles of population, community, and ecosystem ecology.

Journal of the American Water Resources Associatio↗

Phytoplankton primary production in the world's estuarine-coastal ecosystems

Estuaries are biogeochemical hot spots because they receive large inputs of nutrients and organic carbon from land and oceans to support high rates of metabolism and primary production. We synthesize published rates of annual phytoplankton primary production (APPP) in marine ecosystems influenced by connectivity to land – estuaries, bays, lagoons, fjords and inland seas. Review of the scientific literature produced a compilation of 1148 values of APPP derived from monthly incubation assays to measure carbon assimilation or oxygen production. The median value of median APPP measurements in 131 ecosystems is 185 and the mean is 252 g C m −2 yr −1 , but the range is large: from −105 (net pelagic production in the Scheldt Estuary) to 1890 g C m −2 yr −1 (net phytoplankton production in Tamagawa Estuary). APPP varies up to 10-fold within ecosystems and 5-fold from year to year (but we only found eight APPP series longer than a decade so our knowledge of decadal-scale variability is limited). We use studies of individual places to build a conceptual model that integrates the mechanisms generating this large variability: nutrient supply, light limitation by turbidity, grazing by consumers, and physical processes (river inflow, ocean exchange, and inputs of heat, light and wind energy). We consider method as another source of variability because the compilation includes values derived from widely differing protocols. A simulation model shows that different methods reported in the literature can yield up to 3-fold variability depending on incubation protocols and methods for integrating measured rates over time and depth. Although attempts have been made to upscale measures of estuarine-coastal APPP, the empirical record is inadequate for yielding reliable global estimates. The record is deficient in three ways. First, it is highly biased by the large number of measurements made in northern Europe (particularly the Baltic region) and North America. Of the 1148 reported values of APPP, 958 come from sites between 30 and 60° N; we found only 36 for sites south of 20° N. Second, of the 131 ecosystems where APPP has been reported, 37% are based on measurements at only one location during 1 year. The accuracy of these values is unknown but probably low, given the large interannual and spatial variability within ecosystems. Finally, global assessments are confounded by measurements that are not intercomparable because they were made with different methods. Phytoplankton primary production along the continental margins is tightly linked to variability of water quality, biogeochemical processes including ocean–atmosphere CO 2 exchange, and production at higher trophic levels including species we harvest as food. The empirical record has deficiencies that preclude reliable global assessment of this key Earth system process. We face two grand challenges to resolve these deficiencies: (1) organize and fund an international effort to use a common method and measure APPP regularly across a network of coastal sites that are globally representative and sustained over time, and (2) integrate data into a unifying model to explain the wide range of variability across ecosystems and to project responses of APPP to regional manifestations of global change as it continues to unfold.

Biogeosciences↗

Gold deposits in the Xiaoqinling-Xiong'ershan region, Qinling mountains, central China

The gold-rich Xiaoqinling-Xiong'ershan region in eastern Shaanxi and western Henan provinces, central China, lies about 30-50 km inland of the southern margin of the North China craton. More than 100 gold deposits and occurrences are concentrated in the Xiaoqinling (west), Xiaoshan (middle), and Xiong'ershan (east) areas. Late Archean gneiss of the Taihua Group, and Middle Proterozoic metavolcanic rocks of the Xiong'er Group are the main host rocks for the deposits. Mesozoic granitoids (ca. 178-104 Ma) are present in most gold districts, but deposits are typically hosted in the Precambrian basement rocks hundreds of meters to as far as 10 km from the intrusions and related hornfels zones. Deposits in the Xiaoqinling and Xiaoshan areas are best classified as orogenic gold deposits, with ores occurring in a number of distinct belts both in quartz veins and disseminated in altered metamorphic rocks. Alteration assemblages are dominated by quartz, sericite, pyrite, and carbonate minerals. The ore-forming fluids were low salinity, CO2-rich, and characterized by isotopically heavy ??18O. Four deposits (Dongchuang, Wenyu, Yangzhaiyu, and Dahu) in the Xiaoqinling area each contain resources of about 1 Moz Au. Some of the gold deposits in the Xiong'ershan area represent more shallowly emplaced tellurium-enriched orogenic systems, which include resources of approximately 1-1.5 Moz Au at Shanggong and Beiling (or Tantou). Others are epithermal deposits (e.g., Qiyugou and Dianfang) that are hosted in volcanic breccia pipes. Isotopic dates for all gold deposits, although often contradictory, generally cluster between 172-99 Ma and are coeval with emplacement of the post-kinematic granitoids. The gold deposits formed during a period of relaxation of far-field compressional stresses, clearly subsequent to the extensive Paleozoic-early Mesozoic accretion of are terranes and the Yangtze craton onto the southern margin of the North China craton. Hydrothermal and magmatic events occurred locally where extension-related Precambrian basement uplifting took place along the craton margin. Fluids for the orogenic gold deposits in the Xiaoqinling, Xiaoshan, and Xiong'ershan areas may have been released from evolving magmas or resulted from prograde metamorphic reactions within the uplift zones. Alternatively, for the epithermal gold deposits at shallower levels in the Xiong'ershan area, gold-transporting fluids were mainly exsolved from coeval magmas, although meteoric water was also involved in these hydrothermal systems.

Mineralium Deposita↗

Numerical simulation of flow in deep open boreholes in a coastal freshwater lens, Pearl Harbor Aquifer, O‘ahu, Hawai‘i

The Pearl Harbor aquifer in southern O‘ahu is one of the most important sources of freshwater in Hawai‘i. A thick freshwater lens overlays brackish and saltwater in this coastal aquifer. Salinity profiles collected from uncased deep monitor wells (DMWs) commonly are used to monitor freshwater-lens thickness. However, vertical flow in DMWs can cause the measured salinity to differ from salinity in the adjacent aquifer or in an aquifer without a DWM. Substantial borehole flow and displacement of salinity in DMWs over several hundred feet have been observed in the Pearl Harbor aquifer. The objective of this study was to evaluate the effects of borehole flow on measured salinity profiles from DMWs. A numerical modeling approach incorporated aquifer hydraulic characteristics and recharge and withdrawal rates representative of the Pearl Harbor aquifer. Borehole flow caused by vertical hydraulic gradients associated with both the natural regional flow system and groundwater withdrawals was simulated. Model results indicate that, with all other factors being equal, greater withdrawal rates, closer withdrawal locations, or higher hydraulic conductivities of the well cause greater borehole flow and displacement of salinity in the well. Borehole flow caused by the natural groundwater-flow system is five orders of magnitude greater than vertical flow in a homogeneous aquifer, and borehole-flow directions are consistent with the regional flow system: downward flow in inland recharge areas and upward flow in coastal discharge areas. Displacement of salinity inside the DMWs associated with the regional groundwater-flow system ranges from less than 1 to 220 ft, depending on the location and assumed hydraulic conductivity of the well. For example, upward displacements of the 2 percent and 50 percent salinity depths in a well in the coastal discharge part of the flow system are 17 and 4.4 ft, respectively, and the average salinity difference between aquifer and borehole is 0.65 percent seawater salinity. Groundwater withdrawals and drawdowns generally occur at shallow depths in the freshwater system with respect to the depth of the DMW and cause upward flow in the DMW. Simulated groundwater withdrawal of 4.3 million gallons per day that is 100 ft from a DMW causes thirty times more borehole flow than borehole flow that is induced by the regional flow field alone. The displacement of the 2 percent borehole salinity depth increases from 17 to 33 ft, and the average salinity difference between aquifer and borehole is 0.85 percent seawater salinity. Peak borehole flow caused by local groundwater withdrawal near DMWs is directly proportional to the pumping rate in the nearby production well. Increasing groundwater withdrawal to 16.7 million gallons per day increases upward displacement of the 50 percent salinity depth (midpoint of the transition zone) from 4.6 to 77 ft, and the average salinity difference between aquifer and borehole is 1.4 percent seawater salinity. Simulated groundwater withdrawal that is 3,000 ft away from DMWs causes less borehole flow and salinity displacements than nearby withdrawal. Simulated effects of groundwater withdrawal from a horizontal shaft and withdrawal from a vertical well in a homogeneous aquifer were similar. Generally, the 50 percent salinity depths are less affected by borehole flow than the 2 percent salinity depths. Hence, measured salinity profiles are useful for calibration of regional numerical models despite borehole-flow effects. Commonly, a 1 percent error in salinity is acceptable in numerical modeling studies. Incorporation of heterogeneity in the model is necessary to simulate long vertical steps observed in salinity profiles in southern O‘ahu. A thick zone of low aquifer hydraulic conductivity limits exchange of water between aquifer and well and creates a long vertical step in the salinity profile. A heterogeneous basalt-aquifer scenario simulates observed vertical salinity steps and borehole flow that is consistent with measured borehole flow from DMWs in southern O‘ahu. However, inclusion of local-scale heterogeneities in regional models generally is not warranted.

Hawai'i↗

Vanadium

Vanadium is used primarily in the production of steel alloys; as a catalyst for the chemical industry; in the making of ceramics, glasses, and pigments; and in vanadium redox-flow batteries (VRBs) for large-scale storage of electricity. World vanadium resources in 2012 were estimated to be 63 million metric tons, which include about 14 million metric tons of reserves. The majority of the vanadium produced in 2012 was from China, Russia, and South Africa. Vanadium is extracted from several different types of mineral deposits and from fossil fuels. These deposits include vanadiferous titanomagnetite (VTM) deposits, sandstone-hosted vanadium (with or without uranium) deposits (SSV deposits), and vanadium-rich black shales. VTM deposits are the principal source of vanadium and consist of magmatic accumulations of ilmenite and magnetite containing 0.2 to 1 weight percent vanadium pentoxide (V 2 O 5 ). SSV deposits are another important source; these deposits have average ore grades that range from 0.1 to greater than 1 weight percent V 2 O 5 . The United States has been and is currently the main producer of vanadium from SSV deposits, particularly those on the Colorado Plateau. Vanadium-rich black shales occur in marine successions that were deposited in epeiric (inland) seas and on continental margins. Concentrations in these shales regularly exceed 0.18 weight percent V 2 O 5 and can be as high as 1.7 weight percent V 2 O 5 . Small amounts of vanadium have been produced from the Alum Shale in Sweden and from ferrophosphorus slag generated during the reduction of phosphate to elemental phosphorus in ore from shales of the Phosphoria Formation in Idaho and Wyoming. Because vanadium enrichment occurs in beds that are typically only a few meters thick, most of the vanadiferous black shales are not currently economic, although they may become an important resource in the future. Significant amounts of vanadium are recovered as byproducts of petroleum refining, and processing of coal, tar sands, and oil shales may be important future sources. Vanadium occurs in one of four oxidation states in nature: +2, +3, +4, and +5. The V 3+ ion has an octahedral radius that is almost identical to that of (Fe 3+ ) and (Al 3+ ) and, therefore, it substitutes in ferromagnesian minerals. During weathering, much of the vanadium may partition into newly formed clay minerals, and it either remains in the +3 valence state or oxidizes to the +4 valence state, both of which are relatively insoluble. If erosion is insignificant but chemical leaching is intense, the residual material may be enriched in vanadium, as are some bauxites and laterites. During the weathering of igneous, residual, or sedimentary rocks, some vanadium oxidizes to the +5 valence state, especially in the intensive oxidizing conditions that are characteristic of arid climates. The average contents of vanadium in the environment are as follows: soils [10 to 500 parts per million (ppm)]; streams and rivers [0.2 to 2.9 parts per billion (ppb)]; and coastal seawater (0.3 to 2.8 ppb). Concentrations of vanadium in soils (548 to 7,160 ppm) collected near vanadium mines in China, the Czech Republic, and South Africa are many times greater than natural concentrations in soils. Additionally, if deposits contain sulfide minerals such as chalcocite, pyrite, and sphalerite, high levels of acidity may be present if sulfide dissolution is not balanced by the presence of acid-neutralizing carbonate minerals. Some of the vanadium-bearing deposit types, particularly some SSV and black-shale deposits, contain appreciable amounts of carbonate minerals, which lowers the acid-generation potential. Vanadium is a micronutrient with a postulated requirement for humans of less than 10 micrograms per day, which can be met through dietary intake. Primary and secondary drinking water regulations for vanadium are not currently in place in the United States. Vanadium toxicity is thought to result from an intake of more than 10 to 20 milligrams per day. Vanadium is essential for some biological processes and organisms. For example, some nitrogen-fixing bacteria require vanadium for producing enzymes necessary to convert nitrogen from the atmosphere into ammonia, which is a more biologically accessible form of nitrogen.

Professional Paper↗

Glaciers along proposed routes extending the Copper River Highway, Alaska

Three inland highway routes are being considered by the Alaska Department of Transportation and Public Facilities to connect the community of Cordova in southcentral Alaska to a statewide road system. The routes use part of a Copper River and Northwest Railway alignment along the Copper River through mountainous terrain having numerous glaciers. An advance of any of several glaciers could block and destroy the roadway, whereas retreating glaciers expose large quantities of unconsolidated, unvegetated, and commonly ice-rich sediments. The purpose of this study was to map historical locations of glacier termini near these routes and to describe hazards associated with glaciers and seasonal snow. Historical and recent locations of glacier termini along the proposed Copper River Highway routes were determined by reviewing reports and maps and by interpreting aerial photographs. The termini of Childs, Grinnell, Tasnuna, and Woodworth Glaciers were 1 mile or less from a proposed route in the most recently available aerial photography (1978-91); the termini of Allen, Heney, and Schwan Glaciers were 1.5 miles or less from a proposed route. In general, since 1911, most glaciers have slowly retreated, but many glaciers have had occasional advances. Deserted Glacier and one of its tributary glaciers have surge-type medial moraines, indicating potential rapid advances. The terminus of Deserted Glacier was about 2.1 miles from a proposed route in 1978, but showed no evidence of surging. Snow and rock avalanches and snowdrifts are common along the proposed routes and will periodically obstruct the roadway. Floods from ice-dammed lakes also pose a threat. For example, Van Cleve Lake, adjacent to Miles Glacier, is as large as 4.4 square miles and empties about every 6 years. Floods from drainages of Van Cleve Lake have caused the Copper River to rise on the order of 20 feet at Million Dollar Bridge.

Water-Resources Investigations Report↗

Geology and paleontology of the Santa Maria district, California

Stratigraphy, paleontology, and geologic history.-A basement' consisting of igneous rocks of the Jurassic(?) Franciscan formation and sediments of the Upper Jurassic Knoxville formation, and formations of Tertiary and Quaternary age are exposed in the Santa Maria district. The outcrop section, exclusive of the Franciscan, has a maximum thickness of about 10,000 feet, the subsurface section about 27,000 feet. At no locality, however, is either outcrop or subsurface section as thick as the total maxima for the formations. The outcrop Franciscan is made up of altered basalt, gabbro (generally- greatly altered), and minor areas of peridotite and serpentine. The marine Knoxville formation, consisting of shale, thin-bedded calcareous sandstone, and conglomerate, was deposited on the igneous rocks of the Franciscan. It has an outcrop thickness of about 500 feet and a known subsurface thickness of at least 1,250 feet. Aucella cf. A. piochii, which occurs in both outcrop and subsurface sections, suggests late Jurassic age. The Franciscan and Knoxville were probably uplifted during Taliaferro's Diablan orogeny at the close of the Jurassic and formed a land area that presumably endured during early Cretaceous and perhaps during much of late Cretaceous time. If sediments were laid down while 10,000 feet of Upper Cretaceous were deposited in the adjoining San Rafael Mountains, they were eroded following uplift at the close of the Cretaceous, when Reed and Hollister's San Rafael uplift was formed. The district is inferred to have been part of a Franciscan and Knoxville land area at the south border of the San Rafael uplift during all of early Tertiary time. The known history of the district as part of a Tertiary basin began in the early Miocene(?), possibly a little earlier or possibly 3: little later, when the nonmarine sediments of the early Miocene(?) Lospe formation were deposited. The Lospe formation has a maximum thickness of 2, 700 feet, and is made up of coarse-grained reddish sandstone and conglomerate, and greenish sandstone, gypsiferous siltstone, and mudstone. White tuff is a minor but conspicuous constituent. The Lospe formation overlies the Knoxville, or· overlaps it and rests on the Franciscan. The first known Tertiary invasion of the sea took place immediately thereafter in early middle Miocene time, and from then on until approximately the end of the Pliocene the sea occupied continuously at least most of the region. The extent of the Miocene basin is not certainly known, but it was part of an extensive· basin that has been designated the Santa Barbara embayment. The early middle Miocene Point Sal formation is the earliest marine Tertiary formation. It has an outcrop thickness of as much as 1,500 feet and a maximum subsurface thickness of 3,600 feet, and consists of siltstone, mudstone, and thin beds of sandstone. The Point Sal formation overlies the Lospe formation without marked discontinuity, or overlaps it and rests on Knoxville or Franciscan. It contains a large foraminiferal fauna representing the Siphogenerina hughesi zone, or the lower part of Kleinpell's Relizian stage. The Monterey shale overlies the Point Sal formation without noticeable discontinuity, or overlaps all the older sedimentary formations and rests on the Franciscan. It has a maximum outcrop thickness of 2,100 feet, but is as much as 5,000 feet thick in some subsurface sections. The Monterey is divided into three mapped members. The lower member is characterized by phosphatic shale and somewhat porcelaneous shale; the middle member by chert and cherty shale; and the upper member by porcelaneous shale, or by both porcelaneous shale and diatomaceous strata. The lower member contains Foraminifera representing the upper part of Kleinpell's Relizian stage and all of his Luisian stage, the middle member a few species indicating the lower part of his Mohnian stage, and the upper member a fauna representing the upper part of the Bolivina hughesi zone, at the top of the Mohnian, and in part a younger unnamed faunal division; that is, the Monterey is of late middle and late Miocene age. The diverse stratigraphic relations of the Lospe and Point Sal formations and Monterey shale along the borders of the Franciscan rocks forming Point Sal Ridge are inferred to be the result of repeated movements during Miocene time in the area of basement rocks west of the district, presumably extending westward beyond the present coast. Toward the close of the Miocene, low ridges appeared on the floor of the sea, ridges that grew during Pliocene time and were destined to become anticlines. At the same time deformation took place in the northeastern part of the district, and elsewhere on some anticlines, or on other structural highs bounded by faults. In those areas the Sisquoc formation overlies the Monterey shale with marked discordance; elsewhere there is no discordance. Two facies of the Sisquoc formation are mapped: a marginal sandstone facies, designated the Tinaquaic sandstone member, and a basin facies. The Tinaquaic sandstone member is 1,400 feet thick and contains megafossils of middle Pliocene age (also early Pliocene just east of the mapped area). It is unconformable on the Monterey. The basin facies, at least 3,000 feet thick in outcrop sections and 5,000 in some subsurface sections, consists of diatomaceous mudstone, other types of diatomaceous ·strata, somewhat porcelaneous mudstone, and porcelaneous shale-deposits that are ordinarily characteristic of the Monterey shale. Even in areas where the two formations are lithologically indistinguishable and conformable, a field basis for differentiating them has been established. The basin facies of the Sisquoc conformably overlies the Monter.ey in outcrop sections, but in some subsurface sections the formations are unconformable, and on the north limb of the Santa Maria Valley syncline the basin facies of the Sisquoc overlaps the Monterey onto the basement, thus forming the overlap trap for the oil in the Monterey in the Santa Maria Valley field. The lower and middle parts of the basin facies contain Foraminifera of the Bolivina obliqua zone. Kleinpell assigned that zone to the lower part of his late upper Miocene Delmontian stage. Assignment to the upper Delmontian, however, appears to be preferable. The upper few hundred feet of the basin facies contain Foraminifera similar to those in the overlying Foxen mudstone and megafossils of middle Pliocene affinities. The basin facies is therefore considered late upper Miocene to middle Pliocene. The submarine ridges were growing during Pilocene time. The Foxen mudstone is missing on them and in the northeastern part of the district. On the north limb of the Santa Maria Valley syncline, the Foxen overlaps the Sisquoc formation and rests on the basement. In the basins between the submarine ridges, the mudstone, siltstone, and fine-grained sandstone of the Foxen (800 feet thick in outcrop sections and as much as 2, 750 in subsurface sections) overlie conformably the Sisquoc formation. In some of the areas where the Foxen is missing as a lithologic unit, it appears to be represented by condensed deposits of phosphatic pellets, mapped with the underlying or overlying Formation, depending on the matrix, or by a condensed section of fine-grained sand, mapped with lithologically indistinguishable sand in the basal part of the overlying formation. Though foraminifera are abundant in the Foxen mudstone, relatively few species are represented. Those from the lower part of the formation may be of middle Pliocene age. Megafossils from the upper part of the Foxen are considered late Pliocene. The formation is therefore assigned to the middle(?) and upper Pliocene. Movements during Pliocene time in the basement area west of Point Sal Ridge are indicated by the occurrence in the western Casmalia Hills of coarse detritus from the Monterey in the upper part of the Sisquoc formation and in the Foxen mudstone. During late Pliocene time, when the Careaga sandstone was deposited, the Pliocene sea had its greatest extent. The region then flooded may be referred to as the Santa Maria basin. Throughout most of the district two mapped members of the Careaga sandstone are differentiated: the Cebada fine-grained member, which has a maximum outcrop thickness of 1,000 feet, and the Graciosa coarse-grained member, 50 to 425 feet thick. Mild deformation continued during the late Pliocene. The Cebada fine-grained member is missing on the embryonic anticlines, where it is overlapped by the Graciosa coarse-grained member. Both members of the Careaga sandstone contain a large megafauna. The nonmarine Paso Robles formation conformably overlies the Careaga sandstone. The Paso Robles consists chiefly of sand and gravel, but clay, marl, and limestone are the most characteristic constituents. The maximum outcrop thickness is 2,000 feet, the estimated maximum subsurface thickness 4,500 feet. The Paso Robles fresh-water fauna is meager, consisting of a few species similar to living forms. The formation is currently assigned to the interval including late Pliocene and early Pleistocene(?). The age assignment late Pliocene(?) and early Pleistocene, however, may be preferable for the Santa Maria district. Then followed the only period of strong general deformation in the known Tertiary and Pleistocene history of the district. The present structural features of the district were formed at that time, and the submarine ridges appeared as fully formed anticlines. Dating of the deformation is uncertain, because of uncertainty concerning the age of the Paso Robles formation. It is, however, without much doubt of the same age as the welldated strong middle Pleistocene deformation in the Ventura basin. Terrace deposits, laid down on both wave-cut and stream-cut platforms, are rather arbitrarily assigned to the late Pleistocene. The oldest and most extensive terrace deposits, not more than 100 feet thick, are designated the Orcutt sand. The Orcutt sand itself is tilted as much as 12° on the limbs of anticlines and is faulted on the north limb of the Graciosa anticline, indicating renewed growth of the anticlines, presumably in late Pleistocene time. Terrace deposits apparently younger than the Orcutt sand are arched in a low anticline west of lower Foxen Canyon. Structure.-Santa Maria Valley is the boundary between two structural provinces. To the north, valleys and hills are either synclines or anticlines. To the south, on the contrary, major valleys coincide generally with major synclines, and the hills are anticlinal. Santa Maria Valley itself is a syncline. Unlike valleys farther south in the district, however, it lies athwart an older uplift. Also unlike most of the valleys farther south, the axis of the syncline is not in the middle of the valley, but is far to the south near the bordering hills. Indeed, in the western Casmalia Hills, the anticline in the bordering hills is overturned and overrides the axis of the syncline. The major structural features of the district have a general west-northwestward trend parallel to the trend of the basin. Minor westward-trending and northward-trending folds and faults, however, extend across the trend of the major features. The district includes areas of wide, open folds and also areas of narrow, closely spaced, and steeply tilted folds, as well as some major overturned anticlines, most of the latter overturned northward. The closely spaced folds coincide almost invariably with outcrops of the Monterey shale and Point Sal formation. Physiography.-The surface on which the terrace deposits, designated the Orcutt sand, were deposited is extensive but is locally deformed, and only remnants are preserved. Toward the coast it changes from a stream-cut surface to a wave-cut surface. In the coastal area three main marine terraces are recognized: the high terrace (altitude about 800 feet), the intermediate terrace (altitude about 600 feet), and the low terrace (altitude 50 to 125 feet). An indurated layer is present at or near the surface at many localities scattered throughout the district. It is suggested that the indurated layer is an ancient hardpan, the incomplete skeleton of a former soil profile developed on a former surface of less relief than the present surface. Sand dunes extend inland from the coast at the north and south borders of the mapped area. They are classified under three age groups: old, intermediate, and modern. The old dunes, which have a protective cover of natural vegetation and are now inactive, are far more extensive than those of the other two groups. They cover many square miles on a terrace bordering Santa Maria Valley, and extend 20 miles inland. It has not been determined whether their inactivity is due to a cutting off of the supply of sand, or to a climatic change. Occurrence of oil.--Oil has been produced in the Santa Maria district since 1901, the total production to the end of 1947 being 269,657,000 barrels. Throughout the district most of the oil is heavy. Seven producing fields are located in the mapped area. In order of discovery from oldest to youngest the fields are as follows: Orcutt, Lompoc, West Cat Canyon, East Cat Canyon, Casmalia, Gato Ridge, and Santa Maria Valley. The Santa Maria Valley field is the largest in both area and productive capacity. It is the largest overlap field in coastal California, and one of the last major fields found in the State up to 1947. The Las Flores (Monterey) pool is transforming the West Cat Canyon field into one of the major fields of the district. The Monterey shale is the chief oil-bearing formation, and the principal reservoir in the Monterey consists of fractured chert and cherty shale. Sand in the Sisquoc formation is the sole reservoir in the minor East Cat Canyon field and in the Pliocene pool of the West Cat Canyon field. The Point Sal formation yields some oil in the southeastern part of the Santa Maria Valley field (the only part of that field where the formation is present) and a recent well in the Casmalia field is producing a small amount of relatively light oil from the Point Sal. The Lospe formation is also productive in a recent well in the Casmalia field. The Knoxville formation is productive in three areas in the northern part of the Santa Maria Valley field. The Point Sal formation offers the greatest promise for deeper-zone production. The Lospe and Knoxville formations can no longer be ignored in areas where younger marine formations overlap against them. Oil possibilities in undeveloped areas.--Two matters weigh heavily in prospecting in the Santa Maria district: the degree of fracturing of chert and cherty shale in the Monterey shale, and the gravity of the oil. Other things being equal, the productivity of the Monterey varies directly with the amount of fracturing. Very heavy oil, too heavy to produce commercially under present conditions, has befln found ii the Monterey in the northeastern part of the district, where five discoveries (one east of the mapped area) have been made in recent years. Among areas of possible interest, three appear to be favorable for prospecting on the basis of surface geology: an area east of Foxen Canyon, where oil may be trapped in the basal part of the Tinaquaic sandstone member of the Sisquoc formation by westward overlap of successively higher Tinaquaic strata onto the Monterey shale; an area so1tth of the I ... ions Head fault, where oil may be trapped by the fault; and the offshore extension of the north border of Point Sal Ridge, where oil may possibly be trapped in the Monterey by overlap of the Sisquoc formation.

California↗

Floods of November-December 1950 in the Central Valley basin, California

The flood of November-December 1950 in the Central Valley basin was the greatest in most parts of the basin since the turn of the century and probably was exceeded in the lower San Joaquin River basin only by the historic flood of 1862. In respect to monetary loss, the 1950 flood was the most disastrous in the history of the basin. Loss of life was remarkably small when one considers the extensive damage and destruction to homes and other property, which is estimated at 33 million dollars. Outstanding features of the flood were its unprecedented occurrence so early in the winter flood season, its magnitude in respect to both peak and volume in most major tributaries, and the occurrence of a succession of near-peak flows with a period of three weeks. The flood was caused by a series of storms during the period November 16 to December 8, which brought exceptionally warm, moisture-laden air inland against the Sierra Nevada range and caused intense rainfall, instead of snowfall, at unusually high altitudes. Basin-wide totals of rainfall during the period ranged from 30 inches over the Yuba and American River basins to 13 inches over the upper Sacramento and Feather River basins. Based on continuous records of discharge on major tributaries for periods ranging from 22 to 55 years and averaging about 43 years, the 1950 flood peaks were the greatest of record on the American, Cosumnes, Mokelumne, Stanislaus, Tuolumne, Merced, Chowchilla, Fresno, lower San Joaquin, Kings, Kaweah, Tule, and Kern Rivers. Second highest peak of record occurred during the flood of March 1928 on the Yuba, American and Mokelumne Rivers; the flood of Marcn 1940 on Cosumnes River; the flood of January 1911 on the Stanislaus and Tuolumne Rivers; the flood of December 1937 on the Merced, Kings, and Kaweah Rivers; the flood of March 1938 on the Chowchilla, Fresno, and lower San Joaquin Rivers; and the flood of March 1943 on the Tule and Kern Rivers. Peak discharges for 1950 did not exceed previous maxima on Bear, Yuba, Feather, and upper Sacramento Rivers, nor on west side tributaries of lower Sacramento River, Calaveras River, and upper San Joaquin River (above Friant Reservoir). Notable high rates of discharge were 354 cfs per square mile from 39.5 square miles in North Fork of Middle Fork Tule River, 225 cfs per square mile from 198 square miles in Rubicon River, 115 cfs per square mile from 999 square miles in North Fork of American River and 93.7 cfs per square mile from 1,921 square miles in American River at Fair Oaks. This report presents a general description of the 1950 flood, details and estimates of the damage incurred, records of stage and discharge for the period of the flood at 171 stream-gaging stations, records of storage in 14 reservoirs, a summary of peak discharges with comparative data for previous floods at 252 measurement points, and tables showing crest stages along the main stem and major tributary channels of the Sacramento and San Joaquin Rivers. The report also includes a discussion of meteorologic and hydrologic conditions associated with the flood, examples of the flood regulation afforded by storage reservoirs, a brief study of runoff characteristics, and a summary and comparison with previous floods in the Central Valley basin.

California↗