Long-term modeling of soil C erosion and sequestration at the small watershed scale

The soil C balance is determined by the difference between inputs (e.g., plant litter, organic amendments, depositional C) and outputs (e.g., soil respiration, dissolved organic C leaching, and eroded C). There is a need to improve our understanding of whether soil erosion is a sink or a source of atmospheric CO₂. The objective of this paper is to discover the long-term influence of soil erosion on the C cycle of managed watersheds near Coshocton, OH. We hypothesize that the amount of eroded C that is deposited in or out of a watershed compares in magnitude to the soil C changes induced via microbial respiration. We applied the erosion productivity impact calculator (EPIC) model to evaluate the role of erosion–deposition processes on the C balance of three small watersheds (∼1 ha). Experimental records from the USDA North Appalachian Experimental Watershed facility north of Coshocton, OH were used in the study. Soils are predominantly silt loam and have developed from loess-like deposits over residual bedrock. Management practices in the three watersheds have changed over time. Currently, watershed 118 (W118) is under a corn (Zea mays L.)–soybean (Glycine max [L.] Merr.) no till rotation, W128 is under conventional till continuous corn, and W188 is under no till continuous corn. Simulations of a comprehensive set of ecosystem processes including plant growth, runoff, and water erosion were used to quantify sediment C yields. A simulated sediment C yield of 43 ± 22 kg C ha⁻¹ year⁻¹ compared favorably against the observed 31 ± 12 kg C ha⁻¹ year⁻¹ in W118. EPIC overestimated the soil C stock in the top 30-cm soil depth in W118 by 21% of the measured value (36.8 Mg C ha⁻¹). Simulations of soil C stocks in the other two watersheds (42.3 Mg C ha⁻¹ in W128 and 50.4 Mg C ha⁻¹ in W188) were off by <1 Mg C ha⁻¹. Simulated eroded C re-deposited inside (30–212 kg C ha⁻¹ year⁻¹) or outside (73–179 kg C ha⁻¹ year⁻¹) watershed boundaries compared in magnitude to a simulated soil C sequestration rate of 225 kg C ha⁻¹ year⁻¹ and to literature values. An analysis of net ecosystem carbon balance revealed that the watershed currently under a plow till system (W128) was a source of C to the atmosphere while the watersheds currently under a no till system (W118 and W188) behaved as C sinks of atmospheric CO₂. Our results demonstrate a clear need for documenting and modeling the proportion of eroded soil C that is transported outside watershed boundaries and the proportion that evolves as CO₂ to the atmosphere.     

Net CO2 exchange of a subarctic mountain birch ecosystem

Summary  Net ecosystem CO₂ exchange was measured over a mountain birch forest in northern Finland throughout the growing season. The maximal net CO₂ uptake rate of about − 0.5 mg(CO₂) m⁻² s⁻¹ was observed at the end of July. The highest nocturnal respiration rates in early August were 0.2 mg(CO₂) m⁻² s⁻¹. The daily CO₂ balances during the time of maximal photosynthesis were about −15 g(CO₂) m⁻² d⁻¹. The mountain birch forest acted as a net sink of CO₂ from 30 June to 28 August. During that period the net CO₂ balance was −448 g(CO₂)m⁻². The interannual representativeness of the observed balances was studied using a simplified daily balance model, with daily mean global radiation and air temperature as the input parameters. The year-to-year variation in the phenological development was parameterised as a function of the cumulative effective temperature sum. The daily balance model was used for estimating the variability in the seasonal CO₂ balances due to the timing of spring and meteorological factors. The sink term of CO₂ in 1996 was lower than the 15-year mean, mainly due to the relatively late emergence of the leaves. Received October 11, 1999 Revised April 25, 2000     

Simulation and Evaluation of Terrestrial Ecosystem NPP with M-SDGVM over Continental China

Using the regional terrestrial Net Primary Production (NPP) from different observations and models over China, we validated the NPP simulations and explored the relationship between NPP and climate variation at interannual and decadal scales in the Modified Sheffield Dynamic Global Vegetation Model (M-SDGVM) during 1981–2000. M-SDGVM shows agreement with the NPP data from 743 sites under the Global Primary Production Data Initiative (GPPDI). The spatial and the zonal averaged NPP of M-SDGVM agree well with ...     

Turbulent transport of carbon dioxide over a paddy field

