The width and height of FEM were 60 m and 6 m, respectively, and two fluid mediums with the size of 59 m 3.6 m were surrounded by a perfectly matched layer (PML); the two fluid mediums consisted of water and sandy sediment with the heights of 1.3 and 2.3 m and had open boundaries. The sediment was assumed to be sand, but shear properties were not considered. Densities () and sound speeds () of water (upper) and sandy sediment (upper and lower fluid mediums, respectively) are 1000 kg/m3, 1900 kg/m3, 1500 m/s, and 1650 m/s, respectively. The aluminum shell target, whose density and longitudinal and transverse speeds are 2700 kg/m3, 6148.9 m/s, and 3097.3 m/s, respectively, was fully buried in the sediment and its inside was filled with air. The target center was at the depth of 0.35 m from the water–sediment interface and the target radius and thickness were 0.25 m and 0.025 m, respectively. The incident plane wave from the water penetrated the sediment with a subcritical angle (20° herein), and scattering occurred owing to the buried target illuminated by the penetrating wave. The scattered signals corresponding to the frequency range of 100 Hz–15 kHz was simulated using the FEM with a mesh size of the one-twelfth corresponding acoustic wavelengths. Triangular and quadrilateral meshes were used for the two fluid mediums and PML, respectively. All parameters including the model and mesh sizes were determined through a convergence test; the FEM results were observed according to various model and mesh sizes.
To simulate the scattered signal from the realistic target, an efficient model for the buried target scattering was developed using the ray approach. It is noteworthy that the ray approach has been used to model scattered signals from a target in a free field or on the sea bottom by combining wave propagation from the source to the target and from the target to the receiver with a target scattering pattern along scattering angles (i.e., the target scattering amplitude or cross-section) [
12,
13]. However, to model long-range buried target scattering, the ray approach is restricted to low-frequency regions, and the target scattering at high frequencies is simulated by considering the bottom roughness [
18,
19,
20]. In this study, the ray approach was extended to model the target scattering at low- and high-frequency regions (100 Hz to 15 kHz) under a perfectly flat bottom by inspecting the wave propagation of incident waves below the critical angle in the sediment and the corresponding target scattering cross-section, thereby allowing the ray-based unified scheme to model scattered signals from buried targets over a long range.
The wave propagation of subcritical incident waves exhibits the following features: (1) The direction of sound is parallel to the boundary between water and the sediment; (2) the wave propagation speeds in the sediment depend on incident angles below the critical angle; (3) the sound becomes evanescent along the depth, and the decay rate depends on the frequency and incident angle.
Section 2.1 provides an explanation of the features in detail.
Section 2.2 presents an analysis of the target scattering cross-section of the target illuminated by the evanescent wave having the same decaying rate and sound speed as the subcritical penetrating wave (i.e., subcritical scattering from a buried target), based on a comparison with that of a standard plane wave.
Section 2.3 presents a unified ray approach for modeling the scattered signals of the long-range buried target based on conclusions from previous subsections.
2.1. Incident Waves below Critical Angles
A long-range incident wave from the source can be approximated using a plane wave with a low grazing angle that penetrates the sea bottom at a subcritical angle in an ocean environment involving a hard bottom. Because the wave traveling along the boundary must coincide with the horizontal components of waves in water and the sediment (i.e.,
, where
and
are the wave vectors of the first (water) and second (sediment) mediums, respectively, and the subscript
x indicates the horizontal components), the penetrating wave in the sediment propagates in the direction parallel to the boundary regardless of the incident angles below the critical angle
(
and
are sound speeds in water and the sediment, respectively) and becomes evanescent along the depth. The decay rate
is defined as follows:
where
is the pressure at depth
. Subsequently, the decay rate is expressed analytically as
for the penetrating wave below the critical angle, where
is the wavenumber of the sediment (i.e.,
), and
(
) is the incident angle of the wave into the sediment [
21].
The penetrating wave, which propagates horizontally in the sediment, has a speed of
owing to the continuity of the wave on the boundary. Hence, the penetrating wave propagates at a lower speed than that of the sediment, and its propagation speed depends on the incident angle. In a relevant study [
14], when modeling the penetrating wave, it was observed that the wave propagation speeds in the sediment deviated from the measured sediment sound speed for incident waves below critical angles. Subsequently, they were fitted to enhance the agreement between the model results and measured acoustic data. The fitted sound speeds in the sediment were 1650, 1685, and 1720 m/s for grazing angles of 22.0°, 24.8°, and 27.2°, respectively, and the sound speeds of water and the sediment were 1530 and 1720 m/s, respectively. The fitted sound speeds are similar to the horizontally propagating sound speed of
, which are 1650.2, 1685.4, and 1720.2 m/s at 22.0°, 24.8°, and 27.2°, respectively. Therefore, they can be explained simply based on the continuity of the wave on the boundary.
The FEM was used to confirm the features of the penetrating wave, which are shown in
Figure 2. The numerical environment was the same as that in
Figure 1 except for the absence of the target. The penetrating waves were simulated for incident waves with two different incident angles (10° and 20°) at 3 kHz. As expected, the wavefronts in the sediment were vertical, and the sounds propagated horizontally. Each penetrating wave propagated at a distance of the corresponding single wavelength in the sediment during one period, and the propagation speed is calculated by dividing the wavelength in the FEM result with one period. Here, the penetrating wave along the horizontal axis at the fixed depth of 0.02 m was used to measure the wavelength required for calculating the propagation speed. The penetrating wave speeds from the FEM were 1525.6 and 1597.4 m/s for incident angles of 10° and 20° at 3 kHz, respectively. They are similar to those computed using
, which are 1523.1 and 1596.3 m/s at 10° and 20°, respectively.
Meanwhile, the plane wave amplitude diminished as the depth increased. To calculate the decaying rate in the FEM, the penetrating wave along the vertical axis at the fixed horizontal range of 0 m was used. Specifically, pressures of the penetrating wave at different depths from the FEM were used, based on Equation (1). The decaying rates using the FEM results were 4.76 and 2.98 for the incident angles of 10° and 20°, respectively, and they are similar to those obtained from the analytical expression of (4.76 and 2.99 at 10° and 20°, respectively).
2.2. Target Strength Based on Subcritical Incident Waves
Sounds were reradiated from a target with different amplitudes along the scattering angles (
) after the target was illuminated by an incoming wave with angle (
), which is the angle incident to the target. Here, the incident and scattering angles were measured from the horizontal axis passing through the target center. The target scattering cross-section
(or its logarithmic measure, TS) was calculated using scattered waves at positions encompassing the target, as follows [
20]:
where
and
are the amplitudes of the incoming and scattered waves, respectively;
is the distance between the receiver and target. Because the 2-D target scattering was simulated for efficiency, the target scattering cross-section was independent of the azimuthal angle, and the square of the distance in the 3-D target scattering was replaced with the distance.
Generally, when calculating the target scattering cross-section, an incident wave with a planar wavefront, whose amplitude is constant (referred to as standard plane wave herein), illuminates the free-field target with an arbitrary incident angle, and the physical properties (such as the density and sound speed) of the medium surrounding the target is used.
Meanwhile, the long-range buried target is illuminated by the penetrating wave below the critical angle, which induces a different target scattering cross-section from the general one. The penetrating wave (the incoming wave for the buried target) propagates horizontally, and the target scattering cross-section becomes independent of the incident angle of the incoming wave to the target; here, . Furthermore, the incoming wave propagates at a lower speed of than the original sound speed of the sediment surrounding the target, and its amplitude decays along the depth.
When calculating the target scattering cross-section for the long-range buried target, the features of the incoming wave must be considered; hence, the target scattering cross-section is a function of
instead of
. Accordingly, the general definition of the target scattering cross-section is replaced with the following equation:
Here,
is the average amplitude of the incoming wave over the target.
where
is the amplitude of the incoming wave at the depth
,
is the depth of the target center measured from the water–sediment interface, and
is the radius of the target. Because the target scattering cross-section is the ratio between incoming wave intensity and scattered wave power, it is independent on
and
.
is the scattered wave amplitude at a distance of
from the target center when the incoming wave into the free-field target is the evanescent wave having the same decaying rate and propagation speed as the penetrating wave at subcritical angle
. Whereas
can be expressed analytically under the 2-D cylindrical shell illuminated by the standard plane wave, the FEM is essential to computing
for the incoming evanescent wave.
When calculating a scattered signal by target, a scattered amplitude
is used and has a relationship with the scattering cross section as
. Specifically, it is defined for the evanescent wave as follows:
is a scattered signal and its amplitude is . is the wavenumber of medium surrounding the target and is equivalent to the sediment wavenumber for the buried target. and compensate spreading loss and phase change during propagation from the target to receiver, respectively.
In this study, we calculated TSs at low (3 kHz) and high (9 kHz) frequencies for two different incoming waves to the target using the FEM; one for the evanescent wave having the same decaying rate and propagation speed as the penetrating wave at the incident angle of 20°, and the other for the standard plane wave in the front of the target. In the TS calculation, the free-field target was surrounded by the same sediment as that in
Figure 1. The mesh size for each TS calculation was the one-twelfth acoustic wavelength at the corresponding frequency. By comparing the TSs, the TS features of the buried target illuminated by the subcritical penetrating wave were investigated.
Figure 3 shows the TSs, where Equations (2) and (3) were used for the standard plane and evanescent waves, respectively.
(or
) for calculating the TS was computed along
from 0° and 360° at the fixed distance (
m). Whereas the standard waves with an incident angle of zero (
) rendered the TSs symmetric, those for the evanescent waves became asymmetric owing to the depth-dependent amplitude. In particular, larger amplitudes of evanescent waves at the upper section of the target yielded larger TSs between 0° and 180°, which manifested primarily at the higher frequency owing to the larger decay rate. The greater TSs in the upper direction by the evanescent waves facilitated the detection of the buried target. However, as the burial depth increased, the scattered signal was weakened by decreasing
.
Meanwhile, the TSs for the evanescent waves exhibited more complicated patterns owing to the higher wavenumber induced by the lower sound speed, compared with those for the evanescent and standard plane waves propagating with the bottom sound speed. The features of the evanescent waves make the corresponding TSs distinguished from those of the standard waves, and, thus, scattering patterns of buried target illuminated by subcritical penetrating waves are different from the general ones.
2.3. Ray Approach for Modeling Scattered Signals from Long-Range Buried Target
The ray approach has been used to model scattered signals from a target owing to its efficiency. In the ray approach, signals scattered from the target are simulated by combining wave propagations from the source to the target and from the target to the receiver with the target scattering amplitude. Hence, in ray-based modeling, the scattered signals can be analyzed based on their characteristics, such as propagation and scattering, thereby enabling their physical interpretation.
Conventionally, the ray approach has been applied to simulate scattered signals from the target above the sea bottom (proud target) by considering the interference of multiple ray paths from the sea bottom [
12,
13]. In this study, the ray approach was modified to model scattered signals from long-range buried targets at low-to-high frequencies. It is noteworthy that the bottom roughness was considered for modeling the target scattering at high frequencies, which enabled an incident wave to penetrate the bottom and illuminate the target [
18,
19,
20,
22]. However, in the current study, with the assumption of a flat bottom, the buried target scattering at high frequency was simulated by considering the features of the penetrating waves below the critical angle and the corresponding scattering amplitude.
Simulating subcritical scattered signals from a buried target is composed of two steps in the ray approach: (1) The subcritical scattering from the buried target and (2) the sound propagating from source to target and from target to receiver. In the first step, the free-field target is surrounded by sediment and is illuminated by the evanescent wave having the same decaying rate and propagation speed as the penetrating waves at a subcritical angle. Subsequently, a sound is reradiated from the target with different amplitudes according to scattering angles (target scattering amplitude or cross-section) and scattered signals arrive at the receiver via single or multiple paths depending on a sound propagating environment.
In this study, for convenience, a receiver for capturing the scattered signal was installed in the bottom, as shown in
Figure 4a. The dominant scattering signals propagated via direct and bottom reflected paths, where the cylindrical spreading loss in the bottom indicated a constant sound speed. When using the ray approach, the signal scattered to the receiver in the bottom (
) is denoted as follows:
in
based on Equation (4) is
where
and
are the incident wave amplitude and the transmission coefficient from the water to the bottom, respectively.
is the reflection coefficient from the bottom to water. The scattering amplitude
is calculated using scattered signals from the free-field target illuminated by the evanescent wave having the same decaying rate and propagation speed as the subcritical penetrating wave, based on Equation (5).
and
are the scattering angles and propagation distances for the direct and bottom reflected paths, respectively (
Figure 4a).
is expressed differently for forward and backward scatterings:
for forward scattering and
for backward scattering.
A scattered signal was simulated using the ray approach for a receiver at the same depth of the target center (0.35 m) under the numerical environment of
Figure 1. The distance between receiver and target is set to 2 m. Subsequently, the signal was compared with that obtained using the FEM with respect to amplitude (sound pressure level) and phase, as shown in
Figure 4b. The incident wave from the water, whose amplitude was 1 Pa, was incident into the bottom with the angle of 20° (i.e.,
Pa and
°). While overall patterns such as peak positions were similar between the scattered signals, discrepancies were observed owing to approximation in the ray approach. Whereas the FEM fully solved the Helmholtz equation with the given boundary conditions, the ray approach derived the scattered signal approximately by combining wave propagation with target scattering pattern (i.e., target scattering amplitude). Particularly, in the ray approach, the target was treated as a directional point source by non-uniform reradiation from the scattering amplitude, which accounted for the difference of the ray result from the FEM result. The discrepancies were more apparent when the receiver was located closer to the target owing to the complicated near-field target scattering [
5,
23]. However, as shown in
Figure 4b, the ray approach can provide reliable scattering results efficiently in the frequency range of 100 Hz–15 kHz in a unified manner, particularly, for the far-field receiver from the target.
While peak amplitudes were less than those yielded by the shell target illuminated by the standard plane wave owing to the larger decaying rate of penetrating wave at higher frequencies, the strong peaks were still yielded by the MFE from the buried shell target.
So far, bottom attenuation was neglected to focus on the effects of penetrating waves on subcritical target scattering. To inspect the MFE detectability in a more realistic ocean environment, the scattered signals were simulated using the FEM and ray approach while considering the bottom attenuation
dB/wavelength, and they were in good agreement (
Figure 4c). In the ray approach, the scattering amplitude was calculated using the free-field target surrounded by the bottom with attenuation and the additional terms of
,
, and
were added to Equations (5), (7) and (8) to account for loss by the bottom attenuation (
, where
is a wavelength) [
24]. The MFE peaks were observed although the high-frequency peaks diminished owing to the significant bottom attenuation at high frequency.
Finally, the frequency-domain solutions from the FEM and ray approach in
Figure 4b,c were converted to the corresponding time signals using Fourier transform as shown in
Figure 5; 1–14 kHz linear frequency modulated pulse signal was used as a source waveform. While fictitious signals from numerical errors during the conversion were observed in the first and later parts, the time signals from the FEM and ray approach were in good agreement (e.g., the same first arrival time near 1.2 ms).