Palynostratigraphy and palaeoenvironments around the Albian-Cenomanian boundary interval(OAEld),North Bulgaria

The palynological assemblages from two Albian-Cenomanian boundary sections in North Bulgaria are described.The samples analyzed yielded a diverse palynological content including dinoflagellate cysts and miospores.Based on dinocyst nutrient and productivity indices a phase of enhanced nutrient availability and high primary productivity is inferred for the latest Albian interval.The pronounced predominance of peridinioid dinocysts in this interval,namely O.verrucosum,O.scabrosum and especially P.infusorioides is considered to reflect eutrophic conditions.It coincides with the increased phosphorus mass accumulation occurring at the top part of the Upper Albian Dekov Formation.OAE 1d is indicated in the Tolovitsa karst spring section,based on palynofacies dominated by high amounts of granular amorphous organic matter(AOM)related to anoxic environmental conditions.These sections serve as evidence suggesting a relationship between Cretaceous peridinioid cysts(including Palaeohystichphora infusorioides,Ovoidinium verrucosum,O.scabrosum)and anoxic/suboxic conditions and/or high primary productivity.The pollen spectrum inferred relatively stable vegetation patterns of surrounding continental areas during and after the Albian/Cenomanian boundary interval and the times of OAE 1d formation.The hinterland vegetation integrated mainly pteridophyte spores and gymnosperms.The area was part of the Southern Laurasian floral province which was characterized by warm temperate to subtropical humid climate.Angiosperms were still minor part of this vegetation.     

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.     

Precursor times of abnormal <Emphasis Type="Italic">b</Emphasis>-values prior to mainshocks

Seismic observations exhibit the presence of abnormal b-values prior to numerous earthquakes. The time interval from the appearance of abnormal b-values to the occurrence of mainshock is called the precursor time. There are two kinds of precursor times in use: the first one denoted by T is the time interval from the moment when the b-value starts to increase from the normal one to the abnormal one to the occurrence time of the forthcoming mainshock, and the second one denoted by T _p is the time interval from the moment when the abnormal b-value reaches the peak one to the occurrence time of the forthcoming mainshock. Let T* be the waiting time from the moment when the abnormal b-value returned to the normal one to the occurrence time of the forthcoming mainshock. The precursor time, T (usually in days), has been found to be related to the magnitude, M, of the mainshock expected in a linear form as log(T)?=?q?+?rM where q and r are the coefficient and slope, respectively. In this study, the values of T, T _p , and T* of 45 earthquakes with 3?≤?M?≤?9 occurred in various tectonic regions are compiled from or measured from the temporal variations in b-values given in numerous source materials. The relationships of T and T _p , respectively, versus M are inferred from compiled data. The difference between the values of T and T _p decreases with increasing M. In addition, the plots of T*/T versus M, T* versus T, and T* versus T-T* will be made and related equations between two quantities will be inferred from given data.     

Resources and environmental systems management under synchronic interval uncertainties

Resources and environmental systems management (RESM) is challenged by the synchronic effects of interval uncertainties in the related practices. The synchronic interval uncertainties are misrepresented as random variables, fuzzy sets, or interval numbers in conventional RESM programming techniques including stochastic programming. This may lead to ineffectiveness of resources allocation, high costs of recourse measures, increased risks of unreasonable decisions, and decreased optimality of system profits. To fill the gap of few corresponding studies, a synchronic interval linear programming (SILP) method is proposed in this study. The proposition of interval sets and interval functions and coupling them with linear programming models lead to development of an SILP model for RESM. This enables incorporation of interval uncertainties in resource constraints and synchronic interval uncertainties in the programming objective into the optimization process. An analysis of the distribution-independent geometric properties of the feasible regions of SILP models results in proposition of constraint violation likelihoods. The tradeoff between system optimality and constraint violation is analyzed. The overall optimality of SILP systems under synchronic intervalness is quantified through proposition of integrally optimal solutions. Integration of these efforts leads to a violation-constrained interval integral method for optimization of RESM systems under synchronic interval uncertainties. Comparisons with selected existing methods reveal the effectiveness of SILP at eliminating negativity of synchronic intervalness, enabling risk management of and achieving overall optimality of RESM systems, and enhancing the reliability of optimization techniques for RESM problems. The exploited framework for analyzing synchronic interval uncertainties in RESM systems is helpful for addressing synchronisms of other uncertainties such as randomness or fuzziness and avoiding the resultant decision mistakes and disasters due to neglecting them.     

Holocene vegetation succession and responses to climate change in the northern sector of Northeast China