Turbulent fluctuations in CO₂ concentrations over a paddy field are measured by a fastresponse device with an open sensing path. This IR device coupled with a sonic anemometer constitutes an eddy correlation instrument to measure CO₂ fluxes. Three experiments were conducted in the surface layer over paddy 90 cm high. The stability (z – d)/L ranged from -0.14 to 0.20, where L denotes the Monin-Obukhov length.CO₂ power spectra show the range of applicability of the -2/3 power law to be between f = 0.2 and f = 2, where f is the frequency normalized by wind speed and height. The cospectral estimate between CO₂ and vertical component of wind speed ranging from f = 0.005 to f = 2 shows a peak at about f = 0.15 under near-neutral stratification.Hourly means of CO₂ flux measured by the eddy correlation method increase with intensity of net radiation. The maximum value of downward flux of CO₂ rises to 0.6 mg cm^-2 hr^-1 over the paddy field at the stage of ear emergence.Some turbulence statistics relating to the CO₂ transport are evaluated: the correlation coefficient between CO₂ and vertical velocity is about -0.3, and that between CO₂ and humidity attains -0.7 ~ -0.8 under unstable stratification; nondimensional gradients _c for CO₂ and _m for wind speed are 0.89 and 0.99, respectively.     

Continuous measurement of methane flux over a larch forest using a relaxed eddy accumulation method

We measured the methane flux of a forest canopy throughout a year using a relaxed eddy accumulation (REA) method. This sampling system was carefully validated against heat and CO₂ fluxes measured by the eddy covariance method. Although the sampling system was robust, there were large uncertainties in the measured methane fluxes because of the limited precision of the methane gas analyzer. Based on the spectral characteristics of signals from the methane analyzer and the diurnal variations in the standard deviation of the vertical wind velocity, we found the daytime and nighttime precision of half-hourly methane flux measurements to be approximately 1.2 and 0.7?μg?CH₄?m^?2?s^?1, respectively. Additional uncertainties caused by the dilution effect were estimated to affect the accuracy by as much as 0.21?μg?CH₄?m^?2?s^?1 on a half-hourly basis. Diurnal and seasonal variations were observed in the measured fluxes. The biological emission from plant leaves was not observed in our studies, and thus could be negligible at the canopy-scale exchange. The annual methane sink was 835?±?175?mg?CH₄?m^?2?year^?1 (8.35?kg?CH₄?ha^?1?year^?1), which was comparable to the flux range of 379–2,478?mg?CH₄?m^?2?year^?1 previously measured in other Japanese forest soils. This study indicated that the REA method could be a promising technique to measure canopy scale methane fluxes over forests, but further improvement of precision of the analyzer will be required.     

Quantification of ecosystem carbon exchange characteristics in a dominant subtropical evergreen forest ecosystem

CO₂ fluxes were measured continuously for three years (2003?C2005) using the eddy covariance technique for the canopy layer with a height of 27 m above the ground in a dominant subtropical evergreen forest in Dinghushan, South China. By applying gapfilling methods, we quantified the different components of the carbon fluxes (net ecosystem exchange (NEE)), gross primary production (GPP) and ecosystem respiration (R_eco) in order to assess the effects of meteorological variables on these fluxes and the atmospherecanopy interactions on the forest carbon cycle. Our results showed that monthly average daily maximum net CO₂ exchange of the whole ecosystem varied from ?3.79 to ?14.24 ??mol m^?2 s^?1 and was linearly related to photosynthetic active radiation. The Dinghushan forest acted as a net carbon sink of ?488 g C m^?2 y^?1, with a GPP of 1448 g Cm^?2 y^?1, and a R_eco of 961 g C m^?2 y^?1. Using a carboxylase-based model, we compared the predicted fluxes of CO₂ with measurements. GPP was modelled as 1443 g C m^?2 y^?1, and the model inversion results helped to explain ca. 90% of temporal variability of the measured ecosystem fluxes. Contribution of CO₂ fluxes in the subtropical forest in the dry season (October-March) was 62.2% of the annual total from the whole forest ecosystem. On average, 43.3% of the net annual carbon sink occurred between October and December, indicating that this time period is an important stage for uptake of CO₂ by the forest ecosystem from the atmosphere. Carbon uptake in the evergreen forest ecosystem is an indicator of the interaction of between the atmosphere and the canopy, especially in terms of driving climate factors such as temperature and rainfall events. We found that the Dinghushan evergreen forest is acting as a carbon sink almost year-round. The study can improve the evaluation of the net carbon uptake of tropical monsoon evergreen forest ecosystem in south China region under climate change conditions.     

