Heterogeneous distributions of CO2 may be more important for dissolution and karstification in coastal eogenetic limestone than mixing dissolution

Mixing dissolution, a process whereby mixtures of two waters with different chemical compositions drive undersaturation with respect to carbonate minerals, is commonly considered to form cavernous macroporosity (e.g. flank margin caves and banana holes) in eogenetic karst aquifers. On small islands, macroporosity commonly originates when focused dissolution forms globular chambers lacking entrances to the surface, suggesting that dissolution processes are decoupled from surface hydrology. Mixing dissolution has been thought to be the primary dissolution process because meteoric water would equilibrate rapidly with calcium carbonate as it infiltrates through matrix porosity and because pCO₂ was assumed to be homogeneously distributed within the phreatic zone. Here, we report data from two abandoned well fields in an eogenetic karst aquifer on San Salvador Island, Bahamas, that demonstrate pCO₂ in the phreatic zone is distributed heterogeneously. The pCO₂ varied from less than log ?2.0 to more than log ?1.0 atm over distances of less than 30 m, generating dissolution in the subsurface where water flows from regions of low to high pCO₂ and cementation where water flows from regions of high to low pCO_2. Using simple geochemical models, we show dissolution caused by heterogeneously distributed pCO₂ can dissolve 2.5 to 10 times more calcite than the maximum amount possible by mixing of freshwater and seawater. Dissolution resulting from spatial variability in pCO₂ forms isolated, globular chambers lacking initial entrances to the surface, a morphology that is characteristic of flank margin caves and banana holes, both of which have entrances that form by erosion or collapse after cave formation. Our results indicate that heterogeneous pCO₂, rather than mixing dissolution, may be the dominant mechanism for observed spatial distribution of dissolution, cementation and macroporosity generation in eogenetic karst aquifers and for landscape development in these settings. Copyright © 2015 John Wiley & Sons, Ltd.     

Speleothems in a north Cuban cave register sea‐level changes and Pleistocene uplift rates

A flight of marine terraces along the Cuban coast records Quaternary sea‐level highstands and a general slowly uplifting trend during the Pleistocene. U/Th dating of these limestone terraces is difficult because fossil reef corals have been affected by open system conditions. Terrace ages are thus often based on geological and geomorphological observations. In contrast, the minimum age of the terraces can be constrained by dating speleothems from coastal mixing (flank margin) caves formed during past sea‐level highstands and carving the marine limestones. Speleothems in Santa Catalina Cave have ages >360 ka and show various cycles of subaerial–subaqueous corrosion and speleothem growth. This suggests that the cave was carved during the MIS 11 sea‐level highstand or earlier. Some stalagmites grew during MIS 11 through MIS 8 and were submerged twice, once at the end of MIS 11 and then during MIS 9. Phreatic overgrowths (POS) covering the speleothems suggest anchialine conditions in the cave during MIS 5e. Their altitude at 16 m above present sea level indicates a late Pleistocene uplift rate of <0.1 mm/ka, but modelling also shows uplift to have been insignificant over a long timespan during the middle Pleistocene since the cave was carved. Our study shows that some flank margin caves in the region of Matanzas are older than commonly believed (i.e. MIS 11 rather than MIS 5). These caves not only can be preserved but are good markers of interglacial sea‐level highstands, more reliable than marine abrasion surfaces. Copyright © 2018 John Wiley & Sons, Ltd.     

Coastal and inland karst morphologies driven by sea level stands: a GIS based method for their evaluation

