Stochastic modelling of multi-grain equivalent dose (De) distributions: Implications for OSL dating of sediment mixtures

A number of recent optically stimulated luminescence (OSL) studies have cited post-depositional mixing as a dominant source of equivalent dose (D_e) scatter across a range of sedimentary environments, including those previously considered ‘best suited’ for OSL dating. The potentially insidious nature of sediment mixing means that this problem may often only be identifiable by careful statistical analysis of D_e data sets. This study aims to address some of the important issues associated with the characterisation and statistical treatment of mixed D_e distributions at the multi-grain scale of analysis, using simulated D_e data sets produced with a simple stochastic model. Using this Monte Carlo approach we were able to generate theoretical distributions of single-grain D_e values, which were then randomly mixed together to simulate multi-grain aliquot D_e distributions containing a known number of mixing components and known corresponding burial doses. A range of sensitivity tests were undertaken using sediment mixtures with different aged dose components, different numbers of mixing components, and different types of dose component distributions (fully bleached, heterogeneously bleached and significantly overdispersed D_e distributions). The results of our modelling simulations reveal the inherent problems encountered when dating mixed sedimentary samples with multi-grain D_e estimation techniques. ‘Phantom’ dose components (i.e. discrete dose populations that do not correspond to the original single-grain mixing components) are an inevitable consequence of the ‘averaging’ effects of multi-grain D_e analysis, and prevent the correct number of mixing components being identified with the finite mixture model (FMM) for all of the multi-grain mixtures tested. Our findings caution against use of the FMM for multi-grain aliquot D_e data sets, even when the aliquots consist of only a few grains.     

CO2 saturation estimates at Sleipner (North Sea) from seismic tomography and rock physics inversion

CO₂ saturations are estimated at Sleipner using a two-step imaging workflow. The workflow combines seismic tomography (full-waveform inversion) and rock physics inversion and is applied to a two-dimensional seismic line located near the injection point at Sleipner. We use baseline data (1994 vintage, before CO₂ injection) and monitor data that was acquired after 12 years of CO₂ injection (2008 vintage). P-wave velocity models are generated using the Full waveform inversion technology and then, we invert selected rock physics parameters using an rock physics inversion methodology. Full waveform inversion provides high-resolution P-wave velocity models both for baseline and monitor data. The physical relations between rock physics properties and acoustic wave velocities in the Utsira unconsolidated sandstone (reservoir formation) are defined using a dynamic rock physics model based on well-known Biot–Gassmann theories. For data prior to injection, rock frame properties (porosity, bulk and shear dry moduli) are estimated using rock physics inversion that allows deriving physically consistent properties with related uncertainty. We show that the uncertainty related to limited input data (only P-wave velocity) is not an issue because the mean values of parameters are correct. These rock frame properties are then used as a priori constraint in the monitor case. For monitor data, the Full waveform inversion results show nicely resolved thin layers of CO₂–brine saturated sandstones under intra-reservoir shale layers. The CO₂ saturation estimation is carried out by plugging an effective fluid phase in the rock physics model. Calculating the effective fluid bulk modulus of the brine–CO₂ mixture (using Brie equation in our study) is shown to be the key factor to link P-wave velocity to CO₂ saturation. The inversion tests are done with several values of Brie/patchiness exponent and show that the CO₂ saturation estimates are varying between 0.30 and 0.90 depending on the rock physics model and the location in the reservoir. The uncertainty in CO₂ saturation estimation is usually lower than 0.20. When the patchiness exponent is considered as unknown, the inversion is less constrained and we end up with values of exponent varying between 5 and 20 and up to 33 in specific reservoir areas. These estimations tend to show that the CO₂–brine mixing is between uniform and patchy mixing and variable throughout the reservoir.     

The influence of gas bubbles on sediment acoustic properties: in situ, laboratory, and theoretical results from Eckernförde Bay, Baltic sea

