Numerical study of three-dimensional shelf circulation on the Scotian Shelf using a shelf circulation model

A numerical shelf circulation model was developed for the Scotian Shelf, using a nested-grid setup consisting of a three-dimensional baroclinic inner model embedded inside a two-dimensional barotropic outer model. The shelf circulation model is based on the Princeton Ocean Model and driven by three-hourly atmospheric forcing provided by a numerical weather forecast model and by tidal forcing specified at the inner model's open boundaries based on pre-calculated tidal harmonic constants. The outer model simulates the depth-mean circulation forced by wind and atmospheric pressure fields over the northwest Atlantic Ocean with a horizontal resolution of 1/12°. The inner model simulates the three-dimensional circulation over the Gulf of St. Lawrence, the Scotian Shelf, and the adjacent slope with a horizontal resolution of 1/16°. The performance of the shelf circulation model is assessed by comparing model results with oceanographic observations made along the Atlantic coast of Nova Scotia and in the vicinity of Sable Island (on the Scotian Shelf) during two periods: October 2000–March 2001 and April–June 2002. Analysis of model results on Sable Island Bank indicates that tidal currents account for as much as ∼80% of the total variance of near-bottom currents, and currents driven by local winds account for ∼30% of the variance of the non-tidal near-bottom currents. Shelf waves generated remotely by winds and propagating into the region also play an important role in the near-bottom circulation on the bank.     

Hydrodynamics and sediment dynamics in the German Bight. A focus on observations and numerical modelling in the East Frisian Wadden Sea

This work deals with analysis of hydrographic observations and results of numerical simulations. The data base includes acoustic Doppler current profilers (ADCP) observations, continuous measurements on data stations and satellite data originating from the medium resolution imaging spectrometer (MERIS) onboard the European Space Agency (ESA) satellite ENVISAT with a spatial resolution of 300 m. Numerical simulations use nested models with horizontal resolutions ranging from 1 km in the German Bight to 200 m in the East Frisian Wadden Sea coupled with a suspended matter transport model. Modern satellite observations have now a comparable horizontal resolution with high-resolution numerical model of the entire area of the East Frisian Wadden Sea allowing to describe and validate new and so far unknown patterns of sediment distribution. The two data sets are consistent and reveal an oscillatory behaviour of sediment pools to the north of the back-barrier basins and clear propagation patterns of tidally driven suspended particulate matter outflow into the North Sea. The good agreement between observations and simulations is convincing evidence that the model simulates the basic dynamics and sediment transport processes, which motivates its further use in hindcasting, as well as in the initial steps towards forecasting circulation and sediment dynamics in the coastal zone.     

Dynamics of the circulation in the Sea of Marmara: numerical modeling experiments and observations from the Turkish straits system experiment

During September 2008 and February 2009, the NR/V Alliance extensively sampled the waters of the Sea of Marmara within the framework of the Turkish Straits System (TSS) experiment coordinated by the NATO Undersea Research Centre. The observational effort provided an opportunity to set up realistic numerical experiments for modeling the observed variability of the Marmara Sea upper layer circulation at mesoscale resolution over the entire basin during the trial period, complementing relevant features and forcing factors revealed by numerical model results with information acquired from in situ and remote sensing datasets. Numerical model solutions from realistic runs using the Regional Ocean Modeling System (ROMS) produce a general circulation in the Sea of Marmara that is consistent with previous knowledge of the circulation drawn from past hydrographic measurements, with a westward meandering current associated with a recurrent large anticyclone. Additional idealized numerical experiments illuminate the role various dynamics play in determining the Sea of Marmara circulation and pycnocline structure. Both the wind curl and the strait flows are found to strongly influence the strength and location of the main mesoscale features. Large displacements of the pycnocline depth were observed during the sea trials. These displacements can be interpreted as storm-driven upwelling/downwelling dynamics associated with northeasterly winds; however, lateral advection associated with flow from the Straits also played a role in some displacements.     

Numerical simulation of the influence of stationary mesoscale orographic waves on the meridional circulation and ozone fluxes in the middle atmosphere

