Hydrogenization of underground storage of natural gas

The intermittent production of the renewable energy imposes the necessity to temporarily store it. Large amounts of exceeding electricity can be stored in geological strata in the form of hydrogen. The conversion of hydrogen to electricity and vice versa can be performed in electrolyzers and fuel elements by chemical methods. The nowadays technical solution accepted by the European industry consists of injecting small concentrations of hydrogen in the existing storages of natural gas. The progressive development of this technology will finally lead to the creation of underground storages of pure hydrogen. Due to the low viscosity and low density of hydrogen, it is expected that the problem of an unstable displacement, including viscous fingering and gravity overriding, will be more pronounced. Additionally, the injection of hydrogen in geological strata could encounter chemical reactivity induced by various species of microorganisms that consume hydrogen for their metabolism. One of the products of such reactions is methane, produced from Sabatier reaction between H₂ and CO₂. Other hydrogenotrophic reactions could be caused by acetogenic archaea, sulfate-reducing bacteria and iron-reducing bacteria. In the present paper, a mathematical model is presented which is capable to reflect the coupled hydrodynamic and bio-chemical processes in UHS. The model has been numerically implemented by using the open source code DuMuX developed by the University of Stuttgart. The obtained bio-chemical version of DuMuX was used to model the evolution of a hypothetical underground storage of hydrogen. We have revealed that the behavior of an underground hydrogen storage is different than that of a natural gas storage. Both, the hydrodynamic and the bio-chemical effects, contribute to the different characteristics.     

A modified monthly degree‐day model for evaluating glacier runoff changes in China. Part I: model development

The retreat of mountain glaciers and ice caps has dominated the rise in global sea level and is likely to remain an import component of eustatic sea‐level rise in the 21st century. Mountain glaciers are critical in supplying freshwater to populations inhabiting the valleys downstream who heavily rely on glacier runoff, such as arid and semi‐arid regions of western China. Owing to recent climate warming and the consequent rapid retreat of many glaciers, it is essential to evaluate the long‐term change in glacier melt water production, especially when considering the glacier area change. This paper describes the structure, principles and parameters of a modified monthly degree‐day model considering glacier area variation. Water balances in different elevation bands are calculated with full consideration of the monthly precipitation gradient and air temperature lapse rate. The degree‐day factors for ice and snow are tuned by comparing simulated variables to observation data for the same period, such as mass balance, equilibrium line altitude and glacier runoff depth. The glacier area–volume scaling factor is calibrated with the observed glacier area change monitored by remote sensing data of seven sub‐basins of the Tarim interior basin. Based on meteorological data, the glacier area, mass balance and runoff are estimated. The model can be used to evaluate the long‐term changes of melt water in all glacierized basins of western China, especially for those with limited observation data. Copyright © 2011 John Wiley & Sons, Ltd.     

Improving a subgrid runoff parameterization scheme for climate models by the use of high resolution data derived from satellite observations

In this study it is shown that the availability of a very high resolution dataset of land surface characteristics leads to the improvement of a surface runoff parameterization scheme. The improved parameterization scheme was developed for application in global and regional climate models and is a further development of the Arno scheme that is widely used in climate models. Here, surface runoff is computed as infiltration excess from a "bucket" type reservoir which takes the subgrid variability of soil saturation within a model gridbox into account. Instead of prescribing a distribution of subgrid scale soil water capacities as in the original Arno scheme, the array of high resolution soil water capacities taken from a global 1 km dataset of land surface parameters is used to obtain individual fractional saturation curves for each model gridbox. From each saturation curve, the three parameters (a shape parameter describing the shape of the subgrid distribution of soil water capacities, subgrid minimum and maximum soil water capacity) required in the modified formulation of the scheme are derived via optimization. As in the original Arno scheme applied in the ECHAM general circulation model and the REMO regional climate model, topography variations will influence the distribution of saturated subgrid areas within a model gridbox. At most gridboxes the net effect of these changes is such that more runoff is produced for high soil water contents and less runoff for low soil water contents. A validation of simulated discharge computed with a simplified land surface scheme applied to reanalysis data of the European Centre for Medium-Range Weather Forecasts and a hydrological discharge model has shown that these changes lead to a more realistic simulation of the annual cycle of discharge for several catchments. In particular this could be shown for the Yangtze Kiang and Amur catchments where adequate input data are available.     

