On natural and laboratory generated dose response curves for quartz of different grain sizes from Romanian loess

One of the most important foundations of luminescence dating is the assumption that the growth of the luminescence signal in nature can be reproduced under laboratory conditions by performing irradiations with a calibrated beta or gamma source. When optically stimulated luminescence (OSL) of quartz with a dominant fast component is measured using the single aliquot regenerative dose (SAR) protocol, laboratory dose response curves that display continuing growth at high doses are increasingly reported in literature. In this study we investigate fine (4–11 μm) and coarse (63–90 μm) quartz extracted from 25 samples taken from L1, S1 and L2 units from the loess-palaeosol section at Costineşti in Romania. Our results indicate that the growth of the OSL signal in nature does not correspond to the laboratory generated laboratory dose response curve. The growth of the signal in nature is consistent with a single saturating exponential function, with the signal of coarse grains starting to saturate at 100–200 Gy, and for fine grains at 200–300 Gy, respectively. Laboratory dose response curves continue to grow for high doses (>300 Gy) for both quartz fractions. The differences observed between the natural and the laboratory dose response for the two quartz fractions are believed to be a cause for the different chronologies previously reported using the two grain sizes of quartz on Romanian loess. In addition, we have applied the single aliquot regeneration and added dose (SARA) procedure to both fine and coarse grains from the youngest sample. Our findings question the reliability of obtaining high equivalent doses for quartz samples displaying laboratory dose response curves obtained by the SAR protocol for which a single saturating exponential model does not describe the data.     

A Bayesian hierarchical age model for optical dating of single grains of quartz

In optical dating, the last time that a sample of sediment was exposed to sunlight is determined by dividing its equivalent dose (D_e) by the dose rate. For single-grain dating, the sample D_e is based on the statistical analysis of the distribution of D_e values estimated for individual grains, whereas the dose rate is usually determined from measurements of the environmental radioactivity of the bulk sample, together with allowances for radiation sources internal to the grains and cosmic rays. Conventionally, the D_e and dose rate are measured and analysed separately to produce an estimate of the depositional age of a sample, but this approach may result in loss of information because distributions of single-grain D_e values are influenced by several factors. Existing statistical models do not incorporate all the key information contributing to age estimation, such as the pattern and scale of dispersion of single-grain D_e values and dose rates, the associated measurement uncertainties, the effect of natural variability among grains, and the outlier probabilities of D_e and dose rate estimates. Here we propose an empirical Bayesian hierarchical age model (BHAM) for optical dating of quartz samples that incorporates the above information to estimate their depositional ages. The BHAM is based on the implementation of standardised growth curve and L_nT_n methods to integrate information on the full distribution of single-grain D_e values, sources of measurement uncertainty, beta-dose heterogeneity (observed or modelled), and detection of outliers. We present the results of validation tests using data sets of optically stimulated luminescence measurements and dose rates obtained for quartz samples dated previously from Denisova Cave (Russia) and for simulated samples. We conclude that the BHAM represents a robust and flexible approach to dealing with data for single grains of quartz within a Bayesian hierarchical framework and is suitable for application to sediments deposited in a variety of depositional settings.     

Progress and pitfalls of the SAR protocol for the quartz violet stimulated luminescence (VSL) signal: A case study from Sardinia

In the quest for extending the upper age limit of optically stimulated luminescence (OSL) dating of quartz, it was shown that violet stimulated luminescence (VSL) may be a promising candidate. However, difficulties in the application of the single aliquot regenerative dose (SAR) protocol for VSL has been reported in previous studies. In this study, a set of experiments was carried out to investigate the behaviour of the VSL signal under different measurement conditions with the aim to improve the SAR VSL protocol in terms of passing the procedural tests of the SAR protocol, i.e., recycling ratio, recuperation and dose recovery. The validity of an optimised SAR protocol was then evaluated on four coarse-grained quartz samples from Sardinia, for which previously reported OSL and post-IR IRSL ages are available. Our result showed that all measured aliquots meet the recuperation, recycling ratio and dose recovery criteria, indicating that the proposed protocol is suitable for the studied samples. The obtained VSL and reference ages agree within uncertainty (2σ) for most of the samples except one sample with a largest expected equivalent dose of ∼320 Gy, for which the VSL significantly (∼50%) underestimates despite satisfactory dose recovery result. This underestimation is most likely due to trapping sensitivity change between the natural and all subsequent regenerated VSL signals induced by the first preheating in the SAR procedure, which is likely dose dependent. To minimise sensitivity change, the sensitivity-corrected multiple aliquot regenerative dose (SC-MAR) protocol was applied. The SC-MAR protocol yields VSL ages in agreement with reference ages for two samples (1σ); the other two samples show overestimation. We show that the observed age overestimation can be explained by incomplete resetting of the natural signal after 7 days SOL2 exposure used for the SC-MAR procedure. Bleaching experiments confirm variable bleaching behaviour among different samples. However, it is unclear whether this different bleaching behaviour arises from measuring samples with different bleaching histories, mineralogical compositions, or VSL source trap properties.     

