Preliminary study on the origin identification of natural gas by the parameters of light hydrocarbon

The light hydrocarbon composition of 209 natural gas samples and individual light hydrocarbon carbon isotopes of 53 natural gas samples from typical humic-sourced gas and sapropelic-sourced gas in the four basins of China have been determined and analyzed. Some identification parameters for humic-sourced gas and sapropelic-sourced gas are proposed or corrected. The differences of compound-specific δ ¹³C value of individual light hydrocarbon between humic-sourced gas and sapropelic-sourced gas have been founded. The humic-sourced gas has the distribution of δ ¹³C_benzene> ?24‰, δ ¹³C_toluene >?23‰, δ ¹³C_cyclohexane > ?24‰ and δ ¹³C_{methyl cyclohexane}> ?24‰, while the sapropelic-sourced gas has the distribution of δ ¹³C_benzene <?24‰, δ ¹³C_toluene< ?24‰, δ ¹³C_cyclohexane< ?24‰ and δ ¹³C_{methyl cyclohexane}< ?24‰. Among the components of C₇ light hydrocarbon compound, such as normal heptane (nC₇), methyl cyclohexane (MCH) and dimethyl cyclopentane (ΣDMCP), etc, relative contents of nC₇ and MCH are influenced mainly by the source organic matter type of natural gas. Therefore, it is suggested that the gas with relative content of nC₇ of more than 30% and relative content of MCH of less than 70% is sapropelic-sourced gas, while gas with relative content of nC₇ of less than 35% and relative content of MCH of more than 50% is humic-sourced gas. Among components of C_5–7 aliphatics, the gas with relative content of C_5–7 normal alkane of more than 30% is sapropelic-sourced gas, while the gas with relative content of C_5–7 normal alkane of less than 30% is humic-sourced gas. These paremeters have been suggested to identify humic-sourced gas and sapropelic-sourced gas.     

Genetic features of gas-condensate reservoirs in the Kekeya Field,southwest depression in the Tarim Basin,China

Oils, condensates and natural gases in the Kekeya Field, southeast depression of the Tarim Basin were studied for their geochemical characteristics. According to the distribution analysis of the C₂/C₃ values with C₁/C₂ values, C₂/C₃ values with C₁/C₃ values, as well as C₂/C₃ values with dryness index, there are two different types of natural gases in the studied field, which are spatially regularly distributed. One is the oil cracking gas, located on shallow reservoirs over X ₅ ² reservoir, namely Upper oil legs; the other is kerogen cracking gas, located on X ₇ ² reservoirs, X₈ reservoirs and E₂ k reservoirs, namely Lower oil legs. In addition, the distribution patterns of molar concentration of oils and condensates with different carbon numbers of the n-alkanes in the Kekeya Field indicate that the crude oils have experienced several kinds of secondary alterations, which were closely related to the charging of gaseous hydrocarbons after petroleum accumulation. These results indicate that, based on the research of δ ¹³C values of individual hydrocarbons, heptane values and isoheptane values of light hydrocarbons and aromatic maturity parameters for oils, condensates and natural gases, oils and gases were charged at different geological time in the Kekeya Field.     

Influence of aquatic plants on the hydrogen isotope composition of sedimentary long-chain n-alkanes in the Lake Qinghai region,Qinghai-Tibet Plateau

