Distribution and morphology of mud volcanoes and other fluid flow-related lake-bed structures in Lake Baikal, Russia

New high-resolution multibeam bathymetry data recorded in 2009 in the deepest lake in the World, Lake Baikal, Siberia, enabled a better understanding of the morphology of ten known lake-bed structures—the Bolshoy, Malenki, Malyutka and Stari mud volcanoes in the South Baikal Basin, the K1–4 structures in the Selenga delta, and the Novosibirsk and St. Petersburg structures in the Central Baikal Basin—and also the discovery of 29 new lake-bed structures. These new structures are the S1, Tolstiy, mTSG and S2 in the South Baikal Basin, the P1–P4, P6–P19 and K5–K8 in the Selenga delta accommodation zone, and the C1, C3 and C4 edifices in the Central Baikal Basin. In all, 39 positive relief structures were identified and their large-scale distribution mapped. Based on their typical shape, the observation of high-reflectivity areas on side-scan sonar data records, and evidence of feeder channels on subsurface data, these structures can be classified as mud volcanoes. This has already been confirmed in other publications for the Bolshoy, Malenki and K2 structures, by the recovery of mud breccias in sediment cores. Most structures occur on or near faults and have orientations parallel with the major faults and main stress orientations in the basins, suggesting a strong structural control on the formation of the mud volcanoes. Their slopes are generally steeper than 5°, consistent with interpretation as mud cones formed by high-viscosity, stiff mud plugs. Only few structures appear to be characterised by a crater, in which case this apparent crater seems to be formed by the coalescence of several single cones, leaving a depression in the centre. Some structures have a moat, which has probably an erosional origin. Furthermore, three depressions have been found, named P5, P20 and C2, which are suggested to be pockmarks.     

Adaptive volume penalization for ocean modeling

The development of various volume penalization techniques for use in modeling topographical features in the ocean is the focus of this paper. Due to the complicated geometry inherent in ocean boundaries, the stair-step representation used in the majority of current global ocean circulation models causes accuracy and numerical stability problems. Brinkman penalization is the basis for the methods developed here and is a numerical technique used to enforce no-slip boundary conditions through the addition of a term to the governing equations. The second aspect to this proposed approach is that all governing equations are solved on a nonuniform, adaptive grid through the use of the adaptive wavelet collocation method. This method solves the governing equations on temporally and spatially varying meshes, which allows higher effective resolution to be obtained with less computational cost. When penalization methods are coupled with the adaptive wavelet collocation method, the flow near the boundary can be well-resolved. It is especially useful for simulations of boundary currents and tsunamis, where flow near the boundary is important. This paper will give a thorough analysis of these methods applied to the shallow water equations, as well as some preliminary work applying these methods to volume penalization for bathymetry representation for use in either the nonhydrostatic or hydrostatic primitive equations.     

Surface chemistry and reactivity of plant phytoliths in aqueous solutions

In order to better understand the reactivity of plant phytoliths in soil solutions, we determined the solubility, surface properties (electrophoretic mobilities and surface charge) and dissolution kinetics of phytoliths extracted from fresh biomass of representative plant species (larch tree and elm, horsetail, fern, and four grasses) containing significant amount of biogenic silica. The solubility product of larch, horsetail, elm and fern phytoliths is close to that of amorphous silica and soil bamboo phytoliths. Electrophoretic measurements yield isoelectric point pH_IEP = 0.9, 1.1, 2.0 and 2.2 for four grasses, elm, larch and horsetail phytoliths respectively, which is very close to that of quartz or amorphous silica. Surface acid–base titrations allowed generation of a 2-pK surface complexation model (SCM) for larch, elm and horsetail phytoliths. Phytoliths dissolution rates, measured in mixed-flow reactors at far from equilibrium conditions at 1 ≤ pH ≤ 8, were found to be very similar among the species, and close to those of soil bamboo phytoliths. Mechanistic treatment of all plant phytoliths dissolution rates provided three-parameters equation sufficient to describe phytoliths reactivity in aqueous solutions:$R(\text{mol}/{\text{cm}}^2/\text{s})=6?{10}^{?16}?a_{\text{H+}}+5.0?{10}^{?18}+3.5?{10}^{?13}?a_{\text{OH}?}^{0.33}$Alternatively, the dissolution rate dependence on pH can be modeled within the concept of surface coordination theory assuming the rate proportional to concentration of > SiOH₂⁺, > SiOH⁰ and > SiO^? species. In the range of Al concentration from 20 to 5000 ppm in the phytoliths, we have not observed any correlation between their Al content and solubility, surface acid–base properties and dissolution kinetics.It follows from the results of this study that phytoliths dissolution rates exhibit a minimum at pH ～ 3. Mass-normalized dissolution rates are similar among all four types of plant species studied and these rates are an order of magnitude higher than those of typical soil clay minerals. The minimal half life time of larch and horsetail phytoliths in the interstitial soil solution ranges from 10–12 years at pH = 2–3 to < 1 year at pH above 6, comparable with mean residence time of phytoliths in soil from natural observations.     

Numerical modeling of endogenic thermal anomalies on Europa