The authors’ parameterization of the dynamic and thermal action of stationary orographic waves generated by the Earth’s surface relief is included into the model of general circulation of the middle and upper atmosphere. Numerical simulation of the general circulation in the troposphere and stratosphere was performed and the influence of stationary orographic waves propagating upward from the Earth’s surface on the meridional and vertical velocity was studied. It is shown that the allowance for the dynamic and thermal action of these waves in the numerical model leads to changes by up to 20–30% in the meridional circulation and ozone fluxes associated with it at heights of the ozone layer maximum.     

Direct BEM for high-resolution global gravity field modelling

The paper presents a high-resolution global gravity field modelling by the boundary element method (BEM). A direct BEM formulation for the Laplace equation is applied to get a numerical solution of the linearized fixed gravimetric boundary-value problem. The numerical scheme uses the collocation method with linear basis functions. It involves a discretization of the complicated Earth’s surface, which is considered as a fixed boundary. Here 3D positions of collocation points are simulated from the DNSC08 mean sea surface at oceans and from the SRTM30PLUS_V5.0 global topography model added to EGM96 on lands. High-performance computations together with an elimination of the far zones’ interactions allow a very refined integration over the all Earth’s surface with a resolution up to 0.1 deg. Inaccuracy of the approximate coarse solutions used for the elimination of the far zones’ interactions leads to a long-wavelength error surface included in the obtained numerical solution. This paper introduces an iterative procedure how to reduce such long-wavelength error surface. Surface gravity disturbances as oblique derivative boundary conditions are generated from the EGM2008 geopotential model. Numerical experiments demonstrate how the iterative procedure tends to the final numerical solutions that are converging to EGM2008. Finally the input surface gravity disturbances at oceans are replaced by real data obtained from the DNSC08 altimetryderived gravity data. The ITG-GRACE03S satellite geopotential model up to degree 180 is used to eliminate far zones’ interactions. The final high-resolution global gravity field model with the resolution 0.1 deg is compared with EGM2008.     

Electric field effects on ionospheric and thermospheric parameters above the EISCAT station for summer conditions

Numerical calculations of the thermospheric and ionospheric parameters above EISCAT are presented for quiet geomagnetic conditions in summer. The Global Self-consistent Model of the Thermosphere, Ionosphere and Protonosphere (GSM TIP) was used. The numerical results were obtained both with a self-consistent calculation of the electric fields of magnetospheric and dynamo-action origin and with the magnetospheric electric fields only. It was found that the dynamo-electric field has some effect on the ionospheric convection pattern during quiet geomagnetic conditions. It has a marked effect mainly on the zonal neutral wind component above EISCAT (±20m/s at 140 km altitude). We have studied the effects of various field-aligned current (FAC) distributions on thermosphere/ionosphere parameters and we show that a qualitative agreement can be obtained with region-I and -II FAC zones at 75° and 65° geomagnetic latitude, respectively. The maximum FAC intensities have been assumed at 03–21 MLT for both regions with peak values of 2.5 × 10^–7 Am^–2 (region I) and 1.25 × 10^–7 A m^–2 (region II). These results are in agreement with statistical potential distribution and FAC models constructed by use of EISCAT data. The lack of decreased electron density in the night-time sector as observed by the EISCAT radar was found to be due to the spatial distribution of ionospheric convection resulting from electric fields of magnetospheric origin.     

The effect of horizontal resolution on gravity waves simulated by the GFDL “SKYHI” general circulation model

To examine the effects of horizontal resolution on internal gravity waves simulated by the 40-level GFDL SKYHI general circulation model, a comparison is made between the 3° and 1° resolution models during late December. The stratospheric and mesospheric zonal flows in the winter and summer extratropical regions of the 1° model are much weaker and more realistic than the corresponding zonal flows of the 3° model. The weaker flows are consistent with the stronger Eliassen-Palm flux divergence (EPFD).The increase in the magnitude of the EPFD in the winter and summer extratropical mesospheres is due mostly to the increase in the gravity wave vertical momentum flux convergence (VMFC). In the summer extratropical mesosphere, the increase in the resolvable horizontal wavenumbers accounts for most of the increase in the gravity wave VMFC. In the winter extratropical mesosphere, the increase of VMFC associated with large-scale eastward moving components also accounts for part of the increase in the gravity wave VMFC.The gravity waves in the summer and winter mesosphere of the 1° model are associated with a broader frequency-spectral distribution, resulting in a more sporadic time-distribution of their VMFC. This broadening is due not only to the increase in resolvable horizontal wavenumbers but also occurs in the large-scale components owing to wave-wave interactions. It was found that the phase velocity and frequency of resolvable small-scale gravity waves are severely underestimated by finite difference approximations.     

