Diagnostic Study of an Extreme Explosive Cyclone over the Kuroshio/Kuroshio Extension Region

Explosive cyclones(ECs)occur frequently over the Kuroshio/Kuroshio Extension region.The most rapidly intensified EC over the Kuroshio/Kuroshio Extension region during the 42 years(1979-2020)of cold seasons(October-April)was studied to reveal the variations of the key factors at different explosive-developing stages.This EC had weak low-level baroclinicity,mid-level cyclonic-vorticity advection,and strong low-level water vapor convergence at the initial explosive-developing stage.The low-level baroclinicity and mid-level cyclonic-vorticity advection increased substantially during the maximum-deepening-rate stage.The diagnostic analyses using the Zwack-Okossi equation showed that diabatic heating was the main contributor to the initial rapid intensification of this EC.The cyclonic-vorticity advection and warm-air advection enhanced rapidly in the middle and upper troposphere and contributed to the maximum rapid intensification,whereas the diabatic heating weakened slightly in the mid-low troposphere.The relative contribution of the diabatic heating decreased from the initial explosive-developing stage to the maximum-deepening-rate stage due to the enhancement of other factors(the cyclonic-vorticity advection and warm-air advection).Furthermore,the physical factors contributing to this EC varied with the explosive-developing stage.The non-key factors at the initial explosive-developing stage need attention to forecast the rapid intensification.

Seasonal Variation of the Sea Surface Temperature Growth Rate of ENSO

El Ni?o–Southern Oscillation(ENSO) exhibits a distinctive phase-locking characteristic, first expressed during its onset in boreal spring, developing during summer and autumn, reaching its peak towards winter, and decaying over the next spring. Several studies have demonstrated that this feature arises as a result of seasonal variation in the growth rate of ENSO as expressed by the sea surface temperature(SST). The bias towards simulating the phase locking of ENSO by many state-of-the-art climate models is also attributed to the unrealistic depiction of the growth rate. In this study, the seasonal variation of SST growth rate in the Ni?o-3.4 region(5°S–5°N, 120°–170°W) is estimated in detail based on the mixed layer heat budget equation and recharge oscillator model during 1981–2020. It is suggested that the consideration of a variable mixed layer depth is essential to its diagnostic process. The estimated growth rate has a remarkable seasonal cycle with minimum rates occurring in spring and maximum rates evident in autumn. More specifically, the growth rate derived from the meridional advection(surface heat flux) is positive(negative) throughout the year. Vertical diffusion generally makes a negative contribution to the evolution of growth rate and the magnitude of vertical entrainment represents the smallest contributor. Analysis indicates that the zonal advective feedback is regulated by the meridional immigration of the intertropical convergence zone, which approaches its southernmost extent in February and progresses to its northernmost location in September, and dominates the seasonal variation of the SST growth rate.

Application of a Shape-Preserving Advection Scheme to the Moisture Equation in an E-grid Regional Forecast Model

This paper presents a methodology which is very useful to design shape-preserving advection finite difference scheme on general E-grid horizontal arrangement of variables through introducing a two-step shape-preserving positive definite advection scheme in the moisture equation of the LASG-REM (LASG regional E-grid eta-coordinate forecast model). By trial-forecasting six local heavy raincases, the efficiency of the shape-preserving advection scheme in practical application has been examined. The LASG-REM with the shape-preserving advection scheme has a good forecasting ability for local precipitation.

Diagnosis of a Moist Thermodynamic Advection Parameter in Heavy-Rainfall Events

A moist thermodynamic advection parameter, defined as an absolute value of the dot product of hori- zontal gradients of three-dimensional potential temperature advection and general potential temperature, is introduced to diagnose frontal heavy rainfall events in the north of China. It is shown that the parameter is closely related to observed 6-h accumulative surface rainfall and simulated cloud hydrometeors. Since the parameter is capable of describing the typical vertical structural characteristics of dynamic, thermodynamic and water vapor fields above a strong precipitation region near the front surface, it may serve as a physical tracker to detect precipitable weather systems near to a front. A tendency equation of the parameter was derived in Cartesian coordinates and calculated with the simulation output data of a heavy rainfall event. Results revealed that the advection of the parameter by the three-dimensional velocity vector, the covariance of potential temperature advection by local change of the velocity vector and general potential temperature, and the interaction between potential temperature advection and the source or sink of general potential temperature, accounted for local change in the parameter. This indicated that the parameter was determined by a combination of dynamic processes and cloud microphysical processes.

