The radar network deployed in southern France during the first special observing period (SOP 1) of the Hydrological Cycle in the Mediterranean Experiment (HyMeX) was designed to precisely document the 3D structure of moist upstream flow impinging on complex terrain as a function of time, height, and along-barrier distance, and surface rainfall patterns associated with orographic precipitation events. This deployment represents one of the most ambitious field experiments yet, endeavoring to collect high-quality observations of thunderstorms and precipitation systems developing over and in the vicinity of a major mountain chain.
Radar observations collected during HyMeX represent a valuable, and potentially unique, dataset that will be used to improve our knowledge of physical processes at play within coastal orographic heavy precipitating systems and to develop, and evaluate, novel radar-based products for research and operational activities. This article provides a concise description of this radar network and discusses innovative research ideas based upon preliminary analyses of radar observations collected during this field project with emphasis on the synergetic use of dual-polarimetric radar measurements collected at multiple frequencies.
An ambitious radar deployment to collect high-quality observations of heavy precipitation systems developing over and in the vicinity of a coastal mountain chain is discussed.
The accurate prediction of orographic convective precipitation is a major meteorological challenge that depends on a wide range of time and space scales as well as complex processes ranging from moist orographic airflow dynamics to cloud microphysics. Forecasting the location and amount of heavy precipitation is particularly critical in coastal mountainous regions. Heavy precipitation can indeed generate rapid and destructive floods, and thus it represents a major threat to lives and infrastructure, especially since communities living in foothill regions have experienced large population increases. Several important field and research programs, such as the Mesoscale Alpine Programme (MAP; Bougeault et al. 2001), the Intermountain Precipitation Experiment (IPEX; Schultz et al. 2002), and the Improvement of Microphysical Parameterization through Observational Verification Experiment (IMPROVE; Stoelinga et al. 2003), among others, have been designed in the last decade to improve the understanding and prediction of orographically generated precipitation. Although the analysis of datasets collected during these field campaigns has dramatically advanced our understanding of moist dynamics and microphysics over complex terrain (e.g., Medina et al. 2005; Rotunno and Houze 2007), our knowledge of coastal orographic precipitation remains insufficient to satisfy societal demands for precise forecasts and warnings.
A fundamental element that has been missing in past field experiments is nearly simultaneous measurements of the moist marine inflow as a function of time, height, and along-barrier distance with measurements of precipitation over orography in addition to turbulence/microphysical transformations occurring in between. The Hydrological Cycle in the Mediterranean Experiment (HyMeX) field campaign (Ducrocq et al. 2014), which is dedicated to the study of the water cycle in the northwestern Mediterranean basin, intends to achieve this goal and to gather such basic but hard-to-obtain information.
High-resolution dual-Doppler and dual-polarization radar measurements are key to resolving small-scale precipitation and flow structures and to studying the role of complex topography and microphysical processes in precipitation formation and enhancement. During the first special observing period (SOP 1; fall 2012) of HyMeX, which principally focused on investigating the influence of mountain topography on moist dynamics and cloud microphysics, important efforts have thus been made to achieve comprehensive weather radar coverage. In southern France, the HyMeX radar network was specifically designed to collect high-quality observations of thunderstorms and precipitation systems developing over and in the vicinity of the Massif Central mountain chain to document the characteristics of the hydrometeors forming within the moist upstream flow impinging on the Cévennes area.
Observations collected in this framework will provide for the investigation of the microphysics and dynamics of mesoscale convective systems (MCSs) developing in this flood-prone area, and will be used to develop and evaluate innovative radar-based products for numerical weather prediction (NWP) and hydrological applications. A description of this observing network is given herein together with a discussion of potential research studies inferred from the preliminary analysis of radar data collected during HyMeX SOP 1.
THE HYMEX EXPERIMENTAL RADAR NETWORK IN SOUTHERN FRANCE.
The precipitation observing network deployed in southern France during HyMeX (Fig. 1) consists of a combination of disdrometers and weather radars operating at S, C, X, Ku, and W bands (a detailed description of these instruments is in Table ES1, available online at http://dx.doi.org/10.1175/BAMS-D-13-00076.2). All these sensors were deployed within a 90,000-km2 domain encompassing the Rhone Valley and extending northward from the Mediterranean Sea to the Massif Central and eastward from the Massif Central to the French Alps.
The core of the radar network is based on a subsample of the French operational radar network Application Radar à la Météorologie Infrasynoptique (ARAMIS; Fig. 1a). It comprises six radars operating at C and S bands, which cover the French coastal Mediterranean region with an average radar baseline of 150 km, and two X-band polarimetric radars deployed in the southern Alps under the auspices of the Risques Hydrométéorologiques en Territoires de Montagnes et Méditerranéens (RHYTMME) project (Beck and Bousquet 2013). During the experimental phase, additional dual-Doppler and dual-polarization radar coverage was also required to more explicitly resolve small-scale microphysical and dynamical processes in complex terrain as well as to collect high-resolution data in areas not well covered by the operational radars. Four transportable X-band scanning weather radars were thus deployed at various locations of the Cévennes–Vivarais area (Fig. 1b) to obtain hydrometeor type, reflectivity, and Doppler velocity measurements in both volume and range–height indicator (RHI) scanning mode.
The NOXP polarimetric Doppler radar, operated by the U.S. National Severe Storms Laboratory (NSSL), was deployed near the foothills of the Massif Central on the top of Mount Bouquet (600 m MSL). This site, located at an equal distance (∼40 km) between the Bollène and Nîmes operational S-band radars, was chosen to allow for high-resolution multiple-Doppler analysis of radar data over an area that is particularly prone to extreme flash floods (Delrieu et al. 2005).
Another mobile Doppler polarimetric radar, MXPOL, was deployed farther north in the vicinity of a high-resolution Hpiconet network (a dense network of rain gauges and disdrometers over a watershed, indicated by the red rectangle in Fig. 1b). This radar, operated by École Polytechnique Fédérale de Lausanne (EPFL, Switzerland), focused on obtaining microphysical characteristics of clouds and precipitation in terms of drop size and shape distribution with an emphasis on shallow orographic rainbands that develop over the Cévennes.
Two fast-scanning “conventional” radars (X1 and X2) operated by the French Laboratoire de Météorologie Physique (LaMP) were also deployed in the northern part of the domain to study the variability of rain at the precipitation cell scale and at high temporal resolution (<1 min). These two radars were positioned at La Bombine, a site located on the Cévennes ridge at about 1000 m MSL and approximately 10 km away from the observed maximum climatological rainfall, and at Le Chade, within the Hpiconet area, respectively.
The radar network was completed by radars profilers deployed throughout the area to investigate space–time variability of precipitation as a function of along-barrier distance. These include an ensemble of six Micro Rain Radars (MRR; Löffler-Mang et al. 1999) extending from the Mediterranean coast to the Massif Central, to collect high-resolution microphysical and kinematic measurements along a coarse cross-barrier transect on the windward slope (southeastern foothills of the Massif Central), and the S-band, dual-polarimetric Transportable Atmospheric Radar (TARA; Unal 2009; Dufournet and Russchenberg 2011), operated by the Delft University of Technology.
Finally, a network of disdrometers composed of 27 Particle Size and Velocity (PARSIVEL; Löffler-Mang and Joss 2000) disdrometers and two two-dimensional video disdrometers (2DVD; Kruger and Krajewski 2002) was also installed to complement the already dense surface observation network available over the area (Delrieu 2003). The 2DVD provides information about the shape of hydrometeors through recording two orthogonal side views of every particle falling through a 10 × 10 cm2 sampling area. The PARSIVEL is an optical disdrometer (OD) with a sampling area of ∼50 cm2, which is commonly used to investigate the microstructure of rainfall. About half of the disdrometers were concentrated in the northern part of the experimental domain, within the ∼35-km2 Hpiconet area. The remaining disdrometers were deployed at various locations along the MRR transect.
This experimental setup was also complemented by the 95-GHz airborne Radar System Airborne (RASTA) cloud radar (Protat et al. 2009) on board the French Falcon 20 (F20) aircraft. About half of the 20 missions flown by this aircraft were conducted over southern France, allowing for gathering information about the physical and radiative properties of ice particles.
Examples of data collected by these various sensors in a variety of weather situations, including a long-lasting bow echo system (IOP 6, 24 September 2012), a frontal system associated with embedded convection (IOP 8, 29 September 2012), and a fast-moving isolated thunderstorm (IOP 16, 21 October 2012), are presented herein.
HIGH-RESOLUTION OBSERVATIONS OF CLOUD DYNAMICS AND MICROPHYSICS.
An important asset of HyMeX is the ability to gather observations at the microphysical scale from the TARA and RASTA high-resolution radar profilers. These instruments are capable of observing the microphysical properties of precipitation at very high spatial and temporal resolutions—that is, 30 m in range every 3 s for TARA and 120 m every 1.5 s for RASTA—and thus are particularly suited to complement storm-scale and mesoscale observations provided by scanning weather radars. Both instruments have multibeam capability.
TARA (Unal et al. 2012) is an S-band (3.298 GHz) ground-based profiling radar with flexible antenna elevation, based on the frequency-modulated continuous wave (FMCW) principle. The antenna system includes three feeds: a dual-polarized on-focus feed, which is directed along the axis of symmetry of the parabolic reflector (i.e., 45°), and two offset feeds that produce single-polarized beams at 15° angles to the main beam to achieve three-dimensional wind speed measurements (additional information about TARA can be found online in Table ES2 at http://dx.doi.org/10.1175/BAMS-D-13-00076.2).
An example of TARA capabilities is provided in Fig. 2 using data collected during the IOP 8 frontal precipitation event (29 September 2012). The time series of observed reflectivity (Fig. 2a), horizontal wind speed (Fig. 2b), and wind direction (Fig. 2c) allow the evolution in precipitation processes to be observed for a 5-h period (1000–1500 UTC) and provide the ability to highlight differences between precipitation and cloud areas during the passage of a frontal system. TARA observations also provide insights into microphysical processes active within clouds, as illustrated by the differential reflectivity (ZDR) spectrogram shown in Fig. 2d. In this figure, ZDR values (related to particle shape) at each height bin are expressed in terms of Doppler velocity, which provides information about the particle fall velocity and, to a further extent, the particle size and density distribution. From this particular example, one can, for instance, notice 1) an increase of ice particle size in the cloud region (related to the increase of Doppler velocity values), 2) a slight decrease of raindrops size while following toward the ground, and 3) a change in shape of the raindrops from spherical to oblate as the Doppler velocity (and therefore size) increases.
RASTA is a cloud radar that can be either operated from the ground or from a mobile platform. During HyMeX, it was operated on board the French Falcon 20 research aircraft (see Table ES3 in the online supplement at http://dx.doi.org/10.1175/BAMS-D-13-00076.2 for more information about this radar). The airborne version of RASTA has a unique configuration, which consists of six antennas—three looking downward and three looking upward—providing 3D observations of the dynamics and the microphysics of clouds above and below the flight track. Once the cross track, along track, and vertical wind fields are retrieved, one can estimate the terminal fall velocity (Vt), which is then combined with radar reflectivity to compute microphysical properties of observed clouds, such as ice water content (IWC), ice particle density, and effective radius (Delanoë et al. 2007). Examples of RASTA observations obtained during IOP 8 are shown in Fig. 3.
All in all, TARA can monitor the temporal evolution of cloud processes by observing the same atmospheric column at high temporal resolution, while RASTA can be used to provide accurate spatial distribution of such processes over a larger area of interest. This complementarity shall be investigated in the near future by using RASTA data collected during the numerous passages of the F20 aircraft over TARA (e.g., Fig. 3e).
RETRIEVAL OF RAINFALL DSDS.
The disdrometer network was deployed to monitor the variability of drop size distribution (DSD) at ground level as well as to provide relevant information for both polarimetric radar data interpretation and radar rainfall estimation assessment. Collected DSD data can be used to evaluate the performance of radar-based DSD retrieval algorithms as well as to establish local power laws adapted to the Cévennes area for the interpretation of radar observations. This potential is illustrated in Fig. 4, which shows the temporal evolution of the rainfall rate, the total drop concentration, and the median volume diameter as observed by disdrometers along the south–north MRR transect (Fig. 1b) and within the Hpiconet area on 29 September 2012 (IOP 8). The large distances between Candillargues, Alès, and the HPiconet area (Fig. 1b) generate a clearly visible time shift in DSD evolution. The small-scale variability of the DSD is illustrated by the range of values measured by the different disdrometers within the Hpiconet network.
EVALUATION OF HYDROMETEOR CLASSIFICATION ALGORITHMS (HCA).
A major advantage of radar polarimetry is the ability to infer the species of hydrometeors composing precipitating systems. This capability is important to attain many objectives of the HyMeX research program, such as the evaluation of new microphysical schemes developed for high-resolution NWP systems, the investigation of relationships between radar observables and variables of interest (e.g., rain intensity), and between cloud microphysics and lightning activity (see below). In this regard, a particular emphasis was put on the evaluation and improvement of hydrometeor classification algorithms at X,-C-, and S-band frequencies.
Figure 5 presents preliminary results of an ongoing study focusing on evaluating the consistency of microphysical retrievals at C band (Montclar radar) and S band (Nimes radar) by comparing microphysical retrievals within a common sampling area of the two radars. The fuzzy logic algorithm used in this study is based upon the approaches proposed by Al-Sakka et al. (2013), Dolan et al. (2013), Park et al. (2009), and Marzano et al. (2006). It discriminates between six hydrometeor species (rain, wet snow, dry snow, ice, graupel, and hail) at S, C, and X bands. Overall, the fraction of hydrometeor types inferred from the global analysis of 4 h of radar data collected during IOP 8 (Fig. 5b) show good consistency despite the different wavelengths and slightly different beam properties of the two radars (1.1° vs 1.25° beamwidth).The time series of hydrometeor species retrieved from the two radars by steps of 5 min (Figs. 5c,d) are also generally in good agreement, although some dissimilarity could be observed around 1400 UTC. Note that a complete volume scan is composed of three different 5-min cycles that are repeated every 15 min. On the Nimes and Montclar radars, only the highest elevation plan position indicators (PPIs) go up high enough to sample ice over the intercomparison area. These PPIs are only performed during the first 5-min cycle, which explained why ice is only detected every 15 min.
To better quantify the performance of this, or other, HCA, statistical studies and systematic comparisons with in situ [Forward Scattering Spectrometer Probe 100 (FSSP-100), Precipitation Imaging Probe (PIP), 2D stereo (2-DS) probe] data will be conducted using flight-level data collected during all F20 missions flown over France. Detailed radar intercomparisons will also be reinforced and expanded by taking the research polarimetric X-band radars National Oceanic and Atmospheric Administration (NOAA) X-Pol (NOXP) and MXPOL into account in the analysis.
REAL-TIME 3D WIND RETRIEVAL.
The ability to perform multiple-Doppler wind synthesis from operational weather radar systems on a real-time basis was investigated by the French weather service Météo-France in 2006 using a network of six Doppler radars covering the greater Paris area (Bousquet et al. 2007, 2008b). In preparation for the field phase of HyMeX, this analysis was tested over regions of complex terrain of southern France (Beck and Bousquet 2013) before being eventually extended to the entire French radar network, so as to produce a nationwide 3D reflectivity and wind field mosaic (Bousquet and Tabary 2014). During the HyMeX field phase, a special version of this radar mosaic was used to guide the Falcon 20 research aircraft over the Cévennes and Rhone Valley areas. Multiple-Doppler wind (MDW) and reflectivity fields were produced in real time every 15 min over a domain of 200 km × 200 km using the C and S band operational radars shown in Fig. 1a.
Figure 6 presents examples of real-time, high-resolution (1 km2), multiple-Doppler analyses of radar data collected within the bow echo observed during HyMeX IOP 6. At the time of observations, the MCS exhibited many common characteristics of bow echo systems (Fujita 1978; Davis et al. 2004), such as a strong (24 m s–1), 5-km-deep subsiding rear-to-front (RTF) inflow jet, cyclonic and anticyclonic bookend vortices, and strong vertical motion (>12 m–1) along the leading edge, in addition to strong deep convection with reflectivity values of 25 dBZ up to 12 km MSL. In the near future, these high-quality wind fields will be further improved by incorporating research radar data collected during SOP 1 into the analysis. In addition to providing insights about the dynamics of MCSs sampled during HyMeX, MDW will also be used to assess (and eventually correct) numerical model forecasts through identifying possible temporal or spatial phase shifts in model outputs (e.g., Beck et al. 2014; Bousquet et al. 2008a).
The HyMeX field campaign also represents a unique opportunity to compare MDW fields with independent measurements, such as in situ or airborne radar data. With this regard, a comparison between MDW and RASTA winds measured during the IOP 6 flight is presented in Fig. 7. Overall, one can notice a good agreement between the two datasets, although RASTA-derived winds obviously show more details because of the superior resolution of airborne measurements (Δh = 500 m, Δz = 120 m) with respect to ground-based observations (Δh = 1000 m, Δz = 500 m). Such comparisons between MDW- and RASTA-derived winds will be performed for all RASTA flights conducted over France, so as to evaluate ground-based and airborne wind retrievals over the plain and high terrain.
SYNERGY BETWEEN RADAR AND LIGHTNING OBSERVATIONS.
Observations collected during HyMeX also offer a unique opportunity to investigate the relationships between microphysics, dynamics, electrification, and lightning occurrence in storms that were observed in southern France. Indeed, during SOP 1, total lightning activity, including intracloud (IC) and cloud-to-ground (CG) flashes, and electrical properties of storms were recorded with a Lightning Mapping Array (LMA; Rison et al. 1999) composed of 12 stand-alone stations, several operational lightning locating systems [e.g., European Cooperation for Lightning Detection (EUCLID)], and a set of ground-based research instruments (electric field sensors, acoustics sensors, video cameras) deployed throughout the experimental area.
Different types of convection from early to decaying stages were documented over land and over sea with the LMA. An example of concurrent observations within the IOP 6 storm from 0215 to 0230 UTC is shown in Fig. 8. At this time, the lightning activity was distributed more or less along a north–south direction and extended farther north outside the LMA network. It was principally located to the west of significant updrafts, as retrieved from radar observations (Figs. 8b,c). The deepest convective cell was recorded south of the system with lightning observed up to a height of 12 km MSL (Fig. 8e). In this thunderstorm, preliminary analysis of the lightning data suggests that the IC–CG ratio was 110:14 for the period 0225–0230, where 64% of CG flashes were of negative polarity. For the studied core, electrical discharges were recorded mainly in cloud regions with reflectivity above 20 dBZ.
Another interesting example illustrating the potential of polarimetric radar observations for lightning studies and vice versa was sampled by NOXP on 21 October 2012 (Fig. 9). At the time of observations, the radar data indicated an intense cell characterized by significant signal loss in the liquid phase of precipitation at lower elevation angles. Above the melting layer, ZDR presents an interesting positive–negative couplet associated with a decreasing differential phase (ΦDP) in the radial direction. This signature is the result of cross coupling between orthogonally polarized waves when the radar is operating in simultaneous transmission and reception (STAR) mode. Radial streaks in ZDR and ΦDP at heights of 7–10 km have been attributed to ice crystals oriented by electrostatic fields in regions of storm electrification (Ryzhkov and Zrnić 2007; Hubbert et al. 2010). While these depolarization signatures, which can cause bias in ZDR, are usually considered undesirable (Hubbert et al. 2010), they can also be used as an indicator of strong electrification and thus be used as an opportunity to detect charged ice particles in absence of a lightning detection sensor.
As it is believed that lightning frequency is proportional to the product of the downward mass flux of graupels and of the upward mass flux of ice crystals (Deierling et al. 2008), statistical analyses of LMA data combined with ice flux and mass budgets inferred from dual-polarimetric radar observations, MDW analyses, and exploitation of the vertical profile of reflectivity (VPR) polarimetric models (e.g., Kirstetter et al. 2013) will be relied upon to investigate relationships between dynamics, microphysics, and lightning activity. These results will also be used to evaluate the explicit lightning scheme Cloud Electrification and Lightning Scheme (CELLS) recently implemented in the Méso-NH research model (Barthe et al. 2012; Pinty et al. 2013).
TOWARD THE ASSIMILATION OF NEW RADAR OBSERVATIONS.
An expected outcome of the HyMeX project is to thoroughly assess the quality of dual-polarimetric radar measurements as well as to investigate new avenues for their utilization in operational weather forecasts. Specifically, this objective includes the development of data assimilation techniques to ingest dual-polarimetric radar data into convection-permitting NWP models. A necessary step for data assimilation was the development of a flexible, fully featured radar simulator of polarimetric radar variables (e.g., Jung et al. 2008; Pfeifer et al. 2008; Ryzhkov et al. 2011) within the mesoscale, nonhydrostatic atmospheric model Méso-NH (Lafore et al. 1998) to enable direct comparisons between model-simulated and observed polarimetric radar variables. The newly developed polarimetric radar simulator, which is based on the conventional radar simulator of Caumont et al. (2006), calculates electromagnetic wave propagation and scattering at S, C, and X bands and considers beam propagation effects, such as (differential) attenuation and phase shift, and beam refraction and broadening. It is fully consistent with the microphysical parameterizations of the Méso-NH model, which uses a single-moment bulk microphysical scheme governing the equations of the six following water species: vapor, cloud water, liquid water, graupel, snow, and pristine ice. It can be used to simulate all polarimetric observables, such as reflectivity at horizontal polarization, ZDR, ΦDP, and specific differential phase shift (KDP), as well as cross-correlation coefficient and differential backscattering phase, among others.
An example of observed and simulated polarimetric variables is shown in Fig. 10 for the IOP 6 bow echo system. According to observed (Figs. 10a,c,e) and simulated (Figs. 10b,d,f) polarimetric data, the MCS appears reasonably well reproduced by the model although one can note that the convective line is slightly shifted northwestward and that the simulated values are lower than their observed counterparts. These spatiotemporal lags and biases will be considered “innovations” in a data assimilation context, providing needed information to adjust model-state variables and improve forecasts at later times. Future work with these datasets will identify the most useful polarimetric observables for data assimilation purposes and will eventually establish the degree to which their assimilation improves short-term forecasts of heavy precipitation events.
Real-time radar refractivity retrieval.
Radar refractivity measurements consist of using the signal returned by ground clutter to estimate the refractive index, which is related to surface pressure, temperature, and especially humidity (Fabry et al. 1997). These observations represent a unique dataset that can be used to map the distribution of low-level moisture, a key parameter for understanding convection initiation and evolution (Besson et al. 2012; Besson and Parent du Chatelet 2013).
The HyMeX campaign was chosen as a test bed for real-time refractivity retrievals by Météo-France. The refractivity retrieval algorithms of Fabry et al. (1997) and Fabry (2004) were first adapted to the characteristics of French operational radars, all of which are equipped with magnetron transmitters (Parent du Châtelet et al. 2012) and then applied to three operational S-band radars (Nîmes, Bollène, and Opoul) located in southern France (Fig. 1a). Refractivity retrievals were evaluated through comparisons with automatic weather stations (AWS; Fig. 11) as well as NWP outputs using the observation operator described by Caumont et al. (2013). The good agreement between radar-measured, observed, and model-simulated refractivity demonstrate the ability to achieve real-time, accurate, low-level refractivity measurements from noncoherent weather radars, and it definitely paves the way for assimilating radar refractivity measurements in convection-permitting NWP models.
SUMMARY AND OUTLOOK.
The ongoing deployment of operational and research polarimetric radar systems throughout the planet allow for a whole new world of challenges to the international radar community. In this regard, the radar observations collected in southern France during HyMeX SOP 1 represent a valuable, and potentially unique, dataset that can be used to improve our knowledge of physical processes at play within coastal orographic heavy precipitating systems, but also to develop, and to evaluate, novel radar-based products for research and operational activities. For this reason, radar data collected by international teams involved in the project have all been made freely available in the HyMeX database (http://mistrals.sedoo.fr/HyMeX/) for noncommercial research and educational outreach activities. The use and dissemination of these data and products are encouraged by all radar scientists to foster collaborations and discussions between international research teams from all over the world.
The HyMeX 10-yr concerted experimental effort, which has brought together more than 400 scientists from all over the world, will last until 2020 (Drobinski et al. 2014). Although no new field phase is expected at the moment, the HyMeX database will continue to be regularly enhanced with new observations from operational weather services and local field experiments. Collaborations and discussions between scientists will also continue under the auspices of the HyMeX international workshop series, which has been organized every year since 2007.
We wish to thank all the scientists who contributed to the realization of the HyMeX radar component. A special mention is due to the members of Météo-France’s operational observation division, who managed to provide operational data to the HyMeX community—B. Fradon, H. Al-Sakka, A.-A. Boumahmoud, J. Parent-du-Chatelet, and P. Tabary—and to scientists S. Coquillat, M. Hagen, L. Labatut, Y Lemaître, and Y. Pointin and to students E. Fontaine, F. Pantillon, and T. Wisman, who participated in field operations and preliminary data analysis. HyMeX SOP 1 was supported by CNRS, Météo-France, CNES, IRSTEA, and INRA through the large interdisciplinary international program Mediterranean Integrated Studies at Regional and Local Scales (MISTRALS), which is dedicated to the understanding of the Mediterranean basin environmental process (www.mistrals-home.org). The radar component of HyMeX SOP 1 was sponsored by Grants ANR-2011-BS56-027 FLOODSCALE and ANR-11-BS56-0005 IODA-MED.
A supplement to this article is available online (10.1175/BAMS-D-13-00076.2)