Introduction

Functional deterioration of ecosystems due to human activities has raised concerns of an accelerated aggravation of the environment along with climate change (Millennium Ecosystem Assessment, 2005). From the perspective of ecosystem services, the functions of an ecosystem can be categorized into four categories, namely, regulating, supporting, provisioning, and cultural, which can be quantified through economical valuations (de Groot et al., 2002; Liu et al., 2010; Petter et al., 2012). These functions typically reflect the interrelationships and interactions among various ecosystem components such as vegetation, water, soil, atmosphere, and even humans. Especially, the climate regulation function of ecosystems is an important factor for assessing ecosystem services and understanding the feedback from ecosystems to the atmosphere, which in-turn provides insight regarding climate changes (Scholes, 2016).

The climate regulation function of terrestrial ecosystems is mainly characterized by two factors: biogeochemical and biogeophysical processes. Biogeochemical processes describe pathways of carbon exchanges between terrestrial ecosystems and the atmosphere, including carboxylation through photosynthetic activities of vegetated surfaces and carbon emissions via the respiration of vegetation and soil micro-organisms. It is known that an increase in carbon concentration in the atmosphere and the resulting global warming intensify photosynthetic activities, which increases the carbon sink at the surface, consequently resulting in an offsetting effect. Meanwhile, biogeophysical processes are related to thermal energy exchanges between ecosystems and the atmosphere. Incoming solar radiation is transformed and reemitted into the atmosphere as radiation forcing such as sensible heat flux (SH), latent heat flux (LH), and ground heat flux (GH). These fluxes characterize the thermal equilibrium between the land and atmosphere and are determined using the thermal characteristics of land covers and land uses. Of the surface fluxes, SH is the main factor that affects atmospheric changes (Anderson-Teixeira et al., 2012; Feddema et al., 2005; Snyder et al., 2004; West et al., 2011) since it constitutes the majority of the land surface energy budget. LH is one of the major factors involved in the hydrological circulation of terrestrial ecosystems because it is emitted in the form of water vapor. This indicates that terrestrial ecosystems are one of the determinants for partitioning surface water balance which is, consequently, important to quantify the water regulating function of ecosystems. Moreover, LH affects the global water circulation in the entire earth system and the regional variabilities in water regimes. Kim et al. (2016) showed that over 60% of the total precipitation (PRCP) in East Asia consists of water vapors returned to the atmosphere through evapotranspiring processes over the region. This further indicates that LH should be treated as an important negative feedback from the terrestrial ecosystem to climate change, considering the cooling effects from clouding and rainfall events.

Thus, to accurately diagnose and predict the impact of climate change, it is important to understand how changes in the climate regulation function of terrestrial ecosystems affect climate change. However, traditional ecological approaches to quantify the ecosystem function tend to focus only on the carbon storing capacities of different types of vegetation (Castro, 2014; Guimarães et al., 2017; Ma et al., 2019; Petter et al., 2012; Wood et al., 2018; Yu et al., 2002). However, these approaches often underestimate the role of ecosystems in the energy transfer processes at the surface (Anderson-Teixeira et al., 2012; 2013). On the contrary, recent developments in land surface models (LSMs) have enabled researchers to describe phenological variations in vegetated surfaces (Myoung et al., 2011). However, feedback from land surfaces to the atmosphere in the energy, water, and carbon cycles has not been characterized well yet as the quantitative effects of this feedback remain unclear (Hong et al., 2007; Lakshmi et al., 2011; Myoung et al., 2013). It is very challenging to represent this feedback across various spatial scales, such as from the microscopic scale to the continental scale. Nevertheless, meaningful implications regarding the contributions of ecosystems to climate change can be extracted if given with sufficient long-term analyses. This is because the atmospheric changes can be assumed to reflect the feedback from terrestrial ecosystem.

This study aims to investigate the spatio-temporal contributions of terrestrial ecosystem to global warming, quantifying the biogeochemical and biogeophysical factors related to the climate regulation function of terrestrial ecosystems. A sufficiently large study area (East Asia) and long time period (more than half century) were set to extract reasonable trend in climate regulating components of terrestrial ecosystem from a LSM in which a phenology model is implemented.

