Physical and numerical modeling of reflections from interfaces and thin layers

Basic models for understanding the reflection of seismic waves by the seafloor or the sublayering consist of isolated interfaces and isolated thin layers. An isolated interface model is used to demonstrate a reasonable agreement between physical modeling and two numerical solutions of the elastic wave equation. An isolated thin-layer model is used to show that reflections from it can be regarded as a scattered wavefield caused by three secondary sources proportional toDeltaM,-Deltamu, and,Deltarhoat each element of the layer, whereDelta, indicates the difference between the layer and the surroundings,Mis the compressional wave modulus,muis the shear wave modulus, andrhois the density. This viewpoint leads to a simple explanation of the observed variations in amplitude and phase of the reflected waves as the offset distance from source to receiver increases.     

Tsunami detectability using open-ocean bottom pressure fluctuations

Rationale for the measurement of open-ocean tsunami signatures are presented, and available pertinent data are reviewed. Models for tsunami signature and background noise are proposed in order to synthesize an optimum tsunami receiver. Using these models, the minimum tsunami amplitude (in cm) to yield the probability of correct tsunami detectionP_{D} = 0.999and probability of false alarmP_{F} = 10^{-3}is found to be0.718/sqrt{f_{0}}, wheref_{0}is the tsunami dominant frequency (in cycles/h). A realizable receiver is proposed and its performance is evaluated using actual tsunami signatures. It is demonstrated that the detection of a tsunami with an average amplitude as small as 0.7 cm is possible for theP_{D}andP_{F}as above. Simulation results using synthesized background noise are shown. Tidal effects on the receiver performance also are considered and are found negligible for a certain range of the receiver parameters, resulting in a considerable reduction of the signal processing required.     

Properties of superresonant systems of spherical scatterers

