Digital rocks: linking forward modelling to carbonate facies

Forward stratigraphic models usually display sediment types on simulated stratigraphic profiles as ‘facies’ defined only by their depth of deposition. More recently, ‘facies’ have been defined and displayed in terms of the dominant processes of deposition (e.g. in situ growth, pelagic production, turbidite deposition). Standard carbonate facies; that is, the Dunham classification, are defined by rock textures and grain composition that imply that a combination of processes acted together to generate a facies. For example, a bioclastic wackestone is a matrix‐supported rock containing up to 90% matrix and > 10% shelly grains. In terms of modelled processes, the muddy matrix could be generated by: (i) reworking of the shallow platform sediments, (ii) from pelagic deposition, or (iii) in situ production. A traditional depth of deposition process display would not be able to distinguish such a wackestone from any other facies deposited at this water depth and a majority process display would not combine reworked, pelagic muds and in situ contribution in one simulated ‘facies’. This paper introduces a new scheme that enables forward models to output simulated facies defined by a range of values for each of the controlling processes and thereby predicts rock textures within simulated stratigraphies. This approach has been applied to the Jurassic carbonate ramps of the Iberian Basin in northeastern Spain. It is shown to provide more accurate information about the processes that are being simulated, allowing more direct comparisons to be made with the facies observed in the field and providing potential for a more rigorous method for assessing the ‘goodness of fit’ of a simulated stratigraphy.     

Icehouse world sea-level behaviour and resulting stratal patterns in late Visean (Mississippian) carbonate platforms: integration of numerical forward modelling and outcrop studies

Late Visean (Asbian–Brigantian) platform carbonates in the British Isles show a pronounced cyclicity marked by the alternation of mainly subtidal carbonates and subaerial exposure surfaces. Whereas the cyclicity of shallow‐water limestones of this age has been well‐documented, there has been little attempt to understand the controls on larger‐scale patterns such as those recognized in Pennsylvanian successions in the U.S.A. The principal aim of this study is to test two contrasting theories of cycle stacking via numerical forward modelling. Earlier studies of Pennsylvanian‐early Permian platform carbonates in south‐west U.S.A. suggest that cycle stacking patterns were controlled by the interaction of third‐ and fourth‐order sea‐level oscillations, with relatively uniform fourth‐order oscillations altered mainly by the harmonic effect of lower‐order sea‐level changes. An alternative model is based on an insolation curve for the Carboniferous calculated using Milankovitch parameters. This model predicts a considerable variability in levels of solar insolation that would have affected the amplitude of fourth‐order sea‐level changes and cycle composition. Both of these ideas were examined via numerous model runs using CARBONATE 6.0 and a new program, CARBOSMUT. Model results were evaluated through the use of 3 key criteria derived from well‐documented outcropping stratigraphies in the U.K.: (1) cycle stacking patterns and the stratigraphic position of major transgressions, (2) stratigraphic position of major faunal changes, (3) degree of development of subaerial exposure surfaces. Computer simulations and comparison with outcrop data suggest that a model invoking the interaction of relatively uniform fourth‐ (c.100 Kyr) and third‐order sea‐level oscillations is most appropriate for much of the late Visean, with major lowstands occurring at the mid‐Asbian and Asbian–Brigantian boundary. Late Visean cycles are important exploration targets in the Pri‐Caspian Basin, Kazakhstan and understanding the controls on stratal patterns is important as a potential exploration tool.     

Stream-terrace genesis: implications for soil development

