SEARCH · Geology Reports

Results for “Studies in Maine Geology”

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 991 records · Page 55Linked to original sources

USGS1973

Species diversity: Patterns in modern and Miocene foraminifera of the eastern margin of North America

Patterns of foraminiferal species diversity were examined along the eastern margin of North America by utilizing the number of species, S, the information function, H(S), and species equitability, E. The 350 modern samples we studied extended from the Arctic to the Gulf of Mexico at depths ranging from a few meters to more than 5,000 m. In addition, 29 samples from Miocene strata of the Atlantic Coastal Plain and continental shelf were studied. Modern species diversity as measured by S and H(S) generally increases as depth increases and latitude decreases. Some notable exceptions occur, however, which are difficult to explain. For example, species diversity in the Arctic depth interval of 0 to 100 m is as high or higher than that found immediately south of Nova Scotia, in the Gulf of Maine, on Browns and Georges Banks, and even off the Gulf of Mexico deltas. At the moderate depth interval of 100 to 1,000 m, however, the entire margin north of Browns and Georges Banks has lower diversities than that to the south. The highest diversity by far in this depth interval occurs in the northeastern Gulf of Mexico. At the depth interval greater than 1,000 m, the more southern areas studied generally have a higher species diversity than the more northern Cape Cod to Maryland area. An exception to this is the northwestern Gulf of Mexico; this area is also an exception in that species diversity is significantly lower in the deeper waters than in the shallower waters in the same area. The measure of species equitability, E, showed no clear pattern with depth or latitude. This may be so because no simple pattern of species proportions exists or because the sampling was inadequate to measure it. Samples from the Miocene strata show a striking resemblance in species diversity to modern samples at similar depths and latitudes. Our observations indicate that species diversity and equitability have not increased during the last 15 × 10 6 yrs. The fossil and modern data indicate that each environment has its own carrying capacity and that this capacity is reached rather quickly. Although time and environmental stability are undoubtedly important in determining species diversity, as presently defined they are inadequate to explain all observed patterns. Long-term observations in various environments will be required to determine the relative importance of variables that affect species diversity.

Geological Society of America Bulletin↗

USGS1969

Micromineralogy of silver-bearing sphalerite from Flat River, Missouri

Detailed mineralogical and chemical study of sphalerite-rich lead ores from Flat River, Mo., confirms the presence of anomalous amounts of silver in the sphalerite. Although silver is closely associated with chlorine and no silver sulfide minerals were identified, geochemical considerations indicate the silver may be in the form of discrete submicron-size grains of sulfide rather than chloride. However, the close correlation between abundance variations in silver, chlorine, cadmium, iron, and minor lead suggests the possible existence of submicron-size grains of a complex chloride or oxychloride. Results from detailed chemical study of sphalerite show: (1) chlorine is present in areas containing silver and changes in abundance of both elements are similar; (2) the mole ratio of chlorine to silver varies from approximately 3:1 to 6:1 with an approximate average of 5:1 and the total content of each element varies widely between grains and within the same grain; (3) taking total silver in any analysis and calculating the amount of chlorine required for silver chloride leaves an excess of chlorine; (4) in chlorine-deficient areas the amounts of iron and cadmium substituting for zinc in the sphalerite structure vary widely; and (5) in areas containing chlorine, changes in abundance of iron and cadmium vary directly with that of chlorine. The fine-grained dusty opaque inclusions characterize sphalerite from this area; they were identified as mainly galena with lesser amounts of pyrite and cassiterite. Galena is silver-free (<0.02 weight percent Ag) and is essentially pure lead sulfide. Other sulfides identified include pyrite, a cobalt-nickel-iron sulfide of composition Co (sub 0.55) Ni (sub 0.25) Fe (sub 0.20) S (sub 2.00), and an iron-cobalt-nickel sulfide containing 22-28 weight percent Fe, 10-15 weight percent Co, and 3-5 weight percent Ni. Anomalous amounts of tin present as cassiterite are associated with the gangue minerals quartz, potash feldspar, a kaolinite-type clay, and trace amounts of cuprite. Carbonate gangue minerals include calcite, dolomite, ankerite, magnesite, and minor amounts of cerussite. Chemical analyses were made of many minerals and physical and textural relationships examined.

Missouri↗

USGS2017

Noble gas data from Goldfield and Tonopah epithermal Au-Ag deposits, ancestral Cascades Arc, USA: Evidence for a primitive mantle volatile source

The He, Ne, and Ar isotopic composition of fluid inclusions in ore and gangue minerals were analyzed to determine the source of volatiles in the high-grade Goldfield and Tonopah epithermal Au-Ag deposits in southwestern Nevada, USA. Ar and Ne are mainly atmospheric, whereas He has only a minor atmospheric component. Corrected 3 He/ 4 He ratios (with atmospheric He removed) range widely from 0.05 to 35.8 times the air 3 He/ 4 He ratio (R A ), with a median of 1.43 R A . Forty-one percent of measured 3 He/ 4 He ratios are ≥4 R A , corresponding to ≥50% mantle He assuming a mantle ratio of 8 R A . These results suggest that mafic magmas were part of the magmatic-hydrothermal system underlying Goldfield and Tonopah, and that associated mantle-sourced volatiles may have played a role in ore formation. The three highest corrected 3 He/ 4 He ratios of 17.0, 23.7, and 35.8 R A indicate a primitive mantle He source and are the highest yet reported for any epithermal-porphyry system and for the Cascades arc region. Compiled 3 He/ 4 He measurements from epithermal-porphyry systems in subduction-related magmatic arcs around the world (n = 209) display a statistically significant correlation between 3 He/ 4 He and Au-Ag grade. The correlation suggests that conditions which promote higher fluid inclusion 3 He/ 4 He ratios (abundance of mantle volatiles and focused upward volatile transport) have some relation to conditions that promote higher Au-Ag grades (focused flow of metal-bearing fluids and efficient chemical traps). Results of this and previous investigations of He isotopes in epithermal-porphyry systems are consistent with the hypothesis posed in recent studies that mafic magmas serve an important function in the formation of these deposits.

Cascades Arc↗

USGS2010