Acoustic turbidity caused by the presence of gas bubbles in seafloor sediments is a common occurrence worldwide,but is as yet poorly understood. The Coastal Benthic Boundary Layer experiment in the Baltic off northern Germany was planned to better characterize the acoustic response of a bubbly sediment horizon. In this context, in situ measurements of compressional wave speed and attenuation were made over the frequency range of 5–400 kHz in gassy sediments of Eckernförde Bay. Dispersion of compressional speed data was used to determine the upper limit of the frequency of methane bubble resonance at between 20 and 25 kHz. These data, combined with bubble size distributions determined from CT scans of sediments in cores retained at ambient pressure, yield estimates of effective bubble sizes of 0.3–5.0 mm equivalent radius. The highly variable spatial distribution of bubble volume and bubble size distribution is used to reconcile the otherwise contradictory frequency-dependent speed and attenuation data with theory. At acoustic frequencies above resonance (>25 kHz) compressional speed is unaffected by bubbles and scattering from bubbles dominates attenuation. At frequencies below resonance (<1 kHz) ‘compressibility effects’ dominate, speed is much lower (250 m s^-1) than bubble-free sediments, and attenuation is dominated by scattering from impedance contrasts. Between 1.5 and 25 kHz bubble resonance greatly affects speed and attenuation. Compressional speed in gassy sediments (1100–1200 m s^-1) determined at 5–15 kHz is variable and higher than predicted by theory (<250 m s^-1). These higher measured speeds result from two factors: speeds are an average of lower speeds in gassy sediments and higher speeds in bubble-free sediments; and the volume of smaller-sized bubbles which contribute to the lower observed speeds is much lower than total gas volume. The frequency-dependent acoustic propagation is further complicated as the mixture of bubble sizes selectively strips energy near bubble resonance frequencies (very high attenuation) allowing lower and higher frequency energy to propagate. It was also demonstrated that acoustic characterization of gassy sediments can be used to define bubble size distribution and fractional volume.     

Heating of the Lithosphere during Meteorite Cratering

The formation of thermal anomalies around the impact sites of large cosmic bodies on the Earth is studied. The parameters of thermal anomalies are compared for the impacts of bodies of various scales—from one to several hundred kilometers in diameter. The cooling time of the rocks under impact craters of various scales is estimated. The estimates obtained are used to model the input of heat by the impacts of small (less than 500 km in diameter) planetesimals late in the accretion of the Earth. The boundary conditions for calculating the thermal evolution of the early Earth are refined by simultaneously analyzing the sizes of impact thermal anomalies and the model size distributions of projectiles (the mass spectrum of planetesimals).     

1971-2009 年珠穆朗玛峰地区尼泊尔境内气候变化

利用珠穆朗玛峰南坡尼泊尔境内(科西河流域) 的10 个气象站1971-2009 年月平均气温、月平均最高、最低气温和逐月降水资料, 采用线性趋势、Sen 斜率估计、Mann-Kendall 等方法分析区域气候变化状况及其时空特征, 并与珠穆朗玛峰北坡地区气候进行比较, 分析了珠穆朗玛峰地区气候变化的特征与趋势。结果表明:(1) 1971-2009 年间, 珠穆朗玛峰南坡年平均气温为20.0℃, 线性升温率为0.25℃/10a, 与北坡主要受年平均最低气温影响相反, 增幅主要受年平均最高气温升高的影响, 并且在1974 年及1992 年间出现两次显著增温, 增温特别明显的月份为2 月和9 月;(2) 该地区降水变化的局地性较强, 近40 年间年平均降水量为1729.01 mm, 年平均降水量以每年约4.27 mm的线性增幅有所增加, 但并不显著, 且降水月变化和季变化特征均不明显;(3) 由于珠穆朗玛峰南坡受到季风带来暖湿气流和喜马拉雅山阻挡的双重影响, 珠峰南坡的年平均降水量远高于北坡;(4) 珠穆朗玛峰南坡气温变暖的海拔依赖性并不明显, 且南坡地区的变暖趋势并没有北坡变暖趋势明显。     

Use of Surface Waves in Statistical Correlations of Shear Wave Velocity and Penetration Resistance of Chennai Soils

