USGS ScienceSearch

SEARCH · USGS Science

Results for “Nature Based Solutions”

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

Where to restore and conserve? A regional benefit cost analysis of coral reef protection and restoration for coastal flood resilience

Momentum is growing for the management of coral reefs as a strategy to reduce climate risks in tropical coastlines. Yet, quantification of the life-time costs, impacts, and benefits remains limited. This study provides one of the first rigorous, spatially explicit, regional-scale Benefit:Cost Analyses (BCA) for coral reef restoration and was designed to meet the BCA requirements of major hazard mitigation funding and programs. This study simulates coastal flooding using a hydrodynamic model under different scenarios representing current coral reef conditions, reef degradation, and reef restoration. These coastal flood maps are used to estimate socioeconomic damages, which are included in a BCA to assess cost-effectiveness and priority areas for coral reef conservation and restoration. The United States Virgin Islands is used as a case study given recent impacts from storms and their new policy that declares reefs as natural infrastructure. The results show that flood risk across the islands of Saint Croix, Saint John, and Saint Thomas is $51.4 million USD per year. Annually, coral reefs prevent flood damages to 481 people and $43.6 million USD of infrastructure, which represents 87 % of the flood risk. These results identify the communities that could receive the greatest benefits from coral conservation and restoration, which overlap significantly. Coral reef restoration is found to be a cost-effective strategy for flood mitigation although how reefs are restored matters. Hybrid coral reef restoration provides greater flood mitigation benefits than ecological restoration, but its cost-effectiveness is overall lower given its greater costs. Many coastal areas have benefit:cost ratios exceeding 1.5, which complies with typical government agencies’ requirements for accessing hazard mitigation funds. This valuation framework helps provide rigorous regional-scale quantification of nature-based coastal protection solutions for coastal risk management decisions.

Journal of Environmental Management

Integrating refugia modeling and decision science to support climate-adaptive conservation in the Sierra Nevada

To reduce detrimental impacts of anthropogenic change, natural resource managers often look for place-based solutions to minimize biodiversity loss. Climate-change refugia, areas buffered from contemporary climate change, can enable the persistence of valued natural resources and prolong the benefits of conservation action. Here we combine climate-change refugia modeling with structured decision-making to inform conservation decisions for the endangered foothill yellow-legged frog ( Rana boylii ) in the Sierra Nevada region of California, USA. We used an ensemble of species distribution models to identify areas projected to remain suitable into the 2040s and the 2080s under an RCP 8.5 emissions scenario, as well as areas projected to transition to suitable habitat during this time. We integrated these projections with a structured decision-making process to align management strategies with refugia model outcomes for R. boylii in a subset of the study area. Habitat suitability for R. boylii is projected to decline in the study area by over 90% by the 2040s and by a subsequent 15% by the 2080s. Climate-change refugia are projected to occupy ~7% of present-day suitable habitat, with high agreement between GCMs and model timesteps. Areas projected to transition to suitable habitat within the existing R. boylii clade boundaries are negligible. Collectively, climate-change refugia modeling and structured decision-making provide opportunities to improve resource allocation and empower conservation practitioners in climate change adaptation for at-risk species.

California

Methodology for quantifying biodiversity