Gold mineralisation throughout about 45 Ma of Archaean orogenesis: protracted flux of gold in the Golden Mile, Yilgarn craton, Western Australia

The Golden Mile deposit was discovered in 1893 and represents today the largest Archaean orogenic lode gold system in the world (50 M oz produced gold). The Golden Mile deposit comprises three major styles of gold mineralisation: Fimiston, Oroya and Charlotte styles. Fimiston-style lodes formed at 250 to 350 °C and 100 to 200 MPa and are controlled by brittle–ductile fault zones, their subsidiary fault zone and vein networks including breccias and open-cavity-infill textures and hydrothermally altered wall rock. Fimiston lodes were formed late D₁, prior to D₂ regional upright folding. Hydrothermal alteration haloes comprise a progression toward the lode of diminishing chlorite, an increase in sericite and in Fe content of carbonates. Lodes contain siderite, pyrite, native gold, 17 different telluride minerals (Au–Ag tellurides contain ~25% of total gold), tourmaline, haematite, sericite and V-rich muscovite. Oroya-style lodes formed at similar P–T conditions as the Fimiston lodes and are controlled by brittle–ductile shear zones, associated dilational jogs that are particularly well developed at the contact between Paringa Basalt and black shale interflow sedimentary rocks and altered wall rock. The orebodies are characterised by micro-breccias and zones of intense shear zone foliation, very high gold grades (up to 100,000 g/t Au) and the common association of tellurides and vanadian mica (green leader). Oroya lodes crosscut Fimiston lodes and are interpreted to have formed slightly later than Fimiston lodes as part of one evolving hydrothermal system spanning D₁ and D₂ deformation (ca. 2,675–2,660 Ma). Charlotte-style lodes, exemplified by the Mt Charlotte deposit, are controlled by a sheeted vein (stockwork) complex of north-dipping quartz veins and hydrothermally altered wall rock. The Mt Charlotte orebody formed at 120 to 440 °C and 150 to 250 MPa during movement along closely spaced D₄ (2,625 Ma) and reactivated D₂ faults with the quartz granophyre in the Golden Mile Dolerite exerting a strong lithological control on gold mineralisation. Veins consist of quartz–carbonate–minor scheelite, and wall-rock alteration comprises chlorite destruction and growth of ferroan carbonate–sericite–pyrite–native gold. Pyrite–pyrrhotite is zoned on the scale of vein haloes and of the entire mine, giving a vertical temperature gradient of 50–100 °C over 1,000 vertical metres. The structural–hydrothermal model proposed consists of four major stages: (1) D₁ thrusting and formation of Fimiston-style lodes, (2) D₂ reverse faulting and formation of Oroya-style lodes, (3) D₃ faulting and dissecting of Fimiston- and Oroya-style lodes, and (4) D₄ faulting and formation of Mt Charlotte-style sheeted quartz vein system. The giant accumulation of gold in the Golden Mile deposit was formed due to protracted gold mineralisation throughout episodes of an Archaean orogeny that spanned about 45 Ma. Fluid conduits formed early in the tectonic history and persisted throughout orogenesis with the plumbing system showing a rare high degree of focussing, efficiency and duration. In addition to the long-lasting fluid plumbing system, the wide variety of transient structural and geochemical traps, multiple fluid sources and precipitation mechanism contributed towards the richest golden mile in the world.Editorial handling: B. Lehmann     

The geology, petrology, palynology and geochemistry of Permian coal basins in Tanzania: 2. Songwe-Kiwira Coalfield

