A reexamination of the valley wind system in the Alpine Inn Valley with numerical simulations

Summary This paper presents idealized numerical simulations of the valley wind circulation in the Alpine Inn Valley, which are compared with existing data and are used to improve our dynamical understanding of the valley wind. The simulations have been performed with the Penn State/NCAR mesoscale model MM5. They use a high-resolution realistic topography but idealized large-scale conditions without any synoptic forcing to focus on the thermally induced valley wind system. The comparison with the available observations shows that this simplified set-up is sufficient to reproduce the essential features of the valley wind.The results show that the tributaries of the Inn Valley have a considerable impact on the along-valley mass fluxes associated with the valley wind circulation. The upvalley mass flux is found to increase where tributaries enter the Inn Valley from the north, that is, from the direction where the Alpine foreland is located. On the other hand, the upvalley mass flux is reduced at the junctions with southern tributaries because part of the upvalley flow is deflected into these tributaries. For the downvalley flow, the situation is essentially reversed, but the influence of the valley geometry on the flow structure is larger than for the upvalley flow. The most important feature is a lateral valley contraction near the valley exit into the Alpine foreland. It reduces the downvalley mass flux at low levels, so that the wind maximum in the interior of the valley is shifted to a fairly large distance from the ground. North of the valley contraction, however, the downvalley flow strongly accelerates and forms a pronounced low-level jet. A dynamical analysis indicates that this acceleration can be interpreted as a transition from subcritical to supercritical hydraulic flow. Another interesting feature is that the low-level jet maintains its structure for several tenths of kilometres into the Alpine foreland. This appears to be related to the fact that the lateral wind shear on the flanks of the jet is associated with a strong dipole of potential vorticity (PV). Due to the conservation properties of the PV, the downstream advection of the PV dipole leads to the formation of a band-like feature that decays fairly slowly.     

Russia at GHG Market

In the first Kyoto commitment period Russia could be the major supplier for the greenhouse gases (GHG) emissions market. Potential Russian supply depends on the ability of Russia to keep GHG emissions lower than the Kyoto target. In the literature there is no common understanding of the total trading potential of Russia at the international carbon market. In this paper we focus on CO₂ emission, which constituted nearly 80%of Russian GHG emission. We compare different projections of Russian CO₂emission and analyze the most important factors, which predetermine the CO₂emission growth. In a transition economy these factors are: Gross Domestic Product(GDP) dynamic, changes of GDP structure, innovation activity, transformation of export-import flows and response to the market signals. The input-output macroeconomic model with the two different input-output tables representing old and new production technologies has been applied for the analysis to simulate technological innovations and structural changes in the Russian economy during transition period. The Russian supply at the international GHG market without forest sector may be up to 3 billion metric ton of CO₂ equivalent. Earlier actions to reduce CO₂ emission are critical to insure theRussiansupply at the international carbon market. With regard to the current status of the Russian capital market, the forward trading with OECD countries is only the possibility to raise initial investments to roll no-regret and low-cost GHG reduction. This paper discusses uncertainties of RussianCO₂emission dynamics and analyzes the different incentives to lower the emission pathway.     

Revising macroseismic data in Switzerland: The December 20, 1720 earthquake in the region of Lake Constance

Upgrading the Earthquake Catalog of Switzerland (ECOS) included revising the earthquake of 1720. This change has major importance for history and seismology.Although that quake has been the subject of several publications, none was based on critical methods. This re-evaluation of the event is built upon a new and more reliable database established after investigating archives and libraries. Using data from such historical sources, we assigned new site intensities, adopting the criteria established by the European macroseismic scale EMS 98 (Grünthal, 1998).We discovered that the event had been assigned an overestimated intensity, due to interpretation errors in former earthquake catalogs and compilations. We recommend reducing the intensity from I₀= VIII to I₀= VI (EMS 98). The moment magnitude is given as M_W= 4.6. Since the event had been considered the largest for its respective area, downgrading it now will influence the seismic hazard assessment for this region.     

Ambiguous Moment Tensors and Radiation Patterns in Anisotropic Media with Applications to the Modeling of Earthquake Mechanisms in W-Bohemia