Protected and Conserved Areas (PCAs) are usually created and managed for multiple values, with biodiversity conservation being the primary value for most PCAs (see Chapter 1). When using Natural Climate Solutions (NCS) or Ecosystem-based Management (EBM) approaches to climate change mitigation, it is important to ensure that there are no unintended negative consequences on biodiversity and co-benefits for biodiversity are sought. This requires baseline biodiversity information and periodic monitoring. The Convention on Biological Diversity (CBD) has developed a monitoring framework to measure progress (CBD COP 15 2022) towards the GBF. The framework consists of: • headline indicators for national, regional and global monitoring; • global level indicators (collated from yes/no responses in national reports and used to provide a count of the number of countries having undertaken specific activities); • component indicators (which are a list of optional indicators that may apply at global, regional, national and sub-national levels); • complementary indicators (which are a list of optional indicators for thematic or in-depth analysis of each goal and target). Further details of the monitoring framework can be found in Decision 15/5: Monitoring framework for the Kunming-Montreal Global Biodiversity Framework (CBD COP 15, 2022). See Table 5.1 for headline indicators and related targets that they measure. These indicators are used to track national progress towards the GBF. However, those designing sub-national or regional PCA frameworks can use them as well. Biodiversity can be measured at different scales, both spatially and temporally, and at different levels and attributes of biological organisation (Noss, 1990). The Kunming-Montreal Global Biodiversity Framework (GBF) (Convention on Biological Diversity, 2022, December 18) includes goals and targets across scales, with Targets 1 and 3 focused on spatial planning and PCA creation, Targets 2 and 4 focused on restoration and species management to prevent extinction, Target 5 focused on fish stocks, Target 8 on minimising climate change impacts, Target 11 on ecosystem services, and Target 21 on biodiversity information for monitoring the GBF.

IUCN WCPA Protected Area Technical Report Series

The effect of temperature on experimental and natural chemical weathering rates of granitoid rocks

The effects of climatic temperature variations (5-35??C) on chemical weathering are investigated both experimentally using flow-through columns containing fresh and weathered granitoid rocks and for natural granitoid weathering in watersheds based on annual solute discharge. Although experimental Na and Si effluent concentrations are significantly higher in the fresh relative to the weathered granitoids, the proportional increases in concentration with increasing temperature are similar. Si and Na exhibit comparable average apparent activation energies (E(a)) of 56 and 61 kJ/mol, respectively, which are similar to those reported for experimental feldspar dissolution measured over larger temperature ranges. A coupled temperature-precipitation model, using an expanded database for solute discharge fluxes from a global distribution of 86 granitoid watersheds, produces an apparent activation energy for Si (51 kJ/mol), which is also comparable to those derived from the experimental study. This correlation reinforces evidence that temperature does significantly impact natural silicate weathering rates. Effluent K concentrations in the column study are elevated with respect to other cations compared to watershed discharge due to the rapid oxidation/dissolution of biotite. K concentrations are less sensitive to temperature, resulting in a lower average E(a) value (27 kJ/mol) indicative of K loss from lower energy interlayer sites in biotite. At lower temperatures, initial cation release from biotite is significantly faster than cation release from plagioclase. This agrees with reported higher K/Na ratios in cold glacial watersheds relative to warmer temperate environments. Increased release of less radiogenic Sr from plagioclase relative to biotite at increasing temperature produces corresponding decreases in 87Sr/86Sr ratios in the column effluents. A simple mixing calculation using effluent K/Na ratios, Sr concentrations and 87Sr/86Sr ratios for biotite and plagioclase approximates stoichiometric cation ratios from biotite/plagioclase dissolution at warmer temperatures (35??C), but progressively overestimates the relative proportion of biotite with decreasing temperature. Ca, Mg, and Sr concentrations closely correlate, exhibit no consistent trends with temperature, and are controlled by trace amounts of calcite or exchange within weathered biotite. The inability of the watershed model to differentiate a climate signal for such species correlates with the lower temperature dependence observed in the experimental studies.

Geochimica et Cosmochimica Acta

Weathering of slags

Weathering is a natural process causing the transformation of minerals, rocks, and related materials like glass under near-surface conditions. Although metallurgical slags are human-made materials, they also undergo natural weathering processes. As base metal slags weather, the released solutions may contain contaminants that could pose an environmental risk. On the other hand, weathering of ferrous slags is important because of the common reuse of the slags. This chapter summarizes factors that affect rates of slag weathering as well as defining and characterizing the two stages of weathering (incipient and advanced) as they are observed in a wide variety of individual phases and slags. We conclude that slags weather at different scales and rates depending on their chemistry and phase composition as well as the external conditions. Furthermore, weathering of slags is a process that is not perfectly reflected by natural processes because these materials are complex (in terms of texture, composition, and disequilibrium with surficial conditions) and not always analogous to rocks. Therefore, some information in this chapter on weathering is applicable to all slag types, some is relevant to only specific slag types, and some aspects require additional study.