用全球大气环流模式嵌套区域气候模式模拟东亚现代气候

将区域气候模式RegCM2与中国科学院大气物理研究所的9层全球格点大气环流模式IAP AGCM单向嵌套,对东亚现代气候进行数值模拟研究,同时检验和分析该嵌套模式的性能.已完成的10年积分结果表明,单向嵌套RegCM2由于具有较高分辨率和较完善的物理过程,因此对地面气温和降水的空间分布形势和季节变化趋势都有较好的模拟能力,且较与之嵌套的IAP AGCM的模拟效果有较大改善,如在中国区域,它模拟的年均地面气温与实况的空间相关系数由全球环流模式的092提高到094,模拟的年均降水由05提高到07. 这与嵌套RegCM2能模拟出IAP AGCM所不能分辨的中尺度信号有很大关系.     

利用径向基函数RBF解算GRACE全球时变重力场

本文利用GRACE(Gravity Recovery And Climate Experiment)level 1b数据和径向基函数RBF(radial basis function)方法解算了全球时变地球重力场.RBF基函数相比传统球谐(spherical harmonic)基函数,其高度的空域局部特性使得正则化过程易于添加先验协方差信息,从而可能揭示更加准确的重力场信号.本文研究表明,RBF基函数算法在精化现有的GRACE全球时变重力场模型,如提升部分区域信号幅度等方面具有一定优势.本文通过将RBF的尺度因子作为待解参数,基于GRACE卫星的Level 1b数据和变分方程法,成功获取了2009-2010年90阶无约束全球时变重力场RBF模型Hust-IGG03,以及正则化全球时变重力场RBF模型Hust-IGG04.通过与GRACE官方数据处理中心GFZ发布的最新90阶球谐基时变模型RL05a进行对比,结果表明:(1)无约束RBF模型Hust-IGG03和GFZ RL05a在空域和频域表现基本一致;(2)正则化RBF模型Hust-IGG04无需进行后处理滤波已经显示较高信噪比,噪音水平接近于球谐基模型GFZ RL05a经400 km高斯滤波后的效果;(3)HustIGG04相比400 km高斯滤波GFZ RL05a在周年振幅图和趋势图上显示出更多的细节信息,并且呈现出更强的信号幅度,如在格陵兰冰川融化趋势估计上Hust-IGG04比GFZ RL05a提高了24.2%.以上结果均显示RBF方法有助于进一步挖掘GRACE观测值所包含的时变重力场信息.     

Seasonal variation of upper layer circulation in the northern part of the East/Japan Sea

Seasonal variation of upper layer circulation in the northern part of the East/Japan Sea and its mechanism were investigated using empirical orthogonal function (EOF) analysis with satellite sea surface heights over the northern East/Japan Sea and a three-dimensional circulation model. The spatial structure and temporal variation of first EOF mode, which explains about 64% of the total variance, indicate that a large cyclonic circulation in the northern East/Japan Sea shows a semi-annual variation with maximum strength in summer and winter. According to numerical model result, the Liman Cold Current, accepted as a major current in the northern East/Japan Sea, is well mixed vertically by the winter monsoon and the current in the upper layer has a relatively deep structure, with a maximum westward speed of about 20 cm/s in winter. On the other hand, in summer the current has a stronger baroclinic structure of velocity than in winter. Numerical experiments showed that in summer the temporal variation of upper layer circulation is controlled by thermal forcing, such as sea surface heat flux and inflow of heat transport into the East/Japan Sea through the Korea/Tsushima Strait. Moreover, the cyclonic circulation in the upper layer of the northern East/Japan Sea is also generated and strengthened by the positive wind stress curl occupying most of the East/Japan Sea during the winter. The seasonal variation of each forcing that drives the circulation is responsible for the strength or weakness of the upper layer circulation in the northern East/Japan Sea. The contribution of each forcing to the seasonal variation of the upper layer circulation is examined through sensitivity experiments. According to these numerical experiments, the upper layer circulation in the northern East/Japan Sea is strengthened twice a year, in winter and summer, and this semi-annual variation is determined by a combination of wind (winter) and thermal (summer) forcing.     

Numerical analysis of the Black Sea currents and mesoscale eddies in 2006 and 2011

Two prognostic experiments taking into account real atmospheric forcing for 2006 and 2011 were carried out based on the eddy-resolving numerical model with a horizontal resolution of 1.6 km for the Black Sea. The main dynamic features such as the Rim Current, the Sevastopol, and Batumi anticyclones are reproduced in both experiments. The model results are confirmed via observation data. We accomplished the analysis of simulated circulation and energetics. The results demonstrate that both the vertical viscosity and vertical diffusion along with the energy inflow from the wind have been the main contributors to the annual and seasonal budgets of kinetic and potential energies of the Black Sea circulation. It is shown that two regimes of the Black Sea general circulation are implemented depending on a magnitude of wind contribution to the kinetic energy in winter. Intensive mesoscale eddy formation was observed along the Anatolian, Caucasian, and Crimean coasts. The analysis of the Black Sea circulation and eddy energetics allowed us to conclude that the generation and development of the mesoscale coastal eddies is associated with the barotropic instability in case of intensive coastal currents and is associated with both the barotropic and baroclinic instability in case of weak coastal currents.     

First results from ground-based daytime optical investigation of the development of the equatorial ionization anomaly

A meridional scanning OI 630.0-nm dayglow photometer was operated from Ahmedabad (17.2^°N dip lat.) scanning a region towards the south in the upper atmosphere extending over \sim5^° in latitude from 10.2^°N to 15.2^°N dip latitude. From the spatial and temporal variabilities of the dayglow intensity in the scanning region we show for the first time, evidence for the passage of the crest of the equatorial ionization anomaly (EIA) in the daytime by means of a ground-based optical technique. The relationship between the daytime eastward electric field over the dip equator in the same longitude zone as inferred from the equatorial electrojet strength and the evolutionary pattern of EIA is clearly demonstrated. The latter as inferred from the dayglow measurements is shown to be consistent with our present understanding of the electrodynamical processes in the equatorial region. The present results reveal the potential of this ground-based optical technique for the investigation of ionospheric/thermospheric phenomena with unprecedented spatial and temporal resolution.     

Impact of partial steps and momentum advection schemes in a global ocean circulation model at eddy-permitting resolution

Series of sensitivity tests were performed with a z-coordinate, global eddy-permitting (1/4°) ocean/sea-ice model (the ORCA-R025 model configuration developed for the DRAKKAR project) to carefully evaluate the impact of recent state-of-the-art numerical schemes on model solutions. The combination of an energy–enstrophy conserving (EEN) scheme for momentum advection with a partial step (PS) representation of the bottom topography yields significant improvements in the mean circulation. Well known biases in the representation of western boundary currents, such as in the Atlantic the detachment of the Gulf Stream, the path of the North Atlantic Current, the location of the Confluence, and the strength of the Zapiola Eddy in the south Atlantic, are partly corrected. Similar improvements are found in the Pacific, Indian, and Southern Oceans, and characteristics of the mean flow are generally much closer to observations. Comparisons with other state-of-the-art models show that the ORCA-R025 configuration generally performs better at similar resolution. In addition, the model solution is often comparable to solutions obtained at 1/6 or 1/10° resolution in some aspects concerning mean flow patterns and distribution of eddy kinetic energy. Although the reasons for these improvements are not analyzed in detail in this paper, evidence is shown that the combination of EEN with PS reduces numerical noise near the bottom, which is likely to affect current–topography interactions in a systematic way. We conclude that significant corrections of the mean biases presently seen in general circulation model solutions at eddy-permitting resolution can still be expected from the development of numerical methods, which represent an alternative to increasing resolution.     

Numerical calculations of tectonic stress field of Chinese mainland and its neighboring regions and their applications to explanation of seismic activity

IntroductionChinesemainlandliesinthesoutheastofEurasianPlateandisclampedbyindian,PacificandPhilippineSeaPlates.Besides,itislocatedatbetweenthetwobigseismicactivebeltsintheworldatthesametime.Differingfromtheshowoftwobigseismicbeltswhichhaveseismiczonaldistributionalongplateboundary,Chinesemainlandanditsperipheralareahaveapeculiarcharacter,expressingseismicallylargestrength,highfrequency,widedistribution,shallowfocus,seriousdamage,multiplestyles,extremelynonuniformspace-timedistributionandsoo…     

A two-way nested-grid ocean-circulation model for the Meso-American Barrier Reef System

A two-way nested-grid ocean-circulation model is developed for the Meso-American Barrier Reef System (MBRS), using a newly developed two-way interactive nesting technique. The unique feature of this new nesting technique is its use of the semi-prognostic method (Sheng et al. 2001) to exchange information between different grids through the model momentum equations. The nested-grid model for the MBRS has a fine-resolution inner model embedded in a coarse-resolution outer model. The outer model is the western Caribbean Sea model developed by Sheng and Tang (2003), with a horizontal resolution of roughly 19 km. The inner model domain covers the northwest Caribbean Sea (NWCS) between 79°W and 89°W and between 15.5°N and 22°N, with a horizontal resolution of roughly 6 km. The nested-grid ocean model is initialized with the January mean temperature and salinity and forced by the monthly mean COADS (comprehensive ocean-atmosphere data set) wind stress and surface heat flux. The model sea-surface salinity is restored to the monthly mean climatology. The nested-grid model is integrated for 2 years and the second-year model results are presented in this paper. The model-calculated annual-mean near-surface currents over the NWCS agree reasonably well with the time-mean near-surface currents inferred by Fratantoni (2001) from trajectories of the satellite-tracked 15-m drogued drifters in the 1990s. The two-way nested model is also used to quantify the role of local wind stress, local density gradients and boundary forcings of the outer model in driving the annual-mean circulation in the region.Responsible Editor: Phil Dyke     

Assimilation of historical head data to estimate spatial distributions of stream bed and aquifer hydraulic conductivity fields

Management of water resources in alluvial aquifers relies mainly on understanding interactions between hydraulically connected streams and aquifers. Numerical models that simulate this interaction often are used as decision support tools for water resource management. However, the accuracy of numerical predictions relies heavily on unknown system parameters (e.g., streambed conductivity and aquifer hydraulic conductivity), which are spatially heterogeneous and difficult to measure directly. This paper employs an ensemble smoother to invert groundwater level measurements to jointly estimate spatially varying streambed and alluvial aquifer hydraulic conductivity along a 35.6‐km segment of the South Platte River in Northeastern Colorado. The accuracy of the inversion procedure is evaluated using a synthetic experiment and historical groundwater level measurements, with the latter constituting the novelty of this study in the inversion and validation of high‐resolution fields of streambed and aquifer conductivities. Results show that the estimated streambed conductivity field and aquifer conductivity field produce an acceptable agreement between observed and simulated groundwater levels and stream flow rates. The estimated parameter fields are also used to simulate the spatially varying flow exchange between the alluvial aquifer and the stream, which exhibits high spatial variability along the river reach with a maximum average monthly aquifer gain of about 2.3 m³/day and a maximum average monthly aquifer loss of 2.8 m³/day, per unit area of streambed (m²). These results demonstrate that data assimilation inversion provides a reliable and computationally affordable tool to estimate the spatial variability of streambed and aquifer conductivities at high resolution in real‐world systems.     

On the circulation in the East Frisian Wadden Sea: numerical modeling and data analysis

 In this paper we use a combination of numerical modeling and data analysis to gain a better understanding of the major characteristics of the circulation in the East Frisian Wadden Sea. In particular, we concentrate on the asymmetry of the tidal wave and its modulation in the coastal area, which results in a complex pattern of responses to the sea-level forcing from the North Sea. The numerical simulations are based on the 3-D primitive equation General Estuarine Transport Model (GETM) with a horizontal resolution of 200 m and terrain-following vertical coordinates. The model is forced at its open boundaries with sea-level data from an operational model for the German Bight (German Hydrographic Office). The validation data for our model simulations include time series of tidal gauge data and surface currents measured at a pile in the back-barrier basin of the Island Langeoog, as well as several ADCP transects in the Accumer Ee tidal inlet. Circulation and turbulence characteristics are investigated for typical situations driven by spring and neap tides, and the analysis is focused on dominating temporal and spatial patterns. By investigating the response of five back-barrier basins with rather different morphologies to external forcing, an attempt is made to elucidate the dominating physical balances controlling the circulation in the individual sub-basins. It is demonstrated that the friction at the seabed tends to slow down the tidal signal in the shallow water. This leads to the establishment of flood dominance in the shallow sea north of the barrier islands. South of the islands, where the water volume of the channels at low tide is smaller than the tidal prism, the asymmetry of the tidal signal is shifted towards ebb dominance, a feature which is particularly pronounced at spring tide. At the northern open boundary, the tidal wave propagating from west to east generates a sea-level difference of ∼1 m along the boundary, and thereby triggers vigorous alongshore currents. The frictional control in the model is located in the inlets, as well as along the northern boundary. The correlation between velocity and turbulent kinetic energy tends to the establishment of a net southward transport, giving theoretical support to the observed accumulation of sediments on the intertidal flats. Weak turbulence along the northern shores of the barrier islands and the small magnitude of the residual currents there promote accumulation of suspended matter in these areas, although wave action will generally counteract this effect. Received: 29 May 2002 / Accepted: 26 September 2002 Responsible Editor: Jean-Marie Beckers Acknowledgements We are indebted to S. Dick for providing the data from the operational model of BSH and to B. Flemming for the useful discussions. The topography data and Fig. 1 have been prepared in cooperation with F. Meyer. Figure 2 has been prepared by G. Brink-Spalink. We also thank for the comments from an anonymous reviewer which helped to improve our paper.     

基于GRACE卫星RL05数据的南极冰盖质量变化分析

CSR(Centre for Space Research)最近发布了RL05数据,其空间分辨率、精度和周期变化特性等都优于RL04数据.本文采用300 km的扇形滤波及P5M11去相关滤波削弱南北条带等重力场模型误差,并采用Paulson2007模型进行冰川均衡模型改正,利用CSR RL05与RL04数据计算分析了南极2002年到2012年的质量变化序列及其变化趋势的空间分布特性,并选取8个特征点进一步分析了其质量变化序列.同时,对CSR、JPL(Jet Propulsion Laboratory)、GFZ (GeoForschungsZentrum)三个机构发布的RL05数据采用相同的滤波方法进行计算,得到整个南极的质量变化分别为-195.7±20.5 Gt/a、-203.8±23.1 Gt/a、-133.2±29.9 Gt/a,对全球海平面变化的影响分别为0.54±0.06 mm/a、0.56±0.06 mm/a、0.37±0.09 mm/a.     

Numerical Simulation of the Caspian Sea Circulation Using the Marine and Atmospheric Research System

The Marine and Atmospheric Research System (MARS) for the Caspian Sea meteorological characteristics is presented, which is implemented in Zubov State Oceanographic Institute. It includes computation of the atmospheric forcing with the Weather Research and Forecasting model, as well as computation of currents, sea level, temperature, salinity and sea ice with the Institute of Numerical Mathematics Ocean Model and the computation of wind wave parameters using the Russian wind-wave model. The results are presented on verification of the hydrometeocharacteristics simulated with the MARS for the Caspian Sea. As well, the retrospective simulation of the thermohydrodynamic characteristics in this basin is performed with MARS for the ice-free period 2003–2013. The important features of the Caspian Sea circulation are shown.