The hydrogen isotopic composition(δD) of leaf wax long-chain n-alkanes(C27, C29, and C31) from lacustrine sediments has been widely applied to reconstruct terrestrial paleoclimatic and paleohydrological changes. However, few studies have addressed whether the aquatic-derived n-alkanes can affect the δD values of lake sedimentary long-chain n-alkanes, which are usually regarded as a recorder of the terrestrial hydrological signals. Here we systematically investigated δD values of long-chain n-alkanes from modern aquatic plants, both near-shore and off-shore surface sediments, surrounding terrestrial plant litters, as well as river water and lake water in Lake Qinghai and its satellite lakes on the northeastern Qinghai-Tibet Plateau. Our data showed that(i) δD values of long-chain n-alkanes from aquatic plants varied from-184‰ to-132‰ for n-C27, from-183‰ to-138‰ for n-C29, and from-189‰ to-130‰ for n-C31, respectively, with no significant differences among the three n-alkanes homologues;(ii) δD values of long-chain n-alkanes from aquatic plants were generally more positive than those from surrounding terrestrial plants, possibly because that they recorded the D-enrichment of lake water in this semi-arid region;(iii) δD values of long-chain n-alkanes from surface sediments showed significant differences among the three n-alkanes homologues, due to the larger aquatic input of n-C27 to the sedimentary lipid pool than that of n-C31, and(iv) n-C27 δD values of near-shore aquatic plants and near-shore sediments are more negative than those from off-shore as a result of lower δD values of near-shore lake water. Our findings indicate that in this region(i) the offset between sedimentary n-C27 and n-C31 δD values(ΔδDC27-C31) could potentially be used to evaluate if sedimentary long-chain n-alkanes are derived from a single source;(ii) while δD values of n-C27 may be influenced by lake water hydrological changes, sedimentary n-C31 is derived predominantly from terrestrial plants and thus its δD can serve as a relatively reliable indicator for terrestrial paleoclimatic and paleohydrological reconstructions.     

Characterization of the negative carbon isotope shift in segment C2, its global implications as a harbinger of OAE1a

Lower Cretaceous C-isotope records show intermittent negative/positive spikes, and consistent patterns of coeval chemostratigraphic curves thus document shifts that signal simultaneous responses of temporal changes in the global carbon reservoir. The standard pattern registered by the δ ¹³C_org and δ ¹³C_carb in Lower Aptian sediments includes distinct isotope segments C1 to C8 (Menegatti et al., 1998). In the El Pui section, Organyà Basin, Spain, C-isotope segment C2 is the longest interval preceding segments C3–C6 associated with oceanic anoxic event 1a (OAE 1a), and reveals a distinct negative shift of ~1.8‰ to ~2.23‰ defining the C-isotope pattern within that interval. Total inorganic carbon (TIC), total organic carbon (TOC), δ ¹³C_org, microfacies, n-alkanes show no difference before, during, or after the negative inflection. The biomarkers indicate that organic matter (OM) mainly originates from algal/microbial sources because short-chain length homologues (≤nC₁₉) dominate. nC₂₀ through nC₂₅ indicate some contribution from aquatic vegetation, but little from higher plants (>nC₂₅), as also suggested by the terrestrial/aquatic ratio of n-alkanes or (TAR) = [(nC₂₇+nC₂₉+nC₃₁)/(nC₁₅+nC₁₇+nC₁₉)] (averages 0.085). We suggest that conjoint pulses of contemporaneous LIPs (Ontong Java) and massive explosive volcanism in northeast Asia, the Songliao Basin (SB-V), best conform to plausible causes of the negative intra-C2 carbon isotopic excursion (CIE) at that time. Because of its apparent common occurrence the intra-C2 inflection could be a useful marker harbinger to the more pronounced CIE C3, the hallmark of OAE1a.     

The effect of volatile components in oil on evolutionary characteristics of diamondoids during oil thermal pyrolysis

On the basis of the results of simulation experiments, now we better understand the contribution of high carbon number hydrocarbons to diamondoid generation during thermal pyrolysis of crude oil and its sub-fractions(saturated, aromatic, resin, and asphalene fractions). However, little is known about the effect of volatile components in oil on diamondoid generation and diamondoid indices due to the lack of attention to these components in experiments. In this study, the effect of volatile components in oil on diamondoid generation and maturity indices was investigated by the pyrolysis simulation experiments on a normal crude oil from the HD23 well of the Tarim Basin and its residual oil after artificial volatilization, combined with quantitative analysis of diamondoids. The results indicate that the volatile components(≤n C12) in oil have an obvious contribution to the generation of adamantanes, which occurs mainly in the early stage of oil cracking(Easy Ro1.0%), and influences the variations in maturity indices of adamantanes; but they have no obvious effect on the generation and maturity indices of diamantanes. Therefore, some secondary alterations e.g., migration, gas washing, and biodegradation, which may result in the loss of light hydrocarbons in oil under actual geological conditions, could affect the identification of adamantanes generated during the late-stage cracking of crude oil, and further influence the practical application of adamantane indices.     

