A novel non-iterative direct displacement-based seismic design procedure for self-centering buckling-restrained braced frame structures

Self-centering buckling-restrained braces (SCBRBs) were proposed recently to minimize residual deformation of the braces induced by yielding or buckling. Although earthquake resilience of structures equipped with the SCBRBs can be well achieved using displacement based designs (DBDs), previously proposed DBD procedures generally involve iterations. In this study, a novel direct displacement-based design method with a non-iterative procedure, named R–C_R DDBD, is proposed and applied to design of steel braced frame structures with SCBRBs. Unlike previously adopted DBD, the yield displacement does not need to be assumed initially in the proposed procedure. Instead, the yield strength and yield displacement are determined directly by the predetermined objective drift (ratio), using the relation of the strength reduction factor (R) and constant-strength inelastic displacement ratio spectra (C_R spectra), i.e. the R–C_R relation. Since the derived R–C_R relation is independent with the peak ground acceleration of the earthquake records when stiffness and strength degradation are not considered, the proposed procedure can be accurate for any seismic level. The R–C_R DDBD is supposed to begin with the knowledge of the seismic excitation level (according to the structure category, site classification and owner’s requirements) and the corresponding target drift; the end of the design is to obtain the cross sections of main frame members and all the bracing parameters. The result of two 7-story buildings designed according to the R–C_R DDBD procedure demonstrates that this procedure can be effective and fairly simple for practical seismic design.     

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.     

Thermoremanent and chemical magnetization of exsolution products of nanosized titanomagnetites

Based on the theory of two-phase interacting nanoparticles, the formation of thermoremanent and chemical remanent magnetization in nanosized titanomagnetites is modeled. It is shown that the value of thermoremanent magnetization barely depends on the degree of titanomagnetite exsolution whereas, chemical remanent magnetization which emerges during the exsolution increases up to at most the value of thermoremanent magnetization. The values of the ratio of thermoremanent to ideal magnetization, R _t , fall within the limits 0.8 ≤ R _t ≤ 1. The analogous ratio of chemical remanent magnetization to the ideal R _c are below R _t at all stages of the exsolution. Besides, the magnetic interaction between the nanoparticles reduces the values of thermoremanent and chemical magnetization but barely affects the ratio.     

An adaptive modal pushover analysis procedure (<Emphasis Type="Italic">VMPA</Emphasis>-<Emphasis Type="Italic">A</Emphasis>) for buildings subjected to bi-directional ground motions

A new modal pushover analysis procedure (VMPA-A) is developed and implemented in MATLAB code for three-dimensional buildings subjected to bidirectional ground motions. VMPA-A uses stepwise force patterns to represent changes in the dynamic characteristics because of the accumulated structural damages. The hybrid-spectrum concept is introduced to account for the bidirectional ground motion effects. Due to enactments of the equal displacement rule and the secant stiffness-based linearization process, nonlinear analysis is performed for specific displacement targets without stipulation of full modal capacity curves for each mode. Horizontal components of an earthquake record are considered simultaneously, and the consistency between the force and displacement vectors for each mode is provided. These are the main advantages of the proposed procedure against modal pushover analysis (MPA). An existing 21-story reinforced concrete building is analyzed to exemplify VMPA-A. The response parameters such as displacements, story drifts, internal forces, strains, etc. are discussed by comparing the results of VMPA-A with nonlinear time history analyses, which is accepted as the “exact solution”. Though consistent demand estimations are obtained for story drifts, displacements and deformations, some conservative results are obtained for story shears.     

The noble gas isotope geochemical composition of chert at the bottom of Cambrian in Tarim Basin,China

In the Tarim Basin, black shale series at the bottom of Cambrian is one of the important marine facies hydrocarbon source rocks. This research focuses on the analysis of the isotope of noble gas of 11 cherts. The R/R _a ratio of chert in the Keping area is 0.032–0.319, and ⁴⁰Ar/³⁶Ar is 338–430. In Quruqtagh the R/R _a ratio is 0.44–10.21, and ⁴⁰Ar/³⁶Ar is 360–765. The R/R _a ratio of chert increases with ⁴⁰Ar/³⁶Ar from the west to the east accordingly. They have evolved from the crust source area to the mantle source area in a direct proportion. Surplus argon ⁴⁰Ar_E in chert is in direct proportion to the R/R _a ratio, indicating that it has the same origin of excess argon as in fluid and mantle source helium. Comparison of the R/R _a ratios between the west and the east shows that the chert in the eastern part formed from the activity system of the bottom hydrothermal venting driven by the mantle source, where the material and energy of crust and mantle had a strong interaction in exchange; whereas in the western part, chert deposited from the floating of hydrothermal plume undersea bottom, which is far away from the centre of activities of the hydrothermal fluid of ocean bottom. In addition, from noble gas isotope composition of chert, it is suggested that the ocean anoxia incident happened at the black shale of the Cambrian bottom probably because of the large-scaled ocean volcanoes and the following hydrothermal activities.     

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.     

