In order to solve the so-called ＂bull-eye＂ problem caused by using a simple bilinear interpolation as an observational mapping operator in the cost function in the multigrid three-dimensional variational （3DVAR） d...In order to solve the so-called ＂bull-eye＂ problem caused by using a simple bilinear interpolation as an observational mapping operator in the cost function in the multigrid three-dimensional variational （3DVAR） data assimilation scheme, a smoothing term, equivalent to a penalty term, is introduced into the cost function to serve as a means of troubleshooting. A theoretical analysis is first performed to figure out what on earth results in the issue of ＂bull-eye＂, and then the meaning of such smoothing term is elucidated and the uniqueness of solution of the multigrid 3DVAR with the smoothing term added is discussed through the theoretical deduction for one-dimensional （1D） case, and two idealized data assimilation experiments （one- and two-dimensional （2D） cases）. By exploring the relationship between the smoothing term and the recursive filter theoretically and practically, it is revealed why satisfied analysis results can be achieved by using such proposed solution for the issue of the multigrid 3DVAR.展开更多
This work addresses the problem of estimating the states of nonlinear dynamic systems with sparse observations.We present a hybrid three-dimensional variation（3DVar） and particle piltering（PF） method,which combine...This work addresses the problem of estimating the states of nonlinear dynamic systems with sparse observations.We present a hybrid three-dimensional variation（3DVar） and particle piltering（PF） method,which combines the advantages of 3DVar and particle-based filters.By minimizing the cost function,this approach will produce a better proposal distribution of the state.Afterwards the stochastic resampling step in standard PF can be avoided through a deterministic scheme.The simulation results show that the performance of the new method is superior to the traditional ensemble Kalman filtering（EnKF） and the standard PF,especially in highly nonlinear systems.展开更多
In order to improve the efficiency of the Ocean Variational Assimilation System （OVALS）, which has been widely used in various applications, an improved OVALS （OVALS2） is developed based on the recursive filter ...In order to improve the efficiency of the Ocean Variational Assimilation System （OVALS）, which has been widely used in various applications, an improved OVALS （OVALS2） is developed based on the recursive filter （RF） algorithm. The first advantage of OVALS2 is that memory storage can be substantially reduced in practice because it implicitly computes the background error covariance matrix; the second advantage is that there is no inversion of the background error covariance by preconditioning the control variable. For comparing the effectiveness between OVALS2 and OVALS, a set of experiments was implemented by assimilating expendable bathythermograph （XBT） and ARGO data into the Tropical Pacific circulation model. The results show that the efficiency of OVALS2 is much higher than that of OVALS. The computational time and the computer storage in the assimilation process were reduced by 83% and 77%, respectively. Additionally, the corresponding results produced by the RF are almost as good as those obtained by OVALS. These results prove that OVALS2 is suitable for operational numerical oceanic forecasting.展开更多
Oceanic dissolved oxygen(DO)concentration is crucial for assessing the status of marine ecosystems.Against the backdrop of global warming,DO shows a general decrease,posing a threat to the health of marine ecosystems....Oceanic dissolved oxygen(DO)concentration is crucial for assessing the status of marine ecosystems.Against the backdrop of global warming,DO shows a general decrease,posing a threat to the health of marine ecosystems.Therefore,there is an urgent need to develop advanced tools to characterize the spatio-temporal variations of three-dimensional(3D)DO.To address this challenge,this study introduces the Light Gradient Boosting Machine(Light-GBM),combining satellite remote sensing and reanalysis data with Biogeochemical Argo data to accurately reconstruct the 3D DO structure in the Mediterranean Sea from 2010 to 2022.Various environmental parameters are incorporated as inputs,including spatiotemporal features,meteorological characteristics,and ocean color properties.The LightGBM model demonstrates excellent performance on the testing dataset with R^(2) of 0.958.The modeled DO agrees better with in-situ measurements than products from numerical models.Using the Shapley Additive exPlanations method,the contributions of input features are assessed.Sea surface temperatures provide a correlation with DO at the sea surface,while spatial coordinates supplement the view of the ocean interior.Based on the reconstructed 3D DO structure,we identify an oxygen minimum zone in the western Mediterranean that expands continuously,reaching depths of approximately 300–800 m.The western Mediterranean exhibits a significant declining trend.This study enhances marine environmental evidence by proposing a precise and cost-effective approach for reconstructing 3D DO,thereby offering insights into the dynamics of DO variations under changing climatic conditions.展开更多
Satellite altimetry missions at high latitude have opened new avenues for understanding the changes occurring over the ice-covered region.By incorporating Arctic satellite remote sensing data-including sea surface tem...Satellite altimetry missions at high latitude have opened new avenues for understanding the changes occurring over the ice-covered region.By incorporating Arctic satellite remote sensing data-including sea surface temperature(SST),sea surface height anomaly(SSHA),and sea surface salinity(SSS).This study employs a variational method to reconstruct the three-dimensional thermohaline structure of the Arctic Ocean.Compared to the Regional Arctic Reanalysis(RARE),the reconstruction well captures both the horizontal and vertical temperature and salinity structures in the Arctic.It demonstrates superior skill over RARE,when compared with Argo profiles and Ice-Tethered Profiler(ITP)observations.The reconstruction is particularly effective in ice-covered regions,where it more accurately captures the transition from Pacific water to Atlantic water compared to RARE.These findings underscore the potential of applying Arctic satellite data to reconstruct vertical thermohaline structures in the Arctic,particularly in areas due to lack of the subsurface observation reanalysis data exhibit significant biases.As Arctic satellite observations continue to advance,the applications of this method are becoming increasingly promising,which is useful for monitoring the ice-covered region environment and can be applied to oceanographic research.展开更多
A large number of observational analyses have shown that lightning data can be used to indicate areas of deep convection. It is important to assimilate observed lightning data into numerical models, so that more small...A large number of observational analyses have shown that lightning data can be used to indicate areas of deep convection. It is important to assimilate observed lightning data into numerical models, so that more small-scale information can be incorporated to improve the quality of the initial condition and the subsequent forecasts. In this study, the empirical relationship between flash rate, water vapor mixing ratio, and graupel mixing ratio was used to adjust the model relative humidity, which was then assimilated by using the three-dimensional variational data assimilation system of the Weather Research and Forecasting model in cycling mode at 10-min intervals. To find the appropriate assimilation time-window length that yielded significant improvement in both the initial conditions and subsequent forecasts, four experiments with different assimilation time-window lengths were conducted for a squall line case that occurred on 10 July 2007 in North China. It was found that 60 min was the appropriate assimilation time-window length for this case, and longer assimilation window length was unnecessary since no further improvement was present. Forecasts of 1-h accumulated precipitation during the assimilation period and the subsequent 3-h accumulated precipitation were significantly improved compared with the control experiment without lightning data assimilation. The simulated reflectivity was optimal after 30 min of the forecast, it remained optimal during the following 42 min, and the positive effect from lightning data assimilation began to diminish after 72 min of the forecast. Overall,the improvement from lightning data assimilation can be maintained for about 3 h.展开更多
Replicating the chaotic characteristics inherent in nonlinear dynamical systems via machine learning(ML)is a key challenge in this rapidly advancing interdisciplinary field.In this work,we explore the potential of var...Replicating the chaotic characteristics inherent in nonlinear dynamical systems via machine learning(ML)is a key challenge in this rapidly advancing interdisciplinary field.In this work,we explore the potential of variational quantum circuits(VQC)for learning the stochastic properties of classical nonlinear dynamical systems.Specifically,we focus on the one-and two-dimensional logistic maps,which,while simple,remain under-explored in the context of learning dynamical characteristics.Our findings reveal that,even for such simple dynamical systems,accurately replicating longterm characteristics is hindered by a pronounced sensitivity to overfitting.While increasing the parameter complexity of the ML model typically enhances short-term prediction accuracy,it also leads to a degradation in the model’s ability to replicate long-term characteristics,primarily due to the detrimental effects of overfitting on generalization power.By comparing the VQC with two widely recognized classical ML techniques,which are long short-term memory(LSTM)networks for timeseries processing and reservoir computing,we demonstrate that VQC outperforms these methods in terms of replicating long-term characteristics.Our results suggest that for the ML of dynamics,it is demanded to develop more compact and efficient models(such as VQC)rather than more complicated and large-scale ones.展开更多
Plastome variation,including single spontaneous nucleotide substitutions and single insertions/deletions,is the major source of leaf variegation in plants.Additionally,one recent study has showed that a simple plastom...Plastome variation,including single spontaneous nucleotide substitutions and single insertions/deletions,is the major source of leaf variegation in plants.Additionally,one recent study has showed that a simple plastome structural variation,which is induced by one pair of small inverted repeats,can also result in leaf variegation.Here we show a complex plastome structural variation caused by intermolecular and intramolecular recombination across three pairs of small inverted repeats accounts for leaf variegation in a widely cultivated shrub Heptapleurum ellipticum(Araliaceae).This plastome structural variation contains two deletions and two duplications,resulting in dramatic expansion of IRs,substantial contraction of LSC and loss of 11 genes that essential for photosynthesis.Plastome heteroplasmy was detected in both green and albino sectors of variegated leaves.Relative to green sectors,albino sectors in the variegated leaves exhibit significantly reduced expression for the 11 genes lost in the mutated plastome as well as 26 other genes,but significantly increased expression for one gene related to translation apparatus.Optical and transmission electron microscopy observations showed that mesophyll cells of albino sectors possess plastids lacking grana lamellae,which likely carry the mutated plastome and contribute to albinism.In both sectors,the first layer of spongy mesophyll cells beneath the lower epidermis contains normal chloroplasts,suggesting periclinal division of the lower epidermis during development.Our study demonstrates that multiple small repeats can collectively mediate intra-and inter-molecular recombination in plastome and offers a new mechanism accounting for leaf variegation in plants.展开更多
Cirrus clouds play a crucial role in the energy balance of the Earth-atmosphere system.We investigated the spatiotemporal variations of cirrus over the South China Sea(SCS)using satellite data(MOD08,MYD08,CALIPSO)and ...Cirrus clouds play a crucial role in the energy balance of the Earth-atmosphere system.We investigated the spatiotemporal variations of cirrus over the South China Sea(SCS)using satellite data(MOD08,MYD08,CALIPSO)and reanalysis data(MERRA-2)from March 2007 to February 2015(eight years).The horizontal distribution reveals lower cirrus fraction values in the northern SCS and higher values in the southern region,with minima observed in March and April and maxima sequentially occurring in August(northern SCS,NSCS),September(middle SCS,MSCS),and December(southern SCS,SSCS).Vertically,the cirrus fraction peaks in summer and reaches its lowest levels in spring.Opaque cirrus dominates during summer in the NSCS and MSCS,comprising 53.6%and 55.9%,respectively,while the SSCS exhibits a higher frequency of opaque cirrus relative to other cloud types.Subvisible cirrus clouds have the lowest frequency year-round,whereas thin cirrus is most prominent in winter in the NSCS(46.3%)and in spring in the MSCS(45.3%).A case study from September 2021 further explores the influence of ice crystal habits on brightness temperature(BT)over the SCS.Simulations utilizing five ice crystal shapes from the ARTS DDA(Atmospheric Radiative Transfer Simulator Discrete Dipole Approximation)database and the RTTOV 12.4 radiative transfer model reveal that the 8-column-aggregate shape best represents BT in the NSCS and SSCS,while the large-block-aggregate shape performs better in the SSCS.展开更多
Understanding the spatial distributions and corresponding variation mechanisms of key soil nutrients in fragile karst ecosystems can assist in promoting sustainable development.However,due to the implementation of eco...Understanding the spatial distributions and corresponding variation mechanisms of key soil nutrients in fragile karst ecosystems can assist in promoting sustainable development.However,due to the implementation of ecological restoration initiatives such as land-use conversions,novel changes in the spatial characteristics of soil nutrients remain unknown.To address this gap,we explored nutrient variations and the drivers of the variation in the 0–15 cm topsoil layer using a regional-scale sampling method in a typical karst area in northwest Guangxi Zhuang Autonomous Region,Southwest China.Descriptive statistics,geostatistics,and spatial analysis were used to assess the soil nutrient variability.The results indicated that soil organic carbon(SOC),total nitrogen(TN),total phosphorus(TP),and total potassium(TK)concentrations showed moderate variations,with coefficients of variance being 0.60,0.60,0.71,and 0.72,respectively.Moreover,they demonstrated positive spatial autocorrelations,with global Moran's indices being 0.68,0.77,0.64,and 0.68,respectively.However,local Moran's index values were low,indicating large spatial variations in soil nutrients.The best-fitting semi-variogram models for SOC,TN,TP,and TK concentrations were spherical,Gaussian,exponential,and exponential,respectively.According to the classification criteria of the Second National Soil Census in China,SOC and TN concentrations were relatively sufficient,with the proportions of rich and very rich levels being up to 90.9 and 96.0%,respectively.TP concentration was in the mediumdeficient level,with the areas of medium and deficient levels accounting for 33.7 and 30.1%of the total,respectively.TK concentration was deficient,with the cumulative area of extremely deficient,very deficient,and deficient levels accounting for 87.6%of the total area.Consequently,the terrestrial ecosystems in the study area were more vulnerable to soil P and K than soil N deficiencies.Furthermore,variance partitioning analysis of the influencing factors showed that,except for the interactions,the single effect of other soil properties accounted more for soil nutrient variations than spatial and environmental variables.These results will aid in the future management of terrestrial ecosystems.展开更多
The National Geophysical Data Center(NGDC)of the United States has collected aeromagnetic data for input into a series of geomagnetic models to improve model resolution;however,in the Tibetan Plateau region,ground-bas...The National Geophysical Data Center(NGDC)of the United States has collected aeromagnetic data for input into a series of geomagnetic models to improve model resolution;however,in the Tibetan Plateau region,ground-based observations remain insufficient to clearly reflect the characteristics of the region’s lithospheric magnetism.In this study,we evaluate the lithospheric magnetism of the Tibetan Plateau by using a 3D surface spline model based on observations from>200 newly constructed repeat stations(portable stations)to determine the spatial distribution of plateau geomagnetism,as well as its correlation with the tectonic features of the region.We analyze the relationships between M≥5 earthquakes and lithospheric magnetic field variations on the Tibetan Plateau and identify regions susceptible to strong earthquakes.We compare the geomagnetic results with those from an enhanced magnetic model(EMM2015)developed by the NGDC and provide insights into improving lithospheric magnetic field calculations in the Tibetan Plateau region.Further research reveals that these magnetic anomalies exhibit distinct differences from the magnetic-seismic correlation mechanisms observed in other tectonic settings;here,they are governed primarily by the combined effects of compressional magnetism,thermal magnetism,and deep thermal stress.This study provides new evidence of geomagnetic anomalies on the Tibetan Plateau,interprets them physically,and demonstrates their potential for identifying seismic hazard zones on the Plateau.展开更多
Dissecting quantitative traits into Mendelian factors is a great challenge in genetics.Apple fruit storability is a complex trait controlled by multi-genes with unequal effects.We previously identified62 quantitative ...Dissecting quantitative traits into Mendelian factors is a great challenge in genetics.Apple fruit storability is a complex trait controlled by multi-genes with unequal effects.We previously identified62 quantitative trait loci(QTLs)associated with apple fruit storability and genomics-assisted prediction(GAP)models were trained using 56 QTLbased markers.Here,three candidate genes,Md NAC83,Md BPM2,and Md RGLG3,were screened from the regions of QTLs with large G'value and large genetic effects.Both a 216-bp deletion and an SNP934 T/C at the promoter of Md NAC83 were associated with higher Md NAC83expression but an SNP388 G/A at the coding region significantly reduced the activity to activate the expression of the target genes Md ACO1,Md MANA3,and Md XTH28.Md BPM2 and Md RGLG3 participated in the ubiquitination of Md NAC83.SNP657 T/A of Md BPM2 and SNP167C/G of Md RGLG3 caused a reduction in the activity to ubiquitinate Md NAC83.By the addition of functional markers to the Geno Baits SNP array,the prediction accuracy of the updated GAP models increased to 0.7723/0.6231 and 0.5639/0.5345 for flesh firmness/crispness at harvest and flesh firmness/crispness retainability,respectively.The variation network involving eight simple Mendelian variations in six genes helps to gain insight into the molecular quantitative genetics,to improve breeding strategy,and to provide targets for future genome editing.展开更多
According to the Mindlin plate theory and the first-order piston theory,this work obtains accurate closed-form eigensolutions for the flutter problem of three-dimensional(3D)rectangular laminated panels.The governing ...According to the Mindlin plate theory and the first-order piston theory,this work obtains accurate closed-form eigensolutions for the flutter problem of three-dimensional(3D)rectangular laminated panels.The governing differential equations are derived by the Hamilton's variational principle,and then solved by the iterative Separation-of-Variable(i SOV)method,which are applicable to arbitrary combinations of homogeneous Boundary Conditions(BCs).However,only the simply-support,clamped and cantilever panels are considered in this work for the sake of clarity.With the closed-form eigensolutions,the flutter frequency,flutter mode and flutter boundary are presented,and the effect of shear deformation and aerodynamic damping on flutter frequencies is investigated.Besides,the relation between panel energy and the work of aerodynamic load is discussed.The numerical comparisons reveal the following.(A)The flutter eigenvalues obtained by the present method are accurate,validated by the Finite Element Method(FEM)and the Galerkin method.(B)When the span-chord ratio is larger than 3,simplifying a 3D panel to 2D(two-dimensional)panel is reasonable and the relative differences of the flutter points predicted by the two models are less than one percent.(C)The reciprocal relationship between the mechanical energy of the panel and the work done by aerodynamic load is verified by using the present flutter eigenvalues and modes,further indicating the high accuracy of the present solutions.(D)The coupling of shear deformation and aerodynamic damping prevents frequency coalescing.展开更多
In this manuscript,we consider a non-autonomous dynamical system.Using the Carathéodory structure,we define a BS dimension on an arbitrary subset and obtain a Bowen’s equation that illustrates the relation of th...In this manuscript,we consider a non-autonomous dynamical system.Using the Carathéodory structure,we define a BS dimension on an arbitrary subset and obtain a Bowen’s equation that illustrates the relation of the BS dimension to the Pesin-Pitskel topological pressure given by Nazarian[24].Moreover,we establish a variational principle and an inverse variational principle for the BS dimension of non-autonomous dynamical systems.Finally,we also get an analogue of Billingsley’s theorem for the BS dimension of non-autonomous dynamical systems.展开更多
Prevention of biological invasion requires understanding how alien species invade native communities.Although studies have identified mechanisms that underlie plant invasion in some habitats,limited attention has focu...Prevention of biological invasion requires understanding how alien species invade native communities.Although studies have identified mechanisms that underlie plant invasion in some habitats,limited attention has focused on invasion patterns along elevational gradients.In this study,we asked which factors drive the global and regional distribution of the invasive plant Galinsoga quadriradiata along elevational gradients.To answer this question,we examined whether human activities(i.e.,roads)promote G.quadriradiata invasion,how seed dispersal-related traits of G.quadriradiata change along elevation gradients,and whether G.quadriradiata has adapted to high-elevation environments through phenotypic plasticity or genetic variation.On the global scale,we found that human activities and road density positively contribute to the G.quadriradiata expansion in mountainous areas.Field surveys in China revealed significant elevational differences in the seed dispersal traits of G.quadriradiata,with higher-elevation populations exhibiting lower dispersal ability and generally lower genetic diversity.Under common conditions,high-elevation populations showed higher leaf mass ratio but lower root mass ratio and reproductive allocation.This suggests that high-elevation environments create a barrier to dispersal for G.quadriradiata,and that G.quadriradiata has adapted phenotypically to these conditions.Our study indicates that the elevational invasion pattern of G.quadriradiata is shaped by multiple factors,particularly human activities and phenotypic adaptability.In addition,our finding that G.quadriradiata invasion at high elevations is not constrained by low genetic diversity indicates that monitoring and management of G.quadriradiata in mountainous areas should be strengthened.展开更多
Amyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disorder characterized by the progressive degeneration of upper and lower motor neurons in the brainstem and spinal cord,leading to muscle weakness,para...Amyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disorder characterized by the progressive degeneration of upper and lower motor neurons in the brainstem and spinal cord,leading to muscle weakness,paralysis,and respiratory failure (Morgan and Orrell,2016).展开更多
The three-dimensional structure and the seasonal variation of the North Pacific meridional overturning circulation （NPMOC） are analyzed based on the Simple Ocean Data Assimilation data and Argo profiling float data....The three-dimensional structure and the seasonal variation of the North Pacific meridional overturning circulation （NPMOC） are analyzed based on the Simple Ocean Data Assimilation data and Argo profiling float data. The NPMOC displays a multi-cell structure with four cells in the North Pacific altogether. The TC and the STC are a strong clockwise meridional cell in the low latitude ocean and a weaker clockwise meridional cell between 7°N and 18°N, respectively, while the DTC and the subpolar cell are a weaker anticlockwise meridional cell between 3°N and 15°N and a weakest anticlockwise meridional cell between 35°N and 50°N, respectively. The DTC, the TC and the STC are all of very strong seasonal variations. As to the DTC, the southward transport is strongest in fall and weakest in spring. For the TC, the northward transport is strongest in winter and weakest in spring, while the southward transport is strongest in fall and weakest in spring, which is associated with the strong southward fiow of the DTC in fall. As the STC, the northward transport is strongest in winter and weakest in summer, while the southward transport is strongest in summer and weakest in spring. This seasonal difference may be associated with the DTC. The zonal wind stress and the east-west slope of sea level play important roles in the seasonal variations of the TC, the STC and the DTC.展开更多
As in-situ observations are sparse,targeted observations of a specific mesoscale eddy are rare.Therefore,it is difficult to study the three-dimensional structure of moving mesoscale eddies.From April to September 2014...As in-situ observations are sparse,targeted observations of a specific mesoscale eddy are rare.Therefore,it is difficult to study the three-dimensional structure of moving mesoscale eddies.From April to September 2014,an anticyclonic eddy located at 135°E-155°E,26°N-42°N was observed using 17 rapidsampling Argo floats,and the spatiotemporal variations in the three-dimensional structure were studied.The results are as follows:(1)the eddy was identified and tracked using satellite altimeter data.It had a lifetime of 269 days and an average radius of 91.5 km.The lifetime of the eddy can be divided into three phases,i.e.,the initiation,maturity,and termination phases.The depth of its influence reached 1000 m;(2)the Argo profiles were divided into seven periods(approximately 20 days in each)for composite analysis,and the composite Argo profiles and CARS2009(CSIRO Atlas of Regional Seas)climatology data were merged following the data-interpolating variational analysis(DIVA)method to reconstruct the three-dimensional structure.The temperature and salinity anomaly cores of the anticyclonic mesoscale eddy are located from 400 to 600 m.From 800 to 900 m,there is an area of low salinity at the center of the eddy.A high concentration anomaly of dissolved oxygen was located at approximately 250 m;(3)to better understand the features of the eddy and its interaction with the surroundings,we calculated the anomalous velocity of the geostrophic flow and the heat,salt,dissolved oxygen transport anomaly,and discussed the eddy's origin and its adjustments to topography.The maximum heat,salt,and oxygen transport caused by eddy were 9.37×10^11 W,3.08×10^3 kg/s,and 2.70×10^2 kg/s,which all occurred during the termination phase.This study highlights the applicability of using Argo floats to understand the three-dimensional structure thermohaline features of eddies in the North Pacific.展开更多
In this paper,a scheme of dual-Doppler radar wind analysis based on a three-dimensional variational method is proposed and performed in two steps.First,the horizontal wind field is simultaneously recovered through min...In this paper,a scheme of dual-Doppler radar wind analysis based on a three-dimensional variational method is proposed and performed in two steps.First,the horizontal wind field is simultaneously recovered through minimizing a cost function defined as a radial observation term with the standard conjugate gradient method,avoiding a weighting parameter specification step.Compared with conventional dual-Doppler wind synthesis approaches,this variational method minimizes errors caused by interpolation from radar observation to analysis grid in the iterative solution process,which is one of the main sources of errors.Then,through the accelerated Liebmann method,the vertical velocity is further reestimated as an extra step by solving the Poisson equation with impermeable conditions imposed at the ground and near the tropopause.The Poisson equation defined by the second derivative of the vertical velocity is derived from the mass continuity equation.Compared with the method proposed by O’Brien,this method is less sensitive to the uncertainty of the boundary conditions and has better stability and reliability.Furthermore,the method proposed in this paper is applied to Doppler radar observation of a squall line process.It is shown that the retrieved vertical wind profile agrees well with the vertical profile obtained with the velocity–azimuth display(VAD)method,and the retrieved radial velocity as well as the analyzed positive and negative velocity centers and horizontal wind shear of the squall line are in accord with radar observations.There is a good correspondence between the divergence field of the derived wind field and the vertical velocity.And,the horizontal and vertical circulations within and around the squall line,as well as strong updrafts,the associated downdrafts,and associated rear inflow of the bow echo,are analyzed well.It is worth mentioning that the variational method in this paper can be applied to simultaneously synthesize the three-dimensional wind field from multiple-Doppler radar observations.展开更多
In the framework of upper bound theorem of limit analysis, the progressive collapse of shallow rectangular tunnels with double-layer rock mass has been theoretically analyzed based on the three-dimensional (3D) veloci...In the framework of upper bound theorem of limit analysis, the progressive collapse of shallow rectangular tunnels with double-layer rock mass has been theoretically analyzed based on the three-dimensional (3D) velocity discontinuity surfaces. According to the virtual work principle, the difference theorem and the variation method, the collapse surface of double-layer rock mass is determined based on the Hoek-Brown failure criterion. The formula can be degenerated to a single-layer rock collapsing problem when the rock mass is homogeneous. To estimate the validity of the result, the numerical simulation software PLAXIS 3D is used to simulate the collapse of shallow tunnels with double-layer rock mass, and the comparative analysis shows that numerical results are in good agreement with upper-bound solutions. According to the results of parametric analysis, the potential range of collapse of a double-layer rock mass above a shallow cavity decreases with a decrease in A1/A2,σci1/σci2 and σtm1/σtm2 and an increase in B1/B2,γ1/γ2. The range will decrease with a decrease in support pressure q and increase with a decrease in surface overload σs. Therefore, reinforced supporting is beneficial to improve the stability of the cavity during actual construction.展开更多
基金The National Basic Research Program of China under contract No. 2013CB430304the National High-Tech R&D Program of China under contract No. 2013AA09A505the National Natural Science Foundation of China under contract Nos 41030854,40906015,40906016,41106005 and 41176003
文摘In order to solve the so-called ＂bull-eye＂ problem caused by using a simple bilinear interpolation as an observational mapping operator in the cost function in the multigrid three-dimensional variational （3DVAR） data assimilation scheme, a smoothing term, equivalent to a penalty term, is introduced into the cost function to serve as a means of troubleshooting. A theoretical analysis is first performed to figure out what on earth results in the issue of ＂bull-eye＂, and then the meaning of such smoothing term is elucidated and the uniqueness of solution of the multigrid 3DVAR with the smoothing term added is discussed through the theoretical deduction for one-dimensional （1D） case, and two idealized data assimilation experiments （one- and two-dimensional （2D） cases）. By exploring the relationship between the smoothing term and the recursive filter theoretically and practically, it is revealed why satisfied analysis results can be achieved by using such proposed solution for the issue of the multigrid 3DVAR.
基金Project supported by the National Natural Science Foundation of China (Grant No. 41105063)
文摘This work addresses the problem of estimating the states of nonlinear dynamic systems with sparse observations.We present a hybrid three-dimensional variation（3DVar） and particle piltering（PF） method,which combines the advantages of 3DVar and particle-based filters.By minimizing the cost function,this approach will produce a better proposal distribution of the state.Afterwards the stochastic resampling step in standard PF can be avoided through a deterministic scheme.The simulation results show that the performance of the new method is superior to the traditional ensemble Kalman filtering（EnKF） and the standard PF,especially in highly nonlinear systems.
基金supported by the Chinese Academy of Science(Contract No. KZCX2-YW-202)the 973 Pro-gram (Grant No. 2006CB403606)the National Natural Science Foundation of China (Grant Nos. 40606008,40776011)
文摘In order to improve the efficiency of the Ocean Variational Assimilation System （OVALS）, which has been widely used in various applications, an improved OVALS （OVALS2） is developed based on the recursive filter （RF） algorithm. The first advantage of OVALS2 is that memory storage can be substantially reduced in practice because it implicitly computes the background error covariance matrix; the second advantage is that there is no inversion of the background error covariance by preconditioning the control variable. For comparing the effectiveness between OVALS2 and OVALS, a set of experiments was implemented by assimilating expendable bathythermograph （XBT） and ARGO data into the Tropical Pacific circulation model. The results show that the efficiency of OVALS2 is much higher than that of OVALS. The computational time and the computer storage in the assimilation process were reduced by 83% and 77%, respectively. Additionally, the corresponding results produced by the RF are almost as good as those obtained by OVALS. These results prove that OVALS2 is suitable for operational numerical oceanic forecasting.
基金supported by the Central Guiding Local Science and Technology Development Fund of Shandong-Yellow River Basin(No.YDZX2023019)Shandong Natural Science Foundation of China(Nos.ZR2020QF067 and ZR2023QD073)+6 种基金the Discipline Cluster Research Project of Qingdao University“Deep mining and intelligent prediction of multimodal big data for marine ecological disasters”(No.20240604)sourced from the International Argo Program and the national programs that contribute to it(https://argo.ucsd.edu)the CMEMS(http://marine.copernicus.eu/)the CDS(https://cds.climate.copernicus.eu/)the EMODnet(https://www.emodnet-chemistry.eu/)obtained from the ERA5(https://www.ecmwf.int)derived from the Glob Colour Project(http://globcolour.info).
文摘Oceanic dissolved oxygen(DO)concentration is crucial for assessing the status of marine ecosystems.Against the backdrop of global warming,DO shows a general decrease,posing a threat to the health of marine ecosystems.Therefore,there is an urgent need to develop advanced tools to characterize the spatio-temporal variations of three-dimensional(3D)DO.To address this challenge,this study introduces the Light Gradient Boosting Machine(Light-GBM),combining satellite remote sensing and reanalysis data with Biogeochemical Argo data to accurately reconstruct the 3D DO structure in the Mediterranean Sea from 2010 to 2022.Various environmental parameters are incorporated as inputs,including spatiotemporal features,meteorological characteristics,and ocean color properties.The LightGBM model demonstrates excellent performance on the testing dataset with R^(2) of 0.958.The modeled DO agrees better with in-situ measurements than products from numerical models.Using the Shapley Additive exPlanations method,the contributions of input features are assessed.Sea surface temperatures provide a correlation with DO at the sea surface,while spatial coordinates supplement the view of the ocean interior.Based on the reconstructed 3D DO structure,we identify an oxygen minimum zone in the western Mediterranean that expands continuously,reaching depths of approximately 300–800 m.The western Mediterranean exhibits a significant declining trend.This study enhances marine environmental evidence by proposing a precise and cost-effective approach for reconstructing 3D DO,thereby offering insights into the dynamics of DO variations under changing climatic conditions.
基金The National Key R&D Program of China under contract No.2022YFE0106400the China Scholarship Council under contract No.202206710071+2 种基金the Postgraduate Research&Practice Innovation Program of Jiangsu Province under contract No.KYCX23_0657the Special Founds for Creative Research under contract No.2022C61540the Opening Project of the Key Laboratory of Marine Environmental Information Technology under contract No.521037412.
文摘Satellite altimetry missions at high latitude have opened new avenues for understanding the changes occurring over the ice-covered region.By incorporating Arctic satellite remote sensing data-including sea surface temperature(SST),sea surface height anomaly(SSHA),and sea surface salinity(SSS).This study employs a variational method to reconstruct the three-dimensional thermohaline structure of the Arctic Ocean.Compared to the Regional Arctic Reanalysis(RARE),the reconstruction well captures both the horizontal and vertical temperature and salinity structures in the Arctic.It demonstrates superior skill over RARE,when compared with Argo profiles and Ice-Tethered Profiler(ITP)observations.The reconstruction is particularly effective in ice-covered regions,where it more accurately captures the transition from Pacific water to Atlantic water compared to RARE.These findings underscore the potential of applying Arctic satellite data to reconstruct vertical thermohaline structures in the Arctic,particularly in areas due to lack of the subsurface observation reanalysis data exhibit significant biases.As Arctic satellite observations continue to advance,the applications of this method are becoming increasingly promising,which is useful for monitoring the ice-covered region environment and can be applied to oceanographic research.
基金National Key Basic Research and Development(973)Program of China(2014CB441406)National Natural Science Foundation of China(91537209 and 41675001)Basic Research Fund of Chinese Academy of Meteorological Sciences(2016Z002and 2016Y008)
文摘A large number of observational analyses have shown that lightning data can be used to indicate areas of deep convection. It is important to assimilate observed lightning data into numerical models, so that more small-scale information can be incorporated to improve the quality of the initial condition and the subsequent forecasts. In this study, the empirical relationship between flash rate, water vapor mixing ratio, and graupel mixing ratio was used to adjust the model relative humidity, which was then assimilated by using the three-dimensional variational data assimilation system of the Weather Research and Forecasting model in cycling mode at 10-min intervals. To find the appropriate assimilation time-window length that yielded significant improvement in both the initial conditions and subsequent forecasts, four experiments with different assimilation time-window lengths were conducted for a squall line case that occurred on 10 July 2007 in North China. It was found that 60 min was the appropriate assimilation time-window length for this case, and longer assimilation window length was unnecessary since no further improvement was present. Forecasts of 1-h accumulated precipitation during the assimilation period and the subsequent 3-h accumulated precipitation were significantly improved compared with the control experiment without lightning data assimilation. The simulated reflectivity was optimal after 30 min of the forecast, it remained optimal during the following 42 min, and the positive effect from lightning data assimilation began to diminish after 72 min of the forecast. Overall,the improvement from lightning data assimilation can be maintained for about 3 h.
基金Project supported in part by Beijing Natural Science Foundation(Grant No.1232025)Peng Huanwu Visiting Pro-fessor Program,and Academy for Multidisciplinary Studies,Capital Normal University.
文摘Replicating the chaotic characteristics inherent in nonlinear dynamical systems via machine learning(ML)is a key challenge in this rapidly advancing interdisciplinary field.In this work,we explore the potential of variational quantum circuits(VQC)for learning the stochastic properties of classical nonlinear dynamical systems.Specifically,we focus on the one-and two-dimensional logistic maps,which,while simple,remain under-explored in the context of learning dynamical characteristics.Our findings reveal that,even for such simple dynamical systems,accurately replicating longterm characteristics is hindered by a pronounced sensitivity to overfitting.While increasing the parameter complexity of the ML model typically enhances short-term prediction accuracy,it also leads to a degradation in the model’s ability to replicate long-term characteristics,primarily due to the detrimental effects of overfitting on generalization power.By comparing the VQC with two widely recognized classical ML techniques,which are long short-term memory(LSTM)networks for timeseries processing and reservoir computing,we demonstrate that VQC outperforms these methods in terms of replicating long-term characteristics.Our results suggest that for the ML of dynamics,it is demanded to develop more compact and efficient models(such as VQC)rather than more complicated and large-scale ones.
基金supported financially by the National Natural Science Foundation of China(31811530297 and 32170217).
文摘Plastome variation,including single spontaneous nucleotide substitutions and single insertions/deletions,is the major source of leaf variegation in plants.Additionally,one recent study has showed that a simple plastome structural variation,which is induced by one pair of small inverted repeats,can also result in leaf variegation.Here we show a complex plastome structural variation caused by intermolecular and intramolecular recombination across three pairs of small inverted repeats accounts for leaf variegation in a widely cultivated shrub Heptapleurum ellipticum(Araliaceae).This plastome structural variation contains two deletions and two duplications,resulting in dramatic expansion of IRs,substantial contraction of LSC and loss of 11 genes that essential for photosynthesis.Plastome heteroplasmy was detected in both green and albino sectors of variegated leaves.Relative to green sectors,albino sectors in the variegated leaves exhibit significantly reduced expression for the 11 genes lost in the mutated plastome as well as 26 other genes,but significantly increased expression for one gene related to translation apparatus.Optical and transmission electron microscopy observations showed that mesophyll cells of albino sectors possess plastids lacking grana lamellae,which likely carry the mutated plastome and contribute to albinism.In both sectors,the first layer of spongy mesophyll cells beneath the lower epidermis contains normal chloroplasts,suggesting periclinal division of the lower epidermis during development.Our study demonstrates that multiple small repeats can collectively mediate intra-and inter-molecular recombination in plastome and offers a new mechanism accounting for leaf variegation in plants.
基金supported by the National Natural Science Foundation of China(Grant Nos.42027804,41775026,and 41075012)。
文摘Cirrus clouds play a crucial role in the energy balance of the Earth-atmosphere system.We investigated the spatiotemporal variations of cirrus over the South China Sea(SCS)using satellite data(MOD08,MYD08,CALIPSO)and reanalysis data(MERRA-2)from March 2007 to February 2015(eight years).The horizontal distribution reveals lower cirrus fraction values in the northern SCS and higher values in the southern region,with minima observed in March and April and maxima sequentially occurring in August(northern SCS,NSCS),September(middle SCS,MSCS),and December(southern SCS,SSCS).Vertically,the cirrus fraction peaks in summer and reaches its lowest levels in spring.Opaque cirrus dominates during summer in the NSCS and MSCS,comprising 53.6%and 55.9%,respectively,while the SSCS exhibits a higher frequency of opaque cirrus relative to other cloud types.Subvisible cirrus clouds have the lowest frequency year-round,whereas thin cirrus is most prominent in winter in the NSCS(46.3%)and in spring in the MSCS(45.3%).A case study from September 2021 further explores the influence of ice crystal habits on brightness temperature(BT)over the SCS.Simulations utilizing five ice crystal shapes from the ARTS DDA(Atmospheric Radiative Transfer Simulator Discrete Dipole Approximation)database and the RTTOV 12.4 radiative transfer model reveal that the 8-column-aggregate shape best represents BT in the NSCS and SSCS,while the large-block-aggregate shape performs better in the SSCS.
基金supported by the National Natural Science Foundation of China(U2344201 and 42101316)the Natural Science Foundation of Hunan Province,China(2022JJ40866)the Outstanding Youth Project of Education Bureau of Hunan Province,China(20B613)。
文摘Understanding the spatial distributions and corresponding variation mechanisms of key soil nutrients in fragile karst ecosystems can assist in promoting sustainable development.However,due to the implementation of ecological restoration initiatives such as land-use conversions,novel changes in the spatial characteristics of soil nutrients remain unknown.To address this gap,we explored nutrient variations and the drivers of the variation in the 0–15 cm topsoil layer using a regional-scale sampling method in a typical karst area in northwest Guangxi Zhuang Autonomous Region,Southwest China.Descriptive statistics,geostatistics,and spatial analysis were used to assess the soil nutrient variability.The results indicated that soil organic carbon(SOC),total nitrogen(TN),total phosphorus(TP),and total potassium(TK)concentrations showed moderate variations,with coefficients of variance being 0.60,0.60,0.71,and 0.72,respectively.Moreover,they demonstrated positive spatial autocorrelations,with global Moran's indices being 0.68,0.77,0.64,and 0.68,respectively.However,local Moran's index values were low,indicating large spatial variations in soil nutrients.The best-fitting semi-variogram models for SOC,TN,TP,and TK concentrations were spherical,Gaussian,exponential,and exponential,respectively.According to the classification criteria of the Second National Soil Census in China,SOC and TN concentrations were relatively sufficient,with the proportions of rich and very rich levels being up to 90.9 and 96.0%,respectively.TP concentration was in the mediumdeficient level,with the areas of medium and deficient levels accounting for 33.7 and 30.1%of the total,respectively.TK concentration was deficient,with the cumulative area of extremely deficient,very deficient,and deficient levels accounting for 87.6%of the total area.Consequently,the terrestrial ecosystems in the study area were more vulnerable to soil P and K than soil N deficiencies.Furthermore,variance partitioning analysis of the influencing factors showed that,except for the interactions,the single effect of other soil properties accounted more for soil nutrient variations than spatial and environmental variables.These results will aid in the future management of terrestrial ecosystems.
基金supported by the CAS Pioneer Hundred Talents Program and Second Tibetan Plateau Scientific Expedition Research Program(2019QZKK0708)as well as the Basic Research Program of Qinghai Province:Lithospheric Geomagnetic Field of the Qinghai-Tibet Plateau and the Relationship with Strong Earthquakes(2021-ZJ-969Q).
文摘The National Geophysical Data Center(NGDC)of the United States has collected aeromagnetic data for input into a series of geomagnetic models to improve model resolution;however,in the Tibetan Plateau region,ground-based observations remain insufficient to clearly reflect the characteristics of the region’s lithospheric magnetism.In this study,we evaluate the lithospheric magnetism of the Tibetan Plateau by using a 3D surface spline model based on observations from>200 newly constructed repeat stations(portable stations)to determine the spatial distribution of plateau geomagnetism,as well as its correlation with the tectonic features of the region.We analyze the relationships between M≥5 earthquakes and lithospheric magnetic field variations on the Tibetan Plateau and identify regions susceptible to strong earthquakes.We compare the geomagnetic results with those from an enhanced magnetic model(EMM2015)developed by the NGDC and provide insights into improving lithospheric magnetic field calculations in the Tibetan Plateau region.Further research reveals that these magnetic anomalies exhibit distinct differences from the magnetic-seismic correlation mechanisms observed in other tectonic settings;here,they are governed primarily by the combined effects of compressional magnetism,thermal magnetism,and deep thermal stress.This study provides new evidence of geomagnetic anomalies on the Tibetan Plateau,interprets them physically,and demonstrates their potential for identifying seismic hazard zones on the Plateau.
基金financially supported by the National Natural Science Foundation of China(32202431)the National Key R&D Program of China(2022YFD1200503)+2 种基金China Postdoctoral Science Foundation(2022M713408)the Earmarked Fund for CARS-27the Key Research and Development Program of Hebei(21326308D)。
文摘Dissecting quantitative traits into Mendelian factors is a great challenge in genetics.Apple fruit storability is a complex trait controlled by multi-genes with unequal effects.We previously identified62 quantitative trait loci(QTLs)associated with apple fruit storability and genomics-assisted prediction(GAP)models were trained using 56 QTLbased markers.Here,three candidate genes,Md NAC83,Md BPM2,and Md RGLG3,were screened from the regions of QTLs with large G'value and large genetic effects.Both a 216-bp deletion and an SNP934 T/C at the promoter of Md NAC83 were associated with higher Md NAC83expression but an SNP388 G/A at the coding region significantly reduced the activity to activate the expression of the target genes Md ACO1,Md MANA3,and Md XTH28.Md BPM2 and Md RGLG3 participated in the ubiquitination of Md NAC83.SNP657 T/A of Md BPM2 and SNP167C/G of Md RGLG3 caused a reduction in the activity to ubiquitinate Md NAC83.By the addition of functional markers to the Geno Baits SNP array,the prediction accuracy of the updated GAP models increased to 0.7723/0.6231 and 0.5639/0.5345 for flesh firmness/crispness at harvest and flesh firmness/crispness retainability,respectively.The variation network involving eight simple Mendelian variations in six genes helps to gain insight into the molecular quantitative genetics,to improve breeding strategy,and to provide targets for future genome editing.
基金support of the National Natural Science Foundation of China(No.12172023)。
文摘According to the Mindlin plate theory and the first-order piston theory,this work obtains accurate closed-form eigensolutions for the flutter problem of three-dimensional(3D)rectangular laminated panels.The governing differential equations are derived by the Hamilton's variational principle,and then solved by the iterative Separation-of-Variable(i SOV)method,which are applicable to arbitrary combinations of homogeneous Boundary Conditions(BCs).However,only the simply-support,clamped and cantilever panels are considered in this work for the sake of clarity.With the closed-form eigensolutions,the flutter frequency,flutter mode and flutter boundary are presented,and the effect of shear deformation and aerodynamic damping on flutter frequencies is investigated.Besides,the relation between panel energy and the work of aerodynamic load is discussed.The numerical comparisons reveal the following.(A)The flutter eigenvalues obtained by the present method are accurate,validated by the Finite Element Method(FEM)and the Galerkin method.(B)When the span-chord ratio is larger than 3,simplifying a 3D panel to 2D(two-dimensional)panel is reasonable and the relative differences of the flutter points predicted by the two models are less than one percent.(C)The reciprocal relationship between the mechanical energy of the panel and the work done by aerodynamic load is verified by using the present flutter eigenvalues and modes,further indicating the high accuracy of the present solutions.(D)The coupling of shear deformation and aerodynamic damping prevents frequency coalescing.
基金supported by the NSFC(12461012)and the NSF of Chongqing(CSTB2024NSCQ-MSX1246).
文摘In this manuscript,we consider a non-autonomous dynamical system.Using the Carathéodory structure,we define a BS dimension on an arbitrary subset and obtain a Bowen’s equation that illustrates the relation of the BS dimension to the Pesin-Pitskel topological pressure given by Nazarian[24].Moreover,we establish a variational principle and an inverse variational principle for the BS dimension of non-autonomous dynamical systems.Finally,we also get an analogue of Billingsley’s theorem for the BS dimension of non-autonomous dynamical systems.
基金supported by the National Natural Science Foundation of China(32271584 and 31600445)the Natural Science Basic Research Plan in Shaanxi Province of China(2020JM-286)+2 种基金the Fundamental Research Funds for the Central Universities(GK202103072,GK202103073)the National College Students'Innovative Entrepreneurial Training Plan Program(202310718085)Special Research Project in Philosophy and Social Sciences of Shaanxi Province(2022HZ1795).
文摘Prevention of biological invasion requires understanding how alien species invade native communities.Although studies have identified mechanisms that underlie plant invasion in some habitats,limited attention has focused on invasion patterns along elevational gradients.In this study,we asked which factors drive the global and regional distribution of the invasive plant Galinsoga quadriradiata along elevational gradients.To answer this question,we examined whether human activities(i.e.,roads)promote G.quadriradiata invasion,how seed dispersal-related traits of G.quadriradiata change along elevation gradients,and whether G.quadriradiata has adapted to high-elevation environments through phenotypic plasticity or genetic variation.On the global scale,we found that human activities and road density positively contribute to the G.quadriradiata expansion in mountainous areas.Field surveys in China revealed significant elevational differences in the seed dispersal traits of G.quadriradiata,with higher-elevation populations exhibiting lower dispersal ability and generally lower genetic diversity.Under common conditions,high-elevation populations showed higher leaf mass ratio but lower root mass ratio and reproductive allocation.This suggests that high-elevation environments create a barrier to dispersal for G.quadriradiata,and that G.quadriradiata has adapted phenotypically to these conditions.Our study indicates that the elevational invasion pattern of G.quadriradiata is shaped by multiple factors,particularly human activities and phenotypic adaptability.In addition,our finding that G.quadriradiata invasion at high elevations is not constrained by low genetic diversity indicates that monitoring and management of G.quadriradiata in mountainous areas should be strengthened.
文摘Amyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disorder characterized by the progressive degeneration of upper and lower motor neurons in the brainstem and spinal cord,leading to muscle weakness,paralysis,and respiratory failure (Morgan and Orrell,2016).
基金Supported by the National Basic Research Development Program of China(973 Program)under contract Nos 2007CB816002,2007CB816005the innovative key project of Chinese Academy of Sciences under contract No.KZCXZ-YW-201
文摘The three-dimensional structure and the seasonal variation of the North Pacific meridional overturning circulation （NPMOC） are analyzed based on the Simple Ocean Data Assimilation data and Argo profiling float data. The NPMOC displays a multi-cell structure with four cells in the North Pacific altogether. The TC and the STC are a strong clockwise meridional cell in the low latitude ocean and a weaker clockwise meridional cell between 7°N and 18°N, respectively, while the DTC and the subpolar cell are a weaker anticlockwise meridional cell between 3°N and 15°N and a weakest anticlockwise meridional cell between 35°N and 50°N, respectively. The DTC, the TC and the STC are all of very strong seasonal variations. As to the DTC, the southward transport is strongest in fall and weakest in spring. For the TC, the northward transport is strongest in winter and weakest in spring, while the southward transport is strongest in fall and weakest in spring, which is associated with the strong southward fiow of the DTC in fall. As the STC, the northward transport is strongest in winter and weakest in summer, while the southward transport is strongest in summer and weakest in spring. This seasonal difference may be associated with the DTC. The zonal wind stress and the east-west slope of sea level play important roles in the seasonal variations of the TC, the STC and the DTC.
基金Supported by the National Key R&D Program of China(No.2018YFC1406202)the National Natural Science Foundation of China(Nos.41830964,41976188,41605051)。
文摘As in-situ observations are sparse,targeted observations of a specific mesoscale eddy are rare.Therefore,it is difficult to study the three-dimensional structure of moving mesoscale eddies.From April to September 2014,an anticyclonic eddy located at 135°E-155°E,26°N-42°N was observed using 17 rapidsampling Argo floats,and the spatiotemporal variations in the three-dimensional structure were studied.The results are as follows:(1)the eddy was identified and tracked using satellite altimeter data.It had a lifetime of 269 days and an average radius of 91.5 km.The lifetime of the eddy can be divided into three phases,i.e.,the initiation,maturity,and termination phases.The depth of its influence reached 1000 m;(2)the Argo profiles were divided into seven periods(approximately 20 days in each)for composite analysis,and the composite Argo profiles and CARS2009(CSIRO Atlas of Regional Seas)climatology data were merged following the data-interpolating variational analysis(DIVA)method to reconstruct the three-dimensional structure.The temperature and salinity anomaly cores of the anticyclonic mesoscale eddy are located from 400 to 600 m.From 800 to 900 m,there is an area of low salinity at the center of the eddy.A high concentration anomaly of dissolved oxygen was located at approximately 250 m;(3)to better understand the features of the eddy and its interaction with the surroundings,we calculated the anomalous velocity of the geostrophic flow and the heat,salt,dissolved oxygen transport anomaly,and discussed the eddy's origin and its adjustments to topography.The maximum heat,salt,and oxygen transport caused by eddy were 9.37×10^11 W,3.08×10^3 kg/s,and 2.70×10^2 kg/s,which all occurred during the termination phase.This study highlights the applicability of using Argo floats to understand the three-dimensional structure thermohaline features of eddies in the North Pacific.
基金supported by the National Key Research and Development Program of China(Grant No.2019YFC1510400)the National Natural Science Foundation of China(Grant Nos.41975054 and 41930967)the Special Fund for Forecasters of China Meteorological Administration(Grant No.CMAYBY2018-040)。
文摘In this paper,a scheme of dual-Doppler radar wind analysis based on a three-dimensional variational method is proposed and performed in two steps.First,the horizontal wind field is simultaneously recovered through minimizing a cost function defined as a radial observation term with the standard conjugate gradient method,avoiding a weighting parameter specification step.Compared with conventional dual-Doppler wind synthesis approaches,this variational method minimizes errors caused by interpolation from radar observation to analysis grid in the iterative solution process,which is one of the main sources of errors.Then,through the accelerated Liebmann method,the vertical velocity is further reestimated as an extra step by solving the Poisson equation with impermeable conditions imposed at the ground and near the tropopause.The Poisson equation defined by the second derivative of the vertical velocity is derived from the mass continuity equation.Compared with the method proposed by O’Brien,this method is less sensitive to the uncertainty of the boundary conditions and has better stability and reliability.Furthermore,the method proposed in this paper is applied to Doppler radar observation of a squall line process.It is shown that the retrieved vertical wind profile agrees well with the vertical profile obtained with the velocity–azimuth display(VAD)method,and the retrieved radial velocity as well as the analyzed positive and negative velocity centers and horizontal wind shear of the squall line are in accord with radar observations.There is a good correspondence between the divergence field of the derived wind field and the vertical velocity.And,the horizontal and vertical circulations within and around the squall line,as well as strong updrafts,the associated downdrafts,and associated rear inflow of the bow echo,are analyzed well.It is worth mentioning that the variational method in this paper can be applied to simultaneously synthesize the three-dimensional wind field from multiple-Doppler radar observations.
基金Projects(51478477,51878074)supported by the National Natural Science Foundation of ChinaProject(2017-123-033)supported by the Guizhou Provincial Department of Transportation Foundation,ChinaProjects(2018zzts663,2018zzts656)supported by the Fundamental Research Funds for the Central Universities,China
文摘In the framework of upper bound theorem of limit analysis, the progressive collapse of shallow rectangular tunnels with double-layer rock mass has been theoretically analyzed based on the three-dimensional (3D) velocity discontinuity surfaces. According to the virtual work principle, the difference theorem and the variation method, the collapse surface of double-layer rock mass is determined based on the Hoek-Brown failure criterion. The formula can be degenerated to a single-layer rock collapsing problem when the rock mass is homogeneous. To estimate the validity of the result, the numerical simulation software PLAXIS 3D is used to simulate the collapse of shallow tunnels with double-layer rock mass, and the comparative analysis shows that numerical results are in good agreement with upper-bound solutions. According to the results of parametric analysis, the potential range of collapse of a double-layer rock mass above a shallow cavity decreases with a decrease in A1/A2,σci1/σci2 and σtm1/σtm2 and an increase in B1/B2,γ1/γ2. The range will decrease with a decrease in support pressure q and increase with a decrease in surface overload σs. Therefore, reinforced supporting is beneficial to improve the stability of the cavity during actual construction.
