Post-processing of sequential indicator simulation realizations for modeling geologic bodies

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Interplay of Late Cenozoic Siliciclastic Supply and Carbonate Response on the Southeast Florida Platform

High-resolution seismic-reflection data collected along the length of the Caloosahatchee River in southwestern Florida have been correlated to nannofossil biostratigraphy and strontium-isotope chemostratigraphy at six continuously cored boreholes. These data are interpreted to show a major Late Miocene(?) to Early Pliocene fluvial– deltaic depositional system that prograded southward across the carbonate Florida Platform, interrupting nearly continuous carbonate deposition since early in the Cretaceous. Connection of the platform top to a continental source of siliciclastics and significant paleotopography combined to focus accumulation of an immense supply of siliciclastics on the southeastern part of the Florida Platform. The remarkably thick (> 100 m), sand-rich depositional system, which is characterized by clinoformal progradation, filled in deep accommodation, while antecedent paleotopography directed deltaic progradation southward within the middle of the present-day Florida Peninsula. The deltaic depositional system may have prograded about 200 km southward to the middle and upper Florida Keys, where Late Miocene to Pliocene siliciclastics form the foundation of the Quaternary carbonate shelf and shelf margin of the Florida Keys. These far-traveled siliciclastic deposits filled accommodation on the southeastern part of the Florida Platform so that paleobathymetry was sufficiently shallow to allow Quaternary recovery of carbonate sedimentation in the area of southern peninsular Florida and the Florida Keys.

Tectonic geomorphology and active strain of the Northern Apennines mountain front

The Northern Apennines (NA) chain is the expression of the active plate margin between Europe and Adria. Given the low convergence rates and the moderate seismic activity, ambiguities still occur in defining a seismotectonic framework and many different scenarios have been proposed for the mountain front evolution. Differently from older models that indicate the mountain front as an active thrust at the surface, a recently proposed scenario describes the latter as the frontal limb of a long-wavelength fold (> 150 km) formed by a thrust fault tipped around 17 km at depth, and considered as the active subduction boundary. East of Bologna, this frontal limb is remarkably very straight and its surface is riddled with small, but pervasive high- angle normal faults. However, west of Bologna, some recesses are visible along strike of the mountain front: these perturbations seem due to the presence of shorter wavelength (15 to 25 km along strike) structures showing both NE and NW-vergence. The Pleistocene activity of these structures was already suggested, but not quantitative reconstructions are available in literature. This research investigates the tectonic geomorphology of the NA mountain front with the specific aim to quantify active deformations and infer possible deep causes of both short- and long-wavelength structures. This study documents the presence of a network of active extensional faults, in the foothills south and east of Bologna. For these structures, the strain rate has been measured to find a constant throw-to-length relationship and the slip rates have been compared with measured rates of erosion. Fluvial geomorphology and quantitative analysis of the topography document in detail the active tectonics of two growing domal structures (Castelvetro - Vignola foothills and the Ghiardo plateau) embedded in the mountain front west of Bologna. Here, tilting and river incision rates (interpreted as that long-term uplift rates) have been measured respectively at the mountain front and in the Enza and Panaro valleys, using a well defined stratigraphy of Pleistocene to Holocene river terraces and alluvial fan deposits as growth strata, and seismic reflection profiles relationships. The geometry and uplift rates of the anticlines constrain a simple trishear fault propagation folding model that inverts for blind thrust ramp depth, dip, and slip. Topographic swath profiles and the steepness index of river longitudinal profiles that traverse the anti- clines are consistent with stratigraphy, structures, aquifer geometry, and seismic reflection profiles. Available focal mechanisms of earthquakes with magnitude between Mw 4.1 to 5.4, obtained from a dataset of the instrumental seismicity for the last 30 years, evidence a clear vertical separation at around 15 km between shallow extensional and deeper compressional hypocenters along the mountain front and adjacent foothills. In summary, the studied anticlines appear to grow at rates slower than the growing rate of the longer- wavelength structure that defines the mountain front of the NA. The domal structures show evidences of NW-verging deformation and reactivations of older (late Neogene) thrusts. The reconstructed river incision rates together with rates coming from several other rivers along a 250 km wide stretch of the NA mountain front and recently available in the literature, all indicate a general increase from Middle to Late Pleistocene. This suggests focusing of deformation along a deep structure, as confirmed by the deep compressional seismicity. The maximum rate is however not constant along the mountain front, but varies from 0.2 mm/yr in the west to more than 2.2 mm/yr in the eastern sector, suggesting a similar (eastward-increasing) trend of the apenninic subduction.