Methane emissions from a landfill measured by eddy correlation using a fast response diode laser sensor

We describe a fast response methane sensor based on the absorption of radiation generated with a near-infrared InGaAsP diode laser. The sensor uses an open path absorption region 0.5 m long; multiple pass optics provide an optical path of 50 m. High frequency wavelength modulation methods give stable signals with detection sensitivity (S/N=1, 1 Hz bandwidth) for methane of 65 ppb at atmospheric pressure and room temperature. Improvements in the optical stability are expected to lower the current detection limit. We used the new sensor to measure, by eddy correlation, the CH₄ flux from a clay-capped sanitary landfill. Simultaneously we measured the flux of CO₂ and H₂O. From seven half-hourly periods of data collected after a rainstorm on November 23, 1991, the average flux of CH₄ was 17 mmol m^–2 hr^–1 (6400 mg CH₄ m^–2 d^–1) with a coefficient of variation of 25%. This measurement may underrepresent the flux by 15% due to roll-off of the sensor response at high frequency. The landfill was also a source of CO₂ with an average flux of 8.1 mmol m^–2 hr^–1 (8550 mg CO₂ m^–2 d^–1) and a coefficient of variation of 26%. A spectral analysis of the data collected from the CH₄, CO₂, and H₂O sensors showed a strong similarity in the turbulent transfer mechanisms.     

Carbon sequestration and environmental effects of afforestation with Pinus radiata D. Don in the Western Cape, South Africa

A three-step methodology to assess the carbon sequestration and the environmental impact of afforestation projects in the framework of the Flexible Mechanisms of the Kyoto Protocol (Joint Implementation and Clean Development Mechanism) was developed and tested using a dataset collected from the Jonkershoek forest plantation, Western Cape, South Africa, which was established with Pinus radiata in former native fynbos vegetation and indigenous forest. The impact of a change in land use was evaluated for a multifunctional, a production and a non-conversion scenario. First, the carbon balance was modelled with GORCAM and was expressed as (1) C sequestration in tC ha⁻¹ year⁻¹ in soil, litter, and living biomass according to the rules of the first commitment period of the Kyoto Protocol, and (2) CO₂ emission reductions in tC ha⁻¹ year⁻¹, which includes carbon sequestered in the above-mentioned pools and additionally in wood products, as well as emission reductions due to fossil fuel substitution. To estimate forest growth, three data sources were used: (1) inventory data, (2) growth simulation with a process-based model, and (3) yield tables. Second, the effects of land use change were assessed for different project scenarios using a method related to Life Cycle Assessment (LCA). The method uses 17 quantitative indicators to describe the impact of project activities on water, soil, vegetation cover and biodiversity. Indicator scores were calculated by comparing indicator values with reference values, estimated for the climax vegetation. The climax vegetation is the site-specific ecosystem phase with the highest exergy content and the highest exergy flow dissipation capacity. Third, the land use impact per functional unit of 1 tC sequestered was calculated by combining the results of step 1 and step 2. The average baselines to obtain carbon additionality are 476 tC ha⁻¹ for indigenous forest and 32 tC ha⁻¹ for fynbos. Results show that the influence of the growth assessment method on the magnitude of C sequestration and hence on the environmental impact per functional unit is large. When growth rate is assessed with the mechanistic model and with the yield table, it is overestimated in the early years and underestimated in the long term. The main conclusion of the scenario analysis is that the production forest scenario causes higher impacts per functional unit than the multifunctional scenario, but with the latter being less efficient in avoiding CO₂ emissions. The proposed method to assess impacts on diverse components of the ecosystem is able to estimate the general tendency of the adverse and positive effects of each scenario. However, some indicators, more specifically about biodiversity and water balance, could be improved or reinterpreted in light of specific local data about threat to biodiversity and water status.     

Carbon cycle–climate feedback sensitivity to parameter changes of a zero-dimensional terrestrial carbon cycle scheme in a climate model of intermediate complexity