A Two-Step Shape-Preserving Advection Scheme

This paper proposes a new two-step non-oscillatory shape-preserving positive definite finite difference advection transport scheme, which merges the advantages of small dispersion error in the simple first-order upstream scheme and small dissipation error in the simple second-order Lax-Wendroff scheme and is completely different from most of present positive definite advection schemes which are based on revising the upstream scheme results. The proposed scheme is much less time consuming than present shape-preserving or non-oscillatory advection transport schemes and produces results which are comparable to the results obtained from the present more complicated schemes. Elementary tests are also presented to examine the behavior of the scheme.

Improvement of the Semi-Lagrangian Advection Scheme in the GRAPES Model:Theoretical Analysis and Idealized Tests

ABSTRACT The Global/Regional Assimilation and PrEdiction System （GRAPES） is the newgeneration numerical weather predic- tion （NWP） system developed by the China Meteorological Administration. It is a fully compressible non-hydrostatical global/regional unified model that uses a traditional semi-Lagrangian advection scheme with cubic Lagrangian interpola tion （referred to as the SL_CL scheme）. The SL_CL scheme has been used in many operational NWP models, but there are still some deficiencies, such as the damping effects due to the interpolation and the relatively low accuracy. Based on Reich＇s semi-Lagrangian advection scheme （referred to as the R2007 scheme）, the Re_R2007 scheme that uses the low- and high-order B-spline function for interpolation at the departure point, is developed in this paper. One- and two-dimensional idealized tests in the rectangular coordinate system with uniform grid cells were conducted to compare the Re..R2007 scheme and the SL_CL scheme. The numerical results showed that： （1） the damping effects were remarkably reduced with the Re_R2007 scheme; and （2） the normalized errors of the Re_R2007 scheme were about 7.5 and 3 times smaller than those of the SL_CL scheme in one- and two-dimensional tests, respectively, indicating the higher accuracy of the Re..R2007 scheme. Furthermore, two solid-body rotation tests were conducted in the latitude-longitude spherical coordinate system with non uniform grid cells, which also verified the Re_R2007 scheme＇s advantages. Finally, in comparison with other global advection schemes, the Re_R2007 scheme was competitive in terms of accuracy and flow independence. An encouraging possibility for the application of the Re_R2007 scheme to the GRAPES model is provided.

Analytical Solution to the One-Dimensional Advection-Diffusion Equation with Temporally Dependent Coefficients

In a one-dimensional advection-diffusion equation with temporally dependent coefficients three cases may arise: solute dispersion parameter is time dependent while the flow domain transporting the solutes is uniform, the former is uniform and the latter is time dependent and lastly the both parameters are time dependent. In the present work analytical solutions are obtained for the last case, studying the dispersion of continuous input point sources of uniform and increasing nature in an initially solute free semi-infinite domain. The solutions for the first two cases and for uniform dispersion along uniform flow are derived as particular cases. The dispersion parameter is not proportional to the velocity of the flow. The Laplace transformation technique is used. New space and time variables are introduced to get the solutions. The solutions in all possible combinations of increasing/decreasing temporal dependence are compared with each other with the help of graphs. It has been observed that the concentration attenuation with position and time is the fastest in case of decreasing dispersion in accelerating flow field.

Role of Horizontal Density Advection in Seasonal Deepening of the Mixed Layer in the Subtropical Southeast Pacific

