Natural analogues: a potential approach for developing reliable monitoring methods to understand subsurface CO2 migration processes

One possible way of mitigating carbon dioxide (CO₂) emissions from fossil fuel combustion is using carbon dioxide capture and storage (CCS) technology. However, public perception concerning CO₂ storage in the geosphere is generally negative, being particularly motivated by perceived leakage risks. Therefore, a main issue when attempting to gain public acceptance is ensuring provision of appropriate monitoring practices, aimed at providing health, safety and environmental risk assessment, so that potential risks from CO₂ storage are minimized. Naturally occurring CO₂ deposits provide unique natural analogues for evaluating and validating methods used for the detection and monitoring of CO₂ spreading and degassing into the atmosphere. Geological and hydrological structures of the Cheb Basin (NW Bohemia, Czech Republic) represent such a natural analogue for investigating CO₂ leakage and offer a perfect location at which to verify monitoring tools used for direct investigation of processes along preferential migration paths. This shallow basin dating from the Tertiary age is characterized by up to 300?m thick Neogene sediment deposits and several tectonically active faults. The objectives of this paper are to introduce the CO₂ analogues concept to present the Eger Rift as a suitable location for a natural CO₂ analogue site and to demonstrate to what extent such an analogue site should be used (with a case study). The case study presents the results obtained from a joint application of geoelectrical measurements in combination with soil CO₂ concentration and flux determination methods, for the detection and characterization of natural CO₂ releases at gas seeps (as part of a hierarchic monitoring concept). To highlight discharge-controlling structural near surface features was the initial motivation for the application of geoelectrical measurements. Soil-gas concentration and flux measurement techniques are relatively simple to employ and are valuable methods that can be used to monitor seeping CO₂ along preferential pathways. Joint interpretation of both approaches yields a first insight into fluid paths and reveals that the thickness and permeability of site-specific near surface sedimentary deposits have a great influence upon the spatial distribution of the CO₂ degassing pattern at surface level.     

The regional pressure impact of CO<Subscript>2</Subscript> storage: a showcase study from the North German Basin

A regional scale, showcase saline aquifer CO₂ storage model from the North German Basin is presented, predicting the regional pressure impact of a small industrial scale CO₂ storage operation on its surroundings. The intention of the model is to bridge the gap between generic and site-specific, studying the role of fluid flow boundary conditions and petrophysical parameters typically found in the North German Basin. The numerical simulation has been carried out using two different numerical simulators, whose results matched well. The most important system parameters proved to be the model’s hydrological boundary conditions, rock compressibility, and permeability. In open boundary aquifers, injection-induced overpressures dissipate back to hydrostatic level within a few years. If a geological flow barrier is present on at least one side of the aquifer, pressure dissipation is seriously retarded. In fully closed compartments, overpressures can never fully dissipate, but equilibrate to a compartment-wide remnant overpressure. At greater distances to the injection well, maximum fluid pressures are in the range of a few bar only, and reached several years to decades after the end of the actual injection period. This is important in terms of long-term safety and monitoring considerations. Regional pressure increase impacts the storage capacities of neighbouring sites within hydraulically connected units. It can be concluded that storage capacities may be seriously over- or underestimated when the focus is on a single individual storage site. It is thus necessary to assess the joint storage capacities and pressure limitations of potential sites within the same hydraulic unit.     

High resolution sensors in space and time for determination saltation and creep intensity

The measurement of aeolian sand transport rates at high temporal and spatial resolution is crucial for further progress in testing and developing numerical models of sand movement by wind and in the modelling of sand dunes, ripples and so on. This paper reports the development and field testing of two sand transport sensors. The first one, a webcam commonly used with personal computers, is a new device in aeolian research. The webcam frame transfer is triggered by a sonic anemometer every 0·1 second. Consecutive frames are compared and analysed in real‐time by a computer program. Changes in pixel light intensity exceeding a threshold level are recorded and interpreted as grain movements. The second sensor is a small‐sized ‘Saltiphone’‐type device made of simple loudspeakers with a diameter of 15 mm as used in MP3 players. It can be deployed as a 2 × 3 array of six such devices distributed to enable horizontal and vertical spatial sampling of the sand flux. The devices are tested under field conditions. Both signals and the sum of microphone impacts over 15 minutes are compared to data gathered using a Guelph sand trap, and very good agreement is found. Measurements in a wind tunnel using sieved natural sand indicate that the webcam can be used to infer additional information about the grain size. As an application, the fluid and impact thresholds for aeolian sand transport are investigated in field measurements by analysing the onset and breakdown of saltation in gust and lull intervals of rising and falling wind speeds, respectively. In this way, constitutive equations for sand transport in terms of the wind speed can be tested. If viable, they can be employed to infer estimates for the thresholds by minimizing the root‐mean‐square error between measured and calculated transport data. Copyright © 2012 John Wiley & Sons, Ltd.     

