准地转大洋风生环流的格子Boltzmann数值模拟

建立了求解准地转相当正压涡度方程的格子Boltzmann (LB)模型. 该模型将准地转相当正压涡度方程作为一个平流-扩散-反应方程来加以处理,在整体二阶精度下,通过Chapman_Enskog多尺度分析法,可将格子Boltzmann方程还原到相当正压涡度方程. 在不同Reynolds数、不同边界条件以及不同风应力驱动下的数值解表明,该模型正确反映了风生环流的基本结构和不同边界的耗散特征,并得到风生环流的多平衡态解等非线性特征. 此外,不同Rossby变形半径下的实验证明,小Rossby变形半径更容易激发环流的非线性模态. 通过与同等类型有限差方案的比较,表明本文的LB模型具有稳定性好、精度高等优点.     

Forecasting search areas using ensemble ocean circulation modeling

We investigate trajectory forecasting as an application of ocean circulation ensemble modeling. The ensemble simulations are performed weekly, starting with assimilation of data for various variables from multiple sensors on a range of observational platforms. The ensemble is constructed from 100 members, and member no. 1 is designed as a standard (deterministic) simulation, providing us with a benchmark for the study. We demonstrate the value of the ensemble approach by validating simulated trajectories using data from ocean surface drifting buoys. We find that the ensemble average trajectories are generally closer to the observed trajectories than the corresponding results from a deterministic forecast. We also investigate an alternative model in which velocity perturbations are added to the deterministic results and ensemble mean results, by a first-order stochastic process. The parameters of the stochastic model are tuned to match the dispersion of the ensemble approach. Search areas from the stochastic model give a higher hit ratio of the observations than the results based on the ensemble. However, we find that this is a consequence of a positive skew of the area distribution of the convex hulls of the ensemble trajectory end points.     

A new two-way nesting technique for ocean modeling based on the smoothed semi-prognostic method

A new two-way nesting technique is presented for a multiple nested-grid ocean modeling system. The new technique uses the smoothed semi-prognostic (SSP) method to exchange information between the different subcomponents of the nested-grid system. Four versions of the new nesting technique are described, together with conventional one-way nesting. The performance of the different nesting techniques is compared, using two independent nested-grid modeling systems, one for the Scotian Shelf of the northwest Atlantic Ocean and the other for the Meso-American Barrier Reef System of the northwestern Caribbean Sea. Nesting using the semi-prognostic method is shown to effectively prevent unrealistic drift of the inner model, while use of the SSP method avoids unnecessary damping of small scales on the inner model grid. Comparison of the annual-mean flow field with the near-surface currents determined by Fratantoni (in J Geophys Res 106:2977–2996, 2001) from observed trajectories of near-surface drifters demonstrates the overall superiority of the nesting technique based on the SSP method.     

On utility of triangular C-grid type discretization for numerical modeling of large-scale ocean flows

Ocean circulation models based on triangular C-grid discretization are frequently employed to simulate coastal ocean dynamics on unstructured meshes. It is shown that on time and space scales dominated by slow geostrophic dynamics, this discretization tends to exhibit checkerboard noise in the field of horizontal velocity divergence and vertical velocity, respectively. The noise is linked to the geometry of triangular C-grid and is amplified in regimes that are close to geostrophic balance through the particular structure of the Coriolis operator. It can be partly suppressed in some cases but remains a problem in a general case and makes the triangular C-grid a suboptimal choice for large-scale ocean modeling.     

Adaptive split-operator methods for modeling transport phenomena in porous medium systems

Split-operator methods are commonly used to approximate environmental models. These methods facilitate the tailoring of different approximation approaches to different portions of the differential operator and provide a means to split large coupled problems into pieces that are more amenable to parallel computation than the original fully-coupled problem. However, split-operator methods introduce an additional source of approximation error into the solution, which is typically either ignored or controlled heuristically. In this work, we develop two methods to estimate and control the error in split-operator methods, which lead to a dynamic adjustment of the temporal splitting step based upon the error estimators. The proposed methods are shown to yield robust solutions that provide the desired control of error. In addition, for a typical nonlinear reaction problem, the new methods are shown to reduce the solution error by more than two orders of magnitude compared to standard methods for an identical level of computational effort. The algorithms introduced and evaluated have widespread applicability in environmental modeling.     