Biological water-quality assessment of selected streams in the Milwaukee Metropolitan Sewerage District Planning Area of Wisconsin, 2007

Changes in the water quality of stream ecosystems in an urban area may manifest in conspicuous ways, such as in murky or smelly streamwater, or in less conspicuous ways, such as fewer native or pollution-sensitive organisms. In 2004, and again in 2007, the U.S. Geological Survey sampled stream organisms—algae, invertebrates, and fish—in 14 Milwaukee area streams to assess water quality as part of the ongoing Milwaukee Metropolitan Sewerage District (MMSD) Corridor Study. In addition, passive-sampling devices (SPMDs, “semipermeable membrane devices”) were deployed at a subset of sites in order to evaluate the potential exposure of stream organisms to certain toxic chemicals. Results of the 2007 sampling effort are the focus of this report. Results of sampling from 2007 are compared with results from 2004. The water quality of sampled streams was assessed by evaluating biological-assemblage data, metrics computed from assemblage data, and an aggregate bioassessment ranking method that combined data for algae, invertebrates, and fish. These data contain information about the abundance (number) of different species in each group of stream organisms and the balance between species that can or cannot tolerate polluted or disturbed conditions. In 2007, the highest numbers of algal, invertebrate, and fish species were found at the Milwaukee River at Milwaukee, the largest sampled site. Algal results indicated water quality concerns at 10 of the 14 sampled sites due to the occurrence of nuisance algae or low percentages of pollution-sensitive algae. When compared to 2004, total algal biovolume was higher in 2007 at 12 of 14 sites, due mostly to more nuisance green algae from unknown causes. Results of several metrics, including the Hilsenhoff Biotic Index (HBI-10), suggest that invertebrate assemblages in the Little Menomonee River, Underwood Creek, and Honey Creek were poorer quality in 2007 compared to 2004. Six sites received “very poor” quality ratings for fish in 2007, mostly because inadequate numbers of fish were collected at five sites to allow computation of an Index of Biotic Integrity (IBI); this resulted in three additional sites receiving “very poor” ratings compared to 2004. Some signs of potential improvement in the fish assemblage were evident at Lincoln Creek, possibly reflecting delayed effects of the restoration of stream habitat, completed in 2002; however, algae and invertebrates did not show signs of improvement. Aggregate bioassessment rankings across all groups of organisms for 2004 and 2007 indicated that water quality at the two Milwaukee River main stem sites (at Milwaukee and near Cedarburg), Jewel Creek, and the Menomonee River at Menomonee Falls was the least-degraded among all sampled sites. Rankings for Oak Creek and Little Menomonee suggested water quality was worse in 2007 compared to 2004 and placed these two sites together with Kinnickinnic River and Underwood Creek, two concrete-line sites, indicating the most-degraded water quality among all sampled sites. The aggregate ranking for Lincoln Creek in 2007 would have placed it in the most-degraded category but for the positive influence of the fish ranking when compared to poor algal and invertebrate rankings. Potential toxicity due to certain manmade chemicals, such as polycyclic aromatic hydrocarbons (PAHs), was found at all six sites where SPMDs were deployed. As was found in 2004, the highest potential toxicity in 2007 was observed at Lincoln Creek where chemical screening in 2007 also showed the highest total PAHs of all six sites; however, potential toxicity at Little Menomonee River, Honey Creek, and Kinnickinnic River was relatively high compared to Milwaukee River near Cedarburg. Although toxicity and chemical results in 2007 did not agree with aggregate rankings for Lincoln Creek because of fish, nor for Honey Creek, the results did agree with aggregate rankings at four of the six sites. In addition to toxicological and chemical influences, the more urbanized sites have high percentages of impervious surface area, resulting in frequent high stream flows that can adversely affect algal, invertebrate, and fish assemblages. Assessments of the ecological status of different groups of organisms and of potential chemical and physical stressors to organisms are important tools in evaluating streamwater quality.

Wisconsin↗

USGS2013

Oceanographic controls on sedimentary and geochemical facies on the Peru outer shelf and upper slope

Concentrations and characteristics of organic matter in surface sediments deposited under an intense oxygen-minimum zone (OMZ) on the Peru margin were mapped and studied in samples from deck-deployed box cores and push cores acquired by submersible on two east-west transects spanning depths of 75 to 1,000 meters (m) at 12°S and 13.5°S. On the basis of sampling and analyses of the top 1–2 centimeters (cm) of available cores, three main belts of sediments were identified on each transect with increasing depth: (1) muds rich in organic carbon (OC); (2) authigenic phosphatic mineral crusts and pavements; and (3) glaucony facies. Sediments rich in OC on the 12°S transect were mainly located on the outer shelf and upper slope (150–350 m), but they occurred in much shallower water (approximately 100 m) on the 13.5°S transect. The organic matter is almost entirely marine as confirmed by Rock-Eval pyrolysis and isotopic composition of OC. Concentrations of OC are highest (up to 18 percent) in sediments within the OMZ where dissolved oxygen (DO) concentrations are <5 micromoles per kilogram (μM). Even at these low concentrations of DO, however, the surface sediments from within the OMZ are dominantly unlaminated. Concentrations of DO may have the dominant effect on organic matter characteristics, but reworking of fine-grained sediment and organic matter by strong bottom currents with velocities as high as 30 centimeters per second (cm/s) on the slope between 150 and 300 m and redeposition on the seafloor in areas of lower energy and higher DO concentration also exert important controls on OC concentration and degree of oxidation in this region. Phosphate-rich sediments and crusts occurred at depths of about 300 to 550 m on both transects. Nodular crusts of sediment cemented by carbonate-fluorapatite (CFA; phosphorite) or dolomite form within the OMZ. These phosphorite crusts evolve through cementation from light olive-green, stiff but friable, phosphatized claystone “protocrusts” through dense, dark phosphorite crusts, cemented breccias, and pavements. The degree of phosphatization and thickness of the crusts depend on the rates of sediment supply and on the strength and frequency of currents that re-expose crusts on the seafloor. Phosphorite crusts and pavements on the Peru margin can only become buried and incorporated into the geologic record once bottom currents slacken sufficiently to allow fine-grained sediment to accumulate. Glaucony-rich surface sediments, relatively undiluted by other components, were found mainly in deeper water on the 13.5°S transect (750 m to at least 1,067 m). These sediments consist almost entirely of sand-size glaucony pellets. These widespread glaucony sands formed in place and were then concentrated and reworked by strong currents that winnowed away the fine-grained matrix. Although the glaucony occurs in sand-size pellets, the pellets are made up of aggregates of authigenic, platy, micaceous clay minerals. Glaucony is predominantly a potassium (K), sodium (Na), iron (Fe), magnesium (Mg) aluminosilicate with an approximate formula of (K,Na)(Fe 3+ ,Al,Mg) 2 (Si,Al) 4 O 10 (OH) 2 . The glaucony on the 13.5°S transect forms by alteration of one or more original “framework” minerals (carbonate and [or] aluminosilicates) to form pellital aggregates of Fe-, K-, and Mg-rich clay minerals. Because Fe, K, and Mg are derived from seawater, sedimentation rates must be extremely slow in order for the original framework minerals to remain in contact with seawater. The close association of glaucony and phosphorite indicates a delicate balance between the slightly oxidizing conditions at the base of the OMZ that form glaucony and the slightly reducing conditions that mobilize iron and phosphate to form phosphorite.

