Parameter-sparse modification of Fourier methods to analyse the shape of closed contours with application to otolith outlines

-Elliptical Fourier descriptors (EFDs) have been used extensively in shape analysis of closed contours and have a range of marine applications, such as automatic identification of fish species and discrimination between fish stocks based on EFDs of otolith contours. A recent method (the ‘MIRR’ method) transforms the two-dimensional contour to a one-dimensional function by mirroring (reflecting) the lower half of the contour around a vertical axis at the right end of the contour. MIRR then applies the fast Fourier transform (FFT) to the vertical contour points corresponding to equidistant coordinate values along the horizontal axis. MIRR has the advantage of reducing the number of Fourier coefficients to two coefficients per frequency component compared with four EFDs. However, both Fourier methods require several frequency components to reproduce a pure ellipse properly. This paper shows how the methods can be easily modified so that a virtually perfect reproduction of a pure ellipse is obtained with only one frequency component. In addition, real otolith examples for cod (Gadus morhua) and Greenland halibut (Reinhardtius hippoglossoides) are used to demonstrate that the modified methods give better approximations to the large-scale shape of the original contour with fewer coefficients than the traditional Fourier methods, with negligible additional computing time

The Making of the NEAM Tsunami Hazard Model 2018 (NEAMTHM18)

The NEAM Tsunami Hazard Model 2018 (NEAMTHM18) is a probabilistic hazard model for tsunamis generated by earthquakes. It covers the coastlines of the North-eastern Atlantic, the Mediterranean, and connected seas (NEAM). NEAMTHM18 was designed as a three-phase project. The first two phases were dedicated to the model development and hazard calculations, following a formalized decision-making process based on a multiple-expert protocol. The third phase was dedicated to documentation and dissemination. The hazard assessment workflow was structured in Steps and Levels. There are four Steps: Step-1) probabilistic earthquake model; Step-2) tsunami generation and modeling in deep water; Step-3) shoaling and inundation; Step-4) hazard aggregation and uncertainty quantification. Each Step includes a different number of Levels. Level-0 always describes the input data; the other Levels describe the intermediate results needed to proceed from one Step to another. Alternative datasets and models were considered in the implementation. The epistemic hazard uncertainty was quantified through an ensemble modeling technique accounting for alternative models’ weights and yielding a distribution of hazard curves represented by the mean and various percentiles. Hazard curves were calculated at 2,343 Points of Interest (POI) distributed at an average spacing of ∼20 km. Precalculated probability maps for five maximum inundation heights (MIH) and hazard intensity maps for five average return periods (ARP) were produced from hazard curves. In the entire NEAM Region, MIHs of several meters are rare but not impossible. Considering a 2% probability of exceedance in 50 years (ARP≈2,475 years), the POIs with MIH >5 m are fewer than 1% and are all in the Mediterranean on Libya, Egypt, Cyprus, and Greece coasts. In the North-East Atlantic, POIs with MIH >3 m are on the coasts of Mauritania and Gulf of Cadiz. Overall, 30% of the POIs have MIH >1 m. NEAMTHM18 results and documentation are available through the TSUMAPS-NEAM project website (http://www.tsumaps-neam.eu/), featuring an interactive web mapper. Although the NEAMTHM18 cannot substitute in-depth analyses at local scales, it represents the first action to start local and more detailed hazard and risk assessments and contributes to designing evacuation maps for tsunami early warning.publishedVersio

Greenland halibut observed by video in front of survey trawl: behaviour, escapement, and spatial pattern

Abstract Video recordings of Greenland halibut (Reinhardtius hippoglossoides) were made at eight trawl stations in Svalbard waters in August 2002. The recordings were made down to 600 m depth using artificial light. A method for calculating actual fish length from the video image was established and the recordings were analysed with respect to length-dependent behaviour, escapement and spatial pattern. All Greenland halibut observed were either lying on the bottom or swimming in a horizontal position close to the bottom, and there was no tendency to schooling. Individual fish reacted in an ordered way to the approaching trawl and were herded along the ends of the ground-gear. Escapement under the ground-gear was higher for smaller fish, while some larger individuals were apparently able to escape the trawl ahead of the observed region.

The BIG’95 submarine landslide-generated tsunami: A numerical simulation

This article presents a reasonable present-day, sea-level highstand numerical simulation and scenario for a potential tsunami generated by a landslide with the characteristics of the BIG’95 debris flow, which occurred on the Ebro margin in the western Mediterranean Sea in prehistoric times (11,500 cal yr BP). The submarine landslide deposit covers an area of 2200 km2 of the slope and base of slope (200–1800-m water depth), involving a volume of 26 km3. A leapfrog finite difference model, COMCOT (Cornell multigrid coupled tsunami model), is used to simulate the propagation of the debris-flow-generated tsunami and its associated impact on the nearby Balearic Islands and Iberian coastlines. As a requisite of the model, reconstruction of the bathymetry before the landslide occurrence and seafloor variation during landsliding have been developed based on the conceptual and numerical model of Lastras et al. (2005). We have also taken into account all available multibeam bathymetry of the area and high-resolution seismic profiles of the debris flow deposit. The results of the numerical simulation are displayed using plots of snapshots at consecutive times, marigrams of synthetic stations, maximum amplitude plots, and spectral analyses. The obtained outputs show that the nearest shoreline, the Iberian coast, would not be the first one hit by the tsunami. The eastward, outgoing wave would arrive at Eivissa Island 18 min after the triggering of the slide and at Mallorca Island 9 min later, whereas the westward-spreading wave would hit the Iberian Peninsula 54 min after the slide was triggered. This noticeable delay in the arrival times at the peninsula is produced by the asymmetric bathymetry of the Catalano-Balearic Sea and the shoaling effect due to the presence of the wide Ebro continental shelf, which in addition significantly amplifies the tsunami wave (>9 m). The wave amplitudes attain 8 m in Eivissa, and waves up to 3 m high would arrive to Palma Bay. Resonance effects produced in the narrow Santa Ponça Bay in Mallorca Island could produce waves up to 9 m high. A similar event occurring today would have catastrophic consequences, especially in summer when human use of these tourist coasts increases significantly