Summary A series of sensitivity runs have been performed with a coupled climate–carbon cycle model. The climatic component consists of the climate model of intermediate complexity IAP RAS CM. The carbon cycle component is formulated as a simple zero-dimensional model. Its terrestrial part includes gross photosynthesis, and plant and soil respirations, depending on temperature via Q ₁₀-relationships (Lenton, 2000). Oceanic uptake of anthropogenic carbon is formulated is a bi-linear function of tendencies of atmospheric concentration of CO₂ and globally averaged annual mean sea surface temperature. The model is forced by the historical industrial and land use emissions of carbon dioxide for the second half of the 19th and the whole of the 20th centuries, and by the emission scenario SRES A2 for the 21st century. For the standard set of the governing parameters, the model realistically captures the main features of the Earth’s observed carbon cycle. A large number of simulations have been performed, perturbing the governing parameters of the terrestrial carbon cycle model. In addition, the climate part is perturbed, either by zeroing or artificially increasing the climate model sensitivity to the doubling of the atmospheric CO₂ concentration. Performing the above mentioned perturbations, it is possible to mimic most of the range found in the C4MIP simulations. In this way, a wide range of the climate–carbon cycle feedback strengths is obtained, differing even in the sign of the feedback. If the performed simulations are subjected to the constraints of a maximum allowed deviation of the simulated atmospheric CO₂ concentration (pCO_2(a)) from the observed values and correspondence between simulated and observed terrestrial uptakes, it is possible to narrow the corresponding uncertainty range. Among these constraints, considering pCO_2(a) and uptakes are both important. However, the terrestrial uptakes constrain the simulations more effectively than the oceanic ones. These constraints, while useful, are still unable to rule out both extremely strong positive and modest negative climate–carbon cycle feedback.     

Estimating climate change effects on net primary production of rangelands in the United States

The potential effects of climate change on net primary productivity (NPP) of U.S. rangelands were evaluated using estimated climate regimes from the A1B, A2 and B2 global change scenarios imposed on the biogeochemical cycling model, Biome-BGC from 2001 to 2100. Temperature, precipitation, vapor pressure deficit, day length, solar radiation, CO₂ enrichment and nitrogen deposition were evaluated as drivers of NPP. Across all three scenarios, rangeland NPP increased by 0.26 % year^?1 (7 kg C ha^?1 year^?1) but increases were not apparent until after 2030 and significant regional variation in NPP was revealed. The Desert Southwest and Southwest assessment regions exhibited declines in NPP of about 7 % by 2100, while the Northern and Southern Great Plains, Interior West and Eastern Prairies all experienced increases over 25 %. Grasslands dominated by warm season (C4 photosynthetic pathway) species showed the greatest response to temperature while cool season (C3 photosynthetic pathway) dominated regions responded most strongly to CO₂ enrichment. Modeled NPP responses compared favorably with experimental results from CO₂ manipulation experiments and to NPP estimates from the Moderate Resolution Imaging Spectroradiometer (MODIS). Collectively, these results indicate significant and asymmetric changes in NPP for U.S. rangelands may be expected.     

Simulated Climate Change Effects on Year-Round Water Temperatures in Temperate Zone Lakes

A deterministic, one-dimensional model is presented to simulate daily water temperature profiles and associated ice and snow covers for dimictic and polymictic lakes of the temperate zone. The lake parameters required as model input are surface area (A_s), maximum depth (H_MAX), and Secchi depth (z_s), the latter, used as a measure of light attenuation and trophic state. The model is driven by daily weather data and operates year-round over multiple years. The model has been tested with extensive data (over 5,000 temperature points). Standard error between simulated and measured water temperatures is 1.4°C in the open water season and 0.5°C in the ice cover season. The model is applied to simulate the sensitivity of Minnesota lake water temperature characteristics to climate change. The projected climate changes due to a doubling of atmospheric CO₂ are obtained from the output of the Canadian Climate Center General Circulation Model (CCC GCM) and the Goddard Institute of Space Studies General Circulation Model (GISS GCM). Simulated lake temperature characteristics have been plotted in a coordinate system with a lake geometry ratio (A _s ^0.25 /H_MAX) on one axis and Secchi depth on the other. The lake geometry ratio expresses a lake's susceptibility to stratification. By interpolation, the sensitivity of lake temperature characteristics to changes of water depth and Secchi depth under the projected climate scenarios can therefore be obtained. Selected lake temperature characteristics simulated with past climate conditions (1961–1979) and with a projected 2 × CO₂ climate scenario as input are presented herein in graphical form. The simulation results show that under the 2 × CO₂ climate scenario ice formation is delayed and ice cover period is shortened. These changes cause water temperature modifications throughout the year.     

