Evapotranspiration capture and stream depletion due to groundwater pumping under variable boreal climate conditions: Sudogda River Basin,Russia

Groundwater pumping and changes in climate-induced recharge lead to lower groundwater levels and significant changes in the water balance of a catchment. Water previously discharged as evapotranspiration can become a source of pumpage. Neglecting this effect leads to overestimated streamflow depletion. A small river basin (Sudogda River Basin, Russia) with a boreal climate and with long-term records of groundwater head and streamflow rate (showing that the measured stream depletion is less than the pumping rate) was investigated. The role of evapotranspiration in the water balance was analyzed by a hydrogeological model using MODFLOW-2005 with the STR package; the annual variation in recharge was obtained with the codes Surfbal and HYDRUS. The Sudogda River Basin was classified according to landscape and unsaturated-zone texture classes, and for each classified zone, the unsaturated-zone flow simulation was used to calculate the annual recharge dynamics for the observation period. Calibration of the regional flow model was conducted using flow and head observations jointly for two steady-state flow conditions—natural (before pumping started) and stressed (pumping). The simulations showed that pumped water originates from three sources: intercepted baseflow (75% of the annual total pumping rate), the capture of groundwater evapotranspiration discharge plus increased groundwater recharge (17%), and induced stream infiltration (8%). Additionally, multi-year precipitation records were analyzed to detect any long-term recharge and pumping water-budget changes. The results showed that increasing groundwater recharge by natural precipitation leads to (1) decreased intercepted baseflow and induced streamflow infiltration and (2) increased intercepted evapotranspiration discharge, thereby reducing stream depletion.     

Spatiotemporal dynamics: the need for an innovative approach in mountain hazard risk management

Starting with an overview on losses due to mountain hazards in the Russian Federation and the European Alps, the question is raised why a substantial number of events still are recorded—despite considerable efforts in hazard mitigation and risk reduction. The main reason for this paradox lies in a missing dynamic risk-based approach, and it is shown that these dynamics have different roots: firstly, neglecting climate change and systems dynamics, the development of hazard scenarios is based on the static approach of design events. Secondly, due to economic development and population dynamics, the elements at risk exposed are subject to spatial and temporal changes. These issues are discussed with respect to temporal and spatial demands. As a result, it is shown how risk is dynamic on a long-term and short-term scale, which has to be acknowledged in the risk concept if this concept is targeted at a sustainable development of mountain regions. A conceptual model is presented that can be used for dynamical risk assessment, and it is shown by different management strategies how this model may be converted into practice. Furthermore, the interconnectedness and interaction between hazard and risk are addressed in order to enhance prevention, the level of protection and the degree of preparedness.     

Elastic wave‐vector decomposition in heterogeneous anisotropic media

The goal of wave‐mode separation and wave‐vector decomposition is to separate a full elastic wavefield into three wavefields with each corresponding to a different wave mode. This allows elastic reverse‐time migration to handle each wave mode independently. Several of the previously proposed methods to accomplish this task require the knowledge of the polarisation vectors of all three wave modes in a given anisotropic medium. We propose a wave‐vector decomposition method where the wavefield is decomposed in the wavenumber domain via the analytical decomposition operator with improved computational efficiency using low‐rank approximations. The method is applicable for general heterogeneous anisotropic media. To apply the proposed method in low‐symmetry anisotropic media such as orthorhombic, monoclinic, and triclinic, we define the two S modes by sorting them based on their phase velocities (S1 and S2), which are defined everywhere except at the singularities. The singularities can be located using an analytical condition derived from the exact phase‐velocity expressions for S waves. This condition defines a weight function, which can be applied to attenuate the planar artefacts caused by the local discontinuity of polarisation vectors at the singularities. The amplitude information lost because of weighting can be recovered using the technique of local signal–noise orthogonalisation. Numerical examples show that the proposed approach provides an effective decomposition method for all wave modes in heterogeneous, strongly anisotropic media.     

Dynamics of Global Climatic Changes and Possibility of Their Prediction Using Chemical Data from Greenland Ice

Records from the Greenland Ice Sheet Project Two ice core are considered in terms of dynamical systems theory and nonlinear prediction. Dynamical systems theory allows us to reconstruct some properties of a phenomenon based only on past behavior without any mechanistic assumptions or deterministic models. A near-term prediction of temperature, including a mean estimate and confidence interval, is made for 800 years into the future. The prediction suggests that the present short-time global warming trend will continue for at least 200 years and be followed by a reverse in the temperature trend.     