Open-File Report↗

USGS2006

Palila Restoration: Lessons from Long-term Research

BACKGROUND The palila (Loxioides bailleui) is a member of the Hawaiian honeycreeper family of birds (Drepanidinae), which is renowned for the profusion of species - many with bizarre bills and specialized feeding habits - that radiated from a single ancestral type. Most of the 57 or so honeycreeper species are extinct, and the palila is endangered because of its high degree of dependence on the mamane tree (Sophora chrysophylla) (Figure 1) and its restricted distribution on the upper slopes of Mauna Kea (Figure 2). Three decades of research have revealed many important facts about palila, providing the foundation and impetus for conservation programs in the wild and captivity. Additionally, an ambitious public conservation campaign arose due to the land-use conflicts on Mauna Kea. Here we summarize progress in palila conservation biology and outline steps that might overcome the remaining major challenges to its recovery. We also highlight lessons learned from palila research that may help the recovery of other Hawaiian forest birds. Palila and two closely-related species on the tiny islands of Nihoa and Laysan are the last of the seed-eating honeycreeper species in the Hawaiian Islands. About a quarter of the honeycreeper species known from living and fossil specimens had finch-like bills suited mainly for eating seeds and fruits. Because of their dietary specialization, palila are vulnerable to changes in forest size and quality, as was also likely the case for extinct species of seed specialists. Palila and many other forest bird species were once distributed in dry, lowland forests. Fossil records indicate that palila also occurred in the lowlands of O`ahu and Kaua`i until human settlement of those islands. However, because lowland habitats have been highly modified by humans and because mamane occurs today primarily at high elevation, palila are the only native bird species found exclusively in dry, subalpine habitat (2000?2850 m). Similar to other feeding specialists, palila lay few eggs, raise few young each year, and take a relatively long time to complete the nesting cycle. Low rates of reproduction result in low rates of population growth and low potential for recovery from disturbances. Long-term studies of palila offer important insights into the conservation biology of all Hawaiian forest bird species, particularly feeding specialists like the palila. Palila face many challenges common to both generalist and specialist Hawaiian honeycreeper species. Habitat loss and degradation, as well as introduced avian diseases, have reduced their numbers and limited their distribution to a very small portion of their historic range. Introduced mammals prey on palila, while alien insects reduce caterpillars that are particularly important in the diet of nestlings. Securing legal protection and funding for palila restoration has been challenging. Understanding how the palila has avoided extinction can help managers plan its recovery, and better design recovery plans for species with different feeding strategies in other habitats.

Fact Sheet↗

USGS1999

Chapter 5: Cumulate nodules in the Thera pyroclastics

Gabbroic and dioritic nodules are common in some horizons of the Thera pyroclastics. These nodules contain pyroxene as the main ferromagnesian phase in contrast to the early Akrotiri centre, which contain nodules with abundant hornblende (Nicholls 1971a; Pichler & Kussmaul 1972). Calcsilicate nodules and gabbroic nodules have also been recognized in Santorini lavas (Fouqu6 1879; Nicholls 1971a). Previous studies of Santorini volcanic rocks indicate that crystal fractionation has played an important role in generating the compositional range basalt-andesite-dacite-rhyodacite (Nicholls 1971a; Mann 1983; Huijsmans 1985; Barton & Huijsmans 1986). This interpretation is also consistent with much of the petrological data presented in Chapter 4 and the geochemical data and models presented in Chapter 6. However, disequilibrium phenocryst assemblages in some lavas (Huijsmans 1985; Chapter 4), banded pumice clasts in many pyroclastic deposits (Druitt et al. 1989; Chapters 4), isotopic disequilibria in phenocrysts (Pyle et al. 1988; Chapters 6 and 7) and complex zoning patterns in plagioclase crystals (Stamatelopoulou-Seymour et al. 1990) indicate that magma mixing has also occurred. In addition, trace element and isotopic abundances can be interpreted in terms of assimilation of continental crust during fractionation (Barton et al. 1983; Stamatelopoulou-Seymour et al. 1990; Chapter 6). The nodules provide information on petrogenetic processes that have influenced the evolution of Santorini magmas. Their petrology is documented in this chapter and compared with that of associated juvenile ejecta. The nodules are identified as cumulates from andesitic and dacitic magmas, con- firming the importance of fractional crystallization in the petrogenesis of Santorini magmas. They also demonstrate that mafic plutonic rocks can be formed from intermediate to silicic magmas.

Santorini↗

USGS2025

Performance analysis of oil recovery and CO2 retention in a greenfield residual oil zone: CO2-EOR in Tall Cotton Field (Permian Basin, West Texas, USA)