Effects of temperature rise and increase in CO2 concentration on simulated wheat yields in Europe

A crop-growth-simulation model based on SUCROS87 was used to study effects of temperature rise and increase of atmospheric CO₂ concentration on wheat yields in several regions in Europe. The model simulated potential and water-limited crop production (growth with ample supply of nutrients and in the absence of damage by pests, diseases and weeds). Historic daily weather data from 13 sites in Western Europe were used as starting point.For potential production (optimal water) a 3 °C temperature rise led to a yield decline due to a shortening of the growing period on all locations. Doubling of the CO₂ concentration caused an increase in yield of 40% due to higher assimilation rates. It was found that effects of higher temperature and higher CO₂ concentration were nearly additive and the combination of both led to a yield increase of 1–2 ton ha^-1. A very small CO₂-temperature interaction was found: the effect of doubled CO₂ concentration on crop yield was larger at higher temperatures. The inter-annual yield variability was hardly affected.When water was limiting crop-production effects of temperature rise and higher CO₂ levels were different than for the potential production. Rise in temperature led to a smaller yield reduction, doubled CO₂ concentration to a larger yield increase and combination of both led to a large yield increase (3 ton ha^-1) in comparison with yields simulated for the present situation. Both rise in temperature and increase in the CO₂ concentration reduced water requirements of the crop. Water shortages became smaller, leading to a reduction in inter-annual variability. It is concluded that when no major changes in precipitation pattern occur a climate change will not affect wheat yields since negative effects of higher temperatures are compensated by positive effects of CO₂ enrichment.     

How is CO2 affecting yields and technological progress? A statistical analysis

This paper analyzes the impact of climate, crop production technology, and atmospheric carbon dioxide (CO₂) on current and future crop yields. The analysis of crop yields endeavors to advance the literature by estimating the effect of atmospheric CO₂ on observed crop yields. This is done using an econometric model estimated over pooled historical data for 1950–2009 and data from the free air CO₂ enrichment experiments. The main econometric findings are: 1) Yields of C3 crops (soybeans, cotton, and wheat) directly respond to the elevated CO₂, while yields of C4 crops (corn and sorghum) do not, but they are found to indirectly benefit from elevated CO₂ in times and places of drought stress; 2) The effect of technological progress on mean yields is non-linear; 3) Ignoring atmospheric CO₂ in an econometric model of crop yield likely leads to overestimates of the pure effects of technological progress on crop yields of about 51, 15, 17, 9, and 1 % of observed yield gain for cotton, soybeans, wheat, corn and sorghum, respectively; 4) Average climate conditions and climate variability contribute in a statistically significant way to average crop yields and their variability; and 5) The effect of CO₂ fertilization generally outweighs the effect of climate change on mean crop yields in many regions resulting in an increase of 7–22, 4–47, 5–26, 65–96, and 3–35 % for yields of corn, sorghum, soybeans, cotton, and wheat, respectively.     

全球农作物对大气CO2及其倍增的吸收量估算

根据农作物产量资料（ＦＡＯ１９９２年），计算出中国和全球各种作物对ＣＯ２的吸收总量分别为５．５×１０８ｔ／ａＣ和２８．９×１０８ｔ／ａＣ。同时以不同ＣＯ２浓度下小麦、玉米、大豆等全生育期光合速率实验数据直接计算的Ｃ吸收量为对照，与相应的中国产量资料计算结果比较，两者相差２．６％。从而进一步依据作物对ＣＯ２倍增反应诊断实验结果，推算出大气ＣＯ２浓度比目前倍增（７００ｐｐｍ）条件下，中国和全球农作物吸收ＣＯ２总量将增长２１％－２６％，分别为６．６×１０８ｔ／ａ—６．９×１０８ｔ／ａ和３４．１×１０８ｔ／ａ—３６．２×１０８ｔ／ａＣ。研究还表明，单位面积作物年吸Ｃ量全球（３．２ｔ／（ｈｍ２·８））比中国（４．２ｔ／（ｈｍ２·ａ））低２５．４％，而且Ｃ４作物普遍高于同类Ｃ３作物。     

Sensitivity of carbon budget to historical climate variability and atmospheric CO2 concentration in temperate grassland ecosystems in China

