The version 2.1 of PETROMAR-3D model, created in the Center for Marine Meteorology of the Meteorology Institute of Cuba, is presented. This Lagrangian model has been designed to describe the physical processes of mari...The version 2.1 of PETROMAR-3D model, created in the Center for Marine Meteorology of the Meteorology Institute of Cuba, is presented. This Lagrangian model has been designed to describe the physical processes of marine oil spills in the face of multiple scenarios of the marine environment. Although it is applicable to any part of the world, it is mainly designed for inter-American seas. The novelty has been to integrate the processes of drift and weathering into a model, with updated methods that incorporate 3D phenomena, a very favorable situation to achieve an operating system in Cuba and the region for the immediate and medium term. Python was chosen as the programming language because it has advanced libraries for numerical modeling, automation work and other useful tools for pre-and post-processing. By means of adapters, an important number of atmospheric, hydrodynamic and wave models have been considered to create the scenarios efficiently. The modular distribution in which the code has been created facilitates its use for other dispersion analysis and biophysical applications. Finally, a set of simple images are presented, aimed at informing decision-makers in order to mitigate the effects of the spill on the environment.展开更多
The Lagrangian equations of motion of small particle in turbulent boundary layer flows, taking into account the effects of the drug force caused by the wall presence, the Saffman and the Magus lift forces et al., is s...The Lagrangian equations of motion of small particle in turbulent boundary layer flows, taking into account the effects of the drug force caused by the wall presence, the Saffman and the Magus lift forces et al., is studied. Using the spectral method, analytical expression relating to the Lagrangian power spectra of particle velocity to that of the fluid are developed and the results are used to evaluate various responses statistics. In this paper, the results clearly show that the turbulent diffusivity of the particle may be larger than that of the fluid for a period of long-time.展开更多
Spatial heterogeneity or“patchiness”of plankton distributions in the ocean has always been an attractive and challenging scientific issue to oceanographers.We focused on the accumulation and dynamic mechanism of the...Spatial heterogeneity or“patchiness”of plankton distributions in the ocean has always been an attractive and challenging scientific issue to oceanographers.We focused on the accumulation and dynamic mechanism of the Acetes chinensis in the Lianyungang nearshore licensed fishing area.The Lagrangian frame approaches including the Lagrangian coherent structures theory,Lagrangian residual current,and Lagrangian particle-tracking model were applied to find the transport pathways and aggregation characteristics of Acetes chinensis.There exist some material transport pathways for Acetes chinensis passing through the licensed fishing area,and Acetes chinensis is easy to accumulate in the licensed fishing area.The main mechanism forming this distribution pattern is the local circulation induced by the nonlinear interaction of topography and tidal flow.Both the Lagrangian coherent structure analysis and the particle trajectory tracking indicate that Acetes chinensis in the licensed fishing area come from the nearshore estuary.This work contributed to the adjustment of licensed fishing area and the efficient utilization of fishery resources.展开更多
A hybrid Lagrangian - Eulerian （HLE） method is developed for sea ice dynamics, which combines the high computational efficiency of finite difference method （FDM） with the high numerical accuracy of smoothed partic...A hybrid Lagrangian - Eulerian （HLE） method is developed for sea ice dynamics, which combines the high computational efficiency of finite difference method （FDM） with the high numerical accuracy of smoothed particle hydrodynamics （SPH）. In this HLE model, the sea ice cover is represented by a group of Lagrangian ice particles with their own thicknesses and concentrations. These ice variables are interpolated to the Eularian gird nodes using the Gaussian interpolation function. The FDM is used to determine the ice velocities at Eulerian grid nodes, and the velocities of Lagrangian ice particles are interpolated from these grid velocities with the Gaussian function also. The thicknesses and concentrations of ice particles are determined based on their new locations. With the HLE numerical model, the ice ridging process in a rectangular basin is simulated, and the simulated results are validated with the analytical solution. This method is also applied to the simulation of sea ice dynamics in a vortex wind field. At last, this HLE model is applied to the Bohai Sea, and the simulated concentration, thickness and velocity match the satellite images and the field observed data well.展开更多
Wind field simulation in the surface layer is often used to manage natural resources in terms of air quality,gene flow(through pollen drift),and plant disease transmission(spore dispersion).Although Lagrangian stochas...Wind field simulation in the surface layer is often used to manage natural resources in terms of air quality,gene flow(through pollen drift),and plant disease transmission(spore dispersion).Although Lagrangian stochastic(LS)models describe stochastic wind behaviors,such models assume that wind velocities follow Gaussian distributions.However,measured surface-layer wind velocities show a strong skewness and kurtosis.This paper presents an improved model,a non-Gaussian LS model,which incorporates controllable non-Gaussian random variables to simulate the targeted non-Gaussian velocity distribution with more accurate skewness and kurtosis.Wind velocity statistics generated by the non-Gaussian model are evaluated by using the field data from the Cooperative Atmospheric Surface Exchange Study,October 1999 experimental dataset and comparing the data with statistics from the original Gaussian model.Results show that the non-Gaussian model improves the wind trajectory simulation by stably producing precise skewness and kurtosis in simulated wind velocities without sacrificing other features of the traditional Gaussian LS model,such as the accuracy in the mean and variance of simulated velocities.This improvement also leads to better accuracy in friction velocity(i.e.,a coupling of three-dimensional velocities).The model can also accommodate various non-Gaussian wind fields and a wide range of skewness–kurtosis combinations.Moreover,improved skewness and kurtosis in the simulated velocity will result in a significantly different dispersion for wind/particle simulations.Thus,the non-Gaussian model is worth applying to wind field simulation in the surface layer.展开更多
In this paper, a 10-level Quasi-Lagrangian Limited Area Model is used to simulate the process of Meiyu(Baiu) front of 1979. Some physical processes, such as large-scale condensation and cumulus convection, are include...In this paper, a 10-level Quasi-Lagrangian Limited Area Model is used to simulate the process of Meiyu(Baiu) front of 1979. Some physical processes, such as large-scale condensation and cumulus convection, are included in the model. The simulation results are encouraging. 24-h numerical simulation shows that the invading of cold air from North China and rapidly northward moving of warm air from South China can be successfully reproduced. The terrain with a maximum of 4175 m is incorporated in the model. Three different kinds of terrain schemes are tested and the dynamic effect of the Plateau on the process of heavy rainfall is found to be very important.展开更多
The spatial distribution and seasonal variation of the tide-induced Lagrangian Residual Circulations （LRC hereafter）, wind-driven LRC, and the coupling dynamic characteristics were simulated using ECOM, given the He...The spatial distribution and seasonal variation of the tide-induced Lagrangian Residual Circulations （LRC hereafter）, wind-driven LRC, and the coupling dynamic characteristics were simulated using ECOM, given the Hellerman and Rosenstein global monthly-mean wind stresses. The results showed that the tide-induced LRC of the harmonic constituent M2 bears an identical pattern in four seasons in the Bohai Sea： the surface one is weak with random directions; however, there exist a southeast current from the Bohai Strait to the Laizhou bay, and a weakly anticlockwise gyre in the south of the Bohai Strait for the bottom layer LRC. The magnitude of bottom layer tide-induced LRC is larger than the surface one, and moreover, it contributes significantly to the whole LRC in the Bohai Sea. Unlike the identical structure of the tide-induced LRC, the wind driven LRC varies seasonally under the prevailing monsoon. It forms a distinct gyre under the summer and winter monsoons in July and January respectively, but it seems weak and non-directional in April and September.展开更多
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 compre...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.展开更多
A flux-form semi-Lagrangian transport scheme （FFSL） was implemented in a spectral atmospheric GCM developed and used at IAP/LASG. Idealized numerical experiments show that the scheme is good at shape preserving with...A flux-form semi-Lagrangian transport scheme （FFSL） was implemented in a spectral atmospheric GCM developed and used at IAP/LASG. Idealized numerical experiments show that the scheme is good at shape preserving with less dissipation and dispersion, in comparison with other conventional schemes, hnportantly, FFSL can automatically maintain the positive definition of the transported tracers, which was an underlying problem in the previous spectral composite method （SCM）. To comprehensively investigate the impact of FFSL on GCM results, we conducted sensitive experiments. Three main improvements resulted： first, rainfall simulation in both distribution and intensity was notably improved, which led to an improvement in precipitation frequency. Second, the dry bias in the lower troposphere was significantly reduced compared with SCM simulations. Third, according to the Taylor diagram, the FFSL scheme yields simulations that are superior to those using the SCM： a higher correlation between model output and observation data was achieved with the FFSL scheme, especially for humidity in lower troposphere. However, the moist bias in the middle and upper troposphere was more pronounced with the FFSL scheme. This bias led to an over-simulation of precipitable water in comparison with reanalysis data. Possible explanations, as well as solutions, are discussed herein.展开更多
A scale-similarity model of a two-point two-time Lagrangian velocity correlation(LVC) was originally developed for the relative dispersion of tracer particles in isotropic turbulent flows(HE, G. W., JIN, G. D., and ZH...A scale-similarity model of a two-point two-time Lagrangian velocity correlation(LVC) was originally developed for the relative dispersion of tracer particles in isotropic turbulent flows(HE, G. W., JIN, G. D., and ZHAO, X. Scale-similarity model for Lagrangian velocity correlations in isotropic and stationary turbulence. Physical Review E, 80, 066313(2009)). The model can be expressed as a two-point Eulerian space correlation and the dispersion velocity V. The dispersion velocity denotes the rate at which one moving particle departs from another fixed particle. This paper numerically validates the robustness of the scale-similarity model at high Taylor micro-scale Reynolds numbers up to 373, which are much higher than the original values(R_λ = 66, 102). The effect of the Reynolds number on the dispersion velocity in the scale-similarity model is carefully investigated. The results show that the scale-similarity model is more accurate at higher Reynolds numbers because the two-point Lagrangian velocity correlations with different initial spatial separations collapse into a universal form compared with a combination of the initial separation and the temporal separation via the dispersion velocity.Moreover, the dispersion velocity V normalized by the Kolmogorov velocity V_η ≡ η/τ_η in which η and τ_η are the Kolmogorov space and time scales, respectively, scales with the Reynolds number R_λ as V/V_η ∝ R_λ^(1.39) obtained from the numerical data.展开更多
文摘The version 2.1 of PETROMAR-3D model, created in the Center for Marine Meteorology of the Meteorology Institute of Cuba, is presented. This Lagrangian model has been designed to describe the physical processes of marine oil spills in the face of multiple scenarios of the marine environment. Although it is applicable to any part of the world, it is mainly designed for inter-American seas. The novelty has been to integrate the processes of drift and weathering into a model, with updated methods that incorporate 3D phenomena, a very favorable situation to achieve an operating system in Cuba and the region for the immediate and medium term. Python was chosen as the programming language because it has advanced libraries for numerical modeling, automation work and other useful tools for pre-and post-processing. By means of adapters, an important number of atmospheric, hydrodynamic and wave models have been considered to create the scenarios efficiently. The modular distribution in which the code has been created facilitates its use for other dispersion analysis and biophysical applications. Finally, a set of simple images are presented, aimed at informing decision-makers in order to mitigate the effects of the spill on the environment.
文摘The Lagrangian equations of motion of small particle in turbulent boundary layer flows, taking into account the effects of the drug force caused by the wall presence, the Saffman and the Magus lift forces et al., is studied. Using the spectral method, analytical expression relating to the Lagrangian power spectra of particle velocity to that of the fluid are developed and the results are used to evaluate various responses statistics. In this paper, the results clearly show that the turbulent diffusivity of the particle may be larger than that of the fluid for a period of long-time.
基金the National Natural Science Foundation of China(No.31802297)。
文摘Spatial heterogeneity or“patchiness”of plankton distributions in the ocean has always been an attractive and challenging scientific issue to oceanographers.We focused on the accumulation and dynamic mechanism of the Acetes chinensis in the Lianyungang nearshore licensed fishing area.The Lagrangian frame approaches including the Lagrangian coherent structures theory,Lagrangian residual current,and Lagrangian particle-tracking model were applied to find the transport pathways and aggregation characteristics of Acetes chinensis.There exist some material transport pathways for Acetes chinensis passing through the licensed fishing area,and Acetes chinensis is easy to accumulate in the licensed fishing area.The main mechanism forming this distribution pattern is the local circulation induced by the nonlinear interaction of topography and tidal flow.Both the Lagrangian coherent structure analysis and the particle trajectory tracking indicate that Acetes chinensis in the licensed fishing area come from the nearshore estuary.This work contributed to the adjustment of licensed fishing area and the efficient utilization of fishery resources.
基金The study was supported by the National Natural Science Foundation of China under contract No.10772041the State 0ceamic Administration Key Laboratory for Ploar Science of China under contract No.KP 2007004.
文摘A hybrid Lagrangian - Eulerian （HLE） method is developed for sea ice dynamics, which combines the high computational efficiency of finite difference method （FDM） with the high numerical accuracy of smoothed particle hydrodynamics （SPH）. In this HLE model, the sea ice cover is represented by a group of Lagrangian ice particles with their own thicknesses and concentrations. These ice variables are interpolated to the Eularian gird nodes using the Gaussian interpolation function. The FDM is used to determine the ice velocities at Eulerian grid nodes, and the velocities of Lagrangian ice particles are interpolated from these grid velocities with the Gaussian function also. The thicknesses and concentrations of ice particles are determined based on their new locations. With the HLE numerical model, the ice ridging process in a rectangular basin is simulated, and the simulated results are validated with the analytical solution. This method is also applied to the simulation of sea ice dynamics in a vortex wind field. At last, this HLE model is applied to the Bohai Sea, and the simulated concentration, thickness and velocity match the satellite images and the field observed data well.
基金financial support for this research from a USDA-AFRI Foundational Grant (Grant No. 2012-67013-19687)from the Illinois State Water Survey at the University of Illinois at Urbana—Champaign
文摘Wind field simulation in the surface layer is often used to manage natural resources in terms of air quality,gene flow(through pollen drift),and plant disease transmission(spore dispersion).Although Lagrangian stochastic(LS)models describe stochastic wind behaviors,such models assume that wind velocities follow Gaussian distributions.However,measured surface-layer wind velocities show a strong skewness and kurtosis.This paper presents an improved model,a non-Gaussian LS model,which incorporates controllable non-Gaussian random variables to simulate the targeted non-Gaussian velocity distribution with more accurate skewness and kurtosis.Wind velocity statistics generated by the non-Gaussian model are evaluated by using the field data from the Cooperative Atmospheric Surface Exchange Study,October 1999 experimental dataset and comparing the data with statistics from the original Gaussian model.Results show that the non-Gaussian model improves the wind trajectory simulation by stably producing precise skewness and kurtosis in simulated wind velocities without sacrificing other features of the traditional Gaussian LS model,such as the accuracy in the mean and variance of simulated velocities.This improvement also leads to better accuracy in friction velocity(i.e.,a coupling of three-dimensional velocities).The model can also accommodate various non-Gaussian wind fields and a wide range of skewness–kurtosis combinations.Moreover,improved skewness and kurtosis in the simulated velocity will result in a significantly different dispersion for wind/particle simulations.Thus,the non-Gaussian model is worth applying to wind field simulation in the surface layer.
基金This research work is financially supported by the Eighth National Five-Year Scientific Project(85-906-04-02)and Nation Natural Founndation of China(4927520).
文摘In this paper, a 10-level Quasi-Lagrangian Limited Area Model is used to simulate the process of Meiyu(Baiu) front of 1979. Some physical processes, such as large-scale condensation and cumulus convection, are included in the model. The simulation results are encouraging. 24-h numerical simulation shows that the invading of cold air from North China and rapidly northward moving of warm air from South China can be successfully reproduced. The terrain with a maximum of 4175 m is incorporated in the model. Three different kinds of terrain schemes are tested and the dynamic effect of the Plateau on the process of heavy rainfall is found to be very important.
基金N ational N atural Science Foundation of China, N o.40271020 K now ledge Innovation Project of the Institute ofG eographic Sciences and N aturalResources Research,CA S,N o.CX IO G -A 04-09
文摘The spatial distribution and seasonal variation of the tide-induced Lagrangian Residual Circulations （LRC hereafter）, wind-driven LRC, and the coupling dynamic characteristics were simulated using ECOM, given the Hellerman and Rosenstein global monthly-mean wind stresses. The results showed that the tide-induced LRC of the harmonic constituent M2 bears an identical pattern in four seasons in the Bohai Sea： the surface one is weak with random directions; however, there exist a southeast current from the Bohai Strait to the Laizhou bay, and a weakly anticlockwise gyre in the south of the Bohai Strait for the bottom layer LRC. The magnitude of bottom layer tide-induced LRC is larger than the surface one, and moreover, it contributes significantly to the whole LRC in the Bohai Sea. Unlike the identical structure of the tide-induced LRC, the wind driven LRC varies seasonally under the prevailing monsoon. It forms a distinct gyre under the summer and winter monsoons in July and January respectively, but it seems weak and non-directional in April and September.
基金jointly sponsored by the Key Project of the Chinese National Programs for Fundamental Research and Development ("973 Program" Grant No.2013CB430106)+1 种基金the Key Project of the Chinese National Science & Technology Pillar Program during the Twelfth Five-year Plan Period (Grant No.2012BAC22B01)the National Natural Science Foundation of China ( Grant No.41375108)
文摘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.
基金supported by the Chinese Academy of Science Strategic Priority Research Program (Grant No. XDA05110303)"973" Program (Grant Nos. 2010CB950403,2012CB417203,and 2013CB955803)+1 种基金"863" Program(Grant No. 2010AA012305)the National Natural Science Foundation of China (Grant Nos. 40925015,40875034,and 41023002)
文摘A flux-form semi-Lagrangian transport scheme （FFSL） was implemented in a spectral atmospheric GCM developed and used at IAP/LASG. Idealized numerical experiments show that the scheme is good at shape preserving with less dissipation and dispersion, in comparison with other conventional schemes, hnportantly, FFSL can automatically maintain the positive definition of the transported tracers, which was an underlying problem in the previous spectral composite method （SCM）. To comprehensively investigate the impact of FFSL on GCM results, we conducted sensitive experiments. Three main improvements resulted： first, rainfall simulation in both distribution and intensity was notably improved, which led to an improvement in precipitation frequency. Second, the dry bias in the lower troposphere was significantly reduced compared with SCM simulations. Third, according to the Taylor diagram, the FFSL scheme yields simulations that are superior to those using the SCM： a higher correlation between model output and observation data was achieved with the FFSL scheme, especially for humidity in lower troposphere. However, the moist bias in the middle and upper troposphere was more pronounced with the FFSL scheme. This bias led to an over-simulation of precipitable water in comparison with reanalysis data. Possible explanations, as well as solutions, are discussed herein.
基金Project supported by the Science Challenge Program(No.TZ2016001)the National Natural Science Foundation of China(Nos.11472277,11572331,11232011,and 11772337)+1 种基金the Strategic Priority Research Program,Chinese Academy of Sciences(No.XDB22040104)the Key Research Program of Frontier Sciences,Chinese Academy of Sciences(No.QYZDJ-SSW-SYS002)
文摘A scale-similarity model of a two-point two-time Lagrangian velocity correlation(LVC) was originally developed for the relative dispersion of tracer particles in isotropic turbulent flows(HE, G. W., JIN, G. D., and ZHAO, X. Scale-similarity model for Lagrangian velocity correlations in isotropic and stationary turbulence. Physical Review E, 80, 066313(2009)). The model can be expressed as a two-point Eulerian space correlation and the dispersion velocity V. The dispersion velocity denotes the rate at which one moving particle departs from another fixed particle. This paper numerically validates the robustness of the scale-similarity model at high Taylor micro-scale Reynolds numbers up to 373, which are much higher than the original values(R_λ = 66, 102). The effect of the Reynolds number on the dispersion velocity in the scale-similarity model is carefully investigated. The results show that the scale-similarity model is more accurate at higher Reynolds numbers because the two-point Lagrangian velocity correlations with different initial spatial separations collapse into a universal form compared with a combination of the initial separation and the temporal separation via the dispersion velocity.Moreover, the dispersion velocity V normalized by the Kolmogorov velocity V_η ≡ η/τ_η in which η and τ_η are the Kolmogorov space and time scales, respectively, scales with the Reynolds number R_λ as V/V_η ∝ R_λ^(1.39) obtained from the numerical data.
