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.展开更多
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.展开更多
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.展开更多
Under the strategic framework of rural revitalization and agricultural modernization, Xinjiashan Specialty Coffee Base, located in Zaotang Village, Lujiang Town, Longyang District, Baoshan City, has been proactively i...Under the strategic framework of rural revitalization and agricultural modernization, Xinjiashan Specialty Coffee Base, located in Zaotang Village, Lujiang Town, Longyang District, Baoshan City, has been proactively investigating innovative models for agricultural development. Through extensive communication and collaboration, this base has established close partnerships with research institutions including Kunming University of Science and Technology, Baoshan University, and Yunnan Academy of Agricultural Sciences, with a commitment to thoroughly exploring the potential for resource recycling and ecological complementarity. An innovative four-in-one three-dimensional integrated planting system incorporating "coffee, bananas, green manure, and bees" has been implemented. Concurrently, technological and digital management strategies have been comprehensively integrated to improve planting efficiency. Under this model, the proportion of specialty coffee attains 71%, and the per-unit yield is 17% greater than that of the conventional planting model. This approach not only substantially enhances economic returns but also promotes the integrated development of ecological and social benefits, offering a valuable practical example and experiential reference for the specialty and sustainable advancement of the coffee industry in comparable regions.展开更多
Aerial surveys are dynamic and continuous processes,and there are different height distributions of the ground in the measurement area,which leads to problems such as overlapping measurement areas and inaccurate altit...Aerial surveys are dynamic and continuous processes,and there are different height distributions of the ground in the measurement area,which leads to problems such as overlapping measurement areas and inaccurate altitude correction during the survey process.Commonly used terrain correction methods are based on the concept of finite elementization of ground surface radioactive sources,using GPS coordinates,radar altitude,and ground elevation distribution information from aerial surveys,combined with the sourceless efficiency calibration method to construct a response matrix,which is then inverted for surface nuclide content.However,most of the sourceless efficiency calibration methods used are numerical calculations that consider the body detector as a point detector and do not consider the changes in intrinsic detection efficiency under different incident directions of gamma rays.Therefore,when the altitude of the measurement area varies significantly or the flight altitude of the aerial survey is relatively low,such sourceless efficiency calibration method calculations tend to have a large bias,which affects the accuracy of the terrain correction.To address the above problems,this study employs a novel sourceless efficiency calibration method based on the Boolean operation of the ray deposition process and simplifies the traditional body source measurement model to a surface source measurement model to achieve fast and accurate efficiency calibration.Then,through the discretization of the measurement process,the static measurement process is superposed as equivalent to the dynamic measurement process,and the dynamic measurement response matrix is built and optimized based on the calibration method.Finally,the PSO-MLEM algorithm was used to solve the dynamic measurement response matrix to achieve dynamic terrain correction of aerial survey data.Analysis of the Baiyun'ebo test area revealed that,after applying dynamic terrain correction,the inverted anomalies in uranium(eU),thorium(eTh),and potassium(K)concentrations were closer to ground measurements(within 5.72%-30.79%)and exhibited clearer anomaly boundaries compared to traditional height-based corrections.However,owing to the inherent statistical fluctuations and characteristics of matrix inversion,higher measurement values tend to absorb lower ones,potentially enlarging the anomalous regions.Nevertheless,the highanomaly regions after inversion largely coincided with the ground truth validation,demonstrating that the proposed method can effectively correct airborne gamma spectrometry data.展开更多
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.展开更多
Signal transduction in a cell is mostly mediated with biochemical reactions which are noisy and often modeled with chemical master equations or Langevin type of dynamics.Thus stochastic simulation constitutes a major ...Signal transduction in a cell is mostly mediated with biochemical reactions which are noisy and often modeled with chemical master equations or Langevin type of dynamics.Thus stochastic simulation constitutes a major part of computation in cell signaling.Nevertheless,the presence of a wide span of time scales or molecular numbers in various pathways may lead to trouble in computation efficiency or accuracy.To avoid this problem,the commonly employed variational method evolves the whole probability distribution and reduces the stochastic equations to deterministic ones of only a few parameters.However,the design of the left variational basis is essential for its successful application,especially to large networks.In this paper,we extend the conventional polynomial basis to the Fourier and further the Gaussian basis,much facilitating description of multi-peaked or localized non-Gaussian distributions and at the same time avoiding numerical instability and computational complexity frequently encountered with conventional basis.The extension here is demonstrated in several typical biochemical signaling networks and achieves similar accuracy as the benchmark Gillespie algorithm,but with much less running time,which seems to open new opportunities in the variational approach to efficient analysis of noisy dynamics.展开更多
Variational quantum algorithms(VQAs)with random structures have poor trainability due to the exponentially vanishing gradient as the circuit depth and the qubit number increase.This result leads to a general belief th...Variational quantum algorithms(VQAs)with random structures have poor trainability due to the exponentially vanishing gradient as the circuit depth and the qubit number increase.This result leads to a general belief that a deep circuit will not be feasible.In this work,we provide a viable solution to the vanishing gradient problem for deep VQAs with theoretical guarantees.Specifically,we prove that for quantum controlled-layer and quantum residual network(QResNet),architectures,the expectation of the gradient norm can be lower bounded by a value that is independent of the qubit number and the circuit depth.Our results follow from a careful analysis of the gradient behavior on parameter space consisting of rotation angles,as employed in almost all VQAs,instead of relying on impractical 2-design assumptions.We conduct several numerical experiments as verifications,where only our circuits are trainable and converge,while hardware-efficient and random circuits with similar number of parameters in comparison cannot converge.展开更多
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.展开更多
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.展开更多
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.展开更多
The scientific design and preliminary results of the data assimilation component of the Global-Regional Prediction and Assimilation System (GRAPES) recently developed in China Meteorological Administration (CMA) are p...The scientific design and preliminary results of the data assimilation component of the Global-Regional Prediction and Assimilation System (GRAPES) recently developed in China Meteorological Administration (CMA) are presented in this paper. This is a three-dimensional variational (3DVar) assimilation system set up on global and regional grid meshes favorable for direct assimilation of the space-based remote sensing data and matching the frame work of the prediction model GRAPES. The state variables are assumed to decompose balanced and unbalanced components. By introducing a simple transformation from the state variables to the control variables with a recursive or spectral filter, the convergence rate of iteration for minimization of the cost function in 3DVar is greatly accelerated. The definition of dynamical balance depends on the characteristic scale of the circulation considered. The ratio of the balanced to the unbalanced parts is controlled by the prescribed statistics of background errors. Idealized trials produce the same results as the analytic solution. The results of real data case studies show the capability of the system to improve analysis compared to the traditional schemes. Finally, further development of the system is discussed.展开更多
The dynamical constrains in three-dimensional variational data assimilation are discussed when considering the impact of stream divergence and convergence on the pressure and wind fields. For the analysis of severe tr...The dynamical constrains in three-dimensional variational data assimilation are discussed when considering the impact of stream divergence and convergence on the pressure and wind fields. For the analysis of severe tropical cyclone, frontal structures, and other rapidly changing structures, the geostrophic balance and linear balance cannot properly represent the relationship between wind and pressure fields. However, the nonlinear balance incremental equation takes into account the information of flow-dependent background, and makes response to the flow-dependent background covariance in the 3D-Var system. Results indicate that the application of the nonlinear balance equation to 3D-Var system improves the quality of severe tropical cyclone assimilation system, which has some positive effects on intensity prediction of tropical cyclones.展开更多
By introducing the dimensional splitting(DS)method into the multiscale interpolating element-free Galerkin(VMIEFG)method,a dimension-splitting multiscale interpolating element-free Galerkin(DS-VMIEFG)method is propose...By introducing the dimensional splitting(DS)method into the multiscale interpolating element-free Galerkin(VMIEFG)method,a dimension-splitting multiscale interpolating element-free Galerkin(DS-VMIEFG)method is proposed for three-dimensional(3D)singular perturbed convection-diffusion(SPCD)problems.In the DSVMIEFG method,the 3D problem is decomposed into a series of 2D problems by the DS method,and the discrete equations on the 2D splitting surface are obtained by the VMIEFG method.The improved interpolation-type moving least squares(IIMLS)method is used to construct shape functions in the weak form and to combine 2D discrete equations into a global system of discrete equations for the three-dimensional SPCD problems.The solved numerical example verifies the effectiveness of the method in this paper for the 3D SPCD problems.The numerical solution will gradually converge to the analytical solution with the increase in the number of nodes.For extremely small singular diffusion coefficients,the numerical solution will avoid numerical oscillation and has high computational stability.展开更多
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.展开更多
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.展开更多
To address the problem of multi-missile cooperative interception against maneuvering targets at a prespecified impact time and desired Line-of-Sight(LOS)angles in ThreeDimensional(3D)space,this paper proposes a 3D lea...To address the problem of multi-missile cooperative interception against maneuvering targets at a prespecified impact time and desired Line-of-Sight(LOS)angles in ThreeDimensional(3D)space,this paper proposes a 3D leader-following cooperative interception guidance law.First,in the LOS direction of the leader,an impact time-controlled guidance law is derived based on the fixed-time stability theory,which enables the leader to complete the interception task at a prespecified impact time.Next,in the LOS direction of the followers,by introducing a time consensus tracking error function,a fixed-time consensus tracking guidance law is investigated to guarantee the consensus tracking convergence of the time-to-go.Then,in the direction normal to the LOS,by combining the designed global integral sliding mode surface and the second-order Sliding Mode Control(SMC)theory,an innovative 3D LOS-angle-constrained interception guidance law is developed,which eliminates the reaching phase in the traditional sliding mode guidance laws and effectively saves energy consumption.Moreover,it effectively suppresses the chattering phenomenon while avoiding the singularity issue,and compensates for unknown interference caused by target maneuvering online,making it convenient for practical engineering applications.Finally,theoretical proof analysis and multiple sets of numerical simulation results verify the effectiveness,superiority,and robustness of the investigated guidance law.展开更多
基金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.
基金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.
基金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.
文摘Under the strategic framework of rural revitalization and agricultural modernization, Xinjiashan Specialty Coffee Base, located in Zaotang Village, Lujiang Town, Longyang District, Baoshan City, has been proactively investigating innovative models for agricultural development. Through extensive communication and collaboration, this base has established close partnerships with research institutions including Kunming University of Science and Technology, Baoshan University, and Yunnan Academy of Agricultural Sciences, with a commitment to thoroughly exploring the potential for resource recycling and ecological complementarity. An innovative four-in-one three-dimensional integrated planting system incorporating "coffee, bananas, green manure, and bees" has been implemented. Concurrently, technological and digital management strategies have been comprehensively integrated to improve planting efficiency. Under this model, the proportion of specialty coffee attains 71%, and the per-unit yield is 17% greater than that of the conventional planting model. This approach not only substantially enhances economic returns but also promotes the integrated development of ecological and social benefits, offering a valuable practical example and experiential reference for the specialty and sustainable advancement of the coffee industry in comparable regions.
基金supported by the National Key Research and Development Program(No.2022YFC2807400)the National Natural Science Foundation of China(Nos.12265003 and 12205044)。
文摘Aerial surveys are dynamic and continuous processes,and there are different height distributions of the ground in the measurement area,which leads to problems such as overlapping measurement areas and inaccurate altitude correction during the survey process.Commonly used terrain correction methods are based on the concept of finite elementization of ground surface radioactive sources,using GPS coordinates,radar altitude,and ground elevation distribution information from aerial surveys,combined with the sourceless efficiency calibration method to construct a response matrix,which is then inverted for surface nuclide content.However,most of the sourceless efficiency calibration methods used are numerical calculations that consider the body detector as a point detector and do not consider the changes in intrinsic detection efficiency under different incident directions of gamma rays.Therefore,when the altitude of the measurement area varies significantly or the flight altitude of the aerial survey is relatively low,such sourceless efficiency calibration method calculations tend to have a large bias,which affects the accuracy of the terrain correction.To address the above problems,this study employs a novel sourceless efficiency calibration method based on the Boolean operation of the ray deposition process and simplifies the traditional body source measurement model to a surface source measurement model to achieve fast and accurate efficiency calibration.Then,through the discretization of the measurement process,the static measurement process is superposed as equivalent to the dynamic measurement process,and the dynamic measurement response matrix is built and optimized based on the calibration method.Finally,the PSO-MLEM algorithm was used to solve the dynamic measurement response matrix to achieve dynamic terrain correction of aerial survey data.Analysis of the Baiyun'ebo test area revealed that,after applying dynamic terrain correction,the inverted anomalies in uranium(eU),thorium(eTh),and potassium(K)concentrations were closer to ground measurements(within 5.72%-30.79%)and exhibited clearer anomaly boundaries compared to traditional height-based corrections.However,owing to the inherent statistical fluctuations and characteristics of matrix inversion,higher measurement values tend to absorb lower ones,potentially enlarging the anomalous regions.Nevertheless,the highanomaly regions after inversion largely coincided with the ground truth validation,demonstrating that the proposed method can effectively correct airborne gamma spectrometry data.
基金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 by the National Natural Science Foundation of China under Grants No.12375030。
文摘Signal transduction in a cell is mostly mediated with biochemical reactions which are noisy and often modeled with chemical master equations or Langevin type of dynamics.Thus stochastic simulation constitutes a major part of computation in cell signaling.Nevertheless,the presence of a wide span of time scales or molecular numbers in various pathways may lead to trouble in computation efficiency or accuracy.To avoid this problem,the commonly employed variational method evolves the whole probability distribution and reduces the stochastic equations to deterministic ones of only a few parameters.However,the design of the left variational basis is essential for its successful application,especially to large networks.In this paper,we extend the conventional polynomial basis to the Fourier and further the Gaussian basis,much facilitating description of multi-peaked or localized non-Gaussian distributions and at the same time avoiding numerical instability and computational complexity frequently encountered with conventional basis.The extension here is demonstrated in several typical biochemical signaling networks and achieves similar accuracy as the benchmark Gillespie algorithm,but with much less running time,which seems to open new opportunities in the variational approach to efficient analysis of noisy dynamics.
基金The national research foundation of Singapore(NRF-P2024-001).
文摘Variational quantum algorithms(VQAs)with random structures have poor trainability due to the exponentially vanishing gradient as the circuit depth and the qubit number increase.This result leads to a general belief that a deep circuit will not be feasible.In this work,we provide a viable solution to the vanishing gradient problem for deep VQAs with theoretical guarantees.Specifically,we prove that for quantum controlled-layer and quantum residual network(QResNet),architectures,the expectation of the gradient norm can be lower bounded by a value that is independent of the qubit number and the circuit depth.Our results follow from a careful analysis of the gradient behavior on parameter space consisting of rotation angles,as employed in almost all VQAs,instead of relying on impractical 2-design assumptions.We conduct several numerical experiments as verifications,where only our circuits are trainable and converge,while hardware-efficient and random circuits with similar number of parameters in comparison cannot converge.
基金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 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.
基金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.
基金Key Technologies Research and Development Program (Grant No. 2001BA607B and 2001BA607B02)National Natural Science Foundation of China (Grant No. 40518001)
文摘The scientific design and preliminary results of the data assimilation component of the Global-Regional Prediction and Assimilation System (GRAPES) recently developed in China Meteorological Administration (CMA) are presented in this paper. This is a three-dimensional variational (3DVar) assimilation system set up on global and regional grid meshes favorable for direct assimilation of the space-based remote sensing data and matching the frame work of the prediction model GRAPES. The state variables are assumed to decompose balanced and unbalanced components. By introducing a simple transformation from the state variables to the control variables with a recursive or spectral filter, the convergence rate of iteration for minimization of the cost function in 3DVar is greatly accelerated. The definition of dynamical balance depends on the characteristic scale of the circulation considered. The ratio of the balanced to the unbalanced parts is controlled by the prescribed statistics of background errors. Idealized trials produce the same results as the analytic solution. The results of real data case studies show the capability of the system to improve analysis compared to the traditional schemes. Finally, further development of the system is discussed.
基金Support by the National Natural Science Foundation of China under No.40175028 and the Key Technologies R & D Programme under No.2004BA607B.
文摘The dynamical constrains in three-dimensional variational data assimilation are discussed when considering the impact of stream divergence and convergence on the pressure and wind fields. For the analysis of severe tropical cyclone, frontal structures, and other rapidly changing structures, the geostrophic balance and linear balance cannot properly represent the relationship between wind and pressure fields. However, the nonlinear balance incremental equation takes into account the information of flow-dependent background, and makes response to the flow-dependent background covariance in the 3D-Var system. Results indicate that the application of the nonlinear balance equation to 3D-Var system improves the quality of severe tropical cyclone assimilation system, which has some positive effects on intensity prediction of tropical cyclones.
基金supported by the Natural Science Foundation of Zhejiang Province,China(Grant Nos.LY20A010021,LY19A010002,LY20G030025)the Natural Science Founda-tion of Ningbo City,China(Grant Nos.2021J147,2021J235).
文摘By introducing the dimensional splitting(DS)method into the multiscale interpolating element-free Galerkin(VMIEFG)method,a dimension-splitting multiscale interpolating element-free Galerkin(DS-VMIEFG)method is proposed for three-dimensional(3D)singular perturbed convection-diffusion(SPCD)problems.In the DSVMIEFG method,the 3D problem is decomposed into a series of 2D problems by the DS method,and the discrete equations on the 2D splitting surface are obtained by the VMIEFG method.The improved interpolation-type moving least squares(IIMLS)method is used to construct shape functions in the weak form and to combine 2D discrete equations into a global system of discrete equations for the three-dimensional SPCD problems.The solved numerical example verifies the effectiveness of the method in this paper for the 3D SPCD problems.The numerical solution will gradually converge to the analytical solution with the increase in the number of nodes.For extremely small singular diffusion coefficients,the numerical solution will avoid numerical oscillation and has high computational stability.
基金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 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.
文摘To address the problem of multi-missile cooperative interception against maneuvering targets at a prespecified impact time and desired Line-of-Sight(LOS)angles in ThreeDimensional(3D)space,this paper proposes a 3D leader-following cooperative interception guidance law.First,in the LOS direction of the leader,an impact time-controlled guidance law is derived based on the fixed-time stability theory,which enables the leader to complete the interception task at a prespecified impact time.Next,in the LOS direction of the followers,by introducing a time consensus tracking error function,a fixed-time consensus tracking guidance law is investigated to guarantee the consensus tracking convergence of the time-to-go.Then,in the direction normal to the LOS,by combining the designed global integral sliding mode surface and the second-order Sliding Mode Control(SMC)theory,an innovative 3D LOS-angle-constrained interception guidance law is developed,which eliminates the reaching phase in the traditional sliding mode guidance laws and effectively saves energy consumption.Moreover,it effectively suppresses the chattering phenomenon while avoiding the singularity issue,and compensates for unknown interference caused by target maneuvering online,making it convenient for practical engineering applications.Finally,theoretical proof analysis and multiple sets of numerical simulation results verify the effectiveness,superiority,and robustness of the investigated guidance law.