Searching the Rhone delta channel in Lake Geneva since François‐Alphonse Forel

In the late 19th century, F.A. FOREL led investigations of the Rhone River delta area of Lake Geneva that resulted in the dis- covery of a textbook example of a river-fed delta system containing impressive subaquatic channels. Well ahead of the marine counterparts, scientific observations and interpretations of water currents shaping the delta edifice for the first time documented how underflow currents carry cold, suspension-laden waters from the river mouth all the way to the deep basin. These early investigations of the Rhone delta laid the basis for follow-up studies in the 20th and 21th centuries. Sediment coring, water-column measurements, manned submersible diving, seismic reflection profiling and bathymetric sur- veying eventually provided a rich database to unravel the key erosional and depositional processes, further documenting the impact of human-induced changes in the catchment. With the merging of old and new scientific knowledge, today a comprehensive understanding prevails of how a delta changes through time, how its channels are formed, and what potential natural hazards may be related to its evolution. New and efficient bathymetric techniques, paired with novel coring operations, provided a time-series of morphologic evolution showing and quantifying the high dynamics of the delta/channel evolution in an unprecedented temporal and spatial reso- lution. Future investigations will continue to further quantify these dynamic processes and to link the evolution of the subaquatic domain with changes and processes in the catchment and with natural hazards. Its size, easy access, and large variety of states and processes will continue to make the Rhone delta area a perfect ‘laboratory’ in which general processes can be studied that could be upscaled or downscaled to other marine and lacustrine deltas.

The history of subaquatic volcanism recorded in the sediments of Lake Kivu

Subaquatic volcanic activity has been ongoing in Lake Kivu since the early Holocene and has a dynamic effect on the biological productivity in the surface water, and the preservation of carbonate in the deep anoxic water. Groundwater discharge into the lake’s deepwater propels the upward advection of the water column that ultimately supplies nutrients to the surface water for biological production. The amount of nutrients supplied from the deepwater can be increased suddenly by (1) a cold meteorological event that drives deep seasonal mixing resulting in increased nutrients from below and oxygen from above, or (2) subaquatic volcanic activity that induces a buoyant hydrothermal plume, which entrains nutrients from the deepwater and results in anoxia or suboxic conditions in the surface water. Previous sedimentological studies in Lake Kivu have hypothesized that regional climatic changes are responsible for sudden changes in the preservation of carbonates in the Main Basin. Here we reveal that sublacustrine volcanic events most likely induce the abrupt changes to the geochemistry in the sediment in Lake Kivu. An unprecedented look into the sediment stratigraphy and geochemistry from high-resolution seismic-reflection, and 15N-isotope analyses was conducted in the Main Basin. The results reveal that buoyant hydrothermal plumes caused by subaquatic volcanic activity are a possible trigger for increased biological productivity and organic matter preservation, and that ongoing hydrothermal activity increases the alkalinity in the deepwater, leading to carbonate preservation. The onset of carbonate preservation since the 1970s that is currently observed in the sediment could indicate that hydrothermal discharge has recently increased in the lake.

Depositional modes and lake-level variability at Lake Towuti, Indonesia, during the past ~29 kyr BP