Genesis of three distinct types of stream terraces can be understood through application of the concepts of tectonically induced downcutting, base level of erosion, complex response, threshold of critical power, diachronous and synchronous response times, and static and dynamic equilibrium. Climatic and tectonic stream terraces are major terraces below which flights of minor complex-response degradation terraces can form.These three types of terraces can be summarized by describing a downcutting-aggradation-renewed downcutting sequence for streams with gravell bedload. By tectonically induced downcutting, streams degrade to achieve and maintain a dynamic equilibrium longitudinal profile at the base level of erosion. Lateral erosion bevels bedrock beneath active channels to create major straths that are the fundamental tectonic stream-terrace landform. Aggradation events record brief reversals of long-term tectonically induced downcutting because they raise active channels. They may be considered as major (the result of climatic perturbations) or minor (the result of complex-response model types of perturbations). Climatically controlled aggradation followed by degradation leaves an aggradation surface; this type of fill-terrace tread is the fundamental climatic stream-terrace landform. Aggradation surfaces may be buried by subsequent episodes of deposition unless intervening tectonically induced downcutting is sufficient for younger aggradation surfaces to form below older surfaces. Raising of the active channel by either tectonic uplift or by climatically induced aggradation provides the vertical space for degradation terraces to form; first in alluvial fill and then in underlying bedrock along tectonically active streams. These are complex-response terraces because they result from interactions of dependent variables within a given fluvial system. Pauses in degradation to a new base level of erosion, and/or minor episodes of backfilling, lead to formation of complex-response fill-cut and strath, or of fill terraces. Fill-cut terraces are formed in alluvium; they are complex-response terraces because they are higher than the base level of erosion. Good exposures and dating are needed to distinguish static equilibrium complex-response minor strath terraces from dynamic equilibrium tectonic (major) straths. Strath terraces may be regarded as complex-response terraces where degradation rates between times terrace-tread formation exceed the long-term uplift rate for the reach based on ages and positions of tectonic terraces.Late Quaternary global climatic changes control aggradation events and even the times of cutting of major (tectonic) straths, because the base level of erosion can not be attained during times of climatically driven aggradation-degradation events.Most terrace soils form on treads of climatic and complex-response terraces. Aggradation surfaces may provide an ideal flight of terraces on which to study a soils chronosequence. Each aggradation event is recorded by a single relict soil where tectonically induced downcutting is sufficient to provide clear altitudinal separation of the terrace treads. Multiple paleosols are typical of tectonically stable regions where younger aggradation events spread alluvium over treads of older climatic terraces. Pedons on a climatic terrace in a small fluvial system commonly are roughly synchronous - variations of soil properties that can be attributed to temporal differences will be minor compared to altitudinally controlled climatic factors. Climatic terraces of adjacent watersheds also should be roughly synchronous (correlatable) - variations of soil properties that can be attributed to temporal differences will be minor compared to lithologic and climatic factors between different watersheds. Such generalizations may not apply to basins with sufficient relief that geomorphic responses to climatic changes occur at different and overlapping times, and to large rivers whose widely separated reaches are characterized by different response times to climatic perturbations. Soils on climatic terraces of distant watershedswill not be synchronous if their respective aggradation events occur during full-glacial times and interglacial times. Soils on some complex-response terraces may be diachronous within a given fluvial system, and typically are diachronous between watersheds.     

Diatom species variation between lake habitats: implications for interpretation of paleolimnological records

A sample of a sediment record contains diatom species that have grown in disparate habitats and eventually accumulated in a deep part of the lake. The original habitats may differ in substrate, depth location, and availability of resources. Identifying the species characteristic of each habitat should improve our ecological and environmental interpretation of the sediment record by distinguishing habitat specific responses. With this aim, we studied the benthic diatom communities of a deep oligotrophic lake across several habitats. The main source of variation in the diatom composition was the substrate type; particularly, sediment biofilms. Depth was the second factor. The thermocline defined a shift in diatom communities that also included changes in the dominant lifeforms. A third factor was the mesoscale heterogeneity (i.e., rock sides). Although most species were present in many habitats, characteristic species were identified for all the main habitats and used for an improved interpretation of the deep sediment record. Appropriate standardization showed increasing species richness and diversity from epilimnetic epilithic samples to hypolimnetic sediment samples. We estimate that more than 5000 valve counts are required for appropriate comparisons. Consequently, in sediment records with lower counts per sample, one has to amalgamate samples—losing temporal resolution—to achieve reliable analyses of diversity changes over time. Deep sediment samples are representative of the gamma-diversity of the lake diatom metacommunity, which result from the local alpha diversity of the habitats and the beta-diversity of the variability in composition among them. This double source of diversity has to be taken into account when using the sediment record for estimating lake biodiversity changes. On the other hand, we show that an estimation of the spatial (habitat) heterogeneity of a reconstructed environmental variable can be achieved using subsets of species characteristic of each habitat. We demonstrate the procedure by reconstructing the pH fluctuations during the last 200 years in several habitats from a single sediment record. The results are coherent with the expected differences between predominantly trophogenic or tropholithic habitats.     

