Relationship between Net Radiation and Broadband Solar Radiation in the Tibetan Plateau

The characteristics of net radiation (R_n)(0.3--10 μm) in Lhasa and Haibei in the Tibetan Plateau were analyzed based on long-term in-situ measurements of surface radiation data. The monthly average of daily R_n reached a minimum during the winter period followed by an increase until May and then a decline until January. This variation is consistent with solar activity. The annual mean daily total R_n values were 0.92 MJ m^-2 d^-1 and 0.66 MJ m^-2d^-1 in Lhasa and Haibei, respectively. A relationship between R_n and broadband solar radiation (R_s) was demonstrated by a good linear correlation at the two sites. R_n can be an accurate estimate from R_s. The estimated R_n values were similar to the observed values, and the relative deviations between the estimates and measurements of R_n were 2.8% and 3.8% in Lhasa and Haibei, respectively. The application of the R_n estimating model to other locations showed that it could provide acceptable estimated R_n values from the R_s data. Furthermore, we analyzed the influence of clouds on R_n by different clear index (K_s), defined as the ratio of R_s to the extraterrestrial solar irradiance on a horizontal surface. The results indicate that more accurate results are associated with increased cloudy conditions. The influence of the albedo was also considered, but its inclusion in the model resulted in only a slight improvement. Because surface albedo is not usually measured, an expression based solely on global solar radiation could be of more extensive use.     

Longwave incoming radiation in the Tropics: results from field work in three African cities

Summary This study investigates differences in longwave incoming radiation (L↓) within and between three African cities, Dar es Salaam (Tanzania), Ouagadougou (Burkina Faso), and Gaborone (Botswana), during the dry season, and evaluates the performance of a model to simulate these fluxes. In each city, direct observations of L↓, shortwave incoming radiation (K↓), air temperature, air humidity, and total suspended particle (TSP) concentration for three land uses (CBD, green residential, and traditional residential) were taken. The observed L↓ flux decreases with increasing latitude, and temperature becomes an increasingly important factor in governing L↓ variations further from the Equator. Humidity, as well as particle loading, differs significantly between the three cities. Differences between observed and modelled ɛ_sky for rural stations near all cities showed a clear diurnal variation, with maximum differences of 0.08 between day and night. This diurnal difference was incorporated in the model and, for urban areas the model overestimates L↓ by around 25 Wm⁻². However, this model performs equally well regardless of the land use considered in any of the cities. The residual (difference between observed and modelled urban L↓) did not show any correlation with particulate pollution. However, the difference between observed and calculated ɛ_sky is around 0.05 higher in Ouagadougou compared to the other cities, likely due to the heavy dust load observed here. It is concluded that tropical urban longwave radiation is not dramatically different from the mid latitudes.     

Radiation balance in a tropical city — Delhi (India)

