Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Interpretation”

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.

At least 829 records · Page 46Linked to original sources

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↗

Comparison of simple averaging and latent class modeling to estimate the area of land cover in the presence of reference data variability

Estimates of the area or percent area of the land cover classes within a study region are often based on the reference land cover class labels assigned by analysts interpreting satellite imagery and other ancillary spatial data. Different analysts interpreting the same spatial unit will not always agree on the land cover class label that should be assigned. Two approaches for accommodating interpreter variability when estimating the area are simple averaging (SA) and latent class modeling (LCM). This study compares agreement between area estimates obtained from SA and LCM using reference data obtained by seven trained, professional interpreters who independently interpreted an annual time series of land cover reference class labels for 300 sampled Landsat pixels. We also compare the variability of the LCM and SA area estimates over different numbers of interpreters and different subsets of interpreters within each interpreter group size, and examine area estimates of three land cover classes (forest, developed, and wetland) and three change types (forest gain, forest loss, and developed gain). Differences between the area estimates obtained from SA and LCM are most pronounced for the estimates of wetland and the three change types. The percent area estimates of these rare classes were usually greater for LCM compared to SA, with the differences between LCM and SA increasing as the number of interpreters providing the reference data increased. The LCM area estimates generally had larger standard deviations and greater ranges over different subsets of interpreters, indicating greater sensitivity to the selection of the individual interpreters who carried out the reference class labeling.

Washington↗

Geochronologic and geochemical data from metasedimentary and associated rocks in the Lane Mountain area, San Bernardino County, California

Eugeoclinal metasedimentary and metavolcanic rocks in the Lane Mountain area, California, are considered part of the El Paso terrane, which is commonly thought to have been displaced several hundred kilometers (km) southeastward from its place of origin during late Paleozoic truncation of the North American continental margin. Uranium-lead dating of detrital zircons from this area was undertaken to limit the depositional ages of these nearly non-fossiliferous metamorphic rocks. Analysis of detrital zircons from 17 metasedimentary rock samples yielded a composite age distribution that ranges from Archean to Jurassic and has significant peaks at ~2,800 2,400 mega-annum (Ma), 2,100–1,600 Ma, and ~300–200 Ma. The Proterozoic and Archean ages indicate derivation from continental sources in ancestral North America, whereas the late Paleozoic and Mesozoic ages are interpreted as derived from a magmatic arc that began to develop along the continental margin in Permian to Triassic time. The 17 detrital zircon samples are from quartzitic and conglomeratic rocks of the Carbide, Williams Well, and Noble Well formations, which were informally named by T.H. McCulloh in 1960. The zircon data indicate that the oldest rocks in the Carbide formation are quartzites likely correlative with the Ordovician Eureka Quartzite of the Cordilleran miogeocline. These rocks lie structurally above the rest of the Carbide formation, different units of which yielded zircons that indicate maximum depositional ages ranging from middle Paleozoic to Late Triassic. Zircons from the Williams Well and Noble Well formations indicate maximum depositional ages of late Paleozoic and Early Jurassic, respectively. The Noble Well formation is interpreted to correlate with the lithologically similar, Early Jurassic, Fairview Valley Formation of the Black and Quartzite Mountain areas some 60 km to the southwest. The above interpretations depend on the presumption that the detrital zircons in these samples did not undergo extreme, postdepositional lead loss, which would result in misleadingly young ages. Although such lead loss is considered unlikely for these samples, further work could test the validity of this interpretation. Zircons from six additional samples were also analyzed: (1) a quartzite from which all the zircons are interpreted to have formed by Late Jurassic metamorphism; (2) three samples interpreted as albitized igneous rocks of Middle Permian age; and (3) two samples interpreted as fine-grained monzonite to diorite of Late Jurassic age. Both sets of igneous rocks were initially thought to be metasedimentary but were reinterpreted as igneous largely on the basis of the zircon data. Based on the interpretations presented here, this study demonstrates that the depositional, magmatic, and deformational history of the El Paso terrane was longer and more complex than previously thought and will require reevaluation of existing tectonic models involving this terrane.

California↗

Integration of geological remote-sensing techniques in subsurface analysis

Geological remote sensing is defined as the study of the Earth utilizing electromagnetic radiation which is either reflected or emitted from its surface in wavelengths ranging from 0.3 micrometre to 3 metres. The natural surface of the Earth is composed of a diversified combination of surface cover types, and geologists must understand the characteristics of surface cover types to successfully evaluate remotely-sensed data. In some areas landscape surface cover changes throughout the year, and analysis of imagery acquired at different times of year can yield additional geological information. Integration of different scales of analysis allows landscape features to be effectively interpreted. Interpretation of the static elements displayed on imagery is referred to as an image interpretation. Image interpretation is dependent upon: (1) the geologist's understanding of the fundamental aspects of image formation, and (2.) his ability to detect, delineate, and classify image radiometric data; recognize radiometric patterns; and identify landscape surface characteristics as expressed on imagery. A geologic interpretation integrates surface characteristics of the landscape with subsurface geologic relationships. Development of a geologic interpretation from imagery is dependent upon: (1) the geologist's ability to interpret geomorphic processes from their static surface expression as landscape characteristics on imagery, (2) his ability to conceptualize the dynamic processes responsible for the evolution 6f interpreted geologic relationships (his ability to develop geologic models). The integration of geologic remote-sensing techniques in subsurface analysis is illustrated by development of an exploration model for ground water in the Tucson area of Arizona, and by the development of an exploration model for mineralization in southwest Idaho.

