ARMA order model detection using minimum of Kurtosis: application on seismic data

The objective of this study is to find the order and coefficients of non-low-phase causal filters for ARMA (auto regressive moving average) filter model, using the Kurtosis minimization criterion. This method is based on the Kurtosis calculation of the treated sample at the input level and its identification at the output of the ARMA model. For this purpose, the order and coefficients of the AR (auto regressive) part are identified using the Yule-Walker algorithm at order two and then extended to order four. To obtain the MA (moving average) part, the AR components are calculated at first from the ARMA filter by deconvolution. Then, spectrally equivalent and minimum phase (SEMP) MA filter is identified using the Durbin algorithm at second and fourth order. Finally, the correct filter is found when the Kurtosis value of the output ARMA filter reconstituted is the closest to the Kurtosis of introduced signal. The proposed method is then tested on simulated processes and applied to real seismic data to perform blind deconvolution and obtain the reflectivity coefficients of subsoil studied.     

Phytoplankton biomass and size fractions in surface waters of the Australian sector of the southern ocean

We collected surface water along the 142nd E meridian from Tasmania to Antarctica in December 1999. We measured temperature, salinity and total chlorophyll a; additionally, we collected suspended particle size fractions and used fluorometric analysis to determine the quantity of chlorophyll a in each of four cell size classes: picoplankton (<3 μm), two nanoplankton fractions (3–10 μm and 10–20 μm) and microplankton (> 20 μm). Changes in temperature and salinity show that we crossed 6 water masses separated by 5 fronts. We found low abundance (<0.2 mg m⁻³) of chlorophyll in all size classes, with the exception of higher values near the continent (0.2 to 0.4 mg m⁻³). Lowest chlorophyll values (<0.1 mg m⁻³) were found in the Polar Frontal Zone (51° to 54°S). Microplankton made up the largest portion of total chlorophyll throughout most of the region. We conclude that biomass of all phytoplankton fractions, especially pico-and nanoplankton, was constrained by limiting factors, most probably iron, throughout the region and that ecosystem dynamics within a zone are not circumpolar but are regionalized within sectors.     

Parametric vibration of submerged floating tunnel tether under random excitation

For the study of the parametric vibration response of submerged floating tunnel tether under random excitation,a nonlinear random parametric vibration equation of coupled tether and tube of submerged floating tunnel is set up.Subsequently,vibration response of tether in the tether-tube system is analyzed by Monte Carlo method.It may be concluded that when the tube is subjected to zero-mean Gaussian white noise random excitation,the displacement and velocity root mean square responses of tether reach the peak if the circular frequency of tube doubles that of tether;the displacement and velocity root mean square responses of tether increase as the random excitation root mean square increases;owing to the damping force of water,the displacement and velocity root mean square responses of tether decrease rapidly compared with tether in air;increasing the damping of the tether or tube reduces the displacement and velocity root mean square responses of tether;the large-amplitude vibration of tether may be avoided by locating dampers on the tether or tube.     

Study on scattering wave force of horizontal and vertical plate type breakwaters

The interaction between wave and horizontal and vertical plates is investigated by the boundary element method,and the relations of wave exciting force with plate thickness,submergence and length are obtained.It is found that:1) The efficient wave exciting force exists while plate submergence is less than 0.5 m,and the plate is very thin with order O(0.005 m).2) The maximum heave wave exciting force exists,and it is the main factor for surface and submerged horizontal plate while the roll force can be ignored.3) The maximum sway wave exciting force exists,it is the main factor for surface or submerged vertical plate,and the roll force is about 20 times of horizontal plate.     

SAR-imaged spiral eddies in Mutsu Bay and their dynamic and kinematic models

