Abstract

Observational data show that the remotely sensed leaf area index (LAI) has a significant downward trend over the east Tibetan Plateau (TP), while a warming trend is found in the same area. Further analysis indicates that this warming trend mainly results from the nighttime warming. The Single-Column Atmosphere Model (SCAM) version 3.1 developed by the National Center for Atmospheric Research is used to investigate the role of land use change in the TP local climate system and isolate the contribution of land use change to the warming. Two sets of SCAM simulations were performed at the Xinghai station that is located near the center of the TP Sanjiang (three rivers) Nature Reserve where the downward LAI trend is largest. These simulations were forced with the high and low LAIs. The modeling results indicate that, when the LAI changes from high to low, the daytime temperature has a slight decrease, while the nighttime temperature increases significantly, which is consistent with the observations. The modeling results further show that the lower surface roughness length plays a significant role in affecting the nighttime temperature increase.

1. Introduction

It has become apparent that climate change over high-elevation regions is occurring at a faster rate than over low-elevation regions [1]. This high-elevation rate of change often generates stronger disturbances within the climate system than low-elevation changes [2]. The Tibetan Plateau (TP) is an immense upland area ( ), with an average elevation greater than 5,000 m. Over the TP, 46% of the forest cover and 50% of the grassland have been converted to farmland, urban areas, or desert during the last century (X. Liu, personal communication). The TP is the source of snowmelt runoff, supplying water resources to users in China and surrounding areas. The mechanisms influencing the observed rapid climate change and their impacts at this high-elevation region are not well understood. In order to reduce these types of uncertainties and provide society with more reliable projections of future outcomes, detection and analysis of these processes need to be determined.

Many researchers have shown that land use and land cover change has remarkab interacts with climate systems [3, 4]. A statistical analysis of the southern half of the Central Valley was performed by Christy et al. [5] using high-quality temperature observations for the San Joaquin Valley in California, where land use has changed dramatically since the presettlement (circa 1860s), due to extensive agricultural expansion. Their study indicated that the daily maximum near-surface air temperature exhibited a cooling trend over the last century for this region. This finding is in opposition to the observed warming trends that prevail over most other regions (see [6, 7]; many others), suggesting that agricultural activity and associated irrigation processes most likely play a key role in producing this cooling trend. Modeling studies have shown that land use change affects regional climate by altering the components in the surface energy budget [8]. Based on the results from a general circulation model (GCM), Snyder et al. [9] indicate that removing all temperate forest and replacing it with bare soil produce cooling in the winter and spring over tropical and boreal areas due to an increase in the surface albedo, but a warming during the summer due to reduced evapotranspiration. If forced with a moderate amount of carbon dioxide, GCMs could produce similar results for the entire 21st century [10]. These GCM predictions also indicate that the response of land use change could be overridden by a strong atmospheric circulation system such as the Asian Monsoon circulation.