Time-domain modeling of global ocean tides generated by the full lunisolar potential

Traditionally, ocean tides have been modeled in frequency domain with a forcing from selected tidal constituents. It is a natural approach; however, it implicitly neglects non-linearities of ocean dynamics. An alternative approach is time-domain modeling with a forcing given by the full lunisolar potential, i.e., all tidal waves are a priori included. This approach has been applied in several ocean tide models; however, some challenging tasks still remain, for example, assimilation of satellite altimetry data. In this paper, we introduce the assimilative scheme applicable in a time-domain model, which is an alternative to existing techniques used in assimilative ocean tide models. We present results from DEBOT, a global barotropic ocean tide model, which has two modes: DEBOT-h, a purely hydrodynamical mode, and DEBOT-a, an assimilative mode. The accuracy of DEBOT in both modes is assessed through a series of tests against tide gauge data which demonstrate that DEBOT is comparable to state-of-the-art global ocean tide models for major tidal constituents. Furthermore, as signals of all tidal frequencies are included in DEBOT, we also discuss modeling of minor tidal constituents and non-linear compound tides. Our modeling approach can be useful for those applications where the frequency domain approach is not suitable.     

Turbulence structure in the upper ocean: a comparative study of observations and modeling

Observations of turbulent dissipation rates measured by two independent instruments are compared with numerical model runs to investigate the injection of turbulence generated by sea surface gravity waves. The near-surface observations are made by a moored autonomous instrument, fixed at approximately 8 m below the sea surface. The instrument is equipped with shear probes, a high-resolution pressure sensor, and an inertial motion package to measure time series of dissipation rate and nondirectional surface wave energy spectrum. A free-falling profiler is used additionally to collect vertical microstructure profiles in the upper ocean. For the model simulations, we use a one-dimensional mixed layer model based on a k–ε type second moment turbulence closure, which is modified to include the effects of wave breaking and Langmuir cells. The dissipation rates obtained using the modified k–ε model are elevated near the sea surface and in the upper water column, consistent with the measurements, mainly as a result of wave breaking at the surface, and energy drawn from wave field to the mean flow by Stokes drift. The agreement between observed and simulated turbulent quantities is fairly good, especially when the Stokes production is taken into account.     

Potential volume for CO2 deep ocean sequestration: an assessment of the area located on western Pacific Ocean

Captured CO₂ could be deliberately injected into the ocean at great depth, where most of it would remain isolated from the atmosphere for centuries. CO₂ can be transported via pipeline or ship for release in the ocean or on the sea floor. In Taiwan, CO₂ release is preliminarily projected from 2010 to 2030 in an average amount of 6.957 Gt within this duration. If deep sea sequestration for CO₂ can be the possible option in Taiwan, it seems to exists possible potential area delimited between 122.0°E to 122.5°E and 21.8°N to 22.3°N for CO₂ sequestration on account of its isolated and flat topography. Apparently, the area to the southeast of Taiwan is found to reach a depth deeper than −3,000 m and can be taken as a testing area for pilot studies. This study searches the area using the contours from the depth of −4,554 to −5,500 m with 1-m interval; the area, topographic volume, maximum mean height (volume/area), and ocean volume are reported. If the emission rate is kept constantly, for 20-year storage it needs 3 m of thickness reaching the sea ridge at the depth −4,554 m using top-down style; for 100 years of storage it needs 12 m. On the other hand, if it accounts for the bottom the sea floor is taken as the reference and the accumulated CO₂ is stored from the depth at −4,900 m using bottom-up style, it requires about 37 m for the 20-year storage and 61 m for one decade.     

Wave effects in global ocean modeling: parametrizations vs. forcing from a wave model

