Acoustical Measurement and Biot Model for Coral Reef Detection and Quantification

Coral reefs are coastal resources and very useful for marine ecosystems. Nowadays, the existence of coral reefs is seriously threatened due to the activities of blast fishing, coral mining, marine sedimentation, pollution, and global climate change. To determine the existence of coral reefs, it is necessary to study them comprehensively. One method to study a coral reef by using a propagation of sound waves is proposed. In this research, the measurement of reflection coefficient, transmission coefficient, acoustic backscattering, hardness, and roughness of coral reefs has been conducted using acoustic instruments and numerical modeling using Biot theory.The results showed that the quantification of the acoustic backscatter can classify the type of coral reef.


Introduction
As an archipelago country, Indonesia has many types of coral reefs.Coral reefs consist of geological structures of calcium carbonate built over time by tiny living organisms.These structures are found in the marine waters which contain nutrients.Coral reefs had been used for marine tourism, shoreline protection, and fishing ground to the local communities.Coral reefs, which are called rainforests of the sea, are the most biologically diverse ecosystems.They support variety of marine live forms, provide pharmaceutical materials, and generate income from tourism [1].
Coral reefs are valuable marine resources; however, they are in the vulnerable condition [2].In spite of their advantages and in spite of the fact that humans continue to benefit from coral reefs, the marine ecosystem has been subjected to serious anthropogenic threats.Nonsustainable activities and developments have inevitably caused degradation to the coral reef ecosystem.Blast fishing activities, mining of coral reef for building materials, and other threats like sedimentation and pollution in coastal areas are diminishing the ecosystem.
The greatest threats to coral reefs are rising seawater temperature and ocean acidification.With the recent global ocean warming, coral bleaching and mortality have become more frequent.Human induced disturbances, storm, and natural predators also affected the damage of coral reef [2].
However, the lack of accurate and comprehensive coral reef database has always been one of the limiting factors in their conservation and management efforts.The marine scientists and the coastal resource managers need to know the distribution and the status of these coral reefs.

State of the Art and Scientific Contribution.
Detection and mapping of coral reef ecosystem should be at a scale that is adequate to detail information on the distribution of the major coral reef substrates.Unfortunately, the current coral reef research methodologies are still unable to fulfill this need.This is mainly due to the high cost involved and inherent shortcomings of the remote sensing technologies.The conventional coral reef research methods use a diver by applying the line transect, quadrate plot, belt transect, and manta tow survey [3].All of these methods are used to derive both quantitative and qualitative data of coral reefs on a small scale surveyed area.In order to produce contiguous, broad scale resource map, marine scientists have explored the underwater videography and LIDAR survey using airborne and satellite remote sensing [4].These methods yield different degrees of resolution and accuracy depending on the sensor type, cost, time, ground truthing, and postprocessing.LIDAR airborne and satellite remote sensing derived images.Several researchers mapped the coral reef using satellite remote sensing [5].This remote sensing technology was originally designed for terrestrial application.A number of limitations become obvious, when remote sensing is applied to map underwater marine resources.The electromagnetic wave used in satellite remote sensing was more absorbed in ocean environment.
To overcome this problem, the application of underwater acoustic technologies is used to detect and quantify coral reef.Acoustic technologies offer a promising method for underwater remote sensing techniques [6][7][8].The sound waves can propagate deeper compared with electromagnetic wave.This method allows for large area coverage in greater details within a shorter time at a relatively low cost.Acoustic detection using side scan sonar and multibeam echosounder has been successful for water column and seabed imaging [9,10].Bottom backscattering strength was computed for acoustic seabed classification [11,12].
The objective of this research is to develop a methodology for effective ecosystem monitoring and mapping of tropical marine coral reefs using underwater acoustic technology in Seribu Island waters, Indonesia.

Research Method
The research was conducted at Ocean Acoustics and Instrumentation Laboratory, Department of Marine Science and Technology, Faculty of Fisheries and Marine Sciences, Bogor Agricultural University, and Seribu Island waters, Jakarta, Indonesia (Figure 1).