A pattern of slick streaks winding into a spiral, known as a spiral eddy, was identified in 5 images taken by the ERS-1/2 synthetic aperture radar (SAR) in Mutsu Bay (Japan); dynamic and kinematic models of these spiral eddies have been proposed. Common characteristics of the five spiral eddies are: 1) an eddy diameter of about 15 km; 2) their location in the western part of the bay; and 3) their cyclonic direction of rotation. Moreover, the wind conditions over the bay were common: prior to acquiring the images, a strong easterly wind continued blowing for more than one day. The wind field on the bay is known to be orographically steered and has strong windstress vorticity, which generates cyclonic circulation. The diameter and location of the circulation simulated with a numerical ocean model corresponded well to those of the identified spiral eddies. Based on these facts, we propose a dynamic model for the movement of a slick streak, and a kinematic model for the formation of a spiral eddy. We have assumed calm air, a microlayer and seawater with a cyclonic circulation in the dynamic model. The balance of forces is established in the microlayer among the frictional force from the seawater, the frictional force from the calm air, the gravitational force, and the Coriolis force. As a result, the velocity vector of the microlayer deflects slightly towards the center of the cyclonic circulation. We have assumed a point source of the microlayer in the kinematic model. The shapes of a slick streak simulated with the models agree well with the identified patterns in the SAR images.     

Reevaluation of historical ocean heat content variations with time-varying XBT and MBT depth bias corrections

As reported in former studies, temperature observations obtained by expendable bathythermographs (XBTs) and mechanical bathythermographs (MBTs) appear to have positive biases as much as they affect major climate signals. These biases have not been fully taken into account in previous ocean temperature analyses, which have been widely used to detect global warming signals in the oceans. This report proposes a methodology for directly eliminating the biases from the XBT and MBT observations. In the case of XBT observation, assuming that the positive temperature biases mainly originate from greater depths given by conventional XBT fall-rate equations than the truth, a depth bias equation is constructed by fitting depth differences between XBT data and more accurate oceanographic observations to a linear equation of elapsed time. Such depth bias equations are introduced separately for each year and for each probe type. Uncertainty in the gradient of the linear equation is evaluated using a non-parametric test. The typical depth bias is +10 m at 700 m depth on average, which is probably caused by various indeterminable sources of error in the XBT observations as well as a lack of representativeness in the fall-rate equations adopted so far. Depth biases in MBT are fitted to quadratic equations of depth in a similar manner to the XBT method. Correcting the historical XBT and MBT depth biases by these equations allows a historical ocean temperature analysis to be conducted. In comparison with the previous temperature analysis, large differences are found in the present analysis as follows: the duration of large ocean heat content in the 1970s shortens dramatically, and recent ocean cooling becomes insignificant. The result is also in better agreement with tide gauge observations. On leave from the Meteorological Research Institute of the Japan Meteorological Agency.     

Equatorial Pacific gravity lineaments: interpretations with basement topography along seismic reflection lines

The central equatorial Pacific is interesting for studying clues to upper mantle processes, as the region lacks complicating effects of continental remnants or major volcanic plateaus. In particular, the most recently produced maps of the free-air gravity field from satellite altimetry show in greater detail the previously reported lineaments west of the East Pacific Rise (EPR) that are aligned with plate motion over the mantle and originally suggested to have formed from mantle convection rolls. In contrast, the gravity field 600 km or farther west of the EPR reveals lineaments with varied orientations. Some are also parallel with plate motion over the mantle but others are sub-parallel with fracture zones or have other orientations. This region is covered by pelagic sediments reaching ~?500–600 m thickness so bathymetry is not so useful for seeking evidence for plate deformation across the lineaments. We instead use depth to basement from three seismic reflection cruises. In some segments of these seismic data crossing the lineaments, we find that the co-variation between gravity and basement depth is roughly compatible with typical densities of basement rocks (basalt, gabbro or mantle), as expected for some explanations for the lineaments (e.g., mantle convection rolls, viscous asthenospheric inter-fingering or extensional deformation). However, some other lineaments are associated with major changes in basement depth with only subtle changes in the gravity field, suggesting topography that is locally supported by varied crustal thickness. Overall, the multiple gravity lineament orientations suggest that they have multiple origins. In particular, we propose that a further asthenospheric inter-fingering instability mechanism could occur from pressure variations in the asthenosphere arising from regional topography and such a mechanism may explain some obliquely oriented gravity lineaments that have no other obvious origin.     