Open-File Report↗

Sub-crop geologic map of pre-Tertiary rocks in the Yucca Flat and northern Frenchman Flat areas, Nevada Test Site, southern Nevada

This map displays interpreted structural and stratigraphic relations among the Paleozoic and older rocks of the Nevada Test Site region beneath the Miocene volcanic rocks and younger alluvium in the Yucca Flat and northern Frenchman Flat basins. These interpretations are based on a comprehensive examination and review of data for more than 77 drillholes that penetrated part of the pre-Tertiary basement beneath these post-middle Miocene structural basins. Biostratigraphic data from conodont fossils were newly obtained for 31 of these holes, and a thorough review of all prior microfossil paleontologic data is incorporated in the analysis. Subsurface relationships are interpreted in light of a revised regional geologic framework synthesized from detailed geologic mapping in the ranges surrounding Yucca Flat, from comprehensive stratigraphic studies in the region, and from additional detailed field studies on and around the Nevada Test Site. All available data indicate the subsurface geology of Yucca Flat is considerably more complicated than previous interpretations have suggested. The western part of the basin, in particular, is underlain by relics of the eastward-vergent Belted Range thrust system that are folded back toward the west and thrust by local, west-vergent contractional structures of the CP thrust system. Field evidence from the ranges surrounding the north end of Yucca Flat indicate that two significant strike-slip faults track southward beneath the post-middle Miocene basin fill, but their subsurface traces cannot be closely defined from the available evidence. In contrast, the eastern part of the Yucca Flat basin is interpreted to be underlain by a fairly simple north-trending, broad syncline in the pre-Tertiary units. Far fewer data are available for the northern Frenchman Flat basin, but regional analysis indicates the pre- Tertiary structure there should also be relatively simple and not affected by thrusting. This new interpretation has implications for ground water flow through pre-Tertiary rocks beneath the Yucca Flat and northern Frenchman Flat areas, and has consequences for ground water modeling and model validation. Our data indicate that the Mississippian Chainman Shale is not a laterally extensive confining unit in the western part of the basin because it is folded back onto itself by the convergent structures of the Belted Range and CP thrust systems. Early and Middle Paleozoic limestone and dolomite are present beneath most of both basins and, regardless of structural complications, are interpreted to form a laterally continuous and extensive carbonate aquifer. Structural culmination that marks the French Peak accommodation zone along the topographic divide between the two basins provides a lateral pathway through highly fractured rock between the volcanic aquifers of Yucca Flat and the regional carbonate aquifer. This pathway may accelerate the migration of ground-water contaminants introduced by underground nuclear testing toward discharge areas beyond the Nevada Test Site boundaries. Predictive three-dimensional models of hydrostratigraphic units and ground-water flow in the pre-Tertiary rocks of subsurface Yucca Flat are likely to be unrealistic due to the extreme structural complexities. The interpretation of hydrologic and geochemical data obtained from monitoring wells will be difficult to extrapolate through the flow system until more is known about the continuity of hydrostratigraphic units.

Nevada↗

How good are estimates of transmissivity from slug tests in fractured rock?

Slug tests in fractured rock usually are interpreted with models that assume homogeneous formation properties, even though hydraulic properties of fractures can vary by many orders of magnitude over the length of boreholes. To investigate the impact of heterogeneity on the interpretation of slug tests in fractured rock, slug tests were conducted over large intervals of boreholes in crystalline rock in central New Hampshire, and interpreted using a homogeneous model. The results of the slug tests were then compared with estimates of transmissivity from fluid-injection tests conducted over shorter intervals in the same boreholes. The fluid-injection tests showed transmissivity to vary more than six orders of magnitude over the length of the boreholes; however, the sum of the transmissivities from the fluid-injection tests were within an order of magnitude of the transmissivity estimated from the slug tests. Although the two estimates of transmissivity were within an order of magnitude of each other, the water level responses during the slug tests did not exactly match the responses predicted by the homogeneous model. To investigate the effect of heterogeneity on water level responses during slug tests, a Laplace-transform solution was developed for slug tests conducted in boreholes containing multiple fractures with hydraulic properties that vary over the length of the borehole. A comparison of this solution with the homogeneous model shows no difference between the shape of water level responses in a homogeneous formation and a (layered) heterogeneous formation. Furthermore, the transmissivity estimated using a homogeneous model is within an order of magnitude of the prescribed transmissivity in the heterogeneous model. Thus, differences between responses predicted from a homogeneous model and measured water levels during slug tests can be attributed to phenomena such as nonradial flow in the vicinity of the borehole, and not heterogeneous hydraulic properties over the length of the borehole. The experimental results of this investigation show that even when conditions such as nonradial flow are present in the vicinity of the borehole, interpretations of slug tests using a homogeneous model provided order-of-magnitude estimates of transmissivity in the crystalline rock terrane under consideration.Slug tests in fractured rock are usually interpreted with models that assume homogeneous formation properties, although hydraulic properties of fractures can vary by many order of magnitude over the length of the boreholes. To investigate the impact of heterogeneity on the interpretation of slug tests in fractured rocks, slug tests were conducted over large intervals of boreholes in crystalline rock in central New Hampshire. Tests results were interpreted using a homogeneous model. A Laplace-transform solution was developed for slug tests conducted in boreholes containing multiple fractures with hydraulic properties. Finally, a comparison was made between this solution and the homogeneous model.