Complexity vs. Simplicity: Groundwater Model Ranking Using Information Criteria

A groundwater model characterized by a lack of field data about hydraulic model parameters and boundary conditions combined with many observation data sets for calibration purpose was investigated concerning model uncertainty. Seven different conceptual models with a stepwise increase from 0 to 30 adjustable parameters were calibrated using PEST. Residuals, sensitivities, the Akaike information criterion (AIC and AICc), Bayesian information criterion (BIC), and Kashyap's information criterion (KIC) were calculated for a set of seven inverse calibrated models with increasing complexity. Finally, the likelihood of each model was computed. Comparing only residuals of the different conceptual models leads to an overparameterization and certainty loss in the conceptual model approach. The model employing only uncalibrated hydraulic parameters, estimated from sedimentological information, obtained the worst AIC, BIC, and KIC values. Using only sedimentological data to derive hydraulic parameters introduces a systematic error into the simulation results and cannot be recommended for generating a valuable model. For numerical investigations with high numbers of calibration data the BIC and KIC select as optimal a simpler model than the AIC. The model with 15 adjusted parameters was evaluated by AIC as the best option and obtained a likelihood of 98%. The AIC disregards the potential model structure error and the selection of the KIC is, therefore, more appropriate. Sensitivities to piezometric heads were highest for the model with only five adjustable parameters and sensitivity coefficients were directly influenced by the changes in extracted groundwater volumes.     

A Comparative Study and Evaluation of Mixing-Height Estimation Based on Sodar-RASS, Ceilometer Data and Numerical Model Simulations

The unmanned aerial system SUMO (Small Unmanned Meteorological Observer) has been used for the observation of the structure and behaviour of the atmospheric boundary layer above the Advent Valley, Svalbard during a two-week period in early spring 2009. Temperature, humidity and wind profiles measured by the SUMO system have been compared with measurements of a small tethered balloon system that was operated simultaneously. It is shown that both systems complement each other. Above 200?m, the SUMO system outperforms the tethered balloon in terms of flexibility and the ability to penetrate strong inversion layers of the Arctic boundary layer. Below that level, the tethered balloon system provides atmospheric profiles with higher accuracy, mainly due to its ability to operate at very low vertical velocities. For the observational period, a numerical mesoscale model has been run at high resolution and evaluated with SUMO profiles reaching up to a height of 1500?m above the ground. The sensitivity to the choice of atmospheric boundary-layer schemes and horizontal resolution has been investigated. A new scheme especially suited for stable conditions slightly improves the temperature forecast in stable conditions, although all schemes show a warm bias close to the surface and a cold bias above the atmospheric boundary layer. During one cold and cloudless night, the SUMO system could be operated nearly continuously (every 30?C45?minutes). This allowed for a detailed case study addressing the structure and behaviour of the air column within and above Advent Valley and its interaction with the local topography. The SUMO measurements in conjunction with a 10-m meteorological mast enabled the identification of a very stable nocturnal surface layer adjacent to the valley bottom, a stable air column in the valley and a strong inversion layer above the summit height. The results indicate the presence of inertial-gravity waves during the night, a feature not captured by the model.     

Evaluation of the Interpretation of Ceilometer Data with RASS and Radiosonde Data

Since 2006 different remote monitoring methods for determining mixing-layer height have been operated in parallel in Augsburg (Germany). One method is based on the operation of eye-safe commercial mini-lidar systems (ceilometers). The optical backscatter intensities recorded with ceilometers provide information about the range-dependent aerosol concentration; gradient minima within this profile mark the tops of mixed layers. Special software for these ceilometers provides routine retrievals of lower atmospheric layering. A second method, based on sodar observations, detects the height of a turbulent layer characterized by high acoustic backscatter intensities due to thermal fluctuations and a high variance of the vertical velocity component. This information is extended by measurements with a radio-acoustic sounding system (RASS) that directly provides the vertical temperature profile from the detection of acoustic signal propagation and thus temperature inversions that mark atmospheric layers. Ceilometer backscatter information is evaluated by comparison with parallel measurements. Data are presented from 2 years of combined ceilometer and RASS measurements at the same site and from comparison with a nearby (60 km) radiosonde for larger-scale humidity information. This evaluation is designed to ensure mixing-layer height monitoring from ceilometer data more reliable.     

