Modelling the poroelasticity of rocks and ice

For many geophysical and astrophysical applications the relationship between physical bulk properties and porosity is needed. As part of the preparation for the 'Rosetta' comet rendezvous mission, a simple model for the porosity dependence of the elastic properties of granular media, i.e. the elastic moduli and the propagation velocities of elastic waves, has been developed based on textural properties and the contact stiffness of the constituent particles. It is shown that the derived relationships fit very well with sandstone data. The model is also consistent with data for snow and ice and is in agreement with the transformation mechanisms from snow to ice. A short review shows the relevance to another physical bulk property, the thermal conductivity.     

Another Simple Method of Spectral Correction to Obtain Robust Eddy-Covariance Results

The eddy-covariance method is the state-of-the-art tool to measure mass and energy fluxes, though many measuring systems (particularly closed-path systems) show strong spectral attenuation that causes significant underestimates of actual mass fluxes. The standard way to correct these underestimates is to use theoretically derived transfer functions. Practical implementation has shown that the attenuation, especially of closed-path systems, is underestimated by this method, even when the additional tube damping is considered. This paper introduces and tests three alternative site-specific and hardware specific correction methods based on spectral analysis, which typically enhance—additionally to the classical correction—the calculated mass fluxes based on closed-path analysers, as their inevitable tubing and the filters used cause additional flow attenuation. Two of the three methods are based on a site-specific and hardware specific transfer functions, the third uses direct comparison of cospectra. Primarily the methods based on transfer functions proved to be easy-to-handle once established for the specific set-up and measurement site. They represent practical and robust methods to correct for spectral attenuation.     

Seasonal to interannual climate predictability in mid and high northern latitudes in a global coupled model

The upper limit of climate predictability in mid and high northern latitudes on seasonal to interannual time scales is investigated by performing two perfect ensemble experiments with the global coupled atmosphere–ocean–sea ice model ECHAM5/MPI-OM. The ensembles consist of six members and are initialized in January and July from different years of the model’s 300-year control integration. The potential prognostic predictability is analyzed for a set of oceanic and atmospheric climate parameters. The predictability of the atmospheric circulation is small except for southeastern Europe, parts of North America and the North Pacific, where significant predictability occurs with a lead time of up to half a year. The predictability of 2 m air temperature shows a large land–sea contrast with highest predictabilities over the sub polar North Atlantic and North Pacific. A combination of relatively high persistence and advection of sea surface temperature anomalies into these areas leads to large predictability. Air temperature over Europe, parts of North America and Asia shows significant predictability of up to half a year in advance. Over the ice-covered Arctic, air temperature is not predictable at time scales exceeding 2 months. Sea ice thickness is highly predictable in the central Arctic mainly due to persistence and in the Labrador Sea due to dynamics. Surface salinity is highly predictable in the Arctic Ocean, northern North Atlantic and North Pacific for several years in advance. We compare the results to the predictability due to persistence and show the importance of dynamical processes for the predictability.     

Decadal fingerprints of freshwater discharge around Greenland in a multi-model ensemble

The recent increase in the rate of the Greenland ice sheet melting has raised with urgency the question of the impact of such a melting on the climate. As former model projections, based on a coarse representation of the melting, show very different sensitivity to this melting, it seems necessary to consider a multi-model ensemble to tackle this question. Here we use five coupled climate models and one ocean-only model to evaluate the impact of 0.1 Sv (1 Sv = 10⁶ m³/s) of freshwater equally distributed around the coast of Greenland during the historical era 1965–2004. The ocean-only model helps to discriminate between oceanic and coupled responses. In this idealized framework, we find similar fingerprints in the fourth decade of hosing among the models, with a general weakening of the Atlantic Meridional Overturning Circulation (AMOC). Initially, the additional freshwater spreads along the main currents of the subpolar gyre. Part of the anomaly crosses the Atlantic eastward and enters into the Canary Current constituting a freshwater leakage tapping the subpolar gyre system. As a consequence, we show that the AMOC weakening is smaller if the leakage is larger. We argue that the magnitude of the freshwater leakage is related to the asymmetry between the subpolar-subtropical gyres in the control simulations, which may ultimately be a primary cause for the diversity of AMOC responses to the hosing in the multi-model ensemble. Another important fingerprint concerns a warming in the Nordic Seas in response to the re-emergence of Atlantic subsurface waters capped by the freshwater in the subpolar gyre. This subsurface heat anomaly reaches the Arctic where it emerges and induces a positive upper ocean salinity anomaly by introducing more Atlantic waters. We found similar climatic impacts in all the coupled ocean–atmosphere models with an atmospheric cooling of the North Atlantic except in the region around the Nordic Seas and a slight warming south of the equator in the Atlantic. This meridional gradient of temperature is associated with a southward shift of the tropical rains. The free surface models also show similar sea-level fingerprints notably with a comma-shape of high sea-level rise following the Canary Current.     