Flow Diagram of the Research.
The flow diagram of this research in the laboratory is shown in Figure 2. Research output is to develop and apply the methodology of high resolution underwater acoustic technology, calibrating the instrument and measuring acoustic backscattering from coral reef.

Hydroacoustic Data Acquisition, Processing, and Analysis.
Based on underwater acoustic theory, the ability of this method is applied to discriminate coral reef types, to map coral substrates, and to measure acoustic reflectivity of coral reef at Seribu Island waters.The result encourage where the acoustic instrument was able to pick up distinct echo when ensonifying over coral reef type.All ocean surveys were conducted on the research vessel.The ship is outfitted with a suite of underwater acoustics instrument.Matlab software was applied and used in postprocessing and visualization to overlay pings and backscatter returns in both space and time.
Raw data was corrected for ambient noise and interference artifacts.A subset of pings were identified in the acoustic data and concurrent pings are selected in the acoustic data records.Coral reef detection algorithms were developed and used in Matlab to delineate bottom fish and seabed backscatter.
Consecutive echoes representing individual coral are similarly delineated in the hydroacoustic image.In the case of dense coral, tracks are selected on the outer limits of the coral.The average backscattering strength is calculated for the series of echoes.Figure 3 shows a simplified block diagram for coral reef backscattering measurement by quantitative echo sounder.
The backscattered pressure signal from the bottom received by the transducer (Figure 3),  RB [12], is where   is the source pressure level.The bottom echo signal is amplified by the amplifier to give where  RB is the echo amplitude at the preamplifier output,  is the receiving sensitivity of the transducer, and   is the preamplifier gain.The echo amplitude,  RB , is shown from ( 1) and (2) as where  TR =     is the transmitting receiving coefficient.The time varied gain (TVG) amplifier output of  RB corrected for absorption and spreading losses,  TB , is The raw SV value of the bottom echo,  VB , is where   is the multiple echo coefficient and is given by where The above Ψ is the equivalent beam angle of the volume scattering.Substitution of  TB and   into (6) yields This equation is called ring surface scattering model.

Geoacoustic Biot Model.
Coral reefs are porous, with the possibility that the fluid and granular phases will vibrate differently in response to acoustic excitation.Poroelastic theory or Biot theory treats both porosity and elasticity [13,14].
The incident and reflected waves in the seawater and coral will have displacement potential (Figure 4): where   = /  . is the angular frequency, and   is the sound velocity in the seawater.Biot developed a pair of coupled differential equations to describe acoustic wave propagation in an isotropic saturated porous medium with permeability , pore fluid viscosity , and bulk density   : where  is the incremental volume of fluid which enters or leaves the frame and  is the volumetric strain of the frame.Stoll developed the theory for Biot moduli in terms of measurable sediment properties: where   is the grain bulk modulus,   is the modulus of the pore fluid, and  is porosity.The frame bulk modulus was calculated by and the frame shear modulus was calculated by where   and   are the bulk and shear logarithmic decrement.
The boundary conditions are required at a seawater-coral reef interface.
(3) For equilibrium of tangential traction Four linear complex equations can be obtained for acoustic reflection measurement: where the components of {} and {} are given by the physical parameters of seawater and coral.  is the complex amplitude of the reflected wave.The reflection coefficient can be computed from this equation.Physical parameters for coral reef models were shown in Table 1.

