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.     

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)     

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.     

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.     

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.     

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)     

Was a dacite dome implicated in the 9,500 <Emphasis Type="SmallCaps">b.p.</Emphasis> collapse of Mt Ruapehu? A palaeomagnetic investigation

Constraining the process by which volcanoes become unstable is difficult. Several models have been proposed to explain the driving forces which cause volcanic edifices to catastrophically collapse. These include models for destabilisation of volcanic flanks by wedging due to dyke intrusion and the weakening of mechanical properties by pressurisation of pore fluids. It is not known which, if any, of the models are relevant to particular sector collapse events. Recent developments in the palaeomagnetic estimation of emplacement temperatures of volcaniclastic rocks have shown that even relatively low emplacement temperatures can be recorded by volcaniclastics with high fidelity. We have carried out a palaeomagnetic study of emplacement temperatures to investigate the role of igneous activity in the initiation of the 9,500 b.p. Murimotu sector collapse of Mt Ruapehu, New Zealand. This debris avalanche deposit has three fades which are stratigraphically superimposed, and the lowermost fades contains three lithological assemblages representing different segments of the edifice which were transported with little internal mixing within the flow. We have determined that some of the dacite-bearing assemblage 1, fades 1 was hot (∼350 °C) during transport and emplacement, whereas none of the other lithological assemblages of fades contained hot material. Our interpretation is that a dacite dome was active on the ancient Ruapehu edifice immediately prior to the Murimotu sector collapse. The partially cooled carapace of the dome and material shed from this part was incorporated into the avalanche deposit, along with cold lavas and volcaniclastics. We have not found evidence for incorporation of material at or close to magmatic temperatures, at least in the sampled locations. Our palaeomagnetic work allows us to develop a comprehensive, new palaeomagnetic classification of volcaniclastics. Published online: 25 January 2003 Editorial responsibility: D. Dingwell     

The Tolud Burst of Seismicity and the Earthquake of November 30, 2012 (<Emphasis Type="Italic">M</Emphasis><Subscript>C</Subscript> = 5.4, <Emphasis Type="Italic">M</Emphasis><Subscript>W</Subscript> = 4.8) that Accompanied the Start of the 2012‒2013 Tolbachik Eruption

This paper reports a study of the Tolud earthquake sequence; the sequence was a burst of shallow seismicity between November 28 and December 7, 2012; it accompanied the initial phase in the Tolbachik Fissure Eruption of 2012?2013. The largest earthquake (the Tolud earthquake of November 30, 2012, to be referred to as the Tolud Earthquake in what follows, with K_S = 11.3, M_L = 4.9, M_C = 5.4, and M_W = 4.8) is one of the five larger seismic events that have been recorded at depths shallower than 10 km beneath the entire Klyuchevskoi Volcanic Cluster in 1961?2015. It was found that the Tolud earthquake sequence was the foreshock–aftershock process of the Tolud Earthquake. This is one of the larger seismicity episodes ever to have occurred in the volcanic areas of Kamchatka. Data of the Kamchatka seismic stations were used to compute some parameters for the Tolud Earthquake and its largest (M_L = 4.3) aftershock; the parameters include the source parameters and mechanisms, and the moment magnitudes, since no information on these is available at the world seismological data centers. The focal mechanisms for the Tolud Earthquake and for its aftershock are consistent with seismic ruptures at a tension fault in the rift zone. Instrumental data were used to estimate the intensity of shaking due to the Tolud Earthquake. We discuss the sequence of events that was a signature of the time-dependent seismic and volcanic activity that took place in the Tolbachik zone in late November 2012 and terminated in the Tolud burst of seismicity. Based on the current ideas of the tectonics and magma sources for the Tolbachik volcanic zone, we discuss possible causes of these earthquakes.     

Structural analysis and interpretation of the surface deformation associated with the Ning’er,Yunnan Province,China <Emphasis Type="Italic">M</Emphasis><Subscript>S</Subscript>6.4 earthquake of June 3, 2007