However, due to the lack of detailed long-term and densely distributed spatial observations, researchers have not been able to determine how and to what extent these land use changes affect the early and accelerated warming observed on the TP [11]. A high-quality climate model is an efficient tool to advance our understanding of this problem, and observed climate change could be quantitatively interpreted through the model details. The modeling results would further yield new insight into future field experimental design. Therefore, in this study, we use the Single-Column Atmosphere Model (SCAM) version 3.1 coupled with the Community Land Model version 3 (http://www.ccsm.ucar.edu/models/atm-cam/docs/scam/) to investigate the role of land use change in the TP local climate system and isolate the contribution of land use change to the warming of the TP from that of greenhouse gas emissions. The SCAM was developed by the National Center for Atmospheric Research (NCAR). In addition, the mechanisms associated with land use change signals are adequately quantified with the model output and our existing quality-controlled observations.

2. Model, Data, and Methodology

The SCAM used in this study is a one-dimensional time-dependent single-column atmospheric model. The SCAM is embedded in the NCAR Community Atmospheric Model version 3. Since the SCAM includes only one model grid cell, it can be used to effectively study land atmosphere interactions. This grid cell can include up to 10 land use types, giving a benefit to characterize land surface heterogeneity. In this study, the land use types are prescribed in the SCAM to represent land surface features in the TP. The SCAM is also configured with 26 vertical layers ranging from approximately 3 mb to the surface, and it includes all the physical schemes existing in its global counterpart, CAM3. The reanalyzed data by National Centers for Environmental Prediction (NCEP-R1) provide the initial and lateral boundary conditions for SCAM. The latter were updated every 6 hours. The initial conditions include soil moisture and temperature and atmospheric pressure, moisture, temperature, and winds. The 16 km resolution Advanced Very High Resolution Radiometer (AVHRR) leaf area indices (LAIs) were used to examine vegetation changes over the TP and also were inputted into SCAM to more realistically reflect vegetation variations. Station observations were used to analyze climate change over the TP and compare with modeling results.

3. Results

3.1. LAI, Temperature, and Precipitation Analysis

Land use and land cover on the TP have been changed significantly. Such change can be indentified with the AVHRR LAI trends over the period of 1983–2000 for July, August, and September, a season that the LAIs usually reach their maximum values (Figure 1). It is clearly seen that the LAIs have negative trends (blue colors) in the northeast TP, and the maximum trend can reach around −1.0/decade. These trends are calculated through a linear regression method. The region with the negative trends contains the headwaters of the Yellow River (Huang He), the Yangtze River (Chang Jiang), and the Mekong River (Lancang Jiang) that provide a significant amount of water resources for China. This area is also called Sanjiang (three rivers) Headwaters Nature Reserve. In some locations of the southeast TP, LAIs also show positive trends (red colors). Figure 2 shows the time series of the annual maximum LAI at a station called Xinghai that is located near the center of the Sanjiang area. This figure shows that the annual maximum LAI has a downward trend with a value of −0.67/decade over the period of 1983–2000. According to the United States Geological Survey data, Xinghai has a flat terrain and mostly is covered by grassland during the summer. All the following simulations were performed at this station to understand how vegetation affects local climate.

Figure 3(a) shows the time series of the annual precipitation anomalies at the Xinghai station for the period of 1960–2000. It is shown that the trend is 0.5 mm/decade that is quite insignificant. However, the annual mean temperature trend is C/decade (Figure 3(b)). The temperature increase is about C at this station over the study period. When temperatures at 2 AM and 2 PM are examined, it can be seen that the nighttime temperature trend is C/decade (Figure 3(c)) while the daytime temperature trend is only C/decade (Figure 3(d)). Thus, the annual mean temperature increase results from the nighttime warming. Identifying physical factors and processes that contribute to the nighttime warming is the major focus in the following sections, and the SCAM is used for such purpose.

3.2. SCAM Evaluation

Before physical processes that affect temperature change are studied, the SCAM needs to be fully evaluated with observations to ensure that it can reasonably describe climate processes in our study area. SCAM simulations at the Xinghai station were performed over the period of 1951–2000. The first 32-year (1951–1982) simulations were discarded for spinning up the model, and only the simulations over the period of 1983–2000 were analyzed. Figure 4(a) shows that although the SCAM produces a weaker seasonal cycle of the temperature it still generates reasonable temperature simulations at the Xinghai station when compared to the observations. The correlation coefficient between the observation and the simulation is 0.98, the RMSE is C, and the bias is C. The model can also reproduce the observed precipitation, where the correlation coefficient, RMSE, and bias are 0.74, 28 mm/month, and 13 mm/month, respectively (Figure 4(b)). These results show that the SCAM can simulate well for the seasonal cycles of precipitation and temperature at the Xinghai station, indicating that this model can be used to further explore physical processes that affect temperature and precipitation changes in this area.

3.3. Sensitivity Studies for Understanding the Temperature Changes

To examine how vegetation affects the warming in this area, two additional sets of simulations were performed with high and low LAIs for 1983–2000, a period where the remotely sensed LAIs are available. Figure 5 shows seasonal variations of high and low LAIs. The high LAIs are the average over 1983–1985, while the low LAIs are the average over 1998–2000. The largest difference between the high and low LAIs is more than one. The first set of simulations was forced with the high LAIs, where the LAI only varied seasonally but did not have interannual variations. The modeling settings for the second set of simulations were exactly the same as those for the first set of simulations but with the low LAIs. Figure 6(a) shows the 2 m height temperature difference between the high LAI and the low LAI cases averaged over the 18-year period (1983–2000). It is seen that the major changes in temperature (low LAI minus high LAI) occur in the warm season where the LAI difference between the two cases is the largest. During the cold season (from October through next May), there is almost no change for the temperature. Precipitation has very minor changes with a maximum difference of less than 1 mm/month (figure not shown). Figure 6(a) indicates that the lower LAI in this region generates a slight cooling during the summer at noontime with a maximum temperature change of C in July (solid line). All the components in the energy balance equation are examined. The results show that the sensible heat flux has most significant changes that affect 2 m height temperature (Figure 6(b)). It is seen that sensible heat flux decreases by 10 W/m2 at noontime during the summer (solid line). Further analysis indicates that such a decrease results from the lower roughness in the low LAI case. In SCAM, the roughness length for vegetation in this area is 0.06 m, and for bare soil, it is 0.024 m. The low LAI reduces the areal weight of vegetation in the model grid cell, but increases the weight of bare soil when compared to the settings in the high LAI case. Thus, on average, the low LAIs generate a flatter surface than the high LAIs, which decreases the upward sensible heat flux during daytime and slightly lower the near-surface temperature.

Figure 6(a) shows the nighttime temperature difference between the low and high LAI cases (low-high), indicating that the near-surface temperature increases by more than C in July due to changing the LAI from high to low (dashed line). During the evening, the sensible heat flux is downward (negative). The reduced roughness length in the low LAI case suppresses the energy exchanges between the surface and near-surface air. Figure 6(b) shows that the increase in sensible heat flux at midnight is up to 2 Wm2 (dashed line) that leads to more energy remaining in the near-surface air (less downward sensible heat) and increases the temperature there. Detailed calculation of sensible heat flux and how it is related to the roughness length are described by Oleson et al. in [12].

The above modeling results strongly suggest that the warming results from nighttime temperature changes that are consistent with the observations. The decreased vegetation LAI in this area significantly contributes to the nighttime temperature increase. The LAI decrease may not be a controlling factor that produces the observed warming trend, but most likely it reinforces the existing warming trend that is believed to be triggered by the carbon increase in the atmosphere [4].

4. Conclusions and Discussions

In northeast Tibetan Plateau, the remotely sensed LAIs show a significant downward trend, while a warming trend is observed for temperature. Our modeling results indicate that such a vegetation change is an important factor that contributes to the warming in this region. Further analysis indicates that this contribution largely comes from the evening temperature increase. The evening temperature increase results from the more stable near-surface boundary layer that restricts the downward sensible heat flux and keeps more energy in the near-surface air. The deterioration of vegetation leads to a flatter area, reduces the surface roughness, and results in more stable near-surface boundary layer.

This study is not intended to indicate that vegetation deterioration is a determining factor that leads to a warming in this area that is most likely related to global climate change triggered by the atmospheric carbon increase. Instead, the analysis shown in this study describes a fact that the reduced biomass in this area could strengthen the warming trend in this region. This study is based on the modeling results that need to be further verified with more observed evidence. The large- and mesoscale atmospheric information in the SCAM was introduced through the lateral boundary data. Thus, the atmosphere in the model column does not produce any feedback to that outside the column. This drawback could be overcome by applying a regional or global climate model to this region, which could be a future research topic. However, these modeling results give strong insight into future modeling and observational studies for climate change in this high-elevation area.

Acknowledgments

Support for this work was provided by the Utah Agricultural Experiment Station. J. Jin was also supported by K. C. Wong Fellowships, Chinese Academy of Sciences.