Materials|Methods

Land surface model with vegetation dynamics

The interactions between land and atmosphere are essential features to describe the function of terrestrial ecosystems in climate regulation. These interactions are typically implemented as a module in climate and weather forecasting models at the global or regional scale. The first incorporation of land surface processes in a global or regional atmospheric model was a simple bucket model that considered only a one-layer surface (Manabe, 1969). Then, improved LSMs were developed, which considered more complex processes in the energy and water cycles (Deardorff, 1978; Dickinson et al., 1986; Eagleson, 1982; Pan & Mahrt, 1987). Since terrestrial ecosystems have been known to play very important roles in the energy and water cycles, climate models have increasingly introduced vegetation schemes to reflect the physiological and phenological processes of vegetated surfaces. To quantify the emissions of greenhouse gases, it is important to represent the biological mechanism for the carbon cycle within vegetated surfaces via the photosynthetic activities and respiration of vegetation. In the past 20 years, dynamic vegetation models have been developed to simulate the temporal variations in biogeophysical and biogeochemical processes with those in atmospheric features such as temperature, PRCP, wind, and etc. In addition to considering the interactions between land and atmosphere, these models also embody functional relationships between the climate and ecological factors (Daly et al., 2000; Prentice et al., 2000; Running & Couglan, 1988).

The Noah LSM with multi-parameterization options (Noah-MP) was introduced to represent the changes in climate regulation function of terrestrial ecosystems. Noah-MP is an advanced version of the Noah LSM (Niu et al., 2011). The Noah LSM describes land and atmospheric interactions based on thermal energy and water balance and was first developed by Oregon State University (Pan & Mahrt, 1987). This model was extended by incorporating a vegetation transpiration scheme and a surface runoff scheme (Chen et al., 1996). Since then, this model has been widely used in a standalone off-line mode, as well as in a coupled mode with a regional weather forecasting model (Chen & Dudhia, 2001). On the basis of the Noah LSM 3.0, Noah-MP provides additional scheme options: a dynamic vegetation model (Dickinson et al., 1998), a simple groundwater model with a TOPMODEL runoff scheme (Niu et al., 2005), a radiation transfer scheme modified for the three-dimensional structure effect of vegetation canopy (Niu & Yang, 2004), a 3-layer snow model (Yang & Niu, 2003), and a frozen soil scheme for separated permeable and impermeable fractions (Niu & Yang, 2006). Thus, Noah-MP allows a huge number of ensemble simulations to be conducted via scheme combinations. Through various experiments using Noah-MP, Yang et al. (2011) proposed a combination of physical schemes for more accurate simulations at the global and continental scales. Hong et al. (2014) investigated the applicability of genetic algorithms to determine an optimal set of physical schemes for Noah-MP for South Korea, and Hong et al. (2015) successfully applied this approach to East Asia. In this study, the optimal scheme combination of Noah-MP determined by Hong et al. (2015) was used to extract the variables related to the climate regulation function of vegetated surfaces over East Asia.

Data and study domain

To operate Noah-MP, we obtained atmospheric forcing data from the Global Land Data Assimilation System (GLDAS) developed by NASA (Rodell et al., 2004). This data comprises reanalysis data obtained by assimilating atmospheric data from various operational forecasts and observations, for example, from data from the Global Data Assimilation System National Oceanic and Atmospheric Administration of USA. The list of six inputs as atmospheric forcings from GLDAS and outputs as variables in land surface processes is shown in Table 1. The temporal and spatial resolution of the final Noah-MP outputs are 3-hour and 0.25° (approximately 25 km), which are the same as those of GLDAS atmospheric forcing data. The GLDAS data from 1960 onwards are available at ldas.gsfc.nasa.gov.