Residual oil zones (ROZs) can offer significant oil resources via enhanced oil recovery (EOR) as well as subsurface carbon dioxide (CO 2 ) retention during injection. If injected CO 2 is anthropogenic, the ROZs can offer a substantial geologic storage potential. The ROZs below the oil/water contact (OWC) of main pay zones (MPZ) in conventional reservoirs or brownfields, are more commonly developed for CO 2 injection and oil production and reported in the literature. However, CO 2 -EOR in greenfield ROZs, reservoirs without a MPZ present, have rarely been developed for CO 2 -EOR operation. The Tall Cotton Field of West Texas, Permian Basin, which started production in 2015 (Phase 1) and expanded in 2017 (Phase 2) from the San Andres Limestone, is one of the first examples of greenfield ROZs developed for EOR by injecting CO 2 . This paper analyses EOR and CO 2 retention performance of Tall Cotton Field using allocated injection and production data from inverted 5-spot well patterns of Phase-1 and -2 developments. Production and injection data allocated to each of the 28 identified patterns (nine 20-acre patterns for Phase-1, three 20-acre and sixteen 10-acre patterns for Phase-2) were analyzed for historical and forecasted oil recovery using ratio-trend decline analysis, and for CO 2 retention performance of the patterns. The allocated data were further used to calculate injected reservoir pore volume and void replacement ratios (VRR) for the analysis period. Quantitative results indicated that oil recovery factors of the 5-spot patterns varied between 4–10 %, and 5–30 % between the end of injection and the forecast periods, respectively. Storage of CO 2 , on the other hand, increased to a mean value of ∼7130 MMscf per pattern in Phase-1 and to a mean storage of 3700 MMscf per pattern in Phase-2 until the end of injection, followed by a decline after the end of injection and into the forecast period. Resulting CO 2 utilization factors ∼6–50 Mscf/bbl were estimated at the end of injection. Overall, presented results suggested that developing greenfield ROZs for CO 2 -EOR can be as promising as brownfield ROZs and mature MPZs for EOR and underground storage of injected CO 2 . For Tall Cotton Field, results suggest that Phase-2 patterns generally outperformed Phase-1 for oil recovery factors, while Phase-1 performed better in CO 2 retention performance metrics. This is the first study in the literature that reports a detailed CO 2 -EOR performance analysis of a greenfield ROZ in the Permian Basin, which can potentially allow for comparison with MPZs and brownfield ROZs.

Texas↗

USGS1989

Morphology and stratigraphy of small barrier-lagoon systems in Maine

The coast of Maine contains over 200 individual barrier-lagoon systems, most quite small, with an aggregate length of nearly 100 km. Although they represent less than 5% of the tidally influenced coastline of Maine, they are widely distributed and occur in a variety of dynamic regimes and physiographic regions. Their morphology and backbarrier stratigraphy are different from better studied coastal plain systems, and provide important clues to the Holocene evolution of the Maine coast. In a study of geomorphic form and backbarrier stratigraphy, inlet processes and Holocene sea-level rise have been identified as the principal controls on coarse-grained barrier stratigraphy. Barriers in Maine are found in five distinct geomorphic forms, identified herein as: barrier spits, pocket barriers, double tombolos, cuspate barriers and looped barriers. The few long sandy beaches in southwestern Maine are mostly barrier spits. The remainder of the barrier types is composed primarily of gravel or mixed sand and gravel. The barriers protect a variety of backbarrier environments: fresh and brackish ponds, lagoons and fresh- and saltwater marshes. The barriers may or may not have inlets. Normal wave action, coarse-grain size and a deeply embayed coast result in barriers with steep, reflective profiles several meters above MHW. Occasional storm events completely wash over the barriers, building steep, lobate gravel fans along their landward margin. Few, if any, extensive storm layers are recognized as extending into the distal backbarrier environments, however. During sea-level rise and landward barrier retreat, this abrupt, storm-generated transition zone inters the backbarrier sediments. Statistical comparisons of barrier morphology, location and backbarrier environment type with backbarrier stratigraphy show that Holocene backbarrier stratigraphy is best predicted by the modern backbarrier environment type. This, in turn, is influenced most by the absence or presence, and long-term stability or instability of a tidal inlet. Geomorphic barrier form and location in coastal geomorphic compartments show little or no correlation with backbarrier stratigraphy. In contrast to previous classifications of barrier-lagoon systems based primarily on sandy, coastal plain examples, in Maine the shape or origin of the backbarrier system is relatively unimportant. The presence or absence of a tidal inlet is of paramount importance in shaping the Holocene stratigraphy of the backbarrier region.

Marine Geology↗

USGS2004

Twentieth century demographic changes in cirio and cardón in Baja California, Mexico

Aim&ensp; Our purpose was to discern long-term large-scale patterns of survivorship and recruitment of two dominant, charismatic and useful desert plants, cirio [Fouquieriaceae; Fouquieria columnaris (Kell.) Kell. ex Curran] and cardón [Cactaceae; Pachycereus pringlei (S.Watson) Britton & Rose], and to evaluate the effects of physical and human factors. Location&ensp; The study included 77 sites distributed over c. 13,200 km 2 in the northern Vizcaíno Region of the Sonoran Desert, state of Baja California, México. Method&ensp; Censuses used identified individuals ( n = 3780 cirio and 2246 cardón) and were based on repeat photography. Time intervals between censuses ranged from 29 to 96 years, ending in 1996–2002. Earlier repetition (1963) of old photographs at 14 sites allowed analysis of temporal variation. The effect of elapsed time was modelled with Weibull functions for survivorship and linear functions for recruitment. To explain the residual variations, we tested categorical and continuous variables related to substrate, geomorphology, climate and geography. Human impacts were weighed with historical and quantitative analysis. Results&ensp; Elapsed time accounted for 45% of the variation of survivorship in cirio but only 35% in cardón. The fitted Weibull curves were used to estimate longevity: all-size cohorts of 1000 individuals would die out in 388 (cirio) and 723 years (cardón). Recruitment in cirio was linearly related to time ( r 2 = 0.71) and averaged < 0.003 year −1 but was less than mortality at 90% of the sites. For cardón, recruitment averaged > 0.005 year −1 , was linearly related to time ( r 2 = 0.66) and was less than mortality at only 33% of the sites. Vital rates were not strongly correlated between mid- and late-century but were mostly similar to the long-term regional rates. Removing the factor of elapsed time, the residual variations showed that survivorship was greater on sites protected from winds for cirio and was less on very coarse soils for cardón. Recruitment increased with latitude in cirio and was greater on southern exposures for both species. The residual variations were not clearly affected, at a regional scale, by other factors thought to be important, including elevation, distance to the Pacific Ocean, geology, slope gradient, soil stability, older vs. young alluvial soils, and soil Ca : Mg and Na : K ratios. Human impacts have been sporadic and heterogeneous but locally strong; our quantitative indices of accessibility did not show regionally significant effects. Blowdown by hurricanes is a sporadic cause of major mortality for cirio but not necessarily for cardón. Main conclusions&ensp; At our scales, effects of time outweigh those of substrate, and human impact was scant. Large patterns were pronounced: cirio experienced a prolonged, widespread decline in the 1900s, while cardón fluctuated in different directions and degree among local populations. Cirio was found to be inherently much shorter-lived than cardón. We also suggest that recruitment in cirio was low due to long periods with relatively dry winters that affected the entire region, while spatial heterogeneity of cardón recruitment seemed more related to the variation in summer rains.

