P-Vector inverse method evaluated using the modular ocean model (MOM)

Several major inverse methods (Stommel-Schott method, Wunsch method, and Bernoulli method) have been successfully developed to quantitatively estimate the geostrophic velocity at the reference level from hydrographic data. No matter the different appeance, they are based on the same dynamical sophistication: geostrophy, hydrostatic, and potential density (ϱ) conservation (Davis, 1978). The current inverse methods are all based on two conservation principles: potential density and potential vorticity (q=f∂ϱ/∂z) and require β-turning. Thus, two necessary conditions can be incorporated into any inverse methods: (1) non-coincidence of potential density and potential vorticity surfaces and (2) existence of vertical turning of the velocity (β-turning). This can be done using the P-Vector, a unit vector in the direction of ▽ϱ×▽q (Chu, 1994, 1995). The first necessary condition becomes the existence of the P-vector, and the second necessary condition leads to the existence of the P-vector turning in the water column. Along this line, we developed the P-vector inverse emthod with a pre-requirement check-up. The method was verified in this study using the Modular Ocean Model (MOM) from Pacanowskiet al. (1991) version of Bryan-Cox-Semtner ocean general circulation model (OGCM), which is based on the work of Bryan (1969). The statistically steady solutions of temperature and salinity from MOM are used as a “no-error data” set for computing absolute geostrophic velocities by the P-vector inverse method. Circulations are similar between the MOM statistically steady solutions and the P-vector solutions. Furthermore, the quantitative analysis shows that this inverse method has capability of picking up the major signal of the velocity field.     

Methane hydrate dissolution rates in undersaturated seawater under controlled hydrodynamic forcing

The dissolution of in-situ generated methane hydrate in undersaturated, synthetic seawater (S = 35) was investigated in a series of laboratory-based experiments at P-/T-conditions within the hydrate stability field. A controlled flow field was generated across the smooth hydrate surface to test if, in addition to thermodynamic variables, the dissolution rate is influenced by changing hydrodynamic conditions. The dissolution rate was found to be strongly dependent on the friction velocity, showing that hydrate dissolution in undersaturated seawater is a diffusion-controlled process. The experimental data was used to obtain diffusional mass transfer coefficients k_d, which were found to correlate linearly with the friction velocity, u. The resulting k_d/u-correlation allows predicting the flux of methane from natural gas hydrate exposures at the sediment/seawater interface into the bulk water for a variety of natural P, T and flow conditions. It also is a tool for estimating the rate of hydrate regrowth at locations where natural hydrate outcrops at the seafloor persist in contact with undersaturated seawater.     

Behavior of heavy metals and particulate matters in seawater expected from that of radioactive nuclides

Vertical profiles of²¹⁰Pb and²³⁰Th in the deep water were analyzed by using a simple one-dimensional model. Both nuclides are considered to settle down with the velocity of 1×10^–4 cm/sec. The diameter of particle corresponding to the velocity is calculated to be 5, while only about 10 % of the nuclides can be collected on filter paper with a pore size of 0.5. It is supposed that the nuclides exist in particulate materials which is changeable in size. This suggestion is supported from the following evidences. (1) The directly observed behavior of marine snow and the size distribution of particles observed with a coulter counter. (2) The existence of many chemical elements of which residence time is about 150 years. (3) Their possible existence as eutectic solid phases in the seawater. (4) The consistency of the observed accumulation rate of pelagic sediments with that estimated from the settling velocity. (5) The consistency of the decomposition rate of organic matter in the deep water estimated from the oxygen comsumption with that from the settling velocity.     

Sediment velocity estimation using iterative 3-D migrations of short offset seismic reflection data in deep water

In deep ocean settings where water depth greatly exceeds the source-to-receiver length, the geometry is insufficient for accurate determinations of velocity from reflection-moveout. However, velocities are crucial for estimates of physical properties and image processing. Focusing analyses with conventional post-stack two-dimensional migration improves images, but does not produce geologically meaningful velocities except in the special case of a two-dimensional earth. For the more general case of the three-dimensional earth there is no a priori method to determine the degree of geometrical complexity. We present a technique using a short-offset three-dimensional (3-D) data set over the 5 km deep trench west of the Lesser Antilles. These data illustrate highly sensitive post-stack 3-D focusing analyses (± 20 m s^–1 interval velocities), and the relationship of these seismically derived velocities to rock velocities. In our Barbados example we were able to establish the presence of a widespread 80-160 m thick low-velocity zone at and above the main low-angle fault. This observation suggests the water-rich décollement leaks water into the overlying sections. Also evident is a low-velocity section associated with turbidite sands. These results are confirmed with sparse logging data and well samples. Deep-water short offset 3-D experiments provide a potentially effective approach for velocity estimation, replacing the operational complexity of long-offsets with simpler short-offset techniques. In areas of structural complications and abundant diffracted energy, it is a surprisingly accurate method, utilizing the high fidelity 3-D wavefield and the information carried in zero-offset diffraction ellipsoids. The velocity used to properly collapse a diffraction ellipsoid is explicitly the velocity of propagation in the media since the travel path is known exactly. Thus, the derived velocities should closely represent rock velocities, unlike the 2-D case where the propagation geometry is not known.     