One of the main challenges of the Copernicus Marine Service is the implementation of coupled ocean/waves systems that accurately estimate the momentum and energy fluxes provided by the atmosphere to the ocean. This study aims to investigate the impact of forcing the Nucleus for European Modelling of the Ocean (NEMO) ocean model with forecasts from the wave model of Météo-France (MFWAM) to improve classical air-sea flux parametrizations, these latter being mostly driven by the 10-m wind. Three wave-related processes, namely, wave-state-dependent stress, Stokes drift-related effects (Stokes-Coriolis force, Stokes drift advection on tracers and on mass), and wave-state-dependent surface turbulence, are examined at a global scale with a horizontal resolution of 0.25°. Three years of sensitivity simulations (2014–2016) show positive feedback on sea surface temperature (SST) and currents when the wave model is used. A significant reduction in SST bias is observed in the tropical Atlantic Ocean. This is mainly due to the more realistic momentum flux provided by the wave model. In mid-latitudes, the most interesting impact occurs during the summer stratification, when the wind is low and the wave model produces a reduction in the turbulence linked with wave breaking. Magnitudes of the large-scale currents in the equatorial region are also improved by 10% compared to observations. In general, it is shown that using the wave model reduces on average the momentum and energy fluxes to the ocean in tropical regions, but increases them in mid-latitudes. These differences are in the order of 10 to 20% compared with the classical parametrizations found in stand-alone ocean models.     

A high-resolution modeling study of the Western Iberian Margin mean and seasonal upper ocean circulation

The mean seasonal hydrography and circulation of the Western Iberian Margin (WIM) are studied by means of a high-resolution configuration of the Regional Oceanic Modeling System. A comparison of 5-year model averages for January and July with climatological datasets shows a general good agreement in the reproduction of the mean water mass properties and hydrographic distribution. We find that there is a prevailing tendency for slope poleward flow at about 80–100 km offshore at all latitudes from the surface to 1,500 m with strong vertical coupling. This northward flow, which is mainly along slope and amounts up to 8–10 cm s^?1, exhibits several mean flow recirculation regions on its way and evidences of an offshore pathway of poleward flow. Transports at different zonal sections further confirm the poleward flow tendency with two peaks of poleward transport in summer (3–10 Sv) and winter (2–7 Sv). The transport time series emphasize the seasonal character of the alongshore circulation and the interannual intrinsic variability of the circulation, since the forcing fields are climatological. As a conceptual essay with the purpose of assessing the Mediterranean Water flow influence on the WIM mean circulation, a second model configuration is setup, where the Mediterranean outflow into the study domain is removed. We find that there is an attenuation of the mesoscale field, but the slope poleward flow intensifies and remains as a mean dynamical feature closer to the upper slope.     

Improved space–time sea surface salinity mapping in Western Pacific ocean using contingogram modeling

Stochastic Environmental Research and Risk Assessment - In the context of statistical correlation theory and geostatistics, the covariance function has been widely used to characterize the...     

Adaptive multi-fidelity probabilistic collocation-based Kalman filter for subsurface flow data assimilation: numerical modeling and real-world experiment

Stochastic Environmental Research and Risk Assessment - The ensemble Kalman filter (EnKF) has received substantial attention in hydrologic data assimilation due to its ease of implementation. In...     

OpenDA-NEMO framework for ocean data assimilation

Data assimilation methods provide a means to handle the modeling errors and uncertainties in sophisticated ocean models. In this study, we have created an OpenDA-NEMO framework unlocking the data assimilation tools available in OpenDA for use with NEMO models. This includes data assimilation methods, automatic parallelization, and a recently implemented automatic localization algorithm that removes spurious correlations in the model based on uncertainties in the computed Kalman gain matrix. We have set up a twin experiment where we assimilate sea surface height (SSH) satellite measurements. From the experiments, we can conclude that the OpenDA-NEMO framework performs as expected and that the automatic localization significantly improves the performance of the data assimilation algorithm by successfully removing spurious correlations. Based on these results, it looks promising to extend the framework with new kinds of observations and work on improving the computational speed of the automatic localization technique such that it becomes feasible to include large number of observations.     

基于Remez迭代算法求解差分系数的双相介质自适应有限差分正演模拟