Results and Analysis
In order to obtain accurate and precise quantitative measurements of coral reef resources, proper calibration of an acoustic system is necessary when using acoustic techniques such as echo detection and echo integration.Figure 5 shows the acoustic calibration using sphere ball method.The blue color shows the sphere detected on −42.5 dB and this value agreed with the theoretical value obtained by the manufacturer.The result of calibration was shown in Table 2. Acoustic data were collected using four frequencies (10, 50, 100, and 200 kHz) and pulse transmission was simultaneous at all frequencies.All transducers were suspended from a small boat and data were collected while drifting slowly over 15-40 m deep water.During the 6 h of data collection, the boat drifted approximately 3 km.Figure 6 shows the acoustic image of soft coral and hard coral.The amplitude intensity for soft coral ranged from −35.0 to −30.0 dB and that for hard coral ranged from −28.0 to −10.0 dB.
The hardness and roughness of coral reefs were shown in Figure 7.The hardness of soft coral ranged from −33.0 dB to −15.0 dB (Figure 7(a)) and the hardness of hard coral was about −22.0 dB to −15.0 dB.The roughness of soft coral  ranged from −50.0 to −20.0 dB and that of hard coral ranged from −65.0 to −22.0 dB.The harder coral had more energy due to higher acoustic impedance compared to soft coral.However, a very rough and hard surface can scatter much energy acoustically [15,16].
The parameters in Table 1 were used to examine the behavior of the density as a function of frequency.Figure 8 showed a plot of the real parts of the sound speed as a function of frequency from 0.01 Hz to 1000 kHz.At low frequencies, the pore size parameter is smaller than the viscous skin depth and at high frequencies the skin depth is smaller than the pore size [17].
Figures 9 and 10 show the density and attenuation as a function of frequency for Biot model.The bulk and shear moduli of the frame obviously increase the sound speed slightly.The effect of the larger difference in sound speed will be seen in the reflection, transmission, and backscattering strength.
The attenuation due to water viscosity in Biot model varies as frequency squared up to around 1 kHz.The attenuation due to frame moduli varies linearly with frequency.Numerical computation of Biot model confirms that the sound speed, density, and attenuation were frequency dependent.
Figure 11 shows the reflection coefficient as a function of grazing angle for four different frequencies (10, 50, 100, and 200 kHz) for the measurement and Biot model.The measurement of reflection coefficient was conducted by using the average intensity of the coral reef-water interface echo.The biggest difference is near the critical angle (near 20 ∘ to 30 ∘ ) and this is due to the difference in sound speed between coral reef and seawater.The reflection coefficient is strongly correlated with sediment porosity and bulk density.Based on this research, the porosity of the coral reef was estimated using the reflection coefficient.
The general trend for the transmission coefficient is shown in Figure 12.For the Biot model, the procedure used to calculate the transmitting coefficient using the scalar potential and fast wave.This figure means more energy is transmitted to the coral reef than reflected.
Figure 13 showed the backscattering strength as a function of grazing angle.Backscattering is dependent on many confounding parameters such as the coral composition, surface roughness, slope, and spatial heterogeneity of coral.Backscattering versus grazing angle can be exploited to detect differences in coral properties and to solve geoacoustic model of coral properties.Frequency response may indicate coral reef type.The objective of our study was to quantify backscattering for coral versus both frequency and grazing angle and to use parameters of Biot model to identify the type of coral reef.Backscattering strength of hard coral was higher than that of soft coral.These results were agreed on by previous researchers [8,18].

Conclusion
Underwater acoustic remote sensing measurement and Biot model had been used to detect and quantify the coral reefs.The hardness of hard coral was higher than that of soft coral.The roughness of coral contributed to reflectivity and backscattering energy.Numerical computation of Biot model confirms that the sound speed, density, and attenuation are frequency dependent.Acoustic measurement agreed with Biot model in reflection and backscattering value of coral reef.This remote acoustic measurement and Biot model may be useful for estimating coral reef properties when core data are not available.

Figure 2 :Figure 3 :
Figure 2: Flow diagram of underwater acoustic instrument for measuring coral reef.

Figure 12 :
Figure 12: Transmission coefficient using Biot's model as a function of grazing angle at four different frequencies for soft coral (blue line) and hard coral (black line).
Figure 4: Incident, reflected, and transmitted wave in seawater and coral reef.

Table 1 :
Physical parameters for coral reef models.

Table 2 :
Configuration of the four acoustic frequencies.