Anisotropic material properties are usually neglected during inversions for source parameters of earthquakes. In general anisotropic media, however, moment tensors for pure-shear sources can exhibit significant non-double-couple components. Such effects may be erroneously interpreted as an indication for volumetric changes at the source. Here we investigate effects of anisotropy on seismic moment tensors and radiation patterns for pure-shear and tensile-type sources. Anisotropy can significantly influence the interpretation of the source mechanisms. For example, the orientation of the slip within the fault plane may affect the total seismic moment. Also, moment tensors due to pure-shear and tensile faulting can have similar characteristics depending on the orientation of the elastic tensor. Furthermore, the tensile nature of an earthquake can be obscured by near-source anisotropic properties. As an application, we consider effects of inhomogeneous anisotropic properties on the seismic moment tensor and the radiation patterns of a selected type of micro-earthquakes observed in W-Bohemia. The combined effects of near-source and along-path anisotropy cause characteristic amplitude distortions of the P, S1 and S2 waves. However, the modeling suggests that neither homogeneous nor inhomogeneous anisotropic properties alone can explain the observed large non-double-couple components.The results also indicate that a correct analysis of the source mechanism, in principle, is achievable by application of anisotropic moment tensor inversion.     

Mineral dissolution in the Cape Cod aquifer, Massachusetts, USA: I . Reaction stoichiometry and impact of accessory feldspar and glauconite on strontium isotopes, solute concentrations, and REY distribution

To compare relative reaction rates of mineral dissolution in a mineralogically simple groundwater aquifer, we studied the controls on solute concentrations, Sr isotopes, and rare earth element and yttrium (REY) systematics in the Cape Cod aquifer. This aquifer comprises mostly carbonate-free Pleistocene sediments that are about 90% quartz with minor K-feldspar, plagioclase, glauconite, and Fe-oxides. Silica concentrations and pH in the groundwater increase systematically with increasing depth, while Sr isotopic ratios decrease. No clear relationship between ⁸⁷Sr/⁸⁶Sr and Sr concentration is observed. At all depths, the ⁸⁷Sr/⁸⁶Sr ratio of the groundwater is considerably lower than the Sr isotopic ratio of the bulk sediment or its K-feldspar component, but similar to that of a plagioclase-rich accessory separate obtained from the sediment. The Si-⁸⁷Sr/⁸⁶Sr-depth relationships are consistent with dissolution of accessory plagioclase. In addition, solutes such as Sr, Ca, and particularly K show concentration spikes superimposed on their respective general trends. The K-Sr-⁸⁷Sr/⁸⁶Sr systematics suggests that accessory glauconite is another major solute source to Cape Cod groundwater. Although the authigenic glauconite in the Cape Cod sediment is rich in Rb, it is low in in-grown radiogenic ⁸⁷Sr because of its young Pleistocene age. The low ⁸⁷Sr/⁸⁶Sr ratios are consistent with equilibration of glauconite with seawater. The impact of glauconite is inferred to vary due to its variable abundance in the sediments. In the Cape Cod groundwater, the variation of REY concentrations with sampling depth resembles that of K and Rb, but differs from that of Ca and Sr. Shale-normalized REY patterns are light REY depleted, show negative Ce anomalies and super-chondritic Y/Ho ratios, but no Eu anomalies. REY input from feldspar, therefore, is insignificant compared to input from a K-Rb-bearing phase, inferred to be glauconite. These results emphasize that interpretation of groundwater chemistry, even in relatively simple aquifers, may be complicated by solute contributions from “exotic” accessory minerals such as glauconite. To detect such peculiarities, groundwater studies should combine the study of elemental concentration and isotopic composition of several solutes that show different geochemical behavior.     

The effect of dissolved barium on biogeochemical processes at cold seeps

A numerical model was applied to investigate and quantify the biogeochemical processes fueled by the expulsion of barium and methane-rich fluids in the sediments of a giant cold-seep area in the Derugin Basin (Sea of Okhotsk). Geochemical profiles of dissolved Ba²⁺, Sr²⁺, Ca²⁺, SO₄²⁻, HS⁻, DIC, I⁻ and of calcium carbonate (CaCO₃) were fitted numerically to constrain the transport processes and the kinetics of biogeochemical reactions. The model results indicate that the anaerobic oxidation of methane (AOM) is the major process proceeding at a depth-integrated rate of 4.9 μmol cm⁻² a⁻¹, followed by calcium carbonate and strontian barite precipitation/dissolution processes having a total depth-integrated rate of 2.1 μmol cm⁻² a⁻¹. At the low seepage rate prevailing at our study site (0.14 cm a⁻¹) all of the rising barium is consumed by precipitation of barite in the sedimentary column and no benthic barium flux is produced. Numerical experiments were run to investigate the response of this diagenetic environment to variations of hydrological and biogeochemical conditions. Our results show that relatively low rates of fluid flow (<∼5 cm a⁻¹) promote the dispersed precipitation of up to 26 wt% of barite and calcium carbonate throughout the uppermost few meters of the sedimentary column. Distinct and persistent events (several hundreds of years long) of more vigorous fluid flow (from 20-110 cm a⁻¹), instead, result in the formation of barite-carbonate crusts near the sediment surface. Competition between barium and methane for sulfate controls the mineralogy of these sediment precipitates such that at low dissolved methane/barium ratios (<4-11) barite precipitation dominates, while at higher methane/barium ratios sulfate availability is limited by AOM and calcium carbonate prevails. When seepage rates exceed 110 cm a⁻¹, barite precipitation occurs at the seafloor and is so rapid that barite chimneys form in the water column. In the Derugin Basin, spectacular barite constructions up to 20 m high, which cover an area of roughly 22 km² and contain in excess of 5 million tons of barite, are built through this process. In these conditions, our model calculates a flux of barium to the water column of at least 20 μmol cm⁻² a⁻¹. We estimate that a minimum of 0.44 × 10⁶ mol a⁻¹ are added to the bottom waters of the Derugin Basin by cold seep processes, likely affecting the barium cycle in the Sea of Okhotsk.     