Systems of identical precisely spaced bubbles or similar monopole scatterers in water-e.g., inflated balloons or thin-walled shells-insonified at frequenciesomega_{SR}dose to their fundamental radial resonanceomega_{0}(bubble) frequency may themselves display resonance modes or superresonances (SR's) [1]. Ordinary single-bubble resonances magnify the local free-field pressure amplitudep_{1}by a factor(ka)^{-1},abeing the radius andkthe wavenumber in water: for air bubbles or balloons in water, this factor is of the order of 70. Under SR conditions each member of the system amplifies the local free-field amplitude by a further factor of order(ka)^{-1}. Depending upon geometry and other constraints, the pressure fieldP_{SR}on the surface and in the interior of each scatterer will then be in the range of10^{3}p_{1}to5 times 10^{3} p_{1}. This paper investigates the sensitivity of this phenomenon to small departures from the ideal model. In particular, it examines the effect of small differences in scatter positioning and volumes in the context of an SR system consisting of two bubbles/balloons close to the boundary of a thin elastic plate overlying a fluid half-space. It is found that, to observe the SR phenomenon, radii and positions should be controlled to within approximately 1/2 percent.P_{SR}is also sensitive to the angle of incidence of the plane wave train. For the simple system examined here, this sensitivity is considerable for either flexural wave trains or volume acoustic waves incident upon the bubble/ balloon pair (doublet). Practical uses of the phenomenon may range from the design of passive high-Qacoustical filter/amplifiers and acoustical lenses to improved source efficiencies.     

Ocean acoustic field matching: normal mode filtering andnon-Gaussian sources

Two qualitative results concerning statistical sonar signal processing and acoustic field matching are obtained. First, normal-mode field predictions are integrated with statistical signature analysis by constructing a boundary-value problem in the acoustic waveguide. From this construction it is found that the normal-mode filter is the unique acoustic preprocessor which does not confound deterministic waveguide correlation structure with stochastic source covariance structure. Second, the origin of deterministic, Gaussian, and non-Gaussian source signatures is investigated by associating physical parameters with the classical Lindeberg central limit conditions. From construction it is found that there are important objects that are not adequately represented either by infinitesimal points or by infinite surfaces. If receiver resolution is inadequate to resolve source complexity, these objects will exhibit a non-Gaussian acoustic signature via an entirely linear progression from internal excitation, to source radiation, through waveguide propagation, and finally to reception     

中国对虾幼体对运动物体的感觉能力及其与环境照度的关系

甲壳动物具有明显的视觉运动反应。这种对视野内移动目标的反应有两种形式,一种是眼柄的颤动反应(Optokinetic response),另一种为全身性的视觉运动反应(Optomotor response)。对于前一种反应形式国外已有许多研究,尤其对爬行亚目某些种类眼柄颤动反应特点和机制进行了大量实验。后一种反应形式的研究工作较少。有关游泳亚目视觉     

Dissipation observations of drag coefficients over the open ocean

Forty-three open-ocean observations of drag coefficients observed at Argus Island Tower near Bermuda by the dissipation technique resulted in constant drag coefficients for mean horizontal wind velocities between 7.8 and10.4 m .s^{1}in good agreement with the larger near-neutral data set of DeLeonibus and Simpson [1] and the neutral data of Large and Pond [2], both of whom observed10^{3}C_{10} = 1.2whereC_{10}is the drag coefficient at l0 m. Ratios of vertical-to-horizontal wind velocity spectral densities averaged over an inertial subrange of 0.8 to 1.6 Hz ranged from 0.7 to 1.07 in agreement with the Busch and Panofsky [3] result that isotropy is approached only when the observation height is much greater than the Nyquist wavelength.     

Matched-Filter Envelope Detector Deflection Performance for a Correlated Phase Channel

In this paper, the effect of channel phase coherence upon a matched filter envelope detector output is investigated for a pulsed radar or active sonar. A novel model for the correlated channel phases allows the explicit calculation of the loss in detection performance using the deflection criteria. The theoretical model yields good agreement with simulations when the phase correlation coefficients between the first and last pulses are between 0.1 and 1.0. It is shown that a 3-dB loss in performance, as compared to the optimum detector for perfect coherence, requires phase correlation between adjacent pulses ofrho_{i,i+1} = 0.91, 0.96, and 0.96 for 10, 20, and 30 pulses, respectively. On the other hand, the same performance is obtained with a noncoherent combiner of the matched filter pulse returns when correlation between adjacent pulses,rho_{i,i+1} = 0.8, 0.835, and 0.84 for 10, 20, and 30 pulses, respectively. Ifrho_{i,i+1}is smaller than these quantities, one is better off performing noncoherent detection.     

Shallow-Water Acoustic Tomography Performed From a Double-Beamforming Algorithm: Simulation Results

Recent shallow-water experiments in sea channels have been performed using two vertical coplanar densely sampled source and receive arrays. Applying a double-beamforming algorithm on the two arrays both on synthetic numerical simulations and on experimental data sets, we extract efficiently source and receive angles as well as travel times for a large number of acoustic rays that propagate and bounce in the shallow-water waveguide. We then investigate how well sound-speed variations in the waveguide are reconstructed using a ray time-delay tomography based on a Bayesian inversion formulation. We introduce both data and model covariance matrices and we discuss on the synthetic numerical example how to choose the a priori information on the sound-speed covariance matrix. We attribute the partial sound-speed reconstruction to the ray-based tomography and we suggest that finite-frequency effects should be considered as the vertical and horizontal size of the Fresnel zone significantly spreads in the waveguide. Finally, the contribution of a different set of ray angles for tomography goal is also presented.      

Wind-induced drift of objects at sea: The leeway field method

A method for conducting leeway field experiments to establish the drift properties of small objects (0.1-25 m) is described. The objective is to define a standardized and unambiguous procedure for condensing the drift properties down to a set of coefficients that may be incorporated into existing stochastic trajectory forecast models for drifting objects of concern to search and rescue operations and other activities involving vessels lost at sea such as containers with hazardous material.An operational definition of the slip or wind and wave-induced motion of a drifting object relative to the ambient current is proposed. This definition taken together with a strict adherence to a 10 m wind speed allows us to refer unambiguously to the leeway of a drifting object. We recommend that all objects if possible be studied using what we term the direct method, where the object’s leeway is studied directly using an attached current meter.We establish a minimum set of parameters that should be estimated for a drifting object for it to be included in the operational forecast models used for prediction of search areas for drifting objects.We divide drifting objects into four categories, depending on their size. For the smaller objects (less than 0.5 m), an indirect method of measuring the object’s motion relative to the ambient current must be used. For larger objects, direct measurement of the motion through the near-surface water masses is strongly recommended. Larger objects are categorized according to the ability to attach current meters and wind monitoring systems to them.The leeway field method proposed here is illustrated with results from field work where three objects were studied in their distress configuration; a 1:3.3 sized model of a 40-foot Shipping container, a World War II mine and a 220 l (55-gallon) oil drum.     

Inversion in shallow-water using the WKB modal condition

In this paper, we propose an efficient scheme for extracting the complex bottom reflection coefficient (phase and amplitude) in a shallow water waveguide by using the WKB modal condition. The input data are the measured modal wavenumber (k_m) and the modal attenuation (β_m). The main advantages of this scheme are: 1) it is efficient, because there are no replica calculations as in the conventional matched field processing (MFP) scheme, 2) the inverted error induced by the “noise” of the contaminated data can be estimated analytically, and 3) the impact of the environmental (SSP) mismatch can also be estimated analytically. Numerical simulations illustrate that the proposed scheme works well in different scenarios of shallow water waveguides     