Mesoscale circulation along the Sakhalin Island eastern coast

The seasonal and interannual variability of mesoscale circulation along the eastern coast of the Sakhalin Island in the Okhotsk Sea is investigated using the AVISO velocity field and oceanographic data for the period from 1993 to 2016. It is found that mesoscale cyclones with the horizontal dimension of about 100 km occur there predominantly during summer, whereas anticyclones occur predominantly during fall and winter. The cyclones are generated due to a coastal upwelling forced by northward winds and the positive wind stress curl along the Sakhalin coast. The anticyclones are formed due to an inflow of low-salinity Amur River waters from the Sakhalin Gulf intensified by southward winds and the negative wind stress curl in the cold season. The mesoscale cyclones support the high biological productivity at the eastern Sakhalin shelf in July– August.     

Automated placement of supplementary contour lines

Supplementary contour lines are placed between regular contour lines to visualize small but important forms that regular contour lines are unable to show. On topographic maps, typical forms are hillcrests, depressions, saddles, terraces, banks, and levees. No automated method for the selection of supplementary contour lines has been described so far. We document cartographic design principles for the selection of supplementary contour lines for topographic maps, and present an automated method for their placement. Results of the automated method are similar to manually placed supplementary contour lines. Our method helps map authors to create contour line maps that more effectively illustrate relevant small details in maps showing terrain elevation or other scalar fields.     

37 ka以来日本海沉积物有机质碳和氮稳定同位素变化及其古海洋学意义

海洋沉积物有机质碳氮稳定同位素(δ¹³C、δ¹⁵N)广泛用于有机质来源示踪、古生产力和古海洋环境重建。日本海沉积物δ¹³C和δ¹⁵N值一个显著特征是在末次冰盛期(LGM)同步负偏,但是对这一现象产生的原因以及他们的演化过程的认识仍然存在明显不足。在本研究中,我们详细调查了37 ka以来日本海中部LV53-23-1岩心沉积物δ¹³C和δ¹⁵N演化历史。结果显示,沉积物δ¹³C和δ¹⁵N分别介于-26.3‰至-22.5‰和1.6‰至6.1‰,低值出现在LGM(26.5~17 ka)暗色层状泥发育时期,指示较强的陆源输入贡献。在Heinrich冰阶1时期(17~14.5 ka),δ¹³C和δ¹⁵N快速正偏,表明日本海海洋环境发生了明显的转换,对应于对马海峡淹没及对马暖流入侵。14.5 ka之后,沉积物δ¹⁵N值恢复到5‰,与开阔大洋海水硝酸盐的δ¹⁵N值近似。我们采用二端员混合模型粗略地估算了有机质来源的相对贡献。LGM时期陆源有机质贡献介于65%至80%,14.5 ka以后海源有机质贡献介于60%至80%。除了增加的陆源有机质贡献以外,LGM时期沉积物δ¹⁵N亏损还涉及如下过程:(1)较高的含Fe沙尘供给提高日本海表层海洋生物固氮效率;(2)缺氧环境盛行减弱成岩作用对沉积物δ¹⁵N影响。37 ka以来,日本海沉积物δ¹³C和δ¹⁵N变化与有机质来源、营养盐的供给、表层生产力和沉积物氧化还原条件相关,实际受海平面和全球气候制约。     

Effect of pH and organic ligands on the kinetics of smectite dissolution at 25 °C

Forward dissolution rates of Na-Montmorillonite (Wyoming) SWy-2 smectite (Ca_0.06Na_0.56)[Al_3.08Fe(III)_0.38Mg_0.54] [Si_7.93 Al_0.07]O₂₀(OH)₄ were measured at 25 °C in a mixed-flow reactor equipped with interior dialysis compartment (6-8 kDa membrane) as a function of pH (1-12), dissolved carbonate (0.5-10 mM), phosphate (10⁻⁵ to 0.03 M), and nine organic ligands (acetate, oxalate, citrate, EDTA, alginate, glucuronic acid, 3,4-dihydroxybenzoic acid, gluconate, and glucosamine) in the concentration range from 10⁻⁵ to 0.03 M. In organic-free solutions, the Si-based rates decrease with increasing pH at 1 ? pH ? 8 with a slope close to −0.2. At 9 ? pH ? 12, the Si-based rates increase with a slope of ∼0.3. In contrast, non-stoichiometric Mg release weakly depends on pH at 1 ? pH ? 12 and decreases with increasing pH. The empirical expression describing Si-release rates [R, mol/cm²/s] obtained in the present study at 25 °C, I = 0.01 M is given by