Ground Water↗

Borehole geophysical investigation of a formerly used defense site, Machiasport, Maine, 2003-2006

The U.S. Geological Survey, in cooperation with the U.S. Army Corps of Engineers, collected borehole geophysical logs in 18 boreholes and interpreted the data along with logs from 19 additional boreholes as part of an ongoing, collaborative investigation at three environmental restoration sites in Machiasport, Maine. These sites, located on hilltops overlooking the seacoast, formerly were used for military defense. At each of the sites, chlorinated solvents, used as part of defense-site operations, have contaminated the fractured-rock aquifer. Borehole geophysical techniques and hydraulic methods were used to characterize bedrock lithology, fractures, and hydraulic properties. In addition, each geophysical method was evaluated for effectiveness for site characterization and for potential application for further aquifer characterization and (or) evaluation of remediation efforts. Results of borehole geophysical logging indicate the subsurface is highly fractured, metavolcanic, intrusive, metasedimentary bedrock. Selected geophysical logs were cross-plotted to assess correlations between rock properties. These plots included combinations of gamma, acoustic reflectivity, electromagnetic induction conductivity, normal resistivity, and single-point resistance. The combined use of acoustic televiewer (ATV) imaging and natural gamma logs proved to be effective for delineating rock types. Each of the rock units in the study area could be mapped in the boreholes, on the basis of the gamma and ATV reflectivity signatures. The gamma and mean ATV reflectivity data were used along with the other geophysical logs for an integrated interpretation, yielding a determination of quartz monzonite, rhyolite, metasedimentary units, or diabase/gabbro rock types. The interpretation of rock types on the basis of the geophysical logs compared well to drilling logs and geologic mapping. These results may be helpful for refining the geologic framework at depth. A stereoplot of all fractures intersecting the boreholes indicates numerous fractures, a high proportion of steeply dipping fractures, and considerable variation in fracture orientation. Low-dip-angle fractures associated with unloading and exfoliation are also present, especially at a depth of less than 100 feet below the top of casing. These sub-horizontal fractures help to connect the steeply dipping fractures, making this a highly connected fracture network. The high variability in the fracture orientations also increases the connectivity of the fracture network. A preliminary comparison of all fracture data from all the boreholes suggests fracturing decreases with depth. Because all the boreholes were not drilled to the same depth, however, there is a clear sampling bias. Hence, the deepest boreholes are analyzed separately for fracture density. For the deepest boreholes in the study, the intensity of fracturing does not decline significantly with depth. It is possible the fractures observed in these boreholes become progressively tighter or closed with depth, but this is difficult to verify with the borehole methods used in this investigation. The fact that there are more sealed fractures at depth (observed in optical televiewer logs in some of the boreholes) may indicate less opening of the sealed fractures, less water moving through the rock, and less weathering of the fracture infilling minerals. Although the fracture orientation remained fairly constant with depth, differences in the fracture patterns for the three restoration sites indicate the orientation of fractures varies across the study area. The fractures in boreholes on Miller Mountain predominantly strike northwest-southeast, and to a lesser degree they strike northeast. The fractures on or near the summit of Howard Mountain strike predominantly east-west and dip north and south, and the fractures near the Transmitter Site strike northeast-southwest and dip northwest and southeast. The fracture populations for the boreholes on or near the summit of Howard Mountain show more variation than at the other two sites. This variation may be related to the proximity of the fault, which is northeast of the summit of Howard Mountain. In a side-by-side comparison of stereoplots from selected boreholes, there was no clear correspondence between fracture orientation and proximity to the fault. There is, however, a difference in the total populations of fractures for the boreholes on or near the summit of Howard Mountain and the boreholes near the Transmitter Site. Further to the southwest and further away from the fault, the fractures at the Transmitter Site predominantly strike northeast-southwest and northwest-southeast.Heat-pulse flowmeter (HPFM) logging was used to identify transmissive fractures and to estimate the hydraulic properties along the boreholes. Ambient downflow was measured in 13 boreholes and ambient upflow was measured in 9 boreholes. In nine other bedrock boreholes, the HPFM did not detect measurable vertical flow. The observed direction of vertical flow in the boreholes generally was consistent with the conceptual flow model of downward movement in recharge locations and upward flow in discharge locations or at breaks in the slope of land surface. Under low-rate pumping or injection rates [0.25 to 1 gallon per minute (gal/min)], one to three inflow zones were identified in each borehole. Two limitations of HPFM methods are (1) the HPFM can only identify zones within 1.5 to 2 orders of magnitude of the most transmissive zone in each borehole, and (2) the HPFM cannot detect flow rates less than 0.010 + or - 0.005 gal/min, which corresponds to a transmissivity of about 1 foot squared per day (ft2/d). Consequently, the HPFM is considered an effective tool for identifying the most transmissive fractures in a borehole, down to its detection level. Transmissivities below that cut-off must be measured with another method, such as packer testing or fluid-replacement logging. Where sufficient water-level and flowmeter data were available, HPFM results were numerically modeled. For each borehole model, the fracture location and measured flow rates were specified, and the head and transmissivity of each fracture zone were adjusted until a model fit was achieved with the interpreted ambient and stressed flow profiles. The transmissivities calculated by this method are similar to the results of an open-hole slug test; with the added information from the flowmeter, however, the head and transmissivity of discrete zones also can be determined. The discrete-interval transmissivities ranged from 0.16 to 330 ft2/d. The flowmeter-derived open-hole transmissivity, which is the combined total of each of the transmissive zones, ranged from 1 to 511 ft2/d. The whole-well open-hole transmissivity values determined with HPFM methods were compared to the results of open-hole hydraulic tests. Despite the fact that the flowmeter-derived transmissivities consistently were lower than the estimates derived from open-hole hydraulic tests alone, the correlation was very strong (with a coefficient of determination, R2, of 0.9866), indicating the HPFM method provides a reasonable estimate of transmissivities for the most transmissive fractures in the borehole. Geologic framework, fracture characterization, and estimates of hydraulic properties were interpreted together to characterize the fracture network. The data and interpretation presented in this report should provide information useful for site investigators as the conceptual site groundwater flow model is refined. Collectively, the results and the conceptual site model are important for evaluating remediation options and planning or implementing the design of a well field and borehole completions that will be adequate for monitoring flow, remediation efforts, groundwater levels, and (or) water quality. Similar kinds of borehole geophysical logging (specifically the borehole imaging, gamma, fluid logs, and HPFM) should be conducted in any newly installed boreholes and integrated with interpretations of any nearby boreholes. If boreholes are installed close to existing or other new boreholes, cross-hole flowmeter surveys may be appropriate and may help characterize the aquifer properties and connections between the boreholes.