The corridor correction concept

Atmospheric delays are contributors to the GNSS error budget in precise GNSS positioning that can reduce positioning accuracy considerably if not compensated appropriately. Both ionospheric and tropospheric delay corrections can be determined with help of reference stations in active GNSS networks. One approach to interpolate these error terms to the user’s location that is employed in Germany’s SAPOS network is the determination of area correction parameters (ACP, German: “Fl?chenkorrekturparameter—FKP”). A 2D interpolation scheme using data from at least 3 reference stations surrounding the rover is employed. A modification of this method was developed which only makes use of as few as 2 reference stations and provides 1D linear correction parameters along a “corridor” in which the user’s rover is moving. We present the results of a feasibility study portraying results from use of corridor correction parameters for precise RTK-like positioning. The differences to the reference coordinates (3D) attained in average for 1 h of data employing selected network nodes in Germany are between 0.8 and 2.0 cm, which compares well with the traditional area correction method that yields an error of 0.7 up to 1.1 cm.     

青藏高原纳木错湖阶沉积的发生特征及环境指示意义

选取的21个湖相沉积物样品分别采自纳木错湖西岸塔吉古日（剖面Nam Co 1）和纳木错湖东北岸的干玛弄附近（剖面Nam Co 2）。按照其发生特征分别将剖面从顶部到底部划分为12个和9个层次,通过分析各个层次沉积物的粒径分布特征、质地组成特点、总有机碳（TOC）、总无机碳（TIC）、总碳（TC）的百分含量、酸度值（pH）、电导率（EC）、含水量等若干环境指示因子。结果显示:无论在比较两个剖面的发生特征,还是环境指示因子特征,剖面Nam Co 1发育程度较剖面Nam Co 2更加发育;通过对测年数据结果分析,两个剖面均形成于晚更新世,并且通过比较各自剖面在垂直方向上的各项环境指示因子的变化趋势,可知纳木错湖流域在晚更新世（约40 ka BP）阶段发生过3次相对明显的冷暖交替过程.     

Continental rift-setting and evolution of Neoproterozoic Sindreth Basin in NW-India

The Neoproterozoic Sindreth Basin, NW India, and its surrounding area represent a half graben structure situated between the undeformed Malani Igneous Suite (MIS) in the west and a corridor of coeval Cryogenian ductile deformation, anatexis and granite intrusion in the east. The main lithologies observed in the basin are conglomerate, fanglomerate, debris flow and lake deposits derived from a nearby continental provenance, intercalated with concurrent mafic and felsic lava flows. Based on geological traverses across the strike of the basin, we propose a three-fold classification comprising Lower Clastic Unit and an Upper Clastic Unit and a Bimodal (basalt–rhyolite) Volcanic Unit separating the two. Tilting due to basin inversion and faulting has been observed; however, the rocks are unmetamorphosed and show undisturbed primary sedimentary features. The stratigraphic record of the basin is characteristic for deposition and magmatism in a fault-related continental setting. Implications of the findings have been discussed in the context of Neoproterozoic crustal dynamics in NW India. This study provides conclusive evidence for a continental setting for Sindreth Basin evolution and contests the recent models of active subduction setting (either back-arc basin or accretionary sediments over a subduction zone).     

A regional groundwater-flow model for sustainable groundwater-resource management in the south Asian megacity of Dhaka,Bangladesh

The water resources that supply most of the megacities in the world are under increased pressure because of land transformation, population growth, rapid urbanization, and climate-change impacts. Dhaka, in Bangladesh, is one of the largest of 22 growing megacities in the world, and it depends on mainly groundwater for all kinds of water needs. The regional groundwater-flow model MODFLOW-2005 was used to simulate the interaction between aquifers and rivers in steady-state and transient conditions during the period 1981–2013, to assess the impact of development and climate change on the regional groundwater resources. Detailed hydro-stratigraphic units are described according to 150 lithology logs, and a three-dimensional model of the upper 400 m of the Greater Dhaka area was constructed. The results explain how the total abstraction (2.9 million m³/d) in the Dhaka megacity, which has caused regional cones of depression, is balanced by recharge and induced river leakage. The simulated outcome shows the general trend of groundwater flow in the sedimentary Holocene aquifers under a variety of hydrogeological conditions, which will assist in the future development of a rational and sustainable management approach.