A theoretical study of the multigrid three-dimensional variational data assimilation scheme using a simple bilinear interpolation algorithm

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.

Hybrid three-dimensional variation and particle filtering for nonlinear systems

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.

Application of a Recursive Filter to a Three-Dimensional Variational Ocean Data Assimilation System

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.

Three-dimensional limit variation analysis on the ultimate pullout capacity of the anchor cables based on the Hoek-Brown failure criterion

Only simplified two-dimensional model and a single failure mode are adopted to calculate the ultimate pullout capacity(UPC)of anchor cables in most previous research.This study focuses on a more comprehensive combination failure mode that consists of bond failure of an anchorage body and failure of an anchored rock mass.The three-dimensional ultimate pullout capacity of the anchor cables is calculated based on the Hoek-Brown failure criterion and variation analysis method.The numerical solution for the curvilinear function in fracture plane is obtained based on the finite difference theory,which more accurately reflects the failure state of the anchor cable,as opposed to that being assumed in advance.The results reveal that relying solely on a single failure mode for UPC calculations has limitations,as changes in parameter values not only directly impact the UPC value but also can alter the failure model and thus the calculation method.

Construction and application of a three-dimensional vascular variation-based nephrometry scoring system for completely endophytic renal tumors

Background:Completely endophytic renal tumors(CERT)pose significant challenges due to their anatomical complexity and loss of visual clues about tumor location.A facile scoring model based on three-dimensional(3D)reconstructed images will assist in better assessing tumor location and vascular variations.Methods:In this retrospective study,80 patients diagnosed with CERT were included.Forty cases underwent preoperative assessment using 3D reconstructed imaging(3D-Cohort),while the remaining 40 cases were assessed using two-dimensional imaging(2D-Cohort).Vascular variations were evaluated by ascertaining the presence of renal arteries>1,prehilar branching arteries,and arteries anterior to veins.The proposed scoring system,termed RAL,encompassed three critical components:(R)adius(maximal tumor diameter in cm),(A)rtery(occurrence of arterial variations),and(L)ocation relative to the polar line.Comparison of the RAL scoring system was made with established nephrometry scoring systems.Results:A total of 48(60%)patients exhibited at least one vascular variation.In the 2D-Cohort,patients with vascular variations experienced significantly prolonged operation time,increased bleeding volume,and extended warm ischemia time compared with those without vascular variations.Conversely,the presence of vascular vari-ations did not significantly affect operative parameters in the 3D-Cohort.Furthermore,the 2D-Cohort demon-strated a notable decline in both short-and long-term estimated glomerular filtration rate(eGFR)changes com-pared with the 3D-Cohort,a trend consistent across patients with warm ischemia time≥25 min and those with vascular variations.Notably,the 2D-Cohort exhibited a larger margin of normal renal tissue compared with the 3D-Cohort.Elevated RAL scores correlated with larger tumor size,prolonged operation time,extended warm is-chemia time,and substantial postoperative eGFR decrease.The RAL scoring system displayed superior predictive capabilities in assessing postoperative eGFR changes compared with conventional nephrometry scoring systems.Conclusions:Our proposed 3D vascular variation-based nephrometry scoring system offers heightened proficiency in preoperative assessment,precise prediction of surgical complexity,and more accurate evaluation of postoper-ative renal function in CERT patients.

Assimilation of Total Lightning Data Using the Three-Dimensional Variational Method at Convection-Allowing Resolution

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.

Three-dimensional models:from cell culture to Patient-Derived Organoid and its application to future liposarcoma research

Liposarcoma is one of the most common soft tissue sarcomas,however,its occurrence rate is still rare compared to other cancers.Due to its rarity,in vitro experiments are an essential approach to elucidate liposarcoma pathobiology.Conventional cell culture-based research(2D cell culture)is still playing a pivotal role,while several shortcomings have been recently under discussion.In vivo,mouse models are usually adopted for pre-clinical analyses with expectations to overcome the issues of 2D cell culture.However,they do not fully recapitulate human dedifferentiated liposarcoma(DDLPS)characteristics.Therefore,three-dimensional(3D)culture systems have been the recent research focus in the cell biology field with the expectation to overcome at the same time the disadvantages of 2D cell culture and in vivo animal models and fill in the gap between them.Given the liposarcoma rarity,we believe that 3D cell culture techniques,including 3D cell cultures/co-cultures,and Patient-Derived tumor Organoids(PDOs),represent a promising approach to facilitate liposarcoma investigation and elucidate its molecular mechanisms and effective therapy development.In this review,we first provide a general overview of 3D cell cultures compared to 2D cell cultures.We then focus on one of the recent 3D cell culture applications,Patient-Derived Organoids(PDOs),summarizing and discussing several PDO methodologies.Finally,we discuss the current and future applications of PDOs to sarcoma,particularly in the field of liposarcoma.

Role of variation coefficient of stone density in determining success of shock wave lithotripsy in urinary calculi

BACKGROUND Various stone factors can affect the net results of shock wave lithotripsy(SWL).Recently a new factor called variation coefficient of stone density(VCSD)is being considered to have an impact on stone free rates.AIM To assess the role of VCSD in determining success of SWL in urinary calculi.METHODS Charts review was utilized for collection of data variables.The patients were subjected to SWL,using an electromagnetic lithotripter.Mean stone density(MSD),stone heterogeneity index(SHI),and VCSD were calculated by generating regions of interest on computed tomography(CT)images.Role of these factors were determined by applying the relevant statistical tests for continuous and categorical variables and a P value of<0.05 was gauged to be statistically significant.RESULTS There were a total of 407 patients included in the analysis.The mean age of the subjects in this study was 38.89±14.61 years.In total,165 out of the 407 patients could not achieve stone free status.The successful group had a significantly lower stone volume as compared to the unsuccessful group(P<0.0001).Skin to stone distance was not dissimilar among the two groups(P=0.47).MSD was significantly lower in the successful group(P<0.0001).SHI and VCSD were both significantly higher in the successful group(P<0.0001).CONCLUSION VCSD,a useful CT based parameter,can be utilized to gauge stone fragility and hence the prediction of SWL outcomes.