Book chapter

Magnitude and frequency of floods in Nebraska

Estimates of flood characteristics with recurrence intervals up to 100 years can be obtained at most sites in Nebraska by use of techniques presented in this report. Instructions, equations, and graphs are presented to aid the design engineer in estimating the magnitude and frequency of floods. For natural-flow streams, the estimating equations and graphical solutions are based on regional relations between floods of a specific return period and selected basin characteristics. Nebraska was subdivided into five hydrologic regions by means of regressions and residuals from the regressions. The boundaries of Regions 1 and 2 were determined to a great extent by differences in soil type. The divisions between Region 2 and Region 4 were also determined by differences in soil type. Boundaries of Region 3 and Region 5 are along or near basin divides. Flood magnitude and frequency solution diagrams are presented for major controlled streams such as the North Platte, South Platte, Platte, and Republican Rivers. Flood information on small controlled streams is limited to available station data. Flood characteristics are tabulated for 303 gaging stations having 13 or more years of record. These flood characteristics provide the best information on floods at these gaged sites. Observed maximum flood peaks at 303 gaging stations with 13 or more years of record and significant peaks at 57 short-term stations and 31 miscellaneous sites are useful in designing flood-control works for maximum safety from flood damage. Comparison is made with maximum observed floods in the United States.

Nebraska

Nonsteady flow to a well of constant drawdown in an extensive aquifer

A mathematical theory is given for the discharge of a well of constant drawdown, discharging as by natural flow from an effectively infinite aquifer of uniform transmissibility and uniform compressibility. This theory is based on the solution by L. P. Smith of the analogous problem in heat conduction. The mathematical function involved in the solution, which cannot be integrated directly, is evaluated by numerical integration. A table of its values is given for a wide range of its argument. This function is compared with other,asymptotic solutions, and simple, useful approximations are given. Two graphical methods are outlined for determining the coefficients of storage and transmissibility from variations in the rate of discharge of wells flowing at constant drawdown. Data from the Grand Junction, Colo., artesian basin are treated by these methods. In the Grand Junction artesian basin there are about 25 flowing wells ranging in depth from 600 to 1600 ft, most of which obtain water from the Entrada sandstone. A few of the wells obtain water from a sandstone in the overlying Morrison formation and a few tap the underlying Wingate sandstone. The procedure of the tests is outlined, and the “ink-well” mercury gage used to measure the artesian pressures is described. Recovery tests were run on the same wells after the discharge tests. Values of transmissibility obtained from the recovery tests check those obtained by means of the discharge tests.

Colorado

Classification of organic solutes in water by using macroreticular resins

A series of macroreticular resin adsorbents was evaluated for ability to extract and fractionate organic solutes found in natural waters. Studies with organic solute standards and natural water samples lead to the development of dissolved organic carbon (DOC) fractionation analysis whereby the DOC is concentrated and fractionated by macroreticular resins into acid, base, and neutral hydrophobic organic solute fractions and into acid, base, and neutral hydrophilic organic solute fractions. The hydrophobic solutes were sorbed and fractionated on nonionic, nonpolar macroreticular resins, and the hydrophilic solutes were sorbed and fractionated by macroreticular ion-exchange resins. Organic-inorganic complexes as well as organic solutes were isolated by resin sorption. DOC fractionation analysis was developed to serve as a compound classification that fills the gap in organic solute characterization between the organic solute concentration (DOC) and specific compound identification.

Journal of Research of the U.S. Geological Survey

A distribution-based parameterization for improved tomographic imaging of solute plumes

