Temporal and spatial variations of <Emphasis Type="Italic">δ</Emphasis><Superscript>15</Superscript>N and <Emphasis Type="Italic">δ</Emphasis><Superscript>18</Superscript>O for atmospheric N<Subscript>2</Subscript>O above the oceanic surface from Shanghai to Antarctica

During the 22nd Chinese Antarctic Research Expedition (CHINARE-22), the atmospheric gas samples above the oceanic surface and near the surface were collected on the track for the scientific ship “Xuelong” and on Millor Peninsula of eastern Antarctica, respectively, using the Tedlar gas bags. Every day the sampling times were 10:00 and 22:00 (local time), respectively. In the laboratory, high-precision measurement of the isotopic compositions for N₂O in these gas samples was conducted using Thermo Finnigan MAT-253 Isotopic Mass Spectrometer with a fully automated interface for the pre-GC concentration (PreCon) of trace gases. The temporal and spatial variations of δ ¹⁵N and δ ¹⁸O in atmospheric N₂O were analyzed. The mean δ ¹⁵N and δ ¹⁸O-N₂O values above the oceanic surface were (7.21±0.50)‰ and (44.52±0.52)‰, respectively. From 30°N to Antarctica, the δ ¹⁵N (6.05‰–7.88‰) linearly increased with the rate of about 0.01‰ with the latitude while the δ ¹⁸O (43.05‰–48.78‰) showed a large fluctuation. The δ ¹⁵N negatively correlated with air temperature and N₂O concentration, and slightly positively correlated with δ ¹⁸O. The summertime variations of δ ¹⁵N and δ ¹⁸O-N₂O appeared the same trend on Millor Peninsula of eastern Antarctica. They significantly positively correlated with each other and negatively with N₂O concentration. The δ ¹⁵N and δ ¹⁸O-N₂O at different sites averaged (7.42±0.35)‰ and (44.69±0.49)‰, respectively, slightly higher than those above the oceanic surface, significantly higher than those of atmospheric N₂O in the low-latitude regions of Northern Hemisphere. The predominant factors affecting the spatial variations of δ ¹⁵N and δ ¹⁸O values were also discussed. The isotopic data given in this study can help to investigate the global and regional N₂O budgets. Supported by the National Natural Science Foundation of China (Grant Nos. 40676005 and 40406001)     

Temporal variations of <Emphasis Type="Italic">δ</Emphasis><Superscript>18</Superscript>O of atmospheric water vapor at Delingha

Oxygen stable isotope of atmospheric water vapor is widely used to study the modern process of climate. Atmospheric water vapor samples were collected at Dlingha, northeast of Tibetan Plateau during the period from July 2005 to February 2006. The variation of δ ¹⁸O and the relationships between δ ¹⁸O and both the temperature and specific humidity are analyzed in this paper. Results show that the seasonal variation of δ ¹⁸O of atmospheric water vapor at Delingha is remarkable with higher δ ¹⁸O in summer and lower δ ¹⁸O in winter. The temporal variation of vapor δ ¹⁸O shows obvious fluctuations, with magnitude of over 37‰ The daily variation of the δ ¹⁸O is highly correlated with air temperature. The relationship between δ ¹⁸O and atmospheric water vapor content is complex. Study shows that δ ¹⁸O of atmospheric water vapor is positively correlated with specific humidity in winter in seasonal scale and inversely correlated with specific humidity in summer rainy period. The δ ¹⁸O values of atmospheric water vapor are lower than those of precipitation at Delingha, and the average difference is 10.7‰ Variations of δ ¹⁸O of atmospheric water vapor is also found to be affected by precipitation events, The model results show that the precipitation effect could have caused the vapor δ ¹⁸O in the raining season to lower by 7% in average in July and August. Supported by the National Natural Science Foundation of China (Grant Nos. 40671043, 40571039 and 40771048) and National Basic Research Program of China (Grant No. 2005CB422002)     

Atypical soil behavior during the 2011 Tohoku earthquake (<Emphasis Type="Italic">М</Emphasis><Subscript><Emphasis Type="Italic">w</Emphasis></Subscript> = 9)

