Carbonates calibrated against oceanographic parameters along a latitudinal transect in the Gulf of California, Mexico

Trophic resources are an important control governing carbonate production. Though this importance has long been recognized, no calibration exists to quantitatively compare biogenic assemblages within trophic resource fields. This study presents a field calibration of carbonate producers in a range of settings against high‐resolution in situ measurements of nutrients, temperature and salinity. With its latitudinal extent from 30° to 23° N, the Gulf of California, Mexico, spans the warm‐temperate realm and encompasses nutrient regimes from oligo‐mesotrophic in the south to eutrophic in the north. Accordingly, from south to north carbonates are characterized by: (i) coral‐dominated shallow carbonate factories (5–20 m water depth) with average sea‐surface temperatures of 25 °C (min. 18 °C, max. 31 °C), average salinities of 35·06‰ and average chlorophyll a levels, which are a proxy for nutrients, of 0·25 mg Chl a m^?3 (max. 0·48, min. 0·1). (ii) Red algal‐dominated subtidal to inner‐shelf carbonate formation (10–25 m) in the central Gulf of California exhibiting average temperatures of 23 °C (min. 18 °C, max. 30 °C), average salinities of 35·25‰, and average Chl a levels of 0·71 Chl a m^?3 (max. 5·62, min. 0). (iii) Molluskan bryozoan‐rich inner to outer shelf factories in the northern Gulf of California (20–50 m) with average sea surface temperatures of only 20 °C (min. 13 °C, max 29 °C), average salinities of 35·01‰, and average contents of 2·2 mg Chl a m^?3 (max. 8·38, min. 0). By calibrating sedimentological data with in situ measured oceanographic information in different environments, the response of carbonate producers to environmental parameters was established and extrapolated to carbonates on a global scale. The results demonstrate the importance of recognizing and quantifying trophic resources as a dominant control determining the biogenic composition and facies character of both modern and fossil carbonates.     

A sequence stratigraphic model for an intensely bioturbated shallow-marine sandstone: the Bridport Sand Formation, Wessex Basin, UK

The Bridport Sand Formation is an intensely bioturbated sandstone that represents part of a mixed siliciclastic‐carbonate shallow‐marine depositional system. At outcrop and in subsurface cores, conventional facies analysis was combined with ichnofabric analysis to identify facies successions bounded by a hierarchy of key stratigraphic surfaces. The geometry of these surfaces and the lateral relationships between the facies successions that they bound have been constrained locally using 3D seismic data. Facies analysis suggests that the Bridport Sand Formation represents progradation of a low‐energy, siliciclastic shoreface dominated by storm‐event beds reworked by bioturbation. The shoreface sandstones form the upper part of a thick (up to 200 m), steep (2–3°), mud‐dominated slope that extends into the underlying Down Cliff Clay. Clinoform surfaces representing the shoreface‐slope system are grouped into progradational sets. Each set contains clinoform surfaces arranged in a downstepping, offlapping manner that indicates forced‐regressive progradation, which was punctuated by flooding surfaces that are expressed in core and well‐log data. In proximal locations, progradational shoreface sandstones (corresponding to a clinoform set) are truncated by conglomerate lags containing clasts of bored, reworked shoreface sandstones, which are interpreted as marking sequence boundaries. In medial locations, progradational clinoform sets are overlain across an erosion surface by thin (<5 m) bioclastic limestones that record siliciclastic‐sediment starvation during transgression. Near the basin margins, these limestones are locally thick (>10 m) and overlie conglomerate lags at sequence boundaries. Sequence boundaries are thus interpreted as being amalgamated with overlying transgressive surfaces, to form composite erosion surfaces. In distal locations, oolitic ironstones that formed under conditions of extended physical reworking overlie composite sequence boundaries and transgressive surfaces. Over most of the Wessex Basin, clinoform sets (corresponding to high‐frequency sequences) are laterally offset, thus defining a low‐frequency sequence architecture characterized by high net siliciclastic sediment input and low net accommodation. Aggradational stacking of high‐frequency sequences occurs in fault‐bounded depocentres which had higher rates of localized tectonic subsidence.     

Hydrology and hydrogeology of Sakaryabaşı Karstic springs, Çifteler, Turkey

The Sakarya River is one of the largest rivers in Turkey and is fed mainly from Sakaryabaşı springs. The Sakaryabaşı springs are located in the Central Anatolia and issue from confined/semi-confined karst having a thermal component and therefore, having quite different hydrogeological characteristics as compared to the Taurus Karst region, a typical example of the Mediterranean type of karst. The karstic carbonate rocks that form the groundwater reservoir are overlain by a thick semi-pervious overburden of mainly clastics of Neogene age. Tectonics is the major factor controlling the occurrence of the karst springs in the area where topography is rather flat. This study aimed at explaining the occurrence and movement of the karst groundwater within the system by use of hydrogeological, chemical, and isotopic tools. Isotopic composition of the waters revealed that all waters in the region are of meteoric origin and the thermal component is due to deep circulation. The catchment area of the hydrogeological system extends to the south and groundwater movement is towards the outlets, which are in a depression along a major fault. The movement of the groundwater, based on analysis of remotely sensed images, is controlled mainly by structural elements.     