Difference geophysical tomography (e.g. radar, resistivity and seismic) is used increasingly for imaging fluid flow and mass transport associated with natural and engineered hydrologic phenomena, including tracer experiments, in situ remediation and aquifer storage and recovery. Tomographic data are collected over time, inverted and differenced against a background image to produce ‘snapshots’ revealing changes to the system; these snapshots readily provide qualitative information on the location and morphology of plumes of injected tracer, remedial amendment or stored water. In principle, geometric moments (i.e. total mass, centres of mass, spread, etc.) calculated from difference tomograms can provide further quantitative insight into the rates of advection, dispersion and mass transfer; however, recent work has shown that moments calculated from tomograms are commonly biased, as they are strongly affected by the subjective choice of regularization criteria. Conventional approaches to regularization (Tikhonov) and parametrization (image pixels) result in tomograms which are subject to artefacts such as smearing or pixel estimates taking on the sign opposite to that expected for the plume under study. Here, we demonstrate a novel parametrization for imaging plumes associated with hydrologic phenomena. Capitalizing on the mathematical analogy between moment-based descriptors of plumes and the moment-based parameters of probability distributions, we design an inverse problem that (1) is overdetermined and computationally efficient because the image is described by only a few parameters, (2) produces tomograms consistent with expected plume behaviour (e.g. changes of one sign relative to the background image), (3) yields parameter estimates that are readily interpreted for plume morphology and offer direct insight into hydrologic processes and (4) requires comparatively few data to achieve reasonable model estimates. We demonstrate the approach in a series of numerical examples based on straight-ray difference-attenuation radar monitoring of the transport of an ionic tracer, and show that the methodology outlined here is particularly effective when limited data are available.

Geophysical Journal International

Chemistry and transport of soluble humic substances in forested watersheds of the Adirondack Park, New York

Studies were conducted in conjunction with the Integrated Lake-Watershed Acidification Study (ILWAS) to examine the chemistry and leaching patterns of soluble humic substances in forested watersheds of the Adirondack region. During the summer growing season, mean dissolved organic carbon (DOC) concentrations in the ILWAS watersheds ranged from 21–32 mg C l −1 in O/A horizon leachates, from 5–7 mg C l −1 in B horizon leachates, from 2–4 mg C l −1 in groundwater solutions, from 6–8 mg C l −1 in first order streams, from 3–8 mg C l −1 in lake inlets, and from 2–7 mg C l −1 in lake outlets. During the winter, mean DOC concentrations dropped significantly in the upper soil profile. Soil solutions from mixed and coniferous stands contained as much as twice the DOC concentration of lysimeter samples from hardwood stands. Results of DOC fractionation analysis showed that hydrophobia and hydrophilic acids dominate the organic solute composition of natural waters in these watersheds. Charge balance and titration results indicated that the general acid-base characteristics of the dissolved humic mixture in these natural waters can be accounted for by a model organic acid having an average p K a of 3.85, an average charge density of 4–5 μeq mg −1 C at ambient pH, and a total of 6–7 meq COOH per gram carbon.

New York

Soil-solution chemistry in a low-elevation spruce-fir ecosystem, Howland, Maine

Soil solutions were collected monthly by tension and zero-tension lysimeters in a low-elevation red spruce stand in east-central Maine from May 1987 through December 1992. Soil solutions collected by Oa tension lysimeters had higher concentrations of most constituents than the Oa zero-tension lysimeters. In Oa horizon soil solutions growing season concentrations for SO 4 , Ca, and Mg averaged 57, 43, and 30 μmol L −1 in tension lysimeters, and 43, 28, and 19 μmol L −1 in zero-tension lysimeters, respectively. Because tension lysimeters remove water held by the soil at tensions up to 10 kPa, solutions are assumed to have more time to react with the soil compared to freely draining solutions collected by zero-tension lysimeters. Solutions collected in the Bs horizon by both types of collectors were similar which was attributed to the frequency of time periods when the water table was above the Bs lysimeters. Concentrations of SO 4 and NO 3 at this site were lower than concentrations reported for most other eastern U.S. spruce-fir sites, but base cation concentrations fell in the same range. Aluminum concentrations in this study were also lower than reported for other sites in the eastern U.S. and Ca/Al ratios did not suggest inhibition of Ca uptake by roots. Concentrations of SO 4 , Ca, K, and Cl decreased significantly in both the Oa and Bs horizons over the 56-month sampling period, which could reflect decreasing deposition rates for sulfur and base cations, climatic influences, or natural variation. A longer record of measured fluxes will be needed to adequately define temporal trends in solution chemistry and their causes.