The seismogenic fault and the dynamic mechanism of the Ning’er, Yunnan Province M_S6.4 earthquake of June 3, 2007 are studied on the basis of the observation data of the surface fissures, sand blow and water eruption, landslide and collapse associated with the earthquake, incorporating with the data of geologic structures, focal mechanism solutions and aftershock distribution for the earthquake area. The observation of the surface fissures reveals that the Banhai segment of the NW-trending Ning’er fault is dominated by right-lateral strike-slip, while the NNE-trending fault is dominated by left-lateral strike-slip. The seismo-geologic hazards are concentrated mainly within a 330°-extending zone of 13.5 km in length and 4 km in width. The major axis of the isoseismal is also oriented in 330° direction, and the major axis of the seismic intensity VIII area is 13.5 km long. The focal mechanism solutions indicate that the NW-trending nodal plane of the Ning’er M_S6.4 earthquake is dominated by right-lateral slip, while the NE-trending nodal plane is dominated by left-lateral slip. The preferred distribution orientation of the aftershocks of M_S≥2 is 330°, and the focal depths are within the range of 3～12 km, predominantly within 3～10 km. The distribution of the aftershocks is consistent with the distribution zone of the seismo-geologic hazards. All the above-mentioned data indicate that the Banhai segment of the Ning’er fault is the seismogenic fault of this earthquake. Moreover, the driving force of the Ning’er earthquake is discussed in the light of the active block theory. It is believed that the northward pushing of the Indian plate has caused the eastward slipping of the Qinghai-Tibetan Plateau, which has been transformed into the southeastern-southernward squeezing of the southwest Yunnan region. As a result, the NW-trending faults in the vicinity of the Ning’er area are dominated by right-lateral strike-slip, while the NE-trending faults are dominated by left-lateral strike-slip. This tectonic framework might be the main cause of the frequent occurrence of M_S6.0～6.9 earthquakes in the area.     

Automated,high time-resolution measurements of SO<Subscript>2</Subscript> flux at Soufrière Hills Volcano,Montserrat

We report here the first results from an automated, telemetered UV scanning spectrometer system for monitoring SO ₂ emission rates at Soufrière Hills Volcano, Montserrat. Two spectrometers receive light by way of a motor-driven stepping prism and telescope in order to make vertical scans of the volcanic plume. Spectral data from these spectrometers, situated 2,800 m apart and 4,500 m from the volcano, are relayed back to the observatory every 4–5 s via radio modems. A full scan of the plume is accomplished every 1–6 min by the (time-synchronised) spectrometers and a SO ₂ emission rate is calculated using the SO ₂ slant concentrations, scan angles and plume speeds estimated from the wind speed from a telemetered weather station near to the volcano. The plume's position and dimensions are calculated using the angular data from the two spectrometers. The plume height varies significantly diurnally and seasonally and is important in order to minimise the error on SO ₂ emission rates. The new scanning system (Scanspec) provides SO ₂ emission rates from 08:00 to 16:00 h local time every day. Preliminary results highlight a number of features of the SO ₂ time series and plume dynamics and give our first indications of the errors and limits of detection of this system. SO ₂ emission rates vary widely on all time scales (minutes, days, months). This new system has already provided the first real and consistent indication that SO ₂ emission rates vary on a minutes to hours basis, which can be correlated with volcanic activity (for example, rockfall and pyroclastic flow activity). It is anticipated that this system at Soufrière Hills will yield information on shallow processes occurring on short time scales (periods of minutes to hours) as well as deep processes relating to magma supply rates, which will be associated with longer wavelength SO ₂ signals of weeks to months.     

Ionospheric <Emphasis Type="Italic">N</Emphasis><Subscript>e</Subscript> disturbances before 2007 Pu’er,Yunnan, China,earthquake

Selecting three half orbits near the epicenter of Pu’er earthquake, we analyzed the N_e data recorded in their revisited orbits during a year before this earthquake, and extracted N_e precursors. The results show that: ① There are significant seasonal variations of ionospheric N_e in night time, which exhibit different shapes respectively in four seasons; ② There are three main shapes of N_e: single-peak, saddle-shaped and even-shaped, all of which may occur in four seasons, but each season with its typical shape relatively; ③ Spatial images of N_e showed high values near the epicenter in 30 days before the earthquake, and there is a good correlation between anomaly and distribution of earthquake in space and time, which reflects that these spatial anomalies were indeed concerned with the earthquake; ④ There shows a certain similarity of the N_e curves among revisited orbits, which can provide background information for distinguishing and identification of seismic anomaly.     

Reconstruction of the climate in the Tianmu Mountain area,Zhejiang Province,in the last 160 years by σ<Superscript>13</Superscript>C sequence of tree ring α-cellulose