Maine↗

Guidebook to the Gaudalupian symposium

Compared to the Guadalupe Mountains of Texas and New Mexico the depositional environments of the Permian strata of the Glass Mountains (and adjacent Del Norte Mountains) are less well known. In general, the Guadalupian facies in the the Glass and Del Norte mountains changes from predominantly carbonate facies in the northeast to thicker clastic facies in the southwest. Philip B. Kind (1931) originally considered this trend to reflect an uplifted clastic source to the southwest, with carbonate facies developing away from the source area. Ross (1986) interpreted the eastern portion of the Road Canyon and Word formations to consist the shelf, shelf-edge bioherm, and reef facies, and the southwest area to consist of deeper water siliceous shale, clastic limestone, and basinal sandstone facies. Probably the best known controversy in the Glass Mountains involves the depositional environment of the Skinner Ranch Formation (Leonardian according to Ross, 1986; Wolfcampian according to Cooper and Grant, 1972) at its type section on Leonard Mountain. Cooper and Grant (1964) identified in situ patch reefs at the base of the section, which were subsequently interpreted as displaced limestone blocks deposited in a slope environment (Rogers, 1972; Cys and Mazzullo, 1978; Ross, 1986). Later Flores, McMillan, and Watters (1977) interpreted the same units as subtidal and intertidal deposits. The Skinner Ranch Formation illustrates the complexities involved in interpreting the paleogeography of the Glass Mountains. If the Sinner Ranch contains displaced blocks, some eroded from older units, it explains the occurrence of Wolfcampian fossils in the Skinner Ranch (Ross, 1986).The slop facies interpretation also is used to place the shelf edge at that time between Skinner Ranch outcrops at Leonard Mountain and the lagoonal, backreef deposits of the Hess Formation to the east, although most of the actual shelf edge is not preserved (Ross, 1987:30). Similar conflicting interpretations exist in younger rocks in the western facies of the Leonardian Guadalupian to the southwest in the Del Norte Mountains. Ross (1986, 1987) considered the western facies of the Road Canyon and Word formations to be basinal shales and turbidites. Wardlaw et al. (1990) and Rohr et al. (1987) have interpreted this area to be shallow intertidal to lagoonal environments adjacent to an uplifted area to the south. The type section of the Road Canyon Formation is also a subject of disagreement and will be discusses in more detail later.

Texas↗

Understanding the resource potential of natural hydrogen on Earth: Scientific gaps, uncertainties and recommendations

