Assessing the potential of autonomous submarine gliders for ecosystem monitoring across multiple trophic levels (plankton to cetaceans) and pollutants in shallow shelf seas

A combination of scientific, economic, technological and policy drivers is behind a recent upsurge in the use of marine autonomous systems (and accompanying miniaturized sensors) for environmental mapping and monitoring. Increased spatial–temporal resolution and coverage of data, at reduced cost, is particularly vital for effective spatial management of highly dynamic and heterogeneous shelf environments. This proof-of-concept study involves integration of a novel combination of sensors onto buoyancy-driven submarine gliders, in order to assess their suitability for ecosystem monitoring in shelf waters at a variety of trophic levels. Two shallow-water Slocum gliders were equipped with CTD and fluorometer to measure physical properties and chlorophyll, respectively. One glider was also equipped with a single-frequency echosounder to collect information on zooplankton and fish distribution. The other glider carried a Passive Acoustic Monitoring system to detect and record cetacean vocalizations, and a passive sampler to detect chemical contaminants in the water column. The two gliders were deployed together off southwest UK in autumn 2013, and targeted a known tidal-mixing front west of the Isles of Scilly. The gliders’ mission took about 40 days, with each glider travelling distances of >1000 km and undertaking >2500 dives to depths of up to 100 m. Controlling glider flight and alignment of the two glider trajectories proved to be particularly challenging due to strong tidal flows. However, the gliders continued to collect data in poor weather when an accompanying research vessel was unable to operate. In addition, all glider sensors generated useful data, with particularly interesting initial results relating to subsurface chlorophyll maxima and numerous fish/cetacean detections within the water column. The broader implications of this study for marine ecosystem monitoring with submarine gliders are discussed

A two-step approach to model precipitation extremes in California based on max-stable and marginal point processes

In modeling spatial extremes, the dependence structure is classically inferred by assuming that block maxima derive from max-stable processes. Weather stations provide daily records rather than just block maxima. The point process approach for univariate extreme value analysis, which uses more historical data and is preferred by some practitioners, does not adapt easily to the spatial setting. We propose a two-step approach with a composite likelihood that utilizes site-wise daily records in addition to block maxima. The procedure separates the estimation of marginal parameters and dependence parameters into two steps. The first step estimates the marginal parameters with an independence likelihood from the point process approach using daily records. Given the marginal parameter estimates, the second step estimates the dependence parameters with a pairwise likelihood using block maxima. In a simulation study, the two-step approach was found to be more efficient than the pairwise likelihood approach using only block maxima. The method was applied to study the effect of El Ni\~{n}o-Southern Oscillation on extreme precipitation in California with maximum daily winter precipitation from 35 sites over 55 years. Using site-specific generalized extreme value models, the two-step approach led to more sites detected with the El Ni\~{n}o effect, narrower confidence intervals for return levels and tighter confidence regions for risk measures of jointly defined events.Comment: Published at http://dx.doi.org/10.1214/14-AOAS804 in the Annals of Applied Statistics (http://www.imstat.org/aoas/) by the Institute of Mathematical Statistics (http://www.imstat.org

Max-infinitely divisible models and inference for spatial extremes

For many environmental processes, recent studies have shown that the dependence strength is decreasing when quantile levels increase. This implies that the popular max-stable models are inadequate to capture the rate of joint tail decay, and to estimate joint extremal probabilities beyond observed levels. We here develop a more flexible modeling framework based on the class of max-infinitely divisible processes, which extend max-stable processes while retaining dependence properties that are natural for maxima. We propose two parametric constructions for max-infinitely divisible models, which relax the max-stability property but remain close to some popular max-stable models obtained as special cases. The first model considers maxima over a finite, random number of independent observations, while the second model generalizes the spectral representation of max-stable processes. Inference is performed using a pairwise likelihood. We illustrate the benefits of our new modeling framework on Dutch wind gust maxima calculated over different time units. Results strongly suggest that our proposed models outperform other natural models, such as the Student-t copula process and its max-stable limit, even for large block sizes

Quantitative Regular Expressions for Arrhythmia Detection Algorithms

Motivated by the problem of verifying the correctness of arrhythmia-detection algorithms, we present a formalization of these algorithms in the language of Quantitative Regular Expressions. QREs are a flexible formal language for specifying complex numerical queries over data streams, with provable runtime and memory consumption guarantees. The medical-device algorithms of interest include peak detection (where a peak in a cardiac signal indicates a heartbeat) and various discriminators, each of which uses a feature of the cardiac signal to distinguish fatal from non-fatal arrhythmias. Expressing these algorithms' desired output in current temporal logics, and implementing them via monitor synthesis, is cumbersome, error-prone, computationally expensive, and sometimes infeasible. In contrast, we show that a range of peak detectors (in both the time and wavelet domains) and various discriminators at the heart of today's arrhythmia-detection devices are easily expressible in QREs. The fact that one formalism (QREs) is used to describe the desired end-to-end operation of an arrhythmia detector opens the way to formal analysis and rigorous testing of these detectors' correctness and performance. Such analysis could alleviate the regulatory burden on device developers when modifying their algorithms. The performance of the peak-detection QREs is demonstrated by running them on real patient data, on which they yield results on par with those provided by a cardiologist.Comment: CMSB 2017: 15th Conference on Computational Methods for Systems Biolog