Journal of Biogeography↗

USGS2015

California State Waters map series — Offshore of Scott Creek, California

Introduction 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 the 3-nautical-mile limit of 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 subsurface geology. The Offshore of Scott Creek map area is located in central California, on the Pacific Coast about 65 km south of San Francisco and 12 km northwest of Santa Cruz. The onshore part of the map area is sparsely populated; the only cultural center is Davenport, a small community with a population of less than 500. The hilly coastal area is virtually undeveloped, and a large percentage of coastal land is incorporated in open-space trusts. Agricultural land is almost entirely limited to coastal areas between the shoreline and the northwest-trending Santa Cruz Mountains, on Pleistocene alluvial fan deposits and the lowest emergent marine terrace. The Santa Cruz Mountains are part of the northwest-trending Coast Ranges that run roughly parallel to the San Andreas Fault Zone. The map area is cut by the San Gregorio Fault Zone, and it lies a few kilometers southwest of the San Andreas Fault Zone. Regional folding and uplift along the coast has been attributed to a westward bend in the San Andreas Fault Zone and also to right-lateral movement along the San Gregorio Fault Zone. The irregular coastal geomorphology of this area, which consists of low, rocky cliffs and sparse, small pocket beaches backed by low, terraced hills, is partly attributable to this ongoing deformation. The shelf in the map area is underlain by variable amounts (0 to 25 m) of upper Quaternary shelf, nearshore, and fluvial sediments deposited as sea level fluctuated in the late Pleistocene. The northernmost part of the map area is characterized by the presence of uplifted bedrock that has been linked to a local transpressional zone in the San Gregorio Fault Zone. This uplift, coupled with high wave energy, has resulted in little or no sediment cover in this area where exposures of bedrock are present at water depths of as much as 45 m. The thickest deposits of sediment lie offshore of both Davenport and the mouth of Waddell Creek. Coastal sediment transport in the map area is characterized by north-to-south littoral transport of sediment that is derived mainly from streams in the Santa Cruz Mountains and also from local coastal erosion. Shoreline-change studies indicate long-term erosion; within the region between San Francisco and Davenport, the highest long- and short-term coastal-erosion rates occur north of the map area, just north of Point Año Nuevo. During the last approximately 300 years, as much as 18 million cubic yards (14 million cubic meters) of sand-sized sediment has been eroded from the area between Año Nuevo Island and Point Año Nuevo and transported south. Once widened by this pulse of eroded sediment, beaches in the map area are now narrowing as the tail end of this mass of sand progresses farther south. The Offshore of Scott Creek map area lies within the cold-temperate biogeographic zone that is called either the “Oregonian province” or the “northern California ecoregion.” This biogeographic province is maintained by the long-term stability of the southward-flowing California Current, the eastern limb of the North Pacific subtropical gyre that flows from southern British Columbia to Baja California. At its midpoint off central California, the California Current transports subarctic surface (0–500 m deep) waters southward, about 150 to 1,300 km from shore. Seasonal northwesterly winds that are, in part, responsible for the California Current, generate coastal upwelling. The south end of the Oregonian province is at Point Conception (about 320 km south of the map area), although its associated phylogeographic group of marine fauna may extend beyond to the area offshore of Los Angeles in southern California. The ocean off of central California has experienced a warming over the last 50 years that is driving an ecosystem shift away from the productive subarctic regime towards a depopulated subtropical environment. Seafloor habitats in the Offshore of Scott Creek map area, which lie within the Shelf (continental shelf) megahabitat, range from significant rocky outcrops that support kelp-forest communities nearshore to rocky-reef communities in deeper water. Biological productivity resulting from coastal upwelling supports populations of Sooty Shearwater, Western Gull, Common Murre, Cassin’s Auklet, and many other less populous bird species. In addition, an observable recovery of Humpback and Blue Whales has occurred in the area; both species are dependent on coastal upwelling to provide nutrients. The large extent of exposed inner shelf bedrock supports large forests of “bull kelp,” which is well adapted for high-wave-energy environments. The kelp beds are the northernmost known habitat for the population of southern sea otters. Common fish species found in the kelp beds and rocky reefs include lingcod and various species of rockfish and greenling.

California↗

USGS2016

Regional chloride distribution in the Northern Atlantic Coastal Plain aquifer system from Long Island, New York, to North Carolina