Robust tracking of multiple objects in sector-scan sonar imagesequences using optical flow motion estimation

The fast update rate and good performance of new generation electronic sector scanning sonars is now allowing practicable use of temporal information for signal processing tasks such as object classification and motion estimation. Problems remain, however, as objects change appearance, merge, maneuver, move in and out of the field of view, and split due to poor segmentation. This paper presents an approach to the segmentation, two-dimensional motion estimation, and subsequent tracking of multiple objects in sequences of sector scan sonar images. Applications such as ROV obstacle avoidance, visual servoing, and underwater surveillance are relevant. Initially, static and moving objects are distinguished in the sonar image sequence using frequency-domain filtering. Optical flow calculations are then performed on moving objects with significant size to obtain magnitude and direction motion estimates. Matches of these motion estimates, and the future positions they predict, are then used as a basis for identifying corresponding objects in adjacent scans. To enhance robustness, a tracking tree is constructed storing multiple possible correspondences and cumulative confidence values obtained from successive compatibility measures. Deferred decision making is then employed to enable best estimates of object tracks to be updated as subsequent scans produce new information. The method is shown to work well, with good tracking performance when objects merge, split, and change shape. The optical flow is demonstrated to give position prediction errors of between 10 and 50 cm (1%-5% of scan range), with no violation of smoothness assumptions using sample rates between 4 and 1 frames/s     

Channel capacity in bits per joule

Underwater acoustic telemetry has a total input energy constraint, since the energy is stored in the transmitter's batteries. This work is primarily rephrasing the work in channel capacity in terms ofC_{J}bit/J as compared toCbit/s, to emphasize the energy efficiency and to deemphasize the speed of the telemetry. The energy channel capacityC_{J}, for any waveform channel with well-defined capacityCbit/s at signal powerSwatts, is defined asC/Sbit/J. It is shown that for coherent binary frequency shift keying (BFSK) and waveform channels, the supremum ofC_{J}overSis approached asSapproaches zero. For the Gaussian channel the best coding uses narrow bands with the highest S/N.     

Effects of sedimentary interfaces on fracture pattern,linkage, and cluster formation in peritidal carbonate rocks

The permeability of a reservoir is particularly dependent upon the proportion of its fractures that penetrate or are arrested at interfaces such as contacts and discontinuities. Here we report on fracture penetration and fracture arrest in Lower Cretaceous peritidal deposits exposed in the Pizzicoli Quarry, Gargano Promontory, southern Italy. We measured more than 2000 fractures, in the field and using LIDAR data, of which 564 fractures from the field and 518 from LIDAR studies are the focus of this paper. Fracture arrest/deflection and penetration depend much on the effects of peritidal cycle interfaces such as paleosol horizons, laminated carbonate mudstones, and stylonodular horizons. The laminated mudstones have the greatest effect; 63–99% of the fractures are deflected or arrested at such interfaces, whereas 63–90% are deflected/arrested at paleosols, and 20–35% at stylonodular horizons. In the mudstones, many fractures are arrested at thin, internal laminae, such that few penetrate the entire laminated layer, and fewer still the boundaries between the layers. Paleosol interfaces deflect/arrest more than 60% of all fractures. However, when small-offset fractures above and below paleosols are regarded as penetrating, they are evenly spaced (non-clustered), so that fracture-related fluid transport may occur across the entire paleosol. Stylonodular horizons deflect/arrest and split some fractures, but generally have little effect compared with the other types of interfaces. We present three main mechanisms for fracture deflection and/or arrest: (1) the fracture-induced tensile stress ahead of its tip, referred to as the Cook-Gordon debonding mechanism; (2) rotation of the principal stresses at and across the interface, resulting in the formation of stress barriers; and (3) large elastic mismatch (particularly as regards Young’s moduli) between layers across an interface. All these mechanisms are likely to have operated during fracture propagation and arrest in the carbonate rocks of the Pizzicoli Quarry.     

Liquid waveguide capillary cell for the spectrophotometric determination of nanomolar iodate concentrations in marine waters