The mechanisms behind the seasonal deepening of the mixed layer （ML） in the subtropical Southeast Pacific were investigated using the monthly Argo data from 2004 to 2012. The region with a deep ML （more than 175 m） was found in the region of （22°-30°S, 105°-90°W）, reaching its maximum depth （-200 m） near （27°-28°S, 100°W） in September. The relative importance of horizontal density advection in determining the maximum ML location is discussed qualitatively. Downward Ekman pumping is key to determining the eastern boundary of the deep ML region. In addition, zonal density advection by the subtropical countercurrent （STCC） in the subtropical Southwest Pacific determines its western boundary, by carrying lighter water to strengthen the stratification and form a ＂shallow tongue＂ of ML depth to block the westward extension of the deep ML in the STCC region. The temperature advection by the STCC is the main source for large heat loss from the subtropical Southwest Pacific. Finally, the combined effect of net surface heat flux and meridional density advection by the subtropical gyre determines the northern and southern boundaries of the deep ML region： the ocean heat loss at the surface gradually increases from 22~S to 35~S, while the meridional density advection by the subtropical gyre strengthens the strat- ification south of the maximum ML depth and weakens the stratification to the north. The freshwater flux contribution to deepening the ML during austral winter is limited. The results are useful for understanding the role of ocean dynamics in the ML formation in the subtropical Southeast Pacific.

RETROSPECTIVE TIME INTEGRAL SCHEME AND ITS APPLICATIONS TO THE ADVECTION EQUATION

To put more information into a difference scheme of a differential equation for making an accurate prediction, a new kind of time integration scheme, known as the retrospective (RT) scheme, is proposed on the basis of the memorial dynamics. Stability criteria of the scheme for an advection equation in certain conditions are derived mathematically. The computations for the advection equation have been conducted with its RT scheme. It is shown that the accuracy of the scheme is much higher than that of the leapfrog (LF) difference scheme.

Energetics of lateral eddy diffusion/advection: Part III. Energetics of horizontal and isopycnal diffusion/ advection

Gravitational Potential Energy （GPE） change due to horizontal/isopycnal eddy diffusion and advection is examined. Horizontal/isopycnal eddy diffusion is conceptually separated into two steps： stirring and sub scale diffusion. GPE changes associated with these two steps are analyzed. In addition, GPE changes due to stirring and subscale diffusion associated with horizontal/isopycnal advection in the Eulerian coordinates are analyzed. These formulae are applied to the SODA data for the world oceans. Our analysis indicates that horizontal/isopycnal advection in Eulerian coordinates can introduce large artificial diffusion in the model. It is shown that GPE source/sink in isopycnal coordinates is closely linked to physical property distribution, such as temperature, salinity and velocity. In comparison with z-coordinates, GPE source/sink due to stir ring/cabbeling associated with isopycnal diffusion/advection is much smaller. Although isopycnal coordi nates may be a better choice in terms of handling lateral diffusion, advection terms in the traditional Eule rian coordinates can produce artificial source of GPE due to cabbeling associated with advection. Reducing such numerical errors remains a grand challenge.

Revisiting mesoscale eddy genesis mechanism of nonlinear advection in a marginal ice zone

A three-dimensional（3-D） ocean model is coupled with a two-dimensional（2-D） sea ice model, to revisit a nonlinear advection mechanism, one of the most important mesoscale eddy genesis mechanisms in the marginal ice zone. Two-dimensional ocean model simulations suggest nonlinear advection mechanism is more important when the water gets shallower. Instead of considering the ocean as barotropic fluid in the 2-D ocean model, the 3-D ocean model allows the sea ice to affect the current directly in the surface layer via ocean-ice interaction. It is found that both mesoscale eddy and sea surface elevation are sensitive to changes in a water depth in the 3-D simulations. The vertical profile of a current velocity in 3-D experiments suggests that when the water depth gets shallower, the current move faster in each layer, which makes the sea surface elevation be nearly inverse proportional to the water depth with the same wind forcing during the same time. It is also found that because of the vertical motion, the magnitude of variations in the sea surface elevation in the 3-D simulations is very small,being only 1% of the change in the 2-D simulations. And it seems the vertical motion to be the essential reason for the differences between the 3-D and 2-D experiments.

Study of Effects of Beta Term and Nonlinear Advection on the Structure of Tropical Cyclones

The effects of the Beta term on the typhoon structure are examined within the linear framework in terms of an analytical method of 2-D Fourier representation and numerical experiments by a Beta-plane quasi-geostrophic barotropic model. Results show that the joint effects of the difference of Rossby phase velocities and the dispersion of typhoon energy keep the maximum wind velocity reasonably evolving rather than irrestrictively increasing. On the one hand, the nonlinear advection accelerates typhoon vortex damping, and on the other, the high pressure system formed downstream due to energy dispersion makes it easy to maintain.