This model requires land cover and soil texture information over the study area as static inputs. Land cover information, in particular, is used to characterize the biogeochemical features over vegetated covers. Land cover information was from the land use and land cover data used by Kang et al. (2010) (hereafter land cover data from Kongu national university [KLC]), which were generated from the normalized difference vegetation index data of the Moderate Resolution Image Spectro-radiometer (Fig. 1). The land cover classification method used in KLC follows that of the International Geosphere-Biosphere Programme (Friedl et al., 2002). The soil texture data that follow a typical triangular grain-size-based classification were introduced from the Food and Agricultural Organization (FAO, 2002). The study domain was set from −10 to 90° and from 90 to 150° E in latitude and longitude, respectively.

The Noah-MP simulations begin with arbitrary initial conditions regarding the surface status, such as initial soil moisture, SH, LH, and water table depth. Thus, a spin-up process, which is a warm-up simulation for a certain period, is required for the surface conditions to reach equilibrium. With the exception of water table depth, the surface conditions in Noah-MP typically reach equilibrium after an approximately 10-year spin-up (Cai et al., 2014). Thus, in this study, a 60-year simulation was performed from 1951 to 2010 and only 50-year model outputs from 1961 was used for the analyses. The outputs from the initial 10-year period were excluded as a spin-up.

Quantification of climate regulation

Three main components can be considered to quantify the contributions of the surface to climate change. The first is the effect of surface fluxes related to greenhouse gases (ex. CO2, NO2, CH4, etc) absorbed or emitted through biogeochemical processes at the surface. A larger amount of greenhouse gas absorbed by the ecosystem leads to higher functional improvement of the climate regulation. Of the greenhouse gases mentioned above, carbon fluxes at the terrestrial ecosystem can be presented by the Net Ecosystem Exchange (NEE). NEE represents the net amount of emitted or absorbed carbon per unit area and time at the surface. NEE can be converted into radiative forcing for CO2 as follows (Anderson-Teixeira & DeLucia, 2011).

(1)
G CO 2 = E C ( t ) Ddt

GCO2 is the equivalent radiation forcing of CO2 gas at a given time step, and E is the effective radiative efficiency of CO2 (1.4×104 nWm-2ppb-1) (Anderson-Teixeira & DeLucia, 2011; Forster et al., 2007). C(t) represents the additional CO2 concentration in the atmosphere as a rate of the CO2 moles of emitted from the ecosystem surface to the atmosphere. D denotes the CO2 pulse decay over time. In this study, the annual variations in climate regulation were investigated, and thus D becomes about 0.875 (k mol ha-1 yr-1). E is the effective radiative efficiency of CO2 (1.4×104) (Foster et al., 2007).

The second component is net radiation (Rn), which represents the net amount of incoming and outgoing thermal energy fluxes at the surface. This component is a biogeophysical contributor. Rn is determined mainly by albedo, the radiation reflectivity of the surface. Since albedo is determined by the characteristics of vegetation covers, terrestrial ecosystems play an important role in the energy exchange. Rn is quantified by the net sum of incoming and outgoing radiation that are balanced with transferred energy forms such as SH, LH, and GH. The variance of Rn (∆Rn) during a given time period represents the total residual amount in the atmosphere due to energy transfer processes at the surface, which can be categorized as positive or negative feedback. Thus, a higher ∆Rn value indicates a greater degradation of the climate regulation function.

LH can be considered as the third contributor to the climate regulation function. It represents the surface energy flux due to the evapotranspiring process that includes direct evaporation at soil surfaces and transpiration from vegetation. Transpiration is accompanied with photosynthetic activities. Since the water vapors emitted to the atmosphere via this process are resupplied to the land surface and lead to cooling effects, LH can be considered as a positive effect in terms of the climate regulation function of terrestrial ecosystems.

Anderson-Texieira et al. (2012) proposed an equation to quantitatively represent the contribution of terrestrial ecosystems to climate change by combining the three aforementioned components as follows.

(2)
Δ F CR = Δ G CO 2 - Δ R n + Δ LH

Positive values of ∆FCR, the variance of climate regulation function of terrestrial ecosystems, indicate that the ecosystem is functioning to regulate the increase in global warming. The quantitative values of ∆FCR can be numerically calculated when all three components are in the same unit as radiative flux (Wm-2).

Results|Discussion

General trends of climate change over East Asia