Identification of marine natural gases with different origin sources

Kinetic experiments of gas generation for typical samples of marine gas precursors including low-maturity kerogen, residual kerogen and oil as well as dispersed liquid hydrocarbon (DLH) in source rocks were performed by closed system, and the evolution trends of molecular and isotopic compositions of natural gases from different precursors against the maturity (R ₀%) at laboratory conditions were analyzed. Several diagrams of gas origin were calibrated by using the experimental data. A diagram based on the ratio of normal and isomerous butane and pentane (i/nC₄ ? i/nC₅) was proposed and used to identify the origins of the typical marine natural gases in the Sichuan Basin and the Tarim Basin, China. And the maturities of natural gases were estimated by using the statistical relationships between the gaseous molecular carbon isotopic data and maturities (δ¹³C-R ₀%) with different origins. The results indicate that the molecular and isotopic compositions of simulated gases from different precursors are different from each other. For example, the dryness index of the oil-cracking gas is the lowest; the dryness indices of gases from DLH and kerogen in closed system are almost the same; and the dryness index of gases from residual kerogen is extremely high, indicating that the kerogen gases are very dry; the contents of non-hydrocarbon gases in kerogen-cracking gases are far higher than those in oil-cracking and DLH-cracking gases. The molecular carbon isotopes of oil-cracking gases are the lightest, those of kerogen in closed system and GLH-cracking gases are the second lightest, and those of cracking gases from residual kerogen are the heaviest. The calibration results indicate that the diagrams of In(C₁/C₂)-In(C₂/C₃) and δ^{4 3}C₂-δ^{4 3}C₃-In(C₂/C₃) can discriminate primary and secondary cracking gases, but cannot be used to identify gas origin sources, while the diagram of i/nC₄ ? i/nC₅ can differentiate the gases from different precursors. The application results of these diagrams show that gas mixtures extensively exist in China, which involved the gases from multiple precursors and those from different maturity stages. For example, marine gases in the Sichuan Basin involve the mixture of oil-cracking gases and high-over-maturated kerogen gases, while those in the Tarim Basin involve not only the mixture of gases from multiple precursors, but also those from different maturity gases and post-reservoir alternations such as oxidized degradation and gas intrusion processes.     

Classification and origin of natural gases from Lishui Sag,the East China Sea Basin

Natural gases discovered up to now in Lishui Sag, the East China Sea Basin, differ greatly in gaseous compositions, of which hydrocarbon gases amount to 2%–94% while non-hydrocarbon gases are dominated by CO₂. Their hydrocarbon gases, without exception, contain less than 90% of methane and over 10% of C₂ ⁺ heavier hydrocarbons, indicating a wet gas. Carbon isotopic analyses on these hydrocarbon gases showed that δ ¹³C₁, δ ¹³C₂ and δ ¹³C₃ are basically lighter than ?44‰, ?29‰ and ?26‰, respectively. The difference in carbon isotopic values between methane and ethane is great, suggesting a biogenic oil-type gas produced by the mixed organic matter at peak generation. δ ¹³ \(C_{CO_2 } \) values of nonhydrocarbon gases are all heavier than ?10‰, indicating a typical abiogenic gas. The simulation experiment on hydrocarbon generation of organic matter in a closed gold-tube system showed that the proportion of methane in natural gases produced by terrigenous organic matter in the Lingfeng Formation marine deposit is obviously higher than that in natural gases derived from the aquatic and terrigenous mixed organic matter in the Yueguifeng Formation lacustrine deposit, consequently the proportion of heavier hydrocarbons of the former is remarkably lower than that of the latter. Moreover, δ ¹³C₁ values of natural gases produced by terrigenous organic matter in the Lingfeng Formation marine deposit are about 5‰ heavier than those of natural gases derived from the aquatic and terrigenous mixed organic matter in the Yueguifeng Formation lacustrine deposit while δ ¹³C₂ and δ ¹³C₃ values of the former are over 9‰ heavier than those of the latter. Currently the LS36-1 oil-gas pool is the only commercial oil-gas reservoir in Lishui Sag, where carbon isotopic compositions of various hydrocarbon components differ greatly from those of natural gases produced by the Lingfeng Formation organic matter but are very similar to those of natural gases derived from the Yueguifeng Formation organic matter, therefore, natural gases in the LS36-1 oil-gas pool are mainly derived from the Yueguifeng Formation lacustrine source rock rather than the Lingfeng Formation marine or Mingyuefeng Formation coal-measures source rocks.     