Sediment pollen samples from the Huola Basin in the northern sector of northeast China, and surface pollen samples from its environs, were analyzed to reconstruct accurately the historical response of vegetation to climate change since 9100 cal yr BP. Pollen analysis of the Huola Section indicates that vegetation experienced a transformation from early-mid Holocene warm-cold mixed vegetation to late Holocene cold-temperate vegetation. From 9100 to 6000 cal yr BP, the study area was warmer and moister than at present, developing Corylus, Carpinus, Pinus, Picea, Betula and Larix-dominated forests. Two cooling events at 6000–5000 and 3500–2500 cal yr BP led to a decrease in Corylus, Carpinus and other warmth-loving vegetation, whereas cold temperate forests composed of Larix and Betula expanded. After 2500 cal yr BP, Larix and Betula dominated cold-temperate vegetated landscapes. The Holocene warm period in NE China(9100–6000 cal yr BP) suggests that such warming could have resulted in a strengthening of the influence from East Asian Summer Monsoon on northernmost NE China and would have benefited the development of warm-temperate forest vegetation and an improved plant load, which also provides the similarity model for the possible global warming in the future.     

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.     

The maximum expected earthquake magnitudes in different future time intervals of six seismotectonic zones of Iran and its surroundings

One of the crucial components in seismic hazard analysis is the estimation of the maximum earthquake magnitude and associated uncertainty. In the present study, the uncertainty related to the maximum expected magnitude μ is determined in terms of confidence intervals for an imposed level of confidence. Previous work by Salamat et al. (Pure Appl Geophys 174:763-777, 2017) shows the divergence of the confidence interval of the maximum possible magnitude m_max for high levels of confidence in six seismotectonic zones of Iran. In this work, the maximum expected earthquake magnitude μ is calculated in a predefined finite time interval and imposed level of confidence. For this, we use a conceptual model based on a doubly truncated Gutenberg-Richter law for magnitudes with constant b-value and calculate the posterior distribution of μ for the time interval T_f in future. We assume a stationary Poisson process in time and a Gutenberg-Richter relation for magnitudes. The upper bound of the magnitude confidence interval is calculated for different time intervals of 30, 50, and 100 years and imposed levels of confidence α?=?0.5, 0.1, 0.05, and 0.01. The posterior distribution of waiting times T_f to the next earthquake with a given magnitude equal to 6.5, 7.0, and 7.5 are calculated in each zone. In order to find the influence of declustering, we use the original and declustered version of the catalog. The earthquake catalog of the territory of Iran and surroundings are subdivided into six seismotectonic zones Alborz, Azerbaijan, Central Iran, Zagros, Kopet Dagh, and Makran. We assume the maximum possible magnitude m_max?=?8.5 and calculate the upper bound of the confidence interval of μ in each zone. The results indicate that for short time intervals equal to 30 and 50 years and imposed levels of confidence 1???α?=?0.95 and 0.90, the probability distribution of μ is around μ?=?7.16???8.23 in all seismic zones.     

The seismic wave absorption in the crust and upper mantle in the vicinity of the Kislovodsk seismic station

The Q-factor estimates of the Earth’s crust and upper mantle as the functions of frequency (Q(f)) are obtained for the seismic S-waves at frequencies up to ~35 Hz. The estimates are based on the data for ~40 earthquakes recorded by the Kislovodsk seismic station since 2000. The magnitudes of these events are M_W > 3.8, the sources are located in the depth interval from 1 to 165 km, and the epicentral distances range from ~100 to 300 km. The Q-factor estimates are obtained by the methods developed by Aki and Rautian et al., which employ the suppression of the effects of the source radiation spectrum and local site responses in the S-wave spectra by the coda waves measured at a fixed lapse time (time from the first arrival). The radiation pattern effects are cancelled by averaging over many events whose sources are distributed in a wide azimuthal sector centered at the receiving site. The geometrical spreading was specified in the form of a piecewise-continuous function of distance which behaves as 1/R at the distances from 1 to 50 km from the source, has a plateau at 1/50 in the interval from 50–70 km to 130–150 km, and decays as \({\raise0.7ex\hbox{$1$} \!\mathord{\left/ {\vphantom {1 {\sqrt R }}}\right.\kern-\nulldelimiterspace} \!\lower0.7ex\hbox{${\sqrt R }$}}\) beyond 130–150 km. For this geometrical spreading model and some of its modifications, the following Q-factor estimates are obtained: Q(f) ~ 85f^0.9 at the frequencies ranging from ~1 to 20 Hz and Q(f) ~ 75f^1.0 at the frequencies ranging from ~1 to 35 Hz.     

