Different from other domestic and foreign research in which the optimum interpolation(OI) merging algorithm is commonly used,this paper constructs the non-Gaussian model for generalized variational precipitation data ...Different from other domestic and foreign research in which the optimum interpolation(OI) merging algorithm is commonly used,this paper constructs the non-Gaussian model for generalized variational precipitation data merging research based on the non-Gaussianity of precipitation data. For CMORPH data correction,the probability density function( PDF) matching method is adopted,during which the GAMMA function fitting is utilized,and the generalized variational merging based on non-Gaussian model is used to merge corrected CMORPH precipitation data and station ground observation precipitation data. Meanwhile,we carry out an experiment on CMORPH precipitation data correction and the merging of multisource precipitation data based on non-Gaussian model. By measuring the structural similarity between the merged field and the reference field,we get a merging method that can better retain useful " outliers" which represent weather phenomena. The experimental results accord with our expectations.展开更多
The paper discusses the core parameters of the 3 D and 4 D variational merging based on L1 norm regularization,namely optimization characteristic correlation length of background error covariance matrix and regulariza...The paper discusses the core parameters of the 3 D and 4 D variational merging based on L1 norm regularization,namely optimization characteristic correlation length of background error covariance matrix and regularization parameter. Classical 3 D/4 D variational merging is based on the theory that error follows Gaussian distribution. It involves the solution of the objective functional gradient in minimization iteration,which requires the data to have continuity and differentiability. Classic 3 D/4 D-dimensional variational merging method was extended,and L1 norm was used as the constraint coupling to the classical variational merged model. Experiment was carried out by using linear advection-diffusion equation as four-dimensional prediction model,and parameter optimization of this method is discussed. Considering the strong temporal and spatial variation of water vapor,this method is further applied to the precipitable water vapor( PWV) merging by calculating reanalysis data and GNSS retrieval.Parameters were adjusted gradually to analyze the influence of background field on the merging result,and the experiment results show that the mathematical algorithm adopted in this paper is feasible.展开更多
The ground-state properties of a system with a small number of interacting bosons over a wide range of densities are investigated. The system is confined in a two-dimensional isotropic harmonic trap, where the interac...The ground-state properties of a system with a small number of interacting bosons over a wide range of densities are investigated. The system is confined in a two-dimensional isotropic harmonic trap, where the interaction between bosons is treated as a hard-core potential. By using variational Monte Carlo method, we diagonalize the one-body density matrix of the system to obtain the ground-state energy, condensate wavefunction and the condensate fraction. We find that in the dilute limit the depletion of central condensate in the 2D system is larger than in a 3D system for the same interaction strength; however as the density increases, the depletion at the centre of 2D trap will be equal to or even lower than that at the centre of 3D trap, which is in agreement with the anticipated in Thomas-Fermi approximation. In addition, in the 2D system the total condensate depletion is still larger than in a 3D system for the same scattering length.展开更多
To improve traffic performance when on-ramp vehicles merge into the mainstream,a collaborative merging control strategy is proposed to determine the merging sequence and trajectory control of vehicles.Merging trajecto...To improve traffic performance when on-ramp vehicles merge into the mainstream,a collaborative merging control strategy is proposed to determine the merging sequence and trajectory control of vehicles.Merging trajectory planning takes the minimization of vehicle acceleration as the optimization objective.Either the variational method or the quadratic programming method is utilized to determine arrival time,optimal time and control variables for each vehicle.As a supplement,the adaptive cruise control(ACC)model is used to calculate each control variable in each time interval on special occasions.Simulation results show that the cooperative merging control strategy outperforms the optimal control strategy.The root mean square(RMS)of acceleration and the root mean square error(RMSE)of time headway are significantly decreased,with the reductions up to 90.1%and 25.2%,respectively.Under the cooperative control strategy,the difference between the average speed and desired speed consistently approaches zero.In addition,few or no collisions occur.To conclude,the proposed strategy favours the improvements in passenger comfort,traffic efficiency,traffic stability and safety around highway on-ramps.展开更多
A number of useful techniques associated with two-dimensional correlation spectroscopy(2DCOS)to improve its performance and utility have been developed in the last 30years.Evolution of these 2DCOS techniques,including...A number of useful techniques associated with two-dimensional correlation spectroscopy(2DCOS)to improve its performance and utility have been developed in the last 30years.Evolution of these 2DCOS techniques,including some of the very recent developments,is reviewed with examples.Topics include merged or modified asynchronous 2Dcorrelation spectrum,two-dimensional codistribution spectroscopy(2DCDS),Pareto scaling,and null-space projection treatment of spectral dataset.展开更多
This study investigates diurnal variations of precipitation during May–August, 1998–2012, over the steep slopes of the Himalayas and adjacent regions(flat Gangetic Plains–FGP, foothills of the Himalayas–FHH, the s...This study investigates diurnal variations of precipitation during May–August, 1998–2012, over the steep slopes of the Himalayas and adjacent regions(flat Gangetic Plains–FGP, foothills of the Himalayas–FHH, the steep slope of the southern Himalayas–SSSH, and the Himalayas-Tibetan Plateau tableland–HTPT). Diurnal variations are analyzed at the pixel level utilizing collocated TRMM precipitation radar and visible infrared data. The results indicate that rain parameters(including rain frequency, rain rate, and storm top altitude) are predominantly characterized by afternoon maxima and morning minima at HTPT and FGP, whereas, maximum rain parameters at FHH typically occur in the early morning. Rain parameters at SSSH are characterized by double peaks;one in the afternoon and one at midnight. Over HTPT and FGP,convective activity is strongest in the afternoon with the thickest crystallization layer. Over FHH, the vertical structure of precipitation develops most vigorously in the early morning when the most intense collision and growth of precipitation particles occurs. Over SSSH, moist convection is stronger in the afternoon and at midnight with strong mixing of ice and water particles. The results of harmonic analysis show that rain bands move southward from lower elevation of SSSH to FHH with apparent southward propagation of the harmonic phase from midnight to early morning. Moreover, the strongest diurnal harmonic is located at HTPT, having a diurnal harmonic percentage variance of up to 90%. Large-scale atmospheric circulation patterns exhibit obvious diurnal variability and correspond well to the distribution of precipitation.展开更多
Background The full lifespan of long-lived trees includes a seedling phase,during which a seed germinates and grows to a size large enough to be measured in forest inventories.Seedling populations are usually studied ...Background The full lifespan of long-lived trees includes a seedling phase,during which a seed germinates and grows to a size large enough to be measured in forest inventories.Seedling populations are usually studied separately from adult trees,and the seedling lifespan,from seed to sapling,is poorly known.In the 50-ha Barro Colorado forest plot,we started intensive censuses of seeds and seedlings in 1994 in order to merge seedling and adult demography and document complete lifespans.Methods In 17 species abundant in seedling censuses,we subdivided populations into six size classes from seed to 1cm dbh,including seeds plus five seedling stages.The smallest seedling class was subdivided by age.Censuses in two consecutive years provided transition matrices describing the probability that a seedling in one stage moved to another one year later.For each species,we averaged the transition matrix across 25 censuses and used it to project the seedling lifespan,from seed until 1cm dbh or death.Results The predicted mean survival rate of seeds to 1cm dbh varied 1000-fold across species,from 2.9×10^(−6)to 4.4×10^(−3);the median was 2.0×10^(−4).The seedling lifespan,or the average time it takes a seed to grow to 1cm dbh,varied across species from 5.1 to 53.1 years,with a median of 20.3 years.In the median species,the 10%fastest-growing seeds would reach 1cm dbh in 9.0 years,and the slowest 10%in 34.6 years.Conclusions Combining seedling results with our previous study of lifespan after 1cm dbh,we estimate that the focal species have full lifespans varying from 41 years in a gap-demanding pioneer to 320 years in one shade-tolerant species.Lifetime demography can contribute precise survival rates and lifespans to forestry models.展开更多
Since January 2012,the National Satellite Ocean Application Service has released operational wind products from the HY-2A scatterometer(HY2-SCAT),using the maximum-likelihood estimation(MLE) method with a median filte...Since January 2012,the National Satellite Ocean Application Service has released operational wind products from the HY-2A scatterometer(HY2-SCAT),using the maximum-likelihood estimation(MLE) method with a median filter. However,the quality of the winds retrieved from HY2-SCAT depends on the sub-satellite cross-track location,and poor azimuth separation in the nadir region causes particularly low-quality wind products in this region. However,an improved scheme,i.e.,a multiple solution scheme(MSS) with a two-dimensional variational analysis method(2DVAR),has been proposed by the Royal Netherlands Meteorological Institute to overcome such problems. The present study used the MSS in combination with a 2DVAR technique to retrieve wind data from HY2-SCAT observations. The parameter of the empirical probability function,used to indicate the probability of each ambiguous solution being the "true" wind,was estimated based on HY2-SCAT data,and the 2DVAR method used to remove ambiguity in the wind direction. A comparison between MSS and ECMWF winds showed larger deviations at both low wind speeds(below 4 m/s) and high wind speeds(above 17 m/s),whereas the wind direction exhibited lower bias and good stability,even at high wind speeds greater than 24 m/s. The two HY2-SCAT wind data sets,retrieved by the standard MLE and the MSS procedures were compared with buoy observations. The RMS error of wind speed and direction were 1.3 m/s and 17.4°,and 1.3 m/s and 24.0° for the MSS and MLE wind data,respectively,indicating that MSS wind data had better agreement with the buoy data. Furthermore,the distributions of wind fields for a case study of typhoon Soulik were compared,which showed that MSS winds were spatially more consistent and meteorologically better balanced than MLE winds.展开更多
Consider a two-dimensional renewal risk model,in which the claim sizes{Xk;k≥1}form a sequence of i.i.d.copies of a non-negative random vector whose two components are dependent.Suppose that the claim sizes and inter-...Consider a two-dimensional renewal risk model,in which the claim sizes{Xk;k≥1}form a sequence of i.i.d.copies of a non-negative random vector whose two components are dependent.Suppose that the claim sizes and inter-arrival times form a sequence of i.i.d.random pairs,with each pair obeying a dependence structure via the conditional distribution of the inter-arrival time given the subsequent claim size being large.Then a precise large-deviation formula of the aggregate amount of claims is obtained.展开更多
基金Supported by the Science Technology Foundation of State Grid Corporation of ChinaNatural Science Foundation of Anhui Province(1708085QD89)+1 种基金Huaihe River Basin Meteorological Open Research Fund(HRM201407)Shenyang Institute of Atmospheric Environment of China Meteorological Administration Open Fund Project(2016SYIAE14)
文摘Different from other domestic and foreign research in which the optimum interpolation(OI) merging algorithm is commonly used,this paper constructs the non-Gaussian model for generalized variational precipitation data merging research based on the non-Gaussianity of precipitation data. For CMORPH data correction,the probability density function( PDF) matching method is adopted,during which the GAMMA function fitting is utilized,and the generalized variational merging based on non-Gaussian model is used to merge corrected CMORPH precipitation data and station ground observation precipitation data. Meanwhile,we carry out an experiment on CMORPH precipitation data correction and the merging of multisource precipitation data based on non-Gaussian model. By measuring the structural similarity between the merged field and the reference field,we get a merging method that can better retain useful " outliers" which represent weather phenomena. The experimental results accord with our expectations.
基金Supported by Open Foundation Project of Shenyang Institute of Atmospheric Environment,China Meteorological Administration(2016SYIAE14)Natural Science Foundation of Anhui Province,China(1708085QD89)National Natural Science Foundation of China(41805080)
文摘The paper discusses the core parameters of the 3 D and 4 D variational merging based on L1 norm regularization,namely optimization characteristic correlation length of background error covariance matrix and regularization parameter. Classical 3 D/4 D variational merging is based on the theory that error follows Gaussian distribution. It involves the solution of the objective functional gradient in minimization iteration,which requires the data to have continuity and differentiability. Classic 3 D/4 D-dimensional variational merging method was extended,and L1 norm was used as the constraint coupling to the classical variational merged model. Experiment was carried out by using linear advection-diffusion equation as four-dimensional prediction model,and parameter optimization of this method is discussed. Considering the strong temporal and spatial variation of water vapor,this method is further applied to the precipitable water vapor( PWV) merging by calculating reanalysis data and GNSS retrieval.Parameters were adjusted gradually to analyze the influence of background field on the merging result,and the experiment results show that the mathematical algorithm adopted in this paper is feasible.
文摘The ground-state properties of a system with a small number of interacting bosons over a wide range of densities are investigated. The system is confined in a two-dimensional isotropic harmonic trap, where the interaction between bosons is treated as a hard-core potential. By using variational Monte Carlo method, we diagonalize the one-body density matrix of the system to obtain the ground-state energy, condensate wavefunction and the condensate fraction. We find that in the dilute limit the depletion of central condensate in the 2D system is larger than in a 3D system for the same interaction strength; however as the density increases, the depletion at the centre of 2D trap will be equal to or even lower than that at the centre of 3D trap, which is in agreement with the anticipated in Thomas-Fermi approximation. In addition, in the 2D system the total condensate depletion is still larger than in a 3D system for the same scattering length.
基金The Scientific Innovation Research of Graduate Students in Jiangsu Province(No.KYCX17_0145,KYCX17_0141)
文摘To improve traffic performance when on-ramp vehicles merge into the mainstream,a collaborative merging control strategy is proposed to determine the merging sequence and trajectory control of vehicles.Merging trajectory planning takes the minimization of vehicle acceleration as the optimization objective.Either the variational method or the quadratic programming method is utilized to determine arrival time,optimal time and control variables for each vehicle.As a supplement,the adaptive cruise control(ACC)model is used to calculate each control variable in each time interval on special occasions.Simulation results show that the cooperative merging control strategy outperforms the optimal control strategy.The root mean square(RMS)of acceleration and the root mean square error(RMSE)of time headway are significantly decreased,with the reductions up to 90.1%and 25.2%,respectively.Under the cooperative control strategy,the difference between the average speed and desired speed consistently approaches zero.In addition,few or no collisions occur.To conclude,the proposed strategy favours the improvements in passenger comfort,traffic efficiency,traffic stability and safety around highway on-ramps.
文摘A number of useful techniques associated with two-dimensional correlation spectroscopy(2DCOS)to improve its performance and utility have been developed in the last 30years.Evolution of these 2DCOS techniques,including some of the very recent developments,is reviewed with examples.Topics include merged or modified asynchronous 2Dcorrelation spectrum,two-dimensional codistribution spectroscopy(2DCDS),Pareto scaling,and null-space projection treatment of spectral dataset.
基金funded by the National Natural Science Foundation of China (grant no. 41705011, 91837310)the National Key R&D Program of China (2018YFC1506803, 2018YFC1507302, 2018YFC1507200)the Second Tibetan Plateau Scientific Expedition and Research (STEP) program (Grant No. 2019QZKK0104)。
文摘This study investigates diurnal variations of precipitation during May–August, 1998–2012, over the steep slopes of the Himalayas and adjacent regions(flat Gangetic Plains–FGP, foothills of the Himalayas–FHH, the steep slope of the southern Himalayas–SSSH, and the Himalayas-Tibetan Plateau tableland–HTPT). Diurnal variations are analyzed at the pixel level utilizing collocated TRMM precipitation radar and visible infrared data. The results indicate that rain parameters(including rain frequency, rain rate, and storm top altitude) are predominantly characterized by afternoon maxima and morning minima at HTPT and FGP, whereas, maximum rain parameters at FHH typically occur in the early morning. Rain parameters at SSSH are characterized by double peaks;one in the afternoon and one at midnight. Over HTPT and FGP,convective activity is strongest in the afternoon with the thickest crystallization layer. Over FHH, the vertical structure of precipitation develops most vigorously in the early morning when the most intense collision and growth of precipitation particles occurs. Over SSSH, moist convection is stronger in the afternoon and at midnight with strong mixing of ice and water particles. The results of harmonic analysis show that rain bands move southward from lower elevation of SSSH to FHH with apparent southward propagation of the harmonic phase from midnight to early morning. Moreover, the strongest diurnal harmonic is located at HTPT, having a diurnal harmonic percentage variance of up to 90%. Large-scale atmospheric circulation patterns exhibit obvious diurnal variability and correspond well to the distribution of precipitation.
基金funded by the Environmental Seed Arrival and Interspecific Associations in Seedling Sciences Program of the Smithsonian Institutionthe National Science Foundation (DEB-0075102,DEB-0823728,DEB-0640386,DEB-1242622,DEB-1464389)the Andrew Mellon Foundation,The Ohio State University,and Yale University
文摘Background The full lifespan of long-lived trees includes a seedling phase,during which a seed germinates and grows to a size large enough to be measured in forest inventories.Seedling populations are usually studied separately from adult trees,and the seedling lifespan,from seed to sapling,is poorly known.In the 50-ha Barro Colorado forest plot,we started intensive censuses of seeds and seedlings in 1994 in order to merge seedling and adult demography and document complete lifespans.Methods In 17 species abundant in seedling censuses,we subdivided populations into six size classes from seed to 1cm dbh,including seeds plus five seedling stages.The smallest seedling class was subdivided by age.Censuses in two consecutive years provided transition matrices describing the probability that a seedling in one stage moved to another one year later.For each species,we averaged the transition matrix across 25 censuses and used it to project the seedling lifespan,from seed until 1cm dbh or death.Results The predicted mean survival rate of seeds to 1cm dbh varied 1000-fold across species,from 2.9×10^(−6)to 4.4×10^(−3);the median was 2.0×10^(−4).The seedling lifespan,or the average time it takes a seed to grow to 1cm dbh,varied across species from 5.1 to 53.1 years,with a median of 20.3 years.In the median species,the 10%fastest-growing seeds would reach 1cm dbh in 9.0 years,and the slowest 10%in 34.6 years.Conclusions Combining seedling results with our previous study of lifespan after 1cm dbh,we estimate that the focal species have full lifespans varying from 41 years in a gap-demanding pioneer to 320 years in one shade-tolerant species.Lifetime demography can contribute precise survival rates and lifespans to forestry models.
基金Supported by the National High Technology Research and Development Program of China(863 Program)(No.2013AA09A505)the Shandong Joint Fund for Marine Science Research Centers(No.U1406404)+1 种基金the National Natural Science Foundation of China(No.41106152)he National Key Technology R&D Program of China(No.2013BAD13B01)
文摘Since January 2012,the National Satellite Ocean Application Service has released operational wind products from the HY-2A scatterometer(HY2-SCAT),using the maximum-likelihood estimation(MLE) method with a median filter. However,the quality of the winds retrieved from HY2-SCAT depends on the sub-satellite cross-track location,and poor azimuth separation in the nadir region causes particularly low-quality wind products in this region. However,an improved scheme,i.e.,a multiple solution scheme(MSS) with a two-dimensional variational analysis method(2DVAR),has been proposed by the Royal Netherlands Meteorological Institute to overcome such problems. The present study used the MSS in combination with a 2DVAR technique to retrieve wind data from HY2-SCAT observations. The parameter of the empirical probability function,used to indicate the probability of each ambiguous solution being the "true" wind,was estimated based on HY2-SCAT data,and the 2DVAR method used to remove ambiguity in the wind direction. A comparison between MSS and ECMWF winds showed larger deviations at both low wind speeds(below 4 m/s) and high wind speeds(above 17 m/s),whereas the wind direction exhibited lower bias and good stability,even at high wind speeds greater than 24 m/s. The two HY2-SCAT wind data sets,retrieved by the standard MLE and the MSS procedures were compared with buoy observations. The RMS error of wind speed and direction were 1.3 m/s and 17.4°,and 1.3 m/s and 24.0° for the MSS and MLE wind data,respectively,indicating that MSS wind data had better agreement with the buoy data. Furthermore,the distributions of wind fields for a case study of typhoon Soulik were compared,which showed that MSS winds were spatially more consistent and meteorologically better balanced than MLE winds.
基金by the National Social Science Foundation of China(No.20BTJ050).
文摘Consider a two-dimensional renewal risk model,in which the claim sizes{Xk;k≥1}form a sequence of i.i.d.copies of a non-negative random vector whose two components are dependent.Suppose that the claim sizes and inter-arrival times form a sequence of i.i.d.random pairs,with each pair obeying a dependence structure via the conditional distribution of the inter-arrival time given the subsequent claim size being large.Then a precise large-deviation formula of the aggregate amount of claims is obtained.