A new equation of state for high compression

A new equation of state for metals is developed. The starting point for this development is a realistic expression for the free energy, allowance being made not only for the lattice energy but also for the Fermi energy, the exchange and correlation energies of the electrons. This expression is extended by means of three free parameters. The equation of state is successfully transformed in such a way that only macrophysical material parameters occur in it. In this way, the equation can be applied both to experimental high-pressure physics and to geophysics, its advantage resting in the fact that the bonding structure has been taken into account better than in other equations of state.The new equation, along with two other well-proven equations of state, was tried out on isothermal and Hugoniot data of 40 materials. Good agreement with the measured data was found not only for metallic elements, but also for halides and some oxides of geophysical importance. This agreement even proved to be better for most of the materials than that obtained with the other two equations of state.     

Convection currents in the Earth's mantle and the spherical harmonic development of the topography of the earth

Summary A system of simple principles concerning the kinematics of the convection currents in the Earth's mantle was postulated as the starting point. The mantle is supposed to become more and more homogeneous downward. Here it is essential, for the reason of symmetry, that the current system has a degree of regularity as high as possible and is modified only in the upper part, mainly due to the distribution of continents and oceans. The adequate system is chosen by the help of the theory of point groups. Evaluating the distribution to be expected of the coefficients of the spherical harmonic development of the Earth's theoretical topography by means of electronic computers, accordance is found for the extrema to the 31 st order, i.e. as far as they have been determined until now. — The system yields the seismic discontinuities of the Earth's mantle as a secondary product. That does not mean that the system explains them, but it is absolutely consistent with them. That seems to be aquantitative proof of the system regarded as probable on the strength of group-theoretic considerations.Report No. 74, lectured in Potsdam, May 7, 1970.     

Application of Gibbs Energy Minimization to ModelEarly-Diagenetic Solid-Solution Aqueous-Solution EquilibriaInvolving Authigenic Rhodochrosites in Anoxic BalticSea Sediments

The natural early-diagenetic environment ``anoxicporewater – authigenic mineral phases' has beencharacterized in sediment of the Gotland Deep,Baltic Sea, by a closed-system model. Occurrence ofcarbonate precipitates as thin almost pure whitelaminae was considered as a natural experiment forlong-term equilibration between these phases andporewater. Plots of distribution coefficientsindicate that metastable equilibrium exists betweenporewater and the authigenic Ca-rich rhodochrositephases below 7 cm depth. A thermodynamic model ofporewater geochemistry at in situ P = 25 barand T = 5 °C was developed using the Gibbsenergy minimization (GEM) approach. The values of isobaric-isothermal potentials of Mn, Ca, Fe, Mg,Sr, Ba, C, and O, calculated from the porewatercomposition, were used in a new ``dual thermodynamic'calculation approach to estimate solid activitycoefficients of the end-members in the non-idealsolid solution (Mn, Ca, Mg, Sr, Ba, Fe)CO₃,i.e., at full major and minor multi-componentcomplexity. The regular Margules interactionparameters for the composing binaries estimated bythis model were _Mn-Ca = 1.9 ± 0.5,_Mn-Mg = 0.6,_Ca-Mg = 3.7,_Mn-Fe = 0.2,_Ca-Fe = 2.8,_Mn-Sr = 9.7,_Ca-Sr = 2.15,_Mn-Ba = 4.0,_Ca-Ba = 1.4,validating the theoretical predictions given byLippmann in his pioneering 1980's paper. Thestrictly thermodynamic equilibrium model is not onlyable to match both the measured porewater andcarbonate solid-solution composition, but also topredict that the porewater pH, pe, alkalinity, anddissolved Mn, Fe, and S concentrations arecontrolled by the authigenic mineral bufferingassemblage mackinawite-greigite-rhodochrosite. Ourmodel is only compatible with the idea of ACRformation with typical composition (X_Mnbetween 70–75%) in the topmost sediment layerwhich, however, needs a major source ofMn_aq ^II. This is provided by reduction ofparticulate Mn oxides precipitated in significantamounts in the water column upon major inflow eventsin the Baltic Sea. The model enables also to set upscenarios of changing environmental conditions, e.g.,to predict the non-linear response of the carbonatesolid-solution composition to changes in Mn loading,alkalinity and salinity of the sediment-watersystem. The results suggest that the major andespecially minor element contents (Sr, Mg, Ba) inauthigenic carbonates can be applied as anenvironmental paleoproxy.     