A new spectrophotometric method based on a liquid waveguide capillary cell for an enhanced detection was developed to measure nanomolar iodate concentrations. This method has a detection limit and precision of 1–2 nmol/L, which is equivalent to 10% that of conventional methods, a recovery of 97.7%–104.0%, and a working range of 10–120 nmol/L. Water samples were collected from three estuaries and one coastal ocean for testing, and the proposed technique detected as low as 11 nmol/L and 18 nmol/L ...     

Effects of seabed properties on acoustic wave fields in a seismo-acoustic ocean waveguide

Acoustic wave fields in an ocean waveguide with a sediment layer having continuously varying density and sound speed overlying an elastic subbottom are considered in this analysis. The objective of this study is to investigate the effects of seabed acoustic properties, including the density and sound speed of the sediment layer and subbottom, on the characteristics of the wave fields. Examination of the reflection coefficient, wavenumber spectrum, and noise intensity of the sound field through numerical analysis has shown that the variation in the acoustic properties in the sediment layer is an important factor in determining the reflected or noise sound fields. In particular, the sediment thickness-to-wavelength ratio and the types of variation of acoustic properties inside the layer give rise to many characteristics that potentially allow for acoustic inversion of the seabed properties. With regard to the wave-field components in a shallow-water environment, the various types of waves existing in a seismo-acoustic waveguide have been illustrated. The results indicate that the effects of the sediment properties on the wavenumber spectrum are primarily on the continuous and evanescent regimes of the wave field. The noise intensity generated by distributive random monopoles at various depths, together with the effect of refractive sound-speed distribution in the water column, has been obtained and analyzed.     

联合星载光学和SAR影像的漳江口红树林与互花米草遥感监测

文章以福建省漳江口国家级红树林自然保护区为研究区, 提出了一种联合星载光学和合成孔径雷达(SAR)影像的红树林与互花米草分布的自动提取算法, 提取研究区内红树林和互花米草的空间分布。本研究选择2016年、2017年和2018年各一景低潮时的Sentinel-2 A光学影像数据, 获取植被和其他地物的光谱和纹理信息。算法首先基于归一化植被指数、增强型植被指数、地表水分指数以及数字高程模型的相关规则计算红树林和互花米草的潜在分布区; 通过随机森林分类算法区分红树林和互花米草, 2016年、2017年和2018年影像的分类总体精度和Kappa系数分别为98.53%和0.980、96.52%和0.952、98.71%和0.978, 分类效果良好; 使用当年所有Sentinel-1 A/B的SAR影像获得研究区的常年海水分布范围, 使用与海水交界的判据, 实现红树林、互花米草提取范围的优化。研究表明, 研究区2016年、2017年和2018年红树林总面积分别为56.85hm²、59.88hm²和58.61hm², 互花米草总面积分别为109.23hm²、124.00hm²和142.39hm², 与前人的研究成果在红树林和互花米草的空间分布和面积量级上有较好的一致性。     

Hydroacoustic Field of a Point Source in an Inhomogeneous Marine Medium Containing a Cylindrical Body Floating on the Surface

We study the problem of determination of the sound field of a point harmonic source in the coastal zone and the influence of a cylindrical body floating above the source on the sound field formed in the marine medium. A numerical-analytic method is proposed for the determination of the velocity potential. According to this method, the unknown coefficients in the general solution of the problem are determined from the corresponding infinite system of linear algebraic equations by the method of reduction. We present results of numerical calculations for a special case of a waveguide whose parameters are typical of the coastal part of the sea and perform the comparative analysis of the data obtained as a result of variation of the indicated parameters.     

三维弹性侧壁液舱内液体晃动波面的实验研究

通过物理模型实验,对弹性侧壁液舱和刚性液舱内液体晃荡问题进行了研究。由于流固耦合的影响,弹性侧壁液舱内液体晃荡的最低阶固有频率稍小于同尺寸的刚性液舱内液体晃荡的最低阶固有频率。液舱模型处于纵向简谐激励作用下,其中激励频率在最低阶固有频率附近。实验分析两种相对液深比h/L=0.167和h/L=0.333,在二阶模态的次共振和一阶模态的共振状况下,对弹性侧壁液舱与刚性液舱内不同测点的波面、振幅谱和晃动波高进行对比分析。结果表明:在浅液深(h/L=0.167)一阶共振下,流固耦合对波面形态的影响比较明显,弹性侧壁液舱内测点晃动波高明显大于刚性液舱内对应测点波高;而在一般液深(h/L=0.333)一阶共振下,水弹性效应减弱,弹性侧壁液舱与刚性液舱内对应测点处波高差异较小。