Fig. 2 describes the spatio-temporal climate-change trends over East Asia which was from the temperature and PRCP data used for the Noah-MP forcing inputs. Over 50 years, most East Asia regions have exhibited temperature increases, but higher latitudes have shown higher warming trends; regions near the Arctic have shown temperature increases up to 0.05°C per year (Fig. 2A). The temperature variabilities, however, have decreased at high latitudes, indicating that the rates of increase of the annual minimum temperatures are higher than those of the maximum temperatures (Fig. 2B). It is assumed that these increasing trends but decreasing variabilities in temperature at high latitudes are due to higher temperature increases during winter and spring than those during summer and fall (not shown), which reduces inter-seasonal differences in temperature.

Variations in PRCP have been relatively high at 30 to 50° latitudes for 50 years (Fig. 2C): decreasing trends in northern part of China and southern part of Japan and increasing trends in southern part of the Gobi Desert and Korean Peninsula. The variabilities in PRCP, on the contrary, have increased in the Korean Peninsula and southern part of Japan (Fig. 2D). Interestingly, considering the contrasting PRCP trends over these two regions, the PRCP regimes have progressed very differently from each other. For example, the PRCP information can indicate the increasing possibility of extreme events such as heavy rainfall and severe drought in the Korean Peninsula and Japan, respectively. This indication of the variations in extreme events over the two regions must be verified further with more PRCP information at higher spatial and temporal resolutions. The variational trends over the tropics, however, are not significant, considering the very high averages in annual total PRCP.

Variations in quantified climate regulation from terrestrial ecosystem over East Asia

The spatial distributions of the 50-year trends in the quantified climate regulation due to terrestrial ecosystems over East Asia is shown in Fig. 3A. The final outputs were obtained from equation 2 after the quantification of each component simulated using the Noah-MP model. The regions showing an increasing trend of FCR are high latitude regions around Siberia, the Korean Peninsula, Japan Islands, and southeast part of China around the Yangzi River. Meanwhile, those showing a decreasing trend are mainly in the tropics at low latitudes. Interestingly, around the Gobi Desert, the FCR decreases eastward but increases southward.

The significance of the variation in FCR can be examined by calculating it as a percentage of the 50-year average for each region (Fig. 3B). For instance, the increasing trends over Siberian regions are significant, showing increases up to 0.4% per year, which arithmetically indicates an increase over 20% for the past 50 years. Unfortunately, the embedded biogeochemical processes in Noah-MP cannot represent possible offsetting effects from greenhouse gas fluxes via the warming and melting processes of permafrost such as methane. This may result in an overestimation in FCR over such regions. Nevertheless, assuming that the atmospheric variations used for the model inputs are reflected by those in terrestrial ecosystems with climate change, the ecological changes and their impact on the atmosphere have been apparent over the regions. In addition, the differences in the trends FCR between the eastern and southern parts of the Gobi Desert clearly imply the eastward expansion and southward shrinkage of desert areas. This is due to LH effect caused by change in precipitation regimes in the regions. However, it is unclear whether the variations reflect ecological changes over the regions due to very small magnitudes of FCR. The tropics showed high decreasing trends in FCR, but their impacts may not have been significant, comparing to the ordinary fluctuation scales of climate regulation function.

The 50-year trends of the three components, GCO2, LH, and Rn, are shown in Fig. 4. The distribution pattern of GCO2 is very similar to that of the total FCR. This implies that the enhancement in FCR over Siberian regions resulted mainly from biogeochemical processes. Note that some areas in Fig. 4A are blank because some vegetation schemes, such as desert land covers like the Gobi Desert, were not activated in Noah-MP. The contributions from LH to FCR have been relatively low in Siberian regions (Fig. 4B). On the contrary, regions in the Korean Peninsula and around the border between China and Russia around the southeastern part of Siberia show increases in LH. Contributions of variations in net radiation derived mainly from biogeophysical processes have been relatively low over the entire East Asia region (Fig. 4C). Note that the opposite direction of color bar due to the opposite feedback to the atmosphere as the definition of climate regulation mentioned above. Thus, the contributions to the total FCR from the three components can be arranged as follows, in decreasing order: GCO2, LH, and Rn. In other words, warming over East Asia has led to biogeochemical processes at the land surfaces at high latitudes.