The aquifers of the Northern Atlantic Coastal Plain are the principal source of water supply for the region’s nearly 20 million residents. Water quality and water levels in the aquifers, and maintenance of streamflow, are of concern because of the use of this natural resource for water supply and because of the possible effects of climate change and changes in land use on groundwater. The long-term sustainability of this natural resource is a concern at the local community scale, as well as at a regional scale, across state boundaries. In 2010, the U.S. Geological Survey (USGS) began a regional assessment of the Northern Atlantic Coastal Plain aquifers. An important part of this assessment is a regional interpretation of the extent of saltwater and the proximity of saltwater to fresh-groundwater resources and includes samples and published interpretations of chloride concentrations newly available since the last regional chloride assessment in 1989. This updated assessment also includes consideration of chloride samples and refined interpretations that stem from the 1994 discovery of the buried 35 million year old Chesapeake Bay impact structure that has substantially altered the understanding of the hydrogeologic framework and saltwater distribution in eastern Virginia. In this study, the regional area of concern for the chloride samples and interpretations extends from the Fall Line in the west to the outer edge of the Continental Shelf in the east and from the eastern tip of Long Island in the north to about halfway down the North Carolina coast in the south. Discussions of chloride distribution are presented for each of the 10 regional aquifer layers of the Northern Atlantic Coastal Plain, including the offshore extents. Maps of interpreted lines of equal concentration or isochlors were manually prepared for nine of the regional aquifers; a map was not prepared for the surficial regional aquifer. The isochlor interpretations include the offshore extent of the nine regional aquifers and are presented on a 1:2,000,000 scale base map. Vertically, the chloride samples and interpretations range from deepest (oldest) to shallowest (youngest)—Potomac-Patuxent, Potomac-Patapsco, Magothy, Matawan, Monmouth-Mount Laurel, Aquia, Piney Point, Lower Chesapeake, and Upper Chesapeake regional aquifers. The approach of this study maximizes the overall density of chloride information and data by assessing relevant published interpretations, all USGS chloride samples, and all relevant offshore samples in one comprehensive interpretation. Published isochlors, where they were interpreted by regional aquifer, were used as much as possible for this regional isochlor assessment. Publication dates for the isochlors used range from 1982 to 2015, and the scales for the isochlors range from local (county or municipality) to state (sub-regional) to regional. The USGS National Water Information System database provided well sample data for the parts of aquifers that are mainly beneath the land areas and yielded 37,517 water-quality records for 1903 through 2011. Published data reports from four phases of research-related offshore coring (1976, 1993, 1997, 2009) were the main source of water-quality data for the parts of aquifers from the shoreline to the outer edge of the Continental Shelf and yielded samples from multiple depths of each of 13 cores. This study also used interpretations and offshore core data from the last regional chloride assessment (1989) which, in addition to 7 offshore cores, included water-quality data from about 500 wells, and borehole geophysics interpretations from a subset of 11 wells. All published information and data that were used in this study were considered time independent and did not assess the published interpretations or data for temporal trends. The approach used here examined only published interpretations and available chloride data, and did not directly use supplemental techniques that can provide insight into the distribution of saltwater, such as geochemical characterization, borehole geophysical information, and geochronology. Isochlor maps for this study are limited to manual interpretations of the 250-milligram per liter (mg/L) and 10,000-mg/L boundaries developed for 9 of the 10 regional aquifers that constitute the regional hydrogeologic framework of the Northern Atlantic Coastal Plain. For a given aquifer, the approach was to initially consider published isochlor interpretations, where available, then to modify the published interpretations, if necessary, to the extent indicated by the well and core samples. The final step was to interpolate isochlors to the full extent of each aquifer layer in areas with sufficient samples or cited interpretations, or to extrapolate isochlors in areas with no samples or where samples were sparse. The principal limitation of this study is that, because of its regional extent, data and information density can vary greatly, and thus confidence in interpretations can vary widely for onshore and offshore areas across the study area. In areas of sparse data, some samples of elevated chloride could be misinterpreted as being part of a regional elevated chloride trend, and in other cases, an elevated concentration could be misinterpreted as being of only local importance. The interpretive work of this study was applied to a 1:2,000,000 scale base map. Locations of isochlors, wells, cores, political boundaries, and shorelines are meant to be considered approximate. The isochlors presented in this study were manually interpreted for each aquifer unit as a conceptual representation of an equal concentration line approximately in the middle of an aquifer’s thickness. Differences in chloride concentration lines between the top and bottom of an aquifer could be substantial, especially for the thick parts of aquifers, but that information is not presented in this regional assessment. Although additional offshore chloride data are available compared to 27 years ago (1989), the offshore information remains sparse, resulting in less confidence in the offshore interpretations than in the onshore interpretations. Regionally, the 250- and 10,000-mg/L isochlors tend to map progressively eastward from the deepest to the shallowest aquifers across the Northern Atlantic Coastal Plain aquifer system but with some exceptions. The additional data, conceptual understanding, and interpretations in the vicinity of the buried Chesapeake Bay impact structure in eastern Virginia resulted in substantial refinement of isochlors in that area. Overall, the interpretations in this study are updates of the previous regional study from 1989 but do not comprise major differences in interpretation and do not indicate regional movement of the freshwater-saltwater interface since then.

Delaware, Maryland, New Jersey, New York, North Ca↗

USGS2010

Effects of water use and land use on streamflow and aquatic habitat in the Sudbury and Assabet River Basins, Massachusetts

Water withdrawals from surface-water reservoirs and groundwater have affected streamflow in the Sudbury and Assabet River Basins. These effects are particularly evident in the upper Sudbury River Basin, which prompted the need to improve the understanding of water resources and aquatic habitat in these basins. In 2004, the U.S. Geological Survey, in cooperation with the Massachusetts Department of Conservation and Recreation, developed a precipitation-runoff model that uses Hydrologic Simulation Program-FORTRAN (HSPF) to evaluate the effects of water use and projected future water-use and land-use change on streamflow. As part of this study, the aquatic habitat in the basins and the effects of streamflow alteration also were evaluated. Chapter 1 of the report covers the development of the HSPF model that focuses on the upper Sudbury River Basin (106 square miles) but covers the entire Sudbury and Assabet River Basins (339 square miles). The model was calibrated to an 11-year period (1993-2003) using observed or estimated streamflow at four streamgages. The model was then used to simulate long-term (1960-2004) streamflows to evaluate the effects of average 1993-2003 water use and projected 2030 water-use and land-use change over long-term climatic conditions. Simulations indicate that the average 1993-2003 withdrawals most altered streamflow relative to no withdrawals in small headwater subbasins where the ratios of mean annual withdrawals to mean annual streamflow are the highest. The effects of withdrawals are also appreciable in other parts of the upper Sudbury River Basin as a result of the perpetuation of the effects of large withdrawals in upstream reaches or in subbasins that also have a high ratio of withdrawal to streamflow. The simulated effects of potential 2030 water-use and land-use change indicate small decreases in flows as a result of increased water demands, but these flow alterations were offset as a result of decreased evapotranspiration associated with the loss of deep-rooted vegetation. Simulations of reactivating production wells near the north end of Lake Cochituate indicate pumping could substantially affect lake levels and flows at the lake outlet or in nearby reaches in the Sudbury River during periods of low flow, but the effects vary depending on the source of the water to the wells, which is largely unknown. Chapter 2 of the report covers the fish-community assessment and comparison of streamflow-setting standards for protecting aquatic habitat. The fish-community assessment indicates the main stems of the Sudbury and Assabet Rivers are dominated by macrohabitat generalists. Water temperatures recorded in seven free-flowing reaches in the upper Sudbury River Basin at three sites unaffected by withdrawals or impoundments are generally suitable for cold-water fish; however, summer temperatures often rose to a level considered critical to long-term survival of brook trout. At four sites downstream from withdrawals or reservoirs, or both, summer water temperatures were often in the upper critical range for brook trout survival. Physically and statistically based methods for determining streamflows for protecting aquatic habitat were applied at 10 selected riffle sites in the Sudbury and Assabet River Basins. Physically based methods, R2Cross and Wetted-Perimeter, use site-specific physical and hydraulic information and a one-dimensional hydraulics model, HEC-RAS, to determine flows that meet the criteria set forth by the method. The median flow that meets 2-of-3 of the R2Cross hydraulic criteria (percentage of bankfull wetted perimeter, average velocity, and mean depth) ranged from about 0.07 to 0.72 cubic feet per second per square mile (ft3/s/mi2) with an overall median of about 0.24 ft3/s/mi2; the median Wetted-Perimeter target flow ranged from about 0.10 to 0.51 ft3/s/mi2 with an overall median of about 0.25 ft3/s/mi2. Statistically based methods?Tennant, New England Aquatic Base Flow (ABF)