To understand physical mechanisms of generation of abnormally high peak ground acceleration (PGA; >1g) during the Tohoku earthquake, models of nonlinear soil behavior in the strong motion were constructed for 27 KiK-net stations located in the near-fault zones to the south of FKSH17. The method of data processing used was developed by Pavlenko and Irikura, Pure Appl Geophys 160:2365–2379, 2003 and previously applied for studying soil behavior at vertical array sites during the 1995 Kobe (М _w ?=?6.8) and 2000 Tottori (М _w ?=?6.7) earthquakes. During the Tohoku earthquake, we did not observe a widespread nonlinearity of soft soils and reduction at the beginning of strong motion and recovery at the end of strong motion of shear moduli in soil layers, as usually observed during strong earthquakes. Manifestations of soil nonlinearity and reduction of shear moduli during strong motion were observed at sites located close to the source, in coastal areas. At remote sites, where abnormally high PGAs were recorded, shear moduli in soil layers increased and reached their maxima at the moments of the highest intensity of the strong motion, indicating soil hardening. Then, shear moduli reduced with decreasing the intensity of the strong motion. At soft-soil sites, the reduction of shear moduli was accompanied by a step-like decrease of the predominant frequencies of motion. Evidently, the observed soil hardening at the moments of the highest intensity of the strong motion contributed to the occurrence of abnormally high PGA, recorded during the Tohoku earthquake.     

測量大气臭氧垂直分布的逆轉方法<C>——对逆轉方法<A><B>的改进

本文给出了1个新的测量大气臭氧垂直分布的逆转方法,该方法将臭氧层分为6层,各层的臭氧含量,前后两次用图解法由联立方程定量的解出。使用这个方法,我们对在北京观测的10条逆转曲线作了计算,并将这些计算结果作了讨论及初步的误差估计。     

Simulating CO<Subscript>2</Subscript> transport into the ocean from a CO<Subscript>2</Subscript> lake at the seafloor using a <Emphasis Type="Italic">z</Emphasis>- and a <Emphasis Type="Italic">σ</Emphasis>-coordinate model

The ocean takes up approximately 2 GT carbon per year due to the enhanced CO₂ concentrations in the atmosphere. Several options have been suggested in order to reduce the emissions of CO₂ into the atmosphere, and among these are CO₂ storage in the deep ocean. Topographic effects of dissolution and transport from a CO₂ lake located at 3,000-m depth have been studied using the z-coordinate model Massachusetts Institute of Technology general circulation model (MITgcm) and the σ-coordinate model Bergen ocean model (BOM). Both models have been coupled with the general ocean turbulence model (GOTM) in order to account for vertical subgrid processes. The chosen vertical turbulence mixing scheme includes the damping effect from stable stratification on the turbulence intensity. Three different topographic scenarios are presented: a flat bottom and the CO₂ lake placed within a trench with depths of 10 and 20 m. The flat case scenario gives good correlation with previous numerical studies of dissolution from a CO₂ lake. When topography is introduced, it is shown that the z-coordinate model and the σ-coordinate model give different circulation patterns in the trench. This leads to different dissolution rates, 0.1 μmol cm^− 2 s^− 1 for the scenario of a 20-m-deep trench using BOM and 0.005–0.02 μmol cm^− 2 s^− 1 for the same scenario using the MITgcm. The study is also relevant for leakages of CO₂ stored in geological formations and to the ocean.     

Geochemical,isotopic and single crystal <Superscript>40</Superscript>Ar/<Superscript>39</Superscript>Ar age constraints on the evolution of the Cerro Galán ignimbrites