Lake Towuti (2.5°S, 121.5°E) is a long-lived, tectonic lake located on the Island of Sulawesi, Indonesia, and in the center of the Indo-Pacific warm pool (IPWP). Lake Towuti is connected with upstream lakes Matano and Mahalona through the Mahalona River, which constitutes the largest inlet to the lake. The Mahalona River Delta is prograding into Lake Towuti’s deep northern basin thus exerting significant control on depositional processes in the basin. We combine high-resolution seismic reflection and sedimentological datasets from a 19.8-m-long sediment piston core from the distal edge of this delta to characterize fluctuations in deltaic sedimentation during the past ~29 kyr BP and their relation to climatic change. Our datasets reveal that, in the present, sedimentation is strongly influenced by deposition of laterally transported sediments sourced from the Mahalona River Delta. Variations in the amount of laterally transported sediments, as expressed by coarse fraction amounts in pelagic muds and turbidite recurrence rates and cumulative thicknesses, are primarily a function of lake-level induced delta slope instability and delta progradation into the basin. We infer lowest lake-levels between ~29 and 16, a gradual lake level rise between ~16 and 11, and high lake-levels between ~11 and 0 kyr BP. Periods of highest turbidite deposition, ~26 to 24 and ~18 to 16 kyr BP coincide with Heinrich events 2 and 1, respectively. Our lake-level reconstruction therefore supports previous observations based on geochemical hydroclimate proxies of a very dry last glacial and a wet Holocene in the region, and provides new evidence of millennial-scale variations in moisture balance in the IPWP.

Facies investigations on sediment core CRP-3 from the Ross Sea, Antarctica

The Cenozoic Victoria Land Basin (VLB) stratigraphic section penetrated by CRP-3 is mostly of Early Oligocene age. It contains an array of lithofacies comprising fine-grained mudrocks, interlaminated and interbedded mudrocks/sandstones, mud-rich and mud-poor sandstones, conglomerates and diamctites that are together interpreted as the products of shallow marine to possibly non-marine environments of deposition, affected by the periodic advance and retreat of tidewater glaciers. This lithofacies assemblage can be readily rationalised using the facies scheme designed originally for CRP-2/2A, and published previously. The uppermost 330 metres below sea floor (mbsf) shows a cyclical arrangement of lithofacies also similar to that recognised throughout CRP-2/2A, and interpreted to reflect cyclical variations in relative sea-level driven by ice volume fluctuations ('Motif A'). Between 330 and 480 mbsf, a series of less clearly cyclical units, generally fining-upward but nonetheless incorporating a significant subset of the facies assemblage, has been identified and noted in the Initial Report as 'Motif B' Below 480 mbsf, the section is arranged into a repetitive succession of fining-upward units, each of which comprises dolerite clast conglomerate at the base passing upward into relatively thick intervals of sandstones. The cycles present down 480 mbsf are defined as sequences, each interpreted to record cyclical variation of relative sea-level. The thickness distribution of sequences in CRP-3 provides some insights into the geological variables controlling sediment accumulation in the Early Oligocene section. The uppermost part of the section in CRP-3 comprises two or three thick, complete sequences that show a broadly symmetrical arrangement of lithofacies (similar to Sequences 9-11 in CRP-2/2A). This suggests a period of relatively rapid tectonic subsidence, which allowed preservation of the complete facies cycle. Below Sequence 3, however, is a considerable interval of thin, incomplete and erosionally truncated sequences (4-23), which incorporates both the remainder of Motif A sequences and all Motif B sequences recognised. The thinner and more truncated sequences suggest sediment accumulation under conditions of reduced accommodation, and given the lack of evidence for glacial conditions (see Powell et al., this volume) tends to argue for a period of reduced tectonic subsidence. The section below 480 mbsf consists of a series of fining-upward, conglomerate to sandstone intervals which cannot be readily interpreted in terms of relative sea-level change. A relatively mudrock-rich interval above the basal conglomerate/breccia (782-762 mbsf) may record initial flooding of the basin during early rift subsidence. The lithostratigraphy summarised above has been linked to seismic reflection data using depth conversion techniques (Henrys et al., this volume). The three uppermost reflectors ('o', 'p' and 'q') correlate to the package of thick sequences 1-3, and several deeper reflectors can also be correlated to sequence boundaries. The package of thick Sequences 1-3 shows a sheet-like cross-sectional geometry on seismic reflection lines, unlike the similar package recognised in CRP-2/2A.

Sedimentology and geochemistry of Oxygen Isotope Stage 3 and Holocene sediments from Laguna Potrok Aike, Patagonia

