Sequentially sampled gas hydrate water,coupled with pore water and bottom water isotopic and ionic signatures at the Kukuy mud volcano,Lake Baikal: ambiguous deep-rooted source of hydrate-forming water

The isotopic and ionic composition of pure gas hydrate (GH) water was examined for GHs recovered in three gravity cores (165–193 cm length) from the Kukuy K-9 mud volcano (MV) in Lake Baikal. A massive GH sample from core St6GC4 (143–165 cm core depth interval) was dissociated progressively over 6 h in a closed glass chamber, and 11 sequentially collected fractions of dissociated GH water analyzed. Their hydrogen and oxygen isotopic compositions, and the concentrations of Cl^– and HCO₃ ^– remained essentially constant over time, except that the fraction collected during the first 50 minutes deviated partly from this pattern. Fraction #1 had a substantially higher Cl^– concentration, similar to that of pore water sampled immediately above (135–142 cm core depth) the main GH-bearing interval in that core. Like the subsequent fractions, however, the HCO₃ ^– concentration was markedly lower than that of pore water. For the GH water fractions #2 to #11, an essentially constant HCO₃ ^–/Cl^– ratio of 305 differed markedly from downcore pore water HCO₃ ^–/Cl^– ratios of 63–99. Evidently, contamination of the extracted GH water by ambient pore water probably adhered to the massive GH sample was satisfactorily restricted to the initial phase of GH dissociation. The hydrogen and oxygen isotopic composition of hydrate-forming water was estimated using the measured isotopic composition of extracted GH water combined with known isotopic fractionation factors between GH and GH-forming water. Estimated δD of ?126 to ?133‰ and δ¹⁸O of ?15.7 to ?16.7‰ differed partly from the corresponding signatures of ambient pore water (δD of ?123‰, δ¹⁸O of ?15.6‰) and of lake bottom water (δD of ?121‰, δ¹⁸O of ?15.8‰) at the St6GC4 coring site, suggesting that the GH was not formed from those waters. Observations of breccias in that core point to a possible deep-rooted water source, consistent with published thermal measurements for the neighboring Kukuy K-2 MV. By contrast, the pore waters of core St6GC4 and also of the neighboring cores GC2 and GC3 from the Kukuy K-9 MV show neither isotopic nor ionic evidence of such a source (e.g., elevated sulfate concentration). These findings constrain GH formation to earlier times, but a deep-rooted source of hydrate-forming water remains ambiguous. A possible long-term dampening of key deep-water source signatures deserves further attention, notably in terms of diffusion and/or advection, as well as anaerobic oxidation of methane.     

Seasonal Cycle of the Near-Surface Diurnal Wind Field Over the Bay of La Paz,Mexico

The results of numerical simulations of the troposphere over the Bay of La Paz, calculated for the months of January, April, July and October during the period 2006–2010 with the Weather Research and Forecast (WRF v3.5) regional model, are used to describe the seasonal features of the diurnal cycle of planetary boundary-layer winds. Two distinct near-surface diurnal flows with strong seasonal variability were identified: (1) a nocturnal and matutinal breeze directed from the subtropical Pacific Ocean, over the Baja California peninsula and the Bay of La Paz, into the Gulf of California that is associated with the regional sea-surface temperature difference between those two major water bodies; and (2) a mid to late afternoon onshore sea-breeze related to the peninsula’s daily cycle of insolation heating that evolves with counter-clockwise rotation over the Bay of La Paz. The model results reveal the interaction over Baja California of opposing afternoon sea-breeze fronts that originate from the subtropical Pacific Ocean and the Gulf of California, with a convergence line forming over the peaks of the peninsula’s topography and the associated presence of a closed vertical circulation cell over the Bay of La Paz and the adjacent Gulf. The collision of the opposing sea-breeze fronts over the narrow peninsula drives convection that is relatively weak due to the reduced heat source and only appears to produce precipitation sporadically. The spatial structure of the sea-breeze fronts over the Bay of La Paz region is complex due to shoreline curvature and nearby topographic features. A comparison of the numerical results with available meteorological near-surface observations indicates that the modelling methodology adequately reproduced the observed features of the seasonal variability of the local planetary boundary-layer diurnal wind cycle and confirms that the low-level atmospheric circulation over the Bay of La Paz is dominated by kinetic energy in the diurnal band. The strongest (weakest) diurnal flows occur during the summer (winter) in response to the seasonally varying magnitudes of the daily land–sea thermal contrast and the regional subtropical Pacific Ocean–Gulf of California sea-surface temperature difference.     

Differential melt scaling for oblique impacts on terrestrial planets