Maine

Activity product constant of cryolite at 25°C and one atmosphere using selective-ion electrodes to estimate sodium and fluoride activities

The activity product constant of cryolite (Na 3 AlF 6 ) at 25°C and 1 atm total pressure was calculated from data for solutions from which synthetic cryolite or mixtures of cryolite and a solid apparently related to ralstonite had precipitated. The activities of fluoride and of sodium were estimated using specific ion electrodes. The activity of Al 3+ was calculated from aluminum concentrations using equations which involve aluminum species existing in water solutions containing fluoride and hydroxide ions. For the reaction Na 3 AlF 6 ( c ) = 3 Na ( aq ) + + Al 3+ ( aq ) + 6 E - ( aq ) the activity product constant K so is 10 −33.84 , and the standard free energy of formation at 25°C is -745.4±1.0 kcal mole −1 . This latter is identical to a value from the literature which is based on calorimetric data. Although cryolite and related minerals are easily precipitated in the laboratory, their relative scarcity in nature suggests that the effect of influences, perhaps solutes other than those occurring in the structure, may inhibit or prevent the formation of cryolite.

Geochimica et Cosmochimica Acta

A guide to multi-objective optimization for ecological problems with an application to cackling goose management

Choices in ecological research and management are the result of balancing multiple, often competing, objectives. Multi-objective optimization (MOO) is a formal decision-theoretic framework for solving multiple objective problems. MOO is used extensively in other fields including engineering, economics, and operations research. However, its application for solving ecological problems has been sparse, perhaps due to a lack of widespread understanding. Thus, our objective was to provide an accessible primer on MOO, including a review of methods common in other fields, a review of their application in ecology, and a demonstration to an applied resource management problem. A large class of methods for solving MOO problems can be separated into two strategies: modelling preferences pre-optimization (the a priori strategy), or modelling preferences post-optimization (the a posteriori strategy). The a priori strategy requires describing preferences among objectives without knowledge of how preferences affect the resulting decision. In the a posteriori strategy, the decision maker simultaneously considers a set of solutions (the Pareto optimal set) and makes a choice based on the trade-offs observed in the set. We describe several methods for modelling preferences pre-optimization, including: the bounded objective function method, the lexicographic method, and the weighted-sum method. We discuss modelling preferences post-optimization through examination of the Pareto optimal set. We applied each MOO strategy to the natural resource management problem of selecting a population target for cackling goose ( Branta hutchinsii minima ) abundance. Cackling geese provide food security to Native Alaskan subsistence hunters in the goose's nesting area, but depredate crops on private agricultural fields in wintering areas. We developed objective functions to represent the competing objectives related to the cackling goose population target and identified an optimal solution first using the a priori strategy, and then by examining trade-offs in the Pareto set using the a posteriori strategy. We used four approaches for selecting a final solution within the a posteriori strategy; the most common optimal solution, the most robust optimal solution, and two solutions based on maximizing a restricted portion of the Pareto set. We discuss MOO with respect to natural resource management, but MOO is sufficiently general to cover any ecological problem that contains multiple competing objectives that can be quantified using objective functions.

Ecological Modelling

Indirect spectrophotometric determination of traces of bromide in water