Attenuation coefficients of Rayleigh and <Emphasis Type="Italic">Lg</Emphasis> waves

Analysis of the frequency dependence of the attenuation coefficient leads to significant changes in interpretation of seismic attenuation data. Here, several published surface-wave attenuation studies are revisited from a uniform viewpoint of the temporal attenuation coefficient, denoted by χ. Theoretically, χ( f) is expected to be linear in frequency, with a generally non-zero intercept γ?=?χ(0) related to the variations of geometrical spreading, and slope dχ/df = π/Q _e caused by the effective attenuation of the medium. This phenomenological model allows a simple classification of χ( f) dependences as combinations of linear segments within several frequency bands. Such linear patterns are indeed observed for Rayleigh waves at 500–100-s and 100–10-s periods, and also for Lg from ~2 s to ~1.5 Hz. The Lg χ( f) branch overlaps with similar linear branches of body, Pn, and coda waves, which were described earlier and extend to ~100 Hz. For surface waves shorter than ~100 s, γ values recorded in areas of stable and active tectonics are separated by the levels of \(\gamma _{D} \approx 0.2 \times 10^{-3}\) s^???1 (for Rayleigh waves) and 8 ×10^???3 s^???1 (for Lg). The recently recognized discrepancy between the values of Q measured from long-period surface waves and normal-mode oscillations could also be explained by a slight positive bias in the geometrical spreading of surface waves. Similarly to the apparent χ, the corresponding linear variation with frequency is inferred for the intrinsic attenuation coefficient, χ _i , which combines the effects of geometrical spreading and dissipation within the medium. Frequency-dependent rheological or scattering Q is not required for explaining any of the attenuation observations considered in this study. The often-interpreted increase of Q with frequency may be apparent and caused by using the Q-based model of attenuation and following preferred Q( f) dependences while ignoring the true χ( f) trends within the individual frequency bands.     

Shear Wave Velocity Estimates through Combined Use of Passive Techniques in a Tectonically Active Area

We made an attempt to assess the shear wave velocity values V _S and, to a lesser extent, the V _P values from ambient noise recordings in an array configuration. Five array sites were situated in the close proximity to borehole sites. Shear wave velocity profiles were modeled at these five array sites with the aid of two computational techniques, viz. spatial autocorrelation (SPAC) and H/V ellipticity. Out of these five array sites, velocity estimates could be reliably inferred at three locations. The shear wave velocities estimated by these methods are found to be quite consistent with each other. The computed V _S values up to 30 m depth are in the range from 275 to 375 m/s in most of the sites, which implies prevalence of a low velocity zone at some pocket areas. The results were corroborated by evidence of site geology as well as geotechnical information.     

Gamma ray spectra of some lavas from Vesuvius

The authors have carried out the gamma-ray spectrometry of some Vesuvius lavas, by means of a S.E.L.O. equipment which consists of a 100 channels analyzer joined to a photomultiplier tube annexed to aNaI(T1) phosphor. TheU ²³⁸,Th ²³²,Ra ²²⁶ andK ₂ O contents have been determined and the following conclusions have been drawn: theU ²³⁸ seems to be present in the examined lavas and the contents are fairly equal within the errors; also theTh ²³² andK ₂ O contents show no sensible variations in the different lavas. TheRa ²²⁶ content is remarkably higher than that necessary for the equilibrium with the observedU ²³⁸ content and shows a decrease with the age of the lavas. TheRa ²²⁶ content which exceeds that necessary for the equilibrium with the observedU ²³⁸ content decreases with a law similar to the natural radioactive decay ofRa ²²⁶.     

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.     

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.     

Behavior of Parameters of Enhancements in the Nighttime Electron Concentration in the Ionospheric <Emphasis Type="Italic">F</Emphasis>2 Layer

We have analyzed the behavior of the F2 layer parameters during nighttime periods of enhanced electron concentration by the results of vertical sounding of the ionosphere carried out with five-minute periodicity in Almaty (76°55′ E, 43°15′ N) in 2001–2012. The results are obtained within the frameworks of the unified concept of different types of ionospheric plasma disturbances manifested as variations in the height and half-thickness of the layer accompanied by an increase and decrease of N _m F2 at the moments of maximum compression and expansion of the layer. A good correlation is found between height h _Am , which corresponds to the maximum increase, and layer peak height h _m F, while h _Am is always less than h _m F. The difference between h _Am and h _m F linearly increases with increasing h _m F. Whereas the difference is ~38 km for h _m F = 280 km, it is ~54 km for h _m F = 380 km. Additionally, the correlation is good between the increase in the electron concentration in the layer maximum ΔN _m and the maximum enhancement at the fixed height ΔN; the electron concentration enhancement in the layer maximum is about two to three times lower than its maximum enhancement at the fixed height.     

