A piecewise curve-fitting technique for vertical oceanographic profiles and its application to density distribution

A unified method of approximation, extrapolation, and objective layering is offered for processing vertical oceanographic profiles. The method is demonstrated using seawater density and consists of adjustable splitting of each individual profile into N vertical layers based on tentative, piecewise linear homogeneous approximation with specified accuracy and a final fitting of an N-layered analytical model to data. A set of 3N coefficients of the model includes one density value at the sea surface; N−1 depths of layer interfaces; and N pairs of coefficients that describe a profile shape within the n-th layer—an asymptotic density value (a key parameter for extrapolation) and a vertical scale of maximum density variability (related to vertical gradient). Several distinctive characteristics of the technique are: (1) It can be used for the analysis of the vertical structure of individual profiles when N is an unknown parameter, and spatial interpolation when N should be equal for all profiles. (2) A justified downward extrapolation of incomplete data is possible with the model, especially if historical deepwater profiles are available. (3) Layer interfaces, as well as other coefficients, are derived with only one fitting to the entire profile. (4) The technique, using its general formulation, can serve as a parent for developing various types of models. The simpler step-like (with hyperbolic or exponential approximation) and more complicated smooth (continuous in gradient space) models were designed and tested against a large number of density profiles from the Sea of Okhotsk and the Gulf of California. Comparison of parametric, z-levels and isopycnal averaging was done for the region off the northeastern coast of Sakhalin.     

Numerical modeling of the arctic ocean deepwater formation: Part I. Idealized study

In this study the capability of numerical models to transport anomalously dense bottom waters from the areas of their generation on the shelf to deepwater areas is considered. These numerical models have a number of mechanisms that are favorable for this transport, but unrelated directly with the real physics of that motion. A more accurate account of the resupply of oceanic deep waters requires either a very fine resolution for the numerical model or an efficient parametrization. Two such parametrizations are proposed in this study. The first is based on amplified lateral diffusion exchange if slope instability exists. The second is related to an additional procedure allowing dense waters to leak under waters of lesser density along a shelf slope with their displacement to higher layers. The first parametrization is shown to be less efficient because it leads to a significant interaction between anomaly waters and ambient waters and to a quick loss of their original features. The second procedure is more preferable because it leads to a significant resupply of dense deep waters at the expense of dense waters generated at a distance significantly far from the shelf break. With the help of the given parametrization, during a numerical experiment these waters can be transported considerable distances away from their place of generation to the shelf boundary and deepwater areas of the ocean with a slight change in their original density characteristics.     

Study of the Neutral Ekman Flow Using an Algebraic Reynolds Stress Model

A recently developed fully explicit algebraic model of Reynolds stress and turbulent heat flux in a thermally stratified planetary atmospheric boundary layer without stratification has been used for a numerical study of the Ekman turbulent boundary layer over a homogeneous rough surface for different dimensionless surface Rossby numbers. A comparative analysis has been conducted for a closure model of the transport term in the prognostic equation of turbulent kinetic energy dissipation including third-order moments. Dependences of the total wind rotation angle on the Rossby number have been obtained. The calculated vertical profiles of mean velocity, turbulent stress, turbulent kinetic energy, surface-friction velocity, and boundary-layer height agree satisfactorily with observational and earlier obtained LES data.     

黄、东海边界悬浮物输运量的卫星遥感估算

利用2002—2015年的中等分辨率成像光谱仪(Moderate Resolution Imaging Spectroradiometer, MODIS)第4波段反射率产品对黄、东海边界的悬浮物输运情况进行了分析, 同时利用黄、东海2003年现场实测的悬浮物含量数据和这期间准实时的MODIS影像数据, 建立了基于MODIS第4波段遥感反射率的悬浮物含量遥感反演模型。基于该模型, 反演黄、东海表层悬浮物浓度, 并计算悬浮物扩散面积。利用反演得到的表层悬浮物浓度和悬浮物浓度垂向预测模型, 获取垂直方向上的悬浮物浓度, 同时结合水深、扩散面积数据, 计算特定时刻水体中的悬浮物含量。利用插值方法获得整个输运过程中各时刻水体中的悬浮物的含量, 计算各个时刻悬浮物含量的总和作为悬浮物的输送量。对2002年10月至2003年4月黄海向东海的悬浮输送量进行估算, 估算结果为153Mt。     

Large-scale vertical hydroacoustic structure of the North Atlantic and its seasonal variability

