1
|
Pak U, Guo Q, Liu Z, Wang X, Liu Y, Jin G. Spatial Distribution of Pinus koraiensis Trees and Community-Level Spatial Associations in Broad-Leaved Korean Pine Mixed Forests in Northeastern China. PLANTS (BASEL, SWITZERLAND) 2023; 12:2906. [PMID: 37631117 PMCID: PMC10459911 DOI: 10.3390/plants12162906] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/16/2023] [Revised: 07/28/2023] [Accepted: 08/04/2023] [Indexed: 08/27/2023]
Abstract
Investigating the spatial distributions and associations of tree populations provides better insights into the dynamics and processes that shape the forest community. Korean pine (Pinus koraiensis) is one of the most important tree species in broad-leaved Korean pine mixed forests (BKMFs), and little is known about the spatial point patterns of and associations between Korean pine and community-level woody species groups such as coniferous and deciduous trees in different developmental stages. This study investigated the spatial patterns of Korean pine (KP) trees and then analyzed how the spatial associations between KP trees and other tree species at the community level vary in different BKMFs. Extensive data collected from five relatively large sample plots, covering a substantial area within the natural distribution range of KP in northeastern China, were utilized. Uni- and bivariate pair correlation functions and mark correlation functions were applied to analyze spatial distribution patterns and spatial associations. The DBH (diameter at breast height) histogram of KP trees in northeastern China revealed that the regeneration process was very poor in the Changbai Mountain (CBS) plot, while the other four plots exhibited moderate or expanding population structures. KP trees were significantly aggregated at scales up to 10 m under the HPP null model, and the aggregation scales decreased with the increase in size classes. Positive or negative spatial associations were observed among different life stages of KP trees in different plots. The life history stages of the coniferous tree group showed positive spatial associations with KP saplings and juvenile trees at small scales, and spatial independence or negative correlations with larger KP trees at greater scales. All broad-leaved tree groups (canopy, middle, and understory layers) exhibited only slightly positive associations with KP trees at small scales, and dominant negative associations were observed at most scales. Our results demonstrate that mature KP trees have strong importance in the spatial patterns of KP populations, and site heterogeneity, limited seed dispersal, and interspecific competition characterize the spatial patterns of KP trees and community-level spatial associations with respect to KP trees, which can serve as a theoretical basis for the management and restoration of BKMFs in northeastern China.
Collapse
Affiliation(s)
- Unil Pak
- Center for Ecological Research, Northeast Forestry University, Harbin 150040, China; (U.P.); (Q.G.); (Z.L.)
| | - Qingxi Guo
- Center for Ecological Research, Northeast Forestry University, Harbin 150040, China; (U.P.); (Q.G.); (Z.L.)
- Key Laboratory of Sustainable Forest Ecosystem Management-Ministry of Education, Northeast Forestry University, Harbin 150040, China
| | - Zhili Liu
- Center for Ecological Research, Northeast Forestry University, Harbin 150040, China; (U.P.); (Q.G.); (Z.L.)
- Key Laboratory of Sustainable Forest Ecosystem Management-Ministry of Education, Northeast Forestry University, Harbin 150040, China
- Northeast Asia Biodiversity Research Center, Northeast Forestry University, Harbin 150040, China
| | - Xugao Wang
- CAS Key Laboratory of Forest Ecology and Management, Institute of Applied Ecology, Chinese Academy of Sciences, Shenyang 110164, China;
| | - Yankun Liu
- Heilongjiang Forestry Engineering and Environment Institute, Harbin 150040, China;
- Key Laboratory of Forest Ecology and Forestry Ecological Engineering of Heilongjiang Province, Harbin 150040, China
| | - Guangze Jin
- Center for Ecological Research, Northeast Forestry University, Harbin 150040, China; (U.P.); (Q.G.); (Z.L.)
- Key Laboratory of Sustainable Forest Ecosystem Management-Ministry of Education, Northeast Forestry University, Harbin 150040, China
- Northeast Asia Biodiversity Research Center, Northeast Forestry University, Harbin 150040, China
| |
Collapse
|
2
|
Abstract
Satellite-based normalized difference vegetation index (NDVI) time series data are useful for monitoring the changes in vegetation ecosystems in the context of global climate change. However, most of the current NDVI products cannot effectively reconcile high spatial resolution and continuous observations in time. Here, to produce a global-scale, long-term, and high-resolution NDVI database, we developed a simple and new data downscaling approach. The downscaling algorithm considers the pixel-wise ratios of the coefficient of variation (CV) between the coarse- and fine-resolution NDVI data and relative changes in the NDVI against a baseline period. The algorithm successfully created a worldwide monthly NDVI database with 250 m resolution from 1982 to 2018 by translating the fine spatial information from MODIS (Moderate-resolution Imaging Spectroradiometer) data and the long-term temporal information from AVHRR (Advanced Very High Resolution Radiometer) data. We employed the evaluation indices of root mean square error (RMSE), mean absolute error (MAE), and Pearson’s correlation coefficient (Pearson’s R) to assess the accuracy of the downscaled data against the MODIS NDVI. Both the RMSE and MAE values at the regional and global scales are typically between 0 and 0.2, whereas the Pearson’s R values are mostly above 0.7, which implies that the downscaled NDVI product is similar to the MODIS NDVI product. We then used the downscaled data to monitor the NDVI changes in different plant types and places with significant vegetation heterogeneity, as well as to investigate global vegetation trends over the last four decades. The Google Earth Engine platform was used for all the data downscaling processes, and here we provide a code for users to easily acquire data corresponding to any part of the world. The downscaled global-scale NDVI time series has high potential for the monitoring of the long-term temporal and spatial dynamics of terrestrial ecosystems under changing environments.
Collapse
|
3
|
Forest Canopy Height Mapping by Synergizing ICESat-2, Sentinel-1, Sentinel-2 and Topographic Information Based on Machine Learning Methods. REMOTE SENSING 2022. [DOI: 10.3390/rs14020364] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
Spaceborne LiDAR has been widely used to obtain forest canopy heights over large areas, but it is still a challenge to obtain spatio-continuous forest canopy heights with this technology. In order to make up for this deficiency and take advantage of the complementary for multi-source remote sensing data in forest canopy height mapping, a new method to estimate forest canopy height was proposed by synergizing the spaceborne LiDAR (ICESat-2) data, Synthetic Aperture Radar (SAR) data, multi-spectral images, and topographic data considering forest types. In this study, National Geographical Condition Monitoring (NGCM) data was used to extract the distributions of coniferous forest (CF), broadleaf forest (BF), and mixed forest (MF) in Hua’ nan forest area in Heilongjiang Province, China. Accordingly, the forest canopy height estimation models for whole forest (all forests together without distinguishing types, WF), CF, BF, and MF were established, respectively, by Radom Forest (RF) and Gradient Boosting Decision Tree (GBDT). The accuracy for established models and the forest canopy height obtained based on estimation models were validated consequently. The results showed that the forest canopy height estimation models considering forest types had better performance than the model grouping all types of forest together. Compared with GBDT, RF with optimal variables had better performance in forest canopy height estimation with Pearson’s correlation coefficient (R) and the root-mean-squared error (RMSE) values for CF, BF, and MF of 0.72, 0.59, 0.62, and 3.15, 3.37, 3.26 m, respectively. It has been validated that a synergy of ICESat-2 with other remote sensing data can make a crucial contribution to spatio-continuous forest canopy height mapping, especially for areas covered by different types of forest.
Collapse
|
4
|
Huang C, Yang Q, Huang W. Analysis of the Spatial and Temporal Changes of NDVI and Its Driving Factors in the Wei and Jing River Basins. INTERNATIONAL JOURNAL OF ENVIRONMENTAL RESEARCH AND PUBLIC HEALTH 2021; 18:11863. [PMID: 34831620 PMCID: PMC8618191 DOI: 10.3390/ijerph182211863] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 09/23/2021] [Revised: 11/08/2021] [Accepted: 11/08/2021] [Indexed: 12/02/2022]
Abstract
This study aimed to explore the long-term vegetation cover change and its driving factors in the typical watershed of the Yellow River Basin. This research was based on the Google Earth Engine (GEE), a remote sensing cloud platform, and used the Landsat surface reflectance datasets and the Pearson correlation method to analyze the vegetation conditions in the areas above Xianyang on the Wei River and above Zhangjiashan on the Jing River. Random forest and decision tree models were used to analyze the effects of various climatic factors (precipitation, temperature, soil moisture, evapotranspiration, and drought index) on NDVI (normalized difference vegetation index). Then, based on the residual analysis method, the effects of human activities on NDVI were explored. The results showed that: (1) From 1987 to 2018, the NDVI of the two watersheds showed an increasing trend; in particular, after 2008, the average increase rate of NDVI in the growing season (April to September) increased from 0.0032/a and 0.003/a in the base period (1987-2008) to 0.0172/a and 0.01/a in the measurement period (2008-2018), for the Wei and Jing basins, respectively. In addition, the NDVI significantly increased from 21.78% and 31.32% in the baseline period (1987-2008) to 83.76% and 92.40% in the measurement period (2008-2018), respectively. (2) The random forest and classification and regression tree model (CART) can assess the contribution and sensitivity of various climate factors to NDVI. Precipitation, soil moisture, and temperature were found to be the three main factors that affect the NDVI of the study area, and their contributions were 37.05%, 26.42%, and 15.72%, respectively. The changes in precipitation and soil moisture in the entire Jing River Basin and the upper and middle reaches of the Wei River above Xianyang caused significant changes in NDVI. Furthermore, changes in precipitation and temperature led to significant changes in NDVI in the lower reaches of the Wei River. (3) The impact of human activities in the Wei and Jing basins on NDVI has gradually changed from negative to positive, which is mainly due to the implementation of soil and water conservation measures. The proportions of areas with positive effects of human activities were 80.88% and 81.95%, of which the proportions of areas with significant positive effects were 11.63% and 7.76%, respectively. These are mainly distributed in the upper reaches of the Wei River and the western and eastern regions of the Jing River. These areas are the key areas where soil and water conservation measures have been implemented in recent years, and the corresponding land use has transformed from cultivated land to forest and grassland. The negative effects accounted for 1.66% and 0.10% of the area, respectively, and were mainly caused by urban expansion and coal mining.
Collapse
Affiliation(s)
- Chenlu Huang
- College of Tourist (Institute of Human Geography), Xi’an International Studies University, Xi’an 710127, China;
| | - Qinke Yang
- College of Urban and Environment Sciences, Northwest University, Xi’an 710127, China
| | - Weidong Huang
- Hydrology and Water Resources Bureau of Gansu Province, Lanzhou 730000, China;
| |
Collapse
|
5
|
Wang W, Sun J, Zhong Z, Xiao L, Wang Y, Wang H. Relating macrofungal diversity and forest characteristics in boreal forests in China: Conservation effects, inter-forest-type variations, and association decoupling. Ecol Evol 2021; 11:13268-13282. [PMID: 34646468 PMCID: PMC8495802 DOI: 10.1002/ece3.8049] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/26/2021] [Revised: 08/04/2021] [Accepted: 08/09/2021] [Indexed: 11/11/2022] Open
Abstract
QUESTION How conservation and forest type affect macrofungal compositional diversity is not well understood. Even less is known about macrofungal associations with plants, soils, and geoclimatic conditions. LOCATION Southern edge of boreal forest distribution in China, named as Huzhong Nature Reserve. METHODS We surveyed a total of 72 plots for recording macrofungi, plants, and topography in 2015 and measured soil organic carbon, nitrogen, and bulk density. Effects of conservation and forest types on macrofungi and plants were compared, and their associations were decoupled by structural equation modeling (SEM) and redundancy ordination (RDA). RESULTS Conservation and forest type largely shaped macrofungal diversity. Most of the macrofungal traits declined with the conservation intensities or peaked at the middle conservation region. Similarly, 91% of macrofungal traits declined or peaked in the middle succession stage of birch-larch forests. Forest conservation resulted in the observation of sparse, larch-dominant, larger tree forests. Moreover, the soil outside the Reserve had more water, higher fertility, and lower bulk density, showing miscellaneous wood forest preference. There is a complex association between conservation site characteristics, soils, plants, and macrofungi. Variation partitioning showed that soil N was the top-one factor explaining the macrofungal variations (10%). As shown in SEM coefficients, conservation effect to macrofungi (1.1-1.2, p < .05) was like those from soils (1.2-1.6, p < .05), but much larger than the effect from plants (0.01-0.14, p > .10). For all tested macrofungal traits, 89%-97% of their variations were from soils, and 5%-21% were from conservation measures, while plants compensated 1%-10% of these effects. Our survey found a total of 207 macrofungal species, and 65 of them are new updates in this Reserve, indicating data shortage for the macrofungi list here. CONCLUSION Our findings provide new data for the joint conservation of macrofungi and plant communities, highlighting the crucial importance of soil matrix for macrofungal conservation in boreal forests.
Collapse
Affiliation(s)
- Wenjie Wang
- Urban Forests and Wetlands groupNortheast Institute of Geography and AgroecologyChinese Academy of ScienceChangchunChina
- Key Laboratory of Forest Plant EcologyMinistry of EducationNortheast Forestry UniversityHarbinChina
| | - Jingxue Sun
- Key Laboratory of Forest Plant EcologyMinistry of EducationNortheast Forestry UniversityHarbinChina
| | - Zhaoliang Zhong
- Key Laboratory of Forest Plant EcologyMinistry of EducationNortheast Forestry UniversityHarbinChina
| | - Lu Xiao
- Urban Forests and Wetlands groupNortheast Institute of Geography and AgroecologyChinese Academy of ScienceChangchunChina
| | - Yuanyuan Wang
- Urban Forests and Wetlands groupNortheast Institute of Geography and AgroecologyChinese Academy of ScienceChangchunChina
| | - Huimei Wang
- Key Laboratory of Forest Plant EcologyMinistry of EducationNortheast Forestry UniversityHarbinChina
| |
Collapse
|
6
|
Gao X, Huang X, Lo K, Dang Q, Wen R. Vegetation responses to climate change in the Qilian Mountain Nature Reserve, Northwest China. Glob Ecol Conserv 2021. [DOI: 10.1016/j.gecco.2021.e01698] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022] Open
|
7
|
Wu C, Venevsky S, Sitch S, Mercado LM, Huntingford C, Staver AC. Historical and future global burned area with changing climate and human demography. ACTA ACUST UNITED AC 2021. [DOI: 10.1016/j.oneear.2021.03.002] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
|
8
|
Wang M, Venevsky S, Wu C, Berdnikov S, Sorokina V, Kulygin V. Description of local carbon flux from large scale gridded climate data by a dynamic global vegetation model at variable time steps: Example of Euroflux sites. THE SCIENCE OF THE TOTAL ENVIRONMENT 2021; 756:143492. [PMID: 33302082 DOI: 10.1016/j.scitotenv.2020.143492] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/03/2019] [Revised: 10/26/2020] [Accepted: 10/28/2020] [Indexed: 06/12/2023]
Abstract
Dynamic Global Vegetation Models (DGVMs) are commonly used to describe the land biogeochemical processes and regulate carbon and water pools. However, the simulation efficiency and validation of DGVMs are limited to varying temporal and spatial resolutions. Additionally, the uncertainties caused by different interpolation methods used in DGVMs are still not clear. In this study, we employ Socio-Economic and natural Vegetation ExpeRimental (SEVER) DGVM to simulate Net Ecosystem Exchange (NEE) flux with large scale National Centers for Environmental Prediction (NCEP) daily climate data as inputs for the years 1997-2000 at 14 Euroflux sites. It is shown that daily local NEE flux on chosen sites can be reasonably simulated, and daily temperature and shortwave radiation are the most essential inputs for daily NEE simulation compared with precipitation and the ratio of sunshine hours. Different running means (1 to 30 days) methods are analysed for each Euroflux site, and the best results of both averaged regression coefficient and averaged slope of regression are discovered by using 5 days running mean method. SEVER DGVM, driven by linearly interpolated daily climate data is compared at the monthly time step with Lund-Potsdam-Jena (LPJ) DGVM, which combines the linear interpolation of daily temperature with stochastic generation of daily precipitation. The comparison demonstrates that the stochastic generation of daily precipitation provides an acceptable fit to local observed NEE, but with a slight decrease in accuracy. Simulation experiments with SEVER DGVM demonstrate that daily local NEE flux inside a grid cell for a region as large as Europe can be modelled by DGVMs, using only large scale climate data as inputs.
Collapse
Affiliation(s)
- Menghui Wang
- Ministry of Education Key Laboratory for Earth System Modeling, Department of Earth System Science, Tsinghua University, Beijing 100084, China
| | - Sergey Venevsky
- Ministry of Education Key Laboratory for Earth System Modeling, Department of Earth System Science, Tsinghua University, Beijing 100084, China.
| | - Chao Wu
- Ministry of Education Key Laboratory for Earth System Modeling, Department of Earth System Science, Tsinghua University, Beijing 100084, China; Department of Ecology and Evolutionary Biology, Yale University, New Haven, CT 06511, USA.
| | - Sergey Berdnikov
- Federal Research Centre Тhe Southern Scientific Centre of the Russian Academy of Sciences (SSC RAS), Rostov-on-Don 344006, Russia
| | - Vera Sorokina
- Federal Research Centre Тhe Southern Scientific Centre of the Russian Academy of Sciences (SSC RAS), Rostov-on-Don 344006, Russia
| | - Valerii Kulygin
- Federal Research Centre Тhe Southern Scientific Centre of the Russian Academy of Sciences (SSC RAS), Rostov-on-Don 344006, Russia
| |
Collapse
|
9
|
Spatio-Temporal Vegetation Dynamic and Persistence under Climatic and Anthropogenic Factors. REMOTE SENSING 2020. [DOI: 10.3390/rs12162612] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
Land degradation reflected by vegetation is a commonly used practice to monitor desertification. To retrieve important information for ecosystem management accurate assessment of desertification is necessary. The major factors that drive vegetation dynamics in arid and semi-arid regions are climate and anthropogenic activities. Progression of desertification is expected to exacerbate under future climate change scenarios, through precipitation variability, increased drought frequency and persistence of dry conditions. This study examined spatiotemporal vegetation dynamics in arid regions of Sindh, Pakistan, using annual and growing season Normalized Difference Vegetation Index (NDVI) data from 2000 to 2017, and explored the climatic and anthropogenic effects on vegetation. Results showed an overall upward trend (annual 86.71% and growing season 82.7%) and partial downward trend (annual 13.28% and growing season 17.3%) in the study area. NDVI showed the highest significant increase in cropland region during annual, whereas during growing season the highest significant increase was observed in savannas. Overall high consistency in future vegetation trends in arid regions of Sindh province is observed. Stable and steady development region (annual 48.45% and growing 42.80%) dominates the future vegetation trends. Based on the Hurst exponent and vegetation dynamics of the past, improvement in vegetation cover is predicted for a large area (annual 44.49% and growing 30.77%), and a small area is predicted to have decline in vegetation activity (annual 0.09% and growing 3.04%). Results revealed that vegetation growth in the study area is a combined result of climatic and anthropogenic factors; however, in the future multi-controls are expected to have a slightly larger impact on annual positive development than climate whereas positive development in growing season is more likely to continue in future under the control of climate variability.
Collapse
|
10
|
Ge J, Berg B, Xie Z. Climatic seasonality is linked to the occurrence of the mixed evergreen and deciduous broad‐leaved forests in China. Ecosphere 2019. [DOI: 10.1002/ecs2.2862] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022] Open
Affiliation(s)
- Jielin Ge
- State Key Laboratory of Vegetation and Environmental Change Institute of Botany Chinese Academy of Sciences No. 20 Nanxincun Xiangshan Beijing 100093 China
| | - Björn Berg
- Section of Biology University of Gävle Gavle SE‐80176 Sweden
- Department of Forest Sciences University of Helsinki Helsinki Finland
| | - Zongqiang Xie
- State Key Laboratory of Vegetation and Environmental Change Institute of Botany Chinese Academy of Sciences No. 20 Nanxincun Xiangshan Beijing 100093 China
- University of Chinese Academy of Sciences Beijing 100049 China
| |
Collapse
|
11
|
Chu H, Venevsky S, Wu C, Wang M. NDVI-based vegetation dynamics and its response to climate changes at Amur-Heilongjiang River Basin from 1982 to 2015. THE SCIENCE OF THE TOTAL ENVIRONMENT 2019; 650:2051-2062. [PMID: 30290347 DOI: 10.1016/j.scitotenv.2018.09.115] [Citation(s) in RCA: 93] [Impact Index Per Article: 18.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/14/2018] [Revised: 09/08/2018] [Accepted: 09/08/2018] [Indexed: 06/08/2023]
Abstract
Vegetation in Northern Hemisphere, being sensitive to climate change, plays an important role in the carbon cycles between land and the atmosphere. The response of vegetation to climate change was analyzed at pixel, biome and regional scale in Amur-Heilongjiang River Basin (AHRB) for growing season, spring, summer and autumn using Normalized Difference Vegetation Index and gridded climate data for the period 1982-2015. NDVI and climate variables trend detection methods and correlation analysis were applied. The potential impacts of human activities on growing season NDVI dynamics were investigated further using residual trend analysis. Results showed that at river basin scale, growing season vegetation experienced a discontinuous greening trend with two reversals, demonstrating that NDVI initially increased to mid-1990s, then declined to mid-2000s, and finally rebounded to 2015. This may be attributed to the shifting between drought and wet trends, indicating growing season NDVI was mainly regulated by precipitation. Temperature was the dominant factor on affecting spring vegetation growth while autumn NDVI showed negative correlation with precipitation due to the relation of precipitation with sunshine hours available for photosynthesis. The response of vegetation growth to climatic variations varied among vegetation types. Grassland NDVI exhibited positive correlation with precipitation in all time ranges. NDVI of needleleaved forest, broadleaved forest, mixed forest and woodland were positively correlated with temperature in all seasons, while showing significant negative correlation with autumn precipitation. Residual trend analysis revealed that human activities might lead to the vegetation degradation in China farming zone of AHRB. Fires also play an important role in regulating vegetation dynamics in the region. Results of our analysis can be used by national governments from three countries of AHRB in managing and negotiating vegetation resources of the region.
Collapse
Affiliation(s)
- Hongshuai Chu
- Ministry of Education Key Laboratory for Earth System Modeling, Department of Earth System Science, Tsinghua University, Beijing 100084, China
| | - Sergey Venevsky
- Ministry of Education Key Laboratory for Earth System Modeling, Department of Earth System Science, Tsinghua University, Beijing 100084, China.
| | - Chao Wu
- Ministry of Education Key Laboratory for Earth System Modeling, Department of Earth System Science, Tsinghua University, Beijing 100084, China; College of Life and Environment Sciences, University of Exeter, Exeter EX4 4QF, UK
| | - Menghui Wang
- Ministry of Education Key Laboratory for Earth System Modeling, Department of Earth System Science, Tsinghua University, Beijing 100084, China
| |
Collapse
|
12
|
|