Influence of Solar Wind Plasma Parameters on the Intensity of Isolated Magnetospheric Substorms

Parameters of the interplanetary magnetic field and solar wind plasma during periods of 163 isolated substorms have been studied. It is shown that the solar wind velocity V and plasma density N remain approximately constant for at least 3 h before substorm onset Т _o and 1 h after Т _o . On average, the velocity of the solar wind exhibits a stable trend toward anticorrelation with its density over the whole data array. However, the situation is different if the values of V and N are considered with respect to the intensity of substorms observed during that period. With the growth of substorm intensity, quantified as the maximum absolute value of AL index, an increase in both the solar wind plasma velocity and density, at which these substorms appear, is obsreved. It has been found that the magnitude of the solar wind dynamic pressure P is closely related to the magnetosphere energy load defined as averaged values of the Kan–Lee electric field E_KL and Newell parameter dΦ/dt averaged for 1 h interval before Т _o . The growth of the dynamic pressure is accompanied by an increase in the load energy necessary for substorm generation. This interrelation between P and values of EKL and dΦ/dt is absent in other, arbitrarily chosen periods. It is believed that the processes accompanying increasing dynamic pressure of the solar wind result in the formation of magnetosphere conditions that increasingly impede substorm generation. Thus, the larger is P, the more solar wind energy must enter the Earth’s magnetosphere during the period of the growth phase for substorm generation. This energy is later released during the period of the substorm expansion phase and creates even more intense magnetic bays.     

Estimation of key surface parameters in semi-arid region and their impacts on improvement of surface fluxes simulation

Uncertainties in some key parameters in land surface models severely restrict the improvement of model capacity for successful simulation of surface-atmosphere interaction. These key parameters are related to soil moisture and heat transfer and physical processes in the vegetation canopy as well as other important aerodynamic processes. In the present study, measurements of surface-atmosphere interaction at two observation stations that are located in the typical semi-arid region of China, Tongyu Station in Jilin Province and Yuzhong Station in Gansu Province, are combined with the planetary boundary layer theory to estimate the value of two key aerodynamic parameters, i.e., surface roughness length z_0m and excess resistance κB^-1. Multiple parameterization schemes have been used in the study to obtain values for surface roughness length and excess resistance κB^-1 at the two stations. Results indicate that z_0m has distinct seasonal and inter-annual variability. For the type of surface with low-height vegetation, there is a large difference between the default value of z_0m in the land surface model and that obtained from this study. κB^-1 demonstrates a significant diurnal variation and seasonal variability. Using the modified scheme for the estimation of z_0m and κB^-1 in the land surface model, it is found that simulations of sensible heat flux over the semi-arid region have been greatly improved. These results suggest that it is necessary to further evaluate the default values of various parameters used in land surface models based on field measurements. The approach to combine field measurements with atmospheric boundary layer theory to retrieve realistic values for key parameters in land surface models presents a great potential in the improvement of modeling studies of surface-atmosphere interaction.     

Diagnostics of Turbulent Dynamo from the Flux Emergence Rate in Solar Active Regions

Line-of-sight magnetograms acquired by the Helioseismic and Magnetic Imager (HMI) onboard the Solar Dynamic Observatory (SDO) and by the Michelson Doppler Imager (MDI) onboard the Solar and Heliospheric Observatory (SOHO) for 14 emerging ARs were used to study the derivative of the total unsigned flux–the flux emergence rate, R(t). We found that the emergence regime is not universal: each AR displays a unique emergence process. Nevertheless, two types of the emergence process can be identified. First type is a “regular” emergence with quasi-constant behavior of R(t) during a 1–3 day emergence interval with a rather low magnitude of the flux derivative, R_max = (0.57 ± 0.22) × 10²² Mx day^–1. The second type can be described as “accelerated” emergence with a long interval (>1 day) of the rapidly increasing flux derivative R(t) that result in a rather high magnitude of R_max= (0.92 ± 0.29) × 10²² Mx day^–1, which later changes to a very short (about a one third of day) interval of R(t) = const followed by a monotonous decrease of R(t). The first type events might be associated with emergence of a flux tube with a constant amount of flux that rises through the photosphere with a quasi-constant speed. Such events can be explained by the traditional largescale solar dynamo generating the toroidal flux deep in the convective zone. The second-type events can be interpreted as a signature of sub-surface turbulent dynamo action that generates additional magnetic flux (via turbulent motions) as the magnetic structure makes its way up to the solar surface.     