Scientific Investigations Report↗

USGS1983

Stratigraphy of the Silurian outcrop belt on the east side of the Cincinnati Arch in Kentucky, with revisions in the nomenclature

Silurian rocks form a narrow arcuate outcrop belt about 100 mi long on the east side of the Cincinnati Arch in Kentucky. They range from as much as 300 ft thick in the north to a pinchout edge in the south. The nomenclature of this sequence is revised to reflect mappability and lithologic uniformity on the basis of detailed mapping at a scale of 1:24,000 by the U.S. Geological Survey in cooperation with the Kentucky Geological Survey. The Silurian rocks are divided into two parts: the Crab Orchard Group, raised in rank from Crab Orchard Formation and redefined, in the lower part of the Silurian section, and Bisher Dolomite in the upper part of the section. The Crab Orchard Group is subdivided into the Drowning Creek Formation (new name) at the base of the Silurian, overlain by the Alger Shale (adopted herein) south of Fleming County and by the Estill Shale (elevated to formational rank) north of Bath County. The Brassfield Member (reduced in rank from Brassfield Dolomite or Formation) and the Plum Creek Shale and Oldham Members of the former Crab Orchard Formation are included as members of the Drowning Creek; the Lulbegrud Shale, Waco, and Estill Shale Members of the former Crab Orchard Formation are now included in the Alger. The Drowning Creek Formation, 20 to 50 ft thick, is composed mainly of gray fine to coarse-grained dolomite with shale interbeds. The dolomite beds average several inches thick, with bedding surfaces that are locally smooth but generally irregular and are fossiliferous in many places; fossils include brachiopods, crinoid columnals, horn corals, colonial corals, trilobites, pelecypods, and bryozoans. The shale interbeds average several inches thick, except for its Plum Creek Shale Member which is entirely shale and as much as 12 ft thick, and are most abundant in the upper half of the formation. The members of the Drowning Creek intergrade and are indistinguishable in the northern part of the area. The Alger Shale, as much as 170 feet thick, is predominantly grayish-green clay shale with a thin (0.5-3 ft) dolomite member (the Waco, or its northern equivalent, the Dayton Dolomite Member, reduced in rank from Dayton Limestone) near the base. North of Bath County, the Lulbegrud Shale and Dayton Dolomite Members are reassigned to the underlying Drowning Creek Formation, the Estill Shale Member is elevated to formational status, and the Alger is dropped. The Bisher Dolomite, which overlies the Estill Shale in the northernmost part of the Silurian belt, ranges from 0 to 300 ft in thickness and is composed of medium-to coarse-grained, gray, fossiliferous dolomite. The Silurian section overlies Upper Ordovician rocks in apparent conformity, although faunal studies suggest a minor hiatus, and is overlain by Middle to Upper Devonian rocks in a regional angular unconformity that truncates the entire Silurian section at the southwest end of the outcrop belt, where it is nearest the axis of the Cincinnati Arch. All of the units recognized in the Silurian appear to thicken eastward, away from the axis of the arch and towards the Appalachian basin. This, with the presence of isolated remnants of the Brassfield near the axis, suggest that formation of the arch was initiated in Early Silurian time by subsidence of its eastern flank.

Professional Paper↗

USGS1984

Control of barrier island shape by inlet sediment bypassing: East Frisian Islands, West Germany

A study of the East Frisian Islands has shown that the plan form of these islands can be explained by processes of inlet sediment bypassing. This island chain is located on a high wave energy, high tide range shoreline where the average deep-water significant wave height exceeds 1.0 m and the spring tidal range varies from 2.7 m at Juist to 2.9 m at Wangerooge. An abundant sediment supply and a strong eastward component of wave power (4.4 × 10 3 W m −1 ) have caused a persistent eastward growth of the barrier islands. The eastward extension of the barriers has been accommodated more by inlet narrowing, than by inlet migration. It is estimated from morphological evidence that a minimum of 2.7 × 10 5 m 3 of sand is delivered to the inlets each year via the easterly longshore transport system. Much of this sand ultimately bypasses the inlets in the form of large, migrating swash bars. The location where the swash bars attach to the beach is controlled by the amount of overlap of the ebb-tidal delta along the downdrift inlet shoreline. The configuration of the ebb-tidal delta, in turn, is a function of inlet size and position of the main ebb channel. The swash bar welding process has caused preferential beach nourishment and historical shoreline progradation. Along the East Frisian Islands this process has produced barrier islands with humpbacked, bulbous updrift and bulbous downdrift shapes. The model of barrier island development presented in this paper not only explains well the configuration of the German barriers but also the morphology of barriers along many other mixed energy coasts.