This study provides coal quality, petrological, palynological and geochemical (Rock Eval) data on Permian coal seams and associated shales and mudstones of the Karoo Supergroup of the Songwe-Kiwira Coalfield, Tanzania. The coal seams, which have a cumulative thickness of 6.80 m, occur in the shale–coal–sandstone facies of the Mchuchuma Formation of Artinskian to Kungurian(?) age.Coal quality data (calorific values, volatile matter contents) and vitrinite reflectances indicate high volatile C bituminous to high volatile A bituminous coals, having relatively high ash yields (22–49 wt.%) and highly variable sulphur contents (0.17–9.2 wt.%). They could be used to fuel small-scale power generation units thereby providing electricity to nearby towns and villages. Also, the coals could be used as a substitute for wood, which is becoming increasingly scarce. In rural Tanzania, charcoal is still the main energy source for cooking, and wood is used extensively in brick kilns and for making roofing tiles.Petrological analysis indicated that the coals are dominated by dull to banded dull lithotypes, with seams at the base of the Mchuchuma Formation enriched in inertinite macerals (up to 83 vol.%), whereas up-section vitrinite contents increase. Palynological analyses indicated that the assemblage in the lower Mchuchuma Formation (Scheuringipollenites assemblage) is dominated by trilete spores, whereas in the remainder of the section, non-taeniate disaccates dominate (Scheuringipollenites–Protohaploxypinus assemblage). Facies critical macerals suggest for most seams a marsh/wet forest swamp depositional setting, which is consistent with the palynological data.Rock Eval analyses indicate type II/III kerogen, with Tmax (°C) values ranging from 426 to 440, corresponding to the early stage of hydrocarbon generation. Thermal Alteration Indices (2 to 2+) and vitrinite reflectance levels (0.60–0.83 Ro (%) support the Rock Eval maturity assessment, and despite the predominance of terrestrial-derived organic matter, there is evidence of oil generation and expulsion in the form of cavity and fracture filling exsudatinite.     

Gradient in the climate change signal of European discharge predicted by a multi-model ensemble

In order to perform hydrological studies on the PRUDENCE regional climate model (RCM) simulations, a special focus was put on the discharge from large river catchments located in northern and central Europe. The discharge was simulated with a simplified land surface (SL) scheme and the Hydrological Discharge (HD) model. The daily fields of precipitation, 2 m temperature and evapotranspiration from the RCM simulations were used as forcing. Therefore the total catchment water balances are constrained by the hydrological cycle of the different RCMs. The validation of the simulated hydrological cycle from the control simulations shows that the multi-model ensemble mean is closer to the observations than each of the models, especially if different catchments and hydrological variables are considered. Therefore, the multi-model ensemble mean can be used to largely reduce the uncertainty that is introduced by a single RCM. This also provides more confidence in the future projections for the multi-model ensemble means. The scenario simulations predict a gradient in the climate change signal over Northern and Central Europe. Common features are the overall warming and the general increase of evapotranspiration. But while in the northern parts the warming will enhance the hydrological cycle leading to an increased discharge, the large warming, especially in the summer, will slow down the hydrological cycle caused by a drying in the central parts of Europe which is accompanied by a reduction of discharge. The comparison of the changes predicted by the multi-model ensemble mean to the changes predicted by the driving GCM indicates that the RCMs can compensate problems that a driving GCM may have with local scale processes or parameterizations.     

The role of the simulation setup in a long-term high-resolution climate change projection for the southern African region

High-resolution regional climate change simulations have proven to offer an added value compared to available global climate model simulations. However, over many regions of the globe, long-term high-resolution climate change projections are rather sparse. We present a transient high-resolution climate change projection with the regional climate model with the regional climate model REMO over the southern African region, following the SRES A1B emission scenario. The simulation was conducted at 18?km grid spacing for the period from 1960 to 2100, making it to the longest available climate change projection at such a high resolution for the region. In the first part of the study, we focus on the impact of the model setup on the simulated rainfall over the southern African region. In the standard setup, we used the output of the global climate model ECHAM5/MPIOM directly to force REMO. This setup led to a very strong wet bias over the region. Changing it to the double-nesting setup significantly reduced this bias, but a substantial wet bias still persists. The remaining bias could partly be attributed to a warm bias in the SST forcing over the southern Atlantic Ocean. Thus, we applied an SST correction based on the anomaly approach to the data, which led to a further improvement of the rainfall simulation. As the SST bias in the southern Atlantic is a common feature of all global climate models assessed by the IPCC, we recommend the chosen model setup, including the SST correction, as general procedure for dynamical downscaling studies over the southern African region. In the second part, we present the projected spatial and temporal changes of temperature and precipitation, including several rainfall characteristics, over the southern African region. Herby we compare the projections of the high-resolution REMO simulation to those of the forcing regional and global models. We generally find that for temperature the magnitude of the projected changes of the regional model only slightly differs from the GCM projection; however, the spatial patterns are much better resolved in the RCM projections. For precipitation, REMO shows a more intense drying toward the end of the twenty-first century than it is simulated by the global model. This can have a major influence when investigating the impacts of future climate change on a regional or even local scale. In combination with the improved spatial patterns, the application of high-resolution climate change information could therefore improve the results of such applications.     

Extremes and predictability in the European pre-industrial climate of a regional climate model

A high-resolution pre-industrial control simulation with the regional climate model REMO is analyzed in detail for different European subregions. To our knowledge, this is the first long pre-industrial control simulation by a regional climate model as well as at comparable resolution. We assess the ability of the climate model to reproduce the observed climate variability in various parts of the continent. In order to investigate the representation of extreme events in the model under pre-industrial greenhouse gas concentrations, selected seasons are examined with regard to the atmospheric circulation and other climatic characteristics that have contributed to the occurrences. A special focus is dedicated to land-atmosphere interactions. Extreme seasons are simulated by the model under various circumstances, some of them strongly resemble observed periods of extraordinary conditions like the summer 2003 or autumn 2006 in parts of Europe. The regional perspective turns out to be of importance when analyzing events that are constituted by meso-scale atmospheric dynamics. Moreover, the predictability of the European climate on seasonal to decadal time scales is examined by relating the statistics of surface variables to large-scale modes of variability impacting the North Atlantic sector like the Meridional Overturning Circulation, the El Niño Southern Oscillation, and the North Atlantic Oscillation. For this purpose, we introduce a measure of tail dependence that quantifies the correlation between extreme values in two variables that describe the state of the climate system. Significant dependence of extreme events can be detected in various situations.     

Coal mine workers' pneumoconiosis

 The fibrogenic and cytotoxic potential of coal mine dust is independent of the amount of quartz and other inorganic parameters. Results of coal petrographical and organic geochemical investigations of coals and coal mine dust from the Ruhr and Ibbenbüren Region of Germany demonstrate variations of organic dust amount possibly influencing these noxious properties. Coal mine dust of high rank coals is characterized by a pronounced fibrogenic risk. This risk, independent of the quantity of quartz, is probably based on shape variations of different coal macerals. With increasing coalification of the corresponding seam, the vitrinite is enriched in its dust; however, lower concentrations have been determined for inertinite. Vitrinite shows constant shapes and sizes independent of the rank of coal. Inertinite particles with elongated to fibrous shapes tend to larger sizes with increasing coalification. Strikingly, coal mine dust from miners' lungs with high degrees of coal mine workers' pneumoconiosis (CWP) is enriched in inertinite. In contrast, high cytotoxicities in cell tests are known for coal mine dust from low coalified coals. High concentrations of phenolic compounds can be extracted by dichloromethane from low coalified coal mine dust. These compounds, which are characterized by a high water solubility and therefore high bioavailability, explain the high cytotoxicities of coal mine dust. Contamination of dust by diesel emissions in the coal mine can act as additionally supporting parameters for extended cytotoxicities. Received: 18 August 1995 · Accepted: 15 October 1995