A classification based on the number and types of large-scale acoustic waveguides is proposed for the mean seasonal profiles of sound speed propagation. A scheme for North Atlantic zoning, using typical curves of the sound speed vertical distribution, is given. The channel axis's position is shown not to depend on the water mass haline properties, being controlled by the temperature field vertical stratification.Translated by Vladimir A. Puchkin.     

Influence of chemical composition and microstructure on the hygroscopic growth of pyrogenic aerosol

The hygroscopic properties of a model smoke aerosol series including ammonium sulfate, levoglucosan, and oxalic and humic acids were investigated. The condensation growth of particles and the relative humidity of direct and inverse particle-droplet phase transitions were measured with a tandem differential mobility analyzer. The data were obtained in the range of relative humidity from 4 to 95% at 298 K. Experimental growth factors were compared to the results of Koehler theory by using contemporary models of water activity in particles. It was shown that the hygroscopic properties of the aerosols were substantially different from the hygroscopic properties of the corresponding bulk phases. It was found that a complex microstructure of nanoparticles resulted in excess volume absorption of water, in some cases causing their full deliquescence at lower relative humidities than in macrosystems of identical composition. From the measurement data, a parametrization of hygroscopic particle growth was performed. The results of parametrization can be used to construct refined models of a smoke aerosol.     

Assimilating satellite SST/SSH and in-situ T/S profiles with the Localized Weighted Ensemble Kalman Filter

The Localized Weighted Ensemble Kalman Filter(LWEnKF) is a new nonlinear/non-Gaussian data assimilation(DA) method that can effectively alleviate the filter degradation problem faced by particle filtering, and it has great prospects for applications in geophysical models. In terms of operational applications, along-track sea surface height(AT-SSH), swath sea surface temperature(S-SST) and in-situ temperature and salinity(T/S) profiles are assimilated using the LWEnKF in the northern South China ...     

Second-order moment advection scheme applied to Arctic Ocean simulation

We apply the second-order moment (SOM) advection scheme of (Prather, M.J. 1986. Numerical advection by conservation of second-order moments. J. Geophys. Res. 91, 6671–6681.) to the simulation of the large-scale circulation of the Arctic Ocean with a coupled ocean–sea-ice model. Compared to three other advection schemes commonly employed in ocean simulations (centred differences, flux corrected transport, and multidimensional positive definite advection transport), the SOM method helps preserve the vertical structure of Arctic water masses. The depth, thickness and hydrographic properties of the Arctic Surface Water and the Arctic Atlantic Layer are better represented with SOM than with any of the other three advection algorithms. We also present a convenient method for calculating the implicit numerical diffusivity of upstream based schemes, such as the SOM method, and discuss three approaches for improving the monotonicity properties of the SOM algorithm.     

Predicting Bathymetry from the Earth's Gravity Gradient Anomalies

Historically, prediction of ocean floor depth, or bathymetry, has been based on the isostatic modeling and linearized relationships between gravity anomalies and bathymetry. The need for isostatic modeling limits the application of the resulting bathymetry predictions as constraints in geophysical models. An alternative technique making use of the Earth's vertical gravity gradient for predicting bathymetry is explored in this paper. This technique is based on the fact that the observed gravity gradient anomalies result primarily from local mass concentrations on the ocean floor, and that mass compensation by the oceanic crust has an insignificant effect on the gravity gradients, and can be neglected. The resulting bathymetry prediction therefore is independent of isostatic modeling assumptions, allowing it to be used as a constraint on models of lithospheric compensation and for other geodetic and geophysical applications.     

Convection induced by cooling at one side wall in two-dimensional non-rotating fluid —applicability to the deep Pacific circulation—

In recognition that similarity in the density balance leads to resemblance in circulation between the two-dimensional non-rotating and three-dimensional rotating systems which have similar density stratification, we investigate convection induced by cooling at one side wall and heating at the sea surface by using a two-dimensional non-rotating model as idealized representation for the deep Pacific circulation. In the model, various vertical profiles are taken for the side wall cooling, which are assumed to correspond to the density structure of the Anatarctic Circumpolar Current. In a small diffusivity range, two important features are found to be robust against change in the vertical profile of the side wall cooling. One is that the density stratification is horizontally almost uniform. The other is that the balance in the density equation between the vertical advection and the vertical diffusion holds in the interior. Consequently, the vertical density balance, together with the equation of continuity, determines the circulation pattern for the prescribed vertical profile of the side wall cooling. The multi-layered meridional flow, which is expected to exist in the deep Pacific, is shown to form for certain vertical profiles of the side wall cooling.     