A comprehensive scientific research roadmap is essential to bridge knowledge gaps and deepen the understanding of key geological, geochemical, and geophysical aspects of natural hydrogen (H 2 ) as a potential new energy resource. This paper reviews major scientific uncertainties on natural H 2 , suggesting research priorities, as a guide for defining exploration strategies, techniques, and data interpretation. The uncertainties concern all phases of the natural H 2 cycle, from generation (source rocks) through migration (advection and diffusion) and accumulation (reservoir and cap rocks) to the application and interpretation of subsurface and surface geochemical and geophysical exploration techniques. Understanding H 2 sources and generation rates (the amount of H 2 generated by a given volume of rock over time) is crucial for determining whether a geological H 2 system operates as a short-term dynamic system with rapid H 2 production and release, or as a conventional gas system with long-term accumulations, analogous to petroleum reservoirs. Preliminary estimates for serpentinisation, radiolysis, and organic matter degradation suggest that H 2 generation is not inherently fast, especially for non-hydrothermal continental systems (crystalline basement of shields, ophiolites, peridotite massifs, sedimentary basins), and long-term accumulations, like those of fossil natural gas systems, represent the most likely scenario. The mechanisms of H 2 migration through geological formations require application of fundamental principles of fluid-flow physics, distinguishing advection and diffusion, as well as their forms (from gas-phase, bubble flows to aqueous solutions). Additional studies of H 2 accumulation and retention in subsurface reservoirs could improve understanding of mechanisms of H 2 migration by focusing on the rock fluid-bearing properties and the factors affecting H 2 preservation, such as the presence of cap rocks impermeable to H 2 , pressure conditions, residence times, and microbial or abiotic consumption. Advanced techniques, including reservoir modelling, flow simulations, 3D imaging (micro-CT) of H 2 -bearing rocks, and extraction and analysis of gas occluded in rocks, can provide insights into the stability and potential recoverability of H 2 accumulations. The interpretation of surface exploration techniques, including gas geochemistry, geophysics, and remote sensing, long employed in mineral and energy resource exploration, is now being adapted for natural H 2 studies, but challenges remain in the data interpretation. Distinguishing H 2 seepage due to geological degassing from H 2 produced near the surface by modern microbial processes or artificial sources, such as hammering or drilling for soil-gas sampling, drilling into aquifers, and corrosion in boreholes, is an essential step in exploration. The simple detection of H 2 in soils, even in morphological structures like sub-circular depressions or “fairy circles”, cannot be cursorily interpreted as a signal of natural H 2 seepage from a deep source. A holistic geochemical approach, including isotopic analyses of gases associated with H 2 , is recommended to distinguish among the variety of possible H 2 origins. Observations of H 2 in wells should be interrogated to rule out possible artifacts such as corrosion and drill bit metamorphism. The integration of multiple geophysical methods, including seismic, gravimetric, magnetic, and electro-magnetic surveys, is recommended to mitigate interpretation ambiguities regarding the structure of a subsurface H 2 system (source and reservoir rocks, including fluid and gas storage), due to the non-uniqueness of rock-specific physical properties.

Earth-Science Reviews↗

An exploration of open-system uranium-series history of marine terrace corals, Perachora Peninsula, Greece and San Nicolas Island, California, USA

Emergent marine terraces record past periods of relatively high sea level and are common on uplifting coasts worldwide. Such landforms are extensive around the Gulf of Corinth area of Greece and along the coast of California, USA. In the Gulf of Corinth region, marine terraces record Quaternary uplift due to ongoing basin extension on the northern side of the Africa-Eurasia-Aegean Sea plate boundary. Although dating of terrace fossils has been extensive in the Perachora Peninsula part of this region, some apparent ages are much older (∼190 ka to ∼150 ka) and some are much younger (∼91 ka to ∼69 ka) than expected for a marine terrace that has long been hypothesized to be of last-interglacial (∼120 ka) age. New uranium-series analyses of terrace corals from two previously studied localities are reported here. As with previous work, the new results show that these corals have behaved as open systems, with back-calculated, initial 234 U/ 238 U values higher than those found in modern seawater. New and previously published U-series analyses of last-interglacial corals from San Nicolas Island, California (USA) illustrate the possible mechanisms of open-system behavior, which include both “early” alpha-recoil-derived additions of 230 Th and 234 U, and “continuing” alpha-recoil-derived additions of these nuclides. Many corals from the Perachora Peninsula appear to have experienced the “continuing” alpha-recoil-derived open-system history, which explains the apparent older ages, and permits an interpretation that all analyzed corals likely date to the last interglacial period, ∼120 ka. If the estimated age of ∼120 ka reported here is correct, this age, along with considerations of glacial isostatic adjustment (GIA) modeling of past sea level, indicates that the southern side of the Perachora Peninsula has experienced uplift rates ranging from ∼0.23 m/kyr to ∼0.57 m/kyr, increasing to the southeast, in agreement with previous work. Interpretation of the geomorphology of the northern side of the peninsula has been controversial for more than four decades. The “stepped” topography has been interpreted as either a series of successively higher and older marine terraces, as on the southern side, or a marine terrace of one age that has been displaced by a series of shore-parallel, northward-dipping normal faults. In the second interpretation, topographic steps in the landscape are interpreted to be fault scarps, not ancient sea cliffs. Interpretation of the open-system U-series ages indicates that marine deposits at elevations of ∼22 m and ∼44 m are both likely ∼120 ka, supporting the scenario that steps in the landscape are likely fault scarps, not sea cliffs. Using the new U-series ages of ∼120 ka, along with GIA-modeled paleo-sea level confirms another previous study that uplift rates on the northern side of the peninsula are at least as high as 0.39–0.46 m/kyr.