Influence of point-source sediment-supply on modern shelf-slope morphology: implications for interpretation of ancient shelf margins

Present sea-floor bathymetry indicates that the continental-shelf and shelf-break morphology have some unique and predictable characteristics in areas with and without high sediment supply. Using a global bathymetry dataset in open shelf areas in front of rivers that discharge over 25 × 10⁶ tons of sediment per year, five distinct accretionary types of shelf-break are distinguished based on along-shelf gradient variability and inferred shelf-break trajectory. Morphological characteristics of river-mouth shelves (compared with adjacent areas lateral to the immediate fairway of the river) are: (1) an overall lower gradient and greater width, and (2) a relatively high slope gradient/shelf gradient ratio. The exceptions are shelves with active shelf-edge deltas; these are narrower, steeper and have an attenuated shelf break in front of rivers. These observations are at seismic scale and have direct implications for the recognition and positioning of principal cross-shelf, supply fairways on ancient shelves or shelf margins, and therefore the potential by-pass routes for deepwater sands. Higher slope/shelf gradient ratios in areas of actively accreting margins, where the shelf-break is more prominent and easier to recognize on seismic data compared with adjacent areas, predict areas with high sediment supply. Along-strike morphological changes on supply-dominated shelves suggest that identification of the sediment-feed route and depocenter relative to the shelf break during a relative sea level cycle are critical for understanding/predicting the 3-D architecture of the shelf-slope-basin floor clinoform.     

3D Geological modelling and gravity inversion of a structurally complex carbonate area: application for karstified massif localization

This paper proposes a new methodology to improve the location of potential karstified areas by gravity inversion of a 3D geological model. A geological 3D model is built from surface observations, 2D seismic reflection profiles and well data. The reliability of this geological 3D model obtained from integration, interpretation and interpolation of such data is first tested against the structural consistency of the model. Its theoretical gravimetric response is compared to gravity field during the forward problem in order to evaluate the validity/robustness of the geological model. The coherency between the gravity field and the gravimetric response is tested. The litho‐inversion modelling quantifies the distribution of rock density in a probabilistic way, taking into account the geology and physical properties of rocks, while respecting the geological structures represented in the 3D model. The result of the inversion process provides a density distribution within carbonate formations that can be discussed in term of karstification distribution. Thus, lower densities correlate with areas that are strongly karstified. Conversely, higher than mean densities are found in carbonate formations mostly located under marly and impervious formations, preserving carbonate from karstification and paleokarstification.     

Numerical modelling of wavefields in three-dimensional inhomogeneous media

Summary. Numerical modelling is one of the most efficient methods for an investigation of the relationship between structural features and peculiarities of observed wavefields. It is practically the only method for 2-D and 3-D inhomogeneous media.
An algorithm based on ray theory has been developed for calculations of travel times and amplitudes of seismic waves in 3-D inhomogeneous media with curved interfaces. It was applied for numerical modelling of kinematic and dynamic characteristics of seismic waves propagating in laterally inhomogeneous media.
Travel-time and amplitude patterns were studied in the 2-D and 3-D models of a geosyncline, in which velocity distribution was given by an analytical function of the coordinates. For a more complicated model representing a subducting high-velocity lithospheric plate in a transition zone between oceanic and continental upper mantle, the velocity distribution was given by discrete values on a 2-D non-rectangular grid. It was shown that when a source was placed above the lithospheric plate, a shadow zone appeared along a strike of the structure, i.e. in the direction which is perpendicular to a strong lateral velocity gradient. Travel-time residuals were calculated along the seismological profile for a 3-D velocity distribution in the upper mantle beneath Central Asia, obtained as a result of inversion of travel times by the Backus-Gilbert method. They were found to be in a good agreement with the observed data.     