Bayesian analysis of individual and systematic multiplicative errors for estimating ages with stratigraphic constraints in optically stimulated luminescence dating

Many dating techniques include significant error terms which are not independent between samples to date. This is typically the case in Optically Stimulated Luminescence (OSL) dating where the conversion from characteristic equivalent doses to the corresponding ages using the annual dosimetry data includes error terms that are common to all produced datings. Dealing with these errors is essential to estimate ages from a set of datings whose chronological ordering is known. In this work, we propose and we study a Bayesian model to address this problem. For this purpose, we first consider a multivariate model with multiplicative Gaussian errors in a Bayesian framework. This model relates a set of characteristic equivalent doses to the corresponding ages while taking into account for the systematic and non-systematic errors associated to the dosimetry. It thus offers the opportunity to deal properly with stratigraphic constraints within OSL datings, but also with other datings possessing errors which are independent from systematic errors of OSL (e.g. radiocarbon). Then, we use this model to extend an existing Bayesian model for the assessment of characteristic equivalent doses from Single Aliquot and Regenerative (SAR) dose measurements. The overall Bayesian model leads to the joint estimation of all the variables (which include all the dose–response functions and characteristic equivalent doses) of a sequence of, possibly heterogeneous, datings. We also consider a more generic solution consisting in using directly the age model from a set of characteristic equivalent dose estimates and their associated standard errors. We finally give an example of application on a set of five OSL datings with stratigraphic constraints and observe a good adequacy between the two approaches.     

The potential of cryptotephra and OSL dating for refining the chronology of open-air archaeological windblown sand sites: A case study from Mirkowice 33, northwest Poland

The discovery of a cryptotephra (nonvisible volcanic horizon) in a windblown sand archaeological site in Poland highlights how luminescence and tephrostratigraphy may combine to better refine the chronology of such sites. In this study we identify a cryptotephra horizon which on the basis of major and minor element geochemistry and an OSL age of 2.3 ± 0.1 ka is correlated to the Glen Garry tephra. The different methodological strengths of OSL and tephrostratigraphy may be harnessed to counter the limitations of a single approach to produce a more secure chronology. Although in this study the tephra deposition event is shown to post-date the archaeological activity, the methodological approach is clearly demonstrated. Further investigations will reveal if cryptotephra layers are commonly preserved in such environmental settings. If this is so then future applications of this approach may prove to be more widely applicable.     

Modelling dose rate to single grains of quartz in well-sorted sand samples: The dispersion arising from the presence of potassium feldspars and implications for single grain OSL dating

Single grain OSL has become a widely used approach in Quaternary geochronology. However, the origins of D_e distributions and the sources of variation in individual dose estimates are still poorly understood. The amount of scatter in these distributions on top of the known uncertainties in measurement and analysis is defined by overdispersion and this quantity is generally used for weighting individual D_e values to calculate a central equivalent dose. In this study, we address the nature and amount of different sources of dispersion in quartz single grain D_e estimates, by (i) using appropriate statistical tools to characterize D_e populations and (ii) modelling, with a specifically designed Geant4 code, dose rate distributions arising from the presence of potassium feldspar grains in well-sorted sands. The model uses Monte Carlo simulations of beta emissions and interactions in a random close packing of quartz and feldspar spheres representing a sand sample. Based on the simulation results, we explain the discrepancy between intrinsic and natural overdispersion values in a well-bleached sample, thus validating the model. The three parameters having the most influence on dispersion in dose rate distributions, and modelled in this study, appear to be grain size, potassium content and total dose rate.Finally an analysis of measurement uncertainties and other sources of variations in equivalent dose estimates lead us to conclude that all age models (both logged and unlogged) which include an overdispersion value to weight individual D_e values rely mainly on unknown parameters; this ignorance may lead to an inadvertent bias in D_e estimates. Assuming counting statistics make a small contribution to dispersion (as is often the case), we suggest that in some cases it is most appropriate to use unweighted averages of equivalent doses when dividing by commonly measured average dose rates.     

High-resolution dating of Co-rich crusts: A comparative study using the methods of orbital pacing and ~(230)Th_(ex)/~(232)Th dating