California↗

CASERTZ aeromagnetic data reveal late Cenozoic flood basalts(?) in the West Antarctic rift system

The late Cenozoic volcanic and tectonic activity of the enigmatic West Antarctic rift system, the least understood of the great active continental rifts, has been suggested to be plume driven. In 1991-1992, as part of the CASERTZ (Corridor Aerogeophysics of the Southeast Ross Transect Zone) program, an ∼25000 km aeromagnetic survey over the ice-covered Byrd subglacial basin shows magnetic "texture" critical to interpretations of the underlying extended volcanic terrane. The aeromagnetic data reveal numerous semicircular anomalies ∼100-1100 nT in amplitude, interpreted as having volcanic sources at the base of the ice sheet; they are concentrated along north-trending magnetic lineations interpreted as rift fabric. Models constrained by coincident radar ice soundings indicate highly magnetic sources, with a probable high remanent magnetization in the present field direction, strongly suggesting a late Cenozoic age. Magnetic anomalies over exposed late Cenozoic volcanic rocks along part of the rift shoulder and in coastal Marie Byrd Land are similar in form and amplitude. The CASERTZ aeromagnetic results, combined with >100 000 km of widely spaced aeromagnetic profiles, indicate at least 10 6 km 3 of probable late Cenozoic volcanic rock (flood basalt?) in the West Antarctic rift beneath the ice sheet and Ross Ice Shelf. Comparison with other plumes in active rift areas (e.g., Yellowstone and East Africa) indicates that this volume estimate lies in the range of magma generation found in these other low-extension continental rifts.

Geology↗

Minturn and Sangre de Cristo formations of Southern Colorado: A prograding fan delta and alluvial fan sequence shed from the ancestral Rocky Mountains

The Middle Pennsylvania!) Minturn Formation and the Pennsylvanian-Permian Sangre de Cristo Formation of the northern Sangre de Cristo Range form a thick progradational sequence of coarse clastic sediments. These sediments were deposited along the western margin of the central Colorado trough during uplift of the late Paleozoic Uncompahgre highland, a major structural and topographic feature of the Ancestral Rocky Mountains. The Minturn is composed mostly of sandstone and shale deposited by fan deltas that prograded into the central Colorado trough. The Minturn of the northern Sangre de Cristo Range is divisible into a turbidite-bearing facies, a limestone-bearing facies, and a redbed facies. The turbidite-bearing facies is interpreted as deposits of fan deltas that prograded onto the sea bottom below wave base. The limestone-bearing facies is interpreted as deposits of fan deltas that prograded onto a shallow sea bottom above wave base. In this facies, sandstones containing deltaic foresets overlie thin shallow marine limestones and are considered diagnostic of that facies. Both facies contain thick intervals of sandstone and shale interpreted as deltaic and alluvial deposits. Where it onlaps the Uncompahgre highland, the lower part of the Minturn Formation contains quartzose and arkosic redbeds of probable alluvial origin. The continental Sangre de Cristo Formation conformably overlies the Minturn Formation basinward, but it unconformably overlies the Minturn and Precambrian basement near the Uncompahgre highland. The Sangre de Cristo contains two facies. A sandstone facies consists of fining-upward cycles of red conglomeratic sandstone and siltstone interpreted as braided stream deposits on distal parts of alluvial fans. A conglomerate facies (Crestone Conglomerate Member) consists of poorly sorted conglomerates interpreted as debris flow and mudflow deposits and sorted conglomerate and sandstone interpreted as streamflow and sheetflow deposits on proximal parts of alluvial fans. The fan deltas and alluvial fans that deposited the Minturn and Sangre de Cristo formations probably developed in response to faulting and uplift of the Uncompahgre highland. The highland rose at least three times from Middle Pennsylvanian to Early Permian time. Later, during the Laramide orogeny, the faulted boundary between the highland and the central Colorado trough was destroyed by thrusting. The boundary faults were not strike-slip, as shown by the distribution of clasts in conglomerates directly downstream from identical Precambrian rocks that represent part of the highland.

Colorado↗

Comparing landslide inventories: The map depends on the method