Correlations in aftershock and seismicity patterns

Correlations in space and time play a fundamental role in earthquake processes. One direct manifestation of the effects of correlations is the occurrence of aftershocks due to the stress transfer in the vicinity of a main shock. Less obvious and more speculative changes in correlations may occur in the background seismicity before large earthquakes. Using statistical physics it is possible to introduce a measure of spatial correlations through a correlation length. This quantity characterizes how local fluctuations can influence the occurrence of earthquakes over distances comparable with the correlation length. In this work, the physical basis of spatial correlations of earthquakes is discussed in the context of critical phenomena and the percolation problem. The method of two-point correlation function is applied to the seismicity of California. Well defined variations in time of the correlation length are found for aftershock sequences and background seismicity. The scaling properties of our obtained distributions are analyzed with respect to changes in several scaling parameters such as lower magnitude cutoff of earthquakes, the maximum time interval between earthquakes, and the spatial size of the area considered. This scaling behavior can be described in a unified manner by utilizing the multifractal fit. Utilizing the percolation approach the time evolution of clusters of earthquakes is studied with the correlation length defined in terms of the radius of gyration of clusters. This method is applied to the seismicity of California.     

Effects of magnetic interactions in anisotropy of magnetic susceptibility: Models, experiments and implications for igneous rock fabrics quantification

Besides granites of the ilmenite series, in which the anisotropy of magnetic susceptibility (AMS) is mainly controlled by paramagnetic minerals, the AMS of igneous rocks is commonly interpreted as the result of the shape-preferred orientation of unequant ferromagnetic grains. In a few instances, the anisotropy due to the distribution of ferromagnetic grains, irrespective of their shape, has also been proposed as an important AMS source. Former analytical models that consider infinite geometry of identical and uniformly magnetized and coaxial particles confirm that shape fabric may be overcome by dipolar contributions if neighboring grains are close enough to each other to magnetically interact. On these bases we present and experimentally validate a two-grain macroscopic numerical model in which each grain carries its own magnetic anisotropy, volume, orientation and location in space. Compared with analytical predictions and available experiments, our results allow to list and quantify the factors that affect the effects of magnetic interactions. In particular, we discuss the effects of (i) the infinite geometry used in the analytical models, (ii) the intrinsic shape anisotropy of the grains, (iii) the relative orientation in space of the grains, and (iv) the spatial distribution of grains with a particular focus on the inter-grain distance distribution. Using documented case studies, these findings are summarized and discussed in the framework of the generalized total AMS tensor recently introduced by Cañon-Tapia (Cañon-Tapia, E., 2001. Factors affecting the relative importance of shape and distribution anisotropy in rocks: theory and experiments. Tectonophysics, 340, 117–131.). The most important result of our work is that analytical models far overestimate the role of magnetic interaction in rock fabric quantification. Considering natural rocks as an assemblage of interacting and non-interacting grains, and that the effects of interaction are reduced by (i) the finite geometry of the interacting clusters, (ii) the relative orientation between interacting grains, (iii) their heterogeneity in orientation, shape and bulk susceptibility, and (iv) their inter-distance distribution, we reconcile analytical models and experiments with real case studies that minimize the role of magnetic interaction onto the measured AMS. Limitations of our results are discussed and guidelines are provided for the use of AMS in geological interpretation of igneous rock fabrics where magnetic interactions are likely to occur.     

Aftershocks series monitoring of the September 18, 2004 M = 4.6 earthquake at the western Pyrenees: A case of reservoir-triggered seismicity?