Geomagnetically Induced Currents and Space Weather: Pi3 Pulsations and Extreme Values of Time Derivatives of the Geomagnetic Field’s Horizontal Components

The problem of estimating the time derivatives of the horizontal components of the geomagnetic field and forecasting the probability of the occurrence of perturbations that exceed a given threshold level (the over-threshold perturbations) arises in the applications concerned with the geomagnetically induced currents (GICs). In this work, we consider the temporal and spatial structure of the Pi3 pulsations with quasi-periods of 10² to 10³ s during which the auroral and subauroral stations of the IMAGE network record over-threshold values in the derivatives of the meridional (along the longitudinal circle) B_X component and latitudinal (along the latitudinal circle) B_Y component. The extreme |dB_X/dt| values mainly develop against the background of the Pi3 pulsations with a complex frequency content, whereas the extreme |dB_Y/dt| values appear when the buildup (decay) phases of the bay-like disturbance associated with the evolution of a substorm coincide with the respective phases of the field of pulsations. The conditions under which the derivatives |dB_X/dt| and |dB_Y/dt| reach their over-threshold values are studied for subauroral latitudes by the technique of superposed epoch analysis. The extreme values of the derivatives most frequently occur during the main phase of moderate magnetic storms or beyond the storm—during high substorm activity under the conditions of a negative vertical component of the interplanetary magnetic field. The probability of the occurrence of over-threshold values increases at high amplitudes of the Pi3 pulsations and depends on their spectral content. The problem of analyzing and forecasting the over-threshold |dB_Y/dt| perturbations is complicated by the fact that the scale of the perturbations is small along the lines of latitude and large along the meridians. This can result in GIC excitation in the North–South oriented electric power lines by the geomagnetic perturbations localized within a narrow band in longitude which can be missed during the measurements.     

Force reduction factors for steel buildings with welded and post-tensioned connections

The seismic responses of steel buildings with perimeter moment resisting frames (MRF) with welded connections (WC) are estimated and compared to those of similar buildings with semi-rigid post-tensioned connections (PC). The responses are estimated in terms of ductility reduction factors (R _µ,), ductility demands (µ _G ) and force reduction factors (R). Two steel model buildings, which were modeled as complex-3D-MDOF systems, were used in the study. Results indicate that the reduction magnitude of global response parameters is larger than that of local response parameters, contradicting the same reduction implicitly assumed in the static equivalent lateral force procedure, implying that non-conservative design may result. The value of 8 for R, suggested in many codes for ductile steel MRF, and the value of 1 suggested in the well known Newmark and Hall procedure for the ratio of R to µ _G , cannot be justified. The reason for this is that SDOF systems were used to model actual structures, where higher mode effects, energy dissipation and structural overstrength weren’t explicitly considered. The codes should be more transparent regarding the magnitude and the components involved in the force reduction factors. The seismic performance of steel buildings with PC may be superior to that of the buildings with WC, since their force reduction factors are larger and their ductility demands smaller, implying that PC buildings could be designed for smaller lateral seismic forces. The conclusions of this paper are for the particular structural systems and models considered. Much more research is needed to reach more general conclusions.     

Critical frequencies <Emphasis Type="Italic">foF</Emphasis>2 as an indicator of trends in thermospheric dynamics

The time variations in three parameters during the last decades are considered. R(foF2) is the correlation coefficient between the nighttime and daytime values of foF2 for the same day. Stable trends are found for the minimum (R(foF2)(max)) and maximum (R(foF2)(min)) values of R(foF2) during a year. The foF2(night)/foF2(day) ratio demonstrates both, negative and positive trends, and the trend sign depends on the inclination I and declination D of the magnetic field. The correlation coefficient r(h, fo) between foF2 and the 100 hP level in the stratosphere demonstrates a decrease (in the years of maximum and minimum solar activity) from the 1980s to the 1990s. The trends in all three groups of data are considered under the assumption of long-term changes in the circulation in the upper atmosphere.     

Distribution of the field−aligned currents in the ionosphere: dawn–dusk asymmetry and its relation to the asymmetry between the two hemispheres