Catchment environmental change over the 20th Century recorded by sedimentary leaf wax n-alkane δ~(13)C off the Pearl River estuary

The compound-specific stable carbon isotope compositions(δ~(13)C) of leaf wax n-alkanes from two short sediment cores recovered off the Pearl River estuary(PRE) were analyzed to check for their capability of indicating decadal scale catchment environmental change. Sedimentary long-chain n-alkanes exhibited an odd-over-even predominance, with a maximum at n-C_(29) or n-C_(31), indicating their leaf wax origin was from vascular plants. The δ~(13)C values of C_(29) and C_(31) n-alkane in all the sediment samples were in the range of -28.8‰ to -31.2‰, consistent with the C_3 plant-dominated vegetation in the Pearl River catchments. The time series of δ~(13)C records from the two cores were comparable and displayed a decreasing trend from the early 20 th century to the end of the 1970s, followed by a reversal in that change leading to continued increase for ca. 15 years. After being corrected for the effect of atmospheric CO_2 rise and δ~(13)C_(atm) decline, the δ~(13)C_(29) records largely retained their raw changing pattern; the post-1980 increase being more conspicuous. The slightly decreasing trend in corrected δ~(13)C records before around 1980 may have been caused by an increase in precipitation, whereas the subsequent increase of δ~(13)C is likely associated with the observed dry climate and/or intensive anthropogenic deforestation. Our results thus demonstrate that leaf wax n-alkanes buried in the sediments off the PRE may well reflect change in the regional climate and/or human activity in the river catchments over the past century.     

Natural gas geochemistry and its origins in Kuqa depression

According to gas compositional and carbon isotopic measurement of 114 gas samples from the Kuqa depression, accumulation of the natural gases in the depression is dominated by hydrocarbon gases, with high gas dryness (C₁/C_1–4) at the middle and northern parts of the depression and low one towards east and west sides and southern part. The carbon isotopes of methane and its homologues are relatively enriched in ¹³C, and the distributive range of δ ¹³C₁, δ ¹³C₂ and δ ¹³C₃ is ?32‰–?36‰, ?22‰–?24‰ and ?20‰–?22‰, respectively. In general, the carbon isotopes of gaseous alkanes become less negative with the increase of carbon numbers. The δ ¹³ \(C_{CO_2 } \) value is less than ?10‰ in the Kuqa depression, indicating its organogenic origin. The distributive range of ³He/⁴He ratio is within n × 10^?8 and a decrease in ³He/⁴He ratio from north to south in the depression is observed. Based on the geochemical parameters of natural gas above, natural gas in the Kuqa depression is of characteristics of coal-type gas origin. The possible reasons for the partial reversal of stable carbon isotopes of gaseous alkanes involve the mixing of gases from one common source rock with different thermal maturity or from two separated source rock intervals of similar kerogen type, multistages accumulation of natural gas under high-temperature and over-pressure conditions, and sufficiency and diffusion of natural gas.     

Experimental study on the shear behavior of the interface between cushion materials and the concrete raft