High resolution velocity analysis of ocean bottom seismometer data by theτ-p method

A method of high resolution seismic velocity analysis for ocean bottom seismometer (OBS) records is applied to the study of the shallow oceanic crust, especially sedimentary and basement layers. This method is based on the direct-p mapping and the-sum inversion. We use data obtained from a 1989 airgun-OBS experiment in the northern Yamato Basin, Japan Sea and derive P- and S-wave velocity functions that can be compared with the seismic reflection profiles. Using split-spread profile records, we obtain interface dips and true interval velocities from the OBS data. These results show good agreement with the reflection profile records, the acoustic velocities of core samples, and sonic log profiles. We also present a method for estimating errors in the derived velocity functions by calculating covariance of the derived layers' thicknesses. The estimated depth errors are about 150 m at shallow depths, which is close to the seismic wavelength used. The high resolution of this method relies on accurate determination of shot positions by GPS, spatially dense seismic observations, and the use of unsaturated reflected waves arriving after the direct water wave that are observed on low-gain component records.     

AOU as indicator of water-flow direction in the central North Pacific

Salinity, preformed phosphate and AOU distribution on the sigma-t surfaces of 26.8 and 27.3 of the central North Pacific Ocean, as well as the topography of these density surfaces, were studied. The direction of water flow suggested by the AOU distribution on these density surfaces was compared to that indicated by the acceleration potential contours on the _t =125 cl/ton and _t =80cl/ton surfaces drawn by Joseph L.Reid, Jr. The disagreements were explained in terms of mixing and possible gradients of primary production at the sea surface. On the 26.8 sigma-t surface a southward flow connecting the westward flow south of the Aleutian chain and the eastward flow farther south, between 175E and 180 is suggested by the AOU distribution but is not implied by the acceleration potential contours. If the circulation pattern at this density surface is similar to that at the sea surface, this southward flow is very likely to be real.     

基于Argo资料的北太平洋绝对地转流计算

针对海洋实测流速资料极其匮乏的事实, 利用2004 年1 月~2009 年12 月月平均的Argo 温盐格点资料, 结合改进的P-vector 方法重构北太平洋绝对地转流流场。与卫星高度计和实测流速的比较以及相关性分析表明, 重构的绝对地转流是可信的, 可以用来研究不同大尺度环流特征及其动力结构。可以为研究海洋动力过程和气候变化提供一套有用的流场数据。     

Analysis of wake behind a rotating propeller using PIV technique in a cavitation tunnel

A two-frame particle image velocimetry (PIV) technique is used to investigate the wake characteristics behind a marine propeller with 4 blades at high Reynolds number. For each of 9 different blade phases from 0° to 80°, 150 instantaneous velocity fields are measured. They are ensemble averaged to study the spatial evolution of the propeller wake in the region ranging from the trailing edge to one propeller diameter (D) downstream location. The phase-averaged mean velocity shows that the trailing vorticity is related to radial velocity jump, and the viscous wake is affected by boundary layers developed on the blade surfaces and centrifugal force. Both Galilean decomposition method and vortex identification method using swirling strength calculation are very useful for the study of vortex behaviors in the propeller wake region. The slipstream contraction occurs in the near-wake region up to about X/D=0.53 downstream. Thereafter, unstable oscillation occurs because of the reduction of interaction between the tip vortex and the wake sheet behind the maximum contraction point.     

Group velocity and energy transport by Rossby waves

A variational method gives a proof that, for normal modes of Rossby waves in the ocean on a-plane with depth variable in the meridional direction, the time average of the energy flux defined as the product of the pressure and the fluid-particle velocity agrees with the time average of the energy density multiplied by the group velocity. This conclusion holds in considerably general situations and sheds light on the controversial problem of the energy transport by Rossby waves. Moreover, the variational method gives a remarkable relation between the group velocity and the phase velocity. In particular this relation shows that the absolute value of the group velocity never exceeds that of the phase velocity.     