A variety of recent resurfacing features have been observed on Europa, which may produce thermal anomalies detectable by a future mission. However, the likelihood of such a detection depends on their size and lifetimes. The results of this numerical study suggest that the lifetime of a thermal anomaly associated with the emplacement of 100 m of water onto the surface of Europa is several hundred years, and ∼10 years for 10 m of water. If warm ice is emplaced on the surface instead of liquid water, these lifetimes decrease by up to a factor of two. Exploration of model parameters indicates that a thin insulating surface layer can double thermal anomaly lifetimes, anomalies emplaced at a latitude of 80° can remain detectable nearly a factor of two longer than those at equatorial latitudes, and anomalies on the night side can remain detectable for up to ∼20% longer than those on the day side. High temperatures are very short-lived as the surface ice cools very rapidly to below 200 K due to sublimation cooling. Assuming steady-state resurfacing, the number of detectable thermal anomalies associated with the emplacement of 100 m of water would be on the order of 10 if the typical resurfacing area is 15 km². If recent resurfacing is dominated by chaos regions with typical areas of 100 to 1000 km² and lifetimes of 1000 to 4000 years, the number of detectable thermal anomalies would be on the order of 1 to 10.     

On tourmaline

The first part of this paper presents chemical analyses and physical properties for 12 tourmalines of varying color, 10 samples from Mozambique, one from Afghanistan and one from Madagascar. The data are summarized in Tables 1–3. The second part reviews the isomorphism in the tourmaline group on the basis of data gathered from the literature. Plots of c₀ against a₀ for 198 X-rayed tourmalines and of (Fe, Mn, Mg, Ti) against the unit cell volume of 57 analyzed and X-rayed tourmalines supplement Epprecht's review.     

A hydrologic contribution to risk assessment for the Caspian Sea

The Caspian Sea (CS), the world's largest inland sea, may also be considered as large-scale limnic system. Due to strong fluctuations of its water level during the 20th century and the flooding of vast areas in a highly vulnerable coastal zone, economic and environmental risk potentials have to be considered. Since the major water input into the CS is attributed to the Volga river, the understanding of its long-term flow process is necessary for an appropriate risk assessment for the CS and its coastal area. Therefore, a top-down approach based on statistical analyses of long-term Volga flow series is pursued. For the series of annual mean flow (MQ) of the Volga river basin during the 20th century, a complex oscillation pattern was identified. Analyses for multiple gauges in the Volga river basin and Eurasian reference basins revealed that this oscillation pattern resulted from the superposition of oscillations with periods of ∼30 years (MQ) in the western part of the Volga river basin, and ∼14 years (flow volume of snowmelt events) and ∼20 years (flow volume of summer and autumn) in the eastern part of the Volga river basin (Kama river basin). Almost synchronous minima or maxima of these oscillations occurred just in the periods of substantial changes of the Caspian Sea level (CSL). It can thus be assumed that the described mechanism is fundamental for an understanding of the CSL development during the 20th century. Regarding the global climate change, it is still difficult to predict reliably the development of the CSL for the 21st century. Consequently, we suggest an ongoing, interdisciplinary research co-operation among climatology, hydrology, hydraulics, ecology and spatial data management.     

Theory of EM monitoring of sea bottom geothermal areas

Monitoring of geophysical conditions of marine sedimentary basins is necessary for predicting seismic events and for adaptation of geothermal technologies for seismically active (as a rule) sea bottom geothermal areas. These conditions are characterized by seismo-hydro-electromagnetic (EM) geophysical field interaction in the presence of gravity. Based on the main physical principles, geophysical and petrophysical data, we formulate a mathematical model of seismo-hydro-EM interaction in a basin of a marginal sea and calculate the transformation of a seismic excitation in the upper mantle under the central part of the sea of Japan into the low-frequency (0.1 to 10 Hz) EM signals at the top of the sea bottom sedimentary layer, at the sea surface and in the atmosphere up to the lower boundary of the ionosphere. Physics of the EM generation and propagation process is shown including: generation of EM waves in the upper mantle layer M by a seismic wave from under M, spatial modulation of diffusive EM waves by a seismic wave, stopping of the EM wave arrived (before the seismic P wave) from the upper mantle at the top of the sediments because of the high electric conductivity of seawater (3.5 S/m), immediate penetration of the EM wave through the seawater thickness after the delayed seismic P wave shock into the sea bottom, and EM emission from the sea surface into the atmosphere. Let us note that the EM signal in the sea bottom sediments is the first measurable signal of a seismic activation of geological structures beneath the seafloor and this signal is protected by seawater from the influence ionosphere disturbances. Amplitude of the computed magnetic signals (300, 200, 50, and 30 pT at the ocean–atmosphere interface and at the height of 10, 30 and 50 km, respectively), their predominant frequency (0.25 Hz), the delay of the seismic P wave in regard to the magnetic signal for the receivers at the shore (20 s), the amplitude of temperature disturbances in sediments (up to 0.02 K), the parameters of the long (150 km) tsunami wave of a small (up to 20 cm) amplitude far from the shore and other values that characterize the seismo-hydro-EM process are of the orders observed. Recommendations for the EM monitoring of dynamic processes beneath seafloor geothermal areas are given.     

Instrumental Neutron and Photon Activation Analysis in the Geochemical Study of Phonolitic and Trachytic Rocks

More than 130 samples of Tertiary phonolitic, trachytic and trachyandesitic rocks, representing evolved members of strongly and weakly alkaline volcanic rock series of the Ohre (Eger) Rift in the Bohemian Massif, Czech Republic, were analysed by instrumental neutron and photon activation analysis (INAA and IPAA, respectively). Forty-two major and trace elements were determined to explain the origin of such evolved rocks. Specific features of INAA and IPAA of the silicic, highly alkaline rocks are discussed, and results for selected samples are given.