Chinese temperate grasslands play an important role in the terrestrial carbon cycle. Based on the parameterization and validation of Terrestrial Ecosystem Model (TEM, Version 5.0), we analyzed the carbon budgets of Chinese temperate grasslands and their responses to historical atmospheric CO₂ concentration and climate variability during 1951–2007. The results indicated that Chinese temperate grassland acted as a slight carbon sink with annual mean value of 7.3 T?g C, ranging from -80.5 to 79.6 T?g C yr^-1. Our sensitivity experiments further revealed that precipitation variability was the primary factor for decreasing carbon storage. CO₂ fertilization may increase the carbon storage (1.4 %) but cannot offset the proportion caused by climate variability (-15.3 %). Impacts of CO₂ concentration, temperature and precipitation variability on Chinese temperate grassland cannot be simply explained by the sum of the individual effects. Interactions among them increased total carbon storage of 56.6 T?g C which 14.2 T?g C was stored in vegetation and 42.4 T?g C was stored in soil. Besides, different grassland types had different responses to climate change and CO₂ concentration. NPP and R_H of the desert and forest steppes were more sensitive to precipitation variability than temperature variability while the typical steppe responded to temperature variability more sensitively than the desert and forest steppes.     

Simulating the impacts of climate change on cotton production in India

General circulation models (GCMs) project increases in the earth’s surface air temperatures and other climate changes by the mid or late 21st century, and therefore crops such as cotton (Gossypium spp L.) will be grown in a much different environment than today. To understand the implications of climate change on cotton production in India, cotton production to the different scenarios (A2, B2 and A1B) of future climate was simulated using the simulation model Infocrop-cotton. The GCM projections showed a nearly 3.95, 3.20 and 1.85 °C rise in mean temperature of cotton growing regions of India for the A2, B2 and A1B scenarios, respectively. Simulation results using the Infocrop-cotton model indicated that seed cotton yield declined by 477 kg?ha^?1 for the A2 scenario and by 268 kg?ha^?1 for the B2 scenario; while it was non-significant for the A1B scenario. However, it became non-significant under elevated [CO₂] levels across all the scenarios. The yield decline was higher in the northern zone over the southern zone. The impact of climate change on rainfed cotton which covers more than 60 % of the country’s total cotton production area (mostly in the central zone) and is dependent on the monsoons is likely to be minimum, possibly on account of marginal increase in rainfall levels. Results of this assessment suggest that productivity in northern India may marginally decline; while in central and southern India, productivity may either remain the same or increase. At the national level, therefore, cotton production is unlikely to change with climate change. Adaptive measures such as changes in planting time and more responsive cultivars may further boost cotton production in India.     

Climate change,weather variability and corn yield at a higher latitude locale: Southwestern Quebec

Climate change has led to increased temperatures, and simulation models suggest that this should affect crop production in important agricultural regions of the world. Nations at higher latitudes, such as Canada, will be most affected. We studied the relationship between climate variability (temperature and precipitation) and corn yield trends over a period of 33 years for the Monteregie region of south-western Quebec using historical yield and climate records and statistical models. Growing season mean temperature has increased in Monterregie, mainly due to increased September temperature. Precipitation did not show any clear trend over the 33 year period. Yield increased about 118 kg ha⁻¹ year⁻¹ from 1973 to 2005 (under normal weather conditions) due mainly to changes in technology (genetics and management). Two climate variables were strongly associated with corn yield variability: July temperature and May precipitation. These two variables explain more than a half of yield variability associated with climate. In conclusion, July temperatures below normal and May precipitation above normal have negative effects on corn yield, and the growing seasons have warmed, largely due to increases in the September temperature.     

Overestimated Biomass Carbon Pools of the Northern mid- and High Latitude Forests

The biomass carbon (C) stock of forests is one of key parameters for the study of regional and global carbon cycles. Literature reviews shows that inventory-based forest C stocks documented for major countries in the middle and high northern latitudes fall within a narrow range of 36–56 Mg C ha⁻¹ with an overall area-weighted mean of 43.6 Mg C ha⁻¹. These estimates are 0.40 to 0.71 times smaller than those (61–108 Mg C ha⁻¹) used in previous analysis of balancing the global carbon budget. A statistical analysis, using the global forest biomass database, implies that aboveground biomass per hectare is proportional to forest mean height [biomass in Mg/ha = 10.63 (height in m)] in closed-canopy forests in the study regions, indicating that forest height can be a proxy of regional biomass C stocks. The narrow range of C stocks is likely a result of similar forest height across the northern regions. The lower biomass C stock obtained in this study strongly suggests that the role of the northern forests in the global carbon cycle needs to be re-evaluated. Our findings also suggest that regional estimates of biomass could be readily made from the use of satellite methods such as lidar that can measure forest canopy height over large regions.     