Sea level is the base level for groundwater circulation in coastal aquifers. The evolution of karst surface landforms and subsurface drainage systems in these aquifers has been conditioned in geological time by tectonics and glacio‐eustatic sea‐level changes. Present morpho‐structural settings and the type/distribution of karst surface and subsurface forms have developed in different carbonate formations according to differences in lithology, climate and exposure time, all driving the intensity of morphologic and karst processes. The repeated and significant changes of groundwater level linked to ‘sea‐level changes’ have had the most important role in driving the continuous evolution of karstic drainage systems, and has resulted in most cases in a multiphase karst. This study aims at defining a general method for identifying, in karst coastal settings, the elevations of flat or low topographic gradient surfaces (using morphometric analysis of Digital Elevation Models (DEMs) and geographical information systems (GISs), and their comparison with elevations of distinctive karstic levels (passages, lateral solution cavities) observed in vertical shafts and horizontal caves. Of the elevations of flat or low topographic gradient surfaces only those agreeing, within ±10 m or ±20 m, with elevation ranges marked by the high frequency of distinctive karst levels were considered as representative of the more probable past sea‐level stands. The method is applied to a regional coastal carbonate formation in southern Italy, by using a 10 m DEM and information on 140 complex caves and 85 shafts. Of the 15 elevations indicated by DEM analysis [620, 600, 470, 450, 425, 385, 355, 315, 270, 250, 205, 180, 150, 110, and 70 m above sea level (a.s.l.)], 13 match clearly those highlighted by significant frequencies of distinctive karstic levels. These elevations are validated by comparison to the elevation of terraces and karst plains indicated in the literature. Copyright © 2012 John Wiley & Sons, Ltd.     

Formation of phreatic caves in an eogenetic karst aquifer by CO2 enrichment at lower water tables and subsequent flooding by sea level rise

Formation of extensive phreatic caves in eogenetic karst aquifers is widely believed to require mixing of fresh and saltwater. Extensive phreatic caves also occur, however, in eogenetic karst aquifers where fresh and saltwater do not mix, for example in the upper Floridan aquifer. These caves are thought to have formed in their modern settings by dissolution from sinking streams or by convergence of groundwater flow paths on springs. Alternatively, these caves have been hypothesized to have formed at lower water tables during sea level low‐stands. These hypotheses have not previously been tested against one another. Analyzing morphological data and water chemistry from caves in the Suwannee River Basin in north‐central Florida and water chemistry from wells in the central Florida carbonate platform indicates that phreatic caves within the Suwannee River Basin most likely formed at lower water tables during lower sea levels. Consideration of the hydrological and geochemical constraints posed by the upper Floridan aquifer leads to the conclusion that cave formation was most likely driven by dissolution of vadose CO₂ gas into the groundwater. Sea level rise and a wetter climate during the mid‐Holocene lifted the water table above the elevation of the caves and placed the caves tens of meters below the modern water table. When rising water tables reached the land surface, surface streams formed. Incision of surface streams breached the pre‐existing caves to form modern springs, which provide access to the phreatic caves. Phreatic caves in the Suwannee River Basin are thus relict and have no causal relationship with modern surficial drainage systems. Neither mixing dissolution nor sinking streams are necessary to form laterally extensive phreatic caves in eogenetic karst aquifers. Dissolution at water tables, potentially driven by vadose CO₂ gas, offers an underappreciated mechanism to form cavernous porosity in eogenetic carbonate rocks. Copyright © 2012 John Wiley & Sons, Ltd.     

Organic carbon inputs,common ions and degassing: rethinking mixing dissolution in coastal eogenetic carbonate aquifers

Caves deliver freshwater from coastal carbonate landscapes to estuaries but how these caves form and grow remains poorly understood. Models suggest fresh and salt water mixing drives dissolution in eogenetic limestone, but have rarely been validated through sampling of mixing waters. Here we assess controls on carbonate mineral saturation states using new and legacy geochemical data that were collected in vertical profiles through three cenotes and one borehole in the Yucatan Peninsula. Results suggest saturation states are primarily controlled by carbon fluxes rather than mixing. Undersaturation predicted by mixing models that rely on idealized end members is diminished or eliminated when end members are collected from above and below actual mixing zones. Undersaturation due to mixing is limited by CO₂ degassing from fresh water in karst windows, which results in calcite supersaturation. With respect to saline groundwater, controls on capacity for mixing dissolution were more varied. Oxidation of organic carbon increased pCO₂ of saline groundwater in caves (pCO₂ = 10^–2.06 to 10^–0.96 atm) relative to matrix porosity (10^–2.39 atm) and local seawater (10^–3.12 atm). The impact of increased pCO₂ on saturation state, however, depended on the geochemical composition of the saline water and the magnitude of organic carbon oxidation. Carbonate undersaturation due to mixing was limited where gypsum dissolution (Cenote Angelita) or sulfate reduction (Cenote Calica) increased concentrations of common ions (Ca²⁺ or HCO₃^?, respectively). Maximum undersaturation was found to occur in mixtures including saline water that had ion concentrations and ratios similar to seawater, but with moderately elevated pCO₂ (Cenote Eden). Undersaturation, however, was dominated by the initial undersaturation of the saline end member, mixing was irrelevant. Our results add to a growing body of literature that suggests oxidation of organic carbon, and not mixing dissolution, is the dominant control on cave formation and enlargement in coastal eogenetic karst aquifers. Copyright © 2016 John Wiley & Sons, Ltd.     

Assessing groundwater availability of the Maldives under future climate conditions

Groundwater resources of the Republic of the Maldives are threatened by a variety of factors including variable future rainfall patterns, continued population growth and associated pumping demands, rising sea level, and contamination from the land surface. This study assesses changes in groundwater availability due to variable rainfall patterns and sea level rise (SLR) in the coming decades, a key component of water resources management for the country. Using a suite of two‐dimensional density‐dependent groundwater flow models, time‐dependent thickness of the freshwater lens is simulated for a range of island sizes (200 to 1,100 m) during the time period of 2011 to 2050, with recharge to the freshwater lens calculated using rainfall patterns provided by general circulation models for the three distinct geographic regions of the Maldives. The effect of SLR on the freshwater lens is quantified using estimates of shoreline recession and associated decreases in island width. If rainfall is solely considered, groundwater availability is projected to increase, as lens thickness during the 2031–2050 time periods is slightly greater (1–5%) than during the 2011–2030 time period. However, including the impact of SLR indicates an overall decrease in lens thickness, with drastic decreases (60% to 100%) projected for small islands (200 m) and moderate decreases (12% to 14%) expected for 400 m islands, which accommodate one third of the national population. Similar methodologies can be used for other atoll island nations, such as the Republic of Marshall Islands, Federated States of Micronesia, and the Republic of Kiribati. For the Maldives, results from this study can be used in conjunction with population growth estimates to determine the feasibility of including groundwater in water resources planning and management for the country.     

Quantifying threats to groundwater resources in the Republic of Maldives Part I: Future rainfall patterns and sea‐level rise

Reserves of fresh groundwater on atoll islands are extremely fragile due to climatic and anthropogenic stresses. Of major concern is the quantity of water to be available in the coming decades under the influence of variable rainfall patterns, rising sea level, environmental conditions, and expected population growth that depends on groundwater resources. In this study, a 3‐dimensional numerical modelling approach using the SEAWAT modelling code is used to estimate freshwater lens volume fluctuation for 4 representative islands in the Republic of Maldives in response to long‐term changes in rainfall, sea‐level rise (SLR), and anthropogenic stresses such as groundwater pumping and short‐term impacts from tsunami‐induced marine overwash events. This work is divided into 2 papers. This first paper presents numerical model set‐up and calibration, and the effect of future rainfall patterns and SLR on fresh groundwater reserves. The second paper focuses on marine overwash events. The results of simulated future freshwater lens volume presented in the first study contribute to efficient groundwater resources planning and management for the Maldives in the upcoming decades. Freshwater lenses in small atoll islands (area < 0.6 km²) are shown to have a strong variability trends in the upcoming decades with expected reduction in lens volume between 11% and 36% due to SLR. In contrast, freshwater lenses in larger atoll islands (area > 1.0 km²) are shown to have less variability to changing patterns with expected reduction in lens volume between 8% and 26% due to SLR. Study results can provide water resource managers with valuable findings for consideration in water security measures.     

Tidal pumping and biogeochemical processes: Dissolution within the tidal capillary fringe of eogenetic coastal carbonates

Growing evidence suggests microbial respiration of dissolved organic carbon (DOC) may be a principal driver of subsurface dissolution and cave formation in eogenetic carbonate rock. Analyses of samples of vadose zone gasses, and geochemical and hydrological data collected from shallow, uncased wells on San Salvador Island, Bahamas, suggest tidally varying water tables may help fuel microbial respiration and dissolution through oxygenation. Respiration of soil organic carbon transported to water tables generates dysaerobic to anaerobic groundwater, limiting aerobic microbial processes. Positive correlations of carbon dioxide (CO₂), radon-222 (²²²Rn) and water table elevation indicate, however, that tidal pumping of water tables pulls atmospheric air that is rich in oxygen, and low in CO₂ and ²²²Rn, into contact with the tidal capillary fringe during falling tides. Ratios of CO₂ and O₂ in vadose gas relative to the atmosphere indicate this atmospheric oxygen fuels respiration within newly-exposed, wetted bedrock. Deficits of expected CO₂ relative to O₂ concentrations indicate some respired CO₂ is likely removed by carbonate mineral dissolution. Tidal pumping also appears capable of transferring oxygen to the freshwater lens, where it could also contribute to respiration and dissolution; dissolved oxygen concentrations at the water table are at least 5% saturated and decline to anaerobic conditions 1–2 m below. Our results demonstrate how tidal pumping of air to vadose zones can drive mineral dissolution reactions that are focused near water tables and may contribute to the formation of laterally continuous vuggy horizons and potentially caves. © 2020 John Wiley & Sons, Ltd.     

Reconciliation of sea-level observations in the Western North Atlantic during the last glacial cycle

A south to north gradient of increasing marine isotope substage (MIS) 5a (∼80 ka BP) sea level has been recorded across the Caribbean and surrounding region. Relative to present, MIS-5a sea levels range from −19 m to more than +3 m between Barbados, Haiti, the Bahamas, Florida, Bermuda and the US Atlantic Coast. In contrast, no gradient in sea level is observed for the last interglacial period MIS-5e (∼128–118 ka BP) at tectonically stable localities in the same region, with deposits generally lying several metres above present. We demonstrate here that these controversial observations are reconciled by taking into account the isostatic response of the Earth to glacial loading and unloading – a fundamental effect that is commonly overlooked in the interpretation of sea-level observations from different locations to define a ‘global sea-level curve’. Furthermore, the observed gradient can be used to place constraints on Earth rheology and is an important indicator of the behaviour of the North American ice sheets during the last glacial cycle.     

Sources of ground water salinity on islands using 18O, 2H, and 34S

Stable isotopes of 18O and 2H in water, and 34S and 18O in dissolved SO4, are used to verify the interpretation of the chemical evolution and proposed sources of salinity for two islands that have undergone postglacial rebound. Results for delta18O and delta34S in dissolved SO4 on the Gulf Islands, southwest British Columbia, Canada, suggest a three-component mixing between (1) atmospheric SO4 derived largely from recharge of meteoric origin, (2) modern marine SO4 associated with either modern-day salt water intrusion or Pleistocene age sea water, and (3) terrestrial SO4. The age of the marine SO4 is uncertain based on the geochemistry and SO4 isotopes alone. Two options for mixing of saline ground waters are proposed--either between current-day marine SO4 and atmospheric SO4, or between older (Pleistocene age) marine SO4 and atmospheric SO4, delta18O and delta2H compositions are relatively consistent between both islands, with a few samples showing evidence of mixing with water that is a hybrid mixture of Fraser River water and ocean water. The isotopic composition of this hybrid water is approximately delta18O = 10 per thousand. delta18O and delta2H values for many saline ground waters plot close to the global meteoric water line, which is distinctly different from the local meteoric water line. This suggests a meteoric origin for ground waters that is different from the current isotopic composition of meteoric waters. It is proposed these waters may be late Pleistocene in age and were recharged when the island was submerged below sea level and prior to rebound at the end of the last glaciation.     

Predicting Future Groundwater Resources of Coral Atoll Islands

Fresh groundwater reserves on small coral islands are under continual threat of salinization and contamination because of droughts, storm‐surge overwash events, over‐extraction, island community urbanization, and sea level rise. Whereas storm‐surge overwash events can cause sudden groundwater salinization, long‐term changes in rainfall patterns and sea level elevation have the potential of rendering these islands uninhabitable in the coming decades. This study demonstrates the use of a tested freshwater lens thickness simulator to estimate the groundwater resources of a set of atoll islands in the coming decades. The method uses ranges of projected rates of annual rainfall and sea level rise (SLR) to provide a range of probable lens thickness for each island. Projected rainfall is provided by General Circulation Models that accurately replicate the historical rainfall patterns in the geographic region of the islands. Methodology is applied to 68 atoll islands in the Federated States of Micronesia. These islands have widths that range between 150 and 1000 m, and experience annual rainfall rates of between 2.8 and 4.8 m. Results indicate that under average conditions of SLR, beach slope, and rainfall, almost half of the island will experience a 20% decrease in lens thickness by the year 2050. For worst‐case scenarios (high SLR, low rainfall), average decrease in lens thickness is 55%, with almost half of the islands experiencing a decrease of greater than 75% and half of the islands having a lens thickness less than 1.0 m. Small islands (widths less than 400 m) are particularly vulnerable because of shoreline recession. Groundwater on islands in the western region is less vulnerable to SLR because of a projected increase in rainfall during the coming decades. Results indicate the vulnerability of small islands to changing climatic conditions, and can be used for water resources management and community planning. Methodology can be applied to any group of islands as a first approximation of the effect of future climate conditions on groundwater resources. Copyright © 2016 John Wiley & Sons, Ltd.     

Geochemistry of meteoric diagenesis in carbonate islands of the northern Bahamas: 1. Evidence from field studies

Processes driving carbonate diagenesis in islands of the northern Bahamas are investigated using major ion, dissolved oxygen and dissolved organic carbon analyses of water samples from surface and ground waters, and measurements of soil gas P. Meteoric waters equilibrate with aragonite, but reactions are water controlled rather than mineral‐controlled and drive dissolution rather than concurrent precipitation of calcite. Surface runoff waters equilibrate with atmospheric P and rapidly recharge the vadose zone, limiting subaerial bedrock dissolution to only 6·6–15 mg l^?1 Ca. P of soil gas measured in the summer wet season ((7·4 ± 3·7) × 10^?3 atm) is elevated compared with that of the atmosphere, despite the thin skeletal organic nature of the soil and the discontinuous soil cover. Soil waters retained in surface pockets are equilibrated with respect to aragonite and have dissolved 51 ± 19 mg l^?1 Ca. This is substantially less than the 93 ± 18 mg l^?1 Ca in samples from pumping boreholes that sample meteoric waters from the freshwater lens. The high P of the freshwater lens ((16 ± 8·3) × 10^?3 atm for pumping boreholes) suggests that significant additional CO₂ may be derived by oxidation of soil‐ and surface‐derived organic carbon within the lens. The suboxic nature of the majority of the freshwater lens and the observed depletion in sulphate support this suggestion, and indicate that both aerobic and anaerobic oxidation may take place. Shallow lens samples from observation boreholes are calcite supersaturated and have a lower P than deeper lens waters, indicating that CO₂ degasses from the water table, driving precipitation of calcite cements. We suggest that the geochemical evolution of waters in the vadose zone and upper part of the freshwater lens may be determined by the presence of a body of ground air with P controlled by production in the freshwater lens and soil and by degassing to the atmosphere. Copyright © 2007 John Wiley & Sons, Ltd.     

Chronostratigraphy and morphological changes in Cerion land snail shells over the past 130 ka on Long Island,Bahamas

Despite the nearly 600 named species of the land snail Cerion, studies of the geological and paleontological framework of modern species are few. To address this deficiency, the biostratigraphic succession of Cerion was investigated at several areas on Long Island, Bahamas. A chronostratigraphic framework was developed through whole-rock and Cerion land snail aminostratigraphies. About 175 individual Cerion shells from last interglacial and Holocene deposits were age-ranked using stratigraphic position and amino acid racemization (AAR) geochronology. AAR ages were generated using an existing AAR-¹⁴C age model for Cerion from the central Bahamas. The age structure of Cerion fossils in sediments was determined with AAR ages, and the magnitude of “dead carbon” anomalies was evaluated using this chronological approach.Temporal changes in gross shell morphology were examined from four study areas. The last interglacial, marine isotope stage/substage (MIS) 5e (Aminozone E) is characterized by generally large shells and in some cases, bimodal sets of very small (α shells) and very large forms (β shells) coexisting in the same stratigraphic levels (primarily soils), which may encompass the transition from between MIS 5e and 5d/c. Similar bimodality of nearly identical α and β shell forms and sizes is observed at other late MIS 5e sites from the furthest reaches of Great Bahama Bank (including Long, Exumas, Eleuthera, and New Providence Islands). The widespread distribution of α and β forms in soils capping MIS 5e marine and eolian deposits implies that there may have been a synchronous, regional morphological convergence on Great Bahama Bank. None of these forms are observed in Holocene deposits of Aminozone A.The earliest MIS 1 Cerion appear in a oolite deposited 6500 a BP, and are of intermediate size compared to the Pleistocene α and β forms. As MIS 1 progressed, the diversity of shell sizes and shapes increased into modern times. The greater variety of shell forms over the past 1000–2000 a suggests that humans may have played a role in the introduction and redistribution of Cerion across the region. The potential for frequent and widespread human introductions, combined with the propensity of Cerion to hybridise freely may explain the farrago of shell sizes and shapes in the recent snail faunas of Long Island and other Bahama islands.     

Speleogenesis of an exhumed hydrothermal sulphuric acid karst in Cambrian carbonates (Mount San Giovanni,Sardinia)

In the past few years the systematic study of caves intercepted by mine workings in southwest Sardinia has permitted us to observe morphologies due to rare speleogenetic and minerogenetic processes related to ancient hydrothermal activity. These relic morphologies are slowly being overprinted by recent speleogenetic processes that tend to obscure the hypogene origin of these caves. A combined geomorphological and mineralogical investigation has permitted a fairly detailed reconstruction of the various phases of evolution of these caves. Cave formation had already started in Cambrian times, but culminated in the Carboniferous, when most of the large voids still accessible today were formed. A key role in carbonate dissolution was played by sulphuric acid formed by the oxidation of the polymetallic ores present in the rocks since the Cambrian. During the Quaternary a variety of minerals formed inside the caves: calcite and aragonite, that yielded sequences of palaeo‐environmental interest, and also barite, phosgenite, hydrozincite, hemimorphite and many others. These minerals are in part due to a phreatic thermal hypogenic cave forming phase, and in part to later epigene overprinting in an oxidizing environment rich in polymetallic ores. Massive gypsum deposits, elsewhere typical of this kind of caves, are entirely absent due to dissolution during both the phreatic cave formation and the later epigenic stage. Copyright © 2013 John Wiley & Sons, Ltd.     

Vadose CO2 gas drives dissolution at water tables in eogenetic karst aquifers more than mixing dissolution

Most models of cave formation in limestone that remains near its depositional environment and has not been deeply buried (i.e. eogenetic limestone) invoke dissolution from mixing of waters that have different ionic strengths or have equilibrated with calcite at different _pCO₂ values. In eogenetic karst aquifers lacking saline water, mixing of vadose and phreatic waters is thought to form caves. We show here calcite dissolution in a cave in eogenetic limestone occurred due to increases in vadose CO₂ gas concentrations and subsequent dissolution of CO₂ into groundwater, not by mixing dissolution. We collected high‐resolution time series measurements (1 year) of specific conductivity (SpC), temperature, meteorological data, and synoptic water chemical composition from a water table cave in central Florida (Briar Cave). We found SpC, _pCO₂ and calcite undersaturation increased through late summer, when Briar Cave experienced little ventilation by outside air, and decreased through winter, when increased ventilation lowered cave CO_2(g) concentrations. We hypothesize dissolution occurred when water flowed from aquifer regions with low _pCO₂ into the cave, which had elevated _pCO₂. Elevated _pCO₂ would be promoted by fractures connecting the soil to the water table. Simple geochemical models demonstrate that changes in _pCO₂ of less than 1% along flow paths are an order of magnitude more efficient at dissolving limestone than mixing of vadose and phreatic water. We conclude that spatially or temporally variable vadose CO_2(g) concentrations are responsible for cave formation because mixing is too slow to generate observed cave sizes in the time available for formation. While this study emphasized dissolution, gas exchange between the atmosphere and karst aquifer vadose zones that is facilitated by conduits likely exerts important controls on other geochemical processes in limestone critical zones by transporting oxygen deep into vadose zones, creating redox boundaries that would not exist in the absence of caves. Copyright © 2014 John Wiley & Sons, Ltd.     

Granitic coastal geomorphology: applying integrated terrestrial and bathymetric LiDAR with multibeam sonar to examine coastal landscape evolution

Coasts composed of resistant lithologies such as granite are generally highly resistant to erosion. They tend to evolve over multiple sea level cycles with highstands acting to remove subaerially weathered material. This often results in a landscape dominated by plunging cliffs with shore platforms rarely occurring. The long‐term evolution of these landforms means that throughout the Quaternary these coasts have been variably exposed to different sea level elevations which means erosion may have been concentrated at different elevations from today. Investigations of the submarine landscape of granitic coasts have however been hindered by an inability to accurately image the nearshore morphology. Only with the advent of multibeam sonar and aerial laser surveying can topographic data now be seamlessly collected from above and below sea level. This study tests the utility of these techniques and finds that very accurate measurements can be made of the nearshore thereby allowing researchers to study the submarine profile with the same accuracy as the subaerial profile. From a combination of terrestrial and marine LiDAR data with multibeam sonar data, it is found that the morphology of granite domes is virtually unaffected by erosion at sea level. It appears that evolution of these landscapes on the coast is a very slow process with modern sea level acting only to remove subaerially weathered debris. The size and orientation of the joints determines the erosional potential of the granite. Where joints are densely spaced (<2 m apart) or the bedrock is highly weathered can semi‐horizontal surfaces form. Copyright © 2014 John Wiley & Sons, Ltd.     

Constraining Holocene sea levels using U‐Th ages of phreatic overgrowths on speleothems from coastal caves in Mallorca (Western Mediterranean)

Phreatic overgrowths on speleothems (POS) are carbonate formations deposited at the water table of caves in unique karstic coastal settings having morphologies that can be directly related to sea level at the time of formation. The U‐Th ages of calcite and aragonite overgrowths collected from the modern water table in coastal caves on Mallorca (Cova de Cala Varques A and Cova des Pas de Vallgornera) were determined using high‐precision MC‐ICPMS techniques. U‐Th ages indicate that phreatic carbonate deposition occurred between ca 2·8 and at least 0·6 ka BP and are in accord with an archeologically estimated age of 3·7–3·0 ka BP for a drowned prehistoric construction at a depth of 1 m below current sea level in a cave from the same area. Speleothem δ¹³C and δ¹⁸O and chemical composition of cave pools provide supportive evidence that POS reflect mixing between seawater and brackish water table. Copyright © 2010 John Wiley & Sons, Ltd.     

Volcanic geomorphology and tectonics of the Aeolian archipelago (Southern Italy) based on integrated DEM data

The geomorphological and morphometric analysis of the sea floor topography surrounding the Aeolian Islands, South Tyrrhenian Sea, Italy, provides insights into the relationships between the volcanological evolution of the islands and their tectonic features. We constructed geomorphological maps of the submarine portions of the seven large edifices constituting the islands on the basis of a DEM with a 5 m resolution step. These maps include constructional and destructional landforms such as submarine volcanic vents located west of Lipari and north of Alicudi, and hummocky surfaces recognised north of Lipari and Salina. The latter landforms, together with the occurrence of large scars affecting the main edifices on land, suggest that sector collapses affected some islands. Geomorphological data indicate that the location of subaerial and submarine vents is strongly controlled by local tectonic structures striking WNW-ESE (Alicudi-Filicudi sector), NNW-SSE (Salina-Lipari-Vulcano sector) and NE-SW (Panarea-Stromboli sector). The islands can be divided into two groups on the basis of some morphometric parameters: a first group with a pancake-like shape, Dp/D (abrasion platform diameter/basal diameter) higher than 0.40 and H/D (total height/basal diameter) lower than 0.13, and a second group with a conical shape, characterised by Dp/D lower than 0.34 and H/D higher than 0.14. These ratios and other morphometric parameters reflect the different volcanological and structural evolution of the Aeolian Islands. The pancake-like shaped complexes have been created, in addition to their submarine stage, by extrusive and highly explosive activity, whereas the cone-shaped edifices have been characterised by effusive or moderate explosive activity.Editorial responsibility: C Kilburn     

Source speciation resolving hydrochemical complexity of coastal aquifers

There is a growing concern of seawater intrusion to freshwater aquifers due to groundwater overexploitation in the eastern coastal belt of Southern India. The problem becomes complex in the regions where industrial effluents are also contaminating the freshwater aquifers. In order to understand the hydrochemical complexity of the system, topographic elevation, static water level measurements, major ion chemistry, ionic cross plots, water type contours and factor analysis were applied for 144 groundwater samples of shallow and deep sources from Quaternary and Tertiary coastal aquifers, located within the industrial zone of 25 km² area near Cuddalore, Southern India. The ionic cross plots indicates dissolution of halite minerals from marine sources and seawater mixing into inland aquifers up to the level of 9.3%. The factor analysis explains three significant factors totaling 86.3% of cumulative sample variance which includes varying contribution from marine, industrial effluent and freshwater sources.     

Geochemistry of meteoric diagenesis in carbonate islands of the northern Bahamas: 2. Geochemical modelling and budgeting of diagenesis

Geochemical characterization and numerical modelling of surface water and ground water, combined with hydrological observations, provide quantitative estimates of meteoric diagenesis in Pleistocene carbonates of the northern Bahamas. Meteoric waters equilibrate with aragonite, but water‐ rather than mineral‐controlled reactions dominate. Dissolutional lowering of the undifferentiated bedrock surface is an order of magnitude slower than that within soil‐filled topographic hollows, generating small‐scale relief at a rate of 65–140 mm ka^?1 and a distinctive pocketed topography. Oxidation of organic matter within the subsoil and vadose zones generates an average P of 4·0 × 10^?3 atm, which drives dissolution during vadose percolation and/or at the water table. However, these dissolution processes together account for <60% of the average rock‐derived calcium in groundwaters pumped from the freshwater lens. The additional calcium may derive from oxidation of organic carbon within the lens, accounting for the high P of the lens waters. Mixing between meteoric waters of differing chemistry is diagenetically insignificant, but evapotranspiration from the shallow water table is an important drive for subsurface cementation. Porosity generation in the shallow vadose zone averages 1·6–3·2% ka^?1. Phreatic meteoric diagenesis is focused near the water table, where dissolution generates porosity at 1·4–2·8% ka^?1. Maximum dissolution rates, however, are similar to those of evaporation‐driven precipitation, which occludes porosity of 4·0 ± 0·6% ka^?1. This drives porosity inversion, from primary interparticle to secondary mouldic, vug and channel porosity. In the deeper freshwater lens, oxidation of residual organic carbon and reoxidation of reduced sulphur species from deeper anaerobic oxidation of organic carbon may generate porosity up to 0·06% ka^?1. Meteoric diagenesis relies critically on hydrological routing and vadose thickness (controlled by sea level), as well as the geochemical processes active. A thin vadose zone permits direct evaporation from the water table and drives precipitation of meteoric phreatic cements even where mineral stabilization is complete. Copyright © 2007 John Wiley & Sons, Ltd.