Cushion is a layer of granular materials between the raft and the ground. The shear behavior of the interface between the cushion and the raft may influence the seismic performance of the superstructure. In order to quantify such influences, horizontal shear tests on the interfaces between different cushion materials and concrete raft under monotonic and cyclic loading were carried out. The vertical pressure P_v, material type and cushion thickness h_c were taken as variables. Conclusions include: 1) under monotonic loading, P_v is the most significant factor; the shear resistance P_(hmax) increases as P_v increases, but the normalized factor of resistance μ_n has an opposite tendency; 2) for the materials used in this study, μ_n varies from 0.40 to 0.70, the interface friction angle δ_s varies from 20° to 35°, while u_(max) varies from 3 mm to 15 mm; 3) under cyclic loading, the interface behavior can be abstracted as a "three-segment" back-bone curve, the main parameters include μ_n, the displacement u_1 and stiffness K_1 of the elastic stage, the displacement u_2 and stiffness K_2 of the plastic stage; 4) by observation and statistical analysis, the significance of different factors, together with values of K_1, K_2 and μ_n have been obtained.     

Effect of zonal E × B plasma drift on electron density in the low-latitude ionospheric <Emphasis Type="Italic">F</Emphasis> region at a solar activity maximum near vernal equinox

The variations in the density of the ionospheric F2 layer maximum (NmF2) under the action of the zonal plasma drift perpendicularly to the magnetic (B) and electric (E) fields in the direction geomagnetic west-geomagnetic east have been studied using the three-dimensional nonstationary theoretical model of electron and ion densities (N _e and N _i ) and temperatures (T _e and T _i ) in the low-latitude and midlatitude ionospheric F region and plasmasphere. The method of numerical calculations of N _e , N _i , T _e , and T _i , including the advantages of the Lagrangian and Eulerian methods, is used in the model. A dipole approximation of the geomagnetic field (B), taking into account the non-coincidence of the geographic and geomagnetic poles and differences between the positions of the Earth’s and geomagnetic dipole centers, is accepted in the calculations. The calculated NmF2 and altitudes of the F2 layer maximum (hmF2) have been compared with these quantities measured at 16 low-latitude ionospheric sounding stations during the geomagnetically quiet period October 11–12, 1958. This comparison made it possible to correct the input model parameters: the NRLMSISE-00 model [O], the meridional component of the neutral wind velocity according to the HWW90 model, and the meridional component of the equatorial plasma drift due to the electric field specified by the empirical model. It has been indicated that the effect of the zonal E × B plasma drift on NmF2 can be neglected under daytime conditions and changes in NmF2 and hmF2 under the action of this drift are insignificant under nighttime conditions north of 25° and south of ?26° geomagnetic latitude. The effect of the zonal E × B plasma drift on NmF2 and hmF2 is most substantial in the nightside ionosphere approximately from ?20° to 20° geomagnetic latitude, and the neglect of this drift results in an up to 2.4-fold underestimation of NmF2. The found dependence of the effect of the zonal E × B plasma drift on NmF2 and hmF2 on geomagnetic latitude is related to the longitudinal asymmetry of B, asymmetry of the neutral wind about the geomagnetic equator, and changes in the meridional E × B plasma drift at a change in geomagnetic longitude.     

Geochemical characteristics and origin of gases from the Upper,Lower Paleozoic and the Mesozoic reservoirs in the Ordos Basin,China