利用传统有限差分方法对基于Biot理论的双相介质波动方程进行数值求解时,由于慢纵波的存在,数值频散效应较为明显,影响模拟精度.相对于声学近似方程及普通弹性波方程,Biot双相介质波动方程在同等数值求解算法和精度要求条件下,其地震波场正演模拟需要更多的计算时间.本文针对Biot一阶速度-应力方程组发展了一种变阶数优化有限差分数值模拟方法,旨在同时提高其正演模拟的精度和效率.首先结合交错网格差分格式推导Biot方程的数值频散关系式.然后基于Remez迭代算法求取一阶空间偏导数的优化差分系数,并用于Biot方程的交错网格有限差分数值模拟.在此基础上把三类波的平均频散误差参数限制在给定的频散误差阈值和频率范围内,此时优化有限差分算子的长度就能自适应非均匀双相介质模型中的不同速度区间.数值频散曲线分析表明：基于Remez迭代算法的优化有限差分方法相较传统泰勒级数展开方法在大波数范围对频散误差的压制效果更明显;可变阶数的优化有限差分方法能取得与固定阶数优化有限差分方法相近的模拟精度.在均匀介质和河道模型的数值模拟实验中将本文变阶数优化有限差分算法与传统泰勒展开算法、最小二乘优化算法进行比较,进一步证明其在复杂地下介质中的有效性和适用性.

     

信噪比谱约束的自适应谱模拟反褶积方法研究

反褶积是提高地震资料分辨率的重要方法,子波估计与反褶积算子设计是反褶积方法的两个重要方面.本文在分析谱模拟法与自相关法各自特点的基础上,将两种方法的优势相结合,提出了自适应谱模拟方法,提高了谱模拟方法的适应性和准确性.在反褶积算子计算过程中对误差能量进行加权,减少了反褶积过程中处理噪音的产生,同时引入信噪比谱约束,以达到保证信噪比的前提下合理提升地震资料分辨率的目的.     

可控源电磁场三维自适应矢量有限元正演模拟

本文实现了可控源电磁（CSEM）场三维自适应矢量有限元正演算法,该算法采用非结构四面体单元进行三维网格剖分,能够真实模拟地形起伏和复杂电性异常体.采用一次场和二次场分离的方式计算电磁场响应,能够有效解决有限元模拟中的源点奇异性,提高场源附近电磁场数值精度,其中一次场利用CSEM一维正演算法解析求得,二次场采用矢量有限元方法求得.并利用基于后验误差估计的自适应网格细化算法指导网格细化,以减少人为设计网格导致的误差.通过一维和三维模型的数值模拟,验证了本文算法的有效性：一维模型有限元解与解析解吻合得很好,电磁场振幅相对误差在1%左右,相位差整体小于1°;三维模型有限元解与有限体积解吻合得也很好.模拟了一个含三维倾斜板状异常体的可控源电磁响应,表明了本文算法模拟复杂地电结构电磁场的能力和有效性.

     

Coherent variations of the obliquity components in global ice volume and ocean carbon reservoir over the past 5 Ma

The Milankovi theory stresses that the summer insolation in the high northern latitudes that is dominated by the precession cycle controls the glacial/interglacial cycles in global climate change.If the climate system responds linearly to the external insolation forcing,the precession cycle of 23 or 19 ka should dominate the variations in the climatic proxy records.I performed spectral and evolutive cross spectral analyses on the high resolution benthic 18O and 13C records from the South China Sea and the North Atlantic,the proxies of global ice volume and ocean carbon reservoir respectively.I found that the obliquity instead of the eccentricity or the precession is the most marked cycle in the global ice volume and ocean carbon reservoir variations over the past 5 Ma.The analysis further reveals that only at the obliquity band instead of the eccentricity or the precession band does the global ice volume and ocean carbon reservoir display consistently high coherency and stable phase relationship over the past 5 Ma.The consistently positive or near-zero phases of the benthic 18O relative to the benthic13C at the obliquity band suggest that the global carbon cycle is involved in the polar ice sheet growth as an important internal feedback,not a determinative driving factor.The obliquity instead of the precession or the eccentricity takes the dominant role of driving the global climate change during the Pliocene and Pleistocene.