USGS Science⌕ Search

USGS · pp1419

Geology and hydrogeology of the Caribbean Islands aquifer system of the Commonwealth of Puerto Rico and the U.S. Virgin Islands

Abstract

Poorly lithified to unconsolidated carbonate and clastic sedimentary rocks of Tertiary (Oligocene to Pliocene) and Quaternary (Pleistocene to Holocene) age compose the South Coast aquifer and the North Coast limestone aquifer system of Puerto Rico; poorly lithified to unlithified carbonate rocks of late Tertiary (early Miocene to Pliocene) age make up the Kingshill aquifer of St. Croix, U.S. Virgin Islands. The South Coast aquifer, North Coast limestone aquifer system, and Kingshill aquifer are the most areally extensive and function as the major sources of ground water in the U.S. Caribbean Islands Regional Aquifer-System Analysis (CI-RASA) study area. In Puerto Rico's South Coast ground-water province, more than 1,000 meters of clastic and carbonate rocks of Oligocene to Pliocene age infill the South Coast Tertiary Basin. The pattern of lithofacies within this basin appears to have been controlled by changes in base level that were, at times, dominated by tectonic movement (uplift and subsidence), but were also influenced by eustasy. Deposition of the 70-kilometer long and 3- to 8-kilometer wide fan-delta plain that covers much of the South Coast ground-water province occurred largely in response to glacially-induced changes in sea level and climate during the Quaternary period. Tectonic movement played a much less important role during the Quaternary. The North Coast ground-water province of Puerto Rico is underlain by homoclinal coastal plain wedge of carbonate and siliciclastic rocks that infill the North Coast Tertiary Basin and thicken to more than 1,700 meters. A thin basal siliciclastic sequence of late Oligocene age is overlain by a thick section of mostly carbonate rocks of Oligocene to middle Miocene age. Globigerinid limestone of late Miocene to Pliocene age crops out and lies in the shallow subsurface areas of northwestern Puerto Rico. Oligocene to middle Miocene age rocks tentatively can be divided into five depositional sequences and associated systems tracts; these rocks record carbonate and minor siliciclastic deposition that occurred in response to changes in relative sea level. The Cibao Formation represents the most complex of these sequences and contains a varied facies of carbonate, mixed carbonate-siliciclastic, and siliciclastic rocks that reflect differential uplift, subsidence, and transgression of the sea. Uplift, graben formation, and gradual shallowing of the sea are reflected within the bathyal-dominated sedimentary facies of the Kingshill Limestone in St. Croix, U.S. Virgin Islands. Reef-tract limestone beds of Pliocene age were subject to exposure, resubmergence, and meteoric leaching of aragonitic skeletal debris; these beds contain patchy lenses of dolomite that are restricted to a small, structurally-controlled embayment. The South Coast aquifer, the principal water-bearing unit of Puerto Rico's South Coast ground-water province, consists of boulder- to silt-size detritus formed by large and small coalescing fan deltas of Pleistocene to Holocene age. Deep well data indicates that it is possible to vertically separate and group a highly complex and irregular-bedded detrital sequence that underlies distal parts of the fan-delta plain into discrete water-bearing units if correlated with 30- to 40-meter thick, eustatically-controlled depositional cycles. Lithofacies maps show that greatest hydraulic conductivity within the fan-delta plain is generally associated with proximal fan and midfan areas. Distal and interfan areas are least permeable. Alluvial valley aquifers located in the western part of the South Coast ground-water province are important local sources of water supply and appear to contain some of the same physical and hydraulic characteristics as the South Coast aquifer. Older sedimentary rocks within the basin are poor aquifers; conglomeratic beds are well-cemented, and carbonate beds do not contain well-developed solution features, except locally where the beds are overlain by alluvium. Ground-water occurs under unconfined conditions in proximal and midfan areas. Confined conditions within deeper parts of the system and in interfan and some midfan areas are created largely by the intercalated nature of discontinuous fine-grained beds that retard vertical ground-water movement. The development of water resources in southern Puerto Rico has modified the hydrologic system of the South Coast aquifer considerably. Under predevelopment conditions, the South Coast aquifer was recharged in the unconfined, proximal fan and some midfan areas by infrequent rainfall and seepage from streams near the fan apex. Discharge occurred as seabed seepage, baseflow discharge along the lower coastal reach of streams, seepage to coastal wetlands, or evapotranspiration in areas underlain by a shallow water table. Under development conditions, seepage from irrigation canals and areal recharge from furrow irrigation represented a principal mechanism for recharge to the aquifer. Increased ground-water withdrawals in the 1960's and 1970's resulted in declines in the water table to below sea level in some places and intrusion of salt water into the aquifer. By the middle 1980's, a reduction in ground-water withdrawals and a shift from furrow irrigation to drip-irrigation techniques resulted in the recovery of water levels. Under present-day (1986) conditions, regional ground-water flow is coastward but with local movement to some well fields. In addition to the discharge mechanisms described above, ground-water discharges also to coastal canals. The North Coast limestone aquifer system consists of limestone, lesser amounts of dolomite, and minor clastic detritus of Oligocene to Pliocene age that form an unconfined upper aquifer and a confined lower aquifer; these aquifers are separated by a clay, mudstone, and marl confining unit. Topographic relief and incision of carbonate coastal plain rocks by streams are the principal factors controlling the direction of ground-water flow. The North Coast limestone aquifer system is recharged principally by precipitation that enters the upper and lower aquifers where they crop out. Regional groundwater movement from the upper aquifer is to the major rivers, wells, coastal wetlands, coastal, nearshore, and offshore springs, or as seabed seepage. Regional discharge from the lower aquifer is to the major rivers along its unconfined parts or where the confining unit has been breached by streams. Discharge from the lower aquifer also occurs in the San Juan area where the Mucarabones Sand provides an avenue for diffuse upward ground-water flow. Transmissivity within the upper limestone aquifer appears to be largely regulated by the thickness of the freshwater lens. The lens is thickest and transmissivity is greatest in interstream areas that lie in a zone that closely corresponds to the landwardmost extent of the underlying saltwater wedge. Hydraulic conductivity of the upper aquifer generally increases in a coastward direction and reflects lithologic control, karstification in the upper 30 to 100 meters of the section, and enhanced permeability in a zone of freshwater and saltwater mixing. Transmissivity of the lower aquifer is an order of magnitude smaller than that of the upper aquifer; highest transmissivities in the lower aquifer largely correspond to a coarse grainstone-packstone and coral-patch-reef depositional facies contained within the outcropping parts of the Montebello Limestone Member and its subsurface equivalents. Porosity within the North Coast limestone aquifer system is high in grainstone-packstones and low in wackestone and marl. Dolomitized zones and moldic grainstone-packstone strata are the most porous carbonate rocks, but occur in thin beds that usually are only a few meters thick. Processes of karstification that include the development of caverous zones and large vugs, and dissolution along possible regional fracture sets has enhanced permeability within the upper part of the aquifer system. Stratigraphic and lithologic control play an important role controlling permeability within the lower part of the system. The Kingshill aquifer of St. Croix, in large part, is composed of deepwater limestone that contains only microscopic pores and is poorly permeable; however, the upper part of the aquifer, a shallow-water skeletal and reef limestone, is fairly permeable, but restricted in areal extent. Permeability within these uppermost beds of the aquifer has been enhanced by meteoric leaching, dissolution within a mixing zone of saltwater and fresh water, and dolomitization. However, most large-yield wells completed in the Kingshill aquifer are also screened in alluvium that overlies or infills incised channels. The alluvial deposits serve as a temporary storage zone for rainfall, runoff, and ground water slowly entering the Kingshill aquifer.