A rapid, accurate, and sensitive indirect spectrophotometric method for the determination of bromide in natural waters is based on the catalytic effect of bromide on the oxidation of iodine to iodate by potassium permanganate in sulfuric acid solution. The method is applicable to concentrations ranging from 1 to 100 μ g. of bromide per liter, but may be modified to extend the concentration range. Most ions commonly occurring in water do not interfere. The standard deviation is 2.9 at bromide concentrations of 100 μ g. per liter and less at lower concentrations. The determination of bromide in samples containing known added amounts gave values ranging from 99 to 105% of the concentration calculated to be present.

Analytical Chemistry

The 25 October 2010 Mentawai tsunami earthquake, from real-time discriminants, finite-fault rupture, and tsunami excitation

The moment magnitude 7.8 earthquake that struck offshore the Mentawai islands in western Indonesia on 25 October 2010 created a locally large tsunami that caused more than 400 human causalities. We identify this earthquake as a rare slow‐source tsunami earthquake based on: 1) disproportionately large tsunami waves; 2) excessive rupture duration near 125 s; 3) predominantly shallow, near‐trench slip determined through finite‐fault modeling; and 4) deficiencies in energy‐to‐moment and energy‐to‐duration‐cubed ratios, the latter in near‐real time. We detail the real‐time solutions that identified the slow‐nature of this event, and evaluate how regional reductions in crustal rigidity along the shallow trench as determined by reduced rupture velocity contributed to increased slip, causing the 5–9 m local tsunami runup and observed transoceanic wave heights observed 1600 km to the southeast.

Mentawai Islands

Bait preference by the Argentine ant (Hymenoptera: Formicidae) in Haleakala National Park, Hawaii

The Argentine ant, Linepithema humile (Mayr), has proven to be a threat to native arthropod species in Haleakala National Park, Maui, HI, and is also a potential threat to the park's native flora. As it continues to expand its range, an effort has been undertaken to eradicate it, or at the least, control its spread. The 1st part of this effort focused on finding a bait carrier for subsequent toxicant-based control tests. A year-long bait preference test was implemented at each of the ant's 2 infestation sites in Haleakala National Park, in which 6 solid baits and 2 liquid baits were assessed for attractiveness and feasibility for large scale control. At both sites, a toxicant-free formulation of Maxforce, a protein-based granular bait made from ground silkworm, Bombyx mori (L.), pupae, and a 25% sugar water solution were the most attractive baits. Ants took more Maxforce (without toxicant) and sugar water than all other baits, including honey granules and a fish protein bait. Sugar water, however, is difficult to distribute over large natural areas. Maxforce was therefore concluded to be the best bait carrier for toxicant-based control at Haleakala National Park because of its attractiveness and its ease for large scale broadcast dispersal.

Hawaii

Molybdenum isotope fractionation during adsorption to organic matter

Organic matter is of emerging interest as a control on molybdenum (Mo) biogeochemistry, and information on isotope fractionation during adsorption to organic matter can improve interpretations of Mo isotope variations in natural settings. Molybdenum isotope fractionation was investigated during adsorption onto insolubilized humic acid (IHA), a surrogate for organic matter, as a function of time (2–170 h) and pH (2–7). For the time series experiment performed at pH 4.2, the average Mo isotope fractionation between the solution and the IHA (Δ 98 Mo solution-IHA ) was 1.39‰ (± 0.16‰, 2σ, based on 98 Mo/ 95 Mo relative to the NIST 3134 standard) at steady state. For the pH series experiment, Mo adsorption decreased as pH increased from 2.0 to 6.9, and the Δ 98 Mo solution-IHA increased from 0.82‰ to 1.79‰. We also evaluated natural Mo isotope patterns in precipitation, foliage, organic horizon, surface mineral soil, and bedrock from 12 forested sites in the Oregon Coast Range. The average Mo isotope offset observed between precipitation and organic (O) horizon soil was 2.1‰, with light Mo isotopes adsorbing preferentially to organic matter. Fractionation during adsorption to organic matter is similar in magnitude and direction to prior observations of Mo fractionation during adsorption to Fe- and Mn- (oxyhydr)oxides. Our finding that organic matter influences Mo isotope composition has important implications for the role of organic matter as a driver of trace metal retention and isotopic fractionation.

