High-resolution sedimentary record of lead in the coastal mud of the East China Sea in the past one hundred and fifty years

In the East China Sea （ECS）, there are some mud areas, including the south coastal mud area, the north coastal mud area, and the mud area to the southwest of Cheju Island. X-ray fluorescence （XRF） techniques and Thermal Ionization Mass Spectrometry （TIMS） were used to study the high-resolution sedimentary record of Pb concentrations and Pb stable isotopic compositions in the past one hundred and fifty years in the coastal mud of the ECS. Pb concentrations of a ^210Pb dating S5 core in the study area have increased rapidly since 1980, and reached the maximal value with 65.08 μg/g in 2000, corresponding to the fast economic development of China since the implementation of the ＂Reform and Open Policy＂ in 1978; ^206Pb/^207Pb ratios generally had stabilized at 1.195 from 1860 to 1966, and decreased gradually from 1966 to 2000, indicating that the anthropogenic source Pb contribution to the ECS has increased gradually since 1966, especially since 1980. Pb concentrations decreased distinctly from 2000 to 2003 and ^206Pb/^207Pb ratios increased from 2001 to 2003, corresponding closely to the ban of lead gasoline from 2000 in China. From 1950 to 2003, there occurred four distinct decrease events of ^206Pb/^207Pb, possibly responding to the Changjiang River （Yangtze River） catastrophic floods in 1998, 1991, 1981 and 1954; from 1860 to 1966, there were two decrease periods of ^206Pb/^207Pb, which may respond to the catastrophic floods of Changjiang River in 1931 and 1935, and 1870. As a result of the erosion and drowning by the catastrophic floods, the anthropogenic lead accumulated in soil and water environments over a long period of time was brought into the Changjiang River, then part of them was finally transported into the ECS, which leads to changes in Pb stable isotopic composition.     

Temperature–composition–time (T–X–t) data from authigenic K-feldspar: An integrated methodology for dating fluid flow events

Integration of fluid inclusion analysis with high spatial resolution Ar–Ar dating of K-feldspar cements has been used to resolve and reconstruct palaeo-fluid flow. Fluid inclusion analysis allows discrimination of distinct cement phases, thereby identifying discrete episodes of fluid flow. Ar–Ar dating of the same cements via high spatial resolution laserprobe establishes absolute age constraints on the framework previously constructed. Integration of these two datasets yields temperature–composition–time data.     

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.     

Temporal behavior of the background seismicity rate in central Japan, 1998 to mid-2003

We studied the temporal behavior of the background shallow seismicity rate in 700 circular areas across inland Japan. To search for and test the significance of the possible rate changes in background seismicity, we developed an efficient computational method that applies the space–time ETAS model proposed by Ogata in 1998 to the areas. Also, we conducted Monte Carlo tests using a simulated catalog to validate the model we applied. Our first finding was that the activation anomalies were found so frequently that the constant background seismicity hypothesis may not be appropriate and/or the triggered event model with constraints on the parameters may not adequately describe the observed seismicity. However, quiescence occasionally occurs merely by chance. Another outcome of our study was that we could automatically find several anomalous background seismicity rate changes associated with the occurrence of large earthquakes. Very significant seismic activation was found before the M6.1 Mt. Iwate earthquake of 1998. Also, possible seismic quiescence was found in an area 150 km southwest of the focal region of the M7.3 Western Tottori earthquake of 2000. The seismicity rate in the area recovered after the mainshock.     

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.