Experimental Study on Sloshing Characteristics in the Elastic Tank Based on Morlet Wavelet Transform

Hydroelastic effect of sloshing is studied through an experimental investigation. Different excitation frequencies are considered with low-fill-depth and large amplitude. Morlet wavelet transform is introduced to analyze the free surface elevations and sloshing pressures. It focuses on variations and distributions of the wavelet energy in elastic tanks. The evolutions of theoretical and experimental wavelet spectra are discussed and the corresponding Fourier spectrums are compared. Afterwards, average values of the wavelet spectra are extracted to do a quantitative study at various points. From the wavelet analysis, sloshing energies are mainly distributed around the external excitation frequency and expanded to high frequencies under violent condition. In resonance, experimental wavelet energy of the elevation in elastic tanks is obviously less than that in the rigid one; for sloshing pressures, the elastic wavelet energy close to the rigid one and conspicuous impulse is observed. It recommends engineers to concern the primary natural frequency and impulsive peak pressures.     

Dissolved inorganic carbon (DIC) and its carbon isotopic composition in sediment pore waters from the Shenhu area,northern South China Sea

The Shenhu area is one of the most favorable places for the occurrence of gas hydrates in the northern continental slope of the South China Sea. Pore water samples were collected in two piston cores (SH-A and SH-B) from this area, and the concentrations of sulfate and dissolved inorganic carbon (DIC) and its carbon isotopic composition were measured. The data revealed large DIC variations and very negative δ ¹³C-DIC values. Two reaction zones, 0–3 mbsf and below 3 mbsf, are identified in the sediment system. At site SH-A, the upper zone (0–3 mbsf) shows relatively constant sulfate and DIC concentrations and δ ¹³C-DIC values, possibly due to bioturbation and fluid advection. The lower zone (below 3 mbsf) displays good linear gradients for sulfate and DIC concentrations, and δ ¹³C-DIC values. At site SH-B, both zones show linear gradients, but the decreasing gradients for δ ¹³C-DIC and SO₄ ²⁻ in the lower zone below 3 mbsf are greater than those from the upper zone, 0–3 mbsf. The calculated sulfate-methane interface (SMI) depths of the two cores are 10.0 m and 11.1 m, respectively. The depth profiles of both DIC and δ ¹³C-DIC showed similar characteristics as those in other gas hydrate locations in the world oceans, such as the Blake Ridge. Overall, our results indicate an anaerobic methane oxidation (AMO) process in the sediments with large methane flux from depth in the studied area, which might be linked to the formation of gas hydrates in this area.     

Optimization of dissolved urea measurements in coastal waters with the combination of a single reagent and high temperature

Urea is an unstable and intermediate organic nitrogenous compound present in coastal environments and is derived from the excretion of some aquatic organisms and wastewater discharges. Urea plays an important role in the nitrogen cycle, where it is utilized by algae, including diatoms. However, there are very limited relevant data on the production, consumption, and degradation of urea because of the lack of appropriate measurement techniques. The conventional method is based on the formation of a colored product when urea reacts with diacetyl monoxime in a sulfuric acid solution. We examined the optimal conditions for the formation of the colored product; specifically, we evaluated different temperatures (22–80 °C), reaction times, mixing ratios of color reagents, and sample storage times. Application of the single mixed color-developing reagent (COLDER) at 70 °C resulted in the optimal formation of the colored product within a short reaction time of 60 min. This method was then used to measure dissolved urea in different coastal environments. The concentrations detected were as follows: 0.65–0.72, 0.49–0.58, and 1.09–2.28 µM urea-N at coral reef, seagrass, and mangrove sites, respectively. Our results showed high precision (SD = 0.02, CV = 1.2%), a low detection limit (0.03 µM urea-N), and a high recovery rate (94–99%). In summary, this high-temperature procedure for urea measurements should be valuable for obtaining high-precision data that can further the understanding of urea dynamics and its role in coastal ecosystems.