The averaged 50-year trend of FCR for each land cover type is shown in Fig. 5. The trends of GCO2 and LH reflect the regional characteristics of terrestrial ecosystems. Overall, majority of the land covers showed increasing trends of GCO2 and LH, except evergreen broadleaf forests that are mainly distributed in the tropics. “Mixed shrub/grassian,” and “deciduous needleleaf forest” that are the major types of land cover in the Siberian regions showed increases in both components. The opposite phenomena over the “grassian” type of land cover that is mainly distributed around Gobi Desert is also related to the changes in the desert characteristics as mentioned above. Variations in Rn for land cover types have been relatively small, but most of them showed increasing trends, which has negatively affected the total FCR.

In conclusion, Spatio-temporal variations in FCR and its components, GCO2, LH, and Rn, at the land surface were investigated to see their contributions to the past global warming. Those components were forced by atmospheric changes over East Asia in a LSM that represents the interactions between terrestrial ecosystems and the atmosphere. Atmospheric changes according to GLDAS data indicated that the warming trends during the last half century have been stronger at higher latitudes. Strong climatic signals in PRCP variations were also shown in regions between 30° and 50° latitudes. The simulated FCR forced by these atmospheric variations showed significant increasing trends in land covers over high latitudes, and the main contributor was GCO2 via biogeochemical processes. These increases in FCR over the regions suggest the possibility of negative feedback to global warming accompanied by surface ecological changes. In other words, it can be concluded that global warming has stimulated biogeochemical processes at high latitudes, but the terrestrial ecosystems have worked to decrease the warming rate by increasing carbon assimilation. The model also described the variations in FCR around the Gobi Desert, reflecting the possible eastward movement of desert areas.

These results showed that the FCR simulation is sensitive to climate change and the characteristics of land cover types. This indicates that the accuracy of the model can be improved by using accurate information of spatial distributions of land covers, accurate parameterizations representing biogeochemical processes, and application of a long term land cover change model. The current model used in this study cannot reflect ecological successions and spatial changes in climate zones due to temporally fixed land covers. Another limitation of the LSM used in this study is that it cannot represent the possible effects of melted permafrost over high latitudes. Since this study considered only feedback from CO2 gas involved in respiration at the surface, the increase in FCR over high latitudes may be overestimated. The characterization of microbial processes involved in carbon fluxes at the soil surface must be improved.

This study demonstrated how the climate regulation function of terrestrial ecosystems has regionally changed with climate change for the last half century. Climate regulation due to terrestrial ecosystem varies with climatic factors in the atmosphere, and thus better understanding of the feedback from the terrestrial ecosystem to the atmosphere helps develop climate prediction models with higher accuracy. In addition, since regional characteristics of terrestrial ecosystems determine the climate regulation function, the approach used in this study can provide useful information to devise countermeasures for climate change mitigation at the regional scale. The information can be extracted by checking out the regional differences and variations in climate regulation components between the past 50 years and more recent ones in East Asia. For example, predicting changes in climate regulation due to future climate change can be valuable in creating land use policies to improve the climate regulation function and also serve as strategic information for international communities to mitigate climate change.

Acknowledgments

This study was funded by the Assessment of Climate Change Risk for Ecosystems in Korea (NIE-BR-2020-11) from the National Institute of Ecology.

Conflict of Interest

The authors declare that they have no competing interests.

References

1 

Anderson-Teixeira, K.J., & DeLucia, E.H. (2011) The greenhouse gas value of ecosystems Global Change Biology, 17, 425-438 .

2 

Anderson-Teixeira, K.J., Miller, A.D., Mohan, J.E., Hudiburg, T.W., Duval, B.D., & Delucia, E.H. (2013) Altered dynamics of forest recovery under a changing climate Global Change Biology, 19, 2001-2021 .

3 