Energetics of lateral eddy diffusion/advection: Part IV. Energetics of diffusion/advection in sigma coordinates and other coordinates

Gravitational potential energy （GPE） source and sink due to stirring and cabbeling associated with sigma dif fusion/ advection is analyzed. It is shown that GPE source and sink is too big, and they are not closely linked to physical property distribution, such as temperature, salinity and velocity. Although the most frequently quoted advantage of sigma coordinate models are their capability of dealing with topography; the exces sive amount of GPE source and sink due to stirring and cabbeling associated with sigma diffusion/advec tion diagnosed from our analysis raises a very serious question whether the way lateral diffusion/advection simulated in the sigma coordinates model is physically acceptable. GPE source and sink in three coordinates is dramatically different in their magnitude and patterns. Overall, in terms of simulating lateral eddy diffu sion and advection isopycnal coordinates is the best choice and sigma coordinates is the worst. The physical reason of the excessive GPE source and sink in sigma coordinates is further explored in details. However, even in the isopycnal coordinates, simulation based on the Eulerian coordinates can be contaminated by the numerical errors associated with the advection terms.

Implementation of a Conservative Two-step Shape-Preserving Advection Scheme on a Spherical Icosahedral Hexagonal Geodesic Grid

An Eulerian flux-form advection scheme, called the Two-step Shape-Preserving Advection Scheme （TSPAS）, was generalized and implemented on a spherical icosahedral hexagonal grid （also referred to as a geodesic grid） to solve the transport equation. The C grid discretization was used for the spatial discretization. To implement TSPAS on an unstructured grid, the original finite-difference scheme was further generalized. The two-step integration utilizes a combination of two separate schemes （a low-order monotone scheme and a high-order scheme that typically cannot ensure monotonicity） to calculate the fluxes at the cell walls （one scheme corresponds to one cell wall）. The choice between these two schemes for each edge depends on a pre-updated scalar value using slightly increased fluxes. After the determination of an appropriate scheme, the final integration at a target cell is achieved by summing the fluxes that are computed by the different schemes. The conservative and shape-preserving properties of the generalized scheme are demonstrated. Numerical experiments are conducted at several horizontal resolutions. TSPAS is compared with the Flux Corrected Transport （FCT） approach to demonstrate the differences between the two methods, and several transport tests are performed to examine the accuracy, efficiency and robustness of the two schemes.

Consistency Problem with Tracer Advection in the Atmospheric Model GAMIL

The radon transport test, which is a widely used test case for atmospheric transport models, is carried out to evaluate the tracer advection schemes in the Grid-Point Atmospheric Model of IAP-LASG （GAMIL）. Two of the three available schemes in the model are found to be associated with significant biases in the polar regions and in the upper part of the atmosphere, which implies potentially large errors in the simulation of ozone-like tracers. Theoretical analyses show that inconsistency exists between the advection schemes and the discrete continuity equation in the dynamical core of GAMIL and consequently leads to spurious sources and sinks in the tracer transport equation. The impact of this type of inconsistency is demonstrated by idealized tests and identified as the cause of the aforementioned biases. Other potential effects of this inconsistency are also discussed. Results of this study provide some hints for choosing suitable advection schemes in the GAMIL model. At least for the polax-region-concentrated atmospheric components and the closely correlated chemical species, the Flux-Form Semi-Lagrangian advection scheme produces more reasonable simulations of the large-scale transport processes without significantly increasing the computational expense.

Numerical Analysis of Novel Coronavirus (2019-nCov) Pandemic Model with Advection