Variation among fossil chironomid assemblages in surficial sediments of Bodensee-Untersee (SW-Germany): implications for paleolimnological interpretation

In paleolimnology, subfossil head capsules of chironomids play an important role as ecological indicators of lake history. It is important to determine, therefore, whether fossil assemblages are representative of former biocoenoses. There is evidence that headcapsules washed in from other places can make up a significant percentage of the total. As interpretations are usually drawn from the examination of a single core, it is of special interest to know whether a fossil assemblage of a single site properly reflects limnological conditions of the whole lake. This study examined the taxonomic distribution of subfossil chironomids in the surficial sediments of the Bodensee-Untersee, with the aim of assessing the variability in chironomid assemblages. Apparently, most of the head capsules of the profundal fossil assemblages in the Untersee had been washed in from the littoral zone or from the slope. Although the Bodensee-Untersee is a rather large lake, variability is surprisingly low among all samples. Therefore a correct interpretation from a single core may be possible.     

黄土碳酸盐碳同位素环境意义讨论

黄土碳酸盐碳同位素广泛应用于第四纪气候环境变化的研究中,以往的大多数研究中无论是利用钙结核、次生碳酸盐还是成壤碳酸盐,认为其反映了C4植物的丰度。黄土高原碳酸盐碳同位素表现为黄土层高,古土壤层中低,即黄土层中C4植物丰度高于古土壤层。然而,这样的结果和黄土有机碳同位素得到的结果矛盾,有机碳同位素的结果表明温度对C4植物的分布起到了决定性作用。由于有机碳同位素对植物类型的反映更为直接而可靠,因此碳酸盐碳同位素反映C4植物丰度存在疑问。对黄土高原黄土碳酸盐碳同位素的系统概括后认为,第四纪期间黄土碳酸盐碳同位素与C4植物有直接联系,但C4植物丰度不是唯一决定性的因素,碳酸盐碳同位素的指示意义存在复杂性。在黄土高原地区,植被发育程度、与大气CO2交换程度、植被本身的碳同位素值的变化以及原生碳酸盐的影响等因素都会对碳酸盐碳同位素产生影响。由黄土碳酸盐碳同位素的讨论可延伸到不同土壤碳酸盐碳同位素揭示的环境指示意义,不同的土壤环境,其气候条件、植被类型及发育程度、大气CO2的交换情况、微生物的活动及土壤次生碳酸盐受原生碳酸盐溶解的影响等因素都会对碳同位素产生不同程度影响,哪种或哪几种因素产生主要作用,在不同区域土壤环境中是不一样的。具体研究中需确定影响的核心因素,才能确定碳同位素的环境指示意义。     

Water and sediment chemistry of Lake Pumayum Co,South Tibet,China: implications for interpreting sediment carbonate

A combination of water and sediment chemistry was used to investigate carbonate production and preservation in Lake Pumayum Co (altitude 5,030 m a.s.l.), south Tibet, China. We compared the chemical composition of lake water in various parts of the lake with that of input rivers and found that the loss of Ca²⁺ results from calcite sedimentation induced by evaporation and biogenic precipitation. This is supported by evaporation data from the catchment and δ¹⁸O measurements on water. Results suggest that CaCO₃ is the predominant carbonate in this lake. There is a positive correlation in the sediments among concentrations of total inorganic carbon (TIC), Ca, total organic carbon (TOC), and total nitrogen, confirming that most carbonates in sediment are endogenic. The Jiaqu River is the largest inflow to Lake Pumayum Co and has a strong influence on both lake water chemistry and sediment composition. The river and lake bathymetry influence carbonate sedimentation by affecting water flow velocity and growing conditions for macrophytes. Different carbon contents and relationships between TIC and TOC in the two long cores from different depths in the lake reveal that hypolimnetic conditions also influence carbonate precipitation and preservation.