The Ordos Basin, the second largest sedimentary basin in China, contains the broad distribution of natural gas types. So far, several giant gas fields have been discovered in the Upper and Lower Paleozoic in this basin, each having over 1000×10⁸m³ of proven gas reserves, and several gas pools have also been discovered in the Mesozoic. This paper collected the data of natural gases and elucidated the geochemical characteristics of gases from different reservoirs, and then discussed their origin. For hydrocarbons preserved in the Upper Paleozoic, the elevated δ ¹³C values of methane, ethane and propane indicate that the gases would be mainly coal-formed gases; the singular reversal in the stable carbon isotopes of gaseous alkanes suggests the mixed gases from humic sources with different maturity. In the Lower Paleozoic, the δ ¹³C₁ values are mostly similar with those in the Upper Paleozoic, but the δ ¹³C₂ and δ ¹³C₃ values are slightly lighter, suggesting that the gases would be mixing of coal-type gases as a main member and oil-type gases. There are multiple reversals in carbon isotopes for gaseous alkanes, especially abnormal reversal for methane and ethane (i.e. δ ¹³C₁>δ ¹³C₂), inferring that gases would be mixed between high-mature coal-formed gases and oil-type gases. In the Mesozoic, the δ ¹³C values for gaseous alkanes are enriched in ¹²C, indicating that the gases are mainly derived from sapropelic sources; the carbon isotopic reversal for propane and butane in the Mesozoic is caused by microbial oxidation and mixing of gases from sapropelic sources with different maturity. In contrast to the Upper Paleozoic gases, the Mesozoic gases are characterized by heavier carbon isotopes of iso-butane than normal butane, which may be caused by gases generated from different kerogen types. Finally, according to δ ¹³C₁-R ₀ relationship and extremely low total organic carbon contents, the Low Paleozoic gases would not be generated from the Ordovician source as a main gas source, bycontrast, the Upper Paleozoic source as a main gas source is contributed to the Lower Paleozoic gases.     

Detailed Study of Time Variations in the Gutenberg–Richter <Emphasis Type="Italic">b</Emphasis>-value Based on Highly Accurate Seismic Observations at the Garm Prognostic Site,Tajikistan

The time variations in the Gutenberg–Richter b-value are minutely studied based on the data of highly accurate seismological observations at the Garm prognostic site, Tajikistan, where a stationary network of seismic stations of the Complex Seismological Expedition (CSE) of Schmidt Institute of Physics of the Earth (IPE) of the USSR (Russian) Academy of Sciences was in operation from 1955 to 1992. A total of 93035 local earthquakes ranging from 0.0 to 6.3 in the Ml magnitudes are considered. The spatiotemporal fluctuations in the minimal magnitude of completeness of the earthquakes, Mc, are analyzed. The study considers a 25-year interval of the observations at the center of the observation system within which Mc = 0.9. It is shown that in most cases, the b-value and log₁₀E^2/3 experience characteristic time variations before the earthquakes with magnitudes higher than the minimal magnitude of the predicted earthquake (MPE). The 6-year anomaly in the parameters’ b-value, log₁₀E^2/3, and log₁₀N associated with the single strongest earthquake with M = 6.3 that occurred in the observation region on October 26, 1984 is revealed. The inversely proportional relationship is established between the time variations in the b-value and the time variations in the velocities of seismic waves Vp and Vp/Vs. It is shown that the exponent p in the power function which links the time variations of the b-value and log₁₀E^2/3 is higher in the zones of crustal compression than in the zones of extension. It is simultaneously confirmed that the average b-value in the zones of compression is lower than in the zones of extension. It is established that in the case of earthquakes with M ≥ 2.6, the time series of seismic activity log₁₀N_i and the time series of the b-value are highly cross correlated with a coefficient of r ≈ 0.75, whereas in the case of earthquakes with M ≥ 0.9, the coefficient of cross correlation between these time series is close to zero (r ≈ 0.06). The law of variations in the slope of the lines approximating the relationship between the log₁₀N_i time series in the different magnitude ranges (M ≥ Mc_i) and b-value time series is obtained. It is hypothesized that the seismic activity of the earthquakes with high magnitudes can be estimated provided that the parameters of the time series of the b-value and time series of the number of earthquakes logN(ΔM_i) in the range of low magnitudes are known. It is concluded that using the parameter log₁₀N for prognostic estimates of the strong earthquakes only makes sense for earthquakes having moderate and large magnitudes. It is inferred that the time variations in the b-value are predominantly contributed by the time variations of the earthquakes with relatively large magnitudes.     

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.     

Nonadditive entropy <Emphasis Type="Italic">Sq</Emphasis> and nonextensive statistical mechanics: Applications in geophysics and elsewhere