Predicting Bathymetry from the Earth's Gravity Gradient Anomalies

Historically, prediction of ocean floor depth, or bathymetry, has been based on the isostatic modeling and linearized relationships between gravity anomalies and bathymetry. The need for isostatic modeling limits the application of the resulting bathymetry predictions as constraints in geophysical models. An alternative technique making use of the Earth's vertical gravity gradient for predicting bathymetry is explored in this paper. This technique is based on the fact that the observed gravity gradient anomalies result primarily from local mass concentrations on the ocean floor, and that mass compensation by the oceanic crust has an insignificant effect on the gravity gradients, and can be neglected. The resulting bathymetry prediction therefore is independent of isostatic modeling assumptions, allowing it to be used as a constraint on models of lithospheric compensation and for other geodetic and geophysical applications.     

An adjoint method for the assimilation of statistical characteristics into eddy-resolving ocean models

The study investigates perspectives of the parameter estimation problem with the adjoint method in eddy-resolving models. Sensitivity to initial conditions resulting from the chaotic nature of this type of model limits the direct application of the adjoint method by predictability. Prolonging the period of assimilation is accompanied by the appearance of an increasing number of secondary minima of the cost function that prevents the convergence of this method. In the framework of the Lorenz model it is shown that averaged quantities are suitable for describing invariant properties, and that secondary minima are for this type of data transformed into stochastic deviations. An adjoint method suitable for the assimilation of statistical characteristics of data and applicable on time scales beyond the predictability limit is presented. The approach assumes a greater predictability for averaged quantities. The adjoint to a prognostic model for statistical moments is employed for calculating cost function gradients that ignore the fine structure resulting from secondary minima. Coarse resolution versions of eddy-resolving models are used for this purpose. Identical twin experiments are performed with a quasigeostrophic model to evaluate the performance and limitations of this approach in improving models by estimating parameters. The wind stress curl is estimated from a simulated mean stream function. A very simple parameterization scheme for the assimilation of second-order moments is shown to permit the estimation of gradients that perform efficiently in minimizing cost functions.     

Wave forces on, and water-surface fluctuations around a vertical cylinder encircled by a perforated square caisson

The wave forces and moments on and the water surface fluctuations around a vertical circular cylinder encircled by a perforated square caisson were experimentally investigated. The porosity of the outer square caisson was varied from 4.24 to 14.58%. The in-line wave forces on the inner vertical cylinder are influenced by changing the porosity of the outer caisson, whereas the variations in the water surface fluctuations are less influenced in this porosity range. The in-line moment on the vertical cylinder is relatively less sensitive when the porosity is increased from 4.24 to 8.75%, but varies substantially when it is increased from 8.75 to 14.58%. The force and moment ratio (i.e. the ratio of the force or moment on the vertical cylinder, when it is encircled by the perforated caisson to the force or moment on the cylinder without any protection around it) reduces with increased wave height, H, and wave length, L, whereas the wave height ratio (ratio of the wave height at a point in the vicinity of the structure to the incident wave height) is less sensitive for the varying H and L. A new non-dimensional parameter, p^1.5 (D/L)/(H/d), is introduced to predict the in-line force and moment on the inner vertical cylinder, where d is local water depth, D is the diameter of the inner cylinder and p is the porosity of the outer caisson in percentage. Simple predictive equations for forces, moments and water surface fluctuations are provided.     

Multichannel moment detectors for transients in shallow-water noise

This paper presents an evaluation of second, third, and fourth-order moments for the passive detection of transient signals in both simulated Gaussian noise and measured noise. The measured noise was recorded by a vertical array located near the San Diego, CA, harbor and is dominated at low frequencies by ship-generated noise. The detectors assume neither noise nor signal stationarity and can use single or multiple channels of data. Simulation results indicate that the fourth-order moment detector often performs better than the energy detector in the correlated measured noise, with increasing channel contributions to the moment function, resulting in increased gain. The results in simulated Gaussian noise likewise favor the fourth-order moment detector, at least for the signals with significant fourth-order moments, but the ability of the higher order detector to discriminate against correlated noise is evident. Analysis over a 30-min segment of the measured noise with selected signals demonstrates that fourth-order detection gains can be reliably expected as the noise statistics change.     

北部湾南部重力柱状样的MSCL地声学性质测量及分析