Shear wave velocity (V _s) is one of the most important input parameter to represent the stiffness of the soil layers. It is preferable to measure V _s by in situ wave propagation tests, however it is often not economically feasible to perform the tests at all locations. Hence, a reliable correlation between V _s and standard penetration test blow counts (SPT-N) would be a considerable advantage. This paper presents the development of empirical correlations between V _s and SPT-N value for different categories of soil in Chennai city characterized by complex variation of soil conditions. The extensive shear wave velocity measurement was carried out using Multichannel Analysis of Surface Waves (MASW) technique at the sites where the SPT-N values are available. The bender element test is performed to compare the field MASW test results for clayey soils. The correlations between shear wave velocity and SPT-N with and without energy corrections were developed for three categories of soil: all soils, sand and clay. The proposed correlations between uncorrected and energy corrected SPT-N were compared with regression equations proposed by various other investigators and found that the developed correlations exhibit good prediction performance. The proposed uncorrected and energy corrected SPT-N relationships show a slight variation in the statistical analysis indicating that both the uncorrected and energy corrected correlations can predict shear wave velocity with equal accuracy. It is also found that the soil type has a little effect on these correlations below SPT-N value of about 10.     

Age,provenance and tectonic setting of the Canastra and Ibiá Groups (Brasília Belt,Brazil): Implications for the age of a Neoproterozoic glacial event in central Brazil

The Brasília Belt is one of the best preserved Neoproterozoic orogens in Brazil. It comprises a thick Meso–Neoproterozoic sedimentary/metasedimentary pile including the Canastra and Ibiá Groups, which are the object of this study. The Canastra Group constitutes a regressive sedimentary sequence made mainly of greenschist-facies metapelitic and metapsammitic rocks, including phyllite, sandy metarhythmite and quartzite, with minor intercalations of limestone, as well as carbonaceous and carbonatic phyllite. The Ibiá Group is formed of a basal diamictite followed upwards by phyllites and calc-schists. It rests on an erosional unconformity on top of the Canastra Group.A provenance study based on U–Pb zircon geochronology on a selection of seven samples helped to establish the various source areas and maximum depositional ages of the original sediments. In addition, seven new Sm–Nd analyses are presented and discussed together with previously published data.LAM-ICP-MS U–Pb dating of detrital zircon grains indicates a maximum depositional age of the Canastra and Ibiá Groups of ca. 1030 and 640 Ma, respectively. The provenance signature of the Canastra Group comprises a wide range of detrital zircon ages with a significant Paleoproterozoic component (～1.8 and ～2.1 Ga) and an important Mesoproterozoic source (1.1–1.2 Ga), especially for the Paracatu Formation, indicating the São Francisco–Congo Craton as main source. These provenance data, in particular the absence of Neoproterozoic zircon grains, typical of the active margin of the Brasília Belt, allied with the homogeneous Paleoproterozoic T_DM values are consistent with the previous interpretation that the Canastra Group represents a sedimentary sequence deposited on a passive margin setting.Zircon grains from the diamictite of the Ibiá Group yielded ages ranging from 936 to 2500 Ma. In contrast, the overlying calc-phyllite of the Rio Verde Formation reveals a dominant Neoproterozoic provenance pattern with important peaks at 665, 740 and 850 Ma. The São Francisco-Congo Craton and Goiás Magmatic Arc are, most probably, the two main source regions for the Ibiá Group which may represent, therefore, a former fore- or back-arc sedimentary sequence. Tectonically, therefore, the Ibiá Group is equivalent to the Araxá Group exposed in central Goiás and both represent syn-orogenic sedimentary sequences formed with important detrital contributions derived from the Neoproterozoic Goiás Arc.The provenance data presented here indicate that the Cubatão Formation is most possibly representative of a Marinoan or younger glacial event.     

Mg2+ removal in the system Mg2+—amorphous SiO2—H2O by adsorption and Mg-hydroxysilicate precipitation