The dynamics of the system of field-aligned currents (FACs) and closing ionospheric Pedersen currents is considered with the use of original processing methods and the data from four substorms of different types. The total current system comprises of two parts. One is the well-known substorm current wedge (SCW) system, in which the zonal (westward ) current closes FACs in the R1 zone (region). The component 2 consists of two pairs of meridional currents flowing equatorward and poleward in the R1 region and creating regions R0 and R2 (according to the classification of Iijima and Potemra). It is shown that the total FAC of the disturbed magnetosphere–ionosphere system is dominated by the contribution of component 2, which contradicts the original version of the SCW model but is consistent with new data. The quantitative characteristics of the dawn–dusk asymmetry are determined for the FAC distribution in the ionosphere for each substorm. It is shown that the ratio of the average intensities of FACs in the regions R0 and R2 was I_R0/I_R2 ≥ 0.4, which contradicts the popular opinion that there are no FACs in the polar cap. In addition, a relatively rare event of the simultaneous start of the substorm explosive phase and the SSC caused by the dynamic impact of the solar wind when the polar cap expands (rather than compresses as usual) is considered.     

Geochemical characteristics of light hydrocarbons in cracking gases from chloroform bitumen A,crude oil and its fractions

The composition characteristics of light hydrocarbons from crude oil, chloroform bitumen A, saturated hydrocarbon fraction, aromatic hydrocarbon fraction, and asphaltene fraction during cracking have been studied systematically. The results revealed that the content of n-alkanes, branched alkanes and cycloalkanes in light hydrocarbons from the samples gradually decreased as the simulation temperature increased, and finally almost depleted completely, while the abundance of methane, benzene and its homologues increased obviously and became the main products. The ratios of benzene/ n-hexane and toluene/n-heptane can be used as measures for oil cracking levels. Variation characteristics of maturity parameters of light hydrocarbons, for example, iC₄/nC₄, iC₅/nC₅, isoheptane value, 2,2-DMC₄/nC₆, and 2-MC₆+3-MC₆/nC₇ for different samples with increasing pyrolysis temperature, are consistent with those in petroleum reservoirs, indicating that these parameters may be efficient maturity index.     

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.     

The partitioning of radiated energy and the largest aftershock of seismic sequences occurred in the northeastern Italy and western Slovenia

We analyzed the most relevant seismic sequences that occurred from 1977 to 2007 in the Friuli-Venezia Giulia region (northeastern Italy) and western Slovenia. The eight aftershock sequences were triggered by low- to moderate-magnitude earthquakes with mainshock duration magnitude ranging from 3.7 to 5.6. The b-value of the Gutenberg–Richter law varies from 0.8 to 1.1. The modified Omori’s modeling of the sequences evidences values of the p exponent ranging from 0.8 to 1.0. Using the Reasenberg and Jones (Science 243:1173–1176, 1989; Science 265:1251–1252, 1994) approach, we computed the probabilistic estimate of the aftershock rates and the largest aftershock in given time intervals. The difference in magnitude between the mainshock and the largest aftershock is calculated according to the modified Båth law and using an approach that considers the partitioning of the radiated seismic energy between mainshock and aftershocks. The partitioning of the radiated seismic energy appears to play a significant role in the evolution of the sequences. We define the parameter R _ES as the ratio between the radiated seismic energy of the mainshock and the summation of the seismic energy radiated by the aftershocks. The difference in magnitude between the mainshock and the largest aftershock, calculated with the parameter R _ES, agrees well with the observed difference. In most sequences, the parameter R _ES decreases very quickly until the occurrence of the largest aftershock and then becomes constant. By analyzing the values of R _ES during the early hours following the mainshock, we found that the R _ES values after 24 h are well related to the final ones, calculated on the whole sequence, and to the differences in magnitude between the mainshock and the largest aftershock.     

Experimental examination on the heterogeneity parameter <Emphasis Type="Italic">C</Emphasis><Subscript>v</Subscript> of earthquake precursors

Two rock samples with different structures and materials were deformed under a biaxial loading system, and multipoint strain measurements were performed for each sample. The distribution of strain anomalies during the deformation and the instability process were analyzed by using C _v value put forward by WANG Xiao-qing and CHEN Xue-zhong, et al, a parameter to describe the heterogeneous distribution of earthquake precursors, so as to examine the method of C _v value and to explore its physical meaning experimentally. The result shows that the change of C _v value is correlated to the change of deformation characteristics and is an effective parameter to describe the heterogeneity of precursor distribution. C _v value increases firstly and then decreases before the instability, and the instability occurs when C _v value decreases to the level before increasing. This indicates that C _v value may be a useful parameter for earthquake prediction.     