利用钻孔多参数连续记录仪(multi-sensorcorelogger,MSCL),对在南海北部湾南部海区所取得的6个孔重力柱状样进行了测量,获得了连续P波波速、湿密度和孔隙度数据。在室内实验室对柱状样分样之后的沉积物样品进行了孔隙度与湿密度的测定。利用不同的统计回归方法对所获得的6个柱状样的孔隙度与P波波速进行相关分析,并对比室内测量的P波波速、湿密度、孔隙度,建立了基于孔隙度数据对沉积物P波波速进行预测的方法,对MSCL测量方法的优缺点进行了讨论。结果表明,MSCL测量结果与实验室柱状样测定结果较为吻合。样品所在深度对声速的变化影响不大。孔隙度与声速的多项式回归,样条插值回归和GAM模型回归都获得了较高的相关系数,GAM模型可以提供一个较为接近测量值的声速预测方法。MSCL用来测量海底沉积物,可以获得大量的高密度、高精度的沉积物地声学及其他参数数据,但是,如果空气混入沉积物样品中则会导致MSCL测量结果失真。该研究为使用MSCL在区域海底地声学性质、恢复区域海洋沉积历史、海底地声学模型建立等研究提供了参考。     

Probabilistic Models for the Probability of Wave Breaking and Whitecap Coverage Based on Kinematic Breaking Criterion

More and more researches show that neither the critical downward acceleration nor the critical slope of water waves is a universal constant. On the contrary, they vary with particular wave conditions. This fact moders the models either for the probability of wave breaking B or for the whitecap coverage W based on these criteria difficult to apply. In this paper and the one which follows we seek to develop models for the prediction of both B and W based on the kinematical criterion. First, several joint probabihstic distribution functions (PDFs) of wave characteristics are derived, based on which the breaking properties B and W are estimated. The estimation is made on the assumption that a wave breaks ff the horizontal velocity of water particles at its crest exceeds the local wave celerity, and whitecapping occurs in regions of fluid where water particles travel faster than the waves. The consequent B and W depend on wave spectral moments of orders 0 to 4.Then the JONSWAP spectrum is used to represent the fetch-limited sea waves in deep water, so as to relate the probahility of wave breaking and the whitecap coverage with wind parameters. To this end, the time-averaging technique proposed by Glazman (1986) is applied to the estimation of the spectral moments involved, and furthermore, the theoretical models are compared with available observations collected from published literature. From the comparison, the averaging time scale is determined. The final models show that the probability of wave breaking as well as the whitecap coverage depends on the dimensionless fetch. The agreement between these models and the database is reasonable.     

Comparison of the nonlocal transport characteristics of a series of one-dimensional oceanic boundary layer models

Seven one-dimensional oceanic boundary layer models are investigated to assess the possible nonlocal transport characteristics of mass and heat in the upper ocean. The dynamical models have been chosen from the diffusion and bulk types currently in use plus two modifications of the transilient type that have been used extensively for atmospheric work by Stull. The models are forced using wind speed and insolation conditions recorded during a 9-day oceanographic cruise near Bermuda in March 1993 during the decline of the spring bloom. The attenuation of sunlight in the upper ocean is calculated using a full spectral model for downwelling irradiance. The vertical heat transport characteristics are reported and compared. A series of spectral diagnostic tests (Green's function analysis, process spectra, and overall mixing lengths) reveal significant differences in the vertical transport characteristics of the models that are not observed in commonly used diagnostics such as sea surface temperature or mixed layer depth. Age spectra (or modal time since last surface contact) are calculated from Green's functions and reveal how a water mass can be cut off from the surface within a short time period. The large differences in the vertical mixing characteristics of the upper ocean boundary layer models have potential implications on the vertical distributions of short-lived chemical tracers and phytoplankton.     

A temperature gradient column for the rapid measurement of the density of brines

A vertical temperature-gradient column for measuring the density of brine is described. The column is filled with a liquid which is immiscible with the brine and in which a brine droplet comes to rest at a point where its own density matches that of the medium. This method has been used for determination of the density profiles in the interstitial waters of sediment cores from the bottom of the Dead Sea and for other density measurements where rapidity and ease of measurement was required.     

Measurement of water particle velocities in waves

Water surface profiles and horizontal and vertical water particle velocity components have been measured to investigate the properties of intermediate depth waves generated in the laboratory. The data has been compared with linear wave theory. It was found that linear theory predicted the attenuation of velocity field with depth successfully and that it overestimates both components of velocity slightly.