Remote sensing of coccolithophore blooms in selected oceanic regions using the PhytoDOAS method applied to hyper-spectral satellite data

In this study temporal variations of coccolithophore blooms are investigated using satellite data. Eight years (from 2003 to 2010) of data of SCIAMACHY, a hyper-spectral satellite sensor on-board ENVISAT, were processed by the PhytoDOAS method to monitor the biomass of coccolithophores in three selected regions. These regions are characterized by frequent occurrence of large coccolithophore blooms. The retrieval results, shown as monthly mean time series, were compared to related satellite products, including the total surface phytoplankton, i.e. total chlorophyll a (from GlobColour merged data) and the particulate inorganic carbon (from MODIS-Aqua). The inter-annual variations of the phytoplankton bloom cycles and their maximum monthly mean values have been compared in the three selected regions to the variations of the geophysical parameters: sea-surface temperature (SST), mixed-layer depth (MLD) and surface wind-speed, which are known to affect phytoplankton dynamics. For each region, the anomalies and linear trends of the monitored parameters over the period of this study have been computed. The patterns of total phytoplankton biomass and specific dynamics of coccolithophore chlorophyll a in the selected regions are discussed in relation to other studies. The PhytoDOAS results are consistent with the two other ocean color products and support the reported dependencies of coccolithophore biomass dynamics on the compared geophysical variables. This suggests that PhytoDOAS is a valid method for retrieving coccolithophore biomass and for monitoring its bloom developments in the global oceans. Future applications of time series studies using the PhytoDOAS data set are proposed, also using the new upcoming generations of hyper-spectral satellite sensors with improved spatial resolution

Nonnegative unmixing methodology applied on Brillouin Optical Fiber Sensor

International audienceAs a complement to conventional sensors, Distributed Optical Fiber Sensors (DOFS) have gradually played a prominent role in Structural Health Monitoring (SHM) for the last decade. The distributed Brillouin sensor enables to measure strain along kilometers of cable with a spatial resolution of 1 meter. The challenge is to have a centimeter spatial resolution to improve structure defaults detection and localization. A numerical model, based on the sensor physic, is first proposed in order to study Brillouin spectra distortion depending on strain distribution within spatial resolution. Then, based on nonnegative least squares (NNLS) problem, Brillouin spectra are decomposed into several elementary spectra. The estimated central frequencies and maxima permit to estimate a centimeter frequency distribution within the spatial resolution. It has been verified with numerical and experimental data: that method enables to enhance the accuracy and spatial resolution of the sensor from meter to centimeter

Earthquake Arrival Association with Backprojection and Graph Theory

The association of seismic wave arrivals with causative earthquakes becomes progressively more challenging as arrival detection methods become more sensitive, and particularly when earthquake rates are high. For instance, seismic waves arriving across a monitoring network from several sources may overlap in time, false arrivals may be detected, and some arrivals may be of unknown phase (e.g., P- or S-waves). We propose an automated method to associate arrivals with earthquake sources and obtain source locations applicable to such situations. To do so we use a pattern detection metric based on the principle of backprojection to reveal candidate sources, followed by graph-theory-based clustering and an integer linear optimization routine to associate arrivals with the minimum number of sources necessary to explain the data. This method solves for all sources and phase assignments simultaneously, rather than in a sequential greedy procedure as is common in other association routines. We demonstrate our method on both synthetic and real data from the Integrated Plate Boundary Observatory Chile (IPOC) seismic network of northern Chile. For the synthetic tests we report results for cases with varying complexity, including rates of 500 earthquakes/day and 500 false arrivals/station/day, for which we measure true positive detection accuracy of > 95%. For the real data we develop a new catalog between January 1, 2010 - December 31, 2017 containing 817,548 earthquakes, with detection rates on average 279 earthquakes/day, and a magnitude-of-completion of ~M1.8. A subset of detections are identified as sources related to quarry and industrial site activity, and we also detect thousands of foreshocks and aftershocks of the April 1, 2014 Mw 8.2 Iquique earthquake. During the highest rates of aftershock activity, > 600 earthquakes/day are detected in the vicinity of the Iquique earthquake rupture zone

The use of imaging systems to monitor shoreline dynamics

The development of imaging systems is nowadays established as one of the most powerful and reliable tools for monitoring beach morphodynamics. Two different techniques for shoreline detection are presented here and, in one case, applied to the study of beach width oscillations on a sandy beach (Pauanui Beach, New Zealand). Results indicate that images can provide datasets whose length and sample interval are accurate enough to resolve inter-annual and seasonal oscillations, and long-term trends. Similarly, imaging systems can be extremely useful in determining the statistics of rip current occurrence. Further improvements in accuracy and reliability are expected with the recent introduction of digital systems