Explore related subjects

90° N90° S · 180° W ← longitude → 180° E
Source-reported bounding extent: 17.647948051340578° to 18.539512627214105° latitude; -67.29949951171875° to -64.55429077148438° longitude. This indicates report coverage, not an exact sampling location. View area on OpenStreetMap.

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

Robert A. Renken, W. C. Ward, I. P. Gill, Fernando Gómez-Gómez, Jesús Rodríguez-Martínez. 2002. Geology and hydrogeology of the Caribbean Islands aquifer system of the Commonwealth of Puerto Rico and the U.S. Virgin Islands. https://doi.org/10.3133/pp1419

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related USGS reports

Allostratigraphy and formalizing the Cold Harbor Alloformation, Virginia Coastal Plain, Eastern Virginia

Introduction Distinct morphologic features in the Virginia part of the Coastal Plain Province (herein referred to as the Virginia Coastal Plain) were noted in one of the earliest geologic surveys in 1839 (Rogers, 1884). South of the James River, Rogers (1884, p. 252) described a western bench that ranged from “a low rocky coast…of the primeval sea” to an eastern facing scarp east of the City of Suffolk, Virginia. His eastern bench extended from the scarp at Suffolk, which he suggested was a former shoreline, to the present Atlantic Ocean. Notably, the eastern plain was described as less dissected than the western plain. In Virginia, creation and publication of topographic maps by the U. S. Geological Survey began in the mid-1880s. By 1905, many areas in eastern Virginia were mapped as 30-minute and 15-minute series with 10-foot (ft) contour interval accuracy (Virginia Division of Mineral Resources [VDMR] , 2008). After the continued release of these early maps that displayed Virginia Coastal Plain morphology in increasingly greater detail, several geological reports were published where scarps and plains were traced, named, and origins proposed (Shaler, 1890; Shattuck, 1901; Clark and Miller, 1906, 1912; Wentworth, 1930). Shattuck (1901) recognized three continuous scarps in Maryland and most importantly suggested that there was a correlation between surficial deposits (geologic formations) and morphologic terraces; this defined the terrace-formation concept that became the basis of geologic mapping of the Upper Pliocene and Pleistocene sediments in the Virginia Coastal Plain. Clark and Miller (1912) described three Pleistocene formations and supported the terrace-formation concept because the sediments within each formation were lithologically similar to one another, and differentiating them by comparing key vertical sections was not helpful. In Virginia, Wentworth (1930) defined and mapped five Pleistocene terraces on the basis of correlating morphologic terraces with their underlying formation sediments. Similarly, Cooke (1931) identified and traced seven marine shorelines and proposed corresponding terrace names from Maryland to Florida. Cooke (1931) also pointed out that the various elevations of scarps were related to interglacial sea-level highstands whereas prior workers called on vertical land movement of emergence and submergence to create shoreline scarps. Cooke’s (1936) geologic map of eastern South Carolina showing the Pleistocene sediments is also based on the terrace-formation concept. Detailed geologic mapping of the Virginia Coastal Plain by Virginia Division of Mineral Resources (VDMR) and U.S. Geological Survey (USGS) followed lithostratigraphic convention for Early Pliocene and older marine sediments. However, geologic mapping of Upper Pliocene and Pleistocene sediments used the concept of morphologic setting and elevations of scarp toe and flat-plain, or the terrace-formation concept. The 1983 version of the North American Stratigraphic Code (hereafter the Code) (North American Commission on Stratigraphic Nomenclature [NACSN] , 1983) introduced allostratigraphy that gave credence to the terrace-formation concept. According to the current Code (NACSN, 2021; Article 58), “an allostratigraphic unit is a mappable body of rock that is defined and identified on the basis of its bounding discontinuities,” and, “formal allostratigraphic units may be defined or distinguish between different (1) superposed discontinuity-bounded deposits of similar lithology, (2) contiguous discontinuity-bounded deposits of similar lithology, or (3) geographically separated discontinuity-bounded units of similar lithology…The lithology of an allostratigraphic unit plays no part in its definition.” However, definition and placement of bounding discontinuities relies on lateral and vertical knowledge of lithology and is so recognized in Article 58e, “A type locality and type area must be designated; a composite stratotype or a type section and several reference sections are desirable.” We understand that each of the terrace formations of Late Pliocene and Pleistocene age were deposited during a marine transgressive event tied to an interglacial sea-level highstand (scarp above the Poquoson Member of the Tabb Formation may be the result of a still-stand as sea level was falling after deposition of the Lynnhaven Member, refer to fig. 1). The base of each map unit is an unconformity (or is at least in part for the base of the Poquoson) and nomenclature for each unit might be changed from “formation” to “alloformation,” as suggested by Doar and Berquist (2009). Figure 1 shows the relative placement of these mapped and named alloformations as used in the Virginia Coastal Plain and Fall Zone (a zone of transition between the Piedmont and Coastal Plain Provinces) (Berquist, 2021). Except perhaps for the Poquoson, many of these proposed alloformations are unlikely to interfinger or grade laterally seaward into extensive deep marine sediments. After their deposition, erosion during sea-level fall and then rise for the next depositional event would seemingly obliterate any marine connection without a basal unconformity. All these alloformations contain sediments deposited in fluvial, estuarine, and shallow-marine environments and are bound everywhere by unconformities. This paper documents the Cold Harbor Alloformation as a formal stratigraphic unit and the rationale for mapping the sediments using an allostratigraphic classification. The Cold Harbor Alloformation is composed of sediments from a variety of environments such as gravel beds, shoreline beach sands, marsh and estuarine sand and mud, and shallow marine sands. These sediments are found at elevations up to 240 ft and lie unconformably above older sediments and rocks. Our work relied on observations of outcrops and many auger borings to identify basal unconformities.