The celebrated Boltzmann-Gibbs (BG) entropy, S _BG = ?kΣ _i p _i ln p _i , and associated statistical mechanics are essentially based on hypotheses such as ergodicity, i.e., when ensemble averages coincide with time averages. This dynamical simplification occurs in classical systems (and quantum counterparts) whose microscopic evolution is governed by a positive largest Lyapunov exponent (LLE). Under such circumstances, relevant microscopic variables behave, from the probabilistic viewpoint, as (nearly) independent. Many phenomena exist, however, in natural, artificial and social systems (geophysics, astrophysics, biophysics, economics, and others) that violate ergodicity. To cover a (possibly) wide class of such systems, a generalization (nonextensive statistical mechanics) of the BG theory was proposed in 1988. This theory is based on nonadditive entropies such as \(S_q = k\frac{{1 - \sum\nolimits_i {p_i^q } }}{{q - 1}}\left( {S_1 = S_{BG} } \right)\). Here we comment some central aspects of this theory, and briefly review typical predictions, verifications and applications in geophysics and elsewhere, as illustrated through theoretical, experimental, observational, and computational results.     

Earthquake hazard and risk assessment based on Unified Scaling Law for Earthquakes: Greater Caucasus and Crimea

We continue applying the general concept of seismic risk analysis in a number of seismic regions worldwide by constructing regional seismic hazard maps based on morphostructural analysis, pattern recognition, and the Unified Scaling Law for Earthquakes (USLE), which generalizes the Gutenberg-Richter relationship making use of naturally fractal distribution of earthquake sources of different size in a seismic region. The USLE stands for an empirical relationship log₁₀N(M, L)?=?A?+?B·(5 – M)?+?C·log₁₀L, where N(M, L) is the expected annual number of earthquakes of a certain magnitude M within a seismically prone area of linear dimension L. We use parameters A, B, and C of USLE to estimate, first, the expected maximum magnitude in a time interval at seismically prone nodes of the morphostructural scheme of the region under study, then map the corresponding expected ground shaking parameters (e.g., peak ground acceleration, PGA, or macro-seismic intensity). After a rigorous verification against the available seismic evidences in the past (usually, the observed instrumental PGA or the historically reported macro-seismic intensity), such a seismic hazard map is used to generate maps of specific earthquake risks for population, cities, and infrastructures (e.g., those based on census of population, buildings inventory). The methodology of seismic hazard and risk assessment is illustrated by application to the territory of Greater Caucasus and Crimea.     

Geochemical characteristics of pyrolysis gas from epimetamorphic rocks in the northern basement of Songliao Basin,Northeast China

The Xushen gas field, located in the north of Songliao Basin, is a potential giant gas area for China in the future. Its proved reserves have exceeded 1000×10⁸ m³ by the end of 2005. But, the origin of natural gases from the deep strata is still in debating. Epimetamorphic rocks as a potential gas source are widely spreading in the northern basement of Songliao Basin. According to pyrolysis experiments for these rocks in the semi-confined system, gas production and geochemistry of alkane gases are discussed in this paper. The Carboniferous-Permian epimetamorphic rocks were heated from 300°C to 550°C, with temperature interval of 50°C. The gas production was quantified and measured for chemical and carbon isotopic compositions. Results show that δ ¹³C₁ is less than ?20‰, carbon isotope trend of alkane gas is δ ¹³C₁<δ ¹³C₂<δ ¹³C₃ or δ ¹³C₁<δ ¹³C₂>δ ¹³C₃, these features suggest that the gas would be coal-type gas at high-over maturity, not be inorganic gas with reversal trend of gaseous alkanes (δ ¹³C₁>δ ¹³C₂>δ ¹³C₃). These characteristics of carbon isotopes are similar with the natural gas from the basin basement, but disagree with gas from the Xingcheng reservoir. Thus, the mixing gases from the pyrolysis gas with coal-typed gases at high-over maturity or oil-typed gases do not cause the reversal trend of carbon isotopes. The gas generation intensity for epimetamorphic rocks is 3.0×10⁸–23.8×10⁸ m³/km², corresponding to R _o from 2.0% to 3.5% for organic matter.     