Thermo-hydro-mechanical modeling of carbon dioxide injection for enhanced gas-recovery (CO2-EGR): a benchmarking study for code comparison

The objective of this paper was to investigate the THM-coupled responses of the storage formation and caprock, induced by gas production, CO₂-EGR (enhanced gas recovery), and CO₂-storage. A generic 3D planer model (20,000?×?3,000?×?100?m, consisting of 1,200?m overburden, 100?m caprock, 200?m gas reservoir, and 1,500?m base rock) is adopted for the simulation process using the integrated code TOUGH2/EOS7C-FLAC3D and the multi-purpose simulator OpenGeoSys. Both simulators agree that the CO₂-EGR phase under a balanced injection rate (31,500?tons/year) will cause almost no change in the reservoir pressure. The gas recovery rate increases 1.4?% in the 5-year CO₂-EGR phase, and a better EGR effect could be achieved by increasing the distance between injection and production wells (e.g., 5.83?% for 5?km distance, instead of 1.2?km in this study). Under the considered conditions there is no evidence of plastic deformation and both reservoir and caprock behave elastically at all operation stages. The stress path could be predicted analytically and the results show that the isotropic and extensional stress regime will switch to the compressional stress regime, when the pore pressure rises to a specific level. Both simulators agree regarding modification of the reservoir stress state. With further CO₂-injection tension failure in reservoir could occur, but shear failure will never happen under these conditions. Using TOUGH-FLAC, a scenario case is also analyzed with the assumption that the reservoir is naturally fractured. The specific analysis shows that the maximal storage pressure is 13.6?MPa which is determined by the penetration criterion of the caprock.     

Provenance analyses of the heavy-mineral beach sands of the Annaba coast,northeast Algeria,and their consequences for the evaluation of fossil placer deposit

The paper presents the first study of heavy-mineral sand beaches from the Mediterranean coast of Annaba/Algeria. The studied beaches run along the basement outcrops of the Edough massif, which are mainly composed by micaschists, tourmaline-rich quartzo-feldspathic veins, gneisses, skarns and marbles. Sand samples were taken from three localities (Ain Achir, Plage-Militaire and El Nasr). The heavy-mineral fraction comprises between 74 and 91 vol%. The garnets of the beaches are almandine rich and tourmalines vary with respect to their location from schorl to dravite. Tourmaline at Ain Achir and the Plage-Militaire is schorlits, while at El Nasr beach dravite is ubiquitous. The World Shale Average normalised REE of the sands and the basement outcrops reveal: (i) Ain Achir beach: REE pattern of sand and the coastal rocks from the studied beaches reflects a multiple sources; (ii) Plage-Militaire: the sand and the coastal outcrops show similar LREE and a strong enrichment in HREE, suggesting the presence HREE-rich phases found as inclusions in staurolite; (iii) El Nasr: two types of sand patterns are found: one with flat REE pattern similar to the proximal rocks and other one enriched in HREE suggesting a mixed source.     

Late Quaternary glacial and climate history of the Pamir Mountains derived from cosmogenic Be exposure ages

Moraines southwest of Lake Yashilkul, Pamir, Tajikistan, were dated using ¹⁰Be exposure ages of boulder surfaces. We found evidence for (1) an extensive glaciation ∼60,000 yr ago; (2) a less extensive glacial advance, which deposited a characteristic hummocky moraine lobe with exposure ages ranging from ∼11,000 to 47,000 yr, probably deposited at or before 47,000 yr ago; and (3) lateral moraines with exposure ages of ∼40,000 yr, 27,000 yr and 19,000 yr, respectively. Increasing aridity in the Pamir is most likely responsible for the progressively limited extent of the glaciers during the Late Pleistocene.