Due to their slow growth rates, seamount Co-rich crusts are very difficult to date with high resolution and precision. This paper is to test the use of orbital pacing on the growth profile of crusts to determine high-resolution age and growth rate. Crust CB14 from the central Pacific Ocean was selected for this study. We first examined the growth pattern in detail under a reflected-light microscope and ascertained that the growth environment was stable for the sub-layer 1 (0-3 mm). We then used electron mic...     

A study on textural characteristics,heavy mineral distribution and grain-microtextures of recent sediment in the coastal area between the Sarada and Gosthani rivers,east coast of India

The current study aimed to describe textural characteristics, heavy mineral composition, and grain microtextures of the sediment from three micro-environments (foreshore, berm, and dune). A total of forty-one (41) representative surficial sediment samples have been collected from fifteen (15) locations along the beach area between the Sarada and Gosthani rivers on the east coast of India, where the length of the stretch is more than 100 km. The study reveals that most of the coastal sediment is medium to fine sand with relatively high ratios of coarse sand at Yarada beach, and the nature of the sediment is moderately to well sorted. These characteristics indicate a high energy environment. The heavy mineral analysis of the sediment in the current study was done for coarse (+60#) and fine (+230#) size fractions. Studying the weight percentage (WT%) reveals that a high percentage of heavy minerals is associated with fine fractions. Ilmenite, sillimanite, garnet, zircon, and rutile are the major heavy minerals identified in the current investigation. The concentrations of these heavy minerals show great variations from south to north of the study area. From an economic point of view, a considerable amount of heavy minerals (average 48.41%) are present on both sides (north and south) of the Gosthani River mouth. In the Sarada Estuary, the concentration of the economic heavy minerals was found under the minimum economic range. The grain microtextures of the major heavy minerals from the different locations along the study area demonstrate the variation in grain microtextures, which is controlled by the chemical and mechanical processes. These microtextures reflect moderate to high wave energy on the beach area, in addition to high mechanical impact on the grains from the estuary point.     

A comparative study of fuzzy weights of evidence and random forests for mapping mineral prospectivity for skarn-type Fe deposits in the southwestern Fujian metallogenic belt,China

Recent studies have pointed out that the widespread iron deposits in southwestern Fujian metallogenic belt (SFMB) (China) are skarn-type deposits associated with the Yanshanian granites. There is still excellent potential for mineral exploration because large areas in this belt are covered by forest. A new predictive model for mapping skarn-type Fe deposit prospectivity in this belt was developed and focused on in this study, using five criteria as evidence: (1) the contact zones of Yanshanian granites (GRANITE); (2) the contact zones within the late Paleozoic marine sedimentary rocks and the carbonate formations (FORMATION); (3) the NE-NNE-trending faults (FAULT); (4) the zones of skarn alterations (SKARN); and (5) the aeromagnetic anomaly (AEROMAGNETIC). The fuzzy weights of evidence (FWofE) method, developed from the classical weights of evidence (WofE) and based on fuzzy sets and fuzzy probabilities, could provide smaller variances and more accurate posterior probabilities and could effectively minimize the uncertainty caused by omitted or wrongly assigned data and be more flexible than the WofE. It is an efficient and widely used method for mineral potential mapping. Random forests (RF) is a new and useful method for data-driven predictive mapping of mineral prospectivity method, and needs further scrutiny. Both prospectivity results respectively using the FWofE and RF methods reveal that the prediction model for the skarn-type Fe deposits in the SFMB is successful and efficient. Both methods suggested that the GRANITE and FORMATION are the most valuable evidence maps, followed by SKARN, AEROMAGNETIC, and FAULT. This is coincident with the skarn-type Fe deposit mineral model in the SFMB. The unstable performance experienced when FORMATION was omitted might indicate that the highest uncertainty and risk in follow-up exploration is related to the sequences. In addition, the performance of the RF method for the skarn-type Fe deposits prospectivity in the SFMB is better than the FWofE; therefore, it could be used to guide further exploration of skarn-type Fe prospects in the SFMB.     

High-resolution dating of Co-rich crusts: A comparative study using the methods of orbital pacing and <Superscript>230</Superscript>Th<Subscript>ex</Subscript>/<Superscript>232</Superscript>Th dating

Due to their slow growth rates, seamount Co-rich crusts are very difficult to date with high resolution and precision. This paper is to test the use of orbital pacing on the growth profile of crusts to determine high-resolution age and growth rate. Crust CB14 from the central Pacific Ocean was selected for this study. We first examined the growth pattern in detail under a reflected-light microscope and ascertained that the growth environment was stable for the sub-layer 1 (0–3 mm). We then used electron microprobe line-scanning to obtain elemental profiles. The pattern of the power spectrum analysis of the Al-profile revealed that there are significant cycles of 113.9, 87.8, 51.5, 42.2 and 25.8 μm. These cycles correspond to the Milankovitch cycles of 53.1, 41, 24, 19.7 and 12 ka, respectively, and yield the growth rate of about 2.14 mm/Ma and an age of about 1.40 Ma for the boundary between the sub-layer 1 and sub-layer 2. We also used a drilling machine with a numerically controlled drive to obtain high-resolution samples at 0.1mm intervals, and used the ²³⁰Th_ex/²³²Th method to date the samples. For the uppermost 1.3 mm, the growth rate was about 2.15 mm/Ma, and the age for the layer at the depth of 3 mm was about 1.40 Ma, which coincides perfectly with the results obtained from orbital pacing. Thus, it is considered that orbital pacing is a new and effective method to determine the growth rate of the seamount Co-rich crust. This method is applicable for establishing a high-resolution age frame for the crusts of the world’s oceans. Supported by China Ocean Mineral Resources R & P Association (Grant No. DY105-01-01-08) and National Natural Science Foundation of China (Grant Nos. 40106005, 40476050)     

A comparative study of buried structure in soil subjected to blast load using 2D and 3D numerical simulations

The response of underground structures subjected to subsurface blast is an important topic in protective engineering. Due to various constraints, pertinent experimental data are extremely scarce. Adequately detailed numerical simulation thus becomes a desirable alternative. However, the physical processes involved in the explosion and blast wave propagation are very complex, hence a realistic and detailed reproduction of the phenomena would require sophisticated numerical models for the loading and material responses. In this paper, a fully coupled numerical model is used to simulate the response of a buried concrete structure under subsurface blast, with emphasis on the comparative performance of 2D and 3D modeling schemes. The explosive charge, soil medium and the RC structure are all incorporated in a single model system. The SPH (smooth particle hydrodynamics) technique is employed to model the explosive charge and the close-in zones where large deformation takes place, while the normal FEM is used to model the remaining soil region and the buried structure. Results show that the 2D model can provide reasonably accurate results concerning the crater size, blast loading on the structure, and the critical response in the front wall. The response in the remaining part of the structure shows noticeable differences between the 2D and 3D models. Based on the simulation results, the characteristics of the in-structure shock environment are also discussed in terms of the shock response spectra.     

A comprehensive assessment framework for attributing trends in streamflow and groundwater storage to climatic and anthropogenic changes: A case study in the typical semi-arid catchments of southern India

The clearest signs of hydrologic change can be observed from the trends in streamflow and groundwater levels in a catchment. During 1980–2007, significant declines in streamflow (−3.03 mm/year) and groundwater levels (−0.22 m/year) were observed in Himayat Sagar (HS) catchment, India. We examined the degree to which hydrologic changes observed in the HS catchment can be attributed to various internal and external drivers of change (climatic and anthropogenic changes). This study used an investigative approach to attribute hydrologic changes. First, it involves to develop a model and test its ability to predict hydrologic trends in a catchment that has undergone significant changes. Second, it examines the relative importance of different causes of change on the hydrologic response. The analysis was carried out using Modified Soil and Water Assessment Tool (SWAT), a semi-distributed rainfall-runoff model coupled with a lumped groundwater model for each sub- catchment. The model results indicated that the decline in potential evapotranspiration (PET) appears to be partially offset by a significant response to changes in rainfall. Measures that enhance recharge, such as watershed hydrological structures, have had limited success in terms of reducing impacts on the catchment-scale water balance. Groundwater storage has declined at a rate of 5 mm/y due to impact of land use changes and this was replaced by a net addition of 2 mm/y by hydrological structures. The impact of land use change on streamflow is an order of magnitude larger than the impact of hydrological structures and about is 2.5 times higher in terms of groundwater impact. Model results indicate that both exogenous and endogenous changes can have large impacts on catchment hydrology and should be considered together. The proposed comprehensive framework and approach demonstrated here is valuable in attributing trends in streamflow and groundwater levels to catchment climatic and anthropogenic changes.     

南方红壤区长期水土流失综合治理显著减少河流输沙——以鄱阳湖流域赣江上游平江为例

长时间序列水沙数据分析有助于科学评价流域尺度水土保持减流减沙效应,为科学开展区域水土保持成效评估提供实现路径.本文选取南方红壤区水土保持工程最为集中、持续时间最长的典型流域——平江流域,基于1975-2014年的逐日降水量、流量、含沙量数据,综合应用预置白M-K检验法、Theil-Sen趋势度估计法、Pettitt检验...