Impact of Tropical Cyclones on Inhabited Areas of the SWIO Basin at Present and Future Horizons. Part 1: Overview and Observing Component of the Research Project RENOVRISK-CYCLONE

: The international research program “ReNovRisk-CYCLONE” (RNR-CYC, 2017–2021) directly involves 20 partners from 5 countries of the south-west Indian-Ocean. It aims at improving the observation and modelling of tropical cyclones in the south-west Indian Ocean, as well as to foster regional cooperation and improve public policies adapted to present and future tropical cyclones risk in this cyclonic basin. This paper describes the structure and main objectives of this ambitious research project, with emphasis on its observing components, which allowed integrating numbers of innovative atmospheric and oceanic observations (sea-turtle borne and seismic data, unmanned airborne system, ocean gliders), as well as combining standard and original methods (radiosoundings and global navigation satellite system (GNSS) atmospheric soundings, seismic and in-situ swell sampling, drone and satellite imaging) to support research on tropical cyclones from the local to the basin-scale.


Introduction
Due to their highly destructive potential, tropical cyclones (TCs) have long been considered a major risk for populations, territorial economies, and biodiversity. In this regard, predicting their outcome and impacts at present and future times is one of the major concerns of both the Intergovernmental Panel on Climate Change (IPCC) and the World Meteorological Organization (WMO).
As highlighted in the latest reports of WMO's International Workshop on Tropical Cyclones (IWTC), research carried out over the last four years has considerably improved our understanding of TC intensification processes [1,2], as well as TC tracking and intensity forecasting [3,4]. The operational implementation of coupled ocean-atmosphere (OA) numerical weather prediction (NWP) systems by many national weather services has, in particular, played a key role in reducing forecasting errors at all space and time scales [5][6][7][8]. Despite these important advances, additional efforts are still needed to accurately predict and characterize the potential impacts of tropical cyclones on a given territory, especially during landfall. Such efforts include, for instance, the collection of novel atmospheric and oceanic observations, to better constrain (and verify) the performance of coupled NWP systems [9,10], as well as the implementation of wave models and specific microphysical parameterizations to improve roughness, swell, wind speed, and momentum flux representation in TC forecasting systems [11][12][13].
Accurate modelling of OA interactions is particularly crucial in areas such as the tropical south-west Indian Ocean (SWIO) basin (30-90 • E, 0-40 • S), where the atmospheric variability is associated with a particularly strong oceanic response (and vice versa). The SWIO (Figure 1), which contributes to approximately 10-12% of the worldwide cyclonic activity [14][15][16], is indeed widely considered as the cyclonic basin with the highest prevalence of OA interactions [17] due to the unique structure of the thermocline in the Seychelles-Chagos Thermocline Ridge area (55-70 • E, 5-15 • S) [18,19]. Like most TC basins, the SWIO includes many fragile countries, whose economic development and infrastructure, as well as food, medicine and water supply chains, are regularly impacted by tropical cyclones. In very recent years, countries bordering the Mozambique Channel (MC) have indeed been struck by a series of extremely intense and devastating events, whose economic impact will be felt for many years to come. Heavy rains associated with TC DINEO (2017) have caused 700,000 refugees and tens of millions of USD of damage in Mozambique, while the overall cost of TCs ENAWO (2017) and AVA (2018), which affected nearly one million people in Madagascar, was estimated to be more than USD 600 million (about 7% of Madagascar's average annual gross domestic product). These heavy tolls are, however, out of all proportion to those of TCs IDAÏ (considered by the United Nations as the worst natural disaster ever in the MC) and KENNETH (the most intense TC ever reported in the MC [20]), which both made landfall in Mozambique in 2019 [21]. According to the latest economic reports, these two storms have affected a total of nearly 1.7 million people and caused damage and losses estimated at ~USD three billion-plus a further recovery cost estimated at USD 3.4 billion-by the World Bank and Mozambican officials [22].
Given the colossal impact of TCs on the local populations, infrastructure, and economic development of many countries in the SWIO basin, the European Union (EU), together with the Regional Council of Réunion Island and the French State, have designed the transdisciplinary research program "Réunion NOVative research on cyclonic RISKs" (ReNovRisk), to improve the resilience of SWIO countries to TC hazards and mitigate associated economic vulnerability, damages, and risks (e.g., winds, rainfall, landslides, submersion) in inhabited areas. To achieve these objectives, ReNovRisk has been divided into four interlinked research projects, referred to as ReNovRisk-Cyclone, -Erosion, -Impacts and -Transfer, whose overall objectives are described in [23]. The present paper focuses on the Cyclone component of this program, which involves a large international consortium of research institutes, universities, and weather services originating from France (e.g., universities of Réunion Island and Toulouse, Centre National de Recherche Scientifique (CNRS), Météo-France, Institut de Physique de Globe de Paris (IPGP), Institut National de l'Information Géographique et Forestière (IGN), Institut Français de Recherche pour l'Exploitation de la MER (IFREMER)), Mozambique (Eduardo Mondlane and Pemba Unilurio universities, Mozambique Weather Service (INAM)), Madagascar (Institut Supérieur de Technologie de Diego Suarez, Université d'Antananarivo), the Seychelles (Seychelles Meteorological Authority, Seychelles Islands Foundation), Mauritius (Mauritius Oceanography Institute), as well as international institutions such as the European Space Agency (ESA) and WMO, among others.
Through its observing, modelling, climate and outreach components, ReNovRisk-Cyclone (hereafter referred to as RNR-CYC) aims to improve the observation and modelling of TCs, as well as to provide inputs for other components of the global ReNovRisk program focusing on hydrological (RNR-Erosion and RNR-Transfer) and economical In very recent years, countries bordering the Mozambique Channel (MC) have indeed been struck by a series of extremely intense and devastating events, whose economic impact will be felt for many years to come. Heavy rains associated with TC DINEO (2017) have caused 700,000 refugees and tens of millions of USD of damage in Mozambique, while the overall cost of TCs ENAWO (2017) and AVA (2018), which affected nearly one million people in Madagascar, was estimated to be more than USD 600 million (about 7% of Madagascar's average annual gross domestic product). These heavy tolls are, however, out of all proportion to those of TCs IDAÏ (considered by the United Nations as the worst natural disaster ever in the MC) and KENNETH (the most intense TC ever reported in the MC [20]), which both made landfall in Mozambique in 2019 [21]. According to the latest economic reports, these two storms have affected a total of nearly 1.7 million people and caused damage and losses estimated at~USD three billion-plus a further recovery cost estimated at USD 3.4 billion-by the World Bank and Mozambican officials [22].
Given the colossal impact of TCs on the local populations, infrastructure, and economic development of many countries in the SWIO basin, the European Union (EU), together with the Regional Council of Réunion Island and the French State, have designed the transdisciplinary research program "Réunion NOVative research on cyclonic RISKs" (ReNovRisk), to improve the resilience of SWIO countries to TC hazards and mitigate associated economic vulnerability, damages, and risks (e.g., winds, rainfall, landslides, submersion) in inhabited areas. To achieve these objectives, ReNovRisk has been divided into four interlinked research projects, referred to as ReNovRisk-Cyclone, -Erosion, -Impacts and -Transfer, whose overall objectives are described in [23]. The present paper focuses on the Cyclone component of this program, which involves a large international consortium of research institutes, universities, and weather services originating from France (e.g., universities of Réunion Island and Toulouse, Centre National de Recherche Scientifique (CNRS), Météo-France, Institut de Physique de Globe de Paris (IPGP), Institut National de l'Information Géographique et Forestière (IGN), Institut Français de Recherche pour l'Exploitation de la MER (IFREMER)), Mozambique (Eduardo Mondlane and Pemba Unilurio universities, Mozambique Weather Service (INAM)), Madagascar (Institut Supérieur de Technologie de Diego Suarez, Université d'Antananarivo), the Seychelles (Seychelles Meteorological Authority, Seychelles Islands Foundation), Mauritius (Mauritius Oceanography Institute), as well as international institutions such as the European Space Agency (ESA) and WMO, among others.
Through its observing, modelling, climate and outreach components, ReNovRisk-Cyclone (hereafter referred to as RNR-CYC) aims to improve the observation and modelling of TCs, as well as to provide inputs for other components of the global ReNovRisk program focusing on hydrological (RNR-Erosion and RNR-Transfer) and economical (RNR-Impact) consequences of tropical cyclones [23]. Another key objective of RNR-CYC is to foster regional cooperation and improve public policies adapted to present and future TC risks faced by territories bordering the SWIO. The latter is all the more essential because predicted changes in the coupled OA system due to global warming are likely to generate significant modifications of the cyclonic activity in the coming decades. Consequently, regions that are currently spared or moderately affected by TCs, and that often lack experience-based adaptation strategies, may soon have to face potentially increasing TC-related hazards [24,25]. Such changes include, for instance, the widening of the tropical belt resulting from ocean warming [26,27], which has already been shown to induce a poleward migration of TC's lifetime maximum intensity (LMI) in both hemispheres [28][29][30][31], or significant modifications in TC frequency and/or length of the TC season [32].
This paper aims to describing the structure and main objectives of the project RNR-CYC, as well as to present an overview of the main applications and results of its observation component-modelling aspects are presented in more details in the companion paper [33]. This article is organized as follows: Section 2 provides an overview of the four components of RNR-CYC (observation, mesoscale modelling, climate modelling, cooperation, and outreach). Section 3 presents the major achievements of RNR-CYC regarding oceanic and atmospheric observations, while Section 4 concludes and discusses new research topics to be investigated beyond the end of this program.

Structure and Objectives of RNR-CYC
The project RNR-CYC focuses on the meteorological and oceanographic impacts of TCs in the SWIO (Figure 1) at both present and future horizons. It aims, in particular, to better apprehend the impacts of these extreme storms on the main inhabited islands of this oceanic basin by providing innovative modelling and observing products that will also feed the cascade risk analysis tools deployed in the other research components of the global ReNovRisk program [23]. In order to achieve these objectives, RNR-CYC has been divided into four components ( Figure 2): (i) an observation component, to improve both long-term and temporary observations of TCs and their atmospheric and oceanic environments; (ii) a mesoscale modelling component, to improve modelling and short-term forecasting of TCs; (iii) a climate component, to evaluate the consequences of climate change on the variability and structure of TCs at both local and basin scales; and (iv) an outreach component, aimed at improving capacity building in the three aforementioned research areas through strengthening cooperation between SWIO countries. An overview of the structure and objectives of these four components is described hereafter. (RNR-Impact) consequences of tropical cyclones [23]. Another key objective of RNR-CYC is to foster regional cooperation and improve public policies adapted to present and future TC risks faced by territories bordering the SWIO. The latter is all the more essential because predicted changes in the coupled OA system due to global warming are likely to generate significant modifications of the cyclonic activity in the coming decades. Consequently, regions that are currently spared or moderately affected by TCs, and that often lack experience-based adaptation strategies, may soon have to face potentially increasing TC-related hazards [24,25]. Such changes include, for instance, the widening of the tropical belt resulting from ocean warming [26,27], which has already been shown to induce a poleward migration of TC's lifetime maximum intensity (LMI) in both hemispheres [28][29][30][31], or significant modifications in TC frequency and/or length of the TC season [32]. This paper aims to describing the structure and main objectives of the project RNR-CYC, as well as to present an overview of the main applications and results of its observation component-modelling aspects are presented in more details in the companion paper [33]. This article is organized as follows: Section 2 provides an overview of the four components of RNR-CYC (observation, mesoscale modelling, climate modelling, cooperation, and outreach). Section 3 presents the major achievements of RNR-CYC regarding oceanic and atmospheric observations, while Section 4 concludes and discusses new research topics to be investigated beyond the end of this program.

Structure and Objectives of RNR-CYC
The project RNR-CYC focuses on the meteorological and oceanographic impacts of TCs in the SWIO (Figure 1) at both present and future horizons. It aims, in particular, to better apprehend the impacts of these extreme storms on the main inhabited islands of this oceanic basin by providing innovative modelling and observing products that will also feed the cascade risk analysis tools deployed in the other research components of the global ReNovRisk program [23]. In order to achieve these objectives, RNR-CYC has been divided into four components ( Figure 2): (i) an observation component, to improve both long-term and temporary observations of TCs and their atmospheric and oceanic environments; (ii) a mesoscale modelling component, to improve modelling and short-term forecasting of TCs; (iii) a climate component, to evaluate the consequences of climate change on the variability and structure of TCs at both local and basin scales; and (iv) an outreach component, aimed at improving capacity building in the three aforementioned research areas through strengthening cooperation between SWIO countries. An overview of the structure and objectives of these four components is described hereafter.

Observing Component
The observing component of RNR-CYC is aimed at providing additional observations of TCs and their environment by improving regional and local observing capabilities in the SWIO. It has been built around three complementary approaches:

Observing Component
The observing component of RNR-CYC is aimed at providing additional observations of TCs and their environment by improving regional and local observing capabilities in the SWIO. It has been built around three complementary approaches: Atmosphere 2021, 12, 544 5 of 41 (i) A "conventional" approach, based on the reinforcement of regional ground-based meteorological observation facilities and, in particular, of the water vapor Global Navigation Satellite System (GNSS) observation network operated by the International GNSS Service (IGS). Starting in November 2017, 10 new public observation sites (composed of ground-based GNSS receivers and collocated surface weather stations) have been deployed in Madagascar, Eparses Islands, and the Seychelles in the frame of RNR-CYC's sub-program "Indian Ocean GNSS Applications for Meteorology" (IOGA 4 MET [34]) to increase the number of tropospheric GNSS measurements (e.g., zenithal delay, integrated water vapor amounts) and positioning data throughout the western part of the basin (see Section 3.2.2.). The first analysis of GNSS-derived observations collected during RNR-CYC have already demonstrated the benefit of these new permanent stations to investigate the water vapor cycle at diurnal to inter-annual time scales [34,35], but also to provide new and continuous observations to investigate the dynamics of the Earth's crust in this particularly active part of the world [34,36].
(ii) An "experimental" approach, based on the temporary collection of atmospheric and oceanic observations at various locations in the basin. For this purpose, several atmospheric and oceanographic field campaigns have been organized throughout the 3.5-year duration of RNR-CYC, with the goal to provide novel datasets to evaluate numerical model developments and simulations performed in the frame of the project [33]. The main achievements include: The organization of several local observation campaigns to sample sea swell properties during austral winters and summers using acoustic Doppler current profilers (ADCP) and wave gauges deployed near the shore of Réunion Island (see Section 3.1.1.).
(iii) An "exploratory" approach, based on the deployment and evaluation of innovative methods to collect oceanographic observations. A particularly original approach, based on biologging technology, has been evaluated for two years to collect data from sea turtles (ST) equipped with dedicated ARGOS environmental tags in the frame of RNR-CYC's subprogram "Sea Turtle for Ocean Research and Monitoring" (STORM, see Section 3.1.3.). Another original approach, based on the previous work of [37][38][39], was also further investigated to quantify extreme swell phenomena from microseismic noise measurements recorded by ground seismometers (see Section 3.1.2.). The preliminary assessment of terrestrial seismic observations collected in Réunion Island against oceanographic records and offshore wave model data have demonstrated that land-based seismic stations could be particularly useful to observe both austral [40] and cyclonic swell [41] (this Special Issue).

Modelling Component
Protecting life and property requires a precise estimate of the environmental changes associated with the passage of TCs in the vicinity of inhabited areas. The challenge in the face of the cyclonic threat is to simultaneously predict the track and intensity of the storms, but also the consequences resulting from their landfall, or transit near inhabited areas. Hence, damages caused to a given territory, which are essentially related to rainfall intensity, wind strength and sea state (e.g., swell), could significantly differ depending on whether it is affected by a tropical storm (TS), a monsoon depression, or a more or less intense TC. To this end, many operational meteorological services and research centers concerned with TC hazards have made considerable efforts to develop deterministic and ensemble coupled NWP systems providing high spatial resolution forecasts in all TC basins [42][43][44][45][46]). Another fundamental element to improve TC impact predictions, which is partly addressed in the hydrological component (RNR-Erosion) of the global ReNovRisk program, is to also take into account the role of land surface and hydrological processes even after the dissipation of the storm. The recent case of TC IDAÏ, whose catastrophic damage in Beira was mainly caused by flooding occurring days after the storm had passed over the city, is a perfect example of the urgent need to also improve hydrological forecasting.
Improving TC forecasting first and foremost implies a proper representation of the interactions between the storm and the ocean, and vice versa [47]. During the propagation of a TC over an oceanic area, mixing caused by surface winds usually induces a significant drop in surface temperature [48,49] that strongly reduces surface enthalpy and heat fluxes [50,51]. These air-sea fluxes can also be significantly impacted by waves, which redistribute momentum in the near-surface layer and modify the enthalpy fluxes through the emission of sea spray [52][53][54][55][56][57]. In this regard, new parameterizations reproducing the impact of marine aerosols on turbulent heat exchanges have been proposed and validated in recent years [58,59], but are yet to be implemented in atmospheric models to thoroughly evaluate their impact on TC behavior.
Radiative cooling [60,61], evaporation [62,63] and latent heat release [64] have long been recognized to play a key role in the development and intensification of tropical cyclones; therefore, particular attention must also be paid to microphysical schemes implemented in NWP systems. These schemes must, in particular, allow for an efficient representation of the radiative cooling at the top of the storms (which is a constraining criterion of TC intensity) and of the vertical distribution of latent heat (which represents the main source of energy of TCs). They should thus also be able to realistically take into account the role of atmospheric aerosols, that (indirectly) affect the radiation balance by impacting on the radiative and precipitating properties of the clouds. Improving TC forecasting therefore also implies the development of coupled aerosol-microphysical-radiation schemes to be integrated in fully coupled ocean-wave-atmosphere (OWA) models.
An important objective of RNR-CYC was to develop high-resolution OWA and OA modelling systems capable of representing, as exhaustively as possible, the multitude of physical interactions that control the variations of intensity of TCs, as well as their impacts (wind, rain, swell) at the scale of SWIO territories. The main modelling developments made in the frame of this project can be found in [11,65], (this Special Issue, [66]) and in the companion paper [33].

Climate Component
Evaluating the impact of climate change on the frequency and intensity of tropical cyclones is considered as one of the top five issues of concern by the IPCC. Currently, regional and global climate models make it possible to identify the preferred areas of cyclogenesis and occurrence of tropical low-pressure systems at the basin scale, but cannot yet be effectively relied upon to investigate potential changes in their structure and intensification mechanisms. Although it is now widely accepted that the global increase in sea surface temperatures in tropical areas will be a favorable element for TC development, it is not clear how other ingredients involved in their formation and intensification will evolve in the future. In this regard, another important objective of RNR-CYC was to evaluate the global evolution of cyclonic activity in the SWIO, but also to investigate potential structural and intensity changes of TCs resulting from the ongoing modification of their oceanic and atmospheric environments.
The modelling strategy was based on two complementary approaches, allowing to both estimate the evolution of cyclonic activity at the basin scale (i.e., changes in trajectory, intensity, and frequency of TCs at different time scales), and to assess potential structural changes and impacts of TCs at the local scale. This strategy relies on the exploitation of unprecedented high-resolution global climate simulations [32], as well as of mesoscale coupled simulations to estimate the impact of climate change on the intensity, behavior, and consequences of TCs at the scale of a given territory. Examples of results obtained from such high-resolution model runs are discussed in [33,67] (both in this Special Issue).

Regional Cooperation
Another important objective of RNR-CYC is to provide enhanced tools and knowledge to SWIO countries facing cyclonic hazards. Responding to this strong societal issue requires a better structuring of the regional scientific community, as well as significant reinforcement of the cooperation between countries bordering the SWIO basin. While regional collaboration already exists through WMO's regional structures (e.g., Tropical Cyclone Programme (https://public.wmo.int/en/our-mandate/focus-areas/natural-hazards-anddisaster-risk-reduction/tropical-cyclones, accessed on 20 March 2021) and Réunion Island's Regional Specialized Meteorological Center (RSMC La Réunion, http://www.meteo. fr/temps/domtom/La_Reunion/webcmrs9.0/anglais/index.html, accessed on 20 March 2021)) and the Indian Ocean Committee (https://www.commissionoceanindien.org/, accessed on 20 March 2021), interactions remain essentially focused on operational and technical applications, and only modestly promote research development in this area.
One of the strengths of RNR-CYC consists of the implementation of a partnership involving many regional research institutes, universities, and meteorological services that agreed to pool their resources and expertise to strengthen the resilience to TC hazards, and develop public policies better adapted to the risks faced by SWIO territories. RNR-CYC is thus a fundamentally collaborative project that is not only based on the large sharing of data and experiences, but also on training programs and capacity-building initiatives in the fields of observation and forecasting. These actions include, for instance, the organization of forecasting training sessions at RSMC La Réunion, of training courses in climatology and climate change, as well as numerous co-supervised MSc internships based on the analysis of experimental measurements collected in the project.
For instance, the project's field phase, which involved nearly one hundred participants from late January to mid-April 2019, was an opportunity for many students and researchers to initiate themselves to the technique of atmospheric radiosounding (RS). During this 2.5-month campaign,~500 radiosoundings were performed from three experimental sites specifically deployed for this occasion in Maputo (INAM's headquarters, Mozambique), Mayotte (Météo-France weather center, France), and Toamasina (Toamasina international airport, Madagascar). On this occasion, nearly 60 students and academics from Antananarivo (Madagascar), Eduardo Mondlane (Mozambique) and Réunion universities came to Toamasina and Maputo to conduct RS measurements, while many senior forecasters of INAM and Météo-France Mayotte had, for the first time, the opportunity to operate and familiarize themselves with an RS station ( Figure 3).  The involvement of the French consular services in the project also made it possible to communicate widely to the general public, the scientific community, and the media in Mozambique, Madagascar and the Seychelles on the issues of adaptation to climate change and natural hazards. As will be seen in Section 4, these regional collaborations will continue for many years to come through several new research projects initiated by RNR-CYC.

In Situ Swell Observations
One of the main objectives of RNR-CYC is to assess the impact of tropical cyclones on land, including possible submersion resulting from cyclonic swell surge along coastlines. In many tropical islands, the latter are bordered by coral reefs that can serve as a defense against flooding. These reef systems, particularly fringing reefs, protect the coastline by acting as low-pass filters that can reduce the energy of wave flows reaching the coast by up to 98% in the gravitational part of the wave spectrum [68]. The physical processes underlying coral reef coastal protection consist of a complex combination of incident waves, tides, and wind-induced surges [69][70][71].
Ocean wave energy is concentrated in the gravitational frequency band of the wave spectrum, generally between 0.04 Hz and 0.25 Hz. These gravitational waves (GWs) are the main drivers of the hydrodynamics of reef systems, as well as of beaching, runup, and submersion. During and after breaking at the edge of reef systems, GWs are dissipated while low frequency waves (infragravity waves (0.004 < IG (Hz) < 0.04) and very lowfrequency waves (0.001 < VLF (Hz) < 0.004) propagate to the shore. Previous analyses of these propagation and transformation processes across various coral reefs have shown that wave dynamics could vary considerably depending on the characteristics and location of reef systems [71,72].
In order to both quantify the physical processes linked to severe sea states and assess the protective role of reef systems, a cross-shore transect, composed of bottom-fixed wave gauges and ADCPs, was deployed in Réunion Island from February to April in 2019 and 2020. Instruments were installed at the fringing reef of "Trou d'Eau", located along the west coast of the island (Figure 4a). At the near-shore site, reef-base (RS) and reef-flat (RF) stations were deployed through the fringing reef along a cross-shore transect. The RS station was installed at the base of the reef slope at an average depth of 12 m, while three RF stations (RF1, RF2 and RF3 in Figure 4a) were aligned on the flat reef inside the lagoon at a depth of 1 m (Figures 4b,c). RS and RF stations were all equipped with synchronized The involvement of the French consular services in the project also made it possible to communicate widely to the general public, the scientific community, and the media in Mozambique, Madagascar and the Seychelles on the issues of adaptation to climate change and natural hazards. As will be seen in Section 4, these regional collaborations will continue for many years to come through several new research projects initiated by RNR-CYC.

In Situ Swell Observations
One of the main objectives of RNR-CYC is to assess the impact of tropical cyclones on land, including possible submersion resulting from cyclonic swell surge along coastlines. In many tropical islands, the latter are bordered by coral reefs that can serve as a defense against flooding. These reef systems, particularly fringing reefs, protect the coastline by acting as low-pass filters that can reduce the energy of wave flows reaching the coast by up to 98% in the gravitational part of the wave spectrum [68]. The physical processes underlying coral reef coastal protection consist of a complex combination of incident waves, tides, and wind-induced surges [69][70][71].
Ocean wave energy is concentrated in the gravitational frequency band of the wave spectrum, generally between 0.04 Hz and 0.25 Hz. These gravitational waves (GWs) are the main drivers of the hydrodynamics of reef systems, as well as of beaching, runup, and submersion. During and after breaking at the edge of reef systems, GWs are dissipated while low frequency waves (infragravity waves (0.004 < IG (Hz) < 0.04) and very lowfrequency waves (0.001 < VLF (Hz) < 0.004) propagate to the shore. Previous analyses of these propagation and transformation processes across various coral reefs have shown that wave dynamics could vary considerably depending on the characteristics and location of reef systems [71,72].
In order to both quantify the physical processes linked to severe sea states and assess the protective role of reef systems, a cross-shore transect, composed of bottom-fixed wave gauges and ADCPs, was deployed in Réunion Island from February to April in 2019 and 2020. Instruments were installed at the fringing reef of "Trou d'Eau", located along the west coast of the island (Figure 4a). At the near-shore site, reef-base (RS) and reef-flat (RF) stations were deployed through the fringing reef along a cross-shore transect. The RS station was installed at the base of the reef slope at an average depth of 12 m, while three RF stations (RF1, RF2 and RF3 in Figure 4a) were aligned on the flat reef inside the lagoon at a depth of 1 m (Figure 4b,c). RS and RF stations were all equipped with synchronized ocean sensor system instrument (OSSI) wave gauges, enabling continuous recording pressure at a sampling frequency of 10 Hz. The RS station also featured an upward-looking Nortek AQuadopp (AQP) profiler configured to measure current profiles every 20 min, with a 2 Hz hourly burst mode to record wave parameters. In 2020, an ADCP (RDI Sentinel V100) Atmosphere 2021, 12, 544 9 of 41 was also deployed at an ocean offshore (OC) site at a depth of 45 m. This instrument was configured to record incident wave parameters from hourly bursts of 2100 samples at 2 Hz, and current profiles from the bottom to the surface.
here 2021, 12, x FOR PEER REVIEW 9 of 43 ocean sensor system instrument (OSSI) wave gauges, enabling continuous recording pressure at a sampling frequency of 10 Hz. The RS station also featured an upward-looking Nortek AQuadopp (AQP) profiler configured to measure current profiles every 20 min, with a 2 Hz hourly burst mode to record wave parameters. In 2020, an ADCP (RDI Sentinel V100) was also deployed at an ocean offshore (OC) site at a depth of 45 m. This instrument was configured to record incident wave parameters from hourly bursts of 2100 samples at 2 Hz, and current profiles from the bottom to the surface. Observations collected at the offshore stations OC and RS were relied upon to describe and quantify the main properties (height, period, and direction) of the waves impinging on the reef in the GW frequency band (wave periods ranging from 4 to 25 s), while data gathered at RF stations were used to investigate wave transformation (across the fringing reef) and propagation (from outside the lagoon to the shore). Observations collected at OC and RS stations were processed with the RDI software "Velocity" and the PUV method [73], respectively. All OSSI pressure data were corrected from atmospheric mean sea level pressure data recorded at the surface weather station of "La Rivière Des Galets" (located 20 km northward of the experimental site, RDG on Figure 4a) and nonhydrostatic pressure following the linear wave theory [74]. (a) Coastal study area, (b,c) details of the cross-shore transect instrumented during the experiments. The stations are labeled OC for ocean offshore, RS for reef-slope, and RF for reef-flat. In (c), labels OSSI, AQP and RDI ADCP correspond to wave gauges, Nortek AQuadopp current Profiler, and RDI ADCP current profiler, respectively. Instrument deployment dates are indicated in (c).
Observations collected at the offshore stations OC and RS were relied upon to describe and quantify the main properties (height, period, and direction) of the waves impinging on the reef in the GW frequency band (wave periods ranging from 4 to 25 s), while data gathered at RF stations were used to investigate wave transformation (across the fringing reef) and propagation (from outside the lagoon to the shore). Observations collected at OC and RS stations were processed with the RDI software "Velocity" and the PUV method [73], respectively. All OSSI pressure data were corrected from atmospheric mean sea level pressure data recorded at the surface weather station of "La Rivière Des Galets" (located 20 km northward of the experimental site, RDG on Figure 4a) and non-hydrostatic pressure following the linear wave theory [74].
The mean surface elevation was computed from a 20 min moving average of the signal, while the wave spectrum and average wave parameters were obtained through the application of a fast Fourier-transform (FFT). The latter was applied over 4096 data points for incident gravity wave band or swells (0.04 < SW < 0.25 Hz), and 32,768 data values for infragravity (0.004 < IG < 0.04 Hz) and very low frequency (0.001 < VLF < 0.004 Hz) bands.
The recorded signal was not perfectly periodic; therefore, a Hamming window, leading to zero value at the edges, was also applied to mitigate artifacts resulting from leakage issues.
To avoid resulting loss of information at the edges, a 30-min (resp. 3 h) time average was then performed for the incident (resp. IG and VLF) bands. Time series of power spectra density deduced from data collected in 2019 along the cross-shore transect are shown in Figure 5, together with the associated mean power spectra at each station.
Atmosphere 2021, 12, x FOR PEER REVIEW 10 of 43 The mean surface elevation was computed from a 20 min moving average of the signal, while the wave spectrum and average wave parameters were obtained through the application of a fast Fourier-transform (FFT). The latter was applied over 4096 data points for incident gravity wave band or swells (0.04 < SW < 0.25 Hz), and 32,768 data values for infragravity (0.004 < IG < 0.04 Hz) and very low frequency (0.001 < VLF < 0.004 Hz) bands. The recorded signal was not perfectly periodic; therefore, a Hamming window, leading to zero value at the edges, was also applied to mitigate artifacts resulting from leakage issues. To avoid resulting loss of information at the edges, a 30-min (resp. 3 h) time average was then performed for the incident (resp. IG and VLF) bands. Time series of power spectra density deduced from data collected in 2019 along the cross-shore transect are shown in Figure 5, together with the associated mean power spectra at each station.  Figure  4b. The plots are organized in the shoreward direction from the top to the bottom (reef slope RS1 station on the top, and reef flat station RF1 closest to the shore on the bottom). The left panels show the temporal variation of the wave spectrum recorded at each station from February to April 2019. The black solid line indicates the frequency cut at 0.04 Hz between the gravity frequency band (GW) and the infragravity (IG) frequency band. The right panels show the mean power spectral density for each station averaged over the whole period, with the frequencies band GW, IG and VLF (for very low frequency) indicated.
Measurements collected at the reef slope station RS clearly show that the wave energy spectrum is concentrated in the gravity band, with an averaged peak period of 13.5 s over the whole period. After breaking at the reef crest, most of this energy is dissipated (by breaking and/or by friction on the reef bottom), while the remaining part (in IG and VLF bands) is transferred inside the reef system. For the two strong wave events observed in February (likely associated with TC GELENA to the north of Réunion Island) and April 2019 (austral swell event), the total wave energy reduction between the reef slope station RS and the adjacent reef flat station RF3 reaches 98%, a value in good agreement with previous meta-analyses [68]. Within the reef, the transfer of the remaining energy between the three stations was nevertheless slightly different for these two events. In February (TC The plots are organized in the shoreward direction from the top to the bottom (reef slope RS1 station on the top, and reef flat station RF1 closest to the shore on the bottom). The left panels show the temporal variation of the wave spectrum recorded at each station from February to April 2019. The black solid line indicates the frequency cut at 0.04 Hz between the gravity frequency band (GW) and the infragravity (IG) frequency band. The right panels show the mean power spectral density for each station averaged over the whole period, with the frequencies band GW, IG and VLF (for very low frequency) indicated.
Measurements collected at the reef slope station RS clearly show that the wave energy spectrum is concentrated in the gravity band, with an averaged peak period of 13.5 s over the whole period. After breaking at the reef crest, most of this energy is dissipated (by breaking and/or by friction on the reef bottom), while the remaining part (in IG and VLF bands) is transferred inside the reef system. For the two strong wave events observed in February (likely associated with TC GELENA to the north of Réunion Island) and April 2019 (austral swell event), the total wave energy reduction between the reef slope station RS and the adjacent reef flat station RF3 reaches 98%, a value in good agreement with previous meta-analyses [68]. Within the reef, the transfer of the remaining energy between the three stations was nevertheless slightly different for these two events. In February (TC event), the wave energy reduction reached 64% from RF3 to RF2 and 2% from RF2 to RF1, but only 47% from RF3 to RF2 and 15% from RF2 to RF1 in April (austral swell event). The main difference between these two events is related to the period of the incident waves (11.8 s for the TC-related event and 16 s for the austral swell event); this result suggests that the generation and propagation of IG waves inside the reef might be less important for short waves, induced by the wind, than for longer-period swells. This shoreward propagation with less energy reduction is also noticeable in the average spectra, which show a small translation of the peak frequency toward higher frequencies in the IG band, between the RF2 and RF1 station, and a reduction in the peak amplitude.
These preliminary analyses show that Réunion Island's fringing reef plays a strong role in protecting the shore against incident waves. Further investigations will be conducted in order to thoroughly investigate the physical processes and the role of both reef topography and roughness on the dynamics of IG waves. This includes, for instance, the impact of shoreward-propagating IG waves on onshore suspended-sediment transport [75], and the relationship between long wave propagation (and transformation) across the reef and sea water level above the reef-the remaining long waves propagating through the reef have been shown to increase the back-reef set-up and beach runup, which could lead to increased coastal erosion during extreme events such as tropical cyclones or strong storms [76,77].

Ground-Based Swell Observations
The global monitoring of swell activity induced by tropical storms (TS) and TCs is of major interest to quantify the risk associated with extreme swells, but also to validate numerical models used to predict ocean activity. Direct swell observations such as those presented in the previous section are, however, strongly limited by the low number of oceanographic sensors available in this area, as well as by their deployment (and servicing) costs and their vulnerability during tropical cyclones. These limitations motivated the use of indirect observations as alternative and complementary observables to quantify the swell parameters. In this regard, the analysis of wave-induced seismic noise is known to be an interesting substitute for monitoring ocean activity and has been proven to be particularly relevant for assessing the impact of waves on coastal environments [38]. The possibility to derive swell measurements from the seismic noise generated by ocean swell and transmitted to the solid earth as seismic waves recorded by terrestrial seismological instruments [78], is discussed hereafter from data collected in RNR-CYC.
Microseisms recorded by seismic stations worldwide are known to be generated by ocean gravity waves [79]. Seismic energy spectra at terrestrial seismic stations showed two clear peaks in separate frequency bands, known as primary and secondary microseisms (hereafter named PM and SM, respectively), originating from different physical processes involving local or distant sources of ocean wave activity [80]. Primary microseisms (PM) are mostly visible at coastal and island stations and are assumed to be generated by direct interactions of swell-induced pressure variation with the coastal seafloor [81][82][83]. PM noise has the same periods as the ocean swell (i.e., between 8 and 20 s), providing a powerful way to characterize the local impact of swell approaching the shore. On the other hand, SM noise is generally generated in the deep oceans and at larger distances from coastal areas [37,[84][85][86]. It dominates seismic noise at both continental and oceanic stations and exhibits a large peak at half the period of ocean waves (i.e., between 3 and 10 s), generated by the interference of swells of similar periods travelling in opposite directions [79] and generating seismic surface waves travelling horizontally within the solid crust.
In the Indian Ocean, most remote sources of seismic noise are located in the southernmost part of the Austral Ocean basin and are associated with storm systems moving around Antarctica [84,87,88]. Some noise sources may also develop at tropical latitudes in association with tropical cyclones [41] (this Special Issue). Recent seismic deployments on the ocean floor enabled the making of in situ observations of SM underneath TCs in the neighborhood of Réunion Island [37] and confirmed the possibility to track TC and TS from the ocean bottom. Although SM are generally created by distant storms, they can also be generated by the coastal reflection of waves if incident and reflected waves propagate in opposite directions [89,90]. In this latter case, the incoming swell may interfere with its reflected swell, resulting in the generation of standing waves close to coastal areas that oscillate at twice the frequency of the propagating wave [91]. Some studies suggest that local and distant sources of noise in the SM frequency peak may coexist [39,92,93].
Previous seismic analyses conducted in the Pacific and Indian Oceans have already demonstrated that several swell parameters can be derived from the seismic data. This includes: (i) the swell peak period Tp, derived from the dominant frequency of the PM and SM using the power spectral density analysis of the seismic records; (ii) the local or distant wave significant height Hs (in the case of the PM or SM band, respectively), obtained by measuring the microseism amplitudes through hourly Root Mean Square (RMS); and (iii) the wave peak direction Dp for the case of the PM band, or the source direction in the case of the SM band, which can both be deduced from the polarization analysis of the three seismic data components, to determine the dominant direction and strength of the recorded microseismic noise.
To illustrate this innovative approach, we analyse hereafter seismic observations collected during the tropical storm (TS) ELIAKIM that developed in March 2018. This storm has clear signatures, despite not reaching tropical cyclone intensity, and developed at a large distance from Réunion Island, demonstrating the potential of the method in quantifying remote systems. As mapped on Figure  Our analysis was performed from data recorded by seismic stations of the permanent seismic network of the Piton de la Fournaise Volcano Observatory (OVPF, seismic code PF, red triangles on Figure 6b) and from a temporary seismic network deployed in the frame of the RNR program (code ZF, blue triangles on Figure 6b). Seismic data were compared with swell parameters issued from the WaveWatch3 (WW3) model [94,95] at nodes surrounding the island from the global wave model hindcast [96], indicated by stars on Figure 6b.  The hourly RMS amplitudes of the seismic noise recorded at the various seismic stations on the island are shown in Figure 7b,c, for both the SM and the PM bands, respectively. Note that the vertical axes are different, and that the SM amplitude is almost one order of magnitude larger than the PM. The 21 seismic stations analyzed over Réunion island display similar variation patterns. This indicates that this noise is not a purely local source; the storm acting as a distant SM source, or that the source-if local-is larger than the size of the island, which is the case for the swell generating the PM. The observed station-to-station variability in RMS amplitudes suggests the influence of local site effects such as the installation, the coupling of the seismic sensor with the ground, the nature of the bedrock, and the attenuation around the station.
The RMS peaks at seismic station MAT (Figure 7b) are observed on 16 and 19 March. This period corresponds to the maximum intensity of TS ELIAKIM and to the local variation of the swell height, as visible on Hs data issued from the WW3 model at various points around Réunion island ( Figure 8). The maximum Hs observed on 19 March (corresponding to the nodes at the longitude of 55 • E) interacted with the local bathymetry slightly later, which explains the small delay with the PM. A small peak was also observed on 21 March in both SM and PM data with dominating energy at period~10s, likely originating from a distant source and not related to any local swell activity increase. ere 2021, 12, x FOR PEER REVIEW 14 of 43 The hourly RMS amplitudes of the seismic noise recorded at the various seismic stations on the island are shown in Figures 7b,c, for both the SM and the PM bands, respectively. Note that the vertical axes are different, and that the SM amplitude is almost one order of magnitude larger than the PM. The 21 seismic stations analyzed over Réunion island display similar variation patterns. This indicates that this noise is not a purely local source; the storm acting as a distant SM source, or that the source-if local-is larger than the size of the island, which is the case for the swell generating the PM. The observed station-to-station variability in RMS amplitudes suggests the influence of local site effects such as the installation, the coupling of the seismic sensor with the ground, the nature of the bedrock, and the attenuation around the station.
The RMS peaks at seismic station MAT (Figure 7b) are observed on 16 and 19 March. This period corresponds to the maximum intensity of TS ELIAKIM and to the local variation of the swell height, as visible on Hs data issued from the WW3 model at various points around Réunion island ( Figure 8). The maximum Hs observed on 19 March (corresponding to the nodes at the longitude of 55° E) interacted with the local bathymetry slightly later, which explains the small delay with the PM. A small peak was also observed on 21 If one accepts that the SM is generated in the vicinity of the storm center and that the PM is generated in coastal areas closer to the seismic stations, one should observe a delay between the two curves. Such a delay was previously observed for austral swells generated by distant storms [40] and proposed as a precursor for predicting coastal submersion in Réunion island. In the present case, considering the involved distances (TS ELIAKIM was located~800 km of the coast of Réunion Island on 16 March, and 500 km on 19 March) and the involved velocities (3 km s −1 for the surface waves carrying the SM signal, and 50 km h −1 for the long period waves at the surface of the ocean), one should expect a delay of >10 h for each peak. However, Figure 8 shows a delay of only~4 h between the SM and PM. This may indicate that the SM is not generated at the storm center, but likely closer to the island. Alternatively, it may also indicate a slow wind-wave growth that may take from a few hours to few days [97]. The swell-noise amplitude correlation can be used to build a transfer function to translate the amplitude of the seismic noise in terms of swell height. Although such a relation is station-related, and therefore not universal, it was nevertheless shown to provide particularly good results for strong swell events [40,41,81]. March in both SM and PM data with dominating energy at period ~10s, likely originating from a distant source and not related to any local swell activity increase. If one accepts that the SM is generated in the vicinity of the storm center and that the PM is generated in coastal areas closer to the seismic stations, one should observe a delay between the two curves. Such a delay was previously observed for austral swells generated by distant storms [40] and proposed as a precursor for predicting coastal submersion in Réunion island. In the present case, considering the involved distances (TS ELIAKIM was located ~800 km of the coast of Réunion Island on 16 March, and 500 km on 19 March) and the involved velocities (3 km s −1 for the surface waves carrying the SM signal, and ~50 km h −1 for the long period waves at the surface of the ocean), one should expect a delay of >10 h for each peak. However, Figure 8 shows a delay of only ~4 h between the SM and PM. This may indicate that the SM is not generated at the storm center, but likely closer to the island. Alternatively, it may also indicate a slow wind-wave growth that may take from a few hours to few days [97]. The swell-noise amplitude correlation can be used to build a transfer function to translate the amplitude of the seismic noise in terms of swell height. Although such a relation is station-related, and therefore not universal, it was nevertheless shown to provide particularly good results for strong swell events [40,41,81].
The recordings of the three components of the ground motion also make the seismic noise a vectorial observation for measuring the direction and strength of the signal polarization. In the case of the SM band, the polarization is expected to point towards the noise source, i.e., the storm center, whereas it should indicate the very local swell propagation direction in the case of the PM. In this PM case, some angle may exist between the swell propagation direction at a large distance offshore and the swell at the coast, due to refraction of the swell and the coastal bathymetry. In the case of TS ELIAKIM, the polarization was measured on an hourly basis in the SM frequency band. Figure 6b displays the average polarization recorded on 19 March at the various land seismic stations in the 0.1-0.33 Hz frequency band (pink) that can be compared to the azimuth computed from the storm The recordings of the three components of the ground motion also make the seismic noise a vectorial observation for measuring the direction and strength of the signal polarization. In the case of the SM band, the polarization is expected to point towards the noise source, i.e., the storm center, whereas it should indicate the very local swell propagation direction in the case of the PM. In this PM case, some angle may exist between the swell propagation direction at a large distance offshore and the swell at the coast, due to refraction of the swell and the coastal bathymetry. In the case of TS ELIAKIM, the polarization was measured on an hourly basis in the SM frequency band. Figure 6b displays the average polarization recorded on 19 March at the various land seismic stations in the 0.1-0.33 Hz frequency band (pink) that can be compared to the azimuth computed from the storm center position. This map illustrates a very homogeneous orientation of the noise polarization across the island, pointing to the SW (toward the storm center location on that date) and confirming that the SM originated in the vicinity of the storm center.
These results confirm that seismic noise may provide a useful proxy to quantify the swell parameters. In the PM band (10-20 s period), the seismic noise amplitude allows for a robust quantification of the swell height Hs through a transfer function, the spectrogram allows deciphering the swell dominant period Tp, and the polarization analysis allows retrieving the local swell direction Dp. In the SM band (3-10 s period), the amplitude revealed the storm strength and the polarization indicated the storm azimuth. Terrestrial seismic stations therefore provide alternative and complementary observations of both TC and ocean activities. In some cases, the availability of several decades of seismic archives may also provide new opportunities to derive cyclone climatologies [41].

Biologging Observations
Observing the vertical structure of the ocean is essential to improve knowledge of both the coupled OA system and marine ecosystems. In this regard, an increasingly common alternative to gather high-resolution hydrographic profiles in the world's oceans is to rely on animal-borne sensors (a.k.a. biologging) to collect in situ observations in remote and under-instrumented areas. Compared to conventional oceanographic in situ observation approaches (e.g., gliders, ARGO drifters, buoys, research cruises), animal-borne electronic ARGOS tags are relatively inexpensive to operate and can be deployed in remote areas with limited human resources. This approach offers all countries the possibility to actively contribute to the collection of ocean observations; therefore, biologging is expected to grow considerably in the future. The recent decision, in August 2020, of the Global Ocean Observing System (GOOS)'s Executive Committee to create a new observing network exclusively dedicated to animal-borne ocean sensors ("ANIBOS") is definitely in line with this perspective and clearly attests of the immense potential of this approach.
The potential of biologging for sampling the thermal structure of the tropical Indian Ocean was evaluated in the frame of RNR-CYC's sub-program "Sea Turtle for Ocean Research and Monitoring" (STORM). STORM was initiated in January 2019 by Réunion's University and Réunion Island's Sea-Turtle Observatory (Kelonia), with the goal of monitoring the state of the tropical Indian Ocean down to several hundred meters below the surface in near-real-time and at high spatial (<100 m) and temporal (5') resolutions. As of March 2021, 22 animals have been equipped with Temperature-Depth (TD) ARGOS tags before being released from Réunion Island (Figure 9). Note that all animal manipulations were made by accredited ST biologists and that all STs were accidentally captured in fisherman's nets in the vicinity of Réunion Island and brought back to Kelonia's care center to be healed and rehabilitated. During this experiment, two species of late juvenile sea turtles (loggerhead and Olive Ridley) were equipped with Argos TD tags. While some animals stayed in their pelagic habitats, some loggerheads also began their first homing migration to the Oman Gulf, thus enabling the collection of data both in tropical areas, from loggerheads, and subtropical areas, from Olive Ridleys (Figure 9). STs released from Réunion Island principally swam at (or slightly below) the surface (~50% of the time) and near the bottom of the ocean mixed layer (OML, ~25% of the time). They were found to dive up to 100 times a day, sometimes up to 350 m, allowing the collection of numerous hydrographic profiles within and far below the OML [9]. The analysis of data collected during the first year of this experiment has confirmed the great potential of this approach for sampling the vertical structure of the ocean, validating ocean models and spaceborne sensors, as well as to investigate the intra-seasonal variability of the tropical Indian Ocean [9].
Sea turtles are known to generally swim in rings and frontal areas in between, and often travel by moving from one ring to another [98]. These rings and eddies are numerous in the Mozambique Channel [99,100], but less common north of Réunion Island, which During this experiment, two species of late juvenile sea turtles (loggerhead and Olive Ridley) were equipped with Argos TD tags. While some animals stayed in their pelagic habitats, some loggerheads also began their first homing migration to the Oman Gulf, thus enabling the collection of data both in tropical areas, from loggerheads, and subtropical areas, from Olive Ridleys (Figure 9). STs released from Réunion Island principally swam at (or slightly below) the surface (~50% of the time) and near the bottom of the ocean mixed layer (OML,~25% of the time). They were found to dive up to 100 times a day, sometimes up to 350 m, allowing the collection of numerous hydrographic profiles within and far below the OML [9]. The analysis of data collected during the first year of this experiment has confirmed the great potential of this approach for sampling the vertical structure of the ocean, validating ocean models and spaceborne sensors, as well as to investigate the intra-seasonal variability of the tropical Indian Ocean [9].
Sea turtles are known to generally swim in rings and frontal areas in between, and often travel by moving from one ring to another [98]. These rings and eddies are numerous in the Mozambique Channel [99,100], but less common north of Réunion Island, which may explain that loggerhead STs released from there generally tend to follow a straight trajectory (Figure 9)-another explanation (currently under investigation) is that STs rarely feed during their reproductive migration, so as to reach their breeding areas as fast as possible. Nevertheless, according to its circular trajectories, ST Union (Olive Ridley), which moved south-east of Réunion Island in late 2020, apparently travelled in oceanic gyres and eddies (Figure 9).
The surface current analysis of Mercator Ocean's operational model NEMO-PSY4 (a.k.a Glo12, [101]), averaged from 15 October to 15 December, confirmed that ST Union indeed swam in cyclonic and anticyclonic eddies of variable sizes during this period (Figure 10a). Vertical cross-sections of the ocean temperature field across the center of the main eddy (Figure 10b,c) show strong upwellings at the edges of the vortex (while downwelling can be observed at the center), resulting in significant temperature gradients at the surface (Figure 10a). This vertical transport, which brings nutrient-rich waters from the thermocline up to the surface, makes these eddies prime feeding areas for many marine species and are particularly appreciated by sea turtles. Equipping STs with environmental tags is thus an easy and relatively affordable way to sample the properties of these important transient mesoscale features, which will be further investigated in the continuation of the STORM program (see Section 4).    (Figure 11a), and STs India and Tom that were trapped in TC HEROLD (March 2020) during its intensification phase (Figure 11b). Spaceborne oceanic observations are generally unavailable under cloudy conditions; therefore, data collected by these animals also represent a fantastic opportunity to investigate the impact of TCs on the surface and subsurface structure of the ocean.
Sea surface temperature data collected by ST India in TC HEROLD from 14 to 20 March 2020 are shown in Figure 11c. During these six days, the animal remained trapped in the immediate vicinity of the storm center (~30 km from the TC eye) and only moved~50 km in the north/north-west direction. In situ surface temperature observations collected from 14 to 17 March showed SST cooling of~3.5 • C (29 • C to 25.5 • C) in 72 hduring these three days, the animal remained quasi-stationary. As the storm progressively moved south-eastwards, ST India began to slowly move to the north-west, over the area previously affected by the tropical cyclone. Observations collected from 17 to 20 March showed that the temperature surface layer in this area quickly returned to pre-storm conditions, to regain its initial temperature of 29 • C on 19 March. As shown in [33], surface and subsurface observations collected by STs in tropical cyclones can also provide key data to evaluate coupled model forecasts in cyclonic conditions.

Glider Observations
On 22 January 2019, two SlocumG1 gliders operated by CNRS were deployed Réunion Island for a period of two months. The two instruments were programm follow a north-east (glider GLNE) and a north-west (glider GLNW) trajectory to reac northernmost region of the Mascarene Archipelago, where a high probability of T mation was suggested about 400-600 km north of Réunion Island (Figure 12). The acquisition strategy was set in order to complete a saw-tooth navigation pattern, allo the gliders to dive with an angle of 26° between 5 m and 980 m depth (resulting along-track resolution of about 4 km once the profile was normalized on the vertica

Glider Observations
On 22 January 2019, two SlocumG1 gliders operated by CNRS were deployed from Réunion Island for a period of two months. The two instruments were programmed to follow a north-east (glider GLNE) and a north-west (glider GLNW) trajectory to reach the northernmost region of the Mascarene Archipelago, where a high probability of TC formation was suggested about 400-600 km north of Réunion Island ( Figure 12). The data acquisition strategy was set in order to complete a saw-tooth navigation pattern, allowing the gliders to dive with an angle of 26 • between 5 m and 980 m depth (resulting in an along-track resolution of about 4 km once the profile was normalized on the vertical).  The two gliders were programmed to sample the ocean during descending (downcast) and ascending (upcast) profiles. Observations were transmitted in real time by satellite telemetry after each upcast, when directives for modifying the sampling strategy and glider navigation (based on operational ocean model forecasts) were also received. Each glider was equipped with various physical and optical biogeochemical instruments for sampling the ocean temperature, salinity, oxygen, turbidity, and chlorophyll-a concentration at different rates according to depth, as shown in Table 1. All sensors were operational for glider GLNW, but the oxygen and optical sensors were turned off on glider GLNE.
The primary objective of this experiment was to investigate OA kinetic energy exchanges in cyclonic conditions, with emphasis on the fraction of kinetic energy transmitted to the ocean. As shown by [102], who analyzed the vertical structure of the ocean in a tropical cyclone sampled during the CIRENE field phase [103], this kinetic energy is generated on the left side of the TC track in the southern hemisphere, and later consumed by strong vertical mixing resulting from surface water cooling [104]. Due to the presence of strong currents (>0.25 m s −1 ) along the glider trajectories, an important surface drift prevented the full application of the initially planned northward navigation strategy. The navigation parameters were then modified early in the mission to mitigate the battery consumption by reducing the depth range of the glider (because the oil pump that controls the glider buoyancy is the device that consumes the most energy). The gliders profiles were thus reduced to 500 m and 300 m depths (as shown by the colored tracks on Figure  12b,c). Real-time data transmission was also turned-off at some points to reduce the time spent at the surface, when drifting was maximized. The slower-than-expected displacement speed of the gliders prevented them from intercepting the core of TC GELENA, which developed in the area between 5 and 9 February 2019 (Figure 12a). Thanks to realtime trajectory optimization, both systems were nevertheless able to approach relatively close to the storm, as shown by black ellipses displayed in Figure 12a. For instance, glider GLNE was able to move to a distance of about 290 km from the TC center and to collect ocean data in the direct vicinity of this storm. The two gliders were programmed to sample the ocean during descending (downcast) and ascending (upcast) profiles. Observations were transmitted in real time by satellite telemetry after each upcast, when directives for modifying the sampling strategy and glider navigation (based on operational ocean model forecasts) were also received. Each glider was equipped with various physical and optical biogeochemical instruments for sampling the ocean temperature, salinity, oxygen, turbidity, and chlorophyll-a concentration at different rates according to depth, as shown in Table 1. All sensors were operational for glider GLNW, but the oxygen and optical sensors were turned off on glider GLNE.
The primary objective of this experiment was to investigate OA kinetic energy exchanges in cyclonic conditions, with emphasis on the fraction of kinetic energy transmitted to the ocean. As shown by [102], who analyzed the vertical structure of the ocean in a tropical cyclone sampled during the CIRENE field phase [103], this kinetic energy is generated on the left side of the TC track in the southern hemisphere, and later consumed by strong vertical mixing resulting from surface water cooling [104]. Due to the presence of strong currents (>0.25 m s −1 ) along the glider trajectories, an important surface drift prevented the full application of the initially planned northward navigation strategy. The navigation parameters were then modified early in the mission to mitigate the battery consumption by reducing the depth range of the glider (because the oil pump that controls the glider buoyancy is the device that consumes the most energy). The gliders profiles were thus reduced to 500 m and 300 m depths (as shown by the colored tracks on Figure 12b,c). Real-time data transmission was also turned-off at some points to reduce the time spent at the surface, when drifting was maximized. The slower-than-expected displacement speed of the gliders prevented them from intercepting the core of TC GELENA, which developed in the area between 5 and 9 February 2019 (Figure 12a). Thanks to real-time trajectory optimization, both systems were nevertheless able to approach relatively close to the storm, as shown by black ellipses displayed in Figure 12a. For instance, glider GLNE was able to move to a distance of about 290 km from the TC center and to collect ocean data in the direct vicinity of this storm. Vertical profiles of temperature, salinity and density collected by the two gliders between 0 and 150 m and 23 January and 17 February 2019 (along the two northward segments of the tracks) are shown in Figure 13, together with corresponding mean profiles for each parameter. The mixed layer depth (MLD), indicated by a black solid line, is estimated as the depth where the temperature differs by 0.2 • C from a surface reference value of 10 dbar [105,106]. Due to the significant distance separating the two gliders from the TC core, and because the maximum kinetic energy generated in the upper ocean occurred on the left side of the TC track, no significant change in the surface temperature could be noticed during the pre-cyclonic, cyclonic, and post-cyclonic phases (Figure 13a,d,g). A relatively constant temperature of 29 • C was observed between the surface and 30 to 40 m depth, but quickly decreased to reach 22 • C at a depth of 150 m. Salinity measurements (Figure 13b,e,h) indicated a freshening of the water within the top 100 m layer, starting around 6 February (−19.22 • N) for the GLNW track, and from 28 January (−20.15 • N) for the GLNE track. While part of this freshening can be explained by the rainfall generated by TC GELENA between 15 • S and 5 • N, another possible explanation is related to the advection of fresher water originating from the southern branch of the South Equatorial Current, which flows westward with a 20 Sv transport to feed the East Madagascar Current [107]. The velocity field computed from Mercator-Ocean's global reanalysis model PHY-001-030, which indicated the presence of strong westward currents between 28 January and 9 February (not shown), and the vertical extent of this freshwater lens, both seem to support this hypothesis. Density measurements (Figure 13c,f,i) did not seem to be impacted by this freshening and showed a well-mixed layer within the top 40 m of the ocean and a well-stratified layer underneath.
The mixed layer depth shows significant frequency variations in relation to the diurnal cycle, and also appears deeper along the GLNW track until 8 February. GLNW was always located southward of GLNE before this date; therefore, these observations suggest a sloping up of the MLD in the northward direction. As GLNW moved northward of GLNE after 8 February, the MLD measured by GLNE became deeper, thus reinforcing the hypothesis of a northward sloping-up of MLD. A significant rise of the MLD of up to 20 m (on 5 and 6 February for GLNE and GLNW, respectively) could also be observed almost immediately after the TC formation (4 February). The MLD then returned to its initial value between 9 February (for GLNW) and 12 February (for GLNE). The mixed layer depth shows significant frequency variations in relation to the diurnal cycle, and also appears deeper along the GLNW track until 8 February. GLNW was always located southward of GLNE before this date; therefore, these observations suggest a sloping up of the MLD in the northward direction. As GLNW moved northward of GLNE after 8 February, the MLD measured by GLNE became deeper, thus reinforcing the hypothesis of a northward sloping-up of MLD. A significant rise of the MLD of up to 20 m (on 5 and 6 February for GLNE and GLNW, respectively) could also be observed almost immediately after the TC formation (4 February). The MLD then returned to its initial value between 9 February (for GLNW) and 12 February (for GLNE).
The temperature-salinity (TS) diagrams collected along the full glider tracks (23 January to 23 March 2019) are shown in Figure 14. The water masses sampled by the two instruments, known as the Indian Central Water, were formed and subducted in the Subtropical Convergence area of the southern gyre of the Indian Ocean [107,108]. The OML in this water mass was generally between 50 and 150 m in depth, with temperature between 25

UAS Observations
During RNC-CYC field phase, a UAS system called "BOREAL" was operated from Réunion Island to sample the lower part of the atmosphere (Figure 15). The main objective of this experiment was to provide in situ measurements of air-sea interactions in cyclonic and pre-cyclonic conditions to evaluate numerical models and microphysical parameter-

UAS Observations
During RNC-CYC field phase, a UAS system called "BOREAL" was operated from Réunion Island to sample the lower part of the atmosphere (Figure 15). The main objective of this experiment was to provide in situ measurements of air-sea interactions in cyclonic and pre-cyclonic conditions to evaluate numerical models and microphysical parameterizations developed in the frame of the program.
The BOREAL UAS (Boréal SAS, Toulouse, France) is a fixed-wing aircraft with a 4.2 m wingspan propelled by a thermal engine (Figure 15b) that has a maximum takeoff weight of 25 kg and can fly in winds up to 25 m s −1 . For this experiment, a 5 kg scientific payload was developed to study air-sea interactions using an optical particle counter for measuring the aerosol number and size distribution (0.3 < diameter < 3.0 µm; MetOne), a custom-designed multi-hole probe for measurements of turbulence and threedimensional winds, a radar altimeter for wave height and sea state measurements, a broadband shortwave pyranometer for downwelling solar radiation (Licor), an infra-red temperature sensor to measure sea-surface temperature, as well as standard meteorological measurements (atmospheric pressure, temperature and relative humidity). With this payload, the endurance and range of the BOREAL UAS was~7 h, 700 km, and 2000 m above sea level (asl). A live video was also streamed up to 40 km from the ground-station to provide additional safety during low-altitude segments when the UAS flew~40 m asl. The BOREAL UAS is autonomous, although its flight plan could be adapted at any time to accommodate weather conditions or air traffic via a radio or a satellite link. A transponder mounted on its wing also allowed it to integrate into international airspace.
The BOREAL UAS flew over the Indian Ocean in two exclusive zones allocated by the French aviation authority to the southeast and northwest sides of Réunion Island (Figure 15a). These two areas were dedicated for scientific flights over the international waters up to 250 km from the ground station and with a ceiling at 1067 m above sea level. BOREAL UAS operations were conducted from the airfields of Cambaie (when flying towards the north of the island) and Bras Panon (when flying towards the south). In total, 12 scientific flights were carried out between February and March 2019, for a total of 53.7 h of research flights and 5012 km of observations over the ocean. Three flights of more than 200 km from the ground station were carried out, the longest of which lasted 6 h and 18 min and covered 610 km. Examples of observations collected by the BOREAL UAS during this two-month operating period are presented below.  Figure 16 summarizes the BOREAL UAS measurements of aerosol concentrations for particle diameters >0.3 μm and >1.0 μm within 200 m above the ocean surface, compared to horizontal wind speed (Figure 16a) and wave height (Figure 16b) for conditions encountered during the two-month campaign. As expected, there was an increase in particle   (Figure 16b) for conditions encountered during the two-month campaign. As expected, there was an increase in particle concentrations (diameter > 0.3 µm), which is often associated with primary marine aerosol (PMA) emissions, over the observed range of wind speed (2.2 to 13.5 m s −1 ) and wave heights (2 to 3.7 m). However, the range of wind speeds and wave heights encountered in the vicinity of Réunion Island during the RNR-CYC field campaign remained relatively small, partly because the paths of the TCs were never closer than a few hundred kilometers to the ground station during the observed period.  Figure 16 summarizes the BOREAL UAS measurements of aerosol concentratio particle diameters >0.3 μm and >1.0 μm within 200 m above the ocean surface, comp to horizontal wind speed (Figure 16a) and wave height (Figure 16b) for condition countered during the two-month campaign. As expected, there was an increase in pa concentrations (diameter > 0.3 μm), which is often associated with primary marine ae (PMA) emissions, over the observed range of wind speed (2.2 to 13.5 m s −1 ) and heights (2 to 3.7 m). However, the range of wind speeds and wave heights encounter the vicinity of Réunion Island during the RNR-CYC field campaign remained rela small, partly because the paths of the TCs were never closer than a few hundred ki ters to the ground station during the observed period.   (Figure 17a). The signature of this storm could be seen up to several hundred kilometers away in observations of ocean waves, atmospheric turbulence structures, and enhanced PMA emissions. These emissions, which occur over a much larger domain than the cyclonic system, have been shown to modify the TC's track and intensity [11].
As TC JOANINHA moved to the east of Mauritius, the BOREAL UAS flew to the west of Réunion Island to sample the TC's perimeter (see flight track in Figure 17a and wave crests, swells and clouds generated by the cyclone recorded by the on-board camera in Figure 17b). In situ measurements were used, in particular, to assess parameterizations of air-sea interactions simulated with the OWA coupled system developed in RNR-CYC (see [33] for details about this system). Results presented hereafter are derived from a coupled OWA simulation based on the ocean model Coastal and Regional Ocean Community (CROCO, http://www.croco-ocean.org, accessed on 20 March 2021), the wave model Wave Watch 3 (WW3) [94,95], and the atmospheric model Meso-NH (http: //mesonh.aero.obs-mip.fr/, accessed on 20 March 2021), which was developed specifically for comparisons with the BOREAL UAS. Figure 18 presents a vertical profile of aerosol particles conducted during the 25 March flight from near the ocean surface (<50 m asl) to above the marine boundary layer (~1000 m asl). Observed aerosol concentrations (diameter > 0.3 µm) were found to be relatively constant throughout the marine boundary layer (up to~800 m asl) and decreased above the inversion. A surface layer with enhanced emissions, similar to the profiles observed at the surf zone in the frame of project Miriad [109], was also observed at the lowest part of the vertical profile (~50 m asl). Figure 18 presents a vertical profile of aerosol particles conducted during the 25 March flight from near the ocean surface (<50 m asl) to above the marine boundary layer (~1000 m asl). Observed aerosol concentrations (diameter > 0.3 μm) were found to be relatively constant throughout the marine boundary layer (up to ~800 m asl) and decreased above the inversion. A surface layer with enhanced emissions, similar to the profiles observed at the surf zone in the frame of project Miriad [109], was also observed at the lowest part of the vertical profile (~50 m asl). This surface layer, which was not captured in the model simulation ( Figure 18), has not been well-documented so far and has important implications for the transport of aerosols into the well-mixed boundary layer. In addition, the simulations tend to overestimate the aerosol emissions near the ocean surface (i.e., an approximate factor of two enhancement in aerosol concentrations) and show a pronounced vertical gradient throughout the marine boundary layer with an underestimate of aerosol concentrations in the free troposphere. These large differences between observed and simulated vertical profiles of aerosol concentrations suggest that mixing timescales in the marine boundary layer may perhaps be too slow in the model-resulting in an excessive gradient in the aerosol concentration-but could also be related to a possible bias in the CAMS/COPERNICUS analyses used to initialize sea salts in Meso-NH. Although a mass distribution calculation was made to adjust the consistency of aerosol distributions in CAMS and Meso-NH, a small remaining mass difference in the finest modes can also lead to large differences in the number of aerosols, especially in areas with low primary aerosol production, i.e., where aerosol number concentrations mostly arise from initial conditions. This surface layer, which was not captured in the model simulation ( Figure 18), has not been well-documented so far and has important implications for the transport of aerosols into the well-mixed boundary layer. In addition, the simulations tend to overestimate the aerosol emissions near the ocean surface (i.e., an approximate factor of two enhancement in aerosol concentrations) and show a pronounced vertical gradient throughout the marine boundary layer with an underestimate of aerosol concentrations in the free troposphere. These large differences between observed and simulated vertical profiles of aerosol concentrations suggest that mixing timescales in the marine boundary layer may perhaps be too slow in the model-resulting in an excessive gradient in the aerosol concentration-but could also be related to a possible bias in the CAMS/COPERNICUS analyses used to initialize sea salts in Meso-NH. Although a mass distribution calculation was made to adjust the consistency of aerosol distributions in CAMS and Meso-NH, a small remaining mass difference in the finest modes can also lead to large differences in the number of aerosols, especially in areas with low primary aerosol production, i.e., where aerosol number concentrations mostly arise from initial conditions. A few days later, a southern swell event occurred when winds shifted to the north-east direction as TC JOANINHA continued its eastward progression through the Indian Ocean. Flight operations were moved to Bras Panon and the BOREAL UAS flew along a southern curtain for more than 200 km to the east of Réunion Island on 29 March 2019 from 6:33 to 12:04 UTC (Figure 19a). The platform followed a flight plan with ascents from 100 to 1000 A few days later, a southern swell event occurred when winds shifted to the north-east direction as TC JOANINHA continued its eastward progression through the Indian Ocean. Flight operations were moved to Bras Panon and the BOREAL UAS flew along a southern curtain for more than 200 km to the east of Réunion Island on 29 March 2019 from 6:33 to 12:04 UTC (Figure 19a). The platform followed a flight plan with ascents from 100 to 1000 m asl followed by straight-level legs at 400 and 100 m asl for at least 10 km each. Mean and significant wave height observations collected during this flight by the onboard radar altimeter showed that wave height increased by about one meter (from~2.5 to 3.5 m) between locations A and D (Figure 19b). A few days later, a southern swell event occurred when winds shifted to the north-east direction as TC JOANINHA continued its eastward progression through the Indian Ocean. Flight operations were moved to Bras Panon and the BOREAL UAS flew along a southern curtain for more than 200 km to the east of Réunion Island on 29 March 2019 from 6:33 to 12:04 UTC (Figure 19a). The platform followed a flight plan with ascents from 100 to 1000 m asl followed by straight-level legs at 400 and 100 m asl for at least 10 km each. Mean and significant wave height observations collected during this flight by the onboard radar altimeter showed that wave height increased by about one meter (from ~2.5 to 3.5 m) between locations A and D (Figure 19b). A similar increase in wave height was also obtained from the simulations of this case study (Figure 19a). The associated modelled ( Figure 20a) and observed (Figure 20b) spectral density functions provide insight on the composition of the southern swell and show two peaks for both observations and simulations. Azimuth plots from WW3 based on the wave elevation spectrum (Figure 20c) indicate that these peaks corresponded to distinct components of the southern swell; the first component was generated by the northerly winds (red arrow in Figure 20c), and the second component was generated by TC JOAN-INHA to the east (black arrow in Figure 20c). In spite of the observed increase in wave A similar increase in wave height was also obtained from the simulations of this case study (Figure 19a). The associated modelled ( Figure 20a) and observed (Figure 20b) spectral density functions provide insight on the composition of the southern swell and show two peaks for both observations and simulations. Azimuth plots from WW3 based on the wave elevation spectrum (Figure 20c) indicate that these peaks corresponded to distinct components of the southern swell; the first component was generated by the northerly winds (red arrow in Figure 20c), and the second component was generated by TC JOANINHA to the east (black arrow in Figure 20c). In spite of the observed increase in wave height (and relatively constant wind field), a horizontal gradient in the aerosol concentration and size distribution at 100 m asl was not observed during this swell event. These results were predictable because the 1 m increase in wave height is relatively small and not expected to generate a significant difference in PMA emissions. height (and relatively constant wind field), a horizontal gradient in the aerosol concentration and size distribution at 100 m asl was not observed during this swell event. These results were predictable because the 1 m increase in wave height is relatively small and not expected to generate a significant difference in PMA emissions.  Due to the low number of TCs passing in the immediate vicinity of Réunion Island during the TC season 2018-2019, the range of wind speeds and ocean wave conditions were not as large as expected to fully assess parameterizations for PMA emissions. Nevertheless, the datasets collected during this campaign have clearly demonstrated the scientific potential of the BOREAL UAS, and of its associated payload, to collect key observations in remote oceanic areas and in the vicinity to tropical cyclones.

GNSS Observations
A well-known (indirect) atmospheric application of GNSS measurement consists of measuring the atmospheric integrated water vapor (IWV) content from the delay induced by water vapor during the crossing of the earth's atmosphere by the GNSS signal [110]. Due to the fundamental role of water vapor in climate and weather dynamics, tropospheric GNSS measurements have rapidly become one of the main tools used by climatologists to monitor the evolution of the water vapor field at all spatio-temporal scales, but also to improve NWP model forecasts (see [111] for a recent review of current GNSS weather applications). The creation of the International GNSS Service (IGS) network in 1994 [112], which now includes more than 500 stations worldwide, has also enabled GNSS-derived measurements collected in all parts of the world to be widely disseminated to the scientific community. However, with only eight operational stations available prior to RNR-CYC (Figure 21), the density of public GNSS stations in the SWIO was the lowest of all TC basins.
During RNR-CYC, 10 new permanent stations were deployed, including one in the Seychelles (Aldabra), four in the Eparses Islands (three RNR-CYC stations at Juan de Nova, Tromelin and Europa and one shared station installed by OPGP in Glorieuses), and five in Madagascar (Diego-Suarez, Toamasina, Sainte-Marie, Nosy Be and Fort-Dauphin), with the goals of increasing the density of the GNSS network in the SWIO and to provide additional near-real-time IWV measurements at various locations in the basin (Figure 21). Although these stations have not yet been included into the IGS network (pending), most collected GNSS observations are transmitted at hourly time steps to IGN's data center, which already allows for routine real-time GNSS calculations and the wide dissemination of derived meteorological and geophysical products through IGN's permanent GNSS network (RGP, ftp://rgpdata.ign.fr/pub/gnssmayotte, accessed on 20 March 2021).
These new stations have already been shown to represent a fantastic asset for monitoring the spatio-temporal evolution of the water vapor field at local and regional scales [34,35], as well as to evaluate NWP forecasts [34]. A further modelling application is shown in Figure 22, which presents differences between GNSS-derived IWV observations and daily 6 h analyses of AROME-IO operational NWP systems [42] at 00:00, 06:00, 12:00, and 18:00 UTC at Aldabra (ALBR, The Seychelles), Antananarivo (ABPO, Madagascar), Réunion Island (RUN, France) and Sainte-Marie (MASM, Madagascar) during the year 2019. In order to evaluate the performance of the model at various time scales, 30-day (red), 10-day (yellow) and 3-day (black) moving averages were applied to both the model and GNSS data. GNSS observations were processed following the approach proposed by [35]. and daily 6 h analyses of AROME-IO operational NWP systems [42] at 00:00, 06:00, 12:00, and 18:00 UTC at Aldabra (ALBR, The Seychelles), Antananarivo (ABPO, Madagascar), Réunion Island (RUN, France) and Sainte-Marie (MASM, Madagascar) during the year 2019. In order to evaluate the performance of the model at various time scales, 30-day (red), 10-day (yellow) and 3-day (black) moving averages were applied to both the model and GNSS data. GNSS observations were processed following the approach proposed by [35]. , but showed more variations in Aldabra (ALBR, between 0.6 and 1.14 kg m −2 ). Model-observation differences did not show a clear seasonal dependency, except at Aldabra where maximum errors were observed in the middle of the winter, spring, and fall seasons. In the case of this atoll, consisting of land strips of 1-3 km width encircling a lagoon of nearly 30 km × 15 km, the land is considered as submerged by the model, which does not allow taking into account the effect of vegetation and of the diurnal cycle on the atmospheric moisture content. At Aldabra, atmospheric moisture changes in the model thus mostly depend on the variation of ocean temperatures. Interestingly, the sign of the bias error is location-dependent, suggesting that the model is not affected by systematic errors. The associated standard deviation errors also tends to decrease with the length of the smoothing period, because the latter minimized the impact of short-term variations. At short time scales (three days), one can notice that the strongest discrepancies mostly occur during the TC (wet) season (days 1-120 and 330-365), in relation with the passage of tropical cyclones in the vicinity of the GNSS ground-based stations. strips of 1-3 km width encircling a lagoon of nearly 30 km × 15 km, the land is considered as submerged by the model, which does not allow taking into account the effect of vegetation and of the diurnal cycle on the atmospheric moisture content. At Aldabra, atmospheric moisture changes in the model thus mostly depend on the variation of ocean temperatures. Interestingly, the sign of the bias error is location-dependent, suggesting that the model is not affected by systematic errors. The associated standard deviation errors also tends to decrease with the length of the smoothing period, because the latter minimized the impact of short-term variations. At short time scales (three days), one can notice that the strongest discrepancies mostly occur during the TC (wet) season (days 1-120 and 330-365), in relation with the passage of tropical cyclones in the vicinity of the GNSS ground-based stations. As already pointed out by [34], who investigated the seasonal variations of AROMEanalyzed IWV errors at Diego-Suarez (DSUA), model-observation differences seem exacerbated for stations located in the vicinity of complex orography. Except for the particular case of Aldabra, the same behavior can be noticed here for the high-altitude stations of Réunion (600 m asl, REUN) and Antananarivo (2000 m asl, ABPO), which showed higher errors than for the flat island of Sainte-Marie (MASM). These errors might be related to the difficulties of the model to properly capture the modifications of air masses caused by the orography (lifting and subsiding motions) at the local scale due to an insufficient resolution of its topography. A possible way to correct for such errors could be to assimilate GNSS observations (zenithal delays) into the model. While the operational version of the AROME-IO NWP system does not allow for data assimilation, its research configuration includes a 3D-Var scheme that can be used to perform assimilation experiments [34]. Some studies are currently ongoing to determine whether the assimilation of observations collected from these new GNSS stations can help to reduce this moisture bias. As already pointed out by [34], who investigated the seasonal variations of AROMEanalyzed IWV errors at Diego-Suarez (DSUA), model-observation differences seem exacerbated for stations located in the vicinity of complex orography. Except for the particular case of Aldabra, the same behavior can be noticed here for the high-altitude stations of Réunion (600 m asl, REUN) and Antananarivo (2000 m asl, ABPO), which showed higher errors than for the flat island of Sainte-Marie (MASM). These errors might be related to the difficulties of the model to properly capture the modifications of air masses caused by the orography (lifting and subsiding motions) at the local scale due to an insufficient resolution of its topography. A possible way to correct for such errors could be to assimilate GNSS observations (zenithal delays) into the model. While the operational version of the AROME-IO NWP system does not allow for data assimilation, its research configuration includes a 3D-Var scheme that can be used to perform assimilation experiments [34]. Some studies are currently ongoing to determine whether the assimilation of observations collected from these new GNSS stations can help to reduce this moisture bias.

Spaceborne Observations
Visible and infrared satellite observations have long been the main source of knowledge for estimating some of the parameters (e.g., radius of maximum wind, various wind radii) characterizing the structure of tropical cyclones-these parameters are, in particular, inputs of the Dvorak method [113], used by TC forecast centers to produce best-track data and issue TC advisory and forecasts. Since the launch of the first Earth observation satellite in the early 1970s, the quality of ocean surface wind estimates in the vicinity of TCs has constantly been improved.This includes, for instance, the deployment of wind scatterometers [114,115] as well as of multifrequency radiometers [116,117] that both allow for a direct, and more precise, estimate of surface winds under TCs. More recently, a new generation of microwave radiometers has also been put into operation by ESA and the National Aeronautics and Space Administration (NASA) in the frame of the Soil Moisture and Ocean Salinity (SMOS) and Soil Moisture Active Passive (SMAP) space missions, respectively. Thanks to their large coverage and revisit times, these new sensors, which allow for high-resolution (~40 km) surface wind speed sampling in all weather conditions, provide a unique description of TC structure during its whole lifetime [118], and have thus become a key source of information for TC forecasting centers.
In addition to these more or less conventional sensors, a new approach, based on the use of spaceborne synthetic aperture radars (SARs), is also being increasingly used to probe and quantify sea surface properties under extreme wind conditions. Designed in the mid-1990s, SAR systems were initially used for land applications [119], especially for monitoring earth deformation rates ranging from a few millimeters per year (e.g., for glaciers) up to 1 m per h (e.g., earthquakes and landslides). SAR potential for marine applications was already known since the launch of the first SAR satellite [120], but has significantly emerged with the launch of the Canadian Space Agency (CSA)'s RADARSAT-1 (1995) and ESA's Envisat (2002) satellites. Thanks to their unique capability to gather very high-resolution (up to a few meters) surface wind and roughness data in swaths of several hundred kilometers, spaceborne SARs have become key observing systems for monitoring, forecasting, and investigating the properties and evolution of TCs [120,121].
The deployment of a new generation of SAR systems equipped with polarization diversity on-board CSA's RADARSAT-2 (RS2, 2007) and ESA's Sentinel-1 (S1, 2014) satellites has allowed further improvements in the capabilities of these instruments to accurately map the variations of ocean surface winds in TC eyes and eyewalls [122,123]. Acquisitions made by ESA's satellites had, however, never been used to probe tropical cyclones until the implementation of the Satellite Hurricane Observations Campaign (SHOC, 2016-2017), in the Pacific and Atlantic basins. The intercomparison of aircraft reconnaissance wind data against S1 SAR measurements collected in CAT-5 hurricane IRMA (2017) during this experiment definitely demonstrated the capability of these instruments to thoroughly describe TC ocean boundary layer structures at high spatial resolution [124].
The extension of the SHOC program to the SWIO basin (referred to as SHOC-V2) was initiated in 2017, under the frame of RNR-CYC. SAR missions cannot continuously acquire wide swaths of data in high-resolution modes; therefore, a dedicated acquisition procedure was set up to collect data without impacting the operational duty cycle of the satellites. Acquisition requests were passed along to both MDA (MacDonald Dettwiler and Associates, the private company that owns the satellite RS2) and ESA S1 mission planner portals on a 24-to-48 h notice, based on satellite orbit and five-day track forecasts provided by RSMC La Réunion. In order to reduce the workload of participating space agencies, which contributed to this experimental program on a voluntary basis, acquisition requests were also generally limited to storms presenting a threat to SWIO populations.
Between 2017 and 2021, about 150 SAR images were acquired by S1A, S1B and RS2 satellites, allowing the sampling of 20 tropical storms and cyclones over four TC seasons (raw S1 data are available at https://scihub.copernicus.eu, accessed on 20 March 2021). Among these acquisitions, nearly 40% were directly obtained within the eye or the eyewall of the systems ( Table 2, "hits"). SAR images were collected throughout the SWIO basin ( Figure 23) at various storm evolution stages (cyclogenesis, intensification, dissipation). Collected wind data were used for nowcasting and best-track data reanalysis purposes at RSMC La Réunion, model verification [33] and assimilation in NWP systems [125] (this Special Issue).   Table 2. Figure 23. Location of Sentinel-1 (A and B) and RADARSAT-2 swaths listed in Table 2.
Examples of SAR-derived surface wind fields obtained through the application of the retrieval algorithm proposed by [126] are presented hereafter. The first example was taken within TC JOANINHA, which developed in the northeastern part of the Mascarene Archipelago ( Figure 17a) and reached its LMI in the vicinity of Rodrigues Island on 26 March 2019, with 10 min of maximum sustained winds of nearly 60 m s −1 (Figure 24a). The core structure of this system was observed seven times by the SAR systems deployed onboard satellites S1A/S1B and RS2 (Table 2) at various stages of its life cycle. Of particular interest is the image taken by RS2 on 28 March 2019 at 00:52 UTC (Figure 24b), when the system started to experience an eyewall replacement cycle (ERC). The latter, which often occurs when a TC reaches an intensity of 50 m s −1 , has long been recognized as one of the main mechanisms for intense TCs to intensify further [127,128]. During an ERC, outer rainbands strengthen and organize themselves into a ring of thunderstorms that progressively encircle the TC eyewall. The formation of this outer ring, usually referred to as the outer eyewall, eventually cuts off the supply of moisture and angular momentum that maintain convection around the eye, causing a weakening of the system and the eventual dissipation of the inner eyewall. The former ring is then replaced by the outer ring, which gradually contracts and intensifies, often resulting in a more intense TC than prior to the ERC.
The analysis of the SAR-derived surface wind field shown in Figure 24b clearly indicates the presence of two concentric areas of maximum wind located both around the eye and in the east-to-southeastern quadrant of the system core. The inner structure (maximum wind speed of~50.2 m s −1 ) corresponds to the eyewall of the storm, and the outer one (maximum wind speed of~45 m s −1 ) to the cyclonic rainbands progressively wrapping around the eyewall that will eventually replace the inner wall. As mentioned previously, the completion of an ERC is often followed by a re-intensification of the system. However, in the present case, such re-intensification did not occur because TC JOANINHA rapidly encountered a strongly sheared environment that caused its dissipation a couple of days later. The second example is taken from TC IDAÏ (Figure 25), which is considered as the worst natural disaster that has ever affected Mozambique (as well as the surrounding countries of Zimbabwe and Malawi), and the deadliest storm ever recorded in the SWIO. This exceptional TC initiated as a tropical depression along the northern coast of Mozambique on 3 March 2019 (Figure 25a) that initially moved inland in the northwestern part of the country for a few days, with peak winds in the order of 10-15 m s −1 . On 7 March, it made a half-turn near the Mozambique-Malawi border and moved back towards the ocean. After entering the Mozambique Channel, on 9 March, the storm experienced a rapid intensification to reach intense TC intensity (10 min maximum sustained wind >46 m s −1 ), with winds gusts estimated at up to 70 m s −1 on 11 March. IDAÏ then reversed its track back for the second time in the immediate vicinity of the Eparses island of Juan de Nova (located approximately 150 km off the west coast of Madagascar) and began its southwestward progression towards Mozambique. Shortly after this half-turn, it entered into a slight weakening trend following the beginning of its ERC. Right after its completion, IDAÏ immediately re-intensified and reached its LMI on 14 March, with a minimum The second example is taken from TC IDAÏ (Figure 25), which is considered as the worst natural disaster that has ever affected Mozambique (as well as the surrounding countries of Zimbabwe and Malawi), and the deadliest storm ever recorded in the SWIO. This exceptional TC initiated as a tropical depression along the northern coast of Mozambique on 3 March 2019 (Figure 25a) that initially moved inland in the northwestern part of the country for a few days, with peak winds in the order of 10-15 m s −1 . On 7 March, it made a half-turn near the Mozambique-Malawi border and moved back towards the ocean. After entering the Mozambique Channel, on 9 March, the storm experienced a rapid intensification to reach intense TC intensity (10 min maximum sustained wind >46 m s −1 ), with winds gusts estimated at up to 70 m s −1 on 11 March. IDAÏ then reversed its track back for the second time in the immediate vicinity of the Eparses island of Juan de Nova (located approximately 150 km off the west coast of Madagascar) and began its southwestward progression towards Mozambique. Shortly after this half-turn, it entered into a slight weakening trend following the beginning of its ERC. Right after its completion, IDAÏ immediately re-intensified and reached its LMI on 14 March, with a minimum central pressure of 940 hPa and (estimated) 10 min maximum sustained winds of~55 m s −1 . The storm then gradually weakened while progressing towards Beira, where it made landfall on 15 March at the TC stage. Numerous acquisitions were made in this system by S1 satellites, including three images collected directly in the core of the storm (Figure 25b). Two of these hits were made over the ocean by S1A, at the stage of intense TC (11 March at 02:43 UTC and 14 March at 03:08 UTC), while a third one was made at landfall by S1B on 14 March at 16:00 UTC (because 90% of this swath occurred over land, this third acquisition is not exploitable). SAR wind data shown in Figure 25 are the only high-resolution wind observations collected during the oceanic phase of this system-WMO's JDN weather station broke down a week before the passage of the storm over the island and the weather station deployed at Beira in the frame of RNR-CYC was lost during landfall-therefore, these observations represent an invaluable asset for accurately assessing numerical model simulations of this storm [33].
SAR data collected in IDAÏ and other tropical cyclones can also meaningfully complement best-track (BT) data produced by RSMC La Réunion. According to RSMC La Réunion best-track, TC IDAÏ reached its LMI on 14 March 2019 at 00:00 UTC with maximum winds of nearly 55 m s −1 . This value is in good agreement with the maximum intensity measured by S1B a few hours later (~52 m s −1 ). SAR data, however, show that these strong wind values were only observed in the eastern quadrant of the storm's core, but that the average wind speed throughout the eyewall was significantly less intense and mostly comprised between 42 and 45 m s −1 . On 11 March, the agreement between best-track and SAR-derived wind data was less good, with maximum wind speed in the order of 50 m s −1 for the BT and 60 m s −1 for SAR observations. According to SAR observations, the strongest wind speed values were also more or less uniformly distributed throughout the eyewall. Hence, the comparison of the two SAR images suggests that the destructive po- Numerous acquisitions were made in this system by S1 satellites, including three images collected directly in the core of the storm (Figure 25b). Two of these hits were made over the ocean by S1A, at the stage of intense TC (11 March at 02:43 UTC and 14 March at 03:08 UTC), while a third one was made at landfall by S1B on 14 March at 16:00 UTC (because 90% of this swath occurred over land, this third acquisition is not exploitable). SAR wind data shown in Figure 25 are the only high-resolution wind observations collected during the oceanic phase of this system-WMO's JDN weather station broke down a week before the passage of the storm over the island and the weather station deployed at Beira in the frame of RNR-CYC was lost during landfall-therefore, these observations represent an invaluable asset for accurately assessing numerical model simulations of this storm [33].
SAR data collected in IDAÏ and other tropical cyclones can also meaningfully complement best-track (BT) data produced by RSMC La Réunion. According to RSMC La Réunion best-track, TC IDAÏ reached its LMI on 14 March 2019 at 00:00 UTC with maximum winds of nearly 55 m s −1 . This value is in good agreement with the maximum intensity measured by S1B a few hours later (~52 m s −1 ). SAR data, however, show that these strong wind values were only observed in the eastern quadrant of the storm's core, but that the average wind speed throughout the eyewall was significantly less intense and mostly comprised between 42 and 45 m s −1 . On 11 March, the agreement between best-track and SAR-derived wind data was less good, with maximum wind speed in the order of 50 m s −1 for the BT and 60 m s −1 for SAR observations. According to SAR observations, the strongest wind speed values were also more or less uniformly distributed throughout the eyewall. Hence, the comparison of the two SAR images suggests that the destructive potential, and overall global intensity of the system, was probably much greater on 11 March, despite a 10% lower maximum wind speed (50 m s −1 vs. 55 m s −1 ) in BT data.

Conclusions and Perspectives
ReNovRisk-Cyclone (RNR-CYC) is an ambitious international research program aimed at developing regional cooperation in the SWIO with emphasis on the observation and modeling of tropical cyclones at both current and future time horizons. The observation component of RNR-CYC, presented in detail in this paper, has enabled the collection of numerous innovative oceanic (gliders, sea-turtle borne and seismometer-derived data) and atmospheric (UAS and SAR-derived wind data) measurements, together with more conventional observations (GNSS-derived IWV, atmospheric radiosoundings, ADCP, and wave gauge swell measurements) to investigate tropical cyclones and their environment from the local to the basin scale. Its mesoscale modelling and climate components, presented in more detail in the companion paper [33], also allowed for the development of innovative modelling systems, while its outreach component significantly increased regional collaboration between SWIO countries affected by TC hazards. The innovative data collected in the SWIO during RNR-CYC have also proven particularly useful to evaluate the performance of both research and operational TC forecasting models in this undersampled region of the world. SAR measurements, which allow for accurate evaluations of TC forecasts, can also be assimilated into NWP systems [126]; in situ and remote swell measurements can be relied upon to evaluate wave model performance; while biologging measurements are priceless to assess ocean model representation of the ocean mixed layer in cyclonic conditions.
The promising results obtained during the 3.  Tanikely), Seychelles (Aldabra), Mozambique (Ponto di Ouro), and South Africa (iSimanliso). This transdisciplinary project, constructed in close cooperation with ST specialists and oceanographers, will extend the research initiated during RNR-CYC to: (i) investigate Indian Ocean properties and spatio-temporal variability; and (ii) improve knowledge of the ecology of the five species of sea turtles living in the Indian Ocean. A particular emphasis will be put on the observation of mesoscale eddies and coastal currents that develop in the Mozambique Channel (to assess their impact on water mass distribution, transport, and mixing, as well as their role in the dynamics of the Greater Agulhas current system), and on the assimilation of collected ST-borne data in global and limited-area configurations of the ocean model NEMO [129]. This new component of the global STORM program, which is now strongly supported by the EU, UNESCO (through the GOOS/AniBOS), CNES and CNRS, among others, will be conducted for the many years to come, and is likely to be extended to other ocean basins in the mid-term.  Figure 21). Some existing, or soon to be deployed, GNSS stations will also be upgraded with co-located oceanographic sensors to monitor sea level on a regional scale. In the mid-term, these instruments will be integrated into the Global Sea Level Observing System (GLOSS, https://www.gloss-sealevel.org/, accessed on 20 March 2021) network, which monitors sea level rise on a global scale, and into the IGS network, which disseminates GNSS data to the worldwide scientific community. ESPOIRS also includes an ambitious local component aimed at collecting wind and precipitation measurements in tropical cyclones with a transportable polarimetric weather radar to be deployed in Réunion Island (2021), the Seychelles (Mahe, 2022), and Madagascar (Tamatave or Diego Suarez, 2023). This weather radar component will notably focus on the impact of orography on tropical cyclones (or outer cyclonic rainbands in the Seychelles) [130], and will help to develop regional cooperation in the field of atmospheric remote sensing.
The goal of the MAP-IO program (started in January 2021) is to deploy an atmospheric and marine biology observatory in the Indian and Southern Oceans onboard the French vessel Marion Dufresne (Figure 27), operated by the TAAF administration (for its logistical needs in the Eparses Island and French Austral Territories) and IFREMER (for scientific campaigns in the IO). With the deployment of nearly 20 new sensors on-board this RV (e.g., cytometer, titrator, thermosalinograph, cytometer, titrator, NOx, CO, CO2, O3 and CH4 analyzers, GNSS, UV radiometer, CCN, photometer, weather station), observations at the OA interface and within the atmospheric column will be collected over the long term in this region particularly sensitive to the impact of climate change (note that all observations will be transmitted in near-real-time through the RV onboard satellite communication system). The data collected in the Southern Ocean for sea and wind conditions close to those encountered in tropical cyclones will also be used to carry on with the evaluation of OA exchange parameterizations (turbulent fluxes, aerosol, and sea-spray emissions) in extreme conditions [33]. This exploratory research project, which benefits  Figure 21). Some existing, or soon to be deployed, GNSS stations will also be upgraded with co-located oceanographic sensors to monitor sea level on a regional scale. In the mid-term, these instruments will be integrated into the Global Sea Level Observing System (GLOSS, https://www.gloss-sealevel.org/, accessed on 20 March 2021) network, which monitors sea level rise on a global scale, and into the IGS network, which disseminates GNSS data to the worldwide scientific community. ESPOIRS also includes an ambitious local component aimed at collecting wind and precipitation measurements in tropical cyclones with a transportable polarimetric weather radar to be deployed in Réunion Island (2021), the Seychelles (Mahe, 2022), and Madagascar (Tamatave or Diego Suarez, 2023). This weather radar component will notably focus on the impact of orography on tropical cyclones (or outer cyclonic rainbands in the Seychelles) [130], and will help to develop regional cooperation in the field of atmospheric remote sensing.
The goal of the MAP-IO program (started in January 2021) is to deploy an atmospheric and marine biology observatory in the Indian and Southern Oceans onboard the French vessel Marion Dufresne (Figure 27), operated by the TAAF administration (for its logistical needs in the Eparses Island and French Austral Territories) and IFREMER (for scientific campaigns in the IO). With the deployment of nearly 20 new sensors on-board this RV (e.g., cytometer, titrator, thermosalinograph, cytometer, titrator, NO x , CO, CO 2 , O 3 and CH 4 analyzers, GNSS, UV radiometer, CCN, photometer, weather station), observations at the OA interface and within the atmospheric column will be collected over the long term in this region particularly sensitive to the impact of climate change (note that all observations will be transmitted in near-real-time through the RV onboard satellite communication system). The data collected in the Southern Ocean for sea and wind conditions close to those encountered in tropical cyclones will also be used to carry on with the evaluation of OA exchange parameterizations (turbulent fluxes, aerosol, and sea-spray emissions) in extreme conditions [33]. This exploratory research project, which benefits from strong support from the EU, IFREMER, CNRS, TAAF and Réunion University, among others, is likely to be maintained in the long term.
Atmosphere 2021, 12, x FOR PEER REVIEW 37 of 43 from strong support from the EU, IFREMER, CNRS, TAAF and Réunion University, among others, is likely to be maintained in the long term. Funding: This research was principally funded by the European Union, the Regional Council of Réunion Island, and the French State through INTERREG-V Océan Indien 2014-2020 projects "ReN-ovRisk-Cyclones and Climate Change" (TF) and "ReNovRisk Cyclones and Precipitation" (TF), as well as by Université de La Réunion and Centre National de la Recherche Scientifique (CNRS) in the frame of the "Contrat Plan Etat Région" (CPER) "ReNovRisk". The gliders were provided, operated, and co-funded by CNRS DT-INSU (Direction Technique-Institut National des Sciences de l'Univers). The seismic stations deployed in the "Rivière des Pluies" network (code ZF) were provided by the INSU-RESIF/SISMOB instrumental pool. The provision and operating of the BOREAL UAS were co-funded by Météo-France. The installation of GNSS stations on Eparses islands was funded by the Eparses consortium of the Terres Australes et Antarctiques Françaises (TAAF) administration under the frame of project "IOGA4MET-EI". Part of the STORM operations were also co-funded by CNRS-INSU under the frame of the LEFE project "PreSTORM". Access to RADAR-SAT-2 data was supported by Institut Mines Télécom and by public funds (Ministère de l'Education Nationale, de l'Enseignement Supérieur et de la Recherche, FEDER, Région Bretagne, Conseil Général du Finistère, Brest Métropole) received in the framework of the VIGISAT project managed by Groupement Bretagne Télédétection (GIS BreTel -Brittany Remote Sensing).
Institutional Review Board Statement: All sea turtles equipped in this programme were handled by qualified and officially accredited personnel from Kelonia's care center. The used tags meet all the requirements of the international conventions for the protection of sea turtles and were directly purchased from manufacturers specializing in marine biology and biologging. Funding: This research was principally funded by the European Union, the Regional Council of Réunion Island, and the French State through INTERREG-V Océan Indien 2014-2020 projects "ReNovRisk-Cyclones and Climate Change" (TF) and "ReNovRisk Cyclones and Precipitation" (TF), as well as by Université de La Réunion and Centre National de la Recherche Scientifique (CNRS) in the frame of the "Contrat Plan Etat Région" (CPER) "ReNovRisk". The gliders were provided, operated, and co-funded by CNRS DT-INSU (Direction Technique-Institut National des Sciences de l'Univers). The seismic stations deployed in the "Rivière des Pluies" network (code ZF) were provided by the INSU-RESIF/SISMOB instrumental pool. The provision and operating of the BOREAL UAS were co-funded by Météo-France. The installation of GNSS stations on Eparses islands was funded by the Eparses consortium of the Terres Australes et Antarctiques Françaises (TAAF) administration under the frame of project "IOGA4MET-EI". Part of the STORM operations were also co-funded by CNRS-INSU under the frame of the LEFE project "PreSTORM". Access to RADARSAT-2 data was supported by Institut Mines Télécom and by public funds (Ministère de l'Education Nationale, de l'Enseignement Supérieur et de la Recherche, FEDER, Région Bretagne, Conseil Général du Finistère, Brest Métropole) received in the framework of the VIGISAT project managed by Groupement Bretagne Télédétection (GIS BreTel -Brittany Remote Sensing).
Institutional Review Board Statement: All sea turtles equipped in this programme were handled by qualified and officially accredited personnel from Kelonia's care center. The used tags meet all