OSL signal resetting in young deposits determined with a pulsed photon‐stimulated luminescence (PPSL) unit

Optical Stimulated Luminescence (OSL) is a technique that can be used for dating geological materials deposited within the last half‐million years, including sediments transported by air, water or gravity, as well as rocks heated at high temperatures. Recently, several studies have shown that OSL can also provide information on sediment transport. The pulsed photon‐stimulated luminescence (PPSL) unit (also known as a portable OSL reader) developed by the Scottish Universities Environmental Research Centre is an instrument designed to read luminescence signals from bulk (untreated) sediment samples comprising poly‐mineral and poly‐grain fractions. In this contribution, we evaluate the potential of the PPSL unit to assess the degree of OSL signal resetting in 27 young deposits (<2 ka) transported by different geomorphic agents in volcanic, coastal and fluvial depositional settings located in Mexico. Our results are in agreement with previous findings that used the Risø TL/OSL reader, confirming that sediments transported by debrisflows contain the highest inherited luminesce signals. Infrared stimulation (IRSL) values in volcanic ash, lavas, and sand beach and dune deposits exhibit low scatter. However, with blue stimulation (BLSL) these samples reveal a large degree of scattering, attributed to charge transfer in the case of the coastal deposits and to the low sensitivity of quartz in the case of volcanic material. The luminescence signals of fluvial sediments exhibit a highly scattered distribution in both IRSL and BLSL. We conclude that the use of a PPSL unit is a simple approach to assess the degree of OSL signal resetting in deposits sourced from different geological environments. This research contributes to previous studies that have investigated new applications of the PPSL unit to assist in OSL dating of geological materials.     

A robust mesh optimisation method for multiphase porous media flows

Flows of multiple fluid phases are common in many subsurface reservoirs. Numerical simulation of these flows can be challenging and computationally expensive. Dynamic adaptive mesh optimisation and related approaches, such as adaptive grid refinement can increase solution accuracy at reduced computational cost. However, in models or parts of the model domain, where the local Courant number is large, the solution may propagate beyond the region in which the mesh is refined, resulting in reduced solution accuracy, which can never be recovered. A methodology is presented here to modify the mesh within the non-linear solver. The method allows efficient application of dynamic mesh adaptivity techniques even with high Courant numbers. These high Courant numbers may not be desired but a consequence of the heterogeneity of the domain. Therefore, the method presented can be considered as a more robust and accurate version of the standard dynamic mesh adaptivity techniques.     

An algorithm to resolve γ-rays from charged cosmic rays with DAMPE

The DArk Matter Particle Explorer(DAMPE),also known as Wukong in China,which was launched on 2015 December 17,is a new high energy cosmic ray and γ-ray satellite-borne observatory.One of the main scientific goals of DAMPE is to observe Ge V-Te V high energy γ-rays with accurate energy,angular and time resolution,to indirectly search for dark matter particles and for the study of high energy astrophysics. Due to the comparatively higher fluxes of charged cosmic rays with respect to γ-rays,it is challenging to identify γ-rays with sufficiently high efficiency,minimizing the amount of charged cosmic ray contamination. In this work we present a method to identify γ-rays in DAMPE data based on Monte Carlo simulations,using the powerful electromagnetic/hadronic shower discrimination provided by the calorimeter and the veto detection of charged particles provided by the plastic scintillation detector. Monte Carlo simulations show that after this selection the number of electrons and protons that contaminate the selected γ-ray events at～10 Ge V amounts to less than 1% of the selected sample.Finally,we use flight data to verify the effectiveness of the method by highlighting known γ-ray sources in the sky and by reconstructing preliminary light curves of the Geminga pulsar.     

Irrigated area determination: a case study in the Province of San Juan,Argentina

Much of the central-western region of Argentina, where San Juan Province is located, experiences arid to semi-arid climatic conditions with low average annual rainfall accompanied by substantial evapotranspiration. Consequently, a viable crop industry depends to a large extent upon irrigation from major river systems. Increasing demand for water in the lower basin of the San Juan River is emphasizing the need for more accurate estimates of water used for irrigation. Since the water demand for a particular crop is very closely related to crop area, monitoring the area of crop under irrigation is considered a proxy for the amount of water used. Landsat 5 imagery for the growing season, field data and aerial photographs were used to evaluate crop area.     