The impact of urbanization is assessed by comparing values of the radiation parameters at an urban location with those of a rural site. Urban Delhi was divided according to land-use and the effects of urbanization was studied on incoming short-wave (K), albedo, incoming longwave (L, outgoing longwave (L), and net radiation (Q ^*), were individually studied at four representative sites (Rural, Commercial, Residential and Industrial). MaximumK was observed in the rural and commercial areas whereas highL was observed in the commercial and industrial locations. High depletion ofK of the order of 13% was observed for the industrial location in the winter season. An increase ofL in, the industrial location is of the same order as that of the commercial location, i.e. 20%. The residential location shows quite moderate (4.6%) depletion ofK in comparison with other sites.Nomenclature U Urban - R Rural - K Incoming Short-wave radiation - L Incoming Long-wave radiation - L Outgoing Long-wave radiation - Q ^* Net Radiation - Albedo - K ^* Net Short-wave radiation - L ^* Net Long-wave radiation     

用微气象方法估算淮河流域能量平衡(HUBEX/IOP 1998/99)的统计分析和比较研究

In order to study energy and water cycles in the Huaihe River Basin, micrometeorological measurements were carried out in Shouxian County, Anhui Province, during HUBEX/IOP (May to August 1998 and June to July 1999). The employed techniques included Bowen Ratio-Energy Balance (BREB) and Eddy Covariance (EC) methods. In this paper, the basic characteristics of the energy balance components in the district are analyzed. Furthermore, the results are compared with those from other regions of China.The main results are as follows: (1) There was a consistency between the available energy (Rn-G) and the sum of sensible (H) and latent (E) heat fluxes measured by the EC method (H+E)ec, but Ebr was slightly larger (about 10%) than Eec; (2) Most of the net radiation (Rn) was used to evaporate water from the surface. During HUBEX/IOP in 1998 and 1999, the mean daily amounts of Rn were 13.89 MJ m-2 d-1 and 11.83 MJ m-2 d-1, and the mean Bowen Ratios (β) were 0.14 (over ruderal) and 0.06 (over paddy) respectively; (3) The diurnal variation characteristic of β was larger and unsteady at sunrise and sunset, and smaller and steady during the rest of the daytime. Local advection appeared in the afternoon over paddy areas in 1999; (4) In comparison with the results from other regions of China, the nean β was the lowest (0.06) over paddy areas in the Huaihe River Basin and the highest (0.57) during June-August 1998 in Inner Mongolia grassland. The Bowen Ratio β is mainly related to the soil humidity.     

The Theoretical Model of Atmospheric Turbulence Spectrum in Surface Layer

It is shown that the slope of energy spectrum obtained from the velocity solution of Kdv-Burgers equation lies between -5/3 and -2 in the dilogarithmic coordinates paper. The spectrum is very close to one of Kolmogorov’s isotropic turbulence and Frisch’s intermittent turbulence in inertial region. In this paper, the Kdv-Burgers equation to describe atmospheric boundary layer turbulence is obtained. In the equation, the 1 / Re, corresponds to dissipative coefficient v, to dispersive coefficient β, then (v/ 2β)2 corresponds to .We prove that the wave number corresponding to maximum energy spectrum decreases with the decrease of stability (i.e., the increase of in eddy-containing region. And the spectrim amplitude decreases with the increase of (i.e., the decrease of stability). These results are consistent with actual turbulence spectrum of atmospheric surface layer from turbulence data.     

On the annual and monthly mean diurnal variations of diffuse solar radiation at a meteorological station in west Africa

Summary In addition to global solar radiationE ^g , the hourly diffuse componentE ^d incident on a horizontal surface has been measured from February 1993 to January 1995 at a meteorological station in tropical West Africa. The measured diffuse solar irradiance data was corrected for shadow band effects. The monthly mean diurnal variations of diffuse solar irradiance obtained for identical months in the two years have been compared and found to be generally consistent. The corresponding monthly mean hourly values ofE ^d for identical months in 1993 and 1994 agreed to within 9% while yielding correlation coefficients greater than 0.960. In addition, the monthly mean daily totals ofE ^d for identical months were found to agree mostly to within 6% and showed virtually the same annual variations in both years. The monthly mean daily total values of diffuse solar radiation for most months in the two years ranged between 7.94 MJm^–2d^–1 and 10.50 MJm^–2d^–1. The monthly mean of daily hourly maximum values ofE ^d obtained for identical months in the two years have been discussed in relation to the dominant atmospheric conditions during these months. The results been presented here have been compared with those of some investigators within and outside the Africa region.With 8 Figures     

Surface Energy Balance Measurements Above an Exurban Residential Neighbourhood of Kansas City,Missouri

Previous measurements of urban energy balances generally have been limited to densely built, central city sites and older suburban locations with mature tree canopies that are higher than the height of the buildings. In contrast, few data are available for the extensive, open vegetated types typical of low-density residential areas that have been newly converted from rural land use. We made direct measurements of surface energy fluxes using the eddy-covariance technique at Greenwood, a recently developed exurban neighbourhood near Kansas City, Missouri, USA, during an intensive field campaign in August 2004. Energy partitioning was dominated by the latent heat flux under both cloudy and near clear-sky conditions. The mean daytime Bowen ratio (β) values were 0.46, 0.48, and 0.47 respectively for the cloudy, near clear-sky and all-sky conditions. Net radiation (R _n ) increased rapidly from dawn (−34 and −58W m⁻²) during the night to reach a maximum (423 and 630W m⁻²) after midday for cloudy and near clear-sky conditions respectively. Mean daytime values were 253 and 370W m⁻², respectively for the cloudy and near clear-sky conditions, while mean daily values were 114 for cloudy and 171W m⁻² for near clear-sky conditions, respectively. Midday surface albedo values were 0.25 and 0.24 for the cloudy and near clear-sky conditions, respectively. The site exhibited an angular dependence on the solar elevation angle, in contrast to previous observations over urban and suburban areas, but similar to vegetated surfaces. The latent heat flux (Q _E ), sensible heat flux (Q _H ), and the residual heat storage ΔQ _s terms accounted for between 46–58%, 21–23%, and 18–31% of R _n , respectively, for all-sky conditions and time averages. The observed albedo, R _n , and Q _E values are higher than the values that have been reported for suburban areas with high summer evapotranspiration rates in North America. These results suggest that the rapidly growing residential areas at the exurban fringe of large metropolitan areas have a surface energy balance that is more similar to the rural areas from which they were developed than it is to the older suburbs and city centres that make up the urban fabric to which they are being joined.     

Radiation balance over low-lying and mountainous areas in south-west Germany

Summary  Surface radiative fluxes play a major role in the energy exchange process between the atmosphere and earth surface and are thus very crucial to climatic processes within the atmospheric boundary layer. Based on four years REKLIP (REgio-KLIma-Project) data set of measured radiative fluxes and additional supporting meteorological variables, the surface radiation regime for selected lowland site (Bremgarten 212 m a.s.l.) and mountain sites (Geiersnest at 870 m a.s.l.; Feldberg 1489 m a.s.l.) in the southern Upper Rhine valley region (south-west Germany) has been reported. Time series of radiative fluxes and surface albedo showed significant inter-site differences. Possible reasons for the observed differences have been made. Downward atmospheric radiation A _l at the study sites was parameterised in terms of air temperature, vapour pressure and cloud amount, all of which strongly govern the variation of A _l . Effective terrestrial radiation amounted to about 50% of absorbed shortwave radiation at the study sites annually. During clear sky conditions, global solar irradiance G _s constituted about 76.0% of the incident extraterrestrial solar irradiance at Feldberg mountain site but only 68.5% of that at Bremgarten lowland site. Annual cumulative of net radiative flux R _n amounted to 1722 MJm⁻² yr⁻¹ at the lowland site, while that at Geiersnest and Feldberg mountain sites constituted 84% and 73% respectively of the corresponding magnitude for the lowland site. In the same vein, annual mean of radiation efficiency (defined here as R _n /G _s ) amounted to 0.32 in Feldberg, 0.37 in Geiersnest and 0.41 in Bremgarten. Consequently the annual available energy, of which net radiative flux is representative, was smaller at the mountain ous sites relative to the lowland site during the study period. Inter-annual variability of net radiative flux, its constituent variables and derivatives at the study sites were generally below 10%, with longwave fluxes showing the lowest fluctuation. This renders the measured data quite suitable for modelling purposes. In winter, mean daily sums of R _n showed a slow rise with cloud amount N at the lowland site but a sharp rise with N at Feldberg mountain site. In summer however, mean daily sums of R _n declined significantly with N as well as Linke turbidity factor at the study sites. Received June 24, 1999 Revised November 2, 2000     

Atmospheric turbidity of urban and desert areas of the Nile Basin in the aftermath of Mt. Pinatubo’s eruption

Summary  The Linke T _L , ?ngstr?m β and Unsworth-Monteith δ _a turbidity parameters are investigated for two sites in Egypt: Cairo, a densely populated urban area, and Aswan, an arid unpolluted area. These three turbidity parameters are calculated from broadband pyrheliometric measurements recorded hourly over the period 1992–96. Monthly averages of T _L , β and δ _a show relatively flat and identical seasonal variations with a marked main maxima during spring at both sites, due to Khamsin depressions coming from the Great Sahara. A secondary maximum is observed at Aswan in summer, due to dust haze which prevails during that season, and at Cairo in autumn, due to the northern extension of the Sudan monsoon trough, which is accompanied by small scale depressions with dust particles. Annual mean values of T _L , β and δ _a (5.59, 0.250 and 0.372, respectively) at Cairo are larger than at Aswan (3.89, 0.139 and 0.213, respectively). In the same way, the seasonal mean values of T _L , β and δ _a at Cairo are larger than at Aswan. More generally, the monthly and yearly average turbidity values are significantly larger in Cairo than in Aswan for the whole period 1992–96, which is attributable in part to the urbanization/industrialization effect of Cairo. An estimate of the corresponding overburden is obtained by comparison between the present data and older T _L data from 1922–27. It is also shown that turbidity over both sites is largest during 1992, just after the eruption of Mount Pinatubo in 1991. The dependence of β on some meteorological parameters such as wind speed and direction, precipitable water, relative humidity, temperature and visibility, is also analyzed. This reveals in particular that visibility is not a good predictor of turbidity at either site. Conversely, the wind direction and speed have a definite effect on turbidity, and consequently, largest turbidities occur when the wind carries aerosols from the main industrial particle source areas around Cairo. For any season of the year, the average turbidity at the latter site is larger than that at other big cities such as Athens, Rome, and Toronto, but is lower than at Dhahran, Saudi Arabia. Received February 3, 2000 Revised August 25, 2000     

A K-Theory of Dispersion,Settling and Deposition in the Atmospheric Surface Layer

Dispersion estimates with a Gaussian plume model are often incorrect because of particle settling (β), deposition (γ) or the vertical gradient in diffusivity (K _v (z) = K ₀ + μz). These “non-Gaussian” effects, and the interaction between them, can be evaluated with a new Hankel/Fourier method. Due to the deepening of the plume downwind and reduced vertical concentration gradients, these effects become more important at greater distance from the source. They dominate when distance from the source exceeds L _β = K ₀ U/β ², L _γ  = K ₀ U/γ ² and L _μ = K ₀ U/μ ² respectively. In this case, the ratio β/μ plays a central role and when β/μ = 1/2 the effects of settling and K gradient exactly cancel. A general computational method and several specific closed form solutions are given, including a new dispersion relation for the case when all three non-Gaussian effects are strong. A more general result is that surface concentration scales as C(x) ~ γ ⁻² whenever deposition is strong. Categorization of dispersion problems using β/μ, L _γ and L _μ is proposed.     

Aerodynamic Parameters of Regular Arrays of Rectangular Blocks with Various Geometries

The aerodynamic effects of various configurations of an urban array were investigated in a wind-tunnel experiment. Three aerodynamic parameters characterising arrays—the drag coefficient (C _d ), roughness length (z _o) and displacement height (d)—are used for analysis. C _d is based on the direct measurement of the total surface shear using a floating element, and the other two parameters are estimated by logarithmic fitting of the measured wind profile and predetermined total drag force. The configurations of 63 arrays used for measurement were designed to estimate the effects of layout, wind direction and the height variability of the blocks on these parameters for various roughness packing densities. The results are summarised as follows: (1) The estimated C _d and z _o of the staggered arrays peak against the plan area index (λ _p ) and frontal area index (λ _f ), in contrast with values for the square arrays, which are less sensitive to λ _p and λ _f . In addition, the square arrays with a wind direction of 45° have a considerably larger C _d , and the wind direction increases z _o/H by up to a factor of 2. (2) The effect of the non-uniformity of roughness height on z _o is more remarkable when λ _f exceeds 20%, and the discrepancy in z _o is particularly remarkable and exceeds 200%. (3) The effect of the layout of tall blocks on C _d is stronger than that of short blocks. These results indicate that the effects of both wind direction and the non-uniformity of the heights of buildings on urban aerodynamic parameters vary greatly with λ _p and λ _f ; hence, these effects should be taken into account by considering the roughness packing density.     

Micrometeorological Conditions for Dust-Devil Occurrence in the Atacama Desert

We report on field observations in January 2009 (austral summer) of atmospheric dust devils in the northern part of the Atacama Desert in South America (≈20^◦S). An extremely high level of dust-devil activity over the study site has been observed, dependent on local meteorological conditions. We found a high correlation between the dust-devil frequency of occurrence and the Obukhov length scale, L, calculated from meteorological gradient measurements, with a clear tendency for this frequency to increase with decreasing −L. The upper threshold values of −L ≈ 20–30 m, and the 2-m mean wind speed, V ₂ ≈ 8m s⁻¹, for dust-devil occurrence have been found, but the minimal V ₂ threshold was not observed. Parallel routine meteorological measurements enabled us to calculate the main constituents of the surface energy balance, to obtain direct estimates of the surface albedo (α ≈ 0.21 at the solar noon) and to summarize the local conditions.     

A comparison of some aerodynamic resistance methods using measurements over cotton and grass from the 1991 California ozone deposition experiment

Measurements of dry, deposition velocities (V _d ) of O₃ (using the eddy correlation technique) made over a cotton field and senescent grass near Fresno California during July and August 1991 were used to test some dry deposition velocity models. Over the cotton field, the observed maximum daytimeV _d was about 0.8 cm s^–1 and the average nighttime value was about 0.2 cm s^–1. Over the grass, daytime values averaged about 0.2 cm s^–1 and nighttime values about 0.05 cm s^–1. Application of a site-specific model known as ADOM (Acid Deposition and Oxidant Model) over the cotton field generally overestimated the observations except for a few hours in the afternoon when the observations were underestimated The overestimation was attributed to inadequacies in the surface resistance formulation and the underestimation to uncertainties in the aerodynamic formulation. Unlike previous studies which focused on the role of surface resistance, we perform additional tests using a large variety of aerodynamic resistance formulae, in addition to those in ADOM, to determine their influence on the modelledV _d of O₃ over, cotton. Over grass, ADOM considerably overestimated the observations but showed improvement when other surface resistance formulations were applied.     

Atmospheric inorganic nitrogen in marine aerosol and precipitation and its deposition to the North and South Pacific Oceans

Aerosol and rain samples were collected between 48°N and 55°S during the KH-08-2 and MR08-06 cruises conducted over the North and South Pacific Ocean in 2008 and 2009, to estimate dry and wet deposition fluxes of atmospheric inorganic nitrogen (N). Inorganic N in aerosols was composed of ~68% NH₄⁺ and ~32% NO₃^– (median values for all data), with ~81% and ~45% of each species being present on fine mode aerosol, respectively. Concentrations of NH₄⁺ and NO₃⁻ in rainwater ranged from 1.7–55 μmol L⁻¹ and 0.16–18 μmol L⁻¹, respectively, accounting for ~87% by NH₄⁺ and ~13% by NO₃⁻ of total inorganic N (median values for all data). A significant correlation (r = 0.74, p < 0.05, n = 10) between NH₄⁺ and methanesulfonic acid (MSA) was found in rainwater samples collected over the South Pacific, whereas no significant correlations were found between NH₄⁺ and MSA in rainwater collected over the subarctic (r = 0.42, p > 0.1, n = 6) and subtropical (r = 0.33, p > 0.5, n = 6) western North Pacific, suggesting that emissions of ammonia (NH₃) by marine biological activity from the ocean could become a significant source of NH₄⁺ over the South Pacific. While NO₃⁻ was the dominant inorganic N species in dry deposition, inorganic N supplied to surface waters by wet deposition was predominantly by NH₄⁺ (42–99% of the wet deposition fluxes for total inorganic N). We estimated mean total (dry + wet) deposition fluxes of atmospheric total inorganic N in the Pacific Ocean to be 32–64 μmol m⁻² d⁻¹, with 66–99% of this by wet deposition, indicating that wet deposition plays a more important role in the supply of atmospheric inorganic N than dry deposition.     

Variation on PAR to global solar radiation ratio along altitude gradient in Naeba Mountain

Summary We investigated the ratio of photosynthetically active photon flux (Q _p ) to global solar radiation (R _s ) at three sites along different altitudes in Naeba Mountain, Japan at various temporal scales based on 3 years measurement data (1999–2001). The lowest values of the ratio ever reported were found for all sites on both an hourly and a daily scale. A similar slight diurnal pattern was found for all sites based on the monthly mean hourly values of the ratio. However, different sites exhibited different seasonal courses. Statistically significant altitude dependencies were found for the hourly Q _p /R _s under both clear and cloudy weather conditions (t-test, P < 0.001). For clear weather conditions, the hourly Q _p /R _s exhibited an increasing trend with altitude at an average rate of 3.6% per km. The increasing rate was more noticeable below 900 m (8.8% per km) than above (0.7% per km). The inverse trend was found for hourly Q _p /R _s with the altitude under cloudy weather conditions. The hourly Q _p /R _s decreased from 550 m to 1500 m at a rate of 1.8% per km. Again, a major decrease occurred below 900 m, which had the rate of 4.2% per km, compared with 0.2% per km over 900 m. Although the same tendencies were noted for daily Q _p /R _s , under clear sky conditions, they were not as statistically significant as the hourly counterpart (t-test, P < 0.021). The increasing rate of Q _p /R _s at this scale under clear weather conditions was near that of the hourly rate, but below the 900 m rate was reduced to near half of the hourly rate (4.7% per km). And the rate over 900 m increased to 2.7% per km. On the other hand, statistically significant altitude effect was noted for the daily under cloudy weather conditions (t-test, P < 0.002). A rapidly decreasing rate was found for it along the altitude gradient. The reason was due to the large decreasing rate below 900 m (12.4% per km). But again similar to the hourly Q _p /R _s , a very small decreasing rate of daily Q _p /R _s was found over 900 m under cloudy weather conditions. These results suggest the necessity of considering the altitude dependency of Q _p /R _s in future studies.     

A Model to Estimate Global Radiation in Complex Terrain

Global radiation is an important parameter necessary for most ecological models. However, in situ data barely meets the needs of modelling mountainous ecosystems since most field stations are located in flat areas. Consequently, it is usually necessary to extrapolate radiation measurements obtained from an adjacent flat area to the complex terrain of concern. The distribution of radiation in complex terrain depends upon two factors: the local atmospheric conditions, which determine the radiation potentially available to a supposed flat surface in a given location, and the topographic effects on this possible radiation. The latter have been included in detail in most radiation models for complex terrain, but the former are often only simply treated as constant or estimated by over-simplified empirical algorithms. In this paper we propose a novel model that uses a parametric atmospheric model to calculate the potential radiation for a supposed flat surface in a given location, and then account for topographic effects. Direct radiation, diffuse radiation and reflected radiation are calculated separately in the model due to the distinctive characteristics of and the effects by topography. Based on the parametric model, this paper has investigated the relationship between radiation transmittance, clearness indices and altitude under a series of water vapour content and turbidity conditions. This combines three ratios, R _b, R _d, and R _r, defined as the direct radiation, diffuse radiation and reflected radiation received by the arbitrary surface, respectively, to their counterparts in the horizontal surface, to estimate the global radiation for any given location. The model has been validated with data from measurements in National Park Berchtesgaden, Germany, where six measurement sites with various altitudes and topographic characteristics have been deployed. The r ² of modelled and measured hourly global radiation are greater than 0.90 in all six sites, with RMSE varies from 16 to 100 W m⁻². Sensitivity analysis revealed that the model was not sensitive to change in water vapour content, which suggests the possibility to use an exponential algorithm of water vapour content when there is no in situ water vapour content information in complex terrains. The NRMSE was only reduced by 0.04, on average, in five of the six sites when water vapour content information was calculated from the in situ air temperature and relative humidity measurements.     

支持向量机方法在南京太阳总辐射推算中的应用

利用1981、1996和2001年逐日南京站太阳总辐射和日照时数观测资料,建立了基于支持向量机(support vector machine,SVM)方法的太阳总辐射推算模型,预测了1982、1997和2002年的太阳总辐射,并把推算结果和采用线性的气候学方法所得到的推算结果分别与实测值进行对比。采用线性方法得到的1982、1997和2002年的太阳总辐射预测值与实测值间基于1:1线的决定系数(R~2)分别为0.800、0.859和0.838,均方根误差(RMSE)分别为3.250、2.649和2.925 MJ·m~(-2)·d~(-1)。采用SVM方法得到的1982、1997和2002年的R~2分别为0.894、0.938和0.936,RMSE分别为2.353、1.726和1.804 MJ·m~(-2)·d~(-1)。SVM方法得到的太阳总辐射预测值与实测值之间的误差较小,预测精度高于线性方法,更适用于实际太阳总辐射的计算。     

The stability of large—scale horizontal air motion in the non—linear basic zephyr flow under the effect of rossby parameter

The stability of large-scale horizontal motion in the atmosphere is discussed in this paper by using qualitative analysis theory of non-linear ordinary differential equations. Both the non-linear distribution of basic Zephyr flow and the variation of geostrophic vorticity along the latitude (f=f0 + βδy ) are all included in this paper's mathematical model so as to analogue the background field of large-scale horizontal air motion more really in the rotating reference frame of the earth. Some significant results are drawn out from this paper and the conclusions of Li(1986)'s and Wan et al.(1990)'s are extended widely.     

Seasonal variation and meteorological control of CO2 flux in a hilly plantation in the mountain areas of North China

The carbon cycle of terrestrial ecosystems is an important scientific issue in global climate change research.Plantation forest plays an important role in terrestrial carbon budget in China.In this study,eddy covariance flux data measured at Xiaolangdi forest ecosystem research station(XLD) in 2007 and 2008 are used to analyze the seasonal variation and meteorological control of CO2 flux in a 30-yr-old mixed plantation.The plantation forest mainly consists of Quercus variabilis,Platycladus orientalis,and Robinia pseudoacacia.The results show that the seasonal variations of net ecosystem exchange of CO2(NEE),gross primary production(GPP),and ecosystem respiration(Re) display single-peak curves.The maximum of carbon sequestration appears during May and June each year.The relative contribution of carbon release from ecosystem respiration to GPP varied slightly between 2007 and 2008.The relationship between NEE and photosynthetic active radiation(Qp) accords with the rectangular hyperbola model on diurnal scale,and shows a good linear correlation on monthly scale.The ecosystem photosynthetic parameters:the maximum photosynthetic rate(Pmax),the ecosystem photosynthetic photonyield(α),and the daytime ecosystem respiration(Rd) exhibit seasonal variations.Pmax reaches the maximum in August each year,with small interannual difference.The interannual differences of α and Rd are obvious,which is attributed to the changes of meteorological factors,such as solar radiation,vapor pressure deficit(D),precipitation,etc.Parameters Re,GPP,and NEP(net ecosystem production) have obvious exponential relations with temperature on monthly scale.There is a hysteresis in the response of GPP and NEP to temperature,i.e.,the carbon sequestration is not the maximum when the temperature reaches the peak value.The Q10 values were 1.37 and 1.45 in 2007 and 2008,respectively.On monthly scale,Re,GPP,and NEE increase as D increases,but rise slowly and even decrease when D is higher than 1.5 kPa.