Anderson-Teixeira, K.J., Snyder, P.K., Twine, T.E., Cuadra, S.V., Costa, M.H., & DeLucia, E.H. (2012) Climate-regulation services of natural and agricultural ecoregions of the Americas Nature Climate Change, 2, 177-181 .

4 

Cai, X., Yang, Z.L., David, C.H., Niu, G.Y., & Rodell, M. (2014) Hydrological evaluation of the Noah-MP land surface model for the Mississippi River Basin. Journal of Geophysical Research: Atmospheres, 119, 23-38 .

5 

Castro, L.F.R. (2014) Master's thesis Kiel: University of Kiel Ecosystem service assessment for the new vision of the energy project in Samsø (Denmark)- an analysis of the climate regulation service supply

6 

Chen, F., & Dudhia, J. (2001) Coupling an advanced land surface-hydrology model with the Penn State-NCAR MM5 modeling system. Part I: model implementation and sensitivity Monthly Weather Review, 129, 569-585 .

7 

Chen, F., Mitchell, K., Schaake, J., Xue, Y., Pan, H.L., Koren, V., et al. (1996) Modeling of land surface evaporation by four schemes and comparison with FIFE observations. Journal of Geophysical Research: Atmospheres, 101, 7251-7268 .

8 

Daly, C., Bachelet, D., Lenihan, J.M., Neilson, R.P., Parton, W., & Ojima, D. (2000) Dynamic simulation of tree-grass interactions for global change studies Ecological Applications, 10, 449-469 .

9 

Deardorff, J.W. (1978) Efficient prediction of ground surface temperature and moisture, with inclusion of a layer of vegetation Journal of Geophysical Research, 83, 1889-1903 .

10 

de Groot, R.S., Wilson, M.A., & Boumans, R.M.J. (2002) A typology for the classification, description and valuation of ecosystem functions, goods and services Ecological Economics, 41, 393-408 .

11 

Dickinson, R.E., Henderson-Sellers, A., Kennedy, P.J., & Wilson, M.F. (1986) Biosphere-Atmosphere Transfer Scheme (BATS) for the NCAR Community Climate Model Boulder: National Center for Atmospheric Research

12 

Dickinson, R.E., Shaikh, M., Bryant, R., & Graumlich, L. (1998) Interactive canopies for a climate model Journal of Climate, 11, 2823-2836 .

13 

Eagleson, P.S. (1982) Land Surface Processes in Atmospheric General Circulation Models Cambridge: Cambridge University Press

14 

FAO (2003) Digital Soil Map of the World and Derived Soil Properties Rome: FAO, Land and Water Development Division p. Land

15 

Feddema, J.J., Oleson, K.W., Bonan, G.B., Mearns, L.O., Buja, L.E., Meehl, G.A., et al. (2005) The importance of land-cover change in simulating future climates Science, 310, 1674-1678 .

16 

Forster, P., Ramaswamy, V., Artaxo, P., Berntsen, T., Betts, R., Fahey, D.W., et al., Solomon, S.D., Qin, M., Manning, Z., Chen, M., Marquis, K.B., Averyt, M., et al. (Eds.) (2007) Climate Change 2007: The Physical Science Basis New York: Cambridge University Press Changes in atmospheric constituents and radiative forcing, pp. 129-234

17 

Friedl, M.A., McIver, D.K., Hodges, J.C.F., Zhang, X.Y., Muchoney, D., Strahler, A.H., et al. (2002) Global land cover mapping from MODIS: algorithms and early results Remote Sensing of Environment, 83, 287-302 .

18 

Guimarães, H., Braga, R., Mascarenhas, A., & Ramos, T.B. (2017) Indicators of ecosystem services in a military Atlantic forest area, Pernambuco-Brazil Ecological Indicators, 80, 247-257 .

19 

Hong, S., Lakshmi, V., & Small, E.E. (2007) Relationship between vegetation biophysical properties and surface temperature using multisensor satellite data Journal of Climate, 20, 5593-5606 .

20 

Hong, S., Park, S.K., & Yu, X. (2015) Scheme-based optimization of land surface model using a micro-genetic algorithm: assessment of its performance and usability for regional applications SOLA, 11, 129-133 .