The giant ignimbrites that erupted from the Cerro Galán caldera complex in the southern Puna of the high Andean plateau are considered to be linked to crustal and mantle melting as a consequence of delamination of gravitationally unstable thickened crust and mantle lithosphere over a steepening subduction zone. Major and trace element analyses of Cerro Galán ignimbrites (68–71% SiO₂) that include 75 new analyses can be interpreted as reflecting evolution at three crustal levels. AFC modeling and new fractionation corrected δ¹⁸O values from quartz (+7.63–8.85‰) are consistent with the ignimbrite magmas being near 50:50 mixtures of enriched mantle (⁸⁷Sr/⁸⁶Sr ~ 0.7055) and crustal melts (⁸⁷Sr/⁸⁶Sr near 0.715–0.735). Processes at lower crustal levels are predicated on steep heavy REE patterns (Sm/Yb = 4–7), high Sr contents (>250 ppm) and very low Nb/Ta (9-5) ratios, which are attributed to amphibolite partial melts mixing with fractionating mantle basalts to produce hybrid melts that rise leaving a gravitationally unstable garnet-bearing residue. Processes at mid crustal levels create large negative Eu anomalies (Eu/Eu* = 0.45–0.70) and variable trace element enrichment in a crystallizing mush zone with a temperature near 800–850°C. The mush zone is repeatedly recharged from depth and partially evacuated into upper crustal magma chambers at times of regional contraction. Crystallinity differences in the ignimbrites are attributed to biotite, zoned plagioclase and other antecrysts entering higher level chambers where variable amounts of near-eutectic crystallization occurs at temperatures as low as 680°C just preceding eruption. ⁴⁰Ar/³⁹Ar single crystal sanidine weighted mean plateau and isochron ages combined with trace element patterns show that the Galán ignimbrite erupted in more than one batch including a ~ 2.13 Ma intracaldera flow and outflows to the west and north at near 2.09 and 2.06 Ma. Episodic delamination of gravitationally unstable lower crust and mantle lithosphere and injection of basaltic magmas, whose changing chemistry reflects their evolution over a steepening subduction zone, could trigger the eruptions of the Cerro Galán ignimbrites.     

Methods and technique of measurements of ρ<Subscript><Emphasis Type="Italic">k</Emphasis></Subscript> variations by GEOMES-R6 device

Upgrading of the electrical resistance variometer GEOMES-R6 of the Institute of Geophysics of the Poland Academy of Sciences, with enhancement of functional capabilities of the GEOMES assemblage for its application in laboratory measurements on rock samples, has been carried out. The assemblage operation principle, specifications, and theory of operation are given. Results of the assemblage test at measurements on marble sample and concrete model are presented. It is shown that the installation allows long-duration registration of variations of the value ρ _k on high-ohmic samples.     

Temperature evolution from the <Emphasis Type="Italic">δ</Emphasis><Superscript>18</Superscript>O record of Hani peat,Northeast China,in the last 14000 years

From the last deglaciation to the Holocene, the Greenland Ice Core (GISP2) δ ¹⁸O records as well as the records of ice-rafted debris on the surface of the North Atlantic have revealed a succession of sudden cooling events on the centennial to millennial scales. However, the temperature proxy records are rarely studied systematically and directly to ensure that this air temperature cooling pattern simultaneously existed in the East Asian Region, in addition to the repeated pattern occurring in the Greater Atlantic Region. A peat cellulose δ ¹⁸O temperature proxy record proximately existing for 14000 years was picked up from the Hani peat in Jilin Province, China. It suggests by comparison that the sudden cooling events, such as the Older Dryas, Inter-Allerød, Younger Dryas, and nine ice-rafted debris events of the North Atlantic, are almost entirely reiterated in the temperature signals of Hani peat cellulose δ ¹⁸O. These cooling events show that the repeatedly occurring temperature cooling pattern not only appeared in the North Atlantic Region in the high latitudes, but also in the Northwest Pacific Region in the middle latitudes. The climate change events marking the start of the Holocene Epoch, the Holocene Megathermal, the “8.2 kyr” event, the “4.2 kyr” event, the Medieval Warm Period, and the Little Ice Age are further discussed. The sensitivity response of Hani peat cellulose δ ¹⁸O to the land surface temperature and the reason for the age accuracy of peat cellulose ¹⁴C are also discussed based on the characteristics of the peat bog environment.     

Variations of δ<Superscript>18</Superscript>O in Rivers of Crimea in Winter

The results of studies of variations of δ¹⁸O values in river water in Crimea Peninsula in January–February 2015–2017 are given. The variation range of δ¹⁸О in river water over the three years of studies never exceeded 3‰. A tendency toward an increase in δ¹⁸О in the water of the Salgir, Kacha, Al’ma, Bel’bek, Biyuk-Karasu from their sources to mouths was identified and explained by a decrease in evaporation in the mouth areas of the rivers relative to their sources and upper reaches, and the inflow of isotopically light precipitation (rain and snow) into the rivers in their upper reaches. The values of δ¹⁸О in waters of the rivers with regulated stream were found to increase under the effect of the Simferopol Reservoir on the Salgir River and the Izobil’nenskoe Reservoir on the Ulu-Zen’ River. The values of δ¹⁸О in the upper reaches of the large rivers of Kacha and Bel’bek (the northwestern slope of the Crimean Mountains) vary from ?8.7 to ?9.7‰, except for the rivers of Al’ma (?7.7‰) and Kokozka (?10.2‰) because of the different shares of groundwater in the recharge of these rivers.     