Landslide inventory maps are generally prepared by interpreting the geomorphic expression of landsliding on aerial photos, topographic maps, or on the ground. Distinctive landslide geomorphology allows the recognition and mapping of landslides, although there are always landslides that have very subtle expression and are not identified. The difficulties of mapping landslides based on their geomorphic expression are amplified in heavily forested terrain. The ground surface is obscured by tree cover on aerial photographs, and landslide-related features are often hidden. This limitation affects not only aerial photo interpretation, but also interpretation of topographic maps, which are based on aerial photographs. We compared five maps showing landslides in the Laurel Quadrangle in the Santa Cruz Mountains, California. These include a geologic map, a map prepared for the county based on interpretation of aerial photographs, a map prepared by us based on aerial photographs and compilation of previous work, a map of features interpreted from the U.S. Geological Survey 7.5-minute topographic map, and a detailed field-based landslide map. Comparison of these maps shows that the geologic map identifies few landslides, but most landslides on the geologic map are also shown on the other maps. The two maps based mainly on aerial photo interpretation tend to show the larger slides, but there is only about 60 percent correspondence of landslide areas between the two. Comparing the reconnaissance techniques with the much more detailed field mapping shows that the reconnaisance maps emphasize the large slides of bedrock and identify a lower percentage of shallow debris slides and debris flows.

Environmental & Engineering Geoscience↗

Base of principal aquifer for parts of the North Platte, South Platte, and Twin Platte Natural Resources Districts, western Nebraska

Water resources in the North and South Platte River valleys of Nebraska, including the valley of Lodgepole Creek, are critical to the social and economic health of the area, and for the recovery of threatened and endangered species in the Platte River Basin. Groundwater and surface water are heavily used resources, and uses are regulated in the study area. Irrigation is the dominant water use and, in most instances, is supplied by both groundwater and surface-water sources. The U.S. Geological Survey and its partners have collaborated to use airborne geophysical surveys for areas of the North and South Platte River valleys including the valley of Lodgepole Creek in western Nebraska. The objective of the surveys was to map the aquifers and underlying bedrock topography of selected areas to help improve the understanding of groundwater–surface-water relations to guide water-management decisions. This project was a cooperative study involving the North Platte Natural Resources District, the South Platte Natural Resources District, the Twin Platte Natural Resources District, the Conservation and Survey Division of the University of Nebraska-Lincoln, and the Nebraska Environmental Trust. This report presents the interpreted base-of-aquifer surface for part of the area consisting of the North Platte Natural Resources District, the South Platte Natural Resources District, and the Twin Platte Natural Resources District. The interpretations presented herein build on work done by previous researchers from 2008 to 2009 by incorporating additional airborne electromagnetic survey data collected in 2010 and additional test holes from separate, related studies. To make the airborne electromagnetic data useful, numerical inversion was used to convert the measured data into a depth-dependent subsurface resistivity model. An interpretation of the elevation and configuration of the base of aquifer was completed in a geographic information system that provided x, y, and z coordinates. The process of interpretation involved manually picking locations (base-of-aquifer elevations) on the displayed airborne electromagnetic-derived resistivity profile by the project geophysicist, hydrologist, and geologist. These locations, or picks, of the base-of-aquifer elevation (typically the top of the Brule Formation of the White River Group) were then stored in a georeferenced database. The pick was made by comparing the inverted airborne electromagnetic-derived resistivity profile to the lithologic descriptions and borehole geophysical logs from nearby test holes. The database of interpretive picks of the base-of-aquifer elevation was used to create primary input for interpolating the new base-of-aquifer contours. The automatically generated contours were manually adjusted based on the interpreted location of paleovalleys eroded into the base-of-aquifer surface and associated bedrock highs, many of which were unmapped before this study. When contours are overlain by the water-table surface, the saturated thickness of the aquifer can be computed, which allows an estimate of total water in storage. The contours of the base-of-aquifer surface presented in this report may be used as the lower boundary layer in existing and future groundwater-flow models. The integration of geophysical data into the contouring process facilitated a more continuous and spatially comprehensive view of the hydrogeologic framework.

Nebraska↗

Revisions to the Virginia Coastal Plain hydrogeologic framework southwest of the James River

New drilling information reveals that altitudes of some hydrogeologic units of the Virginia Coastal Plain aquifer system differ by as much as 50 feet (ft) from those previously known, namely the Aquia and Potomac aquifers, the Potomac confining zone, and the Nanjemoy-Marlboro and Saint Marys confining units. In addition, the lateral margins of some hydrogeologic units are located as much as several miles from previously estimated locations. The largest revisions to unit margins were for the Aquia aquifer and the Nanjemoy-Marlboro and Saint Marys confining units. Interpretation of new geophysical logs, sediment core, and cuttings as well as revised interpretations to existing data indicate channels and embayments are also preserved on eroded top surfaces of the shallowest hydrogeologic units including the Yorktown confining zone, Yorktown-Eastover aquifer, Saint Marys confining unit, Potomac confining zone, and Potomac aquifer. Enhanced details on the configuration of part of the aquifer system southwest of the James River are provided by sediment cores and cuttings as well as geophysical logs from 36 recently drilled boreholes. These, along with reinterpretation of data from 93 preexisting boreholes, form the basis for revised top-surface altitudes and margins of hydrogeologic units beneath parts of Prince George, Surry, Sussex, Isle of Wight, and Southampton Counties and the cities of Franklin and Suffolk. Groundwater withdrawals in the Virginia Coastal Plain cause widespread water-level declines, create the potential for saltwater intrusion, and contribute to regionwide land subsidence. A description of the aquifer system, termed a hydrogeologic framework, was developed by the U.S. Geological Survey in 2006 and provides information needed to base withdrawal-permitting decisions by the Virginia Department of Environmental Quality. This revision of part of the hydrogeologic framework southwest of the James River is based on interpretations of both new and previously analyzed borehole data. The revision is strictly confined to the study area extent and hydrogeologic units not found within the study area were not revised and are not discussed in this report. The newly determined hydrogeologic-unit altitudes and margins have implications for groundwater-withdrawal permitting. New interpretations have found that the Yorktown Eastover aquifer is absent in the southwestern part of the City of Suffolk, owing to what is most likely an isolated area of sediment-texture facies change. Most notably, the top-surface altitudes of the Aquia and Potomac aquifers have been lowered by as much as 50 ft from previous interpretations. This means that wells previously believed to be screened in the top of the Potomac aquifer could, based on these new interpretations, be screened in the bottom of the Aquia aquifer. These changes to aquifers in which wells are screened means that there is potentially more room in the groundwater withdrawal permitting for the Potomac aquifer, the largest and most productive aquifer in Virginia, and overpumping occurring in the Aquia aquifer.