21 

Hong, S., Yu, X., Park, S.K., Choi, Y.S., & Myoung, B. (2014) Assessing optimal set of implemented physical parameterization schemes in a multi-physics land surface model using genetic algorithm Geoscientific Model Development, 7, 2517-2529 .

22 

Kang, J.H., Suh, M.S., & Kwak, C.H. (2010) Land cover classification over East Asian Region using recent MODIS NDVI data (2006-2008) Atmosphere, 20, 415-426 .

23 

Kim, J.U., Lee, J., Boo, K.O., Shim, S., Kim, J.E., & Byun, Y.H. (2016) Effects of continental evaporation for precipitation over East Asia in the past and the future of HadGEM2-AO climate model Atmosphere, 26, 553-563 .

24 

Lakshmi, V., Hong, S., Small, E.E., & Chen, F. (2011) The influence of the land surface on hydrometeorology and ecology: new advances from modeling and satellite remote sensing Hydrology Research, 42, 95-112 .

25 

Liu, S., Costanza, R., Troy, A., D'Aagostino, J., & Mates, W. (2010) Valuing New Jersey's ecosystem services and natural capital: a spatially explicit benefit transfer approach Environmental Management, 45, 1271-1285 .

26 

Ma, L., Bicking, S., & Müller, F. (2019) Mapping and comparing ecosystem service indicators of global climate regulation in Schleswig-Holstein, Northern Germany The Science of the Total Environment, 648, 1582-1597 .

27 

Manabe, S. (1969) Climate and the ocean circulation, I, The atmospheric circulation and the hydrology of the Earth's surface Monthly Weather Review, 97, 739-774 .

28 

Millennium Ecosystem Assessment (2005) Ecosystems and Human Well-Being. Summary, Our Human Planet: Summary for Decision Makers Washington: Island Press

29 

Myoung, B., Choi, Y., & Park, S.K. (2011) A review on vegetation models and applicability to climate simulations at regional scale Asia-Pacific Journal of Atmospheric Sciences, 47, 463-475 .

30 

Myoung, B., Choi, Y.S., Hong, S., & Park, S.K. (2013) Inter- and intra-annual variability of vegetation in the Northern Hemisphere and its association with precursory meteorological factors Global Biogeochemical Cycles, 27, 31-42 .

31 

Niu, G.Y., & Yang, Z.L. (2004) Effects of vegetation canopy processes on snow surface energy and mass balances. Journal of Geophysical Research: Atmospheres, 109, D23111 .

32 

Niu, G.Y., & Yang, Z.L. (2006) Effects of frozen soil on snowmelt runoff and soil water storage at a continental scale Journal of Hydrometeorology, 7, 937-952 .

33 

Niu, G.Y., Yang, Z.L., Dickinson, R.E., & Gulden, L.E. (2005) A simple TOPMODEL-based runoff parameterization (SIMTOP) for use in global climate models. Journal of Geophysical Research: Atmospheres, 110, D21106 .

34 

Niu, G.Y., Yang, Z.L., Mitchell, K.E., Chen, F., Ek, M.B., Barlage, M., et al. (2011) The community Noah land surface model with multiparameterization options (Noah-MP): 1. Model description and evaluation with local-scale measurements. Journal of Geophysical Research: Atmospheres, 116, D12109 .

35 

Pan, H., & Mahrt, L. (1987) Interaction between soil hydrology and boundary-layer development Boundary-Layer Meteorology, 38, 185-202 .

36 

Petter, M., Mooney, S., Maynard, S.M., Davidson, A., Cox, M., & Horosak, I. (2012) A methodology to map ecosystem functions to support ecosystem services assessments Ecology and Society, 18, 31 .

37 

Prentice, I.C., Heimann, M., & Sitch, S. (2000) The carbon balance of the terrestrial biosphere: ecosystem models and atmospheric observations Ecological Applications, 10, 1553-1573 .

38 

Rodell, M., Houser, P.R., Jambor, U., Gottschalck, J., Mitchell, K., Meng, C.J., et al. (2004) The global land data assimilation system Bulletin of the American Meteorological Society, 85, 381-394 .