Continuous soil CO<Superscript>2</Superscript> and discrete plume SO<Subscript>2</Subscript> measurements at Mt. Etna (Italy) during 1997–2000: a contribution to volcano monitoring

Continuous monitoring of soil CO₂ dynamic concentration (which is proportional to the CO₂ flux through the soil) was carried out at a peripheral site of Mt. Etna during the period November 1997–September 2000 using an automated station. The acquired data were compared with SO₂ flux from the summit craters measured two to three times a week during the same period. The high frequency of data acquisition with both methods allowed us to analyze in detail the time variations of both parameters. Anomalous high values of soil CO₂ dynamic concentration always preceded periods of increased flux of plume SO₂, and these in turn were followed by periods of summit eruptions. The variations were modeled in terms of gas efflux increase due to magma ascent to shallow depth and its consequent depressurization and degassing. This model is supported by data from other geophysical and volcanological parameters. The rates of increase both of soil CO₂ dynamic concentration and of plume SO₂ flux are interpreted to be positively correlated both to the velocity of magma ascent within the volcano and to lava effusion rate once magma is erupted at the surface. Low rates of the increase were recorded before the nine-month-long 1999 subterminal eruption. Higher rates of increase were observed before the violent summit eruption of September-November 1999, and the highest rates were observed during shorter and very frequent spike-like anomalies that preceded the sequence of short-lived but very violent summit eruptions that started in late January 2000 and continued until late June of the same year. Furthermore, the time interval between the peaks of CO₂ and SO₂ in a single sequence of gas anomalies is likely to be controlled by magma ascent velocity.Editorial responsibility: H. Shinohara     

The seasonal behavior of the <Emphasis Type="Italic">B</Emphasis><Subscript><Emphasis Type="Italic">z</Emphasis></Subscript> daily variations according to the monitoring at the Bishkek geodynamic site (1993–1995)

The algorithm for monitoring the magnetotelluric field has been briefly described, and the experimental data, obtained at the Bishkek geodynamic site, have been analyzed. The regularities in the seasonal behavior of the variations in the field vertical magnetic component have been revealed.     

The <Emphasis Type="Italic">M</Emphasis><Subscript><Emphasis Type="Italic">w</Emphasis></Subscript> = 5.8 14 August 2016 middle Sakhalin earthquake on a boundary between Okhotsk and Eurasian (Amurian) plates

An earthquake with the moment magnitude M _w ?=?5.8 occurred in the middle part of the Sakhalin Island, Russian Federation, on 14 August 2016, at 11:17 a.m. UTC. The earthquake source was located west of the Central Sakhalin Fault Zone, which is considered to mark the boundary between the Okhotsk and Eurasian (Amurian) plates. Moment tensor solution of the mainshock as well as the configuration of aftershock cloud suggests that the earthquake was caused by slip on a SW-dipping reverse fault. For the first time for Sakhalin, we have got the felt reports unified in accordance with DYFI. We also analyzed observed PGA values and, based on them, produced shaking maps.     

Evaluation of the dynamic <Emphasis Type="Italic">p</Emphasis>–<Emphasis Type="Italic">y</Emphasis><Subscript>p</Subscript> loops of pile-supported structures on sloping ground

A key issue in the design of pile-supported structures on sloping ground is soil–pile interaction, which becomes more complicated in case of dynamic loading. This study aimed to evaluate the effect of slope on the dynamic behavior of pile-supported structures by performing a series of centrifuge tests. Three models were prepared by varying the slope and soil density of dry sand grounds. The mass supported on 3 by 3 group piles was shaken applying sinusoidal wave with various amplitudes. Test results showed that the location of maximum values and distribution shape of the bending moment below the ground surface varied noticeably with the pile position in the slope case. The relationship between the soil resistance and pile deflection (p–y_p loops) was carefully evaluated by applying the piecewise cubic spline method to fit the measured bending moment curves along piles. It was found that the shape of the p–y_p loops was irregular due to the effect of slope, and immensely influenced by the movement of the unstable zone. In addition, the effect of the pile group in the horizontal case was evaluated by comparing with the previously suggested curves that represent the relationship between the soil resistance and pile–soil relative displacement (p–y curves) to propose the multiplier coefficients.     