East Frisian Islands↗

USGS2010

Vegetation classification and distribution mapping report: Canyon de Chelly National Monument

Executive Summary: The classification and distribution mapping of the vegetation of Canyon de Chelly National Monument (CACH) and surrounding environment was accomplished through a multi-agency effort between 2003 and 2007. The National Park Service’s Southern Colorado Plateau Network facilitated the team that conducted the work, which comprised the U.S. Geological Survey’s Southwest Biological Science Center and Fort Collins Science Center, Navajo Natural Heritage Program, Northern Arizona University, and NatureServe. The project team described 48 plant communities for CACH—35 of which were described from quantitative classification based on field-relevé data collected in 2004. Five additional plant communities were based on field relevés collected in a previous study. The team derived four additional plant communities from field observations during the photointerpretation phase of the project, and field documented them during accuracy assessment. The National Vegetation Classification Standard served as a conceptual framework for assigning these plant communities to the alliance and association level. Ten of the 48 plant communities were designated “park specials”, that is, plant communities with insufficient data to describe them as new alliances or associations. The project team also developed a spatial vegetation map database representing CACH, with three different map-class schemas: base, group, and management map classes. The base map classes represented the finest level of spatial detail. Photointerpreters delineated initial polygons through manual interpretation of 2003/2004 1:12,000-scale true color aerial photography supplemented by occasional computer screen digitizing on a mosaic of digitized aerial photos. These polygons were labeled with base map classes during photointerpretation. Field visits verified interpretation concepts. The vegetation map database includes • ? 53 base map classes, which consist of associations and park specials classified with the quantitative analysis • ? additional associations noted during photointerpretation • ? non-vegetated land cover, such as infrastructure, land use, and geological land cover. The base map classes consist of 4,718 polygons in the project area. A field-based accuracy assessment of the base map classes showed the overall accuracy to be 50.8% The group map classes represent aggregations of the base map classes, approximating the group level of the National Vegetation Classification Standard, Version 2 (Federal Geographic Data Committee 2008). Terrestrial ecological systems, as described by NatureServe (Comer et al. 2003), were used as a first approximation of the group level. The project team identified 16 group map classes in this project. The overall accuracy of the group map classes was determined using the same accuracy assessment data as for the base map classes. The overall accuracy of the group representation of vegetation was 79.9%. In consultation with park staff, the team developed management map classes that consisted of park-defined groupings of base map classes and were intended to represent a balance between maintaining required accuracy and providing a focus on vegetation of particular interest or import to park managers. The 28 management map classes have an overall accuracy of 77.1%. While the main products of this project are the vegetation classification and the vegetation map database, a number of ancillary geographic information system and digital database products were also produced that can be used independently, or to augment the main products. These products include shapefiles of the location of field-collected data and relational databases of field-collected data.

Report↗

USGS2005

Environmental setting, water budget, and stream assessment for the Broad Run watershed, Chester County, Pennsylvania

The Broad Run watershed lies almost entirely in West Bradford Township, Chester County, Pa., and drains 7.08 square miles to the West Branch Brandywine Creek. Because of the potential effect of encroaching development and other stresses on the Broad Run watershed, West Bradford Township, the Chester County Water Resources Authority, and the Chester County Health Department entered into a cooperative study with the U.S. Geological Survey to complete an annual water budget and stream assessment of overall conditions. The annual water budget quantified the basic parameters of the hydrologic cycle for the climatic conditions present from April 1, 2003, to March 31, 2004. These water-budget data identified immediate needs and (or) deficits that were present within the hydrologic cycle during that period, if present; however, an annual water budget encompassing a single year does not identify long-term trends. The stream assessment was conducted in two parts and assessed the overall condition of the watershed, an overall assessment of the fluvial-geomorphic conditions within the watershed and an overall assessment of the stream-quality conditions. The data collected will document present (2004) conditions and identify potential vulnerabilities to future disturbances. For the annual period from April 1, 2003, to March 31, 2004, determination of an annual water budget indicated that of the 67.8 inches of precipitation that fell on the Broad Run watershed, 38.8 inches drained by way of streamflow to the West Branch Brandywine Creek. Of this 38.8 inches of streamflow, local-minimum hydrograph separation techniques determined that 7.30 inches originated from direct runoff and 31.5 inches originated from base flow. The remaining precipitation went into ground-water storage (1.71 inches) and was lost to evapotranspiration (27.3 inches). Ground-water recharge for this period-35.2 inches-was based on these values and an estimated ground-water evapotranspiration rate of 2 inches. Assessment of fluvial-geomorphic conditions included large-scale mapping of stream classes within the Broad Run watershed and in-depth study of three representative stream reaches also within the Broad Run watershed. Based on the total distance of all stream reaches classified within the Broad Run watershed, 61 percent were classified as C-class, 14 percent as E-class, 13 percent as B-class, 5 percent as F-class, 4 percent as undifferentiated B- and F-class, 2 percent as G-class, and less than 1 percent as A-class. The map of stream classes indicates that the Broad Run watershed currently has no large-scale areas of stream impairment and that, generally, the stream is not entrenched and the main branch of the Broad Run has an available, functioning flood plain. Smaller tributary streams, however, showed signs of localized entrenchment due to site-specific influences such as natural stream-channel evolution, localized channelization, localized contraction due to road and driveway crossings, and (or) increased localized runoff. For example, one small reach along a tributary channel was observed to become entrenched due to runoff originating from a new housing development. Entrenched stream reaches are merely located by large-scale mapping and require individual assessment to determine potential causes of entrenchment and (or) future restorative actions. Three in-depth geomorphic study sites showed that the Broad Run can currently be considered graded or in a state of dynamic equilibrium. The sites did, however, identify certain vulnerabilities to future changes within the watershed. These vulnerabilities included disruption of the present sediment supply, including both increases and (or) reductions in the current sediment loads within the Broad Run; increases in both magnitude and duration of storm-water runoff; encroachment of development onto present flood-plain areas, and (or) alterations to riparian zones. Assessment of stream-quality conditions includ

Scientific Investigations Report↗

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↗