Predicting the risk of saltwater contamination of freshwater aquifers during aquifer thermal energy storage

Aquifer thermal energy storage (ATES) is an underground thermal energy storage technology with a large potential to decarbonise the heating and cooling of buildings. ATES installations typically store thermal energy in aquifers that are also exploited for potable water, so a major consideration during development is ensuring that system operation will not lead to groundwater pollution. In this study, the risk of contamination due to upconing of a shallow freshwater/saltwater interface during ATES operation is investigated. Fluid flow, and heat and salt (chloride ion) transport are simulated in a homogeneous aquifer during ATES operation via a well doublet. The impact of geological, hydrological and operational parameters is investigated in a sensitivity analysis. Two new dimensionless numbers are proposed to characterise salt upconing and redistribution during ATES operation and provide a close match to simulated concentrations: C_R,w characterises the contamination risk at the ATES installation, and C_R,d characterises the risk at locations downstream of the ATES installation with respect to background groundwater flow. ATES systems with C_R,w and C_R,d < 10 introduce low risk of contamination in a homogenous aquifer, with chloride concentration at, and downstream of, the ATES system, remaining below the World Health Organisation’s advised limit. ATES installations with C_R,w and C_R,d > 10 cause a rapid increase in aquifer chloride concentration. The results are used to estimate an exclusion distance beyond which ATES system operation will not cause contamination in a homogenous aquifer. The dimensionless parameters proposed allow rapid assessment of the potential for saltwater contamination during ATES operation.

     

Eruption and mass-wasting-induced processes during the late holocene destructive phase of Nevado del Ruiz volcano, Colombia

The November 13, 1985, eruption was characteristic of the Arenas eruptive stage of Nevado del Ruiz, the most recent of a series of twelve eruptive stages that have occurred in the past 11,000 years. Eruptive sequences, deposits and processes similar to that of 1985 have characterized the behavior of Nevado del Ruiz during three major prehistorical and historical eruptive stages: the approximately 3300-3100 yr. B.P. Hedionda, the 16th century Azufrado, and the mid-1800's Lagunillas eruptive stages, that partly destroyed the present Ruiz summit.According to the interpretation of the stratigraphic record of prehistorical eruptions and historical accounts, almost every recent magmatic event was small or short-lived. Nevertheless, rockslide-debris avalanches and catastrophic debris flows were triggered in all the eruptions owing to slope failures related to specific tectonic features of Ruiz volcano and/or to significant interactions between pyroclastic debris and the ice cap. Evidence for headward retreat of avalanche scarps during multiple eruptions reinforce the case that large slope failures can occur repeatedly at a large-volume volcano like Ruiz without reconstruction of the edifice. The latest Ruiz eruptions that involved rockslide-debris avalanches resemble in part the Shiveluch 1964 event for which evidence of lateral blast deposits are lacking, but differ in part from this type because non-eruptive and mass-wasting processes also triggered rockslide-debris avalanches.Many factors render the cluster of domes of the Ruiz summit unstable, including: (1) deeply dissected troughs opened toward the north-northeast (Azufrado), east (Lagunillas), and south (Recio) flanks; (2) strongly hydrothermally altered north and east flanks of the summit; (3) currently glaciated or recently deglaciated, high cliffs; (4) steep unstable margins of the ice cap on the north and east.Thus, in light of its past behavior, a small eruption or an earthquake might trigger catastrophic rockslide-debris avalanches. Furthermore, such avalanches as well as glacial outburst-floods and ice avalanches could induce debris flows by mobilizing weathered, water-saturated, and unconsolidated rocks or deposits.     

天然气水合物稳定带内流体压裂计算的程序耦合方法

海洋天然气水合物稳定带气烟囱结构中存在被水合物充填的裂隙, 表明在自然条件下沉积物中曾发生过流体压裂以及相关的流体流动和水合物形成。在水合物稳定带内实施人为的流体压裂工程, 并联合其他方法(如降压或注热)进行水合物开采, 有望提高开采效率。水合物稳定带内, 无论是自然条件下发生的流体压裂过程, 还是人为实施的流体压裂工程, 都存在水合物反应和沉积物裂隙变形之间的耦合响应。当前, 已有不少数值程序对水合物反应与沉积物弹塑性变形的耦合过程进行了定量研究, 但尚没有数值程序能够计算水合物反应和离散裂隙变形之间的耦合过程。文章将TOUGH+Hydrate程序、IC-FERST和Solidity两者的耦合程序进行了进一步耦合, 为水合物稳定带内的流体压裂计算提供了一种耦合计算方法, 同时通过一个算例初步验证了该耦合计算方法的可行性。验证结果表明, 该耦合计算方法经进一步改进后有望应用于定量研究水合物稳定带内的裂隙变形和水合物反应过程。     