Virginia↗

Formalization of the Elsing Green Alloformation, Virginia Coastal Plain, Eastern Virginia

The Elsing Green Alloformation is herein formalized as a new Pleistocene alloformation in the Virginia Coastal Plain. This allostratigraphic unit is characterized by an unconformity-bounded terrace deposit occupying a morphologic flat with a mapped elevation range of about 26 to 38 feet (7.9 to 11.6 meters); its landward extent is commonly marked at the base of an erosional scarp where the toe is typically at 38 feet (11.6 meters) and is distinguished by its topographic expression and stratigraphic position. Sediments of this unit were previously mapped as part of the Shirley Formation or Sedgefield Member of the Tabb Formation but differ significantly in their bounding scarps and terrace morphology. The type area is within the Tunstall 7.5-minute quadrangle near Elsing Green Plantation, King William County, Virginia. The reference section is provided by boring 2012270-2. Correlation is proposed with other mid-Pleistocene terrace deposits in the mid-Atlantic Coastal Plain, including the younger Lynch Heights Formation of Delaware and the Butlers Bluff Member of the Nassawadox Formation of the Eastern Shore of Virginia. The age of the Elsing Green Alloformation is interpreted to be Marine Isotope Stage 9e (about 330,000 years).

Virginia↗

Resurvey of the Marble Canyon and Bridge Canyon dam sites in Grand Canyon National Park—Changes in sediment storage and evidence supporting the occurrence of bedrock incision through the mid-20th century