Seismic reflection studies in the maar lake Laguna Potrok Aike (51°58? S, 70°23? W) revealed an erosional unconformity associated with a sub-aquatic lake-level terrace at a water depth of 30m. Radiocarbon-dated, multi-proxy sediment studies of a piston core from this location indicate that the sediment below this discontinuity has an age of 45kyr BP (Oxygen Isotope Stage 3), and was deposited during an interval of high lake level. In comparison to the Holocene section, geochemical indicators of this older part of the record either point towards a different sediment source or to a different transport mechanism for Oxygen Isotope Stage 3 sediments. Holocene sedimentation started again before 6790cal. yr BP, providing a sediment record of hydrological variability until the present. Geochemical and isotopic data indicate a fluctuating lake level until 5310cal. yr BP. During the late Holocene the lake level shows a receding tendency. Nevertheless, the lake level did not drop below the 30m terrace to create another unconformity. The geochemical characterization of volcanic ashes reveals evidence for previously unknown explosive activity of the Reclús and Mt. Burney volcanoes during Oxygen Isotope Stage 3.

Sediment wave layer summary in ODP Site 172-1062

Abyssal mud waves (or fine-grained sediment waves) are often cited as evidence for deep current activity because subbottom profiles show that the wave form has migrated with time. The migration history of a fine-grained sediment wave on the Blake-Bahama Outer Ridge (ODP Site 1062) has been studied through the analysis of multiple ODP holes spaced across the wave. Additional information about wave migration patterns comes from 3.5-kHz records and watergun seismic profiles. These data suggest that wave migration has varied during the last not, vert, similar ~10 Myr, although the only sediments sampled are younger than 4.8 Ma. Seismic profiles suggest wave migration was initiated about 8-10 Ma, and wave migration was pronounced from about 5 Ma to about 1 Ma (with an episode of wave reorganization about 4.5 Ma). Analysis of ODP cores suggests that migration rates have been somewhat lower and more variable during the last 1 Myr. Intervals of no wave migration are observed for several time intervals and appear to characterize deglaciations, especially during the last 500 kyr. Comparisons between seismic profiles and the core record show that most of the seismic horizons correlate closely with time horizons, and thus that the seismic profiles give a reasonable representation of sediment wave migration. Models suggest that wave migration is more pronounced during periods of higher bottom current flow and less pronounced during periods of lower current flow. Thus the migration record is consistent with generally higher bottom flow speeds at this site prior to 1 Ma and lower bottom flow speeds after 1 Ma. The Mid-Pleistocene Transition from a dominant climatic periodicity of 40 kyr to a dominant climatic periodicity of 100 kyr starts at about this time, suggesting an overall reduction in bottom flow speed at this site coincident with changing climate patterns. These changes in flow speed could be related to changes in the depth of the Western Boundary Undercurrent as well as to changes in the speed of thermohaline circulation.

Average dip directions from intervals of dipmeter log in drill hole CRP-2A

Bedding dips in the CRP-2A drillhole were determined in two ways (1) analysis of a dipmeter log, and (2) identification of bed boundaries on digital images of the outer core surface. The two methods document the downhole increase in structural dip, to a maximum of 15° in the lowest 150 m of the hole. Dipmeter data, which are azimuthally oriented, indicate a 75° azimuth for structural tilting, in agreement with seismic reflection profiles. Core and log dips indicate that structural dip increases by 5-7° between 325 and 480 mbsf. Both, however, also exhibit high dip inhomogeneity because of depositional (e.g., cross bedding) and post-depositional (e.g., softsediment deformation) processes. This variability adds ambiguity to the search for angular unconformities within the CRP-2A drillhole. Dip directions of different lithologies are generally similar, as are dip directions for the four kinds of systems tracts. Downdip azimuths of sands and muds are slightly different from those of diamicts, possibly reflecting the divergence between ENE offshore dip and ESE glacial advance.

Neogene and Paleocene-Maestrichtian calcareous nannofossils in sediments of the upper continental rise off New Jersey