Two chemical processes can remove Mg²⁺ from suspensions containing amorphous silica (am-SiO₂) at low temperatures: adsorption and precipitation of a Mg-hydroxysilicate resembling sepiolite. Mg²⁺ removal from am-SiO₂ suspensions was investigated, and the relative role of the two removal processes evaluated, as a function of: pH, ionic strength, Mg²⁺ concentration, and temperature.The extent of Mg²⁺ adsorption onto am-SiO₂ decreases with increasing NaCl concentration due to displacement of Mg²⁺ by Na⁺. At NaCl concentrations of 0.05 M and above, adsorption occurs only at pH values above 8.5, where rapid dissolution of am-SiO₂ gives rise to high concentrations of dissolved silica, resulting in supersaturation with respect to sepiolite. Removal of Mg²⁺, at concentrations of 40 to 650 μM, from am-SiO₂ suspensions in 0.70 M NaCl at 25 °C occurs at pH 9.0 and above. Experiments show that under these conditions adsorption and Mg-hydroxysilicate precipitation remove Mg²⁺ at similar rates. For 0.05 M Mg²⁺, at 0.70 M ionic strength and 25 °C, measurable Mg²⁺ removal occurs down to ca. pH 7.5 but is primarily due to Mg-hydroxysilicate precipitation. For the same solution conditions at 5°C, Mg²⁺ removal occurs above pH 8.0 and is primarily due to adsorption.Assuming that increasing pressure does not greatly enhance adsorption, Mg²⁺ adsorption onto am-SiO₂ is an insignificant process in sea water. The surface charge of pristine am-SiO₂ in sea water is primarily controlled by interactions with Na⁺. The principal reaction between Mg²⁺ and am-SiO₂ in marine sediments is sepiolite precipitation.The age distribution of sepiolite in siliceous pelagic sediments is influenced by temperatures of bottom waters and by geothermal gradients.     

Sea-level, sediment supply and coastal changes: Examples from the coast of Ireland

Sediment supply and pre-existing shoreline morphology are crucial factors in controlling coastal changes due to sea-level rise. Using examples from both southeast and northeast Ireland, it can be shown that sea-level change may trigger a sequence of events which leads to both static and dynamic shoreline equilibrium. Cliff erosion and longshore sediment movement in east Co. Wexford has led to injection of sediment onto the shelf, and the growth, under both wave and tide regimes, of linear offshore shoals. These shoals now control the pattern of shoreline erosion and provide a template for possible stepwise evolution of the coast under any future sea-level rise. In contrast, the nearby coast of south Co. Wexford comprises a series of coarse clastic barriers moving monotonously onshore, via overwash processes. Here the behavior of the barrier is conditioned by the antecedent morphology of both the beach face and stream outlet bedforms. Finally, the rock platform coast of Co. Antrim presents a far more resistant shoreline to incident marine processes, yet even here there is strong evidence of present process control over so-called ‘raised’ platforms and embayments. It is concluded that coastal sediment supply and dynamics, together with coastal morphology and its interaction with waves, present a far more complex variety of sea-level indicators than is normally acknowledged.     

Validation and Application of Plants as Biomonitors of Trace Element Atmospheric Pollution – A Co-Ordinated Effort in 14 Countries

The International Atomic Energy Agency (IAEA) organized a co-ordinated research project (CRP) on Validation and application of plants as biomonitors of trace element atmospheric pollution analysed by nuclear and related techniques involving 14 participating countries. The CRPs objective was to identify appropriate bioindicators for local and/or regional application and validate them for general air pollution monitoring. Activities included quantification studies, research into spatial and time resolution for particular organisms, and physiological studies. A number of suitable bioindicators were identified in different parts of the globe and tested during the CRP. Sampling strategies were reviewed and the recommended approach adopted by the group. Appropriate sample preparation procedures were assessed and harmonised to the degree allowed by different geographic and climatic conditions in the participating countries. Two interlaboratory comparison exercises were carried out on lichen and moss materials. Results confirmed definite improvement in analytical performance of the participating laboratories, but also revealed possible inconsistencies due to different sample processing procedures. Several monitoring surveys were carried out and consequently pollution maps drawn for extended areas or countries. Overall results confirmed applicability of lower plants for assessing the degree of atmospheric pollution and provided several countries with effective monitoring tools not used before.