On September 18, 2004, a 4.6 mbLg earthquake was widely felt in the region around Pamplona, at the western Pyrenees. Preliminary locations reported an epicenter less than 20 km ESE of Pamplona and close to the Itoiz reservoir, which started impounding in January 2004. The area apparently lacks of significant seismic activity in recent times. After the main shock, which was preceded by series of foreshocks reaching magnitudes of 3.3 mbLg, a dense temporal network of 13 seismic stations was deployed there to monitor the aftershocks series and to constrain the hypocentral pattern. Aftershock determinations obtained with a double-difference algorithm define a narrow epicentral zone of less than 10 km², ESE–WNW oriented. The events are mainly concentrated between 3 and 9 km depth. Focal solutions were computed for the main event and 12 aftershocks including the highest secondary one of 3.8 mbLg. They show mainly normal faulting with some strike-slip component and one of the nodal planes oriented NW–SE and dipping to the NE. Cross-correlation techniques applied to detect and associate events with similar waveforms, provided up to 33 families relating the 67% of the 326 relocated aftershocks. Families show event clusters grouped by periods and migrating from NW to SE. Interestingly, the narrow epicentral zone inferred here is located less than 4 km away from the 111-m high Itoiz dam. These hypocentral results, and the correlation observed between fluctuations of the reservoir water level and the seismic activity, favour the explanation of this foreshock–aftershock series as a rapid response case of reservoir-triggered seismicity, burst by the first impoundment of the Itoiz reservoir. The region is folded and affected by shallow dipping thrusts, and the Itoiz reservoir is located on the hangingwall of a low angle southward verging thrust, which might be a case sensible to water level fluctuations. However, continued seismic monitoring in the coming years is mandatory in this area to infer more reliable seismotectonic and hazard assessments.     

A comparative error analysis of manual versus automated methods of data acquisition for algebraic strain estimation

This paper presents a statistical analysis of the algebraic strain estimation algorithm of Shimamoto and Ikeda [Shimamoto, T., Ikeda, Y., 1976. A simple algebraic method for strain estimation from deformed eillipsoidal objects: 1. Basic theory. Tectonophysics 36, 315–337]. It is argued that the error in their strain estimation procedure can be quantified using an expected discrepancy measure. An analysis of this measure demonstrates that the error is inversely proportional to the number of clasts used. The paper also examines the role of measurement error, in particular that incurred under (i) a moment based and (ii) manual data acquisition methods. Detailed analysis of these two acquisition methods shows that in both cases, the effect of measurement error on the expected discrepancy is small relative to the effect of the sample size (number of objects). Given their relative speed advantage, this result favours the use of automated measurement methods even if they incur more measurement error on individual objects. Validation of these results is carried out by means of a simulation study, as well as by reference to studies appearing in previous literature. The results are also applied to obtain an upper bound on the error of strain estimation for various studies published in the literature on strain analysis.     

Modèles géostatistiques de concentrations ou de débits le long des cours d'eau

Estimating concentrations or flow rates along a stream network requires specific models. Two classes of models, recently proposed in the literature, are generalized, to the intrinsic case in particular. We present a global construction by ‘streams’, i.e. on the whole set of paths between sources and outlet. Combining stationary or intrinsic one-dimensional random functions leads to stationary or intrinsic models on segments, with discontinuities at the forks. A construction from outlet to sources, leads to stationary or intrinsic models on each stream, without any discontinuity at the forks. The linear variogram is found as a particular case. The extension to the linear model of coregionalization is immediate, allowing a multivariate modelling of concentrations. To cite this article: C. de Fouquet, C. Bernard-Michel, C. R. Geoscience 338 (2006).     

Water temperature behaviour in the River Loire since 1976 and 1881

Analysis of monthly mean river temperatures, recorded on an hourly basis in the middle reaches of the Loire since 1976, allows reconstruction by multiple linear regression of the annual, spring and summer water temperatures from equivalent information on air temperatures and river discharge. Since 1881, the average annual and summer temperatures of the Loire have risen by approximately 0.8?°C, this increase accelerating since the late 1980s due to the rise in air temperature and also to lower discharge rates. In addition, the thermal regime in the Orleans to Blois reach is considerably affected by the inflow of groundwater from the Calcaires de Beauce aquifer, as shown by the summer energy balance. To cite this article: F. Moatar, J. Gailhard, C. R. Geoscience 338 (2006).     

Spatial and temporal variability of Holocene temperature in the North Atlantic region

The early Holocene climate of the North Atlantic region was influenced by two boundary conditions that were fundamentally different from the present: the presence of the decaying Laurentide Ice Sheet (LIS) and higher than present summer solar insolation. In order to assess spatial and temporal patterns of Holocene climate evolution across this region, we collated quantitative paleotemperature records at sub-millennial resolution and synthesized their temporal variability using principal components analysis (PCA). The analysis reveals considerable spatial variability, most notably in the time-transgressive expression of the Holocene thermal maximum (HTM). Most of the region, but especially areas peripheral to the Labrador Sea and hence closest to the locus of LIS disintegration, experienced maximum Holocene temperatures that lagged peak summer insolation by 1000-3000 years. Many sites from the northeastern North Atlantic sector, including the Nordic Seas and Scandinavia, either warmed in phase with maximum summer insolation (11,000-9000 years ago) or were less strongly lagged than the Baffin Bay-Labrador Sea region. These spatially complex patterns of Holocene climate development, which are defined by the PCA, resulted from the interplay between final decay of the LIS and solar insolation forcing.