REE and <Emphasis Type="Italic">ε</Emphasis><Subscript>Nd</Subscript> of clay fractions in sediments from the eastern Pacific Ocean: Evidence for clay sources

Clay fractions in the non-calcareous surface sediments from the eastern Pacific were analyzed for clay minerals, REE and ¹⁴³Nd/¹⁴⁴Nd. Montmorillonite/illite ratio (M/I ratio), total REE contents (ΣREE), LREE/HREE ratio and cerium anomaly (δCe) may effectively indicate the genesis of clay minerals. Clay fractions with M/I ratio <1, δCe >0.85, ΣREE <400 μg/g, LREE/HREE ratio ≈4, and REE patterns similar to those of pelagic sediments are terrigenous and autogenetic mixed clay fractions and contain more autogenetic montmorillonite. Clay fractions with M/I ratio >1, δCe=0.86 to 1.5, ΣREE=200 to 350 μg/g, LREE/HREE ratio ≈6 and REE distribution patterns similar to that of China loess are identified as terrigenous clay fraction. The ¹⁴³Nd/¹⁴⁴Nd ratios or _ɛNd values of clay fractions inherit the features of terrigenous sources of clay minerals. Clay fractions are divided into 4 types according to _ɛNd values. Terrigenous clay minerals of type I with the _ɛNd values of ™8 to ™6 originate mainly from North American fluvial deposits. Those of type II with the ε Nd values of ™9 to ™7 are mainly from the East Asia and North American fluvial deposits. Those of type III with ε Nd values of ™6 to ™3 could come from the central and eastern Pacific volcanic islands. Those of type IV with ε Nd values of ™13 to ™12 may be from East Asia eolian. The terrigenous and autogenetic mixed clay fractions show patchy distributions, indicating that there are volcanic or hot-spot activities in the eastern Pacific plate, while the terrigenous clay fractions cover a large part of the study area, proving that the terrigenous clay minerals are dominant in the eastern Pacific.     

Plasma pressure distribution in the equatorial plane of the Earth’s magnetosphere at geocentric distances of 6–10<Emphasis Type="Italic">R</Emphasis><Subscript><Emphasis Type="Italic">E</Emphasis></Subscript> according to the international THEMIS mission data

The structure of the averaged plasma pressure distribution in the plasma ring around the Earth at geocentric distances of ∼6–10R _E has been determined. The distribution function moments measured on the international THEMIS mission satellites have been used. The plasma pressure distribution in the equatorial plane at 15R _E > X_SM > −15R _E and 15R _E > Y_SM > −15R _E has been statistically studied. The radial dependence of the plasma pressure at the day-night and morning-evening meridians has been analyzed. It has been indicated that the plasma ring around the Earth has a structure, which is close to being azimuthally symmetric. The achieved results have been compared with the pressure distributions obtained previously. It has been indicated that in the overlapping regions, the achieved results agree with the previously obtained data within the pressure determination errors.     

Flow over a rounded backward-facing step,using a <Emphasis Type="Italic">z</Emphasis>-coordinate model and a <Emphasis Type="Italic">σ</Emphasis>-coordinate model

Homogeneous, nonrotating flow over a backward-facing rounded step is simulated using the 2D vertical version of two general circulation models, a z-coordinate model—the Massachusetts Institute of Technology general circulation model (MITgcm)—and a σ-coordinate model—the Bergen Ocean Model (BOM). The backward-facing step is a well-known testcase since it is geometrically simple but still embodies important flow characteristics such as separation point, reattachment length, and recirculation of the flow. The study compares the core of the two models and uses constant eddy viscosities and diffusivities. The Reynolds numbers ranges from 2·10² to 2·10⁶. The results correspond with previously published results having a relatively stationary separation point and a fluctuating reattachment length due to downslope propagating eddies released from the reattachment zone for Reynolds numbers higher than or equal to 2 · 10⁴. For Reynolds number within the laminar regime, the flow is stationary. The discrepancies between the models increase by enhancing Reynolds numbers. The σ-coordinate model experiences a reduction in eddy sizes with increasing resolution and Reynolds numbers in correspondence with published experiments, while the size of the eddies are independent of the Reynolds number using the MITgcm. Due to mixing generated by the staircase topography, the z-coordinate model gives a better convergence of the separation point and reattachment length compared with the BOM; however, this conclusion might change with the inclusion of a relevant turbulence scheme.