Recently,the world is facing the terror of the novel corona-virus,termed as COVID-19.Various health institutes and researchers are continuously striving to control this pandemic.In this article,the SEIAR(susceptible,exposed,infected,symptomatically infected,asymptomatically infected and recovered)infection model of COVID-19 with a constant rate of advection is studied for the disease propagation.A simple model of the disease is extended to an advection model by accommodating the advection process and some appropriate parameters in the system.The continuous model is transposed into a discrete numerical model by discretizing the domains,finitely.To analyze the disease dynamics,a structure preserving non-standard finite difference scheme is designed.Two steady states of the continuous system are described i.e.,virus free steady state and virus existing steady state.Graphical results show that both the steady states of the numerical design coincide with the fixed points of the continuous SEIAR model.Positivity of the state variables is ensured by applying the M-matrix theory.A result for the positivity property is established.For the proposed numerical design,two different types of the stability are investigated.Nonlinear stability and linear stability for the projected scheme is examined by applying some standard results.Von Neuman stability test is applied to ensure linear stability.The reproductive number is described and its pivotal role in stability analysis is also discussed.Consistency and convergence of the numerical model is also studied.Numerical graphs are presented via computer simulations to prove the worth and efficiency of the quarantine factor is explored graphically,which is helpful in controlling the disease dynamics.In the end,the conclusion of the study is also rendered.

Effects of advection,turbulence and radiation on formation of seafog──A theoretical analysis

Ｅｆｆｅｃｔｓｏｆａｄｖｅｃｔｉｏｎ，ｔｕｒｂｕｌｅｎｃｅａｎｄｒａｄｉａｔｉｏｎｏｎｆｏｒｍａｔｉｏｎｏｆｓｅａｆｏｇ──ＡｔｈｅｏｒｅｔｉｃａｌａｎａｌｙｓｉｓＫｅｙｗｏｒｄｓＲｅｌａｔｉｖｅｈｕｍｉｄｉｔｙｅｑｕａｔｉｏｎ，ａｄｖｅｃｔｉｏｎ，...

The effects of horizontal advection on the spring bloom of phytoplankton in the central Southern Huanghai Sea

An algal bloom is defined as a relatively rapid increase in the biomass of phytoplankton in an aquatic system. During 30 March to 24 April 2007, a cruise was conducted in the central Southern Huanghai Sea to investigate the spring bloom processes. The spatial and temporal variations of phytoplankton are discussed based on the in-situ observations and simultaneous remote sensing data. The explosive algal blooming varied quickly in temporal and spatial scales, due to the highly patchy distribution. Data obtained at the 2 anchor stations （BM1 and BM2） were analyzed in the present study. Horizontal advection is speculated to be responsible for the abrupt decrease in the concentration of chlorophyll-a at stations BM1 and BM2. At station BM2, the intermediate high chlorophyll-a concentration, coinciding with the low temperature, was found to be advected from the inshore colder water mass located to the east of the site.

Numerical modeling of contaminant advection impact on hydrodynamic diffusion in a deformable medium

The self-weight of solid waste or machine-rolled compaction can induce or trigger contaminant migration in the landfill.Although the consolidation-induced hydraulic gradient driving solution transport has been extensively investigated,little attention has been paid to ion migration caused by its concentration gradient variation.It is necessary to more precisely predict the multi-stage contaminant transports in deforming porous material.Based on the modified Cam-clay model,the proposed fluid-solid coupled model can simulate the elastoplastic deformation behavior of layered kaolinite and KBr solution transport/sorption,and its modeling results were validated by published laboratory data.The solid-fluid interactions were analyzed by comparing various transport manners of K^(+)and Br^(−)from excess pore pressure generation to dissipation.Results reveal that the consolidation process can accelerate KBr solute advection from the contaminated layer into the uncontaminated layer,and then affects the subsequent diffusion,mechanical dispersion and sorption for K^(+)and Br^(−).The simulations also indicate that consolidation-induced solute transport is time-dependent,and therefore the ion diffusion and mechanical dispersion should receive more attention.

The Influence of Weakly-Nonlinear Vertical Advection on the Wind Field of PBL with Large-Scale Orography

The analytical solutions of the PBL wind distribution under the equilibrium of four forces including both horizontal and vertical advections are obtained in this paper using small parameter method. Utilizing this simple PBL model, we also compute the wind spirals of a certain circular vortex over level ground and over a west-east ranged mountain, respectively. The results show that when there exsits large-scale orographic forcing,vertical advection may exert considerable influence on the wind distribution of PBL, the angle between geostrophic flow and surface wind, and the height of the top of PBL.