Based on cross-dating tree rings from the Tianmu Mountain, Zhejiang Province, the tree rings α-cellulose σ¹³C time series was measured. By spectrum analysis, the variation of tree-ring σ¹³C sequence shows a quasi periodicity of 4.4 years, which is coincident with the standard cycle of El Niño. After removing the long-term decrease trend of the σ¹³C variation related to atmospheric CO₂ concentration, the response of the σ¹³C to climate elements was analyzed using meteorology data from a nearby weather station. The results suggest that there is a distinct relativity between high-frequency variation of tree ring σ¹³C series and seasonal climate parameters, e.g. temperature and precipitation, with a significant time-lag effect. In addition, the high frequency also reflects the strength change of the East Asian Monsoon. The multiple regression method was employed to reconstruct the historical climate, and to analyze the climate change and trend in the last 160 years in the northern Zhejiang Province.     

Strong motion characteristics of the <Emphasis Type="Italic">M</Emphasis><Subscript>w</Subscript> 6.6 Lushan earthquake,Sichuan, China — an insight into the spatial difference of a typical thrust fault earthquake

Near-field strong ground motions are useful for engineering seismology studies and seismic design, but dense observation networks of damaging earthquakes are still rare. In this study, based on the strong-motion data from the M w 6.6 Lushan earthquake, the ground motion parameters in different spatial regions are systematically analyzed, and the contributions from different effects, like the hanging-wall effect, directivity effect, and attenuation effect are separated to the extent possible. Different engineering parameters from the observed ground motions are compared with the local design response spectra and a new attenuation relation of Western China. General results indicate that the high frequency ground motion, like the peak ground acceleration, on two sides of the fault plane is sensitive to the hanging-wall effect, whereas the low frequency ground motion, like the long period spectral acceleration, in the rupture propagation direction is affected by the directivity effect. Moreover, although the M w 6.6 Lushan earthquake is not a large magnitude event, the spatial difference of ground motion is still obvious; thus, for a thrust faulting earthquake, in addition to the hanging effect, the directivity effect should also be considered.     

Altitudinal effect of soil <Emphasis Type="Italic">n</Emphasis>-alkane <Emphasis Type="Italic">δ</Emphasis>D values on the eastern Tibetan Plateau and their increasing isotopic fractionation with altitude

Stable isotope paleoaltimetry has provided unprecedented insights into the topographic histories of many of the world’s highest mountain ranges. However, on the Tibetan Plateau (TP), stable isotopes from paleosols generally yield much higher paleoaltitudes than those based on fossils. It is therefore essential when attempting to interpret accurately this region’s paleoaltitudes that the empirical calibrations of local stable isotopes and the relations between them are established. Additionally, it is vital that careful estimations be made when estimate how different isotopes sourced from different areas may have been influenced by different controls. We present here 29 hydrogen isotopic values for leaf wax-derived n-alkanes (i.e., _δD_wax values, and abundance-weighted average δD values of C₂₉ and C₃₁) in surface soils, as well as the δD values of soil water (δD_sw) samples (totaling 22) from Mount Longmen (LM), on the eastern TP (altitude ~0.8–4.0 km above sea level (asl), a region climatically affected by the East Asian Monsoon (EAM). We compared our results with published data from Mount Gongga (GG). In addition, 47 river water samples, 55 spring water samples, and the daily and monthly summer precipitation records (from May to October, 2015) from two precipitation observation stations were collected along the GG transect for δD analysis. LM soil _δD_wax values showed regional differences and responded strongly to altitude, varying from?160‰ to?219‰, with an altitudinal lapse rate (ALR) of?18‰ km^?1 (R ²=0.83; p<0.0001; n=29). These _δD_wax values appeared more enriched than those from the GG transect by ~40‰. We found that both the climate and moisture sources led to the differences observed in soil _δD_wax values between the LM and GG transects. We found that, as a general rule, ε _wax/rw, ε _wax/p and ε _wax/sw values (i.e., the isotopic fractionation of _δD_wax corresponding to δD_rw, δD_p and δD_sw) increased with increasing altitude along both the LM and GG transects (up to 34‰and 50‰, respectively). Basing its research on a comparative study of _δD_wax, δD_p, δD_rw(δD_springw) and δD_sw, this paper discusses the effects of moisture recycling, glacier-fed meltwater, relative humidity (RH), evapotranspiration (ET), vegetation cover, latitude, topography and/or other factors on ε _wax/p values. Clearly, if ε _wax-p values at higher altitudes are calculated using smaller ε _wax-p values from lower altitudes, the calculated paleowaterδD_p values are going to be more depleted than the actual δD values, and any paleoaltitude would therefore be overestimated.