Geochimica et Cosmochimica Acta

Effects of karst and geologic structure on the circulation of water and permeability in carbonate aquifers

The results of the natural processes caused by solution and leaching of limestone, dolomite, gypsum, salt and other soluble rocks, is known as karst. Development of karst is commonly known as karstification, which may have a pronounced effect on the topography, hydrology and environment, especially where such karst features as sinkholes and vertical solution shafts extend below the land surface and intersect lateral solution passages, cavities, caverns and other karst features in carbonate rocks. Karst features may be divided into two groups: (1) surficial features that do not extend far below the surface; and (2) karst features such as sinkholes that extend below the surface and affect the circulation of water below. The permeability of the most productive carbonate aquifers is due chiefly to enlargement of fractures and other openings by circulation of water. Important controlling factors responsible for the development of karst and permeability in carbonate aquifers include: (1) climate, topography, and presence of soluble rocks; (2) geologic structure; (3) nature of underground circulation; and (4) base level. Another important factor is the condition of the surface of the carbonate rocks at the time they are exposed to meteoric water. A carbonate rock surface, with soil or relatively permeable, less soluble cover, is more favorable for initiation of karstification and solution than bare rocks. Water percolates downward through the cover to the underlying carbonate rocks instead of running off on the surface. Also, the water becomes more corrosive as it percolates through the permeable cover to the underlying carbonate rocks. Where there is no cover or the cover has been removed, the carbonate rocks become case hardened and resistant to erosion. However, in regions underlain not only by carbonate rocks but also by beds of anhydrite, gypsum and salt, such as the Hueco Plateau in southeastern New Mexico, subsurface solution may occur where water without natural acids moves down from bare rock surfaces through cracks to the beds that are more soluble than carbonate rocks. For example, in the area of Carlsbad Caverns in southeastern New Mexico, much of the water responsible for solution that formed the caverns apparently entered the groundwater system through large open fractures and did not form sinkhole topography. East of the Carlsbad Caverns, however, in the Pecos River Valley where the carbonate rocks are overlain by the less soluble Ogallala Formation of Late Tertiary age, solution began along escarpments as the Pecos River and its tributaries cut through the less soluble cover. As these escarpments retreated, sinkholes and other karst features developed. Joints or fractures are essential for initiation of downward percolation of water in compact carbonate rocks such as some Paleozoic limestone in which there is no intergranular permeability. Also joints or fractures and bedding planes may be essential in the initiation of lateral movement of water in the zone of saturation. Where conditions of recharge and discharge are favorable, groundwater may move parallel to the dip. However, the direction of movement of water in most carbonate rocks is not necessarily down dip or parallel to the dip. The general direction of movement of both surface and groundwater may be parallel to the strike in a breached anticline. Faults may restrict the lateral movement of water, especially if water-bearing beds are faulted against relatively impervious beds. Conversely, some fault may serve as avenues through which water may move as, for example, in the Cretaceous Edwards aquifer in the San Antonio area, Texas. Karst aquifers, chiefly carbonate rocks, may be placed in three groups according to water-bearing capacity. Water in aquifers of group 1 occurs chiefly in joints, fractures, and other openings that have not been enlarged by solution. The yield of wells is small. Aquifers in group 2 , with low to intermediate yields, are those in which water occurs in joints and fractures with some cavities and channels enlarged by solution. Aquifers in group 3 are those in which the yield of wells and springs range from intermediate to very large. This group includes five of the most productive aquifers in the U.S.A. The water-bearing beds of all of these productive aquifers, except the Biscayne aquifer in southeastern Florida, contain buried paleokarst in which the permeability has been reactivated and enlarged by the present circulation system.

Journal of Hydrology