Intra-arc opening and closure of a marginal sea: The case of the Guerrero Terrane (SW Mexico)

Abstract The pre-Oligocene structure of southwest Mexico, south of the trans-Mexico volcanic axis, is investigated from Taxco (Guerrero state, abbreviation: Gro) to the Pacific coast. Three volcano-sedimentary units are recognized; from east to west the calc-alkaline Teloloapan, tholeiitic Arcelia and calc-alkaline Zihuatanejo suites. Structural and stratigraphic data show that the Teloloapan volcanic arc, active during ?Late Jurassic and early Cretaceous, was built upon continental basement. The Teloloapan lavas are overlain by the Albian–Cenomanian Morelos platform carbonates and followed by the Upper Cretaceous Mexcala flysch. In contrast, the Arcelia pillow lavas are associated with sandstones and cherts of Albian-?Cenomanian age. The Zihuatanejo arc was also installed upon continental basement and its magmatic activity was in part coeval with Arcelia magmatism. Unlike the almost undeformed Zihuatanejo volcanic rocks, all the other volcanic units are involved in east-vergent thrusting and recumbent folding associated with ductile tectonics, as well as the Late Cretaceous Mexcala flysch overlying the Morelos platform carbonates. Contrasting with previous views, the present results do not support a major mid-Cretaceous thrusting event in the study area. The new geodynamic interpretation proposed here considers that the Arcelia rocks were formed in a marginal basin situated east of the Zihuatanejo arc. Closure of this basin in Paleocene times is responsible for the east vergent thrust tectonics in SW Mexico.     

Designing and testing a network of passive seismic surveying and monitoring in Dehdasht (South Western Iran)

From August 2016 to July 2017, a passive seismic survey was conducted in South Western Iran as a part of a pilot project aimed to improve the imaging in geologically complex areas. Passive seismic methods have shown to be a useful tool to infer the physical properties of the underground geological structures where traditional hydrocarbon exploration methods are challenging. For this purpose, a dense passive seismic network consisting of 119 three-component borehole seismic stations was deployed over an area of 400 km² around the city of Dehdasht. This paper focuses on the details of the network design, which was devoted to high-resolution seismological applications, including local earthquake tomography and seismic attenuation imaging. In this regard, we describe the instrument types and the station installation procedures used to obtain high-quality data that were used to retrieve three-dimensional models of P- and S-wave velocity and P-wave attenuation in the area using tomographic inversion techniques. We also assess the network performance in terms of the seismic ambient noise levels recorded at each station site, and we revise the horizontal orientation of the sensors using surface waves from teleseismic earthquakes.     

A spherical model for computing polarized radiation in Titan's atmosphere

The Huygens descent through Titan's atmosphere in January 2005 will provide invaluable information about Titan's atmospheric composition and aerosol properties. The Descent Imager/Spectral Radiometer (DISR) will perform upward and downward looking radiation observations at various spectral ranges and spatial resolutions. To prepare the DISR data interpretation we have developed a new model for radiation transfer in Titan's atmosphere. The model solves for the full three-dimensional polarized radiation field in spherical geometry. However, the atmosphere itself is assumed to be spherically symmetric. The model is initialized with a fast-to-compute plane–parallel solution based on the doubling and adding algorithm that incorporates a spherical correction for the incoming direct solar beam. The full three-dimensional problem is then solved using the characteristics method combined with the Picard iterative approximation as described in Rozanov et al. (J. Quant. Spectrosc. Radiat. Transfer 69 (2001) 491). Aerosol scattering properties are calculated with a new microphysical model. In this formulation, aerosols are assumed to be fractal aggregates and include methane gas absorption embedded into the extinction coefficient. The resulting radiance of the model atmosphere's internal field is presented for two prescribed DISR wavelengths.