The three-dimensional structure and seasonal variation of the North Pacific meridional overturning circulation

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.

Observed spatiotemporal variation of three-dimensional structure and heat/salt transport of anticyclonic mesoscale eddy in Northwest Pacific

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.

Three-Dimensional Wind Field Retrieved from Dual-Doppler Radar Based on a Variational Method:Refinement of Vertical Velocity Estimates

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.

Limit variation analysis of shallow rectangular tunnels collapsing with double-layer rock mass based on a three-dimensional failure mechanism

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.

Three-Dimensional Models for Analyzing the Cyclic Variations in a Lean Burn CNG Engine

Three-dimensional models, consisting of the flame kernel formation model, flame kernel development model and natural gas single step reaction model, are used to analyze the contribution of cyclic equivalence ratio variations to cyclic variations in the compressed natural gas （CNG） lean burn spark ignition engine. Computational results including the contributions of equivalence ratio cyclic variations to each combustion stage and effects of engine speed to the extent of combustion variations are discussed. It is concluded that the equivalence ratio variations affect mostly the main stage of combustion and hardly influence initial kernel development stage.

Scientific design and preliminary results of three-dimensional variational data assimilation system of GRAPES

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.

Application of Nonlinear Balance Scheme in Three-Dimensional Variational Data Assimilation

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.

Modeling of Geometric Variations Within Three-Dimensional Tolerance Zones

A tolerance modeling method of geometric variations within three-dimensional(3D) tolerance domain is proposed.According to the classification method of feature variations and classes of invariance or symmetry in the new generation Geometrical Product Specifications(GPS) system,the proposed method is based on the Small Displacement Torsor(SDT) concept,and SDT is used to represent the geometric deviations between the nominal feature and the fitting feature. Furthermore,the variation zones of several kinds complex or irregular geometric features are constructed,and the corresponding tolerance modeling methods are discussed respectively. Finally,an example is presented to verify the efficiency of the presented method for modeling 3D dimensional and geometrical tolerances.The proposed modeling method based on the complete mathematical foundation conforms with the design philosophy and principle of the new generation GPS standards system.

Prediction of Salinity Variations in a Tidal Estuary Using Artificial Neural Network and Three-Dimensional Hydrodynamic Models

The simulation of salinity at different locations of a tidal river using physically-based hydrodynamic models is quite cumbersome because it requires many types of data, such as hydrological and hydraulic time series at boundaries, river geometry, and adjusted coefficients. Therefore, an artificial neural network (ANN) technique using a back-propagation neural network (BPNN) and a radial basis function neural network (RBFNN) is adopted as an effective alternative in salinity simulation studies. The present study focuses on comparing the performance of BPNN, RBFNN, and three-dimensional hydrodynamic models as applied to a tidal estuarine system. The observed salinity data sets collected from 18 to 22 May, 16 to 22 October, and 26 to 30 October 2002 (totaling 4320 data points) were used for BPNN and RBFNN model training and for hydrodynamic model calibration. The data sets collected from 30 May to 2 June and 11 to 15 November 2002 (totaling 2592 data points) were adopted for BPNN and RBFNN model verification and for hydrodynamic model verification. The results revealed that the ANN (BPNN and RBFNN) models were capable of predicting the nonlinear time series behavior of salinity to the multiple forcing signals of water stages at different stations and freshwater input at upstream boundaries. The salinity predicted by the ANN models was better than that predicted by the physically based hydrodynamic model. This study suggests that BPNN and RBFNN models are easy-to-use modeling tools for simulating the salinity variation in a tidal estuarine system.

Diagnosis of the Medium-Range Variation of the Subtropical High over the Western Pacific during a Meiyu Process by Three-Dimensional E-P Flux

In this paper, using the daily grid data (2.5 × 2.5) of the ECMWF / WMO, we have computed respectively the three-dimensional wave activity flux in the stages of pre-onset, prevailing and post ending of Meiyu from 1 to 31 July 1982. The potential vorticity field is taken as the physical quantity relating the wave activity flux to the variation of the subtropical high over the Western Pacific. It is found that the three-dimensional wave activity flux is a powerful means for diagnosis of the variation of the subtropical high over the Western Pacific: The region of the subtropical high is just the confluence area of wave energy, whose changes in intensity and range decide the variation of the subtropical high. The confluence of wave energy comes from the monsoon flow in low latitudes, the Meiyu rain belts in middle latitudes and the heating fields on the eastern side of the Qinghai-Xizang Plateau. The relation between these sources and the subtropical high displays the self-adjusting mechanism among members of East-Asia summer monsoon.

A Dimension-Splitting Variational Multiscale Element-Free Galerkin Method for Three-Dimensional Singularly Perturbed Convection-Diffusion Problems

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.

A theory for three-dimensional response of micropolar plates

Through combined applications of the transfer-matrix method and asymptotic expansion technique,we formulate a theory to predict the three-dimensional response of micropolar plates.No ad hoc assumptions regarding through-thickness assumptions of the field variables are made,and the governing equations are two-dimensional,with the displacements and microrotations of the mid-plane as the unknowns.Once the deformation of the mid-plane is solved,a three-dimensional micropolar elastic field within the plate is generated,which is exact up to the second order except in the boundary region close to the plate edge.As an illustrative example,the bending of a clamped infinitely long plate caused by a uniformly distributed transverse force is analyzed and discussed in detail.