Vesteris Seamount: An enigma in the Greenland Basin

Vesteris Seamount is a solitary submarine volcano located at 73°30 N, 9°10W in the Greenland Basin. Steeply rising from a base depth of 3100 m to a minimum depth of ~ 130 m and striking 030°/210°, the feature lies ~ 300 km east of the east Greenland margin on an otherwise nearly flat and featureless seafloor. The main body of the seamount appears to have been formed episodically, the last of which culminated about 110 000 years ago. Subsequent, lower intensity volcanic activity continued sporadically until about 25 000 years ago, as evidenced by ash layers found in cores near the base of the feature. The smoothed surfaces at the summit make it likely that the seamount actually broached the surface during the Weichselian glacial period, between 8000 and 13 000 years ago. Two multibeam bathymetric investigations aboardPFS Polarstern during ARKTIS II/4 (1984) and ARKTIS VII/1 (1990), combined with geologic sampling, single-channel seismic profiling and underwater television coverage, have resulted in a new interpretation of both the morphology and origins of the seamount. Data collected aboardPolarstern from ARKTIS II/4 (1984) have been previously reported by Hempelet al. (1991), however, when combined with the ARKTIS VII/1 (1990) data set, a more detailed interpretation of the morphology and structure was feasible. This included the elongated shape of the feature and showed the existence of several small volcanic cones on the seamount flanks.The U.S. Government right to retain a non-exclusive, royalty free licence in and to any copyright is acknowledged.     

The roughness height and drag law over the water surface based on the hypothesis of local equilibrium

A series of measurements of winds and wind-waves were carried out in wind-wave flumes. A data analysis based on the hypothesis of local equilibrium yielded a new empirical formula on the controversial quantity of roughness heightz ₀ over the water surface: , where the nondimensional roughness height is defined bygz ₀/u _* ² and the wave-wind parameterũ byω _p u _*/g, g being the gravitational acceleration,u _* the friction velocity of air,ω _p the peak frequency of wind-wave spectra. The obtained formula is compared with Charnock's (1955) and Toba's (1979) proposals; is constant in the former and inversely proportional toũ in the latter. As in Toba's, this formula immediately leads to a practically important conclusion that the drag coefficientC _d depends not merely on the usual variableU ₁₀ (wind velocity at 10m height over the water surface), but also on the surface state represented by wind-waves. An explicit expression is provided for the drag coefficient incorporating the wave-wind parameter; it covers the range ofC _d calculated from most of the previous drag formulas, by varying the wave-wind parameter.     

A dynamic equation for the potential energy anomaly for analysing mixing and stratification in estuaries and coastal seas

In this paper, a time-dependent dynamic equation for the potential energy anomaly, , is rigorously derived from dynamic equations for potential temperature and salinity, the continuity equation and the equation of state for sea water. The terms locally changing are (A) the -advection, (B) the depth-mean straining, (C) the non-mean straining, (D) the vertical advection, (E) the vertical mixing, (F) surface and bottom density fluxes, (G) inner sources of density e.g. due to absorption of solar radiation and the non-linearity of the equation of state, and (H) horizontal divergence of horizontal turbulent density fluxes. In order to derive the equation in concise form, a vertical velocity (linearly varying with depth) with respect to depth-proportional vertical coordinates had to be defined. The evaluation of the terms in the -equation is then carried out for a one-dimensional tidal straining study and a two-dimensional estuarine circulation study. Comparisons to empirical estimates for these terms are made for the one-dimensional study. It is concluded that the -equation provides a general reference for empirical bulk parameterisations of stratification and mixing processes in estuaries and coastal seas and that it is a tool for complete analysis of the relevant terms from numerical models.     

Physical dynamics and biological implications of a mesoscale eddy in the lee of Hawai’i: Cyclone Opal observations during E-Flux III