Analytical estimates of melt volumes produced by a given projectile and contained in a given impact crater are derived as a function of impact velocity, impact angle, planetary gravity, target and projectile densities, and specific internal energy of melting. Applications to impact events and impact craters on the Earth, Moon, and Mars are demonstrated and discussed. The most probable oblique impact (45°) produces ～1.6 times less melt volume than a vertical impact, and ～1.6 and 3.7 times more melt volume than impacts with 30° and 15° trajectories, respectively. The melt volume for a particular crater diameter increases with planetary gravity, so a crater on Earth should have more melt than similar-size craters on Mars and the Moon. The melt volume for a particular projectile diameter does not depend on gravity, but has a strong dependence on impact velocity, so the melt generated by a given projectile on the Moon is significantly larger than on Mars. Higher surface temperatures and geothermal gradients increase melt production, as do lower energies of melting. Collectively, the results imply thinner central melt sheets and a smaller proportion of melt particles in impact breccias on the Moon and Mars than on Earth. These effects are illustrated in a comparison of the Chicxulub crater on Earth, linked to the Cretaceous–Tertiary mass extinction, Gusev crater on Mars, where the Mars Exploration Rover Spirit landed, and Tsiolkovsky crater on the Moon. The results are comparable to those obtained from field and spacecraft observations, other analytical expressions, and hydrocode simulations.     

Interface and embedded joint methods for modeling dynamic fracture opening by explosive products

Two computational approaches are proposed in the paper to model dynamic fracture opening by explosive products. The first method assumes that the fractures may be modeled using flow elements embedded along the mesh lines. This method models crack opening in a straightforward way by splitting the nodes of the computational grid. It can account for crack branching; however, the crack directions are constrained by existing mesh faces, which may lead to mesh dependence. Also, the stress in flow elements is calculated explicitly separate from the surrounding solid elements that can impose additional limits on the time step stability condition for explicit integration. The second approach uses embedded flow elements to model the cracks. Typical thickness of the cracks is much smaller than the element size. Therefore, gas pressure in the cracks is assumed to be in stress equilibrium with the element stress. To achieve this, the crack thickness and the state of the gas is updated simultaneously with the state of the solid element which contains the crack. Therefore, the time step is controlled by the explicit solver applied for the solid and does not depend on the thickness of the crack. The main disadvantage of the second approach is due to the complexity of modeling multiple intersecting cracks, which go through the same element. We discuss the areas of possible applications of these 2 methods and the ways to improve and enhance them for future practical applications.     

A semi-analytical approach to a nonlinear stress–strain analysis of buried steel pipelines crossing active faults

In the present paper a semi-analytical methodology for a nonlinear stress–strain analysis of buried steel pipelines at active fault crossings is presented and verified. The developed model introduces a number of critical refinements to the existing methodologies which extend the application area of analytical models in pipeline design. In particular, a strike-slip and normal-slip fault crossings can be analyzed taking into account material and large displacement nonlinearities, nonlinear pipe–soil interaction. The proposed model is verified against the results by other authors and numerical results, obtained with the finite element method.     

Intense diurnal surface currents in the Bay of La Paz,Mexico

Currents in the northern Bay of La Paz were examined using an 8-month Acoustic Doppler Current Profiler (ADCP) record collected in the upper 185 m of the water column during 2007. Flow variability was dominated by tidal motions, which accounted for 43% (33% diurnal, 10% semidiurnal) of the total kinetic energy. The tidal motions had a pronounced vertical structure dominated within a shallow (∼30 m thick) surface layer by intense counterclockwise (CCW) rotary S₁ diurnal radiational currents that were highly coherent with the counterclockwise seabreeze. Motions within the semidiurnal frequency band were primarily associated with significant counterclockwise S₂ radiational tidal currents, which were also coherent with the seabreeze. Both S₁ and S₂ tidal ellipses in the upper layer were aligned perpendicular to the bay entrance with mean semi-major axes of 55 and 20 cm/s, respectively. Below the surface layer, tidal currents decayed rapidly to relatively weak, clockwise rotary barotropic motions. In contrast to those for radiational harmonics, tidal ellipses of the gravitational constituents (M₂, K₁ and O₁) were oriented cross-bay. Energy within the diurnal frequency band in the surface layer was dominated by a coherent component (barotropic, phase-locked baroclinic and radiational), which accounted for roughly 65% (59% from S₁ alone) of the total diurnal kinetic energy. Of the remaining diurnal band energy, 18% was associated with an incoherent baroclinic component and 17% with a background noise component. Below 30 m depth, the corresponding estimates are 40%, 32% and 28%, respectively. The persistent, surface-intensified CCW rotary currents observed at the mooring site are assumed to be forced by strong CCW seabreeze winds in the presence of a “slippery” low-density surface layer. This response may be further augmented by topographic narrowing at the bay entrance and by the close proximity of the diurnal and inertial frequency bands in the region.     

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.