Maestrichtian to Holocene calcareous nannofossils from two closely spaced sites on the upper continental rise some 100 miles (161 km) southeast of Atlantic City, New Jersey, were zoned in order to help date a major canyon-cutting event in the late Miocene and to delineate and correlate other hiatuses with seismic stratigraphy. Mid-middle Eocene through middle Miocene sediments (Zones CP14 to CN6) were not recovered in these holes, but nearly all other zones are accounted for. The Eocene section is described in a companion chapter (Applegate and Wise, 1987, doi:10.2973/dsdp.proc.93.118.1987). Nannofossils are generally sparse and moderately preserved in the clastic sediments of Site 604. Sedimentation rates are extremely high for the upper Pleistocene (201 m/m.y. minimum) above a hiatus calculated to span 0.44 to 1.1 Ma. The associated disconformity is correlated with local seismic reflection Horizon Pr . Sedimentation rates continue to be high (93 m/m.y.) down to a second hiatus in the upper Pliocene dated from about 2.4 to 2.9 (or possibly 3.3) Ma. The disconformity associated with this hiatus is correlated with local seismic reflection Horizon P2 and regional Reflector Blue, which can be interpreted to mark either the onset of Northern Hemisphere continental glaciation or circulation changes associated with the closure of the Central American Seaway. Sedimentation rates in the pre-glacial lower Pliocene are only about a third those in the glacial upper Pliocene. A prominent disconformity in the upper Miocene marks a major lithologic boundary that separates Messinian(?) glauconitic claystones above from lower Tortonian conglomeratic debris flows and turbidites below. The debris flows recovered are assigned to nannofossil Zones CN8a and CN7, but drilling difficulties prevented penetration of the bottom of this sequence some 100 m below the terminal depth of the hole. Correlation of the lower bounding seismic reflector (M2/Merlin?) to a drift sequence drilled on the lower rise at DSDP Site 603, however, predicts that the debris flows began close to the beginning of the late Miocene (upper Zone CN6 time) at about 10.5 Ma. The debris flows represent a major canyon-cutting event that we correlate with the beginning of the particularly severe late Miocene glaciations believed to be associated with the formation of the West Antarctic Ice Sheet. The existence of these spectacular debris flows strongly suggest that the late Miocene glacio-eustatic low stand occurred during Vail Cycle TM3.1 (lower Tortonian) rather than during Vail Cycle TM3.2 (Messinian) as originally published. Beneath a set of coalesced regional disconformities centered upon seismic reflection Horizon Au, coccoliths are abundant and in general are moderately preserved at Site 605 in a 619-m carbonate section extending from the middle Eocene Zone CP13b to the upper Maestrichtian Lithraphidites quadratus Zone. Sedimentation rates are 37 m/m.y. in the Eocene down to a condensed interval near the base (Zone CP9). A disconformity is suspected near the Eocene/Paleocene boundary. Sedimentation rates for the upper Paleocene Zone CP8 are similar to those of the Eocene, but Zones CP7 and CP6 lie within another condensed interval. The highest Paleocene rates are 67 m/m.y. down through Zones CP5 and CP4 to a major disconformity that separates the upper Paleocene from the Danian. This hiatus spans about 2.6 m.y. (upper Zone CP3 to lower Zone CP2) and corresponds to the major sea-level drop at the base of Vail Cycle TE2.1. As the most prominent break in this Paleogene section, it may correspond to seismic reflection Horizon A* of the North American Basin. Sedimentation rates from this point to the Cretaceous/Tertiary boundary drop to 11 m/m.y., still high for a Paleocene DSDP section. No major break in deposition could be detected at the Cretaceous/Tertiary boundary.

Chronostratigraphy of sediments from the Bahamas