Tectonic evolution of the Intra-Pontide suture zone in the Armutlu Peninsula, NW Turkey

The Armutlu Peninsula and adjacent areas in NW Turkey play a critical role in tectonic reconstructions of the southern margin of Eurasia in NW Turkey. This region includes an inferred Intra-Pontide oceanic basin that rifted from Eurasia in Early Mesozoic time and closed by Late Cretaceous time. The Armutlu Peninsula is divisible into two metamorphic units. The first, the Armutlu Metamorphics, comprises a ?Precambrian high-grade metamorphic basement, unconformably overlain by a ?Palaeozoic low-grade, mixed siliciclastic/carbonate/volcanogenic succession, including bimodal volcanics of inferred extensional origin, with a possibly inherited subduction signature. The second unit, the low-grade znik Metamorphics, is interpreted as a Triassic rift infilled with terrigenous, calcareous and volcanogenic lithologies, including basalts of within-plate type. The Triassic rift was unconformably overlain by a subsiding Jurassic–Late Cretaceous (Cenomanian) passive margin including siliciclastic/carbonate turbidites, radiolarian cherts and manganese deposits. The margin later collapsed to form a flexural foredeep associated with the emplacement of ophiolitic rocks in Turonian time. Geochemical evidence from meta-basalt blocks within ophiolite-derived melange suggests a supra-subduction zone origin for the ophiolite. The above major tectonic units of the Armutlu Peninsula were sealed by a Maastrichtian unconformity. Comparative evidence comes from the separate Almacık Flake further east.Considering alternatives, it is concluded that a Mesozoic Intra-Pontide oceanic basin separated Eurasia from a Sakarya microcontinent, with a wider Northern Neotethys to the south. Lateral displacement of exotic terranes along the south-Eurasian continental margin probably also played a role, e.g. during Late Cretaceous suturing, in addition to overthrusting.     

Three-dimensional forward and backward modelling of diapirism: numerical approach and its applicability to the evolution of salt structures in the Pricaspian basin

Numerical studies of ductile deformations induced by salt movements have, until now, been restricted to two-dimensional (2D) modelling of diapirism. This paper suggests a numerical approach to model the evolution of three-dimensional (3D) salt structures toward increasing maturity. This approach is also used here to restore the evolution of salt structures through successive earlier stages. The numerical methodology is applied to study several model examples. We analyse a model of salt diapirs that develop from an initial random perturbation of the interface between salt and its overburden and restore the evolved salt diapirs to their initial stages. We show that the average restoration errors are less than 1%. An evolutionary model of a 2D salt wall loaded by a 2D pile of sediments predicts a decomposition of the salt wall into 3D diapiric structures when the overburden of salt is supplied by 3D synkinematic wedge of sediments. We model salt extrusion feeding a gravity current over the depositional surface and estimate an average rate of extrusion and horizontal velocity of salt spreading. Faulting of the overburden to salt overhangs initiates new secondary diapirs, and we analyse the growth of these secondary diapirs. We also study how lateral flow effects the evolution of salt diapirs. The shape of a salt diapir can be very different if the rate of horizontal flow is much greater than the initial rate of diapiric growth solely due to gravity. We discuss the applicability of the results of the models to the evolution of Late Permian salt structures in the Pricaspian basin (Russia and Kazakhstan). These structures are distinguishable into a variety of styles representing different stages of growth: salt pillows, diapirs, giant salt massifs, 2D diapiric walls and 3D stocks complicated by large overhangs. The different sizes, shapes and maturities of salt structures in different parts of the Pricaspian basin reflect areal differences in salt thickness and loading history. Our results suggest that the numerical methodology can be employed to analyse the evolution of all salt structures that have upbuilt through younger ductile overburdens.