Geochemical study on oil-cracked gases and kerogen-cracked gases (II)—Discrimination methods between oil-cracked gases and kerogen-cracked gases

In the processes of discrimination between oil-cracked gases and kerogen-cracked gases, Behar and Pinzgofer et al.’s results were adopted in the former researches, in which the ratio of C₂/C₃ is basically a constant while the ratio of C₁/C₂ gradually increases in the course of primary cracking of kerogen. Otherwise in the course of secondary cracking of oil, the ratio of C₂/C₃ increases rapidly while C₁/C₂ keeps relatively stable. Our study on analogue experiment shows that, whether it is oil or kerogen, in its process of gas generating by cracking, the ratios of C₂/C₃, C₁/C₂ or C₁/C₃ will all be increased with the growth of thermal conditions. In comparison, the ratio of C₂/C₃, which is affected by genetic type to some comparatively less extent, mainly responds to the maturity of gases, while the value of C₂/C₃ is about 2, and that of C₂/iC₄ is about 10, and the corresponding value of R _o is about 1.5%–1.6%. The influence of gas source on C₂/C₃ is less than that of gas maturity, otherwise C₁/C₂ (or C₁/C₃) is obviously affected by cracking matrices. The ratios of C₁/C₂, C₁/C₃ of oil-cracked gases are less than that of kerogen-cracked gases, under the condition that the ratios of C₂/C₃ are similar in value, so are the value of dryness indexes. There exists wide diffidence between this view and the former discrimination method in theory. The analysis of the spot sample indicates that we can apply the above basic view to dealing efficiently with the problem of the discrimination between oil-cracked gas and kerogen-cracked gas.     

Correlation between foam-bubble size and drag coefficient in hurricane conditions

A recently proposed model of foam impact on the air–sea drag coefficient C _d has been employed for the estimation of the effective foam-bubble radius R _b variation with wind speed U₁₀ in hurricane conditions. The model relates C _d (U₁₀) with the effective roughness length Z _eff (U₁₀) represented as a sum of aerodynamic roughness lengths of the foam-free and foam-covered sea surfaces Z _w (U₁₀) and Z _f (U₁₀) weighted with the foam coverage coefficient α _f (U₁₀). This relation is treated for known phenomenological distributions C _d (U₁₀), Z _w (U₁₀), and α _f (U₁₀) at strong wind speeds as an inverse problem for the effective roughness parameter of foam-covered sea surface Z _f (U₁₀). The present study is aimed at the estimation of the effective roughness of the sea surface assuming that the measurement data for the effective drag coefficient are known. The effective foam-bubble size is found as a function of the wind speed.     

The generalized kinetic equation as a model for the nonlinear transfer in third-generation wave models

An alternative model for the nonlinear interaction term S_nl in spectral wave models, the so called generalized kinetic equation (Janssen J Phys Oceanogr 33(4):863–884, 2003; Annenkov and Shrira J Fluid Mech 561:181–207, 2006b; Gramstad and Stiassnie J Fluid Mech 718:280–303, 2013), is discussed and implemented in the third generation wave model WAVEWATCH-III. The generalized kinetic equation includes the effects of near-resonant nonlinear interactions, and is therefore able, in theory, to describe faster nonlinear evolution than the existing forms of S_nl which are based on the standard Hasselmann kinetic equation (Hasselmann J Fluid Mech 12:481–500, 1962). Numerical simulations with WAVEWATCH have been carried out to thoroughly test the performance of the new form of S_nl, and to compare it to the existing models for S_nl in WAVEWATCH; the DIA and WRT. Some differences between the different models for S_nl are observed. As expected, the DIA is shown to perform less well compared to the exact terms in certain situations, in particular for narrow wave spectra. Also for the case of turning wind significant differences between the different models are observed. Nevertheless, different from the case of unidirectional waves where the generalized kinetic equation represents a obvious improvement to the standard forms of S_nl (Gramstad and Stiassnie 2013), the differences seems to be less pronounced for the more realistic cases considered in this paper.