Continuous cores drilled during the Bahamas Drilling Project (BDP) and the Ocean Drilling Program (ODP) Leg 166 along a transect from the top of Great Bahama Bank to the basin in the Straits of Florida provide a unique data set to test the assumption in seismic stratigraphy that seismic reflections are time lines and, thus, have a chronostratigraphic significance. Seismic reflections that are identified as seismic sequence boundaries (SSBs) were dated by means of biostratigraphy in the five ODP sites and by a combination of biostratigraphy, magnetostratigraphy and Sr isotope stratigraphy in the two BDP sites. The seismic reflection horizons are carried across a variety of facies belts from shallow-water carbonates over slope carbonates to drift deposits in the Straits of Florida. Within this system 17 SSBs were identified and dated. Despite the fact that the seismic reflections cross several facies belts, their ages remain remarkably constant. The average offset in all sites is 0.38 Myr. In no cases do the seismic reflections cut across time lines. The age differences are the combined result of the biostratigraphic sampling frequency, the spacing of marker species that required extrapolation of ages, and the resolution of the seismic data. In addition, uncertainties of age determination in the proximal sites where age-diagnostic fauna are rare add to the age differences between sites. Therefore, it can be concluded that the seismic reflections, which mark the SSBs along the Bahamas Transect, are time lines and can be used as stratigraphic markers. This finding implies that depositional surfaces are preferentially imaged by reflected seismic waves and that an impedance contrast exists across these surfaces. Facies successions across the sequence boundaries indicate that the sequence boundaries coincide with the change of deposition from times of high to low sea level. In the carbonate setting of Great Bahama Bank, sea-level changes produce changes in sediment composition, sedimentation rate and diagenesis from the platform top to the basin. The combination of these factors generates differences in sonic velocity and, thus, in impedance that cause the seismic reflection. The impedance contrasts decrease from the proximal to the distal sites, which is reflected in the seismic data by a decrease of the seismic amplitude in the basinal area.

Wet and dry resistivities along three orthogonal directions of ODP Holes 185-801C and 185-1149B

The results of shore-based three-axis resistivity and X-ray computed tomography (CT) measurements on cube-shaped samples recovered during Leg 185 are presented along with moisture and density, P-wave velocity, resistivity, and X-ray CT measurements on whole-round samples of representative lithologies from Site 1149. These measurements augment the standard suite of physical properties obtained during Leg 185 from the cube samples and samples obtained adjacent to the cut cubes. Both shipboard and shore-based measurements of physical properties provide information that assists in characterizing lithologic units, correlating cored material with downhole logging data, understanding the nature of consolidation, and interpreting seismic reflection profiles.

Geochemistry of evaporite in the Red Sea

One of the major shipboard findings during Leg 23 drilling in the Red Sea was the presence of late Miocene evaporites at Sites 225, 227, and 228. The top of the evaporite sequence correlates with a strong reflector (Reflector S) which has been mapped over much of the Red Sea (Ross et al., 1969, Phillips and Ross, 1970). This indicates that the Red Sea appears to be extent. Miocene sediments, including evaporites, are known from a few outcrops along the coastal plains of the Gulf of Suez to lat 14°N (Sadek, 1959, cited in Friedman, 1972; Heybroek, 1965; Friedman, 1972). Along the length of the Red Sea, the presence of Miocene salt is indicated by seismic reflection studies (Lowell and Genik, 1972) and confirmed by drilling. The recently published data from deep exploratory wells (Ahmed, 1972) demonstrate the great thickness of elastics and evaporites which were deposited in the Red Sea depression during Miocene time. The Red Sea evaporites are of the same age as the evaporites found by deep sea drilling (DSDP Leg 13) in the Mediterranean Sea. Therefore, Reflector S in the Red Sea is comparable to Reflector M in the Mediterranean. It is assumed that during Miocene time a connection between these two basins was established (Coleman, this volume) resulting in a similar origin for the evaporites deposited in the Red Sea and in the Mediterranean Sea. The origin of the Mediterranean evaporites has been discussed in great detail (Hsü et al., 1973; Nesteroff, 1973; Friedman, 1973). The formation of evaporites may be interpreted by three different hypotheses. 1) Evaporation of a shallow restricted shelf sea or lagoon which receives inflows from the open ocean. 2) Evaporation of a deep-water basin which is separated from the open ocean by a shallow sill (Schmalz, 1969). 3) Evaporation of playas or salt lakes which are situated in desiccated deep basins isolated from the open ocean (Hsü et al., 1973). The purpose of this study is to show whether one of these models might apply to the formation and deposition of the Red Sea evaporites. Therefore, a detailed petrographic and geochemical investigation was carried out.