E-Flux III (March 10–28, 2005) was the third and last field experiment of the E-Flux project. The main goal of the project was to investigate the physical, biological and chemical characteristics of mesoscale eddies that form in the lee of Maui and the Island of Hawai’i, focusing on the physical–biogeochemical interactions. The primary focus of E-Flux III was the cyclonic cold-core eddy Opal, which first appeared in the NOAA GOES sea-surface temperature (SST) imagery during the second half of February 2005. During the experiment, Cyclone Opal moved over 160 km, generally southward. Thus, the sampling design had to be constantly adjusted in order to obtain quasi-synoptic observations of the eddy. Analyses of ship transect-depth profiles of CTD, optical and acoustic Doppler current profiler (ADCP) data revealed a well-developed feature characterized by a fairly symmetric circular shape with a radius of about 80 km. Depth profiles of temperature, salinity and density were characterized by an intense doming of isothermal, isohaline and isopycnal surfaces. Isopleths of nutrient concentrations were roughly parallel to isopycnals, indicating the upwelling of deep nutrient-rich water. The deep chlorophyll maximum layer (DCML) shoaled from a depth of about 130 m in the outer regions of the eddy to about 60 m in the center. Chlorophyll concentrations reached their maximum values in Opal's core region (about 40 km in diameter), where nutrients were upwelled into the euphotic layer. ADCP velocity data clearly showed the cyclonic circulation associated with Opal. Vertical sections of tangential velocities were characterized by values that increased linearly with radial distance from near zero close to the center to a maximum of about at roughly 25 km from the center, and then slowly decayed. The vertical extent of the cyclonic circulation was primarily limited to the upper mixed layer, as tangential velocities decayed quite rapidly within a depth range of 90–130 m. Potential vorticity analysis suggests that only a relatively small (about 50 km in diameter) and shallow (to a depth of approximately 70 m) portion of the eddy is isolated from the surrounding waters. Radial movements of water can occur between the center of the eddy and the outer regions along density surfaces within an isopycnal range of σ–t_23.6 () and σ–t_24.4 (). Thus the biogeochemistry of the system might have been greatly influenced by these lateral exchanges of water at depth, especially during Opal's southward migration. While the eddy was translating, deep water in front of the eddy might have been upwelled into the core region, leading to an additional injection of nutrients into the euphotic zone. At the same time, part of the chlorophyll-rich waters in the core region might have remained behind the translating eddy and, thus contributed to the formation of an eddy wake characterized by relatively high chlorophyll concentrations.     

1995与1996年夏季琉球群岛两侧海流

基于1995,1996年夏季日本调查船的观测资料,采用P矢量方法对琉球群岛两侧的海流进行了计算.结果表明:黑潮为琉球群岛以西海域的一支东北向强流,1996年夏季的流速比1995年夏季的强,在深层出现南向逆流.黑潮东、西两侧分别存在一个反气旋式暖涡和一个弱的气旋式冷涡.1995年夏季,琉球群岛以东,从表层至以下层都存在一支沿岸北上的海流,即琉球海流.该海流来自黑潮分支,为本海区的一个主要物理特征.琉球海流以下出现弱的南向流.冲绳岛以东海域,在25°～25°30'N,128°30'～129°10'E附近从表层至700m水深存在一个中尺度的反气旋式暖涡.在温、盐水平分布图上,对应的出现一个较高温、低密水块.1996年夏季,冲绳岛西南海域存在一个中尺度的反气旋式暖涡和一个气旋式冷涡,形成一个偶极子,中间为较强的南向流,该现象为本海区的一个重要物理特征,属首次报道.冲绳岛以东表层主要被南向流控制,琉球海流不明显.200m以深在近岸出现北向流,这表明琉球海流的核心位于次表层.琉球海流的下面出现南向流.计算海区东北部从表层到700m水深出现一个中尺度的反气旋式暖涡,与1995年夏季时比较,其位置向北移动.此外在1996年夏季从近表层到深层,垂直方向和水平方向上的等温线、等盐线波动很大,例如在C断面上冷、暖涡相间出现,且暖     

Irrotational,progressive surface gravity waves near the limiting height

Numerical solutions of irrotational, progressive surface gravity waves in water of a constant depth are obtained by means of an iterative method. Our results suggest that waves with the surface slope angle greater than/6 may exist. The calculated phase velocity of deep water waves near the wave steepness 0.14 is significantly smaller than the value given by the Stokes' fourth approximation.In order to check our method, we apply it to the problem proposed byDavies (1951), which is hypothetical but similar to the present problem, and for which the exact solution is known. In this case our results show good agreement with the exact solution.     

Simulation of Seawater Intrusion Range in Coastal Aquifer Using the FEMWATER Model for Disaster Information