Effects of an Urban Park and Residential Area on the Atmospheric CO2 Concentration and Flux in Seoul, Korea

The CO₂ concentrations and fluxes over an urban forest site (Namsan) and an urban residential region (Boramae) in Seoul, Korea, during the non-growing season (2–4 March 2011), the growing season (10–12 June 2011), and the late-growing season (22–24 September 2011) were analyzed. The CO₂ concentrations of two sites showed nearly the same diurnal variation, with a maximum value occurring during the night and a minimum value occurring during daytime, as well as the same seasonal variation, with a maximum value during the non-growing season (early spring) and a minimum value during the growing season (summer). The CO₂ flux over the urban forest did not show any typical diurnal variation during the non-growing season, but did show diurnal variation with a small positive value during the night and a large negative value during daytime in the growing and late-growing seasons due to photosynthesis in the urban forest. The CO₂ flux over the urban residential region showed a positive daily mean value for all periods, with large values during the non-growing season and small values during the growing season, and it also showed diurnal variation with two maxima at 0600–1000 LST and 1800–2400 LST, and two minima at 0300-0600 LST and 1100-1500 LST, and was strongly correlated with the use of liquefied natural gas for cooking and heating by surrounding houses.     

Characterizing the annual-mean climatic effect of anthropogenic CO<Subscript>2</Subscript> and aerosol emissions in eight coupled atmosphere-ocean GCMs

This work examines the spatial patterns of the transient response of mean annual temperature and precipitation to CO₂ (or CO₂ plus aerosol or aerosol proxy) radiative forcing in eight coupled AOGCMs, generally for the period 1900–2099. Response patterns are characterized using empirical orthogonal functions (EOFs) and the quasi-EOFs of Harvey and Wigley (the first qEOF field, discussed here, is given by the correlation between local year-by-year temperature changes and the global mean temperature change). The first temperature EOF accounts for 80–95% of the space-time variation of the CO₂ run in all of the models, and is almost identical to qEOF₁ of the temperature response or to the temperature change pattern averaged over the last 30 years of the simulations. EOF₁ accounts for 80–95% of the space-time variation in the CO₂+aerosol runs in six of the eight models. The CO₂ response patterns of different models are highly correlated with one another (R ² generally >0.5), and are also highly correlated with the CO₂+aerosol response patterns (R ² 0.85 in all except one model). The difference between CO₂ and CO₂+aerosol runs can be represented by EOF₁ of the year-by-year differences, by qEOF₁ of the year-by-year differences, or by the difference in temperature averaged over the last 30 years of each run. In models where these representations are highly correlated with each other, they are also highly correlated with CO₂ EOF₁. In other cases, aerosol EOF₁ is modestly to highly correlated with control EOF₁ (i.e.: the year-by-year differences between CO₂ and CO₂+aerosol runs are dominated by internal variability), while aerosol qEOF₁ and the 30-year difference are highly correlated with each other. For all models, the decadal mean temperature change can be closely replicated by scaling the CO₂ EOF₁ pattern based on the global mean temperature changes (RMSE for the last decade is <6% of the RMS temperature change for CO₂ runs, <8% for CO₂+aerosol runs). The first EOF of the precipitation response to increasing CO₂ accounts for only 10–30% of the space-time variation, and is generally highly correlated (R ² up to 0.85) with control EOF₁. In all of the models, there is an increase in precipitation in the ITCZ and a decrease in bands at or near 30°S and 30°N. In many models there is an El Niño-like response, including a substantial decrease in precipitation over the Amazon. Global-mean precipitation increases in all models due to CO₂ forcing, but aerosols appear to have a disproportionally large effect in suppressing the increase compared to their effect in suppressing the warming. There is evidence in some models that the non-absorbing aerosols considered here reduce summer monsoon rainfall compared to the changes that would be expected based on the globally averaged effect of aerosols on precipitation. When regional precipitation changes over time are predicted by scaling a fixed precipitation-change pattern with the global mean temperature change, the global mean RMSE in the predicted change in decadal-mean precipitation is 25–35% of the global RMS precipitation changes by the end of the simulation.