Ground‐water velocities within fractures and boreholes, hydraulic heads, and depth profiles of conductivity were measured along a 70 km section of the northeastern coast of the Yucatan Peninsula, Mexico. Hydraulic heads ranged from 40 to 60 cm above mean sea level between 2 and 4 km from the coast. Fluid velocities estimated from point‐dilution tests, in the dual‐porosity rock in a borehole several kilometers from the coast, were 0.021 cm/sec in the fresh‐water lens and 0.082 cm/sec near a fracture in the underlying sea‐water zone. Velocities in large fractures increased from 1 cm/sec 10 kilometers inland to 12 cm/sec near discharge points along the coast. This increase is attributed to the decrease in thickness of the fresh‐water lens. The thickness of the fresh‐water lens is approximately 40% less than the Ghyben‐Herzberg relation predicts for a static system, providing the potential to drive fresh water through fractures into the sea‐water zone below the halocline. Overall, the halocline appears to be in a steady‐state position due to the rapid flow of fresh water and brackish water towards the coast combined with rising sea water in corivectional flow.

Cone-penetrometer testing and computer modeling were utilized to investigate factors controlling fresh-water lens formation at Grand Isle, Louisiana. Measurements of tip resistance, sleeve friction, and electrical conductivity were recorded with depth to permit classification of sediment type and to determine thickness of the fresh-water lens and transition zone. Cone-penetrometer testing provided virtually continuous determinations of change in sediment type and ground-water salinity at a resolution rarely achieved using conventional drilling and water sampling techniques.Three sand bodies are present, each separated by a clay layer. The fresh-water lens is thinner in the center of the island than on the flanks. Fresh-water lens thickness is limited by a clay layer which prohibits downward movement of significant volumes fresh water. The transition zone from fresh water to salt water varies in thickness, being thinnest near the Gulf of Mexico and thickest where silt and clay interfinger with the upper sand.Both the thickness of the fresh-water lens and the shape of the transition zone differ from that predicted by theoretical models. Calibration of STJTRA, a variable-density solute-transport model, indicates that permeability variations are the dominant control on formation of the fresh-water lens. (KEY TERMS: hydrogeology; barrier island; modelling; fresh-water lens; cone penetrometer.)

A temperature spike is reported in the haloclines of three Yucatan sinkholes along a 1 km NW-SE transect from 5 to 4 km inland from the Caribbean coast. The temperature spike decreases in magnitude from 3.5 degrees C to 0.2 degrees C, approaching the coast. The anomaly does not vary diurnally and does not extend down into the underlying sea water. These conditions are inconsistent with explanations such as radiation absorption within the halocline, in situ microbially mediated sulfate reduction within the halocline and the underlying sea water, and sulfide oxidation by photosynthetic purple and green bacteria within the halocline. One explanation consistent with the shape and halocline location of the temperature spike involves a localized sea water convection cell operating near the coast. Cold sea water from the Caribbean Sea enters the coastal limestone at depths of a few hundred meters and heats up because of the geothermal gradient, buoyantly rising in vertical fractures within the unconfined aquifer. Blocked by the less dense fresh water, the movement stops in the halocline where the warm sea water mixes with brackish water. The convection cycle would be completed by the coastward movement of cooling brackish water. The observed temperature anomalies could possibly be snapshots of this warm layer moving coastward.

Geostatistical simulation is used to create a synthetic three‐dimensional geologic environment such that randomly generated values of a physical property have a required covariance and mean. We create a spatially stationary, multivariate Gaussian simulation with a Gaussian covariance. This type of simulation permits modeling of natural variations which may be an important source of noise in geologic investigations, whereas interpolation techniques such as kriging tend to smooth out those variations. Computationally efficient frequency‐domain methods allow the rapid creation of multiple scenarios, each having the required statistical distribution. For the particular problem considered, a given porosity distribution is simulated and then functionally related to physical properties needed for input into geophysical models. Additionally, megascopic cavities having random size, shape, and orientation similar to those observed in karst settings are emplaced in the physical model. Geostatistical simulation plays a particularly important role in exploring the degree to which natural variations in physical parameters lower the resolution of target anomalies that have similar signatures. In the accompanying paper, we discuss the geophysical methods tested and the results of the tests.