The FEMWATER model is a three-dimensional finite element model that can simulate both groundwater flow and mass transport at the same time, and so by using this FEMWATER model seawater intrusion was simulated and the results were compared to observation data. By using the observation data, a conceptual model was made of 3 layers based on the type of soil, and the hydrogeological parameter was selected. There are two types of dispersivities: longitudinal dispersivity and transverse dispersivity. Longitudinal dispersivity was calculated by Xu and Eckstein's method (1995 Xu , M. and Y. Eckstein . 1995 . Use of weighted least-squares method in evaluation of the relationship between dispersivity and field scale . Ground Water 33 ( 6 ): 905 – 908 .[Crossref], [Web of Science ®] , [Google Scholar]) and the transverse dispersivity was put at 3/10 of the longitudinal dispersivity for the simulation. The electrical conductivity suggested by Freeze and Cherry (1979 Freeze , R. A. and J. A. Cherry , J. A. 1979 . Groundwater . Eaglewood Cliffs , NJ : Prentice-Hall . [Google Scholar]) and the correlation with total dissolved solids, TDS, was used to classify the state of the water. Upon analyzing the depths at 15 m, 30 m, 60 m, and 90 m, occurrence of seawater intrusion near the ground surface was limited to the vicinity of the coastline, but as the depths became deeper, seawater intrusion expanded inland. From these results we can come to the conclusion that in order to prevent seawater intrusion there needs to be sufficient preliminary examination.     

An improved “benchmark” method for estimating particle settling velocities from time-series sediment trap fluxes

A new method, based on fitting Fourier series to time-series (TS) data from sediment traps, was developed to estimate the settling velocities (SVs) of sinking particles in the open ocean. This new method was applied to data from MedFlux, as well as from the US JGOFS NABE, EqPac, and ASPS studies. Fluxes of mass and of four chemical tracers, as well as the molar ratios of the latter, were plotted on logarithmic scales; Fourier series were then fit to these data. In each case we determined the most likely settling velocity using a likelihood-based nonlinear fitting algorithm. Variation among estimates using single tracers was significantly less than variation using tracer ratios; we therefore concluded that estimates based on single tracers are to be preferred to estimates based on tracer ratios. Our results also showed no obvious differences among SVs estimated using different single tracers. The best estimate of settling velocity using single-tracer fluxes with good temporal resolution (i.e. for sites with cup rotation times 8.5 days) is 205 m/d, with standard deviation 74 m/d. For MedFlux data alone (which have a resolution of 4–6 days), the estimate is 220±65 m/d. This latter value is within 10% of the estimate of average settling velocity (242±31 m/d) made using MedFlux IRS traps in “settling velocity” mode.     

Error of digitization of multidimensional random hydrophysical fields with power spectra

Exact and approximate formulae are proposed to evaluate the errors of spatial/temporal digitization (aliasing) of then-dimensional random fields withm-power spectra atn=1, , 4 and 1m5, which satisfactorily approximate temperature, velocity, density, and other fields in the sea and atmosphere.Translated by Vladimir A. Puchkin.     

The mechanism of the development of wind-wave spectra

The mechanism of the development of wind-waves will be proposed on the basis of the observed wave spectra in the wind tunnels and at Lake Biwa (Imasato, 1976). It consists of two aspects: One is that the air flow over the wind-waves transfers momentum concentratively to the steepest component waves and the other is that the upper limit of the growth of a wave spectral density is given by the ultimate value in the slope spectral density. The first aspect means that the wave field has the momentum transfer filter on receiving the momentum from the air flow. Wind-waves in the stage of sea-waves receive the necessary amount of momentum by the form drag,e.g. according to the Miles' (1960) inviscid mechanism, through a very narrow frequency region around a dominant spectral peak. On the other hand, wind-waves in the stage of initial-wavelets receive it according to the Miles' (1962a) viscous model through a fairly broad frequency region around the peak. The upper limit ofS _max developing according to viscous mechanism is given byS _max =6.40×10^–4 k _max ^–2cm²s andS _max =2.03C(f _max )^–2cm²s(S _max is the power density of the wave spectral peak with the frequencyf _max ,k _max is the wave number corresponding to the frequencyf _max andC is the phase velocity).From the second aspect, the upper limit of the growth of wave spectral density is given by 33.3f ^–4cm²s in the frequency region of late stage of sea-waves. Therefore, the spectral peak, which has the largest value in the slope spectral density in the component waves of the wave spectrum, rises high over the line 4.15f ^–5cm²s. The energy is transported from the spectral peak to the high frequency part and to the forward face of a wave spectrum by nonlinear wave-wave interaction. This nonlinearity is confirmed by the bispectra calculated from the observed wind-wave data. In the stage of sea-waves, nonlinear rearrangement of the wave energy comes from a narrow momentum transfer filter, and, in the stage of initial-wavelets, it comes mainly from small corrugations and small steepness of the wave field.