Parameter identification, model calibration, and uncertainty quantification are important steps in the model-building process, and are necessary for obtaining credible results and valuable information. Sensitivity ana...Parameter identification, model calibration, and uncertainty quantification are important steps in the model-building process, and are necessary for obtaining credible results and valuable information. Sensitivity analysis of hydrological model is a key step in model uncertainty quantification, which can identify the dominant parameters, reduce the model calibration uncertainty, and enhance the model optimization efficiency. There are, however, some shortcomings in classical approaches, including the long duration of time and high computation cost required to quantitatively assess the sensitivity of a multiple-parameter hydrological model. For this reason, a two-step statistical evaluation framework using global techniques is presented. It is based on （1） a screening method （Morris） for qualitative ranking of parameters, and （2） a variance-based method integrated with a meta-model for quantitative sensitivity analysis, i.e., the Sobol method integrated with the response surface model （RSMSobol）. First, the Morris screening method was used to qualitatively identify the parameters＇ sensitivity, and then ten parameters were selected to quantify the sensitivity indices. Subsequently, the RSMSobol method was used to quantify the sensitivity, i.e., the first-order and total sensitivity indices based on the response surface model （RSM） were calculated. The RSMSobol method can not only quantify the sensitivity, but also reduce the computational cost, with good accuracy compared to the classical approaches. This approach will be effective and reliable in the global sensitivity analysis of a complex large-scale distributed hydrological model.展开更多
The North East China Transect (NECT) is an exemplary region of landmass running along the line of 43°30′ North Latitude and caught between 112° and 130°30′ East Longtitude, and is being sampled by the...The North East China Transect (NECT) is an exemplary region of landmass running along the line of 43°30′ North Latitude and caught between 112° and 130°30′ East Longtitude, and is being sampled by the International Geosphere-Biosphere Program (IGBP) for today’s study on global change. So far, it has already been listed in the first set of IGBP-sponsored transects, becoming a key theatre and "hot spot" for probing terrestrial ecosystems. In terms of biota, the nearly 1,600-kilometre-long tract is located in a mid-latitude zone, featuring a vegetational transition from a temperate timberland of evergreen coniferous forests and broadleaved deciduous forests to a mild-temperate steppe. According to the norms of biome types, it consists of three subzones: meadow steppe, typical steppe and desert steppe in a continuous east-to-west spatial sequence. There are four ecological stations supported by a great number of permanent samples, long-term plots and an enormous build-up of experimental data along展开更多
The article presents a methodology for transitioning to the Amsterdam height system based on utilizing global and regional geoid/quasi-geoid models.The study was conducted for the Polish-Ukrainian cross-border sector ...The article presents a methodology for transitioning to the Amsterdam height system based on utilizing global and regional geoid/quasi-geoid models.The study was conducted for the Polish-Ukrainian cross-border sector and expanded to the entire territory of Poland and Ukraine.The input data comprised two regional and five global geoid/quasi-geoid models.The initial data analysis was conducted for all models relative to GNSS/leveling data in the Baltic height system.The secondary analysis was performed relative to the combined PL-quasi-geoid2021 model and the gravimetric EGG2015 model.Based on the analysis results,a methodology for optimizing heights between regional and global geoid/quasi-geoid models was developed,including the following stages:calculation of conditional global and regional geoid/quasi-geoid heights,calculation of approximately predicted height differences between the conditional regional and global geoid/quasi-geoid,implementation of the refinement(correction)of approximate heights at the regional model,optimization of approximate heights at the regional model,and calculation of a regional combined model in the Amsterdam height system.The developed methodology enables the integration of regional and global geoid/quasi-geoid models into the Amsterdam height system with an accuracy of 1-2 cm by optimizing their heights.The advantage of this methodology is that it requires only a minimal amount of GNSS/leveling data to establish connections between different height systems.展开更多
The kinetic characteristics of plasma-assisted oxidative pyrolysis of ammonia are studied by using the global/fluid models hybrid solution method.Firstly,the stable products of plasma-assisted oxidative pyrolysis of a...The kinetic characteristics of plasma-assisted oxidative pyrolysis of ammonia are studied by using the global/fluid models hybrid solution method.Firstly,the stable products of plasma-assisted oxidative pyrolysis of ammonia are measured.The results show that the consumption of NH_(3)/O_(2)and the production of N_(2)/H_(2)change linearly with the increase of voltage,which indicates the decoupling of nonequilibrium molecular excitation and oxidative pyrolysis of ammonia at low temperatures.Secondly,the detailed reaction kinetics mechanism of ammonia oxidative pyrolysis stimulated by a nanosecond pulse voltage at low pressure and room temperature is established.Based on the reaction path analysis,the simplified mechanism is obtained.The detailed and simplified mechanism simulation results are compared with experimental data to verify the accuracy of the simplified mechanism.Finally,based on the simplified mechanism,the fluid model of ammonia oxidative pyrolysis stimulated by the nanosecond pulse plasma is established to study the pre-sheath/sheath behavior and the resultant consumption and formation of key species.The results show that the generation,development,and propagation of the pre-sheath have a great influence on the formation and consumption of species.The consumption of NH_(3)by the cathode pre-sheath is greater than that by the anode pre-sheath,but the opposite is true for OH and O(1S).However,within the sheath,almost all reactions do not occur.Further,by changing the parameters of nanosecond pulse power supply voltage,it is found that the electron number density,electron current density,and applied peak voltages are not the direct reasons for the structural changes of the sheath and pre-sheath.Furthermore,the discharge interval has little effect on the sheath structure and gas mixture breakdown.The research results of this paper not only help to understand the kinetic promotion of non-equilibrium excitation in the process of oxidative pyrolysis but also help to explore the influence of transport and chemical reaction kinetics on the oxidative pyrolysis of ammonia.展开更多
Radio frequency capacitively coupled plasmas(RF CCPs)operated in Ar/O_(2)gas mixtures which are widely adopted in microelectronics,display,and photovoltaic industry,are investigated based on an equivalent circuit mode...Radio frequency capacitively coupled plasmas(RF CCPs)operated in Ar/O_(2)gas mixtures which are widely adopted in microelectronics,display,and photovoltaic industry,are investigated based on an equivalent circuit model coupled with a global model.This study focuses on the effects of singlet metastable molecule O_(2)(b^(1)∑_(8)^(+)),highly excited Herzberg states O_(2)(A^(3)∑_(u)^(+),A^(3)△_(u),c^(1)∑_(u)^(-)),and the negative ion O_(2)^(-),which are usually neglected in simulation studies.Specifically,their impact on particle densities,electronegativity,electron temperature,voltage drop across the sheath,and absorbed power in the discharge is analyzed.The results indicate that O_(2)(b^(1)∑_(8)^(+))and O_(2)^(-)exhibit relatively high densities in argon-oxygen discharges.While O_(2)(A^(3)∑_(u)^(+),A^(3)△_(u),c^(1)∑_(u)^(-))play a critical role in O_(2)b1S+g production,especially at higher pressure.The inclusion of these particles reduces the electronegativity,electron temperature,and key species densities,especially the O^(-)and O^(*)densities.Moreover,the sheath voltage drop,as well as the inductance and resistance of the plasma bulk are enhanced,while the sheath dissipation power and total absorbed power decrease slightly.With the increasing pressure,the influence of these particles on the discharge properties becomes more significant.The study also explores the generation and loss of main neutral species and charged particles within the pressure range of 20 mTorr-100 mTorr(1 Torr=1.33322×10^(2)Pa),offering insights into essential and non-essential reactions for future low-pressure O_(2)and Ar/O_(2)CCP discharge modeling.展开更多
Accurately simulating mesoscale convective systems(MCSs)is essential for predicting global precipitation patterns and extreme weather events.Despite the ability of advanced models to reproduce MCS climate statistics,c...Accurately simulating mesoscale convective systems(MCSs)is essential for predicting global precipitation patterns and extreme weather events.Despite the ability of advanced models to reproduce MCS climate statistics,capturing extreme storm cases over complex terrain remains challenging.This study utilizes the Global–Regional Integrated Forecast System(GRIST)with variable resolution to simulate an eastward-propagating MCS event.The impact of three microphysics schemes,including two single-moment schemes(WSM6,Lin)and one double-moment scheme(Morrison),on the model sensitivity of MCS precipitation simulations is investigated.The results demonstrate that while all the schemes capture the spatial distribution and temporal variation of MCS precipitation,the Morrison scheme alleviates overestimated precipitation compared to the Lin and WSM6 schemes.The ascending motion gradually becomes weaker in the Morrison scheme during the MCS movement process.Compared to the runs with convection parameterization,the explicit-convection setup at 3.5-km resolution reduces disparities in atmospheric dynamics due to microphysics sensitivity in terms of vertical motions and horizontal kinetic energy at the high-wavenumber regimes.The explicit-convection setup more accurately captures the propagation of both main and secondary precipitation centers during the MCS development,diminishing the differences in both precipitation intensity and propagation features between the Morrison and two single-moment schemes.These findings underscore the importance of microphysics schemes for global nonhydrostatic modeling at the kilometer scale.The role of explicit convection for reducing model uncertainty is also outlined.展开更多
Glacier mass balance is a key indicator of glacier health and climate change sensitivity.Influencing factors include both climatic and nonclimatic elements,forming a complex set of drivers.There is a lack of quantitat...Glacier mass balance is a key indicator of glacier health and climate change sensitivity.Influencing factors include both climatic and nonclimatic elements,forming a complex set of drivers.There is a lack of quantitative analysis of these composite factors,particularly in climate-typical regions like the Tanggula Mountains on the central Tibetan Plateau.We collected data on various factors affecting glacier mass balance from 2000 to 2020,including climate variables,topographic variables,geometric parameters,and glacier dynamics.We utilized linear regression models,ensemble learning models,and Open Global Glacier Model(OGGM)to analyze glacier mass balance changes in the Tanggula Mountains.Results indicate that linear models explain 58%of the variance in glacier mass balance,with seasonal temperature and precipitation having significant impacts.Our findings show that ensemble learning models made the explanations 5.2%more accurate by including the impact of topographic and geometric factors such as the average glacier height,the slope of the glacier tongue,the speed of the ice flow,and the area of the glacier.Interpretable machine learning identified the spatial distribution of positive and negative impacts of these characteristics and the interaction between glacier topography and ice dynamics.Finally,we predicted the responses of glaciers of different sizes to future climate change based on the results of interpretable machine learning.It was found that relatively large glaciers(>1 km~2)are likely to persist until the end of this century under low emission scenarios,whereas small glaciers(<1 km~2)are expected to nearly disappear by 2080 under any emission scenario.Our research provides technical support for improving glacier change modeling and protection on the Tibetan Plateau.展开更多
Motivated by the wide applications of empirical global ocean tide models in Earth science,particularly in shallow waters and coastal regions,we proposed an updated global ocean tide model representing all major diurna...Motivated by the wide applications of empirical global ocean tide models in Earth science,particularly in shallow waters and coastal regions,we proposed an updated global ocean tide model representing all major diurnal and semidiurnal tidal constituents.We integrated the development technique of the Technical University of Denmark(DTU10)model and calculated the residual tides using the intermediary FES2012 model.We utilized all available Topex/Poseidon,Jason-1,Ocean Surface Topography Mission(OSTM)/Jason-2 primary and tandem missions to develop the new model.To reduce the discrepancies between the model and in situ measurements estimated S2 tide constituent,the ERA-Interim model was selected for dry tropospheric correction of altimetry data.The newly developed model with an improved spatial resolution from 1/8°×1/8°to 1/16°×1/16°was assessed and compared against contemporary global ocean tide models using in situ measurements in coastal regions,continental shelf waters,marginal seas,and deep waters.The results demonstrate improvements in shallow waters and coastal regions,particularly in representative coastal re-gions such as the Northwest European Shelf and East Asian Marginal Seas.Moreover,the model performs well in simulating S_(2) and K_(2) tides with root mean square differences of 0.355 and 0.234 cm,respectively,against in situ measurements in deep waters world-wide.Compared with that of DTU10,the root sum square of the new model for the eight main tidal constituents decrease by 8.4%to 0.997 cm.展开更多
ABSTRACT The lAP Dynamic Global Vegetation Model （IAP-DGVM） has been developed to simulate the distribution and structure of global vegetation within the framework of Earth System Models. It incorporates our group...ABSTRACT The lAP Dynamic Global Vegetation Model （IAP-DGVM） has been developed to simulate the distribution and structure of global vegetation within the framework of Earth System Models. It incorporates our group＇s recent developments of major model components such as the shrub sub-model, establishment and competition parameterization schemes, and a process-based fire parameterization of intermediate complexity. The model has 12 plant functional types, including seven tree, two shrub, and three grass types, plus bare soil. Different PFTs are allowed to coexist within a grid cell, and their state variables are updated by various governing equations describing vegetation processes from fine-scale biogeophysics and biogeochemistry, to individual and population dynamics, to large-scale biogeography. Environmental disturbance due to fire not only affects regional vegetation competition, but also influences atmospheric chemistry and aerosol emissions. Simulations under observed atmospheric conditions showed that the model can correctly reproduce the global distribution of trees, shrubs, grasses, and bare soil. The simulated global dominant vegetation types reproduce the transition from forest to grassland （savanna） in the tropical region, and from forest to shrubland in the boreal region, but overestimate the region of temperate forest.展开更多
The capability of an improved Dynamic Global Vegetation Model （DGVM） in reproducing the impact of climate on the terrestrial ecosystem is evaluated. The new model incorporates the Community Land Model- DGVM （CLM3.0...The capability of an improved Dynamic Global Vegetation Model （DGVM） in reproducing the impact of climate on the terrestrial ecosystem is evaluated. The new model incorporates the Community Land Model- DGVM （CLM3.0-DGVM） with a submodel for temperate and boreal shrubs, as well as other revisions such as the ＂two-leaf＂ scheme for photosynthesis and the definition of fractional coverage of plant functional types （PFTs）. Results show that the revised model may correctly reproduce the global distribution of temperate and boreal shrubs, and improves the model performance with more realistic distribution of di？erent vege- tation types. The revised model also correctly reproduces the zonal distributions of vegetation types. In reproducing the dependence of the vegetation distribution on climate conditions, the model shows that the dominant regions for trees, grasses, shrubs, and bare soil are clearly separated by a climate index derived from mean annual precipitation and temperature, in good agreement with the CLM4 surface data. The dominant plant functional type mapping to a two dimensional parameter space of mean annual temperature and precipitation also qualitatively agrees with the results from observations and theoretical ecology studies.展开更多
Heilongjiang province is the largest forest zone in China and the forest coverage rate is 46%. Forests of Heilongjiang province play an important role in the forest ecosystem of China. In this study we investi- gated ...Heilongjiang province is the largest forest zone in China and the forest coverage rate is 46%. Forests of Heilongjiang province play an important role in the forest ecosystem of China. In this study we investi- gated the spatial distribution of forest carbon storage in Heilongjiang province using 3083 plots sampled in 2010. We attempted to fit two global models, ordinary least squares model （OLS）, linear mixed model （LMM）, and a local model, geographically weighted regression model （GWR）, to the relationship between forest carbon content and stand, environment, and climate factors. Five predictors significantly affected forest carbon storage and spatial distribution, viz. average diameter of stand （DBH）, number of trees per hectare （TPH）, elevation （Elev）, slope （Slope） and the product of precipitation and temperature （Rain Temp）. The GWR model outperformed the two global models in both model fitting and prediction because it successfully reduced both spatial auto- correlation and heterogeneity in model residuals. More importantly, the GWR model provided localized model coefficients for each location in the study area, which allowed us to evaluate the influences of local stand conditions and topographic features on tree and stand growth, and forest carbon stock. It also helped us to better understand the impacts of silvi- cultural and management activities on the amount and changes of forest carbon storage across the province. The detailed information can be readily incorporated with the mapping ability of GIS software to provide excellent tools for assessing the distribution and dynamics of the for- est-carbon stock in the next few years.展开更多
With the completion of Chinese BeiDou Navigation Satellite System(BDS),the world has begun to enjoy the Positioning,Navigation,and Timing(PNT)services of four Global Navigation Satellite Systems(GNSS).In order to impr...With the completion of Chinese BeiDou Navigation Satellite System(BDS),the world has begun to enjoy the Positioning,Navigation,and Timing(PNT)services of four Global Navigation Satellite Systems(GNSS).In order to improve the GNSS performance and expand its applications,Low Earth Orbit(LEO)Enhanced Global Navigation Satellite System(LeGNSS)is being vigorously advocated.Combined with high-,medium-,and low-earth orbit satellites,it can improve GNSS performance in terms of orbit determination,Precise Point Positioning(PPP)convergence time,etc.This paper comprehensively reviews the current status of LeGNSS,focusing on analyzing its advantages and challenges for precise orbit and clock determination,PPP convergence,earth rotation parameter estimation,and global ionosphere modeling.Thanks to the fast geometric change brought by LEO satellites,LeGNSS is expected to fundamentally solve the problem of the long convergence time of PPP without any augmentation.The convergence time can be shortened within 1 minute if appropriate LEO constellations are deployed.However,there are still some issues to overcome,such as the optimization of LEO constellation as well as the real time LEO precise orbit and clock determination.展开更多
The interest in the development and improvement of dynamic global vegetation models （DGVMs）, which have the potential to simulate fluxes of carbon, water and nitrogen, along with changes in the vegetation dynamics, ...The interest in the development and improvement of dynamic global vegetation models （DGVMs）, which have the potential to simulate fluxes of carbon, water and nitrogen, along with changes in the vegetation dynamics, within an integrated system, has been increasing. In this paper, some numerical schemes and a higher resolution soil texture dataset were employed to improve the Sheffield Dynamic Global Vegetation Model （SDGVM）. Using eddy covariance-based measurements, we then tested the standard version of the SDGVM and the modified version of the SDGVM. Detailed observations of daily carbon and water fluxes made at the upland oak forest on the Walker Branch Watershed in Tennessee, USA offered a unique opportunity for these comparisons. The results revealed that the modified version of the SDGVM did a reasonable job of simulating the carbon and water flux and the variation of soil water content （SWC）. However, at the end of the growing season, it failed to simulate the effect of the limitations on the soil respiration dynamics and as a result underestimated this respiration. It was also noted that the modified version overestimated the increase in the SWC following summer rainfall, which was attributed to an inadequate representation of the ground water and thermal cycle.展开更多
In the past several decades, dynamic global vegetation models（DGVMs） have been the most widely used and appropriate tool at the global scale to investigate vegetation-climate interactions. At the Institute of Atmosp...In the past several decades, dynamic global vegetation models（DGVMs） have been the most widely used and appropriate tool at the global scale to investigate vegetation-climate interactions. At the Institute of Atmospheric Physics, a new version of DGVM（IAP-DGVM） has been developed and coupled to the Common Land Model（CoLM） within the framework of the Chinese Academy of Sciences＇ Earth System Model（CAS-ESM）. This work reports the performance of IAP-DGVM through comparisons with that of the default DGVM of CoLM（CoLM-DGVM） and observations. With respect to CoLMDGVM, IAP-DGVM simulated fewer tropical trees, more ＂needleleaf evergreen boreal tree＂ and ＂broadleaf deciduous boreal shrub＂, and a better representation of grasses. These contributed to a more realistic vegetation distribution in IAP-DGVM,including spatial patterns, total areas, and compositions. Moreover, IAP-DGVM also produced more accurate carbon fluxes than CoLM-DGVM when compared with observational estimates. Gross primary productivity and net primary production in IAP-DGVM were in better agreement with observations than those of CoLM-DGVM, and the tropical pattern of fire carbon emissions in IAP-DGVM was much more consistent with the observation than that in CoLM-DGVM. The leaf area index simulated by IAP-DGVM was closer to the observation than that of CoLM-DGVM; however, both simulated values about twice as large as in the observation. This evaluation provides valuable information for the application of CAS-ESM, as well as for other model communities in terms of a comparative benchmark.展开更多
A dynamic global vegetation model （DGVM） coupled with a land surface model （LSM） is generally initialized using a spin-up process to derive a physically-consistent initial condition. Spin-up forcing, which is the ...A dynamic global vegetation model （DGVM） coupled with a land surface model （LSM） is generally initialized using a spin-up process to derive a physically-consistent initial condition. Spin-up forcing, which is the atmospheric forcing used to drive the coupled model to equilibrium solutions in the spin-up process, varies across earlier studies. In the present study, the impact of the spin-up forcing in the initialization stage on the fractional coverages （FCs） of plant functional type （PFT） in the subsequent simulation stage are assessed in seven classic climate regions by a modified Community Land Model’s Dynamic Global Vegetation Model （CLM-DGVM）. Results show that the impact of spin-up forcing is considerable in all regions except the tropical rainforest climate region （TR） and the wet temperate climate region （WM）. In the tropical monsoon climate region （TM）, the TR and TM transition region （TR-TM）, the dry temperate climate region （DM）, the highland climate region （H）, and the boreal forest climate region （BF）, where FCs are affected by climate non-negligibly, the discrepancies in initial FCs, which represent long-term cumulative response of vegetation to different climate anomalies, are large. Moreover, the large discrepancies in initial FCs usually decay slowly because there are trees or shrubs in the five regions. The intrinsic growth timescales of FCs for tree PFTs and shrub PFTs are long, and the variation of FCs of tree PFTs or shrub PFTs can affect that of grass PFTs.展开更多
In Dynamic Global Vegetation Models （DGVMs）, the establishment of woody vegetation refers to flowering, fertiliza- tion, seed production, germination, and the growth of tree seedlings. It determines not only the pop...In Dynamic Global Vegetation Models （DGVMs）, the establishment of woody vegetation refers to flowering, fertiliza- tion, seed production, germination, and the growth of tree seedlings. It determines not only the population densities but also other important ecosystem structural variables. In current DGVMs, establishments of woody plant functional types （PFTs） are assumed to be either the same in the same grid cell, or largely stochastic. We investigated the uncertainties in the competition of establishment among coexisting woody PFTs from three aspects： the dependence of PFT establishments on vegetation states; background establishment; and relative establishment potentials of different PFTs. Sensitivity experi- ments showed that the dependence of establishment rate on the fractional coverage of a PFT favored the dominant PFT by increasing its share in establishment. While a small background establishment rate had little impact on equilibrium states of the ecosystem, it did change the timescale required for the establishment of alien species in pre-existing forest due to their disadvantage in seed competition during the early stage of invasion. Meanwhile, establishment purely fiom background （the scheme commonly used in current DGVMs） led to inconsistent behavior in response to the change in PFT specification （e.g., number of PFTs and their specification）. Furthermore, the results also indicated that trade-off between irtdividual growth and reproduction/colonization has significant influences on the competition of establishment. Hence, further development of es- tablishment parameterization in DGVMs is essential in reducing the uncertainties in simulations of both ecosystem structures and successions.展开更多
To investigate the interaction between the tropical Pacific and China seas a variable-grid global ocean circulation model with fine grid covering the area from 20°S to 50°N and from 99° to 150°E is...To investigate the interaction between the tropical Pacific and China seas a variable-grid global ocean circulation model with fine grid covering the area from 20°S to 50°N and from 99° to 150°E is developed. Numerical computation of the annually cyclic circulation fields is performed. The results of the annual mean zonal currents and deep to abyssal western boundary currents in the equatorial Pacific Ocean are reported. The North Equatorial Current,the North Equatorial Countercurrent, the South Equatorial Current and the Equatorial Undercurrent are fairly well simulated. The model well reproduces the northward flowing abyssal western boundary current.From the model results a lower deep western boundary current east of the Bismarck-Solomon-New Hebrides Island chain at depths around 2 000 m has been found. The model results also show that the currents in the equatorial Pacific Ocean have multi-layer structures both in zonal currents and western boundary currents, indicating that the global ocean overturning thermohaline circulation appears of multi-layer pattern.展开更多
On the basis of the newly released International Terrestrial Reference Frame（ITRF2008） by the International Earth Rotation Service （IERS）, a new global plate model ITRF2008 plate for the major plates is establishe...On the basis of the newly released International Terrestrial Reference Frame（ITRF2008） by the International Earth Rotation Service （IERS）, a new global plate model ITRF2008 plate for the major plates is established. This ITRF2008-derived model is analyzed in comparison with NNR-NUVEL1A model, which is mainly based on geological and geophysical data. The Eurasia and Paeifi6 plates display obvious differences in terms of the velocity fields derived from the two plate motion models. Plate acceleration is also introduced to characterize the differences of the two velocity fields which obtained from ITRF2008-plate and NNR-NUVEL1A models for major individual plates. The results show that the Africa, South America and Eurasia plates are undergoing acceleration, while the North America and Australia plates are in the state of deceleration motion展开更多
Using wave measurements from the EMFISIS instrument onboard Van Allen Probes,we investigate statistically the spatial distributions of the intensity of plasmaspheric hiss waves.To reproduce these empirical results,we ...Using wave measurements from the EMFISIS instrument onboard Van Allen Probes,we investigate statistically the spatial distributions of the intensity of plasmaspheric hiss waves.To reproduce these empirical results,we establish a fitting model that is a thirdorder polynomial function of L-shell,magnetic local time(MLT),magnetic latitude(MLAT),and AE*.Quantitative comparisons indicate that the model’s fitting functions can reflect favorably the major empirical features of the global distribution of hiss wave intensity,including substorm dependence and the MLT asymmetry.Our results therefore provide a useful analytic model that can be readily employed in future simulations of global radiation belt electron dynamics under the impact of plasmaspheric hiss waves in geospace.展开更多
Tropical mountainous areas not only provide substantial carbon storage and play an important role in global biological diversity, but also provide basic livelihood for a large number of poor ethnic minorities. However...Tropical mountainous areas not only provide substantial carbon storage and play an important role in global biological diversity, but also provide basic livelihood for a large number of poor ethnic minorities. However, there is no unified and explicit definition for mountainous areas. The local elevation range(LER) is a crucial structural parameter for delineating mountainous areas. However, current LER products are limited by the subjective selection of an optimum statistical window or coarser spatial resolution of topographical data. In this study, we presented an approach using thresholds for three topographic parameters, elevation, slope, and LER, derived from the Advanced Spaceborne Thermal Emission and Reflection Radiometer Global Digital Elevation Model(ASTER GDEM) to redelineate the vast mountainous areas of mainland Southeast Asia(MSEA). The mean change-point analysis method was applied to determine the optimum statistical window of the 1 arc second(approximately 30 m)-resolution GDEM LER. The results showed that: First, the optimum statistical window is 38 × 38 cell units(width × height) in a rectangular neighborhood, or an area of about 1.30 km^2 for calculating GDEM LER in MSEA. Second, the LER of more than 80% of the area ranges from 30 m to 499 m in MSEA. The LERs in the northern and northwestern MSEA are greater than their counterparts in the south and east. Third, the area of the re-delineated mountainous areas was 83.52 × 10~4 km^2, about 38.71% of the total area. Spatially, the mountainous areas are mainly distributed in the north and northeast of MSEA. The re-delineated 30-m resolution map of the mountainous areas will serve as a topographical dataset for monitoring mountainrelated land surface changes in MSEA. The parameter-modified mountain extraction procedure can be expanded to delineate global mountainous areas.展开更多
基金supported by the National Natural Science Foundation of China (Grant No. 41271003)the National Basic Research Program of China (Grants No. 2010CB428403 and 2010CB951103)
文摘Parameter identification, model calibration, and uncertainty quantification are important steps in the model-building process, and are necessary for obtaining credible results and valuable information. Sensitivity analysis of hydrological model is a key step in model uncertainty quantification, which can identify the dominant parameters, reduce the model calibration uncertainty, and enhance the model optimization efficiency. There are, however, some shortcomings in classical approaches, including the long duration of time and high computation cost required to quantitatively assess the sensitivity of a multiple-parameter hydrological model. For this reason, a two-step statistical evaluation framework using global techniques is presented. It is based on （1） a screening method （Morris） for qualitative ranking of parameters, and （2） a variance-based method integrated with a meta-model for quantitative sensitivity analysis, i.e., the Sobol method integrated with the response surface model （RSMSobol）. First, the Morris screening method was used to qualitatively identify the parameters＇ sensitivity, and then ten parameters were selected to quantify the sensitivity indices. Subsequently, the RSMSobol method was used to quantify the sensitivity, i.e., the first-order and total sensitivity indices based on the response surface model （RSM） were calculated. The RSMSobol method can not only quantify the sensitivity, but also reduce the computational cost, with good accuracy compared to the classical approaches. This approach will be effective and reliable in the global sensitivity analysis of a complex large-scale distributed hydrological model.
基金This research work is jointly subsidied by two research projects entitled"A Predictive Study on the Changing Trend of Life-supporting Environment in China Over the Next 20-50 YearsA Modeling Research on the Responses of China's Terrestrial Ecosys
文摘The North East China Transect (NECT) is an exemplary region of landmass running along the line of 43°30′ North Latitude and caught between 112° and 130°30′ East Longtitude, and is being sampled by the International Geosphere-Biosphere Program (IGBP) for today’s study on global change. So far, it has already been listed in the first set of IGBP-sponsored transects, becoming a key theatre and "hot spot" for probing terrestrial ecosystems. In terms of biota, the nearly 1,600-kilometre-long tract is located in a mid-latitude zone, featuring a vegetational transition from a temperate timberland of evergreen coniferous forests and broadleaved deciduous forests to a mild-temperate steppe. According to the norms of biome types, it consists of three subzones: meadow steppe, typical steppe and desert steppe in a continuous east-to-west spatial sequence. There are four ecological stations supported by a great number of permanent samples, long-term plots and an enormous build-up of experimental data along
文摘The article presents a methodology for transitioning to the Amsterdam height system based on utilizing global and regional geoid/quasi-geoid models.The study was conducted for the Polish-Ukrainian cross-border sector and expanded to the entire territory of Poland and Ukraine.The input data comprised two regional and five global geoid/quasi-geoid models.The initial data analysis was conducted for all models relative to GNSS/leveling data in the Baltic height system.The secondary analysis was performed relative to the combined PL-quasi-geoid2021 model and the gravimetric EGG2015 model.Based on the analysis results,a methodology for optimizing heights between regional and global geoid/quasi-geoid models was developed,including the following stages:calculation of conditional global and regional geoid/quasi-geoid heights,calculation of approximately predicted height differences between the conditional regional and global geoid/quasi-geoid,implementation of the refinement(correction)of approximate heights at the regional model,optimization of approximate heights at the regional model,and calculation of a regional combined model in the Amsterdam height system.The developed methodology enables the integration of regional and global geoid/quasi-geoid models into the Amsterdam height system with an accuracy of 1-2 cm by optimizing their heights.The advantage of this methodology is that it requires only a minimal amount of GNSS/leveling data to establish connections between different height systems.
基金Fundamental Research Funds for the Central Universities(M23JBZY00050)National Natural Science Foundation of China(22278032)。
文摘The kinetic characteristics of plasma-assisted oxidative pyrolysis of ammonia are studied by using the global/fluid models hybrid solution method.Firstly,the stable products of plasma-assisted oxidative pyrolysis of ammonia are measured.The results show that the consumption of NH_(3)/O_(2)and the production of N_(2)/H_(2)change linearly with the increase of voltage,which indicates the decoupling of nonequilibrium molecular excitation and oxidative pyrolysis of ammonia at low temperatures.Secondly,the detailed reaction kinetics mechanism of ammonia oxidative pyrolysis stimulated by a nanosecond pulse voltage at low pressure and room temperature is established.Based on the reaction path analysis,the simplified mechanism is obtained.The detailed and simplified mechanism simulation results are compared with experimental data to verify the accuracy of the simplified mechanism.Finally,based on the simplified mechanism,the fluid model of ammonia oxidative pyrolysis stimulated by the nanosecond pulse plasma is established to study the pre-sheath/sheath behavior and the resultant consumption and formation of key species.The results show that the generation,development,and propagation of the pre-sheath have a great influence on the formation and consumption of species.The consumption of NH_(3)by the cathode pre-sheath is greater than that by the anode pre-sheath,but the opposite is true for OH and O(1S).However,within the sheath,almost all reactions do not occur.Further,by changing the parameters of nanosecond pulse power supply voltage,it is found that the electron number density,electron current density,and applied peak voltages are not the direct reasons for the structural changes of the sheath and pre-sheath.Furthermore,the discharge interval has little effect on the sheath structure and gas mixture breakdown.The research results of this paper not only help to understand the kinetic promotion of non-equilibrium excitation in the process of oxidative pyrolysis but also help to explore the influence of transport and chemical reaction kinetics on the oxidative pyrolysis of ammonia.
基金supported by the National Natural Science Foundation of China(Grant Nos.12020101005,12475202,12347131,and 12405289).
文摘Radio frequency capacitively coupled plasmas(RF CCPs)operated in Ar/O_(2)gas mixtures which are widely adopted in microelectronics,display,and photovoltaic industry,are investigated based on an equivalent circuit model coupled with a global model.This study focuses on the effects of singlet metastable molecule O_(2)(b^(1)∑_(8)^(+)),highly excited Herzberg states O_(2)(A^(3)∑_(u)^(+),A^(3)△_(u),c^(1)∑_(u)^(-)),and the negative ion O_(2)^(-),which are usually neglected in simulation studies.Specifically,their impact on particle densities,electronegativity,electron temperature,voltage drop across the sheath,and absorbed power in the discharge is analyzed.The results indicate that O_(2)(b^(1)∑_(8)^(+))and O_(2)^(-)exhibit relatively high densities in argon-oxygen discharges.While O_(2)(A^(3)∑_(u)^(+),A^(3)△_(u),c^(1)∑_(u)^(-))play a critical role in O_(2)b1S+g production,especially at higher pressure.The inclusion of these particles reduces the electronegativity,electron temperature,and key species densities,especially the O^(-)and O^(*)densities.Moreover,the sheath voltage drop,as well as the inductance and resistance of the plasma bulk are enhanced,while the sheath dissipation power and total absorbed power decrease slightly.With the increasing pressure,the influence of these particles on the discharge properties becomes more significant.The study also explores the generation and loss of main neutral species and charged particles within the pressure range of 20 mTorr-100 mTorr(1 Torr=1.33322×10^(2)Pa),offering insights into essential and non-essential reactions for future low-pressure O_(2)and Ar/O_(2)CCP discharge modeling.
基金supported by the National Natural Science Foundation of China(Grant No.42305169)the Basic Research Fund of CAMS(Grant No.2023Y001)the National Key Scientific and Technological Infrastructure project“Earth System Numerical Simulation Facility”(Earth Lab)。
文摘Accurately simulating mesoscale convective systems(MCSs)is essential for predicting global precipitation patterns and extreme weather events.Despite the ability of advanced models to reproduce MCS climate statistics,capturing extreme storm cases over complex terrain remains challenging.This study utilizes the Global–Regional Integrated Forecast System(GRIST)with variable resolution to simulate an eastward-propagating MCS event.The impact of three microphysics schemes,including two single-moment schemes(WSM6,Lin)and one double-moment scheme(Morrison),on the model sensitivity of MCS precipitation simulations is investigated.The results demonstrate that while all the schemes capture the spatial distribution and temporal variation of MCS precipitation,the Morrison scheme alleviates overestimated precipitation compared to the Lin and WSM6 schemes.The ascending motion gradually becomes weaker in the Morrison scheme during the MCS movement process.Compared to the runs with convection parameterization,the explicit-convection setup at 3.5-km resolution reduces disparities in atmospheric dynamics due to microphysics sensitivity in terms of vertical motions and horizontal kinetic energy at the high-wavenumber regimes.The explicit-convection setup more accurately captures the propagation of both main and secondary precipitation centers during the MCS development,diminishing the differences in both precipitation intensity and propagation features between the Morrison and two single-moment schemes.These findings underscore the importance of microphysics schemes for global nonhydrostatic modeling at the kilometer scale.The role of explicit convection for reducing model uncertainty is also outlined.
基金funding from the National Key Research and Development Program of China(2023YFC3206300)the Gansu Provincial Science and Technology Program(22ZD6FA005)+2 种基金the Gansu Youth Science and Technology Fund(E4310103)the Gansu Postdoctoral Science Foundation(E339880112)the Tibet Science and Technology Program(XZ202301ZY0001G and XZ202401JD0007)。
文摘Glacier mass balance is a key indicator of glacier health and climate change sensitivity.Influencing factors include both climatic and nonclimatic elements,forming a complex set of drivers.There is a lack of quantitative analysis of these composite factors,particularly in climate-typical regions like the Tanggula Mountains on the central Tibetan Plateau.We collected data on various factors affecting glacier mass balance from 2000 to 2020,including climate variables,topographic variables,geometric parameters,and glacier dynamics.We utilized linear regression models,ensemble learning models,and Open Global Glacier Model(OGGM)to analyze glacier mass balance changes in the Tanggula Mountains.Results indicate that linear models explain 58%of the variance in glacier mass balance,with seasonal temperature and precipitation having significant impacts.Our findings show that ensemble learning models made the explanations 5.2%more accurate by including the impact of topographic and geometric factors such as the average glacier height,the slope of the glacier tongue,the speed of the ice flow,and the area of the glacier.Interpretable machine learning identified the spatial distribution of positive and negative impacts of these characteristics and the interaction between glacier topography and ice dynamics.Finally,we predicted the responses of glaciers of different sizes to future climate change based on the results of interpretable machine learning.It was found that relatively large glaciers(>1 km~2)are likely to persist until the end of this century under low emission scenarios,whereas small glaciers(<1 km~2)are expected to nearly disappear by 2080 under any emission scenario.Our research provides technical support for improving glacier change modeling and protection on the Tibetan Plateau.
基金supported by the open funding of the Technology Innovation Center for South China Sea Re-mote Sensing,Surveying and Mapping Collaborative Ap-plication,Ministry of Natural Resources,P.R.China(No.RSSMCA-2024-B001)the National Natural Science Foundation of China(Nos.T2261149752 and 42476172).
文摘Motivated by the wide applications of empirical global ocean tide models in Earth science,particularly in shallow waters and coastal regions,we proposed an updated global ocean tide model representing all major diurnal and semidiurnal tidal constituents.We integrated the development technique of the Technical University of Denmark(DTU10)model and calculated the residual tides using the intermediary FES2012 model.We utilized all available Topex/Poseidon,Jason-1,Ocean Surface Topography Mission(OSTM)/Jason-2 primary and tandem missions to develop the new model.To reduce the discrepancies between the model and in situ measurements estimated S2 tide constituent,the ERA-Interim model was selected for dry tropospheric correction of altimetry data.The newly developed model with an improved spatial resolution from 1/8°×1/8°to 1/16°×1/16°was assessed and compared against contemporary global ocean tide models using in situ measurements in coastal regions,continental shelf waters,marginal seas,and deep waters.The results demonstrate improvements in shallow waters and coastal regions,particularly in representative coastal re-gions such as the Northwest European Shelf and East Asian Marginal Seas.Moreover,the model performs well in simulating S_(2) and K_(2) tides with root mean square differences of 0.355 and 0.234 cm,respectively,against in situ measurements in deep waters world-wide.Compared with that of DTU10,the root sum square of the new model for the eight main tidal constituents decrease by 8.4%to 0.997 cm.
基金supported by the Chinese Academy of Sciences Strategic Priority Research Program (Grant No. XDA05110103)the State Key Project for Basic Research Program of China (Grant No. 2010CB951801)
文摘ABSTRACT The lAP Dynamic Global Vegetation Model （IAP-DGVM） has been developed to simulate the distribution and structure of global vegetation within the framework of Earth System Models. It incorporates our group＇s recent developments of major model components such as the shrub sub-model, establishment and competition parameterization schemes, and a process-based fire parameterization of intermediate complexity. The model has 12 plant functional types, including seven tree, two shrub, and three grass types, plus bare soil. Different PFTs are allowed to coexist within a grid cell, and their state variables are updated by various governing equations describing vegetation processes from fine-scale biogeophysics and biogeochemistry, to individual and population dynamics, to large-scale biogeography. Environmental disturbance due to fire not only affects regional vegetation competition, but also influences atmospheric chemistry and aerosol emissions. Simulations under observed atmospheric conditions showed that the model can correctly reproduce the global distribution of trees, shrubs, grasses, and bare soil. The simulated global dominant vegetation types reproduce the transition from forest to grassland （savanna） in the tropical region, and from forest to shrubland in the boreal region, but overestimate the region of temperate forest.
基金supported by Chinese Academy of Sciences (KZCX2-YW-219, 100 Tal-ents Program)Ministry of Science and Technology of China (2009CB421406)
文摘The capability of an improved Dynamic Global Vegetation Model （DGVM） in reproducing the impact of climate on the terrestrial ecosystem is evaluated. The new model incorporates the Community Land Model- DGVM （CLM3.0-DGVM） with a submodel for temperate and boreal shrubs, as well as other revisions such as the ＂two-leaf＂ scheme for photosynthesis and the definition of fractional coverage of plant functional types （PFTs）. Results show that the revised model may correctly reproduce the global distribution of temperate and boreal shrubs, and improves the model performance with more realistic distribution of di？erent vege- tation types. The revised model also correctly reproduces the zonal distributions of vegetation types. In reproducing the dependence of the vegetation distribution on climate conditions, the model shows that the dominant regions for trees, grasses, shrubs, and bare soil are clearly separated by a climate index derived from mean annual precipitation and temperature, in good agreement with the CLM4 surface data. The dominant plant functional type mapping to a two dimensional parameter space of mean annual temperature and precipitation also qualitatively agrees with the results from observations and theoretical ecology studies.
基金financially supported by the Scientific Research Funds for Forestry Public Welfare of China(Granted No.201004026)the Program for Changjiang Scholars and Innovative Research Team in University(IRT1054)
文摘Heilongjiang province is the largest forest zone in China and the forest coverage rate is 46%. Forests of Heilongjiang province play an important role in the forest ecosystem of China. In this study we investi- gated the spatial distribution of forest carbon storage in Heilongjiang province using 3083 plots sampled in 2010. We attempted to fit two global models, ordinary least squares model （OLS）, linear mixed model （LMM）, and a local model, geographically weighted regression model （GWR）, to the relationship between forest carbon content and stand, environment, and climate factors. Five predictors significantly affected forest carbon storage and spatial distribution, viz. average diameter of stand （DBH）, number of trees per hectare （TPH）, elevation （Elev）, slope （Slope） and the product of precipitation and temperature （Rain Temp）. The GWR model outperformed the two global models in both model fitting and prediction because it successfully reduced both spatial auto- correlation and heterogeneity in model residuals. More importantly, the GWR model provided localized model coefficients for each location in the study area, which allowed us to evaluate the influences of local stand conditions and topographic features on tree and stand growth, and forest carbon stock. It also helped us to better understand the impacts of silvi- cultural and management activities on the amount and changes of forest carbon storage across the province. The detailed information can be readily incorporated with the mapping ability of GIS software to provide excellent tools for assessing the distribution and dynamics of the for- est-carbon stock in the next few years.
基金the National Natural Science Funds of China[grant numbers 41874030,42074026]Natural Science Funds of Shanghai[grant number 21ZR1465600]+3 种基金the Program of Shanghai Academic Research Leader[grant number 20XD1423800]the Innovation Program of Shanghai Municipal Education Commission[grant number 2021-01-07-00-07-E00095]the“Shuguang Program”supported by Shanghai Education Development Foundation and Shanghai Municipal Education Commission[grant number 20SG18]the Scientific and Technological Innovation Plan from Shanghai Science and Technology Committee[grant numbers 20511103302,20511103402 and 20511103702].
文摘With the completion of Chinese BeiDou Navigation Satellite System(BDS),the world has begun to enjoy the Positioning,Navigation,and Timing(PNT)services of four Global Navigation Satellite Systems(GNSS).In order to improve the GNSS performance and expand its applications,Low Earth Orbit(LEO)Enhanced Global Navigation Satellite System(LeGNSS)is being vigorously advocated.Combined with high-,medium-,and low-earth orbit satellites,it can improve GNSS performance in terms of orbit determination,Precise Point Positioning(PPP)convergence time,etc.This paper comprehensively reviews the current status of LeGNSS,focusing on analyzing its advantages and challenges for precise orbit and clock determination,PPP convergence,earth rotation parameter estimation,and global ionosphere modeling.Thanks to the fast geometric change brought by LEO satellites,LeGNSS is expected to fundamentally solve the problem of the long convergence time of PPP without any augmentation.The convergence time can be shortened within 1 minute if appropriate LEO constellations are deployed.However,there are still some issues to overcome,such as the optimization of LEO constellation as well as the real time LEO precise orbit and clock determination.
基金This paper is partly supported by the Chinese Academy of Sciences International Partnership Creative Group ＂The Climate System Model Development and Application Studies＂, the 973 project under Grant No. 2005CB321703 the Fund for Innovative Research Groups with Grant No. 40221503+2 种基金the National Natural Science Foundation of China under Grant Nos. 40225013the NSFC project with Grant No. 40233031 The participation of Paul J. Hanson in this work was supported by the U.S. Department of Energy （D0E）, 0ffice of Science, Biological and Environmental Research （BER）, as a part of the Program for Ecosystem Research （PER）. The data from the Walker Branch AmeriFlux tower site （Kell Wilson and Dennis Baldocchi） was developed with funding from the D0E, 0ffice of Science （BER） as a part of its Terrestrial Carbon Processes （TCP） program and from NASA/GEWEX.
文摘The interest in the development and improvement of dynamic global vegetation models （DGVMs）, which have the potential to simulate fluxes of carbon, water and nitrogen, along with changes in the vegetation dynamics, within an integrated system, has been increasing. In this paper, some numerical schemes and a higher resolution soil texture dataset were employed to improve the Sheffield Dynamic Global Vegetation Model （SDGVM）. Using eddy covariance-based measurements, we then tested the standard version of the SDGVM and the modified version of the SDGVM. Detailed observations of daily carbon and water fluxes made at the upland oak forest on the Walker Branch Watershed in Tennessee, USA offered a unique opportunity for these comparisons. The results revealed that the modified version of the SDGVM did a reasonable job of simulating the carbon and water flux and the variation of soil water content （SWC）. However, at the end of the growing season, it failed to simulate the effect of the limitations on the soil respiration dynamics and as a result underestimated this respiration. It was also noted that the modified version overestimated the increase in the SWC following summer rainfall, which was attributed to an inadequate representation of the ground water and thermal cycle.
基金supported by the National Major Research High Performance Computing Program of China(Grant No.2016YFB02008)the National Natural Science Foundation of China(Grant Number 41705070)supported by the National Natural Science Foundation of China(Grant Numbers 41475099 and 41305096)
文摘In the past several decades, dynamic global vegetation models（DGVMs） have been the most widely used and appropriate tool at the global scale to investigate vegetation-climate interactions. At the Institute of Atmospheric Physics, a new version of DGVM（IAP-DGVM） has been developed and coupled to the Common Land Model（CoLM） within the framework of the Chinese Academy of Sciences＇ Earth System Model（CAS-ESM）. This work reports the performance of IAP-DGVM through comparisons with that of the default DGVM of CoLM（CoLM-DGVM） and observations. With respect to CoLMDGVM, IAP-DGVM simulated fewer tropical trees, more ＂needleleaf evergreen boreal tree＂ and ＂broadleaf deciduous boreal shrub＂, and a better representation of grasses. These contributed to a more realistic vegetation distribution in IAP-DGVM,including spatial patterns, total areas, and compositions. Moreover, IAP-DGVM also produced more accurate carbon fluxes than CoLM-DGVM when compared with observational estimates. Gross primary productivity and net primary production in IAP-DGVM were in better agreement with observations than those of CoLM-DGVM, and the tropical pattern of fire carbon emissions in IAP-DGVM was much more consistent with the observation than that in CoLM-DGVM. The leaf area index simulated by IAP-DGVM was closer to the observation than that of CoLM-DGVM; however, both simulated values about twice as large as in the observation. This evaluation provides valuable information for the application of CAS-ESM, as well as for other model communities in terms of a comparative benchmark.
基金supported by the Chinese Academy of Sciences under Grant No.KZCX2-YW-219State Key Project for Basic Research Program of China(973)under Grant No.2010CB951801Key Program of National Natural Science Foundation under Grant No.40830103
文摘A dynamic global vegetation model （DGVM） coupled with a land surface model （LSM） is generally initialized using a spin-up process to derive a physically-consistent initial condition. Spin-up forcing, which is the atmospheric forcing used to drive the coupled model to equilibrium solutions in the spin-up process, varies across earlier studies. In the present study, the impact of the spin-up forcing in the initialization stage on the fractional coverages （FCs） of plant functional type （PFT） in the subsequent simulation stage are assessed in seven classic climate regions by a modified Community Land Model’s Dynamic Global Vegetation Model （CLM-DGVM）. Results show that the impact of spin-up forcing is considerable in all regions except the tropical rainforest climate region （TR） and the wet temperate climate region （WM）. In the tropical monsoon climate region （TM）, the TR and TM transition region （TR-TM）, the dry temperate climate region （DM）, the highland climate region （H）, and the boreal forest climate region （BF）, where FCs are affected by climate non-negligibly, the discrepancies in initial FCs, which represent long-term cumulative response of vegetation to different climate anomalies, are large. Moreover, the large discrepancies in initial FCs usually decay slowly because there are trees or shrubs in the five regions. The intrinsic growth timescales of FCs for tree PFTs and shrub PFTs are long, and the variation of FCs of tree PFTs or shrub PFTs can affect that of grass PFTs.
基金supported by the Strategic Priority Research Program of the Chinese Academy of Sciences(Grant No.XDA05110103)the State Key Project for Basic Research Program of China(Grant No.2010CB951801)the National High Technology Research and Development Program of China(863 Program)(Grant No.2009AA122105)
文摘In Dynamic Global Vegetation Models （DGVMs）, the establishment of woody vegetation refers to flowering, fertiliza- tion, seed production, germination, and the growth of tree seedlings. It determines not only the population densities but also other important ecosystem structural variables. In current DGVMs, establishments of woody plant functional types （PFTs） are assumed to be either the same in the same grid cell, or largely stochastic. We investigated the uncertainties in the competition of establishment among coexisting woody PFTs from three aspects： the dependence of PFT establishments on vegetation states; background establishment; and relative establishment potentials of different PFTs. Sensitivity experi- ments showed that the dependence of establishment rate on the fractional coverage of a PFT favored the dominant PFT by increasing its share in establishment. While a small background establishment rate had little impact on equilibrium states of the ecosystem, it did change the timescale required for the establishment of alien species in pre-existing forest due to their disadvantage in seed competition during the early stage of invasion. Meanwhile, establishment purely fiom background （the scheme commonly used in current DGVMs） led to inconsistent behavior in response to the change in PFT specification （e.g., number of PFTs and their specification）. Furthermore, the results also indicated that trade-off between irtdividual growth and reproduction/colonization has significant influences on the competition of establishment. Hence, further development of es- tablishment parameterization in DGVMs is essential in reducing the uncertainties in simulations of both ecosystem structures and successions.
基金This study is supported by the National Natural Sci-ence Foundation of China under contract No.40136010the Major State Basic Research Program of China under contract No.G1999043808the Youth Fund of National“863”Project of China under contract No.2002AA639350.
文摘To investigate the interaction between the tropical Pacific and China seas a variable-grid global ocean circulation model with fine grid covering the area from 20°S to 50°N and from 99° to 150°E is developed. Numerical computation of the annually cyclic circulation fields is performed. The results of the annual mean zonal currents and deep to abyssal western boundary currents in the equatorial Pacific Ocean are reported. The North Equatorial Current,the North Equatorial Countercurrent, the South Equatorial Current and the Equatorial Undercurrent are fairly well simulated. The model well reproduces the northward flowing abyssal western boundary current.From the model results a lower deep western boundary current east of the Bismarck-Solomon-New Hebrides Island chain at depths around 2 000 m has been found. The model results also show that the currents in the equatorial Pacific Ocean have multi-layer structures both in zonal currents and western boundary currents, indicating that the global ocean overturning thermohaline circulation appears of multi-layer pattern.
基金supported by the International Science & Technology Collaborative Program of China (2010DFB20190)Natural Science Foundation of China (41174004)the State Key Development Program of Basic Research of China (2008CB425705)
文摘On the basis of the newly released International Terrestrial Reference Frame（ITRF2008） by the International Earth Rotation Service （IERS）, a new global plate model ITRF2008 plate for the major plates is established. This ITRF2008-derived model is analyzed in comparison with NNR-NUVEL1A model, which is mainly based on geological and geophysical data. The Eurasia and Paeifi6 plates display obvious differences in terms of the velocity fields derived from the two plate motion models. Plate acceleration is also introduced to characterize the differences of the two velocity fields which obtained from ITRF2008-plate and NNR-NUVEL1A models for major individual plates. The results show that the Africa, South America and Eurasia plates are undergoing acceleration, while the North America and Australia plates are in the state of deceleration motion
基金supported by the B-type Strategic Priority Program of the Chinese Academy of Sciences (Grant No. XDB41000000)the NSFC grants 41674163, 41974186, 41704162, 41904144, and 41904143+1 种基金the pre-research projects on Civil Aerospace Technologies No. D020308, D020104 and D020303funded by China National Space Administration。
文摘Using wave measurements from the EMFISIS instrument onboard Van Allen Probes,we investigate statistically the spatial distributions of the intensity of plasmaspheric hiss waves.To reproduce these empirical results,we establish a fitting model that is a thirdorder polynomial function of L-shell,magnetic local time(MLT),magnetic latitude(MLAT),and AE*.Quantitative comparisons indicate that the model’s fitting functions can reflect favorably the major empirical features of the global distribution of hiss wave intensity,including substorm dependence and the MLT asymmetry.Our results therefore provide a useful analytic model that can be readily employed in future simulations of global radiation belt electron dynamics under the impact of plasmaspheric hiss waves in geospace.
基金supported by the Strategic Priority Research Program of Chinese Academy of Sciences (Grant No. XDA20010203)
文摘Tropical mountainous areas not only provide substantial carbon storage and play an important role in global biological diversity, but also provide basic livelihood for a large number of poor ethnic minorities. However, there is no unified and explicit definition for mountainous areas. The local elevation range(LER) is a crucial structural parameter for delineating mountainous areas. However, current LER products are limited by the subjective selection of an optimum statistical window or coarser spatial resolution of topographical data. In this study, we presented an approach using thresholds for three topographic parameters, elevation, slope, and LER, derived from the Advanced Spaceborne Thermal Emission and Reflection Radiometer Global Digital Elevation Model(ASTER GDEM) to redelineate the vast mountainous areas of mainland Southeast Asia(MSEA). The mean change-point analysis method was applied to determine the optimum statistical window of the 1 arc second(approximately 30 m)-resolution GDEM LER. The results showed that: First, the optimum statistical window is 38 × 38 cell units(width × height) in a rectangular neighborhood, or an area of about 1.30 km^2 for calculating GDEM LER in MSEA. Second, the LER of more than 80% of the area ranges from 30 m to 499 m in MSEA. The LERs in the northern and northwestern MSEA are greater than their counterparts in the south and east. Third, the area of the re-delineated mountainous areas was 83.52 × 10~4 km^2, about 38.71% of the total area. Spatially, the mountainous areas are mainly distributed in the north and northeast of MSEA. The re-delineated 30-m resolution map of the mountainous areas will serve as a topographical dataset for monitoring mountainrelated land surface changes in MSEA. The parameter-modified mountain extraction procedure can be expanded to delineate global mountainous areas.