Enhanced nonlinear static analysis with the drift pushover procedure for tall buildings

In this study a new method for nonlinear static analysis based on the relative displacements of stories is proposed that is able to be implemented in a single stage analysis and considers the effects of an arbitrary number of higher modes. The method is called the extended drift pushover analysis procedure (EDPA). To define the lateral load pattern, values of the relative displacements of stories are calculated using the elastic modal analysis and the modal combination factors introduced. For determining the combination factors, six different approaches are examined. Buildings evaluated in this study consist of four special steel moment-resisting frames with 10–30 stories. Responses including relative displacements of stories, story shear forces and rotation of plastic hinges in each story are calculated using the proposed approaches in addition to modal pushover analysis and nonlinear dynamic time history analyses. The nonlinear dynamic analysis is implemented using ten consistent earthquake records that have been scaled with regard to ASCE7-10. Distribution of response errors of story shears and plastic hinge rotations show that a major part of error corresponds to the second half of the buildings studied. Thus, the mentioned responses are corrected systematically. The final results of this study show that implementing the EDPA procedure using the third approach of this research is able to effectively overcome the limitations of both the traditional and the modal pushover analyses methods and predict the seismic demands of tall buildings with good accuracy.     

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.     

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.     

Shape and position of Earth’s bow shock near-lunar orbit based on ARTEMIS data

Earth’s bow shock is the result of interaction between the supersonic solar wind and Earth’s magnetopause. However, data limitations mean the model of the shape and position of the bow shock are based largely on near-Earth satellite data. The model of the bow shock in the distant magnetotail and other factors that affect the bow shock, such as the interplanetary magnetic field (IMF) B_y, remain unclear. Here, based on the bow shock crossings of ARTEMIS from January 2011 to January 2015, new coefficients of the tail-flaring angle α of the Chao model (one of the most accurate models currently available) were obtained by fitting data from the middle-distance magnetotail (near-lunar orbit, geocentric distance -20R_E>X>-50R_E). In addition, the effects of the IMF B_y on the flaring angle α were analyzed. Our results showed that: (1) the new fitting coefficients of the Chao model in the middle-distance magnetotail are more consistent with the observed results; (2) the tail-flaring angle α of the bow shock increases as the absolute value of the IMF B_y increases. Moreover, positive IMF B_y has a greater effect than negative IMF B_y on flaring angle. These results provide a reference for bow shock modeling that includes the IMF B_y.     

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.     

Seasonal variation of carbon exchange of typical forest ecosystems along the eastern forest transect in China

The long-term and continuous carbon fluxes of Changbaishan temperate mixed forest (CBS), Qianyanzhou subtropical evergreen coniferous forest (QYZ), Dinghushan subtropical evergreen mixed forest (DHS) and Xishuangbana tropical rainforest (XSBN) have been measured with eddy covariance techniques. In 2003, different responses of carbon exchange to the environment appeared across the four ecosystems. At CBS, the carbon exchange was mainly determined by radiation and temperature. 0°C and 10°C were two important temperature thresholds; the former determined the length of the growing season and the latter affected the magnitude of carbon exchange. The maximum net ecosystem exchange (N _EE) of CBS occurred in early summer because maximum ecosystem photosynthesis (G _PP) occurred earlier than maximum ecosystem respiration (R _e). During summer, QYZ experienced severe drought and N _EE decreased significantly mainly as a result of the depression of G _PP. At DHS and XSBN, N _EE was higher in the drought season than the wet season, especially the conversion between carbon sink and source occurring during the transition season at XSBN. During the wet season, increased fog and humid weather resulted from the plentiful rainfall, the ecosystem G _PP was dispressed. The Q ₁₀ and annual respiration of XSBN were the highest among the four ecosystems, while the average daily respiration of CBS during the growing season was the highest. Annual N _EE of CBS, QYZ, DHS and XSBN were 181.5, 360.9, 536.2 and ?320.0 g·C·m^?2·a^?1, respectively. From CBS to DHS, the temperature and precipitation increased with the decrease in latitude. The ratio of N _EE/R _e increased with latitude, while R _e/G _PP, ecosystem light use efficiency (L _UE), precipitation use efficiency and average daily G _PP decreased gradually. However, XSBN usually escaped such latitude trend probably because of the influence of the south-west monsoon climate which does not affect the other ecosystems. Long-term measurement and more research were necessary to understand the adaptation of forest ecosystems to climate change and to evaluate the ecosystem carbon balance due to the complexity of structure and function of forest ecosystems.