39 

Running, S.W., & Coughlan, J.C. (1988) A general model of forest ecosystem processes for regional applications I. Hydrologic balance, canopy gas exchange and primary production processes Ecological Modelling, 42, 125-154 .

40 

Scholes, R.J. (2016) Climate change and ecosystem services WIREs Climate Change, 7, 537-550 .

41 

Snyder, P.K., Delire, C., & Foley, J.A. (2004) Evaluating the influence of different vegetation biomes on the global climate Climate Dynamics, 23, 279-302 .

42 

West, P.C., Narisma, G.T., Barford, C.C., Kucharik, C.J., & Foley, J.A. (2011) An alternative approach for quantifying climate regulation by ecosystems Frontiers in Ecology and the Environment, 9, 126-133 .

43 

Wood, S.L.R., Jones, S.K., Johnson, J.A., Brauman, K.A., Chaplin-Kramer, R., Fremier, A., et al. (2018) Distilling the role of ecosystem services in the Sustainable Development Goals Ecosystem Services, 29, 70-82 .

44 

Yang, Z.L., & Niu, G.Y. (2003) The Versatile Integrator of Surface and Atmosphere processes: Part 1. Model description Global and Planetary Change, 38, 175-189 .

45 

Yang, Z.L., Niu, G.Y., Mitchell, K.E., Chen, F., Ek, M.B., Barlage, M., et al. (2011) The community Noah land surface model with multiparameterization options (Noah-MP): 2. Evaluation over global river basins. Journal of Geophysical Research: Atmospheres, 116, D12110 .

46 

Yu, X.X., Qin, Y.S., Chen, L.H., & Liu, S. (2002) The forest ecosystem services and their valuation of Beijing mountain areas Acta Ecologica Sinica, 22, 783-786 .

Figures and Table
새창으로 보기
Fig. 1

Land use and land cover map from Kongju National University, reinterpreted based on the United States Geological Survey land cover classification criteria (Kang et al., 2010).

PNIE-1-058-f1.jpg
새창으로 보기
Fig. 2

Spatial distributions of (A) Ta trend, (B) trend of Ta variability (10-year running standard deviations), (C) PRCP trend as percentage of local PRCP averages, and (D) trend of PRCP variability in percentage (10-year running standard deviations) for recent 50 years (1961-2010). PRCP, precipitation.

PNIE-1-058-f2.jpg
새창으로 보기
Fig. 3

Spatial distributions of the simulated FCR trend of the terrestrial ecosystem quantified (A) in the energy flux unit and (B) in percentages of local FCR averages for recent 50 years.

PNIE-1-058-f3.jpg
새창으로 보기
Fig. 4

Spatial distributions of the FCR components: (A) carbon flux (GCO2), (B) moisture flux (LH), and (C) net radiation flux (Rn). LH, latent heat flux.

PNIE-1-058-f4.jpg
새창으로 보기
Fig. 5

Trend of FCR components: (A) carbon flux (GCO2), (B) moisture flux (LH), and (C) net radiation flux (Rn) for each land cover type in whisker box plots with minimum, lower quartile, median, upper quartile, and maximum (extracted only for grids showing 0.05 P-value and lower). LH, latent heat flux.

PNIE-1-058-f5.jpg
새창으로 보기
Table 1

List of inputs and outputs of Noah-MP model

Inputs (atmospheric forcings) Outputs (variables in land surface processes)
Precipitation
Incoming solar radiations
Air temperature (2 m height from the surface)
Humidity (2 m height from the surface)
Wind speed (10 m height from the surface)
Surface pressure
Energy cycle
 Soil heat flux, Latent heat flux, ground heat flux
Water cycle
 Evapotranspiration, surface and subsurface runoff
Geophysical features
 Surface temperature, soil temperature and moisture at various layers,
  snow cover and melting rate
Plant physiological features
 GPP, NPP, leaf area index, vegetation fraction
[i]

Noah-MP, Noah land surface model with multi-parameterization options; GPP, gross primary productivity; NPP, net primary productivity.