Seasonal variation of soil respiration under different land use/land cover in arid region

Through combining the soil respiration with the main environmental factors under the planting shelterbelt (Populus woodland) and the natural desert vegetation (Tamarix ramosissima+Phragmites communis community and Haloxylon ammodendron community) in the western Junngar Basin, the difference in soil respiration under different land use/land cover types and the responses of soil respiration to temperature and soil moisture were analyzed. Results showed that the rate of soil respiration increased with temperature. During the daytime, the maximum soil respiration rate occurred at 18:00 for the Populus woodland, 12:00 for T. ramosissima+Ph. communis community, and 14:00 for H. ammodendron community, while the minimum rate all occurred at 8:00. The soil respiration, with the maximum rate in June and July and then declining from August, exhibited a similar trend to the near-surface temperature from May to October. During the growing season, the mean soil respiration rates and seasonal variation differed among the land use/land cover types, and followed the order of Populus woodland > T. ramosissima+Ph. communis community > H. ammodendron community. The difference in the soil respiration rate among different land use/land cover types was significant. The soil respiration of Pouplus woodland was significantly correlated with the near-surface temperature and soil temperature at 10 cm depth (P<0.01) in an exponential manner. The soil respiration of T. ramosissima+Ph. communis and H. ammodendron communities were all linearly correlated with the near-surface temperature and soil surface temperature (P<0.01). Based on the near-surface temperature, the calculated Q₁₀ of Populus woodland, T. ramosissima+Ph. communis community and H. ammodendron community were 1.48, 1.59 and 1.63, respectively. The integrated soil respiration of the three land use/land cover types showed a significant correlation with the soil moisture at 0–5 cm, 5–15 cm and 0–15 cm depths (P<0.01). The quadratic model could best describe the relationship between soil respiration and soil moisture at 0–5 cm depth (P<0.01).     

Spatial variability and long-term change in pollen diversity in Nam Co catchment (central Tibetan Plateau): Implications for alpine vegetation restoration from a paleoecological perspective

Pollen diversity offers abundant clues into the floristic diversity and history of vegetation change. Few palynological studies investigated modern pollen diversity or the past floristic diversity on the Tibetan Plateau (TP). Based on modern pollen assemblages from 37 topsoils and 63 surface lake sediments in the Nam Co catchment on the central TP, this study quantitatively explored spatial distribution of modern pollen diversity using Shannon-Wiener index (H) and palynological richness (E(T _n ), n=600). Pollen diversity indices showed spatial variability among vegetation types, reflecting the differences in terrestrial floristic diversity in the lake catchment. Their values were high in the southeastern region of the lake catchment which is covered by alpine steppe, while values were low for alpine meadow and marsh meadow. The pollen diversity in lacustrine pollen assemblage could be an effective proxy to document past floristic diversity. The past floristic diversity in the lake catchment, recovered from a fossil pollen record of NMLC-1, showed a long-term change of ascending overlaid by several rapid diversity changes during the last 8400 years due to the downward shift of altitudinal vegetation belt driven by a general climatic cooling. The results imply that under the environmental challenge of climate warming and vegetation degradation, alpine vegetation restoration in the Nam Co catchment and the central TP should pay attention to altitudinal vegetation belt and zonal vegetation of alpine steppe, and use the long-term change of floristic diversity as a historical analogue.     

Limits of transversely isotropic elastic parameters for the existence of classical Rayleigh waves

Solutions of P-SV equations of motion in a homogeneous transversely isotropic elastic layer contain a factor exp(±ν _j z), where z is the vertical coordinate and j?=?1, 2. For computing Rayleigh wave dispersion in a multi-layered half space, ν _j is computed at each layer. For a given phase velocity (c), ν _j becomes complex depending on the transversely isotropic parameters. When ν _j is complex, classical Rayleigh waves do not exist and generalised Rayleigh waves propagate along a path inclined to the interface. We use transversely isotropic parameters as α _H , β _V , ξ, ? and η and find their limits beyond which ν _j becomes complex. It is seen that ν _j depends on ? and η, but does not depend on ξ. The complex ν _j occurs when ? is small and η is large. For a given c/β _V , the region of complex ν _j in a ? -η plane increases with the increase of α _H /β _V . Further, for a given α _H /β _V , the complex region of ν _j increases significantly with the decrease of c/β _V . This study is useful to compute dispersion parameters of Rayleigh waves in a layered medium.