Virginia↗

Surface-geophysical investigation of the University of Connecticut landfill, Storrs, Connecticut

A surface-geophysical investigation of the former landfill area at the University of Connecticut in Storrs, Conn. was conducted as part of a preliminary hydrogeologic assessment of the contamination of soil, surface water, and ground water at the site. Geophysical data were used to help determine the dominant direction of fracture strike; subsurface structure of the landfill; locations of possible leachate plumes, fracture zones or conductive lithologic layers; and the location and number of chemical waste-disposal pits. Azimuthal square-array direct-current (dc) resistivity, two-dimensional (2D) dc-resistivity, inductive terrain conductivity, and ground-penetrating radar (GPR) were the methods used to characterize the landfill area.The dominant strike direction of bedrock fractures interpreted from azimuthal square-array resistivity data is north, ranging from 285 to 30 degrees east of True North. These results complement local geologic maps that identify bedrock foliation and fractures that strike approximately north-south and dip 30 to 40 degrees west.The subsurface structure of the landfill was imaged with 2D dc-resistivity profiling data, which were used to interpret a landfill thickness of 10 to 15 meters. Orientation of the landfill trash disposal trenches were detected by azimuthal square-array resistivity soundings; the dimension and the orientation of the trenches were verified by aerial photographs.Inductive terrain conductivity and 2D dc-resistivity profiling detected conductive anomalies that were interpreted as possible leachate plumes near two surface-water discharge areas. The conductive anomaly to the north of the landfill is interpreted to be a shallow leachate plume and dissipates to almost background levels 45 meters north of the landfill. The anomaly to the southwest is interpreted to extend vertically through the overburden and into the shallow bedrock and laterally along the intermittent drainage to Eagleville Brook, terminating 140 meters south of the landfill. Inductive terrain conductivity and 2D dc-resistivity profiling also detected two dipping, sheet-like conductive features that extend vertically into the bedrock. These features were interpreted either as fracture zones filled with conductive fluids or conductive lithologic layers between more resistive layers. One dipping conductive feature was detected south of the landfill, and the other feature was detected to the west of the former chemical waste-disposal pits. Both anomalies strike approximately north-south and dip about 30 degrees to the west.GPR was used unsuccessfully to locate the former chemical waste-disposal pits. Although the entire overburden and the upper few meters of bedrock were imaged, no anomalous features were detected with GPR that could be correlated with the pits. It is possible that the area surveyed by GPR was entirely backfilled after the soil was removed from the site and that the outline of the former chemical waste-disposal pits no longer exists.

Connecticut↗

Use of low-altitude aerial photography to identify submersed aquatic macrophytes

The feasibility of using low-altitude aerial photography to identify beds of submersed macrophytes is demonstrated. True color aerial photos and collateral ground survey information for submersed aquatic macrophyte beds at 10 sites in the St.Clair-Detroit River system were obtained in September 1978. Using the photos and collateral ground survey information, a dichotomous key was developed for the identification of six classes - beds of five genera of macrophytes and one substrate type. A test was prepared to determine how accurately photo interpreters could identify the six classes. The test required an interpreter to examine an unlabeled, outlined area on photographs and identify it using the key. Six interpreters were tested. One pair of interpreters was trained in the interpretation of a variety of aerial photos, a second pair had field experience in the collection and identification of submersed macrophytes in the river system, and a third pair had neither training in the interpretation of aerial photos nor field experience. The criteria that we developed were applied equally well by the interpretors, regardless of their training or experience. Overall accuracy (i.e., omission errors) of all six classes combined was 68% correct, whereas, overall accuracy of individual classes ranged from 50 to 100% correct. Mapping accuracy (i.e. omission and commission errors) of individual classes ranged from 36 to 75%. Although the key developed for this study has only limited application outside the context of the data and sites examined in this study, it is concluded that low-altitude aerial photography, together with limited amounts of collateral ground survey information, can be used to economically identify beds of submersed macrophytes in the St. Clair-Detroit River system and other similar water bodies.

Conference Paper↗