The Bureau of Reclamation developed an extensive plan for a network of dams, water tunnels, and hydropower plants in and around Grand Canyon, Arizona, in the 1940s through 1960s. The two largest of these planned dams were the Marble Canyon and Bridge Canyon Dams on the Colorado River. Though these dams were ultimately never built, Reclamation conducted extensive topographic, bathymetric, and subsurface exploration work at the sites proposed for these dams in the 1940s and 1950s. Resurveys of these dam sites were conducted between 1998 and 2021 to determine the changes in sediment storage at these dam sites caused by the upstream construction and operation of Glen Canyon Dam and by the recession of Lake Mead, the reservoir impounded by Hoover Dam. The resurveys of the Marble Canyon dam sites indicate that the post-1950s changes in sediment storage at these dam sites are broadly consistent with flux-based estimates of voluminous sand erosion from Marble Canyon since the 1963 closure of Glen Canyon Dam. These resurveys also suggest that the pre-dam longitudinal variation in sediment thickness over bedrock played a key role in determining the locations of the sand erosion induced by Glen Canyon Dam; more sand eroded from locations where more sand was present in the 1950s. The resurvey of the Bridge Canyon dam sites indicates that the Colorado River’s incision of the Lake Mead delta is regulated both by bed-sediment grain size and downstream hydraulic controls. Finally, analyses of bed-sediment thickness and sedimentological data at the dam sites, and observations of bed scour and gravel transport, suggest that sufficient bedrock was exposed to allow bedrock incision during commonly recurring pre-dam snowmelt floods that entrained small boulders into transport.

Arizona↗