# Dissolved Radiotracers and Numerical Modeling in North European Continental Shelf Dispersion Studies (1982-2016): Databases, Methods and Applications

Pascal Bailly Du Bois, Franck Dumas, Claire Voiseux, Mehdi Morillon, Pierre-Emmanuel Oms, Luc Solier

► To cite this version:

Pascal Bailly Du Bois, Franck Dumas, Claire Voiseux, Mehdi Morillon, Pierre-Emmanuel Oms, et al. Dissolved Radiotracers and Numerical Modeling in North European Continental Shelf Dispersion Studies (1982-2016): Databases, Methods and Applications. Water, 2020, 12 (6), pp.1667. 10.3390/w12061667. [https://hal.science/hal-02866940v1](https://hal.science/hal-02866940v1)

HAL Id: hal-02866940

Submitted on 12 Jun 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
<figure><img src="image_2.png" /></figure>

Review

# Dissolved Radiotracers and Numerical Modeling in North European Continental Shelf Dispersion Studies (1982–2016): Databases, Methods and Applications

Pascal Bailly du Bois <fnref n="1" />*, Franck Dumas <fnref n="2" />, Claire Voiseux <fnref n="1" />, Mehdi Morillon <fnref n="1" />, Pierre-Emmanuel Oms <fnref n="3" /> and Luc Solier <fnref n="1" />

Received: 2 October 2019; Accepted: 5 June 2020; Published: 10 June 2020

**Abstract:** Significant amounts of anthropogenic radionuclides were introduced in ocean waters following nuclear atmospheric tests and development of the nuclear industry. Dispersion of artificial dissolved radionuclides has been extensively measured for decades over the North-European continental shelf. In this area, the radionuclide measurement and release fluxes databases provided here between 1982 and 2016 represent an exceptional opportunity to validate dispersion hydrodynamic models. This work gives accessibility to these data in a comprehensive database. The MARS hydrodynamic model has been applied at different scales to reproduce the measured dispersion in realistic conditions. Specific methods have been developed to obtain qualitative and quantitative results and perform model/measurement comparisons. Model validation concerns short to large scales with dedicated surveys following the dispersion: it was performed within a two- and three-dimensional framework and from minutes and hours following a release up to several years. Results are presented concerning the dispersion of radionuclides in marine systems deduced from standalone measurements, or according to model comparisons. It allows characterizing dispersion over the continental shelf, pathways, transit times, budgets and source terms. This review presents the main approaches developed and types of information derived from studies of artificial radiotracers using observations, hydrodynamic models or a combination of the two, based primarily on the new featured datasets.

**Keywords:** radionuclide; tracer; data collection; antimony 125 (¹²⁵Sb), tritium (³H), dispersion; modeling; English Channel; North Sea; Biscay Bay

## 1. Introduction

Significant amounts of anthropogenic radionuclides have been introduced in ocean waters since 1945. The main origins were the fallout from the atmospheric nuclear tests that occurred before 1980, controlled releases from the nuclear industry and the accidental releases resulting from the Chernobyl (1986) and Fukushima Dai-ichi (2011) nuclear power plants. Among these sources, in

<footnote n="1">¹ Laboratoire de Radioécologie de Cherbourg, IRSN-LRC, rue Max Pol Fouchet B.P. 10, 50130 Cherbourg en Cotentin, France; claire voiseux@irsn.fr (C.V.); mehdi.morillon@gmail.com (M.M.); luc.solier@irsn.fr (L.S.)</footnote>

<footnote n="2">² Service Hydrographique et Océanographique de la Marine-SHOM, 29200 Brest, France
franck.dumas@shom.fr</footnote>

<footnote n="3">³ Centre d'Etudes et de Valorisation des Algues CEVA, 83 Presqu'île de Pen Lan, 22610 Pleubian, France;
omspe@hotmail.fr</footnote>

<footnote>* Correspondence: pascal.bailly-du-bois@irsn.fr; Tel.: +33-2-33-01-41-05</footnote>
Europe, the releases from nuclear fuel reprocessing plants from Sellafield and La Hague were the most important.

Oceanographic sampling campaigns made it possible to measure the dispersion of the artificial radionuclides remaining dissolved in seawater over the past 30 years. These works have demonstrated the interest of dissolved radionuclides as tools for oceanography [1–3]. Over the European continental shelf, extensive in-situ measurements have been performed between 1950 and 2000 by English, German, Belgian and French institutes [4,5] that allowed to draw up the general circulation pathways and water masses transit times in the Irish Sea, the North Sea, the English Channel (Figure 1) and the Arctic ocean [6–11]. The Cherbourg Radioecology Laboratory (IRSN-LRC) has contributed to these studies at the scale of the English Channel and the North Sea since 1988 [12–21]. These works have been ongoing since 1994 in the North-East Atlantic waters, the Celtic Sea, Irish Sea and the Bay of Biscay. During the 1990s, it appears that dissolved radionuclide measurements represent exceptional tools to test and validate marine hydrodynamic models applied to represent short to long term dispersion processes.

This work presents a review and an update of results from the main dissolved radionuclides that have been measured extensively by the IRSN-LRC, i.e., ³H, ¹²⁵Sb, ¹³⁷Cs, ¹³⁴Cs, ¹⁰⁶Ru and ⁶⁰Co. Other radionuclides could be used as oceanographic tracers, such as ⁹⁰Sr, ⁹⁹Tc, ¹²⁹I, ²³⁸, ²³⁹+²⁴⁰Pu. Radionuclides with a half-life lower than one year has been sparsely detected as ¹¹⁰mAg, ⁵⁴Mn, ⁵⁸Co and ¹³¹I.

Realistic simulations of the dispersion of soluble substances in the marine environment are essential for management of the marine ecosystem. Such simulations were applied to study the fate of chronic or accidental releases into the sea; they may also be used to feed ecological models that encompass exchanges between the different compartments of the environment: seawater, living organisms and sediments. Such tools are particularly relevant in seas that are subject to strong anthropogenic pressures, such as the macro-tidal seas of north-western Europe.

Various methods have been tested for calculating the behavior of water masses, their advection and their dispersion. Models commonly simulate currents at different resolutions (from a tenth of a meter to kilometers), spatial (1 km–1000 km) and temporal coverage (from hour to decades). While these models produce an accurate representation of tidal levels and associated currents, greater requirements are needed to simulate advection of soluble substances over periods longer than the tidal cycle. The models’ ability to reproduce dispersion under realistic conditions of release, wind and tide over several days, weeks or years is a sensitive criterion for assessing their reliability.

Validation of hydrodynamic models applied for realistic simulation of mid- to long term dispersion in seawater requires field data of comparable parameters and coverage. The ideal tracer must have a conservative behavior in the water mass; that is to say, neither fixed by the environmental compartments (sediment, living species) nor modified during its stay in seawater and when subsequently diluted. It must be measurable several hundred or thousand kilometers from its input point (meaning even at very low levels of concentration). The discharge conditions and flow must be well known and it must have few properly controlled origins. The radioactive decay is easily accounted by hydrodynamic models.

Some artificial (issued from nuclear industry) radionuclides released by nuclear plants fully meet these specifications if their half-life is long enough compared to the transit-times in the studied area (from weeks to years). Among them, ¹²⁵Sb and ³H as HTO have proven to be conserved in seawater over years at the scale of European waters. These radionuclides can be measured at very low concentrations, up to the Atlantic seawater background concentrations (concentrations 40,000 times lower than natural radioactivity for gamma-emitters).

Collaboration between Ifremer-DYNECO-PHYSED oceanographic physicians and IRSN-LRC marine radio ecologists contributes to improve the marine hydrodynamic models used [22–29]. Further studies have associated systematically in-situ measurements with model simulations in order to improve the efficiency of measurements to check model’s precision and obtain more reliable and versatile models, applicable to all coastal seas of the European continental shelf. Model/measurements comparisons have been performed at the scale of the English Channel and the
North Sea with targeted two-dimensional (2D) residual and instantaneous models [11,22–24]. This was done in the Bay of Biscay [26] and the Pacific [27] with three-dimensional (3D) models.

This study also included high resolution 2D and 3D model/measurements comparisons at short scale close to a release outfall [25].

<figure><img src="image_3.png" /><figcaption>Figure 1. General circulation of water masses in the north-west of Europe, adapted from [11].</figcaption></figure>

The purpose of this work is to present a review and an update of field data and methods applicable to validate and improve dispersion models in European macro-tidal seas at all scales. The homogeneous database thus gathered can be used by modelers to test the reliability of their models against appropriate data. It complements the existing databases as the IAEA Marine Information System (MARIS [30]), the World Ocean Database 2013 [31,32], IFREMER-SISMER [33] or BODC [34].

This database has an historical value as releases from nuclear plants have significantly decrease for most gamma emitters (two orders of magnitudes for <sup>125</sup>Sb, <sup>106</sup>Ru, <sup>137</sup>Cs and <sup>60</sup>Co presented here); it represents a huge amount of work (80 oceanographic campaigns) and the dispersion plumes measured before 2000 will be difficult to obtain in the future. Moreover, in a context of a decreasing number of marine radio ecologists, it is valuable to share these data with a larger scientific community. We are convinced that they have not given all the knowledge they could provide.

This work is mainly focused on the dispersion of the Orano recycling plant located at La Hague in the mid English Channel. It represents the main source-term of dissolved radionuclides in this sea. Consequently, its dispersion concerns mostly the English Channel and Southern North Sea.
Nevertheless, other radiotracers source terms have been accounted as Sellafield releases in the Irish Sea or nuclear power plant releases reaching the Bay of Biscay through rivers. It could represent a metrological challenge, but tracers exist to extent the work elsewhere in the world.

This work addresses the following successive aspects: the background section describes successively the radionuclides measured, sampling and measurement methods; the measurement and release database achieved; a short description of the models used and the different model/measurement comparisons performed.

The application section focuses on the main features that could be retained from the point of view of radiotracers, hydrodynamic models and methods, from the short scale in the vicinity of an outfall, to the large scale in the English Channel, North Sea, the Biscay Bay and in the Pacific (Fukushima accident). The last section presents perspectives of applications in other areas or oceanographic domains.

## 2. Background

### 2.1. Dissolved Radionuclides as Oceanographic Tracers

To study the dispersion process over short to long periods or carry out and interpret repeated measurements for varying conditions of release or forcing, it is necessary to use tracers which fulfill the following characteristics:

1. It must originate from one or a small number of clearly identified release points;

2. The release conditions must be precisely known (time, fluxes);

3. Labeling in seawater must be significant, in particular in relation to the pre-existing background level. Labeling concerns seawaters where a significant concentration of radionuclide could be measured.

4. The tracer must be soluble and not fix onto living organisms or sediments over time (i.e., the stable element is conservative in seawater);

5. It must be possible to measure the tracer after dilution in the sea over hours, weeks, months or years after its release.

Natural tracers, such as copper, iron, nutrients, cannot generally be used because of the multiplicity of their source terms and the complexity of the phenomena governing their production and fate in the marine environment.

Even if the release conditions are known for artificial tracers, there are often many release points for each one. Such tracers are often involved in geochemical and biogeochemical processes, and therefore their conservative behavior in the marine environment are not guaranteed as shown in Section 3.2.4. for <sup>60</sup>Co and <sup>106</sup>Ru.

Radioactive tracers generally meet criteria 1 and 2. As regards criteria 3, 4 and 5, some radionuclides exhibit long term conservative behaviors in seawater, such as <sup>125</sup>Sb, <sup>99</sup>Tc and <sup>3</sup>H (tritium) [35], and to a lesser extent <sup>90</sup>Sr, <sup>137</sup>Cs and <sup>134</sup>Cs. Between 1970 and 1995 extensive measurements of these radionuclides were carried in the seas of north-western Europe [5,6,11,14–21,24,36]. Reductions in fluxes released during the period from 1980 to 2000 have led to significant decreases in concentrations in the marine system. Out of the radionuclides mentioned above, only tritium released from nuclear fuel re-processing plants has not undergone a reduction since 1980 (Figure 2); nevertheless, it fully satisfies the five criteria, having a radioactive half-life of 12.4 year.
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Annual liquid releases from Sellafield and La Hague nuclear reprocessing plants for ¹²⁵Sb, ³H, ¹³⁷Cs, ¹⁰⁶Ru and ⁶⁰Co.</figcaption></figure>

An additional interest of artificial radionuclides as environmental tools is that there is no risk of samples contamination during in-board sampling process, which is generally not the case for other chemicals that require clean conditions to prevent pollutions.

The radionuclides measured in this study that have the most conservative behavior in seawater are mainly antimony 125 (¹²⁵Sb) and tritium (³H) in the form of tritiated water HTO (Table 1). The radioactive emission determines the technic that could be applied for each radionuclide.

<table><caption>**Table 1.** Characteristics of the radionuclides used in this work.</caption><thead><tr><th>Radionuclide</th><th>Tritium<br>³H</th><th>Cesium<br>137<br><sup>137</sup>Cs</th><th>Cesium<br>134<br><sup>134</sup>Cs</th><th>Antimony<br>125<br><sup>125</sup>Sb</th><th>Ruthenium<br>106<br><sup>106</sup>Ru</th><th>Cobalt<br>60<br><sup>60</sup>Co</th></tr></thead><tbody><tr><td>Radioactive decay (year)</td><td>12.3</td><td>30.2</td><td>2.1</td><td>2.8</td><td>1</td><td>5.3</td></tr><tr><td>Principal radioactive<br>emission</td><td>β</td><td>γ</td><td>γ</td><td>γ</td><td>γ</td><td>γ</td></tr><tr><td>Conservative behavior *</td><td>100%</td><td>83%-86%</td><td>83%-86%</td><td>98%</td><td>19%-26%</td><td>8%-14%</td></tr></tbody></table>

### 2.1.1. Gamma Emitters

Gamma emitter's (¹³⁷Cs, ¹³⁴Cs, ¹⁰⁶Ru and ⁶⁰Co in Table 1) concern radionuclides which the main radioactive emission is a high energy γ photon. Thus, they can be detected all together by their

<footnote>*: Percentage of radionuclide quantities remaining dissolved in seawater at the scale of the English Channel (from [20]).</footnote>
specific energy emission without radiochemistry for the isolation of each element. Specific methods adapted to the dilution conditions observed in the marine environment and the on-board treatment of samples was developed in order to measure artificial activities of the order of 0.3 Bq·m⁻³ [37,38]. This is in contrast with the natural radioactivity of seawater, which is of the order of 12,000 Bq·m⁻³ (mainly ⁴⁰K).

This work focuses mainly on ¹²⁵Sb and ³H which are conservative in seawater; nevertheless, the database provides measurements of the radionuclides, which are measured together (¹⁰⁶Ru, ¹³⁷Cs, ¹³⁴Cs and ⁶⁰Co). They are useful to investigate other process as exchanges between seawater and sediments. Methods to measure ¹²⁵Sb, ¹⁰⁶Ru, ¹³⁷Cs, ¹³⁴Cs and ⁶⁰Co have been described in length in [37,38]. Analyses were carried out at the French Navy-Groupe d'études atomiques in Cherbourg (GEA), and at the IRSN-LRC laboratory in Cherbourg.

### 2.1.2. Tritium

³H is present in all nuclear plants' liquid releases. The two main sources of tritium are the nuclear fuel reprocessing plants at Sellafield into the Irish Sea that have been active since 1952, and from La Hague into the English Channel since 1966. Liquid releases from each nuclear power plant are two orders of magnitude lower than the fluxes issued from reprocessing plants. Nevertheless, these releases must also be accounted for all along the European coasts and rivers, particularly away from the reprocessing plants influences.

Due to the small amount of seawater required for direct measurement by liquid scintillation (8 mL) a large amount of samples can be processed that can be used to have a high frequency 2D picture or even a 3D picture of the tritium dispersion at a short scale close to the La Hague outfall (Sections 2.4.1. and 2.4.2. [26,39]). The power plants labeling could be used in areas out of theses influences, with corresponding levels lower than 1 Bq·L⁻¹, in the Bay of Biscay [26] or in the Mediterranean Sea.

In the releases performed by nuclear reprocessing plants and French nuclear power plants, tritium is the form of the tritiated water molecule HTO; thus, it has a strict conservative behavior in seawater. This is not always the case: for example, the releases performed by the Amersham plant in Cardiff in the Bristol Channel concern organic compounds tagged with tritium. Organically bounded tritium (OBT) has a specific behavior in marine environment with strong assimilation by living species [40]. Due to their locations (mainly the Bristol Channel and Rhone River), these kinds of releases do not influence the results presented here.

The most used method to measure tritium is liquid scintillation, which allows up to 10,000 samples per year for one analyzer. With a detection limit around 1 Bq·L⁻¹, this method could only be relevant within the plume of reprocessing plant releases (eastern English Channel, Irish Sea, rivers). Extra suited methods exists for tackling lower levels of concentration, which requires sampling one liter of seawater. Low level ³H measurements were performed with two methods:

(i) Electrolytic enrichment of water samples [41]. The detection limit reached 0.01 Bq·L⁻¹.

(ii) ³He regrowth and measurement by mass spectrometry [42]. The detection limit could reach 0.001 Bq·L⁻¹.

## 2.2. Database

### 2.2.1. Radionuclide Measurements in Seawater

The measurement database provided includes all oceanographic campaigns performed by IRSN-LRC (CEA-LRM before 2002) between 1982 and 2016 (Table A1). Data are provided for ³H, ¹²⁵Sb, ¹⁰⁶Ru, ¹³⁷Cs, ¹³⁴Cs and ⁶⁰Co. The database concerns 80 oceanographic campaigns and totals 39,642 sampling locations at sea. Data concerning the coastal monitoring station the closest to the La Hague outfall were added (744 measurements between 1984 and 2018).

Table A1 in Appendix A lists the different campaigns, the number of measurements and maximum concentrations measured. Figure 3 shows the locations of all samples obtained. The
measurement database is available in [43] as a supporting material for this work. Part of these data was already given available in [5] (2010 ¹³⁷Cs measurements), in [26,44] (14,494 tritium measurements), and sparsely in previous publications. The new database encompasses 47,387 measurements; more than 60% of these data were previously unpublished.

As said in the introduction, existing databases include measurements of dissolved radioactivity in the European seas. The most important are the IAEA Marine Information System (MARIS [30]), the World Ocean Database 2013 [31,32], the IFREMER-SISMER [33] or the BODC [34]. The data provided here complete them in different ways.

¹³⁷Cs and ¹³⁴Cs are commonly measured in oceanographic studies, but ¹²⁵Sb, ¹⁰⁶Ru and ⁶⁰Co are sparsely detected. ³H was used as a tracer in open oceans but to a lesser extent in coastal waters.

Due to their conservative behavior in seawater, ¹²⁵Sb and ³H are choice tracers to track the dispersion of industrial releases.

The data collection was designed to appraise the dispersion of French nuclear plants releases from short to large spatial and temporal scales through repeated oceanographic campaigns since 1982. The La Hague controlled radioactive releases in seawater are the more important in Europe; the corresponding radiotracer measurement database allows a complete case study of its dispersion in marine systems.

To our knowledge, no other in-situ measurement database allows to follow individual radionuclides releases in the marine environment at short scale.

The provided database is homogeneous with the release data; it could be included further in other existing databases.

<figure><img src="image_5.png" /><figcaption>Figure 3. Location of radionuclide measurements samples obtained between 1982 and 2016.</figcaption></figure>
### 2.2.2. Radioactive releases

Releases in seawater from La Hague reprocessing plant was transmitted by Orano company for each individual release since 1982 up to 2018. Quantity released, date-time of the beginning and end of each release are given. Prior to 1982, only annual releases were available. Part of these data was already given available in [26,45] (7840 tritium releases).

The new database encompasses 22,183 individual releases for tritium, ¹³⁷Cs, ¹³⁴Cs, ¹²⁵Sb, ¹⁰⁶Ru and ⁶⁰Co (110,915 values). In total, 93% of these data are unpublished as part of a comprehensive accessible database until now.

The Sellafield reprocessing plant annual release data have been published by [46] and in MAFF, RIFE and Sellafield Ltd. reports [47,48]. ¹²⁵Sb data have been available only since 1978 with an information gap between 2001 and 2012.

Annual releases from British and French nuclear power plants are issued from annual reports from MAFF, RIFE and EDF reports [47,48].

The releases database is available in [49] as a supplementary materials for this work. Figure 2 presents the annual fluxes released by Sellafield and La Hague reprocessing plants for ¹²⁵Sb, ¹³⁷Cs, ¹⁰⁶Ru and ⁶⁰Co. It shows that the main release period for ¹²⁵Sb, ³H, ¹³⁷Cs, ¹⁰⁶Ru and ⁶⁰Co occurs before 1992, with a maximum for ¹³⁷Cs between 1969 and 1983 at Sellafield (1–5 PBq·y⁻¹), and for ¹⁰⁶Ru between 1969 and 1991 at La Hague (close to 1 PBq·y⁻¹). Tritium releases (in blue) show a more homogeneous temporal evolution, with fluxes in the same magnitude as Sellafield since 1969 (1–3 PBq·y⁻¹) and since 1993 from La Hague (around 10 PBq·y⁻¹). It depends mainly from the quantities of nuclear fuel processed by the plants.

Existing datasets concern mainly annual releases from selected nuclear plants. This work provides a compilation of all the release data in one accessible database for French and British nuclear plants in European seas. The knowledge of each individual release is essential to perform precise comparisons between measurements and hydrodynamic model simulations and evaluate the dispersion in the vicinity of an outfall. The detail of the La Hague controlled radioactive releases represent an unparalleled dataset in this perspective.

## 2.3. *Hydrodynamic Models*

The different models applied are based on the Model for Applications at Regional Scale (MARS) developed by the French Ifremer institute since 1987 [50]. This model was built under various assumptions presented below: using the non-stationary Saint-Venant equations (i.e., 2D); the primitive 3D equations; the "Lagrangian barycentric" method in 2D to filter out the tidal signal.

### 2.3.1. 2D Modeling

Numerous modeling studies [20,22–25] have demonstrated that models using two-dimensional horizontal approximation (i.e., shallow-water equations) are able to simulate a satisfactory dissolved-substance transport in the non-stratified area. These equations were solved using the finite-difference MARS2D model [50].

In the largest domain the baroclinic effects can be neglected (so that 2D models can relevantly be used) as long as they do not influence that much the targeted area. Naturally, they play a major role over the shelf of the Bay of Biscay or next to the mouth of large estuaries (Loire, Gironde, Rhine etc.). In the eastern English Channel and southern North Sea, away from rivers, the ocean is rarely thermally stratified nor stratified in terms of salinity and temperature; this is due to the strong tidally induced mixing.

### 2.3.2. 3D Modeling

The three-dimensional MARS3D model is described in detail in [50]. The model uses a 3D finite difference scheme, applying the Boussinesq approximation and hydrostaticity to resolve primitive equations. The model involves a nesting strategy (example Figure 4), starting from a broad region
covering the entire North-West European continental shelf (with a 5.6 km grid resolution) down to a
detailed domain covering a few tens of km (with a 5–100 m resolution).

The bathymetry at the grid nodes of the different models is estimated with the method
described in [51] from various data sources [25,52].

<figure><img src="image_6.png" /><figcaption>**Figure 4.** Example of model nesting for the La Hague Cape area, Ranks 0–2 are two-dimensional (2D), Rank3 is 2D or three-dimensional (3D), Rank4 is 3D, adapted from [26].</figcaption></figure>

2.3.3. 2D "Lagrangian Barycentric" Method

A tidal residual model was designed to reproduce transport, dilution and decay phenomena over long time scales (ranging from a week to several years) and extensive spatial scales (from 30 to 1000 km). It is based on the hypothesis that the water column is homogeneous and that barotropic phenomena will prevail over dynamic baroclinic phenomena. This model has been extensively described and applied in [18,20,25,52–54]. The interest of applying "Lagrangian barycentric" currents is to investigate the residual circulation occurring at long time scales (more than a tidal cycle of 12.42 h), and to simulate dispersion very quickly at the scale of the English Channel and the North Sea (computation time about 1000 times lower than with a similar 2D model).

## 2.4. *Sampling Close to an Outfall*

2.4.1. Model Assisted Sampling

In the vicinity of a punctual source-term (represented by the location of the end of the release
pipe), the sampling strategy must be adapted to catch the dispersion of a rapidly moving narrow
plume exhibiting short scale features.

During the first hours following release, it is necessary to have good knowledge of the time
schedule of the release and where the plume will be located close to the outfall. The dynamics of the
currents and duration of releases close to the La Hague Cape impose precise positioning in time (less
than 15 min) and space (100 m), to ensure that the sampling sections of the ship's track encompass
the plume area (Figure 5 [25]). The releases schedules were transmitted to the vessel by the nuclear
plant unit (Service de Protection Radiologique-SPR) of the ORANO company. The MARS2D model
(see Section 2.3.1., rank 3 in Figure 4) is used on-board to get some forecasts of the dispersion ahead
of the release and allow precise positioning of the vessel during the hours following the start of the
release.

The procedure made possible to sample at high frequency (30 s) an average of 300 surface
stations to follow the horizontal dispersion of a given release (see for example the figure in Section
3.1.3.). Varying hydrodynamic and meteorological conditions were investigated.

<figure><img src="image_7.png" /><figcaption>**Figure 5.** Sampling locations during the DISVER 2011 campaign (13,400 in-depth samples between 5 and 11 April 2011) (a) 0–2 hours after a release (b) 3–8 hours after a release.</figcaption></figure>

2.4.2. High Frequency In-Depth Sampling

Sampling at depth during oceanographic campaigns requires usually bottles deployed from the surface down to the sea floor on a manual or monitored sampling system, such as a CTD rosette. This requires the ship stops at stations during sampling event (i.e., move the sampling system downward, close bottles and move upward at around 1 m·s⁻¹ speed). Overall, it takes about 15–30 min for each station in shallow waters (20–50 m). In the working area, for security reasons the ship cannot be stopped close to the coast or in areas with strong currents. Furthermore, the plume identification requires much higher sampling rates. In the context of the DISVER experiment (Section 3.1.2. [55]) for studying vertical dispersion close to an outfall, we aim to validate three-dimensional hydrodynamic models by using in-situ in-depth radiotracer measurements; this objective requires high-frequency sampling (one every 30 s) at 10 depth levels simultaneously (thus giving 1800 samples per hour).

The area off the Cap de La Hague represents one of the most difficult areas to investigate for such studies, since it combines strong currents (up to 5 m·s⁻¹), proximity of the coast (less than one km) and a complex topography with many rocky shoals and deeps varying from 25 to 90 m over a few kilometers. For operational and security reasons, the ship must remain under way normally at speeds from 1 to 5 m·s⁻¹ (10 knots) while sampling. A devoted system was developed to perform safe sampling in such rough environment [39].

It comprises three main components:

* A sampling line designed to sample 10 depths simultaneously down to 65 m;

* A deep towed depressor (known as Dynalest), which maintains the line at depth close to the seabed;

* An automatic high-frequency sampler with volume, flux and depth control to get 1800 samples per hour.
To compare several campaigns with variable source terms, measured concentrations could be normalized by the cumulated released flux corresponding to the plume targeted. They become dilution coefficients (DC) or dilution factors with:

$$DC = \frac{Bq \cdot L^{-1}_M}{Bq \cdot L^{-1}_R} \qquad (5)$$

Bq·L⁻¹M: measured concentrations

Bq·L⁻¹R: released flux

This dilution coefficient was applied for model/measurement comparisons at a short scale in the vicinity of an outfall, when individual release plumes could be distinguished [25]. At larger scales, individual releases are mixed together in seawater. A more representative approach is to compare the measured concentration (Bq·m⁻³) with the mean released flux in Bq·s⁻¹. Such dilution coefficients (concentrations corresponding to a given released flux) were applied at the scale of the whole English Channel and the whole North Sea (see Section 3.2.3., [20,24]).

For 3D model/measurement comparisons, the methods described previously could be applied by accounting for the water mass stratification of the radionuclide concentrations. In case of studies in the near outfall field, particular attention should be given to obtain the same reference of x, y, z and t between measurements and simulations: an error of minutes or hundred meters could result in a concentration variation larger than an order of magnitude.

### 2.5.2. Radionuclides Inventories

When campaigns had a sufficient extent and sampling rate, localized radioactivity measurements were interpolated over the whole studied area at the nodes of a regular grid. A well-suited interpolation method is “kriging” [57]. This allows the visualization and comparison of the distribution plumes measured and simulated (see Section 3.2.1.).

If the radionuclide distribution in seawater is not stratified (no variation from surface to sea bottom), the total inventory of radionuclides present could be simply calculated by the multiplication of the surface concentration C by the volume of seawater of each mesh (C × dx × dy × depth). In case of stratification, in-depth measurements must be performed in order to associate the right concentration to each level of the water column [26].

Quantities of measured radionuclides can be straightforwardly compared to the simulated ones at different scales [18,20,24,26,27] or compared to the known releases in order to exhibit transit times or conservative behavior [18]. Inverse calculation has also been applied to estimate the source term (e.g., in the Fukushima case) [58].

### 2.6. Scales for Model/Measurement Comparisons

Different time and space scales are presented in this work; the distinction of the different scales is based on the dispersion characteristics in a macro-tidal context. Super tidal scale (minutes to hours, 100 m–10 km) concerns dispersion close to an outfall when the labeled plume is not vertically averaged. The tidal scale (hours to days, 1–30 km) concerns the dispersion in the vicinity of an outfall when the labeled plumes associated to each specific release could be distinguished by surface measurements. Large scale (from week to years, 30–1000 km) concerns the dispersion when individual plumes are mixed together. This is discussed further in Section 4.3.

## 3. Applications

The obtained database [43] contains data from 1982 up to 2016, with all measurements acquired by the IRSN-LRC laboratory during oceanographic campaigns. It concerns mainly ¹²⁵Sb, ³H, ¹⁰⁶Ru, ¹³⁷Cs, ¹³⁴Cs and ⁶⁰Co.

Such data would have been meaningless without the corresponding known releases that explain the observed labeling (cf. 2.2.2. [49]).
## 3.1. La Hague Cape, Short Scale/High Resolution Studies

### 3.1.1. The La Hague Cape Main Characteristics for Dissolved Radionuclide Dispersion

The North Cotentin includes several nuclear facilities, such as the nuclear power plant at Flamanville Cape, the building of nuclear submarines at Cherbourg and the ANDRA low-level waste deposit at La Hague. The most important in term of liquid releases at sea is the Orano reprocessing plant, of which the outfall is located next to the La Hague Cape (Figure 5). The La Hague Cape forms a physical boundary between the Normandy-Brittany Gulf in the south-west, and the mid-English Channel towards the east. Because of the coastal morphology, the tidal wave coming from the Atlantic is partly blocked in the west-facing bay formed by the Normandy-Brittany Gulf. This embayment is characterized by very large tidal ranges (reaching more than 14 m near the Mont St Michel during spring tides). The Cape de la Hague works as a bottleneck for the water masses involved during the emptying and filling of this bay twice a day. This explains why the tidal currents close to the cape are among the strongest in Europe (they can reach 5 m·s⁻¹ during spring tide), with highly variable tidal range around the cape (from 11 m in the south of the Cape down to 6 m in the north [59]). This area is also characterized by diverse topography exhibiting pronounced gradients (depths from 20 to 100 m), many islands, numerous bays and shallow coves.

A tidal residual currents [52] divergence zone close to the release outfall divides waters flowing into the Normandy-Brittany Gulf from the waters forming part of the general flow from west to east up the Channel and towards the Straits of Dover. As a result, small differences in the release conditions can lead to opposite directions of spreading in the medium-term [25,60]. As a consequence, simulation of the dispersion in the area is challenging for the numerical models' dispersion capability.

### 3.1.2. 3D Dispersion: Super Tidal Time Scale

The La Hague outfall is located on the sea floor, two kilometers off the coast (Figure 5). The vertical dispersion of the plume has been investigated during the DISVER project. In total, 19,000 in-depth samples were taken during three campaigns of four-five days [55]. Figure 5 shows an example of the sampling strategy with plume transects at different distances from the outfall. The transects are slices across the axis of the plume propagation from the surface down to 25 m depth; they are spaced at time intervals of 5 to 10 min. The database encompasses 137 vertical transects up to 65 m depth with 100-200 individual measurements per transect. This provides a rather good picture of the vertical structure of the plume along each of these transects. Figure 6 shows an example of model/measurement comparison during two series of transects at 850-1300 m and 3200-4200 m downstream the outfall. These figures show an accurate representation of the plume in space and time. Figure 6a exhibits unexpected highly complex and variable structures resulting from intense turbulent mixing with eddies of about 100 m that extend to the width of the plume rapidly. The 3D model is not able to reproduce the turbulent mixing at this scale. Further than 3000 m downstream (Figure 6b), the plume reaches the sea surface with a more homogeneous shape that is better reproduced by the model.

The results obtained from all transects have been integrated by the normalization of the concentration measured with the corresponding release. Results are presented in Figure 7 with the measured and simulated extension and dilution of the plume with the distance during a constant release (steady state situation). It shows that if the model does not reproduce the instantaneous turbulence, on average, the shape and dilution coefficient are given the correct order of magnitude. Details of the dilution variations and model/measurement comparisons are given in [55].
<figure><img src="image_8.png" /></figure>
<figure><img src="image_9.png" /><figcaption>**Figure 6.** Examples of measurement/model comparison of vertical transects of the plume between 850–1300 m distance from the outfall (a, 7/10/2010), and 3200–4200 m from the outfall (b, 6/4/2011), adapted from [61]. Time and spatial scales (m) are the same for measurements and simulations. Insets: maps of surface sampling locations, adapted from [55].</figcaption></figure>
<figure><img src="image_10.png" /><figcaption>**Figure 7.** Variations with distance of the outfall of measured and simulated average dilution factors, adapted from [55].</figcaption></figure>

As the plume shape is more complex to the north of the La Hague Cape (Figure 5b), a similar model/measurement comparison is difficult to perform. Figure 8 shows the vertical distribution of concentrations 7 km downstream and five-six hours after a release. In this area, the tritium plume crossed the La Hague trough, of which the depths range from 80 to 100 m. These observations result from current shear between surface and bottom waters: the labeling is first transported in surface waters where the currents are larger.

<figure><img src="image_11.png" /><figcaption>**Figure 8.** Vertical concentrations of the plume five and six hours after a release and at 7 km distance.</figcaption></figure>

Left: distribution of surface sampling locations during the survey the 6/4/2011 (red dots).

Background: simulation of the expected plume at 12:15. Right: interpolation of in-depth measurements obtained during the survey at 11:58 and 13:10 the 6/4/11 (back and forth).

Numbers: concentrations measured at each in-depth location.

### 3.1.3. 2D Dispersion: Tidal Time Scale

As shown in Figure 7, the vertical homogenization of the plume is reached at 3 km downstream the outfall, which corresponds to one hour after the release. Except at particular locations where
strong bathymetry changes occurs (Figure 8), in the eastern English Channel, surface measurements of radionuclide concentrations are representative of the whole water column labeling. Investigations have been performed to test the model capability to simulate the dispersion of the La Hague plant within 1 to 48 h after the beginning of the release (DISPRO project). Beyond 48 h, it is no longer possible to distinguish two daily consecutive releases from another, because dilution and stirring mix them together. Within 48 h following a release, about 200 transects were done with sampling every 30 seconds (i.e., transects crossing a dispersion plume that could be traced back to a known release, and into which the width of the plume could be assessed). These transects represent roughly 3000 individual measurements, which encompass all tidal conditions (from neap to spring tides). An example of the tracking of one release together with a model/measurement comparison is given in Figure 9. Comparisons have been performed by accounting for the measured and simulated concentrations at a sampling location, and other criteria described in Section 2.5.1. Additional comparison criteria. Details of the results obtained are in [26]; Table 2 presents the main results of model/measurement comparisons for short scales model/measurements comparisons. For example, at short scales next to an outfall, the deviation between measured and simulated concentrations (expressed as mean dilution coefficient) is lower than 71% for a 95% confidence interval (Table 2).

<figure><img src="image_12.png" /><figcaption>**Figure 9.** Example of a comparison between the measured and simulated dispersion during 4 h after release, from 16:22 to 19:45 adapted from [25].</figcaption></figure>
<table><caption>**Table 2.** Characteristics of the model used to simulate soluble radioactive release dispersion in seawaters at a short scale next to an outfall around the Hague Cape.</caption><thead><tr><th colspan="2">Dispersion Characteristics 1h to 48 h after Release</th></tr></thead><tbody><tr><td>Geographical boundaries of the model</td><td>49°17′ N–49°55′ N; 2°26′ W–1°31′ W</td></tr><tr><td>Hydrodynamic model</td><td>Two-dimensional: Vertically-averaged velocities and concentrations, 110 m mesh size, 20 s time step</td></tr><tr><td>Average discrepancy between the mean dilution coefficients measured and simulated in the plumes</td><td>9% (–66% → 70%)</td></tr><tr><td>Average discrepancy per transect between simulated and measured maximum dilution coefficients</td><td>3% (–72% → 73%)</td></tr><tr><td>Average measured/simulated plume-width discrepancy</td><td>–6% (–73% → 65%)</td></tr><tr><td>Average discrepancy between measured and simulated plume-position, as a function of distance from the outlet point</td><td>–1% (–22% → 22%)</td></tr></tbody></table>

Figures in brackets indicate the 95% confidence interval.

## 3.2. Large Scale Model/Measurement Comparisons: Multi Tidal Time Scales

### 3.2.1. Individual Measurements

When a sufficient amount of measurements was obtained during oceanic campaigns, we could draw up radionuclide concentration maps at the scale of each survey. Campaigns performed repeatedly covered the English Channel and the North Sea between 1988 and 1996, the north-west Atlantic and Irish Sea between 1994 and 1996, the North Atlantic between 1997 and 2004 and the Bay of Biscay between 2009 and 2016. Hydrodynamic models accounting for all source terms and oceanic forcing (tide, meteorological forcing and large-scale circulation and hydrology) allow comparing of the measured and simulated concentrations at the date of each campaign. Figure 10 gives an example of such a comparison in 1994 for ¹²⁵Sb and 2016 for ³H.

<figure><img src="image_13.png" /><figcaption>**Figure 10.** (a) Measured and simulated ¹²⁵Sb in the English Channel, September 1994, adapted from [24]; measured (b) and simulated (c) surface ³H in the Bay of Biscay, spring 2016, adapted from [26].</figcaption></figure>
Computation of the measured and simulated concentrations at each sampling location and time gives statistics of the model representativeness. For example, in the English Channel and the North Sea, application of a residual Lagrangian model gives a mean difference between the 1400 individual values calculated and measured between 1988 and 1994 of 54% [24]. The deviation at a short scale next to an outfall is lower than 70% (Table 2). At the scale of the Biscay Bay with a 3D model the deviation is 21% [62].

### 3.2.2. Water Masses Labeling

In case of multi-tracer studies, it is possible to associate a specific labeling of the different water masses investigated and to draw up a "picture" of the different plumes associated to the different water masses. This means that a set of initial average concentrations is affected at each water mass entering the studied area. Then, at each sampling location in this area, the resolution of an equations system can determine the contribution of each water mass that fits the measured concentrations of tracers. This implies that the number of unknown (water mass to account) does not exceed the number of independent tracers measured.

This method was particularly fruitful in 1988 in the North Sea, where it was possible to distinguish and map the distribution of the four main water masses entering the North Sea by using the simultaneous measurements of three radionuclides ($^{125}$Sb, $^{137}$Cs, $^{134}$Cs) and salinity (Figure 11) [17]. As compared with another campaign performed two years before, the contribution of the direct fallout after the Chernobyl accident and the rate of renewing of North Sea waters were also calculated.

<figure><img src="image_14.png" /><figcaption>**Figure 11.** Percentage origins of water masses entering the North Sea in 1988 deduced from dissolved radionuclides measurements (a); $^{134}$Cs remaining fallout from the Chernobyl accident and rate of renewing of water masses in two years (b), adapted from [17].</figcaption></figure>

A similar approach was applied in the Bay of Biscay in 2016, by accounting for the different $^{3}$H labeling of waters coming from the Loire and Gironde rivers. Therefore, it is possible to compare
directly the extent of the plumes deduced from measurements and simulated separately as shown in
Figure 12 [40].

<figure><img src="image_15.png" /><figcaption>**Figure 12.** (a) Measured and (b) simulated freshwater contribution of the Loire River relatively to the Gironde River, adapted from [26] in spring 2016.</figcaption></figure>

3.2.3. Integration of Normalized Contributions

In a similar way as presented at the end of Section 2.5.1., it is possible to compare the plume distributions obtained from different campaigns by normalizing the measured concentrations with the corresponding fluxes of releases. This method allows the average dispersion characteristics and dilution coefficient of the considered source term to be mapped. It has been applied on a short scale in 3D (Figure 7), in the English Channel [20] and the North Sea [24], as presented in Figures 13 and 14. It shows that, on average, the Lagrangian residual models are able to properly catch the dispersion process at that large scale. This result was obtained after adapting the wind stress drag, which is one of the main hydrodynamic model calibration parameters.

<figure><img src="image_16.png" /><figcaption>**Figure 13.** Normalized measured and Lagrangian residual simulated distribution of $^{125}$Sb between 1983 and 1994 for a constant release of $10^6 \text{ Bq} \cdot \text{s}^{-1}$, adapted from [20].</figcaption></figure>
<figure><img src="image_17.png" /><figcaption>**Figure 14.** Normalized measured and simulated distribution of <sup>125</sup>Sb between 1988 and 1994 for a constant release of 10<sup>6</sup> Bq·s<sup>-1</sup>, adapted from [24].</figcaption></figure>

### 3.2.4. Inventories

As presented in Section 2.5.2., inventories of radionuclides quantities in different areas and water masses allow for the assimilation of numerous measurements in integrated quantities to account for the labeled depth and distribution of samples. The results could be compared with the known releases supposed to contribute to the labeling, and to model results in the same areas and water masses. Comparison with releases without simulation have been performed in the North Sea by slicing the southern North Sea into boxes where measured radionuclide quantities were compared with the La Hague releases [18] (Figure 15). From this correspondence it was possible to highlight the average transit time (one year from La Hague to the Skagerrak) and fluxes through the Dover Strait (97,000–195,000 m<sup>3</sup>·s<sup>-1</sup>), as well as estimations of losses for non-conservative radionuclides [18].

In the English Channel and the southern North Sea, the equilibrium between the measured and simulated quantities of radionuclides was an efficient tool to tune the wind friction that provide the best reliability of the Lagrangian residual model [24].

The loss of non-conservative radionuclides from seawater to sediments and living species (around 15% for <sup>134</sup>Cs, 75% for <sup>106</sup>Ru and 85% for <sup>60</sup>Co Table 3, Figure 16) was assessed from the measured quantity comparisons with the corresponding releases at the scale of the English Channel. It is then possible to check the environmental models accounting for the fluxes between the different environmental compartments. Such balanced budget exhibited an unknown source term at this scale, which represented twice the expected quantities measured for <sup>137</sup>Cs at the Channel scale (Table 3).

<figure><img src="image_18.png" /></figure>
<figure><img src="image_19.png" /><figcaption>**Figure 16.** Average impact of the La Hague reprocessing plant in seawater corresponding to a constant release of 1 MBq·s⁻¹; (a) ¹²⁵Sb, (b) ¹³⁴Cs, (c) ¹³⁷Cs, (d) ³H, (e) ¹⁰⁶Ru, (f) ⁶⁰Co, adapted from [20].</figcaption></figure>

<table><caption>**Table 3.** Comparison of measured quantities of radioactivity in the English Channel with correspondent releases from La Hague from 1983 to 1994, adapted from [20]. Antimony-125 is used as reference for the calculation; the background concentrations in Atlantic surface waters (¹³⁷Cs, ³H) and Chernobyl fallout (¹³⁷Cs, ¹³⁴Cs) has been deduced.</caption><thead><tr><th rowspan="3">Volume:<br>Equivalent Release Duration:</th><th colspan="5">Whole English Channel<br>4702 Km³<br>32 Weeks–7.3 Months</th><th colspan="5">Eastern English Channel<br>1576 Km³<br>25 Weeks–5.7 Months</th></tr><tr><th>¹³⁷Cs</th><th>¹³⁴Cs</th><th>¹⁰⁶Ru</th><th>¹²⁵Sb</th><th>⁶⁰Co</th><th>³H</th><th>¹³⁷Cs</th><th>¹³⁴Cs</th><th>¹⁰⁶Ru</th><th>¹²⁵Sb</th><th>⁶⁰Co</th><th>³H</th></tr><tr><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th></tr></thead><tbody><tr><th>Number of campaigns</th><td>3</td><td>2</td><td>3</td><td>3</td><td>3</td><td>1</td><td>5</td><td>3</td><td>5</td><td>5</td><td>2</td><td>1</td></tr><tr><th>Fraction of the La Hague<br>release</th><td>233%</td><td>86%</td><td>26%</td><td>98%</td><td>14%</td><td>103%</td><td>139%</td><td>83%</td><td>19%</td><td>98%</td><td>8%</td><td>121%</td></tr></tbody></table>

We hypothesized [18] that the ¹³⁷Cs excess in the English Channel resulted mainly from the influence of Sellafield releases. It was afterwards confirmed [11]; this evidences a pathway from the Irish Sea through the St George Channel in the Celtic Sea, and then a seasonal input in the English Channel as shown in Figure 1. The flux that reaches the English Channel represents around 1% of Sellafield ¹³⁷Cs. Figure 17 shows the English Channel areas that are the most impacted by the Sellafield releases. This corresponds to locations where sediments act as a delayed secondary source term for the water column, as demonstrated in the Irish Sea [63]. It concerns mostly the Hurd deep in the English Channel center and the French coastal areas.
<figure><img src="image_20.png" /><figcaption>**Figure 17.** Calculated from Figure 16 ($^{137}$Cs–$^{134}$Cs): $^{137}$Cs not coming from La Hague (average during 1986–1994).</figcaption></figure>

In the Celtic Sea, comparisons between inventories and releases from La Hague and Sellafield plants give an estimate of the fluxes and residence times of water masses in the Celtic Sea and North Est Atlantic approaches, with labeling from Sellafield releases along the western Ireland coasts [11] Figure 1.

At a global scale, a comparison of the measured ³H inventories and known inputs have been applied at the world oceans scale (Figure 18 [64]). It provides fluxes, estimations of residence time and values for the past and future concentration of ³H in Atlantic waters entering the European waters (78 Bq·m⁻³ in 2016).

<figure><img src="image_21.png" /><figcaption>**Figure 18.** Inventory of anthropogenic tritium by unit of surface area (kBq·m⁻²), adapted from [65].</figcaption></figure>

This value was applied to compare the measured and simulated radionuclide inventories in the Bay of Biscay. At this scale, the concentration simulated by the MARS3D model was compared to the measured ones by accounting for the variation of the ³H concentrations with depth [26]. It provides an estimation of the residence time (around one year for the continental shelf of the Bay of Biscay), and pathways of waters coming, respectively, from the Loire and Gironde rivers [26].
Estimation of radionuclides inventories was also applied during the accidental context
following the nuclear accident at Fukushima Dai-ichi power plant. Many seawater measurements of
¹³⁷Cs were provided continuously by Japanese authorities (Figure 19).

They enable the calculation of successive inventories of ¹³⁷Cs quantities in the surrounding Pacific waters and to estimate the marine source term coming at sea (Figure 20) [58]. This first estimation of 27 PBq (12 PBq–41 PBq) appeared to be very high at the time of publication, but later estimations were in the same order of magnitude [65–67]. The rapid decrease of ¹³⁷Cs quantities during the weeks following the accident reveals a very rapid environmental half-life that has been reproduced by MARS3D simulation (Figure 20 [27]).

<figure><img src="image_22.png" /><figcaption>Figure 19. Concentrations of ¹³⁷Cs in seawater between 11 April and 11 July 2011, adapted from [59].</figcaption></figure>
<figure><img src="image_23.png" /><figcaption>**Figure 20.** Time evolution of measured and simulated <sup>137</sup>Cs inventories present in the Fukushima area, adapted from [27].</figcaption></figure>

## 4. Discussion

This work gives accessibility to a comprehensive database, including in-situ radionuclide measurements performed during oceanographic campaigns and fluxes of dissolved radioactive release between 1982 and 2018 over the North-Western European continental shelf.

The results presented show applications of in-situ radionuclide measurements coupled with hydrodynamic models. They demonstrate how they improve the knowledge of dispersion in seawater. Applied alone, in situ data provide general pathways, transit time and dilution rate at different scales. Coupled with numerical hydrodynamic models, they improve and confirm their reliability with strong constraints. Physical measurements (as current, sea level or drifter measurements) generally represent short term parameters. Dissolved radionuclides complement physical data with information associated from short to long term transport and dispersion in seawater that could not be directly measured in currents.

Use of radiotracers for model validation was initiated by Salomon and Guéguéniat in 1988 with residual Lagrangian models [68]. This pioneering work was applied at different scales with space and time adaptations to obtain results comparable to measurements and simulations.

Artificial dissolved radionuclide tracers are particularly interesting because they fulfill a detection limit close to the background level, and contain few well known source terms and low risk of pollution during sampling and treatment. These advantages allow the detection of labeled plumes from short (hours, 100 m) to large scale (years, 1000 km). It counterbalances the efforts given to radionuclide measurements such as sampling large volumes for gamma emitters and measurement constraints.

### 4.1. Radiotracers

The use of artificial radionuclides such as oceanographic tracers depends from the controlled liquid releases performed by nuclear industry. These releases represent opportunities of existing labeling that could be used to follow water masses. They are issued mainly from European reprocessing plants. The fluxes of gamma emitters have decreased by two orders of magnitudes during the 1990s after application of a more efficient industrial purification process to diminish the environmental risk. It results in that the concentrations measured in the vicinity of the La Hague plant in 2018 are close to the detection limit and lower than the ones measured at 1500 km distance along the Norwegian coasts in 1988 (less than 2 Bq·m<sup>-3</sup>). Other radionuclides could still be applied as long term tracers, such as tritium [11,26,64], <sup>99</sup>Tc [69] <sup>129</sup>I [70,71] or <sup>236</sup>U [72–74].
³H is present in all nuclear plant releases and, in the form of tritiated water HTO, it could not be retained at the source by chemical process. At the nearest monitoring station to the main outfall in the marine system at La Hague Cape, concentrations in seawater are on average 12 Bq·m⁻³ (2017–2019, O. Connan pers. comm.). This concentration is three orders of magnitude lower from the OMS maximum recommended value of 10,000 Bq·L⁻¹ for drinking water. The ³H labeling from the nuclear plant will remain significant as long as the nuclear plants operate, and thus, they remain a relevant long term oceanographic tracer if we have the capability to measure low levels (from the oceanic background of 0.07 Bq·L⁻¹ [64]). The mapping of the dispersion of the Loire and Gironde waters at the scale of the Bay of Biscay demonstrates this approach [26]. Limitations for radiotracers measurement concern mainly the available means (low level counters, ³He degassing and mass spectrometry) and delay between sampling and measurement (from days to years). As an example, the ³He ingrowth method for ³H measurement exists in few laboratories around the world.

## 4.2. Hydrodynamic Models

Oceanographic hydrodynamic models have been continuously improved following the computing capabilities enhancement. It is possible to represent the whole continental shelf in 3D with a resolution of around one kilometer. This calculation efficiency must not mask the requested improvements concerning the knowledge and representation of physical process. The representativeness and precision of forcing the parameters determines the model reliability.

It concerns open boundaries limit conditions (stratification, currents, tides), and local data such as bathymetry precision, meteorological effects, bottom nature and associated friction. A better understanding of the surface and bottom drag coefficients remains an ongoing concern in research [75], as they require significant improvements.

A demonstration that a given hydrodynamic model is reliable does not prove that other similar models are at the same level. Long term advection and dispersion are particularly discriminant, as it represents only a few percent of the instantaneous currents. The radiotracer databases will remain useful to test the next hydrodynamic model generation.

Nevertheless, the level of realism reached between measurements and simulations give way to hydrodynamic model application at the scale of the European shelf or elsewhere in the world as shown in the north-western Pacific [27].

For short scale studies in 3D, MARS3D is not able to represent the real turbulence process (Figure 6). These phenomena are outside the initial scope of this model (i.e., regional scale); it must be improved by applying other numerical schemes or discard assumptions as the hydrostatic one. The computing capabilities remain a limit to simulate together short scales (meters) and longer ones (tens of kilometers); however, methods exist to overcome these limits (AGRIF) [76].

## 4.3. Scales for Model/Measurement Comparisons

This work shows that different approaches must be applied to measure the dispersion from short to medium or large scales. An operational definition to distinguish short scale from longer ones in tidal seas is the capability, or lack thereof, to distinguish each individual release plume. This depends from the diffusion process and the frequency of the main forcing parameter and releases. This results in that in the center of the English Channel, short scale studies concern fewer than four tidal cycles (two days).

The short scale requires models simulating rapid plume displacements (30 km in 3 h) with a high resolution and short time step (lower than 100 m and 20 s). A difference of 10 minutes or 200 m could change the punctual concentration close to the plume of one order of magnitude. If we consider that the model simulates the tidal propagation over thousands of kilometers before calculating the local tide, this implies strong constraints on the numerical scheme and forcing parameters (open boundaries conditions and drag coefficients). The dispersion of individual releases could be measured, but the rapid changes during dispersion do not allow drawing 3D or 2D maps that describe the plume dispersion several times per hour. 2D and one-dimensional (1D) comparisons were performed for 3D and 2D models comparisons, respectively.
The more the model/measurement comparison is done far from a source term, the more the plumes are smoothed and the lower the deviations between the punctual measurement and equivalent simulated values. Larger scales could use, when tidal effects are smoothed, instantaneous or Lagrangian residual models with a mesh size larger than one kilometer and time steps higher than 200 s. Exhaustive measurements could be performed that allow calculating inventories with possibility to compare measured, simulated and released radionuclides quantities.

Long term-large scale models obtain good model/measurement scores without under- or overestimation of radionuclides concentrations if they account well for the residual effects of the different influences (tidal, meteorological, frictions). These influences can be adjusted such as the wind drag coefficient obtained in [24]. The determination of detailed physical phenomena hidden behind these calibration parameters is an open research field.

## 4.4. *Sampling Tools*

The measurement of radionuclide dispersion led to the development of original systems to sample seawater in highly dynamic environments. High frequency samplers allow the simultaneous collection of 10 samples each 30 s while recording the time, location, volume, flow and depth [39]. The different parts of the developed tools may be used separately or together and adapted for other purposes.

Volumes collected could be extended or sampling depth increased that allow continuous in-depth sampling of any dissolved substances with a vessel sailing. 3D sampling is not limited to seawater and could concern larvae or particulate matter with continuous sampling and measurement of vertical suspended particulate matter (SPM) profiles.

This enhances the efficiency of oceanographic campaigns. As an example, the towed depressor Dynalest was applied as a support for instruments during sailing: an acoustic current doppler profiler (ADCP) was deployed to measure currents from 3 m depth up to the bottom without being perturbed by surface waves with a ship sailing normally in the Alderney race (0.5–4 m.s⁻¹).

Variability of plume dispersion made short scales studies impossible without accurate simulation tools to forecast its location. Sampling assisted by model simulation has demonstrated its efficiency in the vicinity of an outfall in a tidal environment (Figure 5 and Figure 9, [25,39]).

## 4.5. *Radionuclides Inventories*

By accounting for all available measurements associated with their representative area and depth, it is possible to map the radionuclide distributions and calculate radionuclides quantities existing in specific areas. Such inventories give integrated values with a better precision than individual concentrations [17,20].

Marine dissolved radionuclides give a rare example in environmental studies where balance budgets could be calculated between measured, simulated and releases inventories with precision better than 10% [25,26,20]. Comparison between radionuclides inventories quantities existing in seawater with the known sources results in determinations of environmental fluxes and transit times [18,26]. It highlights unknown sources and water masses pathways [11,20,58]. Inventories could also be compared with the simulation results and they represent robust check values to improve hydrodynamic models [20,24,27].

The balance obtained between released and simulated quantities of conservative radionuclides gives access to the loss of less conservative radionuclides. It is possible to quantify the fluxes between the different environmental compartments (i.e., seawater, living species and sediments) [18,20]. Radionuclide inventories are then powerful tools to check integrated environmental models accounting for stocks and fluxes of radionuclides at the scale of a continental sea as the English Channel.

## 4.6. *Normalization*
Normalization by a parameter representative of an environmental or anthropogenic forcing could be applied to long term time series or variable spatial distributions. This method made possible, by accounting for adequate time or space scale, to compare highly variable temporal or spatial changes of the environment as a result of variations depending on the season, meteorological effects, tides and releases. Normalization mitigates stochastics effects and highlight average characteristics of the environment difficult to quantify. It gives useful results with the measurements alone, but can also be compared with simulations. Others examples are: 3D short scale dispersion characteristics (Figure 7); average distributions and inventories in the English Channel (Figure 16 and Table 3) and North Sea (Figures 14 and 15); shape of the English Channel pathway in the southern North Sea [18]. Similar methods could be applied in areas where deterministic models fail to represent the water masses circulation as in the Kuroshio gyres in the Northern Pacific [77].

### 4.7. Perspectives

Tracers and methods applied here are potentially applicable elsewhere if the radionuclide labeling is sufficient to distinguish different water masses. This is obvious in the case in the Irish Sea, where the other main artificial source is the Sellafield reprocessing plant. The work performed in the Bay of Biscay demonstrates that nuclear power plant releases of tritium are sufficient to label water masses at a large scale if low level measurement methods are available [26]. This was also possible after precise evaluation of the North Atlantic surface water background concentration [64].

The possible application of ³H as water masses tracer potentially concern all seas reached by nuclear power plant releases. Maps of power plant locations with the corresponding tritium releases (Figure 18 [64,78]) suggest the Mediterranean Sea, North-West Atlantic, Arctic Sea, Western and Eastern Pacific, Black Sea, Baltic Sea, China Sea, Japan Sea, Gulf of Mexico and Arabian Sea.

Radionuclide dispersion models were applied to simulate realistic behavior of controlled industrial releases from 1984 up to 2016 in the English Channel (Figure 21). Such results could be compared with monitoring measurements and give a way to understand and interpret them. It was applied by IRSN to optimize the location of monitoring stations. Another application is the forecasting of accidental situations as it was performed after shipwrecks in the English Channel (Ievoli Sun in 2000, Ece in 2006) or after the Fukushima accident [27]. An operational tool is in course of development at IRSN to forecast the consequences of accidental situations (STERNE tool [79]).
<figure><img src="image_24.png" /><figcaption>**Figure 21.** MARS3D simulation of the dispersion of tritium industrial releases in seawater over the European continental shelf.</figcaption></figure>

Knowledge of dissolved transport (currents, advection and dispersion) is the basis of other oceanographic studies. It made possible pollutant transfer studies between seawater and living species [56,80–82] or particulate matter and physical transport of sediments [83–85]. Radio-ecological models accounting together for biological and hydro sedimentary process are accessible. The first tests have also shown the possibility to follow the sea/atmosphere exchanges by using ³H in water vapor [86].

Radionuclides behavior could be studied at the earth scale [65], as a tool for testing climate changes models (renewing of ocean water masses with ³H, water cycle, carbon cycle with ¹⁴C). Radionuclides were widely applied during the 1980s, but are not outdated environmental tracers.

Our wishes are that the provided database will contribute to these works and will be complemented.

**Supplementary Materials:** The following are available online at: Bailly du Bois, P. IRSN measurements of dissolved radioactivity in seawater, 1982–2016, 39642 stations. Available online: https://doi.pangaea.de/10.1594/PANGAEA.906541 (dataset in review, waiting for publication acceptation) (accessed on 02 October 2019). Bailly du Bois, P. IRSN measurements of dissolved radioactivity in seawater, 1982–2016, including liquid releases database issued from British and French nuclear plants. Available online: https://doi.pangaea.de/10.1594/PANGAEA.906749 (dataset in review, waiting for publication acceptation) (accessed on 2 October 2019).

**Author Contributions:** Conceptualization and Methodology, Formal analysis, P.B.d.B.; Data curation and databases, P.B.d.B., C.V. and P.-E.O.; Funding acquisition, P.B.d.B.; Investigation, P.B.d.B., F.D. and P.-E.O.; Project administration, P.B.d.B.; Software, M.M. and P.B.d.B.; Supervision, P.B.d.B.; Validation, P.B.d.B. and P.-E.O. Visualization, P.B.d.B. Writing—original draft preparation, P.B.d.B. Writing—review and editing, P.B.d.B. and F.D. Oceanographic campaigns, P.B.d.B., C.V., M.M., F.D. and P.-E.O. Radionuclide measurements,
L.S., P.B.d.B., C.L. and C.V. Model development, F.D. and P.B.d.B. Technical development, P.B.d.B. All authors have read and agreed to the published version of the manuscript.

**Funding:** Most of the oceanographic campaign uses vessels and crew provided by the French oceanographic fleet operated by CNRS/INSU (Centre National de la Recherche Scientifique, France), Ifremer-Genavir (Institut Français de Recherche pour l'Exploitation de la Mer, France) and GEA-Navy (Groupe d'études atomiques de la Marine nationale). British "Cirolana" and German "Validivia" oceanographic vessels participate also to the sampling collection. This research was partly funded by Orano Company for oceanographic data acquired between 1990 and 2011.

**Acknowledgments:** It is not possible to mention all people that have contributed to the collection of samples during oceanographic campaigns. In particular, we want to thank: the crews from the CNRS/INSU and Genavir vessels, Pierre Guéguéniat, Dominique Boust, Rémi Gandon, René Léon, Marianne Lamotte, Bruno Fiévet, Olivier Connan, Pascal Morin, Bernard Le Cann, Louis Marié, Robert Lafite, Kins Leonard, Peter Kershaw and Jürgen Hermann. The Orano La Hague SPR team must also be thanked for providing detailed liquid radioactive releases during and after oceanographic campaigns. Lastly, we thank the editors and anonymous reviewers for their fruitful comments during the review process.

**Conflicts of Interest:** The authors certify that they have no conflicts of interest to declare.
Appendix A

<table><caption>**Table A1.** List of oceanographic campaigns accounted in the database. Positions are referenced in WGS84 in decimal degrees (N, E: +; S, W:-). Radionuclide concentrations are in Bq·m⁻³ Abbreviations: LH: La Hague Cape; EC: English Channel; CI: Channel Islands; CS: Celtic Sea; IS: Irish Sea; NS: North Sea; MS: Mediterranean Sea; BB: Biscay Bay; NA: North Atlantic.</caption><thead><tr><th>Campaign name</th><th>Year</th><th>Area abr.</th><th>Stat. nb.Lon.</th><th>Min.Lon.</th><th>Max.Lat.</th><th>min.Lat. max.</th><th>H3</th><th>Cs137</th><th>Cs134</th><th>Sb125</th><th>Ru106</th><th>Co60</th></tr><tr><td colspan="14" style="text-align:center;">Nb. Max. meas. Nb. Max. meas. Nb. Max. meas. Nb. Max. meas. Nb. Max. meas. Nb. Max. meas.</td></tr></thead><tbody><tr><td>Cirolana0582</td><td>1982</td><td>EC, CS, NS</td><td>42</td><td>-7.0</td><td>2.0</td><td>49.5</td><td>58.7</td><td>39</td><td>642</td><td></td><td></td><td></td><td></td></tr><tr><td>Cirolana0583</td><td>1983</td><td>EC, CS, IS</td><td>54</td><td>-7.5</td><td>-1.9</td><td>49.3</td><td>54.7</td><td>44</td><td>13283</td><td>22</td><td>688</td><td>14</td><td>459</td><td>36</td><td>7511</td><td>1</td><td>21.8</td></tr><tr><td>Clione0983</td><td>1983</td><td>EC, CS</td><td>36</td><td>-8.5</td><td>1.3</td><td>48.6</td><td>51.6</td><td>35</td><td>37.0</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>GEA0983</td><td>1983</td><td>EC, CI</td><td>12</td><td>-2.5</td><td>-2.2</td><td>49.3</td><td>49.6</td><td>12</td><td>19.2</td><td>6</td><td>3.0</td><td>9</td><td>23.8</td><td>9</td><td>69</td><td></td><td></td></tr><tr><td>Thalassa1283</td><td>1983</td><td>EC, BB</td><td>56</td><td>-4.8</td><td>0.3</td><td>47.4</td><td>50.0</td><td>55</td><td>59.2</td><td>31</td><td>10.0</td><td>45</td><td>149</td><td>44</td><td>207</td><td>8</td><td>44.4</td></tr><tr><td>Pluteus0583</td><td>1983</td><td>EC, CI, NS</td><td>119</td><td>-2.7</td><td>2.9</td><td>48.9</td><td>51.3</td><td>115</td><td>69.2</td><td>21</td><td>15.5</td><td>105</td><td>1668</td><td>79</td><td>3256</td><td></td><td></td></tr><tr><td>Pluteus0686</td><td>1986</td><td>EC, CI</td><td>185</td><td>-4.8</td><td>2.0</td><td>48.3</td><td>51.3</td><td>184</td><td>34.5</td><td>141</td><td>6.7</td><td>163</td><td>211</td><td>134</td><td>686</td><td>29</td><td>13.0</td></tr><tr><td>Ferry1986</td><td>1986</td><td>EC, CS</td><td>31</td><td>-8.3</td><td>-4.0</td><td>48.8</td><td>51.8</td><td>31</td><td>15.4</td><td>7</td><td>5.3</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Pluteus1286</td><td>1986</td><td>EC</td><td>47</td><td>-1.8</td><td>-1.4</td><td>49.7</td><td>50.3</td><td>47</td><td>18.6</td><td>27</td><td>5.2</td><td>44</td><td>89.7</td><td>38</td><td>180</td><td>27</td><td>8.8</td></tr><tr><td>GEA1986</td><td>1986</td><td>EC, CS</td><td>97</td><td>-6.0</td><td>1.9</td><td>48.7</td><td>51.2</td><td>97</td><td>53.3</td><td>80</td><td>21.8</td><td>45</td><td>46.3</td><td>34</td><td>315</td><td>13</td><td>8.9</td></tr><tr><td>Nies0187</td><td>1987</td><td>NS</td><td>29</td><td>1.5</td><td>9.2</td><td>51.0</td><td>55.0</td><td>29</td><td>84.0</td><td>23</td><td>13.0</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Aurélia0387</td><td>1987</td><td>NS</td><td>31</td><td>3.0</td><td>4.7</td><td>51.7</td><td>53.0</td><td>31</td><td>15.5</td><td>31</td><td>5.7</td><td>30</td><td>57.1</td><td>29</td><td>213</td><td></td><td></td></tr><tr><td>Sépia1987</td><td>1987</td><td>EC, NS</td><td>14</td><td>1.3</td><td>2.3</td><td>50.3</td><td>51.2</td><td>14</td><td>11.2</td><td>4</td><td>1.8</td><td>14</td><td>48.5</td><td>14</td><td>121</td><td>1</td><td>0.4</td></tr><tr><td>Manche0687</td><td>1987</td><td>EC</td><td>50</td><td>-1.9</td><td>0.3</td><td>49.6</td><td>50.8</td><td>50</td><td>10.0</td><td>10</td><td>3.0</td><td>50</td><td>84.1</td><td>40</td><td>495</td><td>7</td><td>10.4</td></tr><tr><td>Luctor0687</td><td>1987</td><td>NS</td><td>63</td><td>2.0</td><td>3.9</td><td>51.1</td><td>51.9</td><td>63</td><td>31.2</td><td>19</td><td>4.0</td><td>61</td><td>62.2</td><td>19</td><td>120</td><td></td><td></td></tr><tr><td>Pluteus0887</td><td>1987</td><td>EC, CI, NS</td><td>165</td><td>-4.0</td><td>7.9</td><td>48.7</td><td>54.2</td><td>164</td><td>90.5</td><td>84</td><td>39.6</td><td>152</td><td>70.1</td><td>139</td><td>180</td><td>14</td><td>2.5</td></tr><tr><td>Aurélia1187</td><td>1987</td><td>NS</td><td>6</td><td>3.5</td><td>4.6</td><td>52.1</td><td>52.7</td><td>6</td><td>12.8</td><td>4</td><td>1.6</td><td>4</td><td>57.4</td><td>6</td><td>104</td><td></td><td></td></tr><tr><td>GEA1987</td><td>1987</td><td>EC</td><td>16</td><td>-2.0</td><td>2.5</td><td>49.5</td><td>51.2</td><td>16</td><td>11.8</td><td>11</td><td>2.6</td><td>16</td><td>44.8</td><td>13</td><td>143</td><td>5</td><td>2.2</td></tr><tr><td>Aurélia0488</td><td>1988</td><td>NS</td><td>9</td><td>3.5</td><td>4.6</td><td>52.1</td><td>52.9</td><td>9</td><td>12.6</td><td></td><td></td><td>6</td><td>58.0</td><td></td><td></td><td></td><td></td></tr><tr><td>Pluteus0488</td><td>1988</td><td>EC</td><td>44</td><td>-1.7</td><td>1.6</td><td>49.5</td><td>51.1</td><td>44</td><td>11.1</td><td>1</td><td>1.1</td><td>43</td><td>43.6</td><td>28</td><td>116</td><td></td><td></td></tr><tr><td>Pluteus0688</td><td>1988</td><td>EC</td><td>105</td><td>-5.1</td><td>1.6</td><td>49.4</td><td>51.1</td><td>104</td><td>12.2</td><td></td><td></td><td>81</td><td>112</td><td>63</td><td>222</td><td></td><td></td></tr><tr><td>Pluteus0788</td><td>1988</td><td>EC, CI</td><td>103</td><td>-4.0</td><td>-1.1</td><td>48.6</td><td>49.8</td><td>103</td><td>12.5</td><td>4</td><td>1.4</td><td>94</td><td>464</td><td>42</td><td>324</td><td>6</td><td>2.7</td></tr><tr><td>GEA0788</td><td>1988</td><td>EC</td><td>105</td><td>-5.0</td><td>1.6</td><td>49.4</td><td>51.1</td><td>104</td><td>12.2</td><td>21</td><td>1.9</td><td>81</td><td>180</td><td>64</td><td>222</td><td>5</td><td>2.6</td></tr><tr><td>Tramanor1</td><td>1988</td><td>EC, NS</td><td>256</td><td>-2.0</td><td>9.0</td><td>49.7</td><td>60.2</td><td>242</td><td>70.7</td><td>192</td><td>13.3</td><td>186</td><td>48.9</td><td>107</td><td>189</td><td>2</td><td>2.0</td></tr><tr><td>Gedynor</td><td>1989</td><td>EC, NS</td><td>127</td><td>-4.0</td><td>8.5</td><td>49.2</td><td>58.9</td><td>119</td><td>38.8</td><td>37</td><td>4.3</td><td>89</td><td>57.1</td><td>33</td><td>151</td><td>1</td><td>1.6</td></tr><tr><td>RadeBrest0689</td><td>1989</td><td>CI</td><td>23</td><td>-4.6</td><td>-4.3</td><td>48.3</td><td>48.4</td><td>23</td><td>5.9</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>RadeBrest0789</td><td>1989</td><td>CI</td><td>23</td><td>-4.6</td><td>-4.3</td><td>48.3</td><td>48.4</td><td>23</td><td>6.3</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>GEA0989</td><td>1989</td><td>EC, CI</td><td>88</td><td>-2.6</td><td>-0.9</td><td>49.4</td><td>50.2</td><td>87</td><td>31.5</td><td>57</td><td>7.8</td><td>88</td><td>107.4</td><td>80</td><td>355</td><td>23</td><td>15.9</td></tr><tr><td>Pluteus0989</td><td>1989</td><td>EC, CI</td><td>124</td><td>-2.8</td><td>-1.0</td><td>49.2</td><td>50.3</td><td>124</td><td>37.0</td><td>71</td><td>4.8</td><td>123</td><td>203.0</td><td>100</td><td>985</td><td>53</td><td>13.9</td></tr></tbody></table>
<table><tbody><tr><td>Disver2008</td><td>2008</td><td>LH, EC, CI</td><td>796</td><td>-2.0</td><td>-2.0</td><td>49.7</td><td>49.7</td><td>784</td><td>190660</td></tr><tr><td>Aspex2009</td><td>2009</td><td>BB</td><td>36</td><td>-6.0</td><td>-1.5</td><td>44.0</td><td>47.8</td><td>34</td><td>336</td></tr><tr><td>Ovide2010</td><td>2010</td><td>NA</td><td>1</td><td>-19.1</td><td>-19.1</td><td>45.8</td><td>45.8</td><td>1</td><td>76</td></tr><tr><td>Aspex2010</td><td>2010</td><td>BB</td><td>62</td><td>-6.4</td><td>-1.6</td><td>44.0</td><td>48.3</td><td>62</td><td>266</td></tr><tr><td>Disver2010</td><td>2010</td><td>LH</td><td>5556</td><td>-2.0</td><td>-2.0</td><td>49.7</td><td>49.7</td><td>5314</td><td>7029776</td></tr><tr><td>Disver2011</td><td>2011</td><td>LH</td><td>12498</td><td>-2.0</td><td>-1.9</td><td>49.7</td><td>49.8</td><td>12498</td><td>8739130</td></tr><tr><td>Aspex2011</td><td>2011</td><td>BB</td><td>74</td><td>-6.0</td><td>-1.5</td><td>44.0</td><td>47.8</td><td>62</td><td>211</td></tr><tr><td>Trimadu2013</td><td>2013</td><td>NA</td><td>18</td><td>-13.2</td><td>55.5</td><td>-32.4</td><td>34.9</td><td>17</td><td>80</td><td>12</td><td>1.1</td></tr><tr><td>Traces2014</td><td>2014</td><td>EC, CI</td><td>1205</td><td>-3.0</td><td>-1.6</td><td>48.6</td><td>49.8</td><td>1170</td><td>125159</td></tr><tr><td>Traces2015</td><td>2015</td><td>EC, CI</td><td>547</td><td>-3.0</td><td>-1.7</td><td>48.8</td><td>49.8</td><td>546</td><td>38320</td></tr><tr><td>Dynsedim2016</td><td>2016</td><td>BB</td><td>31</td><td>-5.3</td><td>-2.6</td><td>45.6</td><td>47.9</td><td>30</td><td>640</td></tr><tr><td>Pelgas2016</td><td>2016</td><td>BB</td><td>130</td><td>-5.8</td><td>-1.3</td><td>43.7</td><td>47.9</td><td>125</td><td>583</td></tr><tr><td>Plume2016</td><td>2016</td><td>BB</td><td>254</td><td>-5.2</td><td>-1.1</td><td>44.7</td><td>48.6</td><td>177</td><td>19072</td><td>12</td><td>1.3</td></tr><tr><td>Goury (time series)</td><td>1984-2018</td><td>LH, EC</td><td>744</td><td>-1.9</td><td>-1.9</td><td>49.7</td><td>49.7</td><td>266</td><td>41999</td><td>615</td><td>79.9</td><td>260</td><td>62.3</td><td>498</td><td>223.1</td><td>484</td><td>1745</td><td>398</td><td>50.7</td></tr><tr><td>Nb. total (campaigns):</td><td></td><td></td><td>39642</td><td></td><td></td><td></td><td></td><td>35663</td><td></td><td>3492</td><td></td><td>1295</td><td></td><td>2242</td><td></td><td>1606</td><td></td><td>568</td><td></td></tr><tr><td>Nb. Total (all):</td><td>80</td><td></td><td>40386</td><td></td><td></td><td></td><td></td><td>35929</td><td></td><td>4107</td><td></td><td>1555</td><td></td><td>2740</td><td></td><td>2090</td><td></td><td>966</td><td></td></tr></tbody></table>
References

1. Guary, J.C.; Guéguéniat, P.; Pentreath, J. *Radionucléides : A Tool for Oceanography*; Elsevier Applied Science: Amsterdam, The Netherlands, 1988; p. 461.

2. Kershaw, J.P.; Woodhead, D.S. *Radionuclides in the Study of Marine Processes*; Elsevier Applied Science: Amsterdam, The Netherlands, 1991; p. 393.

3. Germain, P.; Guary, J.C.; Guéguéniat, P.; Métivier, H.E. Radionuclides in the Oceans. In *Proceedings of the RADOC 96–97, Part 1, Inventories, Behaviour and Processes, Cherbourg-Octeville, France, 7–11 October 1996; Radioprotection numéro spécial, colloques Volume 32, C2; Les Editions de Physique: Les Ulis, France 1997*; p. 422.

4. Herrmann, J.; Kershaw, P.; Bailly du Bois, P.; Guéguéniat, P. The distribution of artificial radionuclides in the English Channel, southern North Sea, Skagerrak and Kattegat, 1990–1993. *J. Mar. Syst.* **1995**, *6*, 427–456, doi:10.1016/0924-7963(95)00026-1.

5. Povinec, P.P.; Bailly du Bois, P.; Kershaw, P.; Nies, H.; Scotto, P. Temporal and spatial trends in the distribution of 137Cs in surface waters of Northern European Seas—A record of 40 years of investigations. *Deep. Sea Res. Part II Top. Stud. Oceanogr.* **2003**, *50*, 2785–2801, doi:10.1016/s0967-0645(03)00148-6.

6. Kershaw, P.; McCubbin, D.; Leonard, K. Continuing contamination of north Atlantic and Arctic waters by sellafield radionuclides. *Sci. Total. Environ.* **1999**, *237*, 1119–1132.

7. Kautsky, H. Investigations on the distribution of 137Cs, 134Cs and 90Sr and the water mass transport times in the Northern North Atlantic and the North Sea. *Ocean Dyn.* **1987**, *40*, 49–69.

8. Dahlgaard, H. Transfer of European coastal pollution to the arctic: Radioactive tracers. *Mar. Pollut. Bull.* **1995**, *31*, 33–37.

9. Mitchell, P.I.; Holm, E.; Dahlgaard, H.; Boust, D.; Leonard, K.S.; Papucci, C.; Salbu, B.; Christensen, G.; Strand, P.; Sánchez-Cabeza, J.A.; et al. *Radioecological Assessment of the Consequences of Contamination of Arctic Waters: Modelling the Key Processes Controlling Radionuclide Behaviour Under Extreme Conditions (ARMARA)*; EC Nuclear Fission Safety Programme 1995-99 Contract No. F14P-CT95-0035 Final Report; Department of Experimental Physics University College Dublin: Dublin, Ireland, 1999; p. 151.

10. Nies, H.; Albrecht, H.; Herrmann, J. In *Radionuclides in Water and Suspended Particulate Matter from the North Sea, Proceedings of the Radionuclides in the Study of Marine Processes, Norwich, UK, 10–13 September 1991*; Kershaw, J.P., Woodhead, D.S., Eds.; Elsevier Applied Science: Amsterdam, The Netherlands, 1991; pp. 24–36.

11. Bailly du Bois, P.; Germain, P.; Rozet, M.; Solier, L. In *Water masses circulation and residence time in the celtic sea and English channel approaches, characterisation based on radionuclides labelling from industrial releases, Proceedings of the International Conference on Radioactivity in Environment, Monaco, Principality of Monaco, 1–5 September 2002*; Borretzen, P., Jolle, T. Strand, P., Eds.; Norwegian Radiation Protection Authority: Østerås, Norway; pp. 395–399.

12. Guéguéniat, P.; Salomon, J.; Wartel, M.; Cabioch, L.; Fraizier, A. Transfer pathways and transit time of dissolved matter in the eastern English channel indicated by Space-Time radiotracers measurement and hydrodynamic modelling. *Estuarine Coast. Shelf Sci.* **1993**, *36*, 477–494, doi:10.1006/ecss.1993.1029.

13. Guéguéniat, P.; Bailly du Bois, P.; Gandon, R.; Salomon, J.; Baron, Y.; Leon, R. Spatial and temporal distribution (1987–91) of 125Sb used to trace pathways and transit times of waters entering the north sea from the English channel. *Estuarine Coast. Shelf Sci.* **1994**, *39*, 59–74, doi:10.1006/ecss.1994.1049.

14. Guéguéniat, P.; Bailly du Bois, P.; Salomon, J.C. Transfer pathways and transit time of dissolved matter in the english channel, north sea and arctic indicated by Space-Time radiotracers measurement. *Estuarine, Coast. Shelf Sci.* **1993**, *36*, 477–494.

15. Guéguéniat, P.; Herrmann, J.; Kershaw, P.; Bailly du Bois, P.; Baron, Y. Artificial radioactivity in the English Channel and the North Sea. In *Radionuclides in the Oceans, Inputs and Inventories, RADOC 96–97; Guéguéniat, P., Germain, P., Métivier, H., Eds.*; Les Editions de Physique: Les Ulis, France, 1997; pp. 121–154.

16. Guéguéniat, P.; Kershaw, P.; Hermann, J.; Bailly du Bois, P. New estimation of La Hague contribution to the artificial radioactivity of Norwegian waters (1992–1995) and Barents Sea (1992–1997). *Sci. Total. Environ.* **1997**, *202*, 249–266, doi:10.1016/s0048-9697(97)00120-4.
17. Bailly du Bois, P.; Guéguéniat, P.; Gandon, R.; Leon, R.; Baron, Y. Percentage contribution of inputs from the Atlantic, Irish sea, English channel and baltic into the north sea during 1988: A Tracer-Based evaluation using artificial radionuclides. *Neth. J. Sea Res.* **1993**, *31*, 1–17, doi:10.1016/0077-7579(93)90012-h.

18. Bailly du Bois, P.; Salomon, J.; Gandon, R.; Guéguéniat, P. A quantitative estimate of English channel water fluxes into the north sea from 1987 to 1992 based on radiotracer distribution. *J. Mar. Syst.* **1995**, *6*, 457–481, doi:10.1016/0924-7963(95)00018-k.

19. Bailly du Bois, P.; Rozet, M.; Thoral, K.; Salomon, J.C. In *Improving knowledge of Water-Mass circulation in the english channel using radioactive tracers, Proceedings of the Radioprotection-Colloques, April 1997, Numéro spécial "Radionuclides in the Oceans", RADOC 96-97, Proceedings Part 1 "Inventories, Behaviour and Processes", Cherbourg-Octeville, France, 7-11 October 1996*; Germain, P., Guary, J.C., Guéguéniat, P., Métivier, H., Eds.; Les Editions de Physique: Les Ulis, France; 1997, Volume 32, pp. 63–69.

20. Bailly du Bois, P.; Guéguéniat, P. Quantitative assessment of dissolved radiotracers in the English channel: Sources, average impact of la Hague reprocessing plant and conservative behaviour (1983, 1986, 1988, 1994). *Cont. Shelf Res.* **1999**, *19*, 1977–2002, doi:10.1016/s0278-4343(99)00049-7.

21. Guéguéniat, P.; Bailly du Bois, P.; Salomon, J.; Masson, M.; Cabioch, L. Fluxmanche radiotracers measurements: A contribution to the dynamics of the English channel and north sea. *J. Mar. Syst.* **1995**, *6*, 483–494, doi:10.1016/0924-7963(95)00019-l.

22. Salomon, J.C.; Guegueniat, P.; Breton, M. In *Mathematical model of 125Sb transport and dispersion in the channel, Proceedings of the Radionuclides in the Study of Marine Processes, Norwich, UK, 10–13 September 1991*; Kershaw J.P., Woodhead D.S., Eds.; Elsevier Applied Science: London, UK, 1991; pp. 74–83.

23. Salomon, J.C.; Breton, M.; Fraizier, A.; Bailly du Bois, P.; Guéguéniat, P. A Semi-Analytic Mathematical Model for Dissolved Radionuclides Dispersion in the Channel Isles Region. *Radioprotect. Colloq.* **1997**, *32*, 375–380.

24. Bailly du Bois, P.; Dumas, F. Fast hydrodynamic model for medium-and Long-Term dispersion in seawater in the English channel and southern north sea, qualitative and quantitative validation by radionuclide tracers. *Ocean Model.* **2005**, *9*, 169–210, doi:10.1016/j.ocemod.2004.07.004.

25. Bailly du Bois, P.; Dumas, F.; Solier, L.; Voiseux, C. In-Situ database toolbox for Short-Term dispersion model validation in Macro-Tidal seas, application for 2D-Model. *Cont. Shelf Res.* **2012**, *36*, 63–82, doi:10.1016/j.csr.2012.01.011

26. Oms, P.-E. Transferts Multi-Échelles des Apports Continentaux dans le golfe de gascogne. Ph.D. Thesis, Université de Bretagne Occidentale: Brest, France, June 2019.

27. Bailly du Bois, P.; Garreau, P.; Laguionie, P.; Korsakissok, I. Comparison between modelling and measurement of marine dispersion, environmental Half-Time and 137Cs inventories after the Fukushima Daiichi accident. *Ocean Dyn.* **2014**, *64*, 361–383, doi:10.1007/s10236-013-0682-5.

28. Garreau, P.; Bailly du Bois, P. Transportation of radionuclides in celtic sea a possible mechanisms. *Radioprotect. Colloq.* **1997**, *32*, 381–385.

29. Perianez, R. Modelling the tidal dispersion of 137Cs and 239240Pu in the English channel. *J. Environ. Radioact.* **2000**, *49*, 259–277.

30. IAEA. MARIS. Availabe online: https://maris.iaea.org/Home.aspx (accessed on 23 September 2019).

31. Boyer, T.P.; Antonov, J.I.; Baranova, O.K.; Garcia, H.; Grodsky, A.; Johnson, D.; Locarnini, R.A.; Mishonov, A.; O'Brien, T.; Paver, C., et al. *World Ocean Database 2013*; NOAA: Silver Spring, MD, USA, 2013; p. 15.

32. NOAA. World Ocean Database. Availabe online: https://www.nodc.noaa.gov/OC5/WOD/pr_wod.html (accessed on 23 September 2019).

33. Ifremer. SISMER. Availabe online: https://data.ifremer.fr/SISMER (accessed on 23 September 2019).

34. British Oceanographic Data Centre (BODC). Availabe online: https://www.bodc.ac.uk/ (accessed on 23 September 2019).

35. Gandon, R.; Bailly du Bois, P.; Baron, E.Y. Caractère conservatif de l'antimoine 125 dans les eaux marines soumises à l'influence des rejets de l'usine de retraitement des combustibles irradiés de La Hague. *Radioprotection* **1998**, *33*, 457–482, doi:10.1051/radiopro:1998116.

36. Kautsky, H. Artificial radioactivity in the North Sea and the northern north Atlantic during the years 1977 to 1986; INIS-MF-11049; International Atomic Energy Agency (IAEA): Hamburg, Germany, 1988, 471–102.

37. Gandon, R.; Guéguéniat, P. Preconcentration 125Sb of onto MnO2 from Seawater Samples for Gamma-Ray Spectrometric Analysis. *Radiochim. Acta* **1992**, *57*, 159–164, doi:10.1524/ract.1992.57.23.159.
38. Gandon, R.; Boust, D.; Bedue, O. Ruthenium complexes originating from the purex process: Coprecipitation with copper ferrocyanides via ruthenocyanide formation. *Radiochim. Acta* **1993**, *61*, 41–45.

39. Bailly du Bois, P.; Pouderoux, B.; Dumas, F. System for High-Frequency simultaneous water sampling at several depths during sailing. *Ocean Eng.* **2014**, *91*, 281–289, doi:10.1016/j.oceaneng.2014.09.022.

40. McCubbin, D.; Leonard, K.S.; A.; Bailey, T.; Williams, J.; Tossell, P. Incorporation of organic tritium (3H) by marine organisms and sediment in the severn estuary/Bristol channel (UK). *Mar. Pollut. Bull.* **2001**, *42*, 2852–2863.

41. Ostlund, H.G.; Dorsey, H.G. *Rapid Electrolytic Enrichment and Hydrogen Gas Proportional Counting of Tritium*; Slovenske Pedagogicke Nakladatelstvo: Bratislava, Czechoslovakia, 1977.

42. Clarke, W.; Jenkins, W.; Top, Z. Determination of tritium by mass spectrometric measurement of 3He. *Int. J. Appl. Radiat. Isot.* **1976**, *27*, 515–522, doi:10.1016/0020-708x(76)90082-x.

43. Bailly du Bois, P. IRSN Measurements of Dissolved Radioactivity in Seawater, 1982–2016, 39642 Stations. PANGAEA Data Publisher: Bremerhaven, Germany, 2011. Availabe online: https://doi.pangaea.de/10.1594/PANGAEA.906541 (accessed on 2 October 2019).

44. Bailly du Bois, P.; Dumas, F.; Solier, L.; Voiseux, C. Tritium sampled and measured in surface water along cruise tracks. PANGAEA Data Publisher: Bremerhaven, Germany, 2011. Availabe online: https://doi.pangaea.de/10.1594/PANGAEA.762261 (accessed on 02 October 2019).

45. Bailly du Bois, P.; Dumas, F.; Solier, L.; Voiseux, C. Controlled tritium liquid releases from Areva-NC reprocessing plant. PANGAEA Data Publisher: Bremerhaven, Germany, 2011. Availabe online: https://doi.pangaea.de/10.1594/PANGAEA.762428 (accessed on 02 October 2019).

46. Gray, J.; Jones, S.R.; Smith, A.D. Discharges to the environment from the Sellafield site, 1951–1992. *J. Radiol. Prot.* **1995**, *15*, 99–131, doi:10.1088/0952-4746/15/2/001.

47. RIFE. RIFE-23 Radioactivity in Food and the Environment, 2017; Environment Agency, Food Standards Agency, Food Standards Scotland, Natural Resources Wales, Northern Ireland Environment Agency, Scottish Environment Protection Agency Preston: Lancashire, UK, 2018; p. 260.

48. Sellafield, L. *Monitoring Our Environment, Discharges and Environmental Monitoring Annual Report 2017*; Sellafield Ltd.: Sellafield Site, Seascale, Cumbria, UK, 2018; p. 260.

49. Bailly du Bois, P. IRSN Measurements of Dissolved Radioactivity in Seawater, 1982–2016, Including Liquid Releases Database Issued from British and French Nuclear Plants. PANGAEA Data Publisher: Bremerhaven, Germany, 2011. Availabe online: https://doi.pangaea.de/10.1594/PANGAEA.906749 (accessed on 2 October 2019).

50. Lazure, P.; Dumas, F. An External–Internal mode coupling for a 3D hydrodynamical model for applications at regional scale (MARS). *Adv. Water Resour.* **2008**, *31*, 233–250, doi:10.1016/j.advwatres.2007.06.010.

51. Bailly du Bois, P. Automatic calculation of bathymetry for coastal hydrodynamic models. *Comput. Geosci.* **2011**, *37*, 1303–1310, doi:10.1016/j.cageo.2010.11.018.

52. Salomon, J.C.; Breton, M. An atlas of Long-Term currents in the Channel. *Oceanol. Acta* **1993**, *16*, 439–448.

53. Salomon, J.; Breton, M.; Guéguéniat, P. A 2D long term Advection–Dispersion model for the Channel and southern North Sea Part B: Transit time and transfer function from Cap de La Hague. *J. Mar. Syst.* **1995**, *6*, 515–527, doi:10.1016/0924-7963(95)00021-g.

54. Orbi, A.; Salomon, J.C. Dynamique de marée dans le Golfe Normand-Breton. *Oceanol. Acta* **1988**, *11* 55–64.

55. Bailly du Bois, P.; Morillon, M.; Solier, L. *Mesure et Modélisation de la Dispersion Verticale dans le raz Blanchard (projet DISVER)*; IRSN/Prp-Env/SERIS2014-21; IRSN: Cadarache, France, 2014; p. 56.

56. Fiévet, B.; Bailly du Bois, P.; Laguionie, P.; Morillon, M.; Arnaud, M.; Cunin, P. A dual pathways transfer model to account for changes in the radioactive caesium level in demersal and pelagic fish after the Fukushima Dai-ichi nuclear power plant accident. *PLoS ONE* **2017**, *12*, e0172442, doi:10.1371/journal.pone.0172442.

57. Matheron, G. Les Concepts de base et L'Evolution de la geostatistique miniere. *Adv. Geostat. Min. Ind.* **1976**, *3–10*, doi:10.1007/978-94-010-1470-0_1.

58. Bailly du Bois, P.; Laguionie, P.; Boust, D.; Korsakissok, I.; Didier, D.; Fiévet, B. Estimation of marine Source-Term following Fukushima Dai-Ichi accident. *J. Environ. Radioact.* **2012**, *114*, 2–9, doi:10.1016/j.jenvrad.2011.11.015.
59. Bailly du Bois, P.; Dumas, F.; Morillon, M.; Furgerot, L.; Voiseux, C.; Poizot, E.; Méar, Y.; Bennis, A.C. The Alderney Race: General hydrodynamic and particular features. Philosophical Transactions of the Royal Society of London A, accepted March 2020. (in press).

60. Bailly du Bois, P.; Dumas, F. Dissolved radionuclide Measurements Used for Qualitative and Quantitative Calibration of Hydrodynamic Models in the English Channel and the North Sea; Validation of "TRANSMER" Model. In Proceedings of the 34th International Liege Colloquium on Ocean Hydrodynamics, Tracer Methods in Geophysical Fluid Dynamics, Liege, Belgium, 6–10 May 2002; p. 7.

61. Bailly du Bois, P. *Dispersion des Radionucléides dans les mers du Nord-Ouest de l'Europe : Observations et Modélisation*, Collection HDR ed.; IRSN, Ed.; IRSN: BP 17-92262 Fontenay -aux-Roses Cedex, France, 2014; p. 329.

62. Oms, P.E.; Bailly du Bois, P.; Dumas, F.; Lazure, P.; Morillon, M.; Voiseux, C.; Le Corre, C.; Solier, L.; Maire, D.; Caillaud, M. Tritium as an original continental runoffs tracer in the bay of biscay : Measurements and modelling. In Proceedings of the ISOBAY, Anglet, France, 7 June 2018.

63. Hunt, J.; Leonard, K.; Hughes, L. Artificial radionuclides in the Irish Sea from sellafield: Remobilisation revisited. *J. Radiol. Prot.* **2013**, *33*, 261–279, doi:10.1088/0952-4746/33/2/261.

64. Oms, P.-E.; Bailly du Bois, P.; Dumas, F.; Lazure, P.; Morillon, M.; Voiseux, C.; Le Corre, C.; Cossonnet, C.; Solier, L.; Morin, P. Inventory and distribution of tritium in the oceans in 2016. *Sci. Total. Environ.* **2019**, *656*, 1289–1303, doi:10.1016/j.scitotenv.2018.11.448.

65. Charette, M.A.; Breier, C.F.; Henderson, P.B.; Pike, S.M.; Rypina, I.I.; Jayne, S.R.; Buesseler, K.O. Radium-Based estimates of cesium isotope transport and total direct ocean discharges from the fukushima nuclear power plant accident. *Biogeosciences* **2013**, *10*, 2159–2167, doi:10.5194/bg-10-2159-2013.

66. Buesseler, K. Fukushima and Ocean Radioactivity. *Oceanography* **2014**, *27*, 92–105.

67. Aoyama, M.; Hamajima, Y.; Inomata, Y.; Kumamoto, Y.; Oka, E.; Tsubono, T.; Tsumune, D. Radiocaesium derived from the TEPCO Fukushima accident in the North Pacific Ocean: Surface transport processes until 2017. *J. Environ. Radioact.* **2018**, *189*, 93–102, doi:10.1016/j.jenvrad.2018.03.014.

68. Salomon, J.C.; Guéguéniat, P.; Orbi, A.; Baron, Y. In *A lagrangian model for long term tidally induced transport and mixing, Verification by artificial radionuclide concentrations, Proceedings of the Radionucléides: A Tool for Oceanography, Cherbourg, France, 1–5 June 1987*; Guary, J.C., Guéguéniat, P., Pentreath, R.J., Eds. Elsevier Applied Science: London, UK, 1988, pp. 384–394.

69. Simonsen, M.; Saetra, Ø.; Isachsen, P.E.; Lind, O.C.; Skjerdal, H.K.; Salbu, B.; Heldal, H.E.; Gwynn, J.P. The impact of tidal and mesoscale eddy advection on the long term dispersion of 99 Tc from Sellafield. *J. Environ. Radioact.* **2017**, *177*, 100–112, doi:10.1016/j.jenvrad.2017.06.002.

70. Alfimov, V.; Aldahan, A.; Possnert, G. Water masses and 129I distribution in the Nordic Seas. *Nucl. Instruments Methods Phys. Res. Sect. Beam Interactions Mater. Atoms* **2013**, *294*, 542–546, doi:10.1016/j.nimb.2012.07.042.

71. Daraoui, A.; Tosch, L.; Michel, R.; Goroncy, I.; Nies, H.; Synal, H.-A.; Alfimov, V.; Walther, C.; Gorny, M.; Herrmann, J. Iodine-129, Iodine-127 and Cesium-137 in seawater from the North Sea and the Baltic Sea. *J. Environ. Radioact.* **2016**, *162*, 289–299, doi:10.1016/j.jenvrad.2016.06.006.

72. Christl, M.; Casacuberta, N.; Lachner, J.; Herrmann, J.; Synal, H.-A. Anthropogenic 236U in the North Sea–A Closer Look into a Source Region. *Environ. Sci. Technol.* **2017**, *51*, 12146–12153, doi:10.1021/acs.est.7b03168.

73. Castrillejo, M.; Casacuberta, N.; Christl, M.; Vockenhuber, C.; Synal, H.-A.; García-Ibáñez, M.I.; Lherminier, P.; Sarthou, G.; Garcia-Orellana, J.; El Zrelli, R. Tracing water masses with 129I and 236U in the subpolar North Atlantic along the GEOTRACES GA01 section. *Biogeosci. Discuss.* **2018**, *2018*, 1–28, doi:10.5194/bg-2018-228.

74. Christl, M.; Casacuberta, N.; Vockenhuber, C.; Elsässer, C.; Bailly du Bois, P.; Herrmann, J.; Synal, H.-A. Reconstruction of the 236U input function for the northeast Atlantic Ocean—Implications for 129I/236U and 236U/238U-based tracer ages. *J. Geophys. Res. Ocean.* **2015**, doi:10.1002/2015JC011116.

75. Bennis, A.C.; Furgerot, L.; Bailly du Bois, P.; Dumas, F.; Odaka, T.; Lathuilière, C.; Filipot, J.F. Numerical modelling of Three-Dimensional Wave-Current interactions in extreme hydrodynamic conditions: Application to Alderney Race. *Appl. Ocean. Res.* **2019**, *95*, 102021.

76. Debreu, L.; Vouland, C.; Blayo, E. AGRIF: Adaptive grid refinement in Fortran. *Comput. Geosci.* **2008**, *34*, 8–13, doi:10.1016/j.cageo.2007.01.009.
77. Science-Council-of-Japan. *A Review of the Model Comparison of Transportation and Deposition of Radioactive Materials Released to the Environment as a Result of the Tokyo Electric Power Company's Fukushima Daiichi Nuclear Power Plant Accident*; Science Council of Japan: Tokyo, Japan, 2014; p. 111.

78. Usman, S.M.; Jocher, G.; Dye, S.; McDonough, W.F.; Learned, J. AGM2015: Antineutrino Global Map 2015. *Sci. Rep.* **2015**, *5*, 13945, doi:10.1038/srep13945.

79. Duffa, C.; Bailly du Bois, P.; Caillaud, M.; Charmasson, S.; Couvez, C.; Didier, D.; Dumas, F.; Fiévet, B.; Morillon, M.; Renaud, P.; et al. Development of emergency response tools for accidental radiological contamination of French coastal areas. *J. Environ. Radioact.* **2016**, *151*, 487–494, doi:10.1016/j.jenvrad.2015.04.019.

80. Fiévet, B.; Plet, D. Estimating biological Half-Lifes of radionuclides in marine compartments from environmental time-series measurements. *J. Environ. Radioact.* **2003**, *65*, 91–107, doi:10.1016/s0265-931x(02)00090-5.

81. Fiévet, B.; Voiseux, C.; Rozet, M.; Masson, M.; Bailly du Bois, P. Transfer of radiocarbon liquid releases from the AREVA La Hague spent fuel reprocessing plant in the English Channel. *J. Environ. Radioact.* **2006**, *90*, 173–196, doi:10.1016/j.jenvrad.2006.06.014.

82. Fiévet, B.; Pommier, J.; Voiseux, C.; Bailly du Bois, P.; Laguionie, P.; Cossonnet, C.; Solier, L. Transfer of tritium released into the marine environment by french nuclear facilities bordering the English channel. *Environ. Sci. Technol.* **2013**, *47*, 6696–6703, doi:10.1021/es400896t.

83. Blanpain, O.; Bailly du Bois, P.; Cugier, P.; Lafite, R.; Lunven, M.; Dupont, J.; Le Gall, E.; Legrand, J.; Pichavant, P. Dynamic sediment profile imaging (DySPI): A new field method for the study of dynamic processes at the sediment-water interface. *Limnol. Oceanogr. Methods* **2009**, *7*, 8–20, doi:10.4319/lom.2009.7.8.

84. Blanpain, O. Dynamique Sédimentaire Multiclasse: De L'étude des Processus À La Modélisation en Manche. PhD Thesis, Rouen University, Mont-Saint-Aignan, France, October 2009.

85. Rivier, A.; Le Hir, P.; Bailly du Bois, P.; Laguionie, P.; Morillon, M. Numerical modelling of heterogeneous sediment transport: New insights for particulate radionuclide transport and deposition. In Proceedings of the Coastal Dynamics, Helsingør, Denmark, 12–16 June 2017; pp. 1767–1778.

86. Bacon, G. *Etude du Transfert Entre L'eau et L'atmosphère D'un Rejet Marin de Tritium (HTO) en Zone Côtière. Etude Préliminaire*; Rapport 3ème année école d'ingénieur; CESI: Caen, France; 2011.

<figure><img src="image_25.png" /></figure>