引潮力对海洋环流模式的影响

The eight main tidal constituents have been implemented in the global ocean general circulation model with approximate 1° horizontal resolution.Compared with the observation data,the patterns of the tidal amplitudes and phases had been simulated fairly well.The responses of mean circulation,temperature and salinity are further investigated in the global sense.When implementing the tidal forcing,wind-driven circulations are reduced,especially those in coastal regions.It is also found that the upper cell transport of the Atlantic meridional overturning circulation(AMOC) reduces significantly,while its deep cell transport is slightly enhanced from 9×106m3/s to 10×106 m3/s.The changes of circulations are all related to the increase of a bottom friction and a vertical viscosity due to the tidal forcing.The temperature and salinity of the model are also significantly affected by the tidal forcing through the enhanced bottom friction,mixing and the changes in mean circulation.The largest changes occur in the coastal regions,where the water is cooled and freshened.In the open ocean,the changes are divided into three layers:cooled and freshened on the surface and below 3 000 m,and warmed and salted in the middle in the open ocean.In the upper two layers,the changes are mainly caused by the enhanced mixing,as warm and salty water sinks and cold and fresh water rises;whereas in the deep layer,the enhancement of the deep overturning circulation accounts for the cold and fresh changes in the deep ocean.     

Frequency management support for remote sea-state sensing using the JINDALEE skywave radar

Successful operation of a skywave (over-the-horizon, OTH) radar in a remote sea-state sensing mode is critically dependent upon the application of comprehensive frequency management techniques. In addition to the problem of selecting a frequency yielding an adequate signal-to-noise ratio in the geographical area under investigation, attention must be paid to minimization of ionospheric multimode and other phenomena capable of distorting or convoluting the sea backscatter spectrum. This paper describes the manner in which these problems have been addressed in the JINDALEE skywave radar, and relates the practical difficulties inherent in the task of frequency management in support of an OTH radar involved in sea-state sensing. Measurement techniques include backscatter and oblique-incidence sounding, HF spectral surveillance, and a low-powered frequency-agile "miniradar" capable of operating in either a conventional backscatter or alternate oblique-incidence mode. In addition to providing the main radar with real-time frequency management advice, a principal emphasis has been the acquisition of a synoptic data base suitable for off-line statistical analysis.     

A spectral barotropic model of the wind-driven world ocean

A global spectral barotropic ocean model is introduced to describe the depth-averaged flow. The equations are based on vorticity and divergence (instead of horizontal momentum); continents exert a nearly infinite drag on the fluid. The coding follows that of spectral atmospheric general circulation models using triangular truncation and implicit time integration to provide a first step for seamless coupling to spectral atmospheric global circulation models and an efficient method for filtering of ocean wave dynamics. Five experiments demonstrate the model performance: (i) Bounded by an idealized basin geometry and driven by a zonally uniform wind stress, the ocean circulation shows close similarity with Munk’s analytical solution. (ii) With a real land–sea mask the model is capable of reproducing the spin-up, location and magnitudes of depth-averaged barotropic ocean currents. (iii) The ocean wave-dynamics of equatorial waves, excited by a height perturbation at the equator, shows wave dispersion and reflection at eastern and western coastal boundaries. (iv) The model reproduces propagation times of observed surface gravity waves in the Pacific with real bathymetry. (v) Advection of tracers can be simulated reasonably by the spectral method or a semi-Langrangian transport scheme. This spectral barotropic model may serve as a first step towards an intermediate complexity spectral atmosphere–ocean model for studying atmosphere–ocean interactions in idealized setups and long term climate variability beyond millennia.     

A spectral barotropic model of the wind-driven world ocean

A global spectral barotropic ocean model is introduced to describe the depth-averaged flow. The equations are based on vorticity and divergence (instead of horizontal momentum); continents exert a nearly infinite drag on the fluid. The coding follows that of spectral atmospheric general circulation models using triangular truncation and implicit time integration to provide a first step for seamless coupling to spectral atmospheric global circulation models and an efficient method for filtering of ocean wave dynamics. Five experiments demonstrate the model performance: (i) Bounded by an idealized basin geometry and driven by a zonally uniform wind stress, the ocean circulation shows close similarity with Munk’s analytical solution. (ii) With a real land–sea mask the model is capable of reproducing the spin-up, location and magnitudes of depth-averaged barotropic ocean currents. (iii) The ocean wave-dynamics of equatorial waves, excited by a height perturbation at the equator, shows wave dispersion and reflection at eastern and western coastal boundaries. (iv) The model reproduces propagation times of observed surface gravity waves in the Pacific with real bathymetry. (v) Advection of tracers can be simulated reasonably by the spectral method or a semi-Langrangian transport scheme. This spectral barotropic model may serve as a first step towards an intermediate complexity spectral atmosphere–ocean model for studying atmosphere–ocean interactions in idealized setups and long term climate variability beyond millennia.     

A spectral model for unidirectional nonlinear wave propagation over arbitrary depths

A weakly-nonlinear and dispersive wave equation recently developed by the authors is used for formulating a spectral-type unidirectional wave propagation model describing spectral transformations of narrow-band waves travelling over arbitrary depths. The essential characteristics of the model equation are recapitulated first and then the spectral domain representation in terms of spatially varying harmonic amplitudes is presented. The resulting evolution equations are used to simulate the experiments concerning harmonic generation in shallow water and nonlinear random wave transformations over a submerged bar. Furthermore, the spectral model predictions are compared with the field measurements in nearshore with satisfactory results.     

A spectral model of underwater irradiance in the Black Sea

We develop a semiempirical spectral model of penetrating irradiance taking into account the biooptical characteristics of the Black Sea. The evaluation of the contributions of the principal optically active components to the total absorption of light in the sea shows that, in the short-wave range (400–500 nm), light is mainly absorbed by the dissolved organic matter (41–77%). The contribution of phytoplankton to the total absorption attains its maximum values (26–37%) in the abyssal part of the sea in the period of spring blooming of diatoms. In the coastal waters, the absorption of light by suspended nonalgae particles in summer is almost twice as intense (20–30%) as in the open sea (8–13%). The analysis of the sensitivity of our model shows that the absorption of light by dissolved organic matter is more significant for the estimation of the photosynthetically active radiation in the Black Sea than the concentration of pigments and backscattering of light by suspended particles. The comparison of the results of model computations with the data of measurements of the underwater irradiance reveals high accuracy of the proposed model.     

A model for Doppler peak spectral shift for low grazing angle sea scatter

A model is formulated for Doppler spectral characteristics of radar sea scatter for low grazing angles, and is compared with previous radar measurements reported in the literature. The Doppler model is based upon the two-scale model for radar scatter, with scatterer motions hypothesized as due to the orbital wave velocity of the large-scale waves, Stokes and wind drift currents, and the phase velocity of the small-scale Bragg scatterers. Expressions for Doppler shifts due to these motions are derived, and are given as a function of wave height, wave period, and wind speed. Although this model appears to account for the peak Doppler shift of the sea-scatter Doppler spectrum for vertical polarization, it is insufficient to describe horizontally and cross-polarized data, which have larger mean Doppler shifts. However, these two cases are found to scale very closely with the nearly simultaneous vertically polarized data for the variety of environmental conditions reported. Implications of the extension of these results to higher-angle remote-sensing applications are discussed.     

A model for numerical simulation of nonstationary sonar reverberation using linear spectral prediction

An innovative approach to the numerical generation of nonstationery reverberation time series is presented and demonstrated. The computer simulated reverberation time series are of high quality, in that they are accurate representations of those which would result from an actual sonar system (transmit/receive and horizontal/ vertical beampatterns; pulse type, shape, length, and power; frequency and sampling rate), platform (speed and depth), and environment (wind speed and direction, backscattering strengths, and propagation loss). Volume, surface, and/or bottom reverberation as seen by a multiple beam sonar on a moving platform is generated. The approach utilizes recent developments in linear spectral prediction research in which the spectra of stochastic processes are modeled as rational functions and algorithms are used to efficiently compute optimal estimates of coefficients which specify the spectra. A two-fold sequence is formulated; first, the expected reverberation spectra for all beams are predicted and, second, the stochastic time series are generated from the expected spectra. The expected spectra are predicted using a numerical implementation, referred to as the REVSPEC (reverberation spectrum) model, of a general formulation of Faure, Ol'shevskii, and Middleton. Given the spectra, the Levinson-Durbin method is used to solve the Yule-Walker equations of the autoregressive formulation of linear spectral prediction. The numerical implementation of the approach, referred to as the REVSIM (reverberation simulation) model, produces nonstationary coherent multiple-beam reverberation time series. The formulation of the REVSIM model is presented and typical results given. A comparison is made between the simulation outputs of the REVSIM model and those of the REVGEN (reverberation generator) model, a standard well-accepted time series simulation model, to demonstrate the validity of the new approach.     

A global heat and freshwater forcing dataset for ocean models

A global dataset based on the ECMWF Re-Analyses (ERA) is presented that can be used as surface boundary conditions for ocean models with sea-ice components. The definition of these conditions is based on bulk formulae. To study the mean ocean circulation, a mean annual cycle on a daily basis was constructed from ERA for all relevant parameters including wind stress. Continental runoff is considered by using information about the catchment areas of the rivers and about the main drainage basins. The bulk formulae were extended by using sea ice concentration.To estimate meridional heat transports (MHT) and to avoid any drift in ocean model simulations, the heat and fresh water budgets have been closed by applying an inverse procedure to fine-tune the fluxes towards observed transports. To improve the MHTs on the Southern Hemisphere the winds and the short wave radiation at southern higher latitudes should be corrected. Furthermore, tests were performed concerning short wave radiation which was increased in the tropics and decreased in the subsidence zones.The heat and fresh water fluxes are assessed by using a scheme of Macdonald and Wunsch based on hydrographic sections. The net heat fluxes of ERA and of the forcing dataset are consistent with the heat flux divergences and convergences estimated by this scheme except for parts of the South Atlantic and the Indian Ocean sector of the Southern Ocean where none of these datasets is consistent with these estimates. In the subtropical South Indian Ocean the forcing dataset is consistent with these estimates while ERA are not. The flux components of ERA and the forcing dataset were compared to several observational datasets (SRB, SOC, HOAPS, GPCP, and CMAP). For each component, at least one of these datasets (especially HOAPS) supports the effects of the inverse procedure and the bulk formulae almost globally with some regional exceptions: short wave radiation in the tropical oceans and the subtropical North Atlantic, latent heat flux at higher latitudes, and precipitation in the northern North Atlantic.Comparisons to the NCAR/NCEP Re-Analyses (NRA) (versions 1 and 2) and the ECHAM model in place of ERA lead to similar results. In the North Atlantic the net heat fluxes of the model based datasets approach the hydrographic estimate with increasing resolution. Applied to any ocean/sea-ice model and compared to ERA, the forcing dataset would induce only a relative small net sea-surface buoyancy loss.A comparison of the forcing dataset to measurements made using one buoy deployed in the western Pacific warm pool and five buoys deployed in the subduction region of the Northeast Atlantic shows that at the site of the first buoy the net heat fluxes of the forcing dataset are in poorer agreement than those of ERA. At the sites of two subduction buoys both datasets show the same level of agreement within the error bars specified. At the sites of the three remaining subduction buoys the forcing dataset shows a marginal improvement on ERA.     

Atmospheric forcing of the eastern tropical Pacific: A review

The increase in marine, land surface, atmospheric and satellite data during recent decades has led to an improved understanding of the air–sea interaction processes in the eastern tropical Pacific. This is also thanks to extensive diagnoses from conceptual and coupled ocean–atmosphere numerical models. In this paper, mean fields of atmospheric variables, such as incoming solar radiation, sea level pressure, winds, wind stress curl, precipitation, evaporation, and surface energy fluxes, are derived from global atmospheric data sets in order to examine the dominant features of the low level atmospheric circulations of the region. The seasonal march of the atmospheric circulations is presented to depict the role of radiative forcing on atmospheric perturbations, especially those dominating the atmosphere at low levels.In the tropics, the trade winds constitute an important north–south energy and moisture exchange mechanism (as part of the low level branch of the Hadley circulation), that determines to a large extent the precipitation distribution in the region, i.e., that associated with the Inter-Tropical Convergence Zone (ITCZ). Monsoonal circulations also play an important role in determining the warm season precipitation distribution over the eastern tropical Pacific through a large variety of air–sea–land interaction mechanisms. Westward traveling waves, tropical cyclones, low latitude cold air intrusions, and other synoptic and mesoscale perturbations associated with the ITCZ are also important elements that modulate the annual rainfall cycle. The low-level jets of the Gulf of California, the Intra-Americas Sea (Gulf of Mexico and Caribbean Sea) and Chocó, Colombia are prominent features of the eastern tropical Pacific low-level circulations related to sub-regional and regional scale precipitation patterns. Observations show that the Intra-Americas Low-Level Jet intensity varies with El Niño/Southern Oscillation (ENSO) phases, however its origin and role in the westward propagation and development of disturbances that may hit the eastern tropical Pacific, such as easterly waves and tropical cyclones, are still unclear. Changes in the intensity of the trade winds in the Caribbean Sea and the Gulf of Mexico (associated with eastern tropical Pacific wind jets) exert an important control on precipitation by means of wind–topography interactions. Gaps in the mountains of southern Mexico and Central America allow strong wind jets to pass over the continent imprinting a unique signal in sea surface temperatures and ocean dynamics of the eastern tropical Pacific.The warm pools of the Americas constitute an important source of moisture for the North American Monsoon System. The northeastern tropical Pacific is a region of intense cyclogenetic activity, just west of the coast of Mesoamerica. Over the oceanic regions, large-scale properties of key variables such as precipitation, moisture, surface energy fluxes and wind stress curl are still uncertain, which inhibits a more comprehensive view of the region and stresses the importance of regional field experiments. Progress has been substantial in the understanding of the ocean and atmospheric dynamics of the eastern tropical Pacific, however, recent observational evidence such as that of a shallow meridional circulation cell in that region, in contrast to the classic concept of the Hadley-type deep meridional circulation, suggests that more in situ observations to validate theories are still necessary.This paper is part of a comprehensive review of the oceanography of the eastern tropical Pacific Ocean.     

Influence of direct sulfate-aerosol radiative forcing on the results of numerical experiments with a climate model of intermediate complexity

The climate model of intermediate complexity developed at the Institute of Atmospheric Physics of the Russian Academy of Sciences (IAP RAS CM) is extended by a block for the direct anthropogenic sulfate-aerosol (SA) radiative forcing. Numerical experiments have been performed with prescribed scenarios of the greenhouse and anthropogenic sulfate radiative forcings from observational estimates for the 19th and 20th centuries and from SRES scenarios A1B, A2, and B1 for the 21st century. The globally averaged direct anthropogenic SA radiative forcing F _ASA by the end of the 20th century relative to the preindustrial state is ?0.34 W/m², lying within the uncertainty range of the corresponding present-day estimates. The absolute value of F _ASA is the largest in Europe, North America, and southeastern Asia. A general increase in direct radiative forcing in the numerical experiments that have been performed continues until the mid-21st century. With both the greenhouse and the sulfate loadings included, the global climate warming in the model is 1.5–2.8 K by the end of the 21st century relative to the late 20th century, depending on the scenario, and 2.1–3.4 K relative to the preindustrial period. The sulfate aerosol reduces global warming by 0.1–0.4 K in different periods depending on the scenario. The largest slowdown (>1.5 K) occurs over land at middle and high latitudes in the Northern Hemisphere in the mid-21st century for scenario A2. The IAP RAS CM response to the greenhouse and the aerosol forcing is not additive.     

An unstructured spectral/hp element model for enhanced Boussinesq-type equations

A spectral/hp element method for solving enhanced Boussinesq-type equations in two horizontal dimensions is introduced. The numerical model is based on the discontinuous Galerkin method on unstructured meshes with expansions of arbitrary order. Numerical computations are used to illustrate that the computational efficiency of the model increases with increasing (i) expansion polynomial order, (ii) integration time and (iii) relative depth. Thus, the spectral/hp element technique appears to offers potentially significant savings in computational time for a fixed numerical error, compared to low-order numerical methods, for large-scale and long-time simulations of dispersive wave propagation. The practical applicability of the model is illustrated by several test cases.     

Stationary model of the upper ocean layer in terms of Ekman's theory

The impact of a downward heat flux on the depth of the upper ocean friction layer is investigated in terms of Ekman's theory for steady drift currents.Translated by V. Puchkin.     

Improved explicit approximation of linear dispersion relationship for gravity waves: A discussion

Recently, an accurate explicit approximation to linear dispersion relationship is proposed based on Eckart's explicit relationship (Beji, 2013). The author has nicely improved Eckart's explicit dispersion relationship by introducing an empirical correction function. The resulting expression is valid for the entire range of relative water depths and accurate to within 0.044%.In this discussion, the proposed expression by the author is simplified and improved to an accuracy of 0.019%. Moreover, a near exact solution with 0.001% accuracy is also given.     

全球海洋环流模式中上层海洋对表面强迫的响应和调整 I. 年际变率

利用一个较高分辨率的全球海洋环流模式在COADS1945～1993年逐月平均资料的强迫下对海温和环流场进行了模拟试验,研究了全球热带海洋(主要是热带太平洋)海温和环流场的年际变化特征及模式ENSO冷暖事件演变的控制机理.结果表明,模式成功地再现了和观测一致的海温和环流的年际变化以及ENSO演变特征.其中热带印度洋年际SST变率的主要模态表现为与ENSO相联系的海盆尺度的一致性增暖或变冷现象,次级模态为热带印度洋偶极子模态;热带大西洋的SST年际变率表现为类ENSO的年际振荡现象.在热带太平洋,SST年际变化主要表现为ENSO型,环流的年际变率表现为与ENSO相对应的热带海洋质量循环圈的年际异常.对应于暖(冷)事件,前期赤道海洋垂直环流圈显示出减弱(增强)的特征.其中南赤道流异常的位相较Nino3区海温总体要超前5个月左右的时间;赤道上翻流异常的位相在表层要超前4个月,并随时间由上至下扩展;赤道潜流的异常则显示出东传特征,其中最早的较为显著的异常发生ENSO成熟前3个月180°附近.在模式ENSO冷暖事件的演变过程中,次表层海温异常沿赤道的东传起了关键作用,模式的ENSO模态主要表现为"时滞振子"模态.