## CloudSat-inferred vertical structure of precipitation over the Antarctic continent

## F. Lemonnier , J-B. Madeleine , C. Claud , C. Palerme , C. Genthon , T. 1 1 1 2 1 L'Ecuyer and N. B. Wood 3 3

1 Sorbonne Universit´, Ecole normale sup´rieure, PSL Research University, Ecole e ´ e ´ polytechnique, CNRS, Laboratoire de M´ et´orologie dynamique, LMD/IPSL, F-75005 Paris, e France 2 Development Centre for Weather Forecasting, Norwegian Meteorological Institute, Oslo, Norway 3 Department of Atmospheric and Oceanic Sciences, University of WisconsinMadison, Madison, Wisconsin, USA

## Key Points:

- · Regridding of the CloudSat observations reveals the 3D structure of precipitation over Antarctica
- · Distributions vary substantially by type of region: ice sheet, peninsula or ice shelves
- · Precipitation rate - temperature relationship is explained by large-scale orographic precipitation

## Abstract

Current global warming is causing significant changes in snowfall in polar regions, directly impacting the mass balance of the ice caps. The only water supply in Antarctica, precipitation, is poorly estimated from surface measurements. The onboard cloudprofiling radar of the CloudSat satellite provided the first real opportunity to estimate precipitation at continental scale. Based on CloudSat observations, we propose to explore the vertical structure of precipitation in Antarctica over the 2007-2010 period. A first division of this dataset following a topographical approach (continent versus peripheral regions, with a 2250 m topographical criterion) shows a high precipitation rate (275 mm.yr -1 at 1200 meters above ground level) with low relative seasonal variation ( ± 11%) over the peripheral areas. Over the plateau, the precipitation rate is low (34 mm.yr -1 at 1200 m.a.g.l.) with a much larger relative seasonal variation ( ± 143%). A second study that follows a geographical division highlights the average vertical structure of precipitation and temperature depending on the regions and their interactions with topography. In particular, over ice-shelves, we see a strong dependence of the distribution of precipitation on the sea-ice coverage. Finally, the relationship between precipitation and temperature is analyzed and compared with a simple analytical relationship. This study highlights that precipitation is largely dependent on the advection of air masses along the topographic slopes with an average vertical wind of 0.02 m.s -1 . This provides new diagnostics to evaluate climate models with a three-dimensional approach of the atmospheric structure of precipitation.

## 1 Introduction

Predicting the mass balance of ice sheets is a major challenge in the context of sea level rise. This surface mass balance depends on the relative magnitudes of precipitation versus sublimation/evaporation, meltwater run-off and blowing snow (Eisen et al., 2008). On an unfriendly continent for field observations, satellites are keys for observing precipitation. Based on these spatial remote sensing observations, recent studies observe an increasing ice loss in West Antarctica due to the acceleration of glacier and basal melting (Shepherd et al., 2018, 2012; Pritchard et al., 2012). This does not seem to be restricted to the West alone, and Rignot et al. (2019) shows that ice loss is more significant than its gain at the scale of the entire continent.

According to the climate models of the 5 th Climate Model Intercomparison Project, precipitation over Antarctic would increase from 5.5% to 24.5% between 1986-2005 and 2080-2099. This increase would have a significant impact on sea level (Church et al., 2013). However, the present-day averaged Antarctic snowfall rates predicted by these models range from 158 to 354 mm.yr -1 while the first model-independent climatology of current Antarctic precipitation yields a value of 172 mm.yr -1 over the 2007 - 2010 period (Palerme et al., 2017). The latter value was recently re-evaluated at 160 mm/year (Palerme et al., 2019) over the 2007 - 2010 period. It is worth noting that Milani et al. (2018) have also produced a climatology that also includes the Southern Ocean.

Apart from these surface precipitation estimates, the vertical structure of precipitation was poorly known until recently, but work by (Grazioli, Madeleine, et al., 2017) compared vertical year-averaged profile of precipitation observed with ground radars with simulated snowfall profiles at Dumont d'Urville station. It showed surface rates in good agreement with climate models but significant biases in the vertical structure of precipitation. Using long periods of ground-based observations to characterize local snowfall events showed also a significant amount of virgas in precipitation events (Maahn et al., 2014; Dur´n-Alarc´n et al., 2019). a o A poor representation of the vertical structure of precipitation reveals deficiencies of process representation in models including microphysical processes.

The thermal structure of the atmosphere over Antarctica, which is deeply involved in the origin of precipitation, is unique. The center of the continent being a cold pole, the structure of the Antarctic circumpolar flow is of baroclinical type. Low pressure systems are thus generated by horizontal temperature gradients in the troposphere and grow through baroclinic instabilities. Low pressure systems bring strong winds to the coastal areas when pressure gradients increase as polar lows over the ocean are moving southward and encounter areas of higher pressure, such as the semi-permanent continental anticyclone (Bromwich &amp; Parish, 1998; King &amp; Turner, 2007; Van de Berg et al., 2007). These winds are the moisture vectors that control large scale precipitation processes. Indeed, winds move air masses against the Antarctic continent and lead to encounters of moist and dry, hot and cold air masses, thereby creating ideal conditions for the formation of precipitation. Over the vertical axis, precipitation undergoes various changes before it reaches the ground. The origin of precipitation is described by cloud microphysics. Microphysical processes such as ice crystal nucleation are initiated in clouds, grow by diffusion at the expense of the supercooled droplets or by collision with other crystals, and fall to the ground (Findeisen et al., 2015; Dye et al., 1974). During its downfall, precipitation can reach drier air masses and sublimate, which can significantly reduce the precipitation rates. On the contrary, precipitation can reach saturated air masses, interact with potential clouds at lower altitude, aggregate and increase its mass flow. It is precisely all these processes and vertical evolution that define surface precipitation.

Precipitation in Antarctica is challenging to study, mostly due to geographic characteristics. Indeed, ground-based measurements are sparse and challenging in Antarctica and the size of this continent does not allow to cover and study the whole distribution, frequency and rate of precipitation. In coastal areas, it is influenced by synoptic conditions such as oceanic fronts (Bromwich, 1988) and it is difficult to distinguish between precipitation and blowing snow caused by strong katabatic winds. Furthermore, re-sublimation processes of the precipitation in the lower layer of the atmosphere have been observed at Dumont d'Urville with a micro-rain radar, and result in a decrease in the precipitation rates at the surface (Grazioli, Madeleine, et al., 2017). At Princess Elisabeth and Dumont d'Urville, the study of Dur´n-Alarc´n et al. (2019) showed that in a o about a third of precipitation event cases, they are virgas that do not reach the surface. There can also be advection of very large amounts of moisture, providing a high annual precipitation contribution over the continent (Gorodetskaya et al., 2014). According to ground observations over peripheral areas and plateau, the ice crystals that are at the origin of the precipitation are thinner than 100 µ m and aggregate to constitute coarser particles (Lachlan-Cope et al., 2001; Lawson et al., 2006). A recent study shows that the occurrence of ice crystals at Dumont d'Urville is similar to that of aggregates at low altitude ( &lt; 500 m above the surface) (Grazioli, Genthon, et al., 2017). Over the continental plateau ( &gt; 2250 m), most of the precipitation is driven by a few fronts while the remaining annual precipitation rate is in the form of 'Diamond Dust' (thin ice crystals) under clear-sky conditions (Bromwich, 1988; Fujita &amp; Abe, 2006; Turner et al., 2019).

The first real opportunity to assess precipitation in polar regions from a spaceborne radar platform appeared with the cloud-profiling radar (CPR) on CloudSat satellite (Stephens &amp; Ellis, 2008; Liu et al., 2008). CloudSat provided day- and night-time solid precipitation observations from August 2006 to April 2011, which led to the first multi-year, modelindependent climatology of Antarctic precipitation, with an average of 171 mm w.e year -1 over the Antarctic ice sheet, north of 82 o S (Palerme et al., 2014). The first hundreds of meters over the surface are contaminated by radar wave reflections from the surface, which is called ground-clutter. A study of Maahn et al. (2014) investigated whether a thinner ground clutter layer could increase the quality of the precipitation measured at Princess Elisabeth, but the data are really usable only from 1200 m.a.g.l. The climatology of Palerme et al. (2014) is based on the 2C-SNOW-PROFILE product (Wood, 2011),, which provides, on average, an uncertainty on single snowfall retrievals ranging between 1.5 and 2.5 times the snowfall rate (Wood, 2011; Palerme et al., 2014). However, recent stud-

ies based on comparisons with ground-based radars reassessed this range of uncertainty for a few snowfall events at two Antarctic stations (Lemonnier et al., 2019; Souverijns et al., 2018), and provided uncertainties ranging from 13 to 22% for the CloudSat snowfall observations. Palerme et al. (2019) have identified some CloudSat observations likely contaminated by the ground clutter over mountainous areas, and have produced a new climatology with a mean snowfall rate of 160 mm.yr -1 .

In this study, we computed the first 3-D multi-year climatology of Antarctic precipitation north of 82 ◦ S from spaceborne remote sensing observations (section 2). As a first step, we characterize the general horizontal and vertical structures of the precipitation averaged in time over the continent (section 3). Then, a comparison between CloudSat snowfall retrievals and concurrent 3D meteorological parameters obtained from the European Center for Medium-Range Weather Forecast (ECMWF) operational weather analysis at each vertical level is conducted in section 4. This allows us to determine which physical or dynamical process predominates on the high continental plateau and on the peripheral areas. At last, we propose a new way of evaluating modelled precipitation using this new dataset, and especially its vertical structure, in atmospheric models such as global GCMs or Antarctic regional models (section 5).

## 2 Data and method

Here, the climatology produced by Palerme et al. (2019) is extended to all vertical levels up to 10000 m. The first 4 bins above the surface are excluded due to potential ground clutter contamination (Wood, 2011). We used the fifth release of 2C-SNOWPROFILE CloudSat product (Wood, 2011; Wood et al., 2014) that gives an estimate of the liquid-equivalent snowfall rate and its Snow Retrieval Statue (SRS) which indicates if some errors occurred in the retrieval process. The retrieval of this product is based on a relationship between reflectivity and snowfall that takes into account a priori estimates of snow particle size distribution and microphysical and scattering properties at each vertical bin (Rodgers, 2000). Then this product is averaged within each cell of a grid over the Antarctic continent from about 1200 m.a.g.l. (meters above ground level) up to 10000 m.a.g.l. (5 th to 45 th bins above the surface) with a 240 m vertical resolution. The same extraction process is applied to the associated ECMWF-AUX operational weather analysis temperatures, humidities and pressures in order to obtain a climatology based on the same sampling.

The low temporal sampling of CloudSat is a source of uncertainty. Following Souverijns et al. (2018), this new CloudSat climatology is processed over a grid of 1 ◦ of latitude by 2 ◦ of longitude. This grid is able to accurately represent a snowfall climatology and is the best balance between satellite omissions and commissions, according to data from three ground radars.

2C-SNOW-PROFILE provides a snowfall uncertainty between 1.5 and 2.5 times the solid precipitation rate for a single measurement. The average over the whole continent over the 2007-2010 period reduces this uncertainty but some uncertainty remains due to systematic errors in algorithm assumptions that are difficult to assess (Palerme et al., 2014). Lemonnier et al. (2019) reassessed this range of uncertainties by a shorttime (a few seconds) and short-spatial scale comparison of CloudSat measurements with 24 GHz Micro Rain ground radars at different locations in East Antarctica (a coastal and a high altitude station) with two different kinds of weather conditions, leading to new values of [-13%, +22%] and providing confidence over the full height of the CloudSat retrievals.

The temperature associated with each precipitation profile is obtained from ECMWFAUX operational weather analysis. The ECMWF-AUX data set is an intermediate product that contains the set of ECMWF state variable data interpolated to each CloudSat

cloud profiling radar bin. Besides the snowfall rate climatology computed on a continental grid at each vertical level of the satellite, the pre-gridded dataset of the 2C-SNOWPROFILE product (snowfall retrieval and geolocation fields) as well as the ECMWFAUX product were extracted in order to evaluate the relationship between precipitation rates and temperatures (see section 4).

## 3 General structure of Antarctic precipitation

Figure 1 shows precipitation maps over Antarctica at four different CloudSat vertical levels. As highlighted in previous studies (Palerme et al., 2014; Milani et al., 2018), there is a dichotomy in precipitation rate between the high plateau of the continent where the snowfall rate is excessively low and the coastal areas. It is also important to note extremely high precipitation rates located west of the Antarctic peninsula and over mount Vinson. This is due to high orography and moisture flux from the Southern Ocean, producing more clouds and precipitation (Nicolas &amp; Bromwich, 2011). Precipitation over the ice-shelves appears to be very low. According to the study of Knuth et al. (2010) carried out over the Ross ice-shelf, it even appears that the blowing and drifting snows prevails over the accumulation of precipitation. Thus, precipitation is regionally very variable in Antarctica. Figure 1 shows that precipitation decreases with altitude over the continent, which is not always the case over the ocean. Indeed, over all oceanic regions south of 60 ◦ S except the eastern part of the peninsula and above the ice-shelves, the precipitation rate increases between the 1200 m.a.g.l. (fig. 1.a) and 2160 m.a.g.l. (fig. 1.b) vertical bins before decreasing further at altitude (fig. 1.c and 1.d). This could be explained by a local formation of snow in regions of shallow convection, as highlighted by Kulie et al. (2016) and Milani et al. (2018), who looked at the differences between the 3 rd and the 6 th vertical bins.

The thin black solid line in fig. 2 delineates regions that are 2250 m above sea-level. This altitude separates peripheral areas where precipitation is mostly controlled by large scale oceanic fronts, from the plateau where snowfall is rare and depends on a few events which are large enough to bring moisture onto the Antarctic plateau.

In order to have a better understanding of the horizontal and vertical structures of snowfall, the Antarctic continent is also divided into several geographic regions. Figure 2 presents these different regions for the following geographic areas :

- · A The first studied area is the East Antarctic, located between 0 ◦ and 150 ◦ of longitude. This area is defined mostly as an homogeneous ice-sheet with a huge plateau at an altitude of 3000 m above sea level.
- · B The second studied area is the West Antarctic, lower in altitude than the East part of the continent. This region has been subdivided in two parts, one comprising mountain ranges including the mount Vinson (B1) and the other part comprising the massive Ross and Ronne-Filchner ice-shelves (B2).
- · C The third region of interest is the peninsula, which is very mountainous and where precipitation is exceptionally high.

## 3.1 From coasts to high continental plateau

By dividing Antarctica according to topographical information (2250 m in surface elevation, see Palerme et al. (2014)), we distinguish two different types of climate. The peripheral areas includes the Peninsula (surface elevation &lt; 2250 m), West Antarctica and eastern coasts. The plateau (surface elevation &gt; 2250 m) includes essentially East Antarctica. Following Palerme et al. (2017), fig. 3 summarizes the seasonal evolution of the precipitation over the Antarctic continent. The seasonal variability of snowfall across the Antarctic continent is mainly influenced by snowfall in peripheral regions, with the periphery of the continent receiving the vast majority of precipitation (Bromwich, 1988;

Figure 1. Mean annual snowfall rate (mm water equivalent / year - mm.yr -1 ) from the 2CSNOW-PROFILE product over the 2007-2010 period a) at 1200 m.a.g.l. (5 th bin); b) at 2160 m.a.g.l. (8 th bin); c) at 3120 m.a.g.l. (12 th bin) and d) at 5040 m.a.g.l. (20 th bin).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00457v2-figure-6c7447ab-e38e-4e13-8fad-c37521ce2442.jpg)

Genthon et al., 2009). This figure corresponds to the 5 th CloudSat vertical bin at about 1200 m.a.g.l. Over the Antarctic continent and especially over peripheral areas (altitude &lt; 2250 m), the maximum snowfall rate is observed in March-April-May (MAM) then decreases until December-January-February (DJF) where the snowfall rate is minimum. Over the plateau, this variability differs with a maximum in precipitation rate over the DJF period and a minimum in June-July-August (JJA). Error bars represent the CPR measurement uncertainties given by Lemonnier et al. (2019), obtained by a comparison of four precipitation events with ground instruments at two different locations. The seasonal evolution of precipitation has been studied for the other vertical levels of the product, and no difference in seasonal variability from one level to another is observed. As already mentioned by Palerme et al. (2017), the seasonal evolution of mean precipitation in Antarctica is in accordance with ERA-Interim and climate models for the peripheral regions. This is not the case for precipitation above the plateau. A more recent study by Lenaerts et al. (2018) focusing on ground accumulation presents similar results, although it focuses on surface snow accumulation, where other processes such as surface sublimation or blowing snow are at work.

Fig. 4 shows another aspect of the vertical structure of precipitation in Antarctica above the ground clutter layer represented by annual average precipitation profiles. Each

Figure 2. Digital Elevation Map of the Antarctic Ice Sheet (Liu et al., 2015) with four frames corresponding to the four studied areas. A in magenta is the East Antarctic continent; B in red is the West Antarctic, which has been subdivided into continental region (B.1) and ice-shelves (B.2); C in green is the Peninsula. Thin black solid line is the iso-altitude 2250 m and separates the coastal areas from the plateau.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00457v2-figure-7e92cdac-6852-4d6a-8a3d-1424bce9516d.jpg)

Figure 3. Seasonal variability of snowfall north of 82 ◦ during the period 2007-2010 for CloudSat in mm/yr at about 1200 m.a.g.l. over the entire continent in red, over the peripheral areas in green and over the high continental plateau in blue. Errorbars represent uncertainties as calculated by Lemonnier et al. (2019) and extrapolated to the entire continent.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00457v2-figure-fc65f599-d020-4944-818b-327e92e117ad.jpg)

profile is averaged over the 2007-2010 period, then area-averaged. The standard deviation on fig. 4.a is calculated from the time dimension (in months) and on fig. 4.b it is

calculated from the gridded spatial dimension. As observed in this figure, monthly temporal variability is low over the Antarctic continent. We note that the time standard deviation is continuous over the entire precipitation profile. The spatial variability is important, whether on the plateau or over peripheral areas. This spatial standard deviation decreases very rapidly with altitude. Thus, the precipitation rates recorded at high altitudes appear to be regionally homogeneous. Precipitation consistently increases from 7000 m in altitude until 1200 m, reaching a rate of approximately 160 mm.yr -1 . At this altitude, the precipitation gradient should become negative and the profile should reach a maximum in precipitation before inverting in the lowest layers over the peripheral areas (in green), according to several studies (Grazioli, Madeleine, et al., 2017; Dur´n-Alarc´n a o et al., 2019). The blue line presents the profile over the plateau. It is characterized by a very small amount of precipitation all along the profile and a relatively large dispersion : 30 mm.yr -1 with a σ -value of almost 50 % in the lowest bin. It suggests a high variability over the plateau. In comparison with the plateau, the green line presents the peripheral profile of precipitation, with higher precipitation values and a lower variability. Fig. 4 confirms that precipitation in Antarctica mostly occurs over the coasts. As highlighted by Grazioli, Madeleine, et al. (2017), there is an inversion of the precipitation profile in the lowest layers ( &lt; 1000 m) over the Dumont d'Urville station due to snow sublimation by katabatic winds. However, CloudSat does not see this low-level sublimation inversion layer because of ground clutter (Lemonnier et al., 2019; Palerme et al., 2019).

Figure 4. a) Averaged vertical profiles of precipitation over the 2007-2010 period of CloudSat observations in solid lines, filled areas are corresponding to the temporal σ standard deviations (1) over the entire continent in red. (2) over the plateau in blue. (3) over the peripheral areas in green. b) Same results as a) but the filled areas correspond to the spatial σ standard deviations (1) over the plateau in blue. (2) over the peripheral areas in green. Over the peripheral areas, the first bin value is 275 mm.yr -1 ± 11%, the precipitation rate at the first bin over the Antartic plateau is 34 mm.yr -1 ± 143% and the precipitation for the entire continent north of 82 ◦ is 163 mm.yr -1 ± 13%.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00457v2-figure-e1bdd3c2-4b1a-43a7-b4a7-dac1a04a81cd.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00457v2-figure-0ddde21e-0b37-4dd4-a886-733150864ef7.jpg)

## 3.2 From geographical areas

The re-gridded precipitation dataset is used to evaluate seasonal variabilities as well as vertical variabilities. This study also takes into account different regions of interest with the aim of providing useful diagnostics for climate models. Precipitation in each geographical region is longitudinally averaged at each CloudSat vertical bin for the entire period of observation. The same computation is applied to the atmospheric temperature from ECMWF operational weather analysis and digital elevation model. The results are presented in the first column of Fig. 5. The baroclinic structure of tempera-

ture, with the isotherms bending down towards the pole, is visible on each latitudinally averaged regions of fig. 5. Then we present in the second column the averaged snowfall rate at the fifth CloudSat level (at 1200 m.a.g.l., same studied altitude above ground level than Palerme et al. (2014)) in blue and its associated temperature in red over both the continent and the oceans (up to 60 ◦ ). The filled areas around the means in blue and red are the corresponding σ -standard deviations of the precipitation rate and temperature. The standard deviations are calculated over the spatial dimension for each region of interest.

Fig. 5.a represents precipitation structure over East Antarctica, the most homogeneous region of the Southern continent, and highlights an obstacle to the progression of precipitation flows with topography. The iso-precipitation rates are following isotherms over the ocean but are no longer parallel over the continent with a faster decrease in precipitation isolines. The precipitation rate drops from 600 mm.yr -1 to less than 10 mm.yr -1 along 10 degrees of latitude and 3000 meters of altitude difference. Fig. 5.b, which only considers retrievals at the first available bin over the surface, shows with more accuracy this latitudinal evolution. The precipitation rate over the plateau is very low and its associated standard deviation is relatively important (almost 100%). From 75 S and north-◦ ward, at the edge of the plateau, precipitation rate increases quickly before reaching a maximum value of ∼ 600 mm.yr -1 over the ocean at latitude 66 ◦ S. The temperature along the slope follows a similar evolution. Over the Southern Ocean, precipitation and temperature stabilize and reveal a more homogeneous oceanic climate.

Regarding the western part of the continent, there are some similarities with the eastern continent, such as the low precipitation rate on the plateau on fig. 5.c. The average precipitation rate at low altitude is lower than in the East, at 500 mm.yr -1 . The transition from continental to oceanic precipitation rates seems to be more gradual. The precipitation rate decreases from 500 mm.yr -1 to about 150 mm.yr -1 at an average elevation of 1500 m without stabilizing above the plateau, probably due to local orographic effects that characterize this region. The sea-ice cover is vast in this region of Antarctica, although the seasonal variability of its area is high, this precipitation pattern may be a persistent signal of winter precipitation behaviour. Indeed, surface humidity fluxes can be reduced in winter by the sea ice cover, preventing the atmosphere from collecting moisture by evaporation. On fig. 5.d, similarly to the eastern part of the continent, precipitation rate and temperature are increasing along the slope to the ocean. Gradient changes in the precipitation curve at 73 ◦ S and 78 S are caused by topography. ◦

Over the western ice-shelves on fig. 5.e, precipitation rate ranges from 100 to 500 mm.yr -1 at the first level. This could be due to a lack of moisture input to the atmosphere due to the sea ice cover and ice-shelves. Indeed, there is a slight correlation between sea ice concentration and low level clouds (Schweiger et al., 2008). When crossing the ice-shelves at 75 ◦ S, the precipitation increases quickly with latitude and altitude. Fig. 6 presents the seasonal variation in precipitation over this region. The reanalyses of ERA-Interim are used to present sea-ice cover fraction over the 2007-2010 period, which is shown in this figure using blue bars and ranging from 0 to 1 (0% to 100% of sea ice cover fraction). In winter (fig. 6.b), the almost complete sea ice coverage shows a maximum precipitation at 60 ◦ S which decreases rapidly above the ice-shelf following isotherms. In summer (fig. 6.a), the average sea ice coverage is less extensive but still over 80%. Also, there is higher averaged temperatures. Therefore, there is stronger precipitation rates closer to the coasts, with values ranging from 300 to 500 mm.yr -1 in the vicinity of the coasts at 75 ◦ S. Fig. 5.f shows a slow evolution of precipitation on ice-shelf with a very low standard deviation (less than 25% the value itself), while the temperature decreases linearly with latitude.

The peninsula is a particular region of interest, with a mountain ridge extending across the circumpolar atmospheric circulation. Over this region, precipitation is mostly driven by orogenic effects. Fig. 5.h shows a high precipitation rate with a large standard

9

Figure 5. First column presents the zonal averaged precipitation rate (shaded colors) and its corresponding atmospheric temperature (contours) of each region, as follows : a) for East continent (A in fig.2), c) for West continent (B.1 in fig.2), e) for West ice-shelves (B.2 in fig.2) and g) for Peninsula (C in fig.2). Black filled areas correspond to the zonal-average elevation of each area. Second column presents the fifth CloudSat level (1200 m.a.g.l.) of precipitation along the latitude (in blue) with its σ -standard deviation calculated over the spatial dimension. In red the associated temperature and its σ -standard deviation are represented. -10-

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00457v2-figure-77570d8c-c954-4ce8-a816-312e151aa798.jpg)

Figure 6. Zonal mean precipitation and its corresponding atmospheric temperature (contours) of the west Antarctic ice-shelves during the averaged summer and winter. The average period are a) December-January-February and b) June-July-August over the 2007-2010 period. The zonal averaged fractional sea-ice coverage obtained with ERA-Interim at a resolution of 0.75 ◦ is represented by the blue bars with a proportion ranging from 0 to 1.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00457v2-figure-50f92b45-4ce7-4f87-abf3-707c87dbe4e3.jpg)

deviation. The maximum snowfall rate is observed at the end of the peninsula with very high rates reaching 1000 mm.yr -1 at 65 ◦ S as seen on fig. 5.g. At 65 S the precipitation ◦ reaches its maximum and there is an abrupt decrease in temperature. At this latitude, the topographic gradient of the peninsula is very intense. Over the peninsula, from 64 ◦ S to 74 ◦ S, precipitation rate is high and strongly affected by topography. It ranges from 500 to 700 mm.yr -1 before decreasing from 75 ◦ S to 82 ◦ S with precipitation rates lower than 200 mm.yr -1 . The topographic obstacle at 75 ◦ S with an average elevation of nearly 1000 m, but locally reaching 2500 m, results in two precipitation regimes, intense over the peninsula and low on the continental side. To better understand these results, we have performed a high-resolution climatology (see fig. 7) of 0.1 ◦ in longitude by 0.1 ◦ in latitude on the first available level of CloudSat by averaging measured precipitations rates and corresponding temperatures in each cell as we did with a 2 ◦ in longitude by 1 ◦ in latitude grid for the whole continent (see section 2). This figure clearly shows the effect of topography on precipitation and temperature fields. Both temperature and precipitation are higher west of the peninsula, with very high precipitation rates reaching 1000 mm.yr -1 observed along the west coast of the peninsula. On the east side, the temperature is lower due to the Larsen C ice-shelf. As seen in fig. 7.b, the temperature above the peninsula is much lower than above the ocean with values ranging from -15 to -25 ◦ C. This highresolution view provides a better understanding of the precipitation and temperature variations observed in fig. 5.g and 5.h. Indeed, they show an average structure of warmer and wetter regions in the west of the peninsula with colder and drier regions in the east of the peninsula, thus increasing the variability of temperature and precipitation.

Fig. 8.a summarizes the precipitation structure over the entire continent. It clearly shows that precipitation over ocean follows isotherms and thus evolves with temperature. Then precipitation isolines depart from isotherms over the topographic slope : precipitation rates are drastically decreasing when confronting the plateau. Indeed, the air masses rising along the topographical slope are now far from oceanic moisture sources, thus explaining the divergence of isotherms and precipitation contours. The average evolution of the topographical slope of Antarctica rises from sea level to 2000 metres above sea level from 66 ◦ S to 78 ◦ S. This topographical evolution is associated with an evolution of precipitation from 700 mm.yr -1 to 150 mm.yr -1 . On the plateau, precipitation can reach values as low as 10 mm.yr -1 . Fig. 8.b shows the evolution of the zonal mean

Figure 7. Averaged annual precipitation and temperature over the 2007-2010 period at 1200 m above ground level. The grid resolution is 0.1 ◦ in longitude by 0.1 ◦ in latitude. Precipitation is obtained from the 2C-SNOW-PROFILE product presented in section 2 and temperature is obtained from the ECMWF-AUX product. The grey grid represents the resolution of the climatology presented in this study. CloudSat passes are not always located in exactly the same place. Thus the passes plotted on this figure with the red dots at the edge of the trackes may contain only a few strong precipitation events with a very low number of satellite overpasses ( &lt; 10 passes). At this resolution and according to Souverijns et al. (2018), there are more omissions than commissions of precipitation events by the CPR and few overflights. The average precipitation rate will thus could be very high.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00457v2-figure-310ee588-98b7-4e12-bdb0-8a84067a1a0c.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00457v2-figure-9fa361ed-5377-448f-9362-c3c52742b18f.jpg)

precipitation rate at the fifth CloudSat bin. Over the plateau, precipitation rate is low but its standard deviation is relatively high, indicating that precipitation is not homogeneous, but brought by local events that transport moisture from the Southern Ocean. Along the topographic slope between coasts and plateau, precipitation and temperature are increasing, then precipitation reaches a maximum of 600 mm.yr -1 at 66 ◦ S.

Figure 8. a) Latitudinally average of the precipitation structure and its corresponding atmospheric temperature (contours) for the entire continent. b) First CloudSat level of precipitation along the latitude (in blue) with its σ -spatial variation, and associated temperature (in red) with its σ -spatial standard deviation.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00457v2-figure-5f526d23-ca2a-48d2-924f-6cb9630c9063.jpg)

## 4 Precipitation distribution over the Antarctic continent

In order to have a better representation of the precipitation over the ice-cap and to develop a new diagnostic tool for models, we looked at the distribution of the 2C-SNOWPROFILE and ECMWF-AUX products raw data over the entire 3D dataset south of 60 S. ◦ To do so, we have arranged these points in temperature bins of linear size and precipitation rate bins of logarithmic size. The study of the general structure of precipitation performed in section 3 showed that the main climatological differences are between the oceanic border regions and the high plateau. In this section we therefore keep distinguishing the peripheral precipitation from the precipitation over the plateau.

## 4.1 Evolution of precipitation distributions with altitude

Each vertical CloudSat level above the surface has been studied separately. Fig. 9 presents histograms of the distribution of precipitation rate observations evolving with altitude. In both areas, the precipitation distribution of the first level is truncated at 0.01 mm/hr, while for the other levels, there is no truncation. However, we observe precipitation rates below 0.01 mm/hr in the 2C-SNOW-PROFILE product. The Snow Retrieval Status values for these measurements are greater than 3, indicating that an error in the retrieval process could occur. It would therefore seem that ground clutter, in addition to generating excessive precipitation rates, does not allow sufficient confidence to be placed in small rates of near-surface precipitation. It is worth noting that we calculated the average precipitation rate at the second available level of CloudSat (at 1440 m.a.g.l.), by taking into account and by excluding the points recorded under 0.01 mm/hr to make sure that the ground clutter does not affect the averaged precipitation rate. This has no impact on the averaged snowfall rate at 1440 m.a.g.l. Doing so, we indeed verified that the missing observations do not affect the average precipitation rate at 1200 m.a.g.l., neither for peripheral regions nor on the plateau.

Figure 9. Histogram plots of snowfall rates in mm.hr -1 and vertical evolutions ( a) for the peripheral areas ( &lt; 2250 m) and b) for the plateau ( &gt; 2250 m) at different altitudes above ground level. Black dashed lines are passing through maxima of each vertical bin distribution.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00457v2-figure-15e1ce68-1d52-43b6-866d-601b059d10f4.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00457v2-figure-dc5c24cd-68d0-494f-8b3a-a06ebef0db00.jpg)

We observe that the evolution of the number of observations and measured precipitation rates as a function of altitude above the coasts and plateau is similar, except for the first level under investigation. For coasts, distribution maxima at 1680 m.a.g.l. (7 th bin) and 2160 m.a.g.l. (9 th bin) do not vary in precipitation rate (0.0221 mm.hr -1 ) but decrease from 14% in records numbers. It then reaches 0.0148 mm.hr -1 at 3120 m.a.g.l. (13 th bin), 0.0099 mm.hr -1 at 4080 m.a.g.l. (17 th bin) and 0.0055 mm.hr -1 at 5040 m.a.g.l. (21 th bin). Along these levels, there appears to be a consistency in the logarithmic de-

cay of the peaks. Above 6000 m.a.g.l., the maximum precipitation rate stabilizes around a value of 0.0030 mm.hr -1 . Over the plateau, the position of the maximum in the lower vertical bins is at a lower precipitation rate than on the coasts. Apart from the truncated level, there is a decrease from 0.0099 mm.hr -1 to 0.0081 mm.hr -1 between 1680 m.a.g.l. (7 th bin) and 2160 m.a.g.l. (9 th bin). A thousand meters above at 4080 m.a.g.l. (17 th bin) the peak is located at 0.0030 mm.hr -1 . Finally, at higher altitudes, the maximum remains in the same precipitation rate, centered on 0.0025 mm.hr -1 . For each region, the positions of the maxima at all altitudes are listed in Appendix A.

## 4.2 Precipitation at saturation by forced lifting

CloudSat observations are also presented on scatterplots in fig. 10 at different altitudes above ground level. The colorbar indicates the relative number of observations for a given precipitation rate and temperature. Over the peripheral regions, there is a significant spread both in precipitation rate and temperature. Over the plateau, the spread in precipitation is smaller and the distribution is centered at lower temperatures. In both regions, there seems to be a relationship between precipitation and temperature, and the spread in precipitation rate decreases when altitude increases.

To further understand these distributions, we can derive an analytical relationship between precipitation and temperature by assuming that precipitation over the Antarctic ice cap is mostly controlled by the topographic lifting of oceanic air parcels that are initially close to saturation. We first assume that all the horizontal motion is converted into vertical motion. Then, assuming a mean horizontal wind velocity over a given slope, we can deduce the vertical speed of the oceanic air parcels. Then, this vertical wind speed can be used to lift air parcels along the moist adiabatic lapse rate. At each level, all the water vapor entering the parcel will be in excess of saturation and said to be equal to the precipitation mass flux. By integrating over the vertical direction, the resulting precipitation rate can be written as:

$$P _ { r } = - \frac { w } { \rho _ { w a t e r } } \int ^ { z } \rho _ { a t m } \frac { L q _ { s a t } ( T, p ) } { R _ { v a p } T ^ { 2 } } \Gamma _ { s a t } d z$$

where w is the vertical wind speed, ρ water is the water density and ρ atm is the air density, L is the latent heat of sublimation, q sat is the humidity at saturation, R vap is the specific gas constant for wet air and Γ sat is the moist adiabatic lapse rate. The demonstration is in the Appendix B. This integral equation is resolved by using 240m ∆ z corresponding to the CloudSat bin vertical dimensions. The precipitation P r expressed by equation 1 is given in water equivalent mm.hr -1 . It is represented by the black curves with markers on fig. 10. The values used for w are 0.0001 m.s -1 , 0.001 m.s -1 , 0.01 m.s -1 , 0.1 m.s -1 and 1 m.s -1 . They are respectively represented in fig. 10 by black dashed lines with diamond, circle, triangle, square and star markers. The rows are corresponding to CloudSat vertical bins above the surface, at respectively 1200, 2160, 3120 and 5040 m.a.g.l. which are corresponding to the 5 th , 9 th , 13 th and 21 th vertical bins above the surface. The white solid line σ limit defining the standard deviation of the population distribution provides additional information on precipitation. In both cases, the spread in precipitation is greater when T increases and follows the trend with the triangles markers with a vertical velocity of 0.01 m/s. This vertical velocity w can be assessed by a simple relationship between the slope dz/dx and horizontal wind velocity. For example, based on fig. 8, a horizontal wind u = 5 m.s -1 (this averaged wind velocity is in agreement with general circulation models, such as the IPSL-CM model in our case) advected over the average slope of the Antarctic ice cap, which extends over 10 degrees of latitude and 2000 m of vertical elevation gain, implies a vertical motion of u.dz/dx = 0 01 m.s . -1 . A much lower vertical wind value implies an air advection along a gentler slope. A vertical speed of 0.1 m.s -1 and above can be explained by sharp topographical obstacles,

such as mountain ranges in West Antarctica and along the peninsula, or stronger largescale winds. The higher the vertical velocity, the higher the precipitation rate.

At 1200 m.a.g.l. over peripheral areas on fig. 10.a, most of the precipitation records are located between -32 ◦ C and -14 C, and between 0.006 mm.hr ◦ -1 and 0.25 mm.hr -1 . The distribution is bounded between circle and square markers. When continuing to increase in altitude, the distribution follows the orientation of the analytical relationship. At the highest level considered in the fig. 10.g, the spread is located between -55 ◦ C and -31 ◦ C in temperature and ranges from 0.001 mm.hr -1 to 0.005 mm.hr -1 in precipitation rate. Throughout the ascent along the CloudSat bins, the distribution evolves just over the line with the triangles markers, which means a vertical speed of about 0.02 m.s -1 . This is again consistent with an orographic precipitation on the margins of the ice-sheet. The table A2 shows the precipitation and temperature locations of these distributions for levels from 5 th to 25 th CloudSat vertical bins. Above the plateau, on fig. 10.b at 1200 m.a.g.l. the distribution is located between -43 ◦ C and -27 C. The lower boundary of the ◦ distribution is slightly above the theoretical triangles markers line. At this altitude, distribution is highly impacted by ground clutter. At a higher altitude of 2160 m.a.g.l. fig. 10.d, this spread ranges from 0.002 mm.hr -1 to 0.05 mm.hr -1 with a temperature interval of -49 ◦ C to -30 ◦ C. As it evolves with altitude, the distribution continues to follow the line with the triangles markers, but only slightly above it. The table A3 shows the precipitation and temperature locations of these distributions for levels from 5 th to 25 th CloudSat vertical bin.

Fig. 11 summarizes distribution of precipitation and temperature summed for all vertical levels. It shows that the distribution over peripheral regions reaches higher temperatures, so its precipitation dispersion is larger at these temperatures. Over peripheral regions, the density-plot is defined for large-scale vertical velocities w ranging from 0.0001 (circle markers) to 0.1 m.s -1 (square markers), while over the continental plateau, the distribution is bounded by triangle markers and square markers. Some observations reach very high precipitation rates, and are higher than the analytical relationship showed by square markers with a vertical velocity of 0.1 m.s -1 . This is either due to high largescale winds analogous to extreme events or high slope analogous to local topographical obstacles along the peninsula and mountain ranges. We calculated a precipitation rate with a wind speed of 0.2 m.s -1 . This agrees either with a slope of 0.2%, corresponding to the slope of the East Antarctic ice cap and a strong large-scale wind blowing at 100 m.s -1 or with a mean large-scale wind of 5 m.s -1 and a steep slope of 4%. This covers many possible combinations of strong large scale winds and slopes ranging from the slope of the East Antarctic ice cap to much steeper slopes. Observations on the plateau and coasts that exceed the resulting analytical relationship with a wind speed of 0.2 m.s -1 are showed on fig. 12. In order to verify whether these measurements are made over areas of high topographic gradient, we have reported these measurements on a high-resolution topographic map Greene et al. (2017) and Howat et al. (2019) in figure 13.

Fig. 13 shows a high-resolution (200 m) topographic map processed by Greene et al. (2017); Howat et al. (2019) over 4 areas where these very high precipitation rates have been recorded by CloudSat. These measurements suggest a vertical velocity of advection at a large scale greater than 0.1 m.s -1 , possible only in regions with a high topographic gradient. The slope of the region of Terre Ad´lie showed on fig. e 13.a respects the topographic criterion causing strong vertical winds. Fig. 13.b presenting Ellsworth Land and Mount Vinson massif satisfies the criterion of a high topographic gradient. Indeed, the Mount Vinson massif is a mountain range and the highest point in Antarctica. The peninsula on fig. 13.c, as previously discussed, is a topographic barrier that crosses and obstructs areas with major circumpolar air currents. And the fourth region on fig. 13.d, Dronning Maud, is also respecting the topographic criterion causing strong vertical winds with the presence of a mountainous barrier a few hundred kilometres from the coast. On the four cases, the red dots corresponding to high precipitation measurements

Figure 10. First column presents scatter plots of precipitation in mm/hr and temperature in ◦ C at different altitudes over the peripheral areas. Second column presents the same results over plateau area. The dashed black lines are the assumptions of theoretical precipitation rates calculated using equation 1 for vertical wind w values of 0.0001 m.s -1 (diamond markers), 0.001 m.s -1 (circle markers), 0.01 m.s -1 (triangle markers), 0.1 m.s -1 (square markers) and 1 m.s -1 (star markers). Colorbars are relative amounts of observations per CloudSat bin. White contours represent the σ standard deviation of the distributions. -16-

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00457v2-figure-5681f39d-5197-48dc-94e9-a32a50a2d1d5.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00457v2-figure-c7fd3017-f934-4087-b7b4-e94a0a432f2c.jpg)

Figure 11. a) Scatter plots of precipitation in mm.hr -1 and temperature in ◦ C at all altitudes over the peripheral areas. a) Same result over plateau area. The dashed black lines are the assumptions of theoretical precipitation rates calculated using equation 1 for vertical wind w values of 0.0001 m.s -1 (diamond markers), 0.001 m.s -1 (circle markers), 0.01 m.s -1 (triangle markers), 0.1 m.s -1 (square markers) and 1 m.s -1 (star markers). Colorbars are relative amounts of observations per CloudSat bin. White contours represent the σ standard deviation of the distributions. Considered levels are 1200, 2160, 3120 and 5040 m.a.g.l.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00457v2-figure-93a7bdb8-484f-410c-8f2f-4df97be0589d.jpg)

Figure 12. Observations on the plateau and coasts that exceed the resulting analytical relationship with a wind speed of 0.2 m.s -1 (red scattered points). The dashed black line represents the analytical relationship with a vertical speed of 0.2 m.s -1 . There are 1251 records shown here.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00457v2-figure-c59a0d13-fe6a-4fb3-b77a-bb139ac7c59c.jpg)

follow each others, they correspond to the same satellite track at a given time. These are rare occasional events of massive precipitation that contribute significantly to the accumulation of snow on the continent.

Fig. 14 summarizes the behaviour of precipitation over the Antarctic continent. Precipitation over the peripheral areas occurs with higher temperature by forced lifting of air masses along the topographic slope. The different topographical slopes as well as the strength of the large-scale horizontal wind thus generates different vertical wind velocities. According to equation 1, the variability of these winds consequently generates a wide spread in precipitation rates, identified by the white dashed line contours. Then over the plateau, the available water quantities, temperatures, low slopes and low largescale horizontal winds cause precipitation to evolve following Clausius-Clapeyron rela-

Figure 13. High resolution (200 m) topographic map (Greene et al., 2017; Howat et al., 2019) over Terre Ad´lie, Ellsworth Land, the peninsula and Dronning Maud Land. The red points core respond to the measurements greater than the 0.2 m.s -1 vertical velocity precipitation hypothesis presented in fig. 12.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00457v2-figure-83e177b9-7db4-4e1d-b421-d31e74ed6185.jpg)

tion with a small spread in precipitation rate. This is represented on fig. 14 by the white solid contours.

## 5 Conclusion

Precipitation is mostly considered as a surface variable, climatologies typically only reporting the 2D horizontal distribution at the surface. The CloudSat radar dataset now allows to explore the 3D structure of precipitation, and this study provides insight into the origins of precipitation and its evolution along the vertical dimension over Antarctica. The 2C-SNOW-PROFILE product has been explored on the three spatial dimensions and the temporal dimension from 2007 to 2010 over the Antarctic continent. This 3D-dataset is computed following Palerme et al. (2014, 2019) with a horizontal resolution of 1 ◦ of latitude by 2 ◦ of longitude grid and a 240 m vertical resolution in order to optimally represent the Southern polar climate. It has been studied over several different regions, geographically divided into four regions (East continent, West continent, West ice-shelves and Peninsula), and topographically divided into two regions: the peripheral areas (where topography &lt; 2250 m) and the Antarctic plateau (where topography &gt; 2250 m). The distinction between peripheral areas and continental plateau revealed that most of the precipitation is located over the peripheral areas with low relative seasonal variability, in contrast with the plateau where it is relatively high. The study of the four geographical regions revealed many differences between the Peninsula and the western iceshelves as compared to the continental eastern and western parts. Indeed, the Peninsula

Figure 14. Scatter plots of precipitation in mm.hr -1 and temperature in ◦ C at all altitudes over the Antarctic continent. The colorbar indicates the relative number of observations for a given precipitation rate and temperature. The dashed black lines are the assumptions of theoretical precipitation rates calculated with the moisture convergence equation for vertical wind w values of 0.0001 m.s -1 (diamond markers), 0.001 m.s -1 (circle markers), 0.01 m.s -1 (triangle markers) 0.1 m.s -1 (square markers) and 1 m.s -1 (star markers). Dotted white line represents the 1σ distribution of precipitation over the peripheral areas and solid white line represents the 1σ distribution of precipitation over the plateau.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00457v2-figure-297f74b3-8685-44ab-b673-0fc3498493fc.jpg)

is an area where precipitation is strong and mainly driven by local orography and snowfall rates are low over the ice-shelves, while precipitation is increasing along the slope to the ocean over continental eastern and western regions. There also seems to be a strong dichotomy in the behaviour of precipitation over the ice-shelves. This could be related to the sea ice cover, which is able to stratify the atmosphere by cooling it, thus limiting the mixing and thus precipitation.

The pre-gridded CloudSat product is then statistically studied in order to evaluate the distribution of snowfall rates and associated temperatures over the peripheral areas and the plateau at each CloudSat vertical level. This shows that precipitation, whatever the altitude, is mainly driven by large-scale convergence of moist fluxes over the topography. This has been confirmed by comparing the observed precipitation/temperature distributions with an analytical relationship of precipitation as a function of temperature and a vertical advection velocity that is directly dependent on large scale horizontal wind strength and slope. Over the peripheral regions, the dispersion of precipitation rates as a function of temperature is large. This may mean that both topographical slopes and vertical winds are highly diversified. Over the continental plateau, the dispersion is reduced and results from gentler slopes than in peripheral regions, as well as slower

vertical uplifts. Thus, the precipitation dispersion at a given temperature of these distributions is justified by varying degrees of horizontal advection on variable slopes. The averaged observed precipitation distribution follows the analytical curve corresponding to a vertical wind of 0.01 m.s -1 , which is typical of the values found in general circulation models, for example in our case the IPSL-CM. This new study of the CloudSat precipitation product provides new and innovative tools to evaluate climate models with a three-dimensional view of the atmospheric structure of precipitation. On the one hand, the use of average continental precipitation profiles and zonal averages are diagnostic tools that are easily comparable to models. On the other hand, precipitation rate and temperature scatter-plots require high-frequency model outputs, but lead to a powerful way of evaluating transport, thermodynamics and microphysical processes in further detail.

## Appendix A Location of the maximum precipitation rate as a function of altitude

Selected bin values for precipitation rate and temperature as presented in figures 9, 11, 12 and 10.

## Temperature:

[-60.0000 -59.1429 -58.2857 -57.4286 -56.5714 -55.7143 -54.8571 -54.0000 -53.1429 -52.2857 -51.4286 -50.5714 -49.7143 -48.8571 -48.0000 -47.1429 -46.2857 -45.4286 -44.5714 -43.7143 -42.8571 -42.0000 -41.1429 -40.2857 -39.4286 -38.5714 -37.7143 -36.8571 -36.0000 -35.1429 -34.2857 -33.4286 -32.5714 -31.7143 -30.8571 -30.0000 -29.1429 -28.2857 -27.4286 -26.5714 -25.7143 -24.8571 -24.0000 -23.1426 -22.2857 -21.4286 -20.5714 -19.7143 -18.8571 -18.0000 -17.1429 -16.2857 -15.4286 -14.5714 -13.7143 -12.8571 -12.0000 -11.1429 -10.2857 -9.4286 -8.5714 -7.7143 -6.8571 -6.0000 -5.1429 -4.2857 -3.4286 -2.5714 -1.7143 -0.8571 0.0000]

Precipitation rate: [0.0001 0.0001 0.0001 0.0002 0.0002 0.0003 0.0003 0.0004 0.0005 0.0006 0.0007 0.0009 0.0011 0.0013 0.0016 0.0020 0.0025 0.0030 0.0037 0.0045 0.0055 0.0067 0.0081 0.0099 0.0122 0.0148 0.0181 0.0221 0.0270 0.0330 0.0403 0.0493 0.0602 0.0735 0.0898 0.1097 0.1340 0.1636 0.1998 0.2441 0.2981 0.3641 0.4447 0.5432 0.6634 0.8103 0.9897 1.2088 1.4765 1.8034 2.2026 2.6903 3.2850 4.0139 4.9021 5.9874 7.3130 8.9321]

The number of decimals for these bins after the decimal point was chosen according to the sensitivity of the CPR measurements.

## Appendix B Precipitation at saturation by forced lifting demonstation

All the air advected into the pole would condense by forced lifting due to the topographical landmass of the continent. We assume a purely vertical motion initiated from a horizontal advection at the origin of large-scale precipitation. The vertical motion of the air mass implies a cooling and a variation in the saturation vapour pressure, and therefore in the precipitation rate. We use moisture flow convergence equations and test several vertical wind values in order to explain precipitation observations. This is based on the net moisture balance in a steady state saturated atmosphere and can be derived from:

$$\frac { D q } { D t } = S$$

where D/Dt is the lagrangian derivative:

$$\frac { D } { D t } = \frac { \partial } { \partial t } + u \frac { \partial } { \partial x } + v \frac { \partial } { \partial y } + w \frac { \partial } { \partial z }$$

Table A1. Location of the maximum precipitation rate as a function of altitude, as presented in figure 9. Values for all CloudSat bins ranging from 1200 m.a.g.l. (5 th bin) to 6000 m.a.g.l. (25 th bin) above the surface are presented. The format is as follows: 'precipitation rate in mm.hr -1 - number of observations'.

| Altitude above ground level   | Coastal areas   | Continental plateau   |
|-------------------------------|-----------------|-----------------------|
| 1200 m                        | 0.0122 - 377278 | 0.0099 - 321044       |
| 1440 m                        | 0.0181 - 377234 | 0.0122 - 260956       |
| 1680 m                        | 0.0221 - 351080 | 0.0099 - 212493       |
| 1920 m                        | 0.0221 - 325865 | 0.0099 - 184575       |
| 2160 m                        | 0.0221 - 305071 | 0.0081 - 163445       |
| 2400 m                        | 0.0181 - 288393 | 0.0067 - 146641       |
| 2640 m                        | 0.0181 - 272086 | 0.0055 - 132604       |
| 2880 m                        | 0.0181 - 257160 | 0.0055 - 121014       |
| 3120 m                        | 0.0148 - 243191 | 0.0045 - 110872       |
| 3360 m                        | 0.0121 - 229763 | 0.0037 - 100530       |
| 3600 m                        | 0.0121 - 216356 | 0.0030 - 91289        |
| 3840 m                        | 0.0099 - 201847 | 0.0030 - 83873        |
| 4080 m                        | 0.0099 - 187903 | 0.0030 - 75265        |
| 4320 m                        | 0.0081 - 175297 | 0.0025 - 66611        |
| 4560 m                        | 0.0081 - 164391 | 0.0025 - 57850        |
| 4800 m                        | 0.0067 - 153230 | 0.0025 - 49477        |
| 5040 m                        | 0.0055 - 142044 | 0.0025 - 40760        |
| 5280 m                        | 0.0045 - 131296 | 0.0025 - 32798        |
| 5520 m                        | 0.0037 - 121678 | 0.0025 - 25271        |
| 5760 m                        | 0.0030 - 111769 | 0.0025 - 19258        |
| 6000 m                        | 0.0030 - 103270 | 0.0025 - 14080        |

Table A2. Location of the σ deviation of the maximum precipitation rate as a function of altitude for the peripheral areas. Values for all CloudSat bins ranging from 1200 m.a.g.l. (5 th bin) to 6000 m.a.g.l. (25 th bin) above the surface are presented. The format is as follows: 'maximum observation concentration'[extrema extrema]. The bold values indicate that the maximum distribution is overlaid with the 1σ distribution boundary.

| Altitude above ground level   |   Precipitation rate in mm.hr - 1 | Precipitation rate in mm.hr - 1   | Precipitation rate in mm.hr - 1   | Temperature in ◦ C              |
|-------------------------------|-----------------------------------|-----------------------------------|-----------------------------------|---------------------------------|
| 1200 m                        |                            0.0148 | [0.0055                           | - 0.2981]                         | -25.7143 [-32.5714 - -14.5714]  |
| 1440 m                        |                            0.0181 | [0.0055                           | - 0.2981]                         | -25.7143 [-33.4286 - -14.5714]  |
| 1680 m                        |                            0.027  | [0.0037                           | - 0.3641]                         | -25.7143 [-35.1429 - -14.5714]  |
| 1920 m                        |                            0.027  | [0.0030                           | - 0.2981]                         | -26.5714 [-36.0000 - -16.2857]  |
| 2160 m                        |                            0.0221 | [0.0020                           | - 0.2981]                         | -28.2857 [-17.1429 - -36.8571]  |
| 2400 m                        |                            0.0221 | [0.0020                           | - 0.2441]                         | -30.0000 [-38.5714 - -18.0000]  |
| 2640 m                        |                            0.027  | [0.0013                           | - 0.1998]                         | -28.2857 [-40.2857 - -18.8571]  |
| 2880 m                        |                            0.027  | [0.0013                           | - 0.1636]                         | -28.2857 [-41.1429 - -20.5714]  |
| 3120 m                        |                            0.0221 | [0.0013                           | - 0.1636]                         | -29.1429 [-42.0000 - -22.2857]  |
| 3360 m                        |                            0.0221 | [0.0013                           | - 0.1339]                         | -30.0000 [-43.7143 - -23.1429]  |
| 3600 m                        |                            0.0181 | [0.0013                           | - 0.1339]                         | -31.7143 [-45.4286 - -24.0000]  |
| 3840 m                        |                            0.0122 | [0.0016                           | - 0.1097]                         | -35.1429 [-48.0000 - -24.8571]  |
| 4080 m                        |                            0.0148 | [0.0016                           | - 0.0898]                         | -36.0000 [-48.8571 - -26.5714]  |
| 4320 m                        |                            0.0099 | [0.0013                           | - 0.0735]                         | -38.5714 [-49.7143 - -28.2857]  |
| 4560 m                        |                            0.0099 | [0.0013                           | - 0.0602]                         | -37.7143 [-50.5714 - -29.1429]  |
| 4800 m                        |                            0.0067 | [0.0013                           | - 0.0602]                         | -41.1429 [ -41.1429 - -30.0000] |
| 5040 m                        |                            0.0099 | [0.0013                           | - 0.0493]                         | -40.2857 [-54.8571 - -30.8571]  |
| 5280 m                        |                            0.0055 | [0.0013                           | - 0.0404]                         | -43.7143 [-55.7143 - -33.4286]  |
| 5520 m                        |                            0.0055 | [0.0013                           | - 0.0330]                         | -45.4286 [-56.5714 - -35.1429]  |
| 5760 m                        |                            0.0037 | [0.0013                           | - 0.0270]                         | -48.0000 [-56.5714 - -35.1429]  |
| 6000 m                        |                            0.003  | [0.0013                           | - 0.0270]                         | -48.0000 [-57.4286 - -36.0000]  |

Table A3. Location of the σ deviation of the maximum precipitation rate as a function of altitude for the plateau. Values for all CloudSat bins ranging from 1200 m.a.g.l. (5 th bin) to 6000 m.a.g.l. (25 th bin) above the surface are presented. The format is as follows: 'maximum observation concentration'[extrema extrema].

| Altitude above ground level   |   Precipitation rate in mm.hr | - 1               |   Precipitation rate in mm.hr | Temperature in ◦ C    |
|-------------------------------|-------------------------------|-------------------|-------------------------------|-----------------------|
| 1200 m                        |                        0.0099 | [0.0081 - 0.0493] |                      -36.8571 | [-43.7143 - -27.4286] |
| 1440 m                        |                        0.0122 | [0.0055 - 0.0602] |                      -37.7143 | [-45.4286 - -27.4286] |
| 1680 m                        |                        0.0099 | [0.0037 - 0.0602] |                      -38.5714 | [-46.2857 - -28.2857] |
| 1920 m                        |                        0.0099 | [0.0030 - 0.0493] |                      -38.5714 | [-47.1429 - -28.2857] |
| 2160 m                        |                        0.0081 | [0.0020 - 0.0493] |                      -40.2857 | [-48.8571 - -30.0000] |
| 2400 m                        |                        0.0099 | [0.0016 - 0.0403] |                      -40.2857 | [-49.7143 - -31.7143] |
| 2640 m                        |                        0.0067 | [0.0013 - 0.0330] |                      -42      | [-50.5714 - -32.5714] |
| 2880 m                        |                        0.0055 | [0.0013 - 0.0330] |                      -43.7143 | [-52.2857 - -34.2857] |
| 3120 m                        |                        0.0055 | [0.0013 - 0.0270] |                      -45.4286 | [-54.0000 - -35.1429] |
| 3360 m                        |                        0.0037 | [0.0013 - 0.0221] |                      -48      | [-54.8571 - -36.8571] |
| 3600 m                        |                        0.0037 | [0.0013 - 0.0221] |                      -48.8571 | [-55.7143 - -38.5714] |
| 3840 m                        |                        0.003  | [0.0016 - 0.0181] |                      -49.7143 | [-56.5714 - -39.4286] |
| 4080 m                        |                        0.003  | [0.0011 - 0.0181] |                      -51.4286 | [-58.2857 - -41.1429] |
| 4320 m                        |                        0.0025 | [0.0011 - 0.0148] |                      -53.1429 | [-59.1429 - -42.8571] |
| 4560 m                        |                        0.003  | [0.0011 - 0.0122] |                      -53.1429 | [-50.5714 - -29.1429] |
| 4800 m                        |                        0.0025 | [0.0011 - 0.0122] |                      -54.8571 | [-60.0000 - -45.4286] |
| 5040 m                        |                        0.0025 | [0.0013 - 0.0122] |                      -56.5714 | [-60.0000 - -46.2857] |
| 5280 m                        |                        0.0025 | [0.0013 - 0.0099] |                      -56.5714 | [-60.0000 - -46.2857] |
| 5520 m                        |                        0.0025 | [0.0011 - 0.0099] |                      -59.1429 | [-60.0000 - -48.0000] |
| 5760 m                        |                        0.0025 | [0.0011 - 0.0099] |                      -59.1429 | [-60.0000 - -49.7143] |
| 6000 m                        |                        0.0025 | [0.0011 - 0.0081] |                      -55.7143 | [-60.0000 - -49.7143] |

with u , v and w , representing the standard three-dimensional wind components, and q is the specific humidity. S represents the budget of water vapor, which is defined by the difference between sources and sinks of water vapor following air parcel motion. S typically takes the form Ev -Pr , where Ev is the evaporation rate into the air parcel and Pr is the precipitation rate.

$$\frac { \partial q } { \partial t } + u \frac { \partial q } { \partial x } + v \frac { \partial q } { \partial y } + w \frac { \partial q } { \partial z } = E v - P r$$

If we assume that the lifted parcels are brought to saturation, Ev = 0 kg.kg -1 .s -1 , q = q sat ( T, p ) and in steady state, ∂q/∂t = 0 kg.kg -1 .s -1 so:

$$\vec { u }. \vec { \nabla } q _ { s a t } ( T, p ) = - P r$$

Assuming a purely vertical motion, we finally reach the form:

$$w \frac { \partial q _ { s a t } ( T, p ) } { \partial z } = - P r$$

We can therefore assess the precipitation rate by integrating a vertical flux over the column observed by CloudSat at a given constant w :

$$P _ { r } = - w \int ^ { z } \rho _ { a t m } \frac { \partial q _ { s a t } ( T, p ) } { \partial z } d z$$

The term ∂q sat ( T, p /∂z ) can be decomposed and then described from Clausius Clapeyron's equation as follow:

$$\frac { d e _ { s a t } } { e _ { s a t } } = \frac { L } { R _ { v a p } T ^ { 2 } } d T$$

with q sat ( T, p ) = 0 622 . e sat /p . At a constant pressure:

$$\frac { \partial q _ { s a t } } { \partial z } = \frac { \partial q _ { s a t } ( T, p ) } { \partial T }. \frac { \partial T } { \partial z } = \frac { L q _ { s a t } ( T, p ) } { R _ { v a p } T ^ { 2 } } \frac { \partial T } { \partial z }$$

with the latent heat of sublimation L and the specific gas constant for wet air R vap . Since we assume that we are at saturation, ∂T / ∂z is the moist adiabatic lapse rate Γ sat with a value of -6.5 K.km -1 . The definitive form of the precipitation equation is:

$$P _ { r } = - w \int ^ { z } \rho _ { a t m } \frac { L q _ { s a t } ( T, p ) } { R _ { v a p } T ^ { 2 } } \Gamma _ { s a t } d z$$

## Acknowledgments

This work was supported by the French National Research Agency (Grant number : ANR15-CE01-0003). CloudSat data is freely available via the CloudSat Data Processing Center (http://www.cloudsat.cira.colostate.edu/). The authors thank Karine Marquois and the IT department of the Laboratoire de M´t´orologie Dynamique for the informatics e e support to generate the 3D climatology. The authors thank Anna-Lea Albright, Margaux Valls, Luca Montabone and M´ elanie Thiriet for the proofreading.

## References

Bromwich, D. H. (1988). Snowfall in high southern latitudes. Reviews of Geophysics , 26 (1), 149-168.

Bromwich, D. H., &amp; Parish, T. R. (1998). Meteorology of the antarctic. In Meteorology of the southern hemisphere (pp. 175-200). Springer.

Church, J. A., Clark, P. U., Cazenave, A., Gregory, J. M., Jevrejeva, S., Levermann,

A., . . . others (2013). Sea-level rise by 2100. Science , 342 (6165), 1445-1445. Dur´ an-Alarc´n, C., Boudevillain, B., Genthon, C., Grazioli, J., Souverijns, N., van o Lipzig, N. P. M., . . . Berne, A. (2019). The vertical structure of precipitation at two stations in east antarctica derived from micro rain radars. The Cryosphere , 13 (1), 247-264. Retrieved from https://www.the-cryosphere .net/13/247/2019/ doi: 10.5194/tc-13-247-2019

Dye, J. E., Knight, C. A., Toutenhoofd, V., &amp; Cannon, T. W. (1974). The mechanism of precipitation formation in northeastern colorado cumulus iii. coordinated microphysical and radar observations and summary. Journal of the Atmospheric Sciences , 31 (8), 2152-2159.

Eisen, O., Frezzotti, M., Genthon, C., Isaksson, E., Magand, O., van den Broeke, M. R., . . . others (2008). Ground-based measurements of spatial and temporal variability of snow accumulation in east antarctica. Reviews of Geophysics , 46 (2).

Findeisen, W., Volken, E., Giesche, A. M., &amp; Br¨nnimann, S. o (2015). Colloidal meteorological processes in the formation of precipitation. Meteorologische Zeitschrift , 24 (4), 443-454.

Fujita, K., &amp; Abe, O. (2006). Stable isotopes in daily precipitation at dome fuji, east antarctica. Geophysical research letters , 33 (18).

Genthon, C., Krinner, G., &amp; Castebrunet, H. (2009). Antarctic precipitation and climate-change predictions: horizontal resolution and margin vs plateau issues. Annals of Glaciology , 50 (50), 55-60.

Gorodetskaya, I. V., Tsukernik, M., Claes, K., Ralph, M. F., Neff, W. D., &amp;

Van Lipzig, N. P. (2014). The role of atmospheric rivers in anomalous snow accumulation in east antarctica. Geophysical Research Letters , 41 (17), 61996206.

Grazioli, J., Genthon, C., Boudevillain, B., Duran-Alarcon, C., Del Guasta, M., Jean-Baptiste, M., &amp; Berne, A. (2017). dumont d'urville, ad´lie land, east antarctica. e The Cryosphere , 11 (4), 1797.

Measurements of precipitation in Grazioli, J., Madeleine, J.-B., Gall´e, H., Forbes, R. M., Genthon, C., Krinner, G., &amp; e Berne, A. (2017). Katabatic winds diminish precipitation contribution to the antarctic ice mass balance. Proceedings of the National Academy of Sciences , 114 (41), 10858-10863.

Greene, C. A., Gwyther, D. E., &amp; Blankenship, D. D. (2017). Antarctic mapping tools for matlab. Computers &amp; Geosciences , 104 , 151-157.

Howat, I. M., Porter, C., Smith, B. E., Noh, M.-J., &amp; Morin, P. (2019). The reference elevation model of antarctica. The Cryosphere , 13 (2), 665-674.

King, J. C., &amp; Turner, J. (2007). Antarctic meteorology and climatology . Cambridge University Press.

Knuth, S. L., Tripoli, G. J., Thom, J. E., &amp; Weidner, G. A. (2010). The influence of blowing snow and precipitation on snow depth change across the ross ice shelf and ross sea regions of antarctica. Journal of Applied Meteorology and Climatology , 49 (6), 1306-1321.

Kulie, M. S., Milani, L., Wood, N. B., Tushaus, S. A., Bennartz, R., &amp; LEcuyer, A shallow cumuliform snowfall census using spaceborne radar.

T. S. (2016). Journal of Hydrometeorology , 17 (4), 1261-1279.

Lachlan-Cope, T., Ladkin, R., Turner, J., &amp; Davison, P. (2001). Observations of cloud and precipitation particles on the avery plateau, antarctic peninsula. Antarctic Science , 13 (3), 339-348.

Lawson, R. P., Baker, B. A., Zmarzly, P., OConnor, D., Mo, Q., Gayet, J.-F., &amp; Shcherbakov, V. (2006). Microphysical and optical properties of atmospheric ice crystals at south pole station. Journal of Applied Meteorology and Climatology , 45 (11), 1505-1524.

Lemonnier, F., Madeleine, J., Claud, C., Genthon, C., Dur´n-Alarc´n, C., Palerme, a o C., . . . others (2019). Evaluation of cloudsat snowfall rate profiles by a comparison with in-situ micro rain radars observations in east antarctica. The Cryosphere, https://doi. org/10.5194/tc-2018-236 .

Lenaerts, J. T., Fyke, J., &amp; Medley, B. (2018). The signature of ozone depletion in recent antarctic precipitation change: A study with the community earth system model. Geophysical Research Letters , 45 (23), 12-931.

Liu, C., Zipser, E. J., Cecil, D. J., Nesbitt, S. W., &amp; Sherwood, S. (2008). A cloud and precipitation feature database from nine years of trmm observations. Journal of Applied Meteorology and Climatology , 47 (10), 2712-2728.

Maahn, M., Burgard, C., Crewell, S., Gorodetskaya, I. V., Kneifel, S., Lhermitte, S., . . . van Lipzig, N. P. (2014). How does the spaceborne radar blind zone affect derived surface snowfall statistics in polar regions? Journal of Geophysical Research: Atmospheres , 119 (24), 13-604.

Milani, L., Kulie, M. S., Casella, D., Dietrich, S., LEcuyer, T. S., Panegrossi, G., . . . Wood, N. B. (2018). Cloudsat snowfall estimates over antarctica and the southern ocean: An assessment of independent retrieval methodologies and multi-year snowfall analysis. Atmospheric research , 213 , 121-135.

Nicolas, J. P., &amp; Bromwich, D. H. (2011). Climate of west antarctica and influence of marine air intrusions. Journal of Climate , 24 (1), 49-67.

Palerme, C., Claud, C., Wood, N., L'Ecuyer, T., &amp; Genthon, C. (2019). How does ground clutter affect cloudsat snowfall retrievals over ice sheets? IEEE Geoscience And Remote Sensing Letters, VOL 16, 342-346 .

Palerme, C., Genthon, C., Claud, C., Kay, J. E., Wood, N. B., &amp; L Ecuyer, T.

(2017). Evaluation of current and projected antarctic precipitation in cmip5 models. Climate Dynamics , 48 (1-2), 225-239.

Palerme, C., Kay, J., Genthon, C., L'Ecuyer, T., Wood, N., &amp; Claud, C. (2014). How much snow falls on the antarctic ice sheet? The Cryosphere , 8 (4), 15771587.

Pritchard, H., Ligtenberg, S., Fricker, H., Vaughan, D., Van den Broeke, M., &amp; Padman, L. (2012). Antarctic ice-sheet loss driven by basal melting of ice shelves. Nature , 484 (7395), 502.

Rignot, E., Mouginot, J., Scheuchl, B., van den Broeke, M., van Wessem, M. J., &amp; Morlighem, M. (2019). Four decades of antarctic ice sheet mass balance from 1979-2017. Proceedings of the National Academy of Sciences , 116 (4), 1095-1103.

Rodgers, C. D. (2000). Inverse methods for atmospheric sounding: theory and practice (Vol. 2). World scientific.

Schweiger, A. J., Lindsay, R. W., Vavrus, S., &amp; Francis, J. A. (2008). Relationships between arctic sea ice and clouds during autumn. Journal of Climate , 21 (18), 4799-4810.

Shepherd, A., Ivins, E., Rignot, E., Smith, B., Van Den Broeke, M., Velicogna, I., . . . others (2018). Mass balance of the antarctic ice sheet from 1992 to 2017. Nature , 558 , 219-222.

Shepherd, A., Ivins, E. R., Geruo, A., Barletta, V. R., Bentley, M. J., Bettadpur, S., . . . others (2012). A reconciled estimate of ice-sheet mass balance. Science , 338 (6111), 1183-1189.

Souverijns, N., Gossart, A., Lhermitte, S., Gorodetskaya, I. V., Grazioli, J.,

Berne, A., . . . van Lipzig, N. P. M. (2018). Evaluation of the cloudsat surface snowfall product over antarctica using ground-based precipitation radars. The Cryosphere , 12 (12), 3775-3789. Retrieved from https://

www.the-cryosphere.net/12/3775/2018/ doi: 10.5194/tc-12-3775-2018 Stephens, G. L., &amp; Ellis, T. D. (2008). Controls of global-mean precipitation increases in global warming gcm experiments. Journal of Climate , 21 (23), 6141-

6155. Turner, J., Phillips, T., Thamban, M., Rahaman, W., Marshall, G. J., Wille, J. D., . . . others (2019). The dominant role of extreme precipitation events in antarctic snowfall variability. Geophysical Research Letters , 46 (6), 3502-3511. Van de Berg, W., Van den Broeke, M., &amp; Van Meijgaard, E. (2007). Heat budget of the east antarctic lower atmosphere derived from a regional atmospheric climate model. Journal of Geophysical Research: Atmospheres , 112 (D23). Wood, N. B. (2011). Estimation of snow microphysical properties with application to millimeter-wavelength radar retrievals for snowfall rate (phd thesis). Wood, N. B., L'Ecuyer, T. S., Heymsfield, A. J., Stephens, G. L., Hudak, D. R., &amp; Rodriguez, P. (2014). Estimating snow microphysical properties using collocated multisensor observations. Journal of Geophysical Research: Atmospheres , 119 (14), 8941-8961.
Florentin Lemonnier Jean-Baptiste Madeleine Chantal Claud Cyril Palerme Christophe Genthon Tristan L'Ecuyer Norman B. Wood
CloudSat-inferred vertical structure of precipitation over the Antarctic continent
Current global warming is causing significant changes in snowfall in polar
regions, directly impacting the mass balance of the ice caps. The only water
supply in Antarctica, precipitation, is poorly estimated from surface
measurements. The onboard cloud-profiling radar of the CloudSat satellite
provided the first real opportunity to estimate precipitation at continental
scale. Based on CloudSat observations, we propose to explore the vertical
structure of precipitation in Antarctica over the 2007-2010 period. A first
division of this dataset following a topographical approach (continent versus
peripheral regions, with a 2250m topographical criterion) shows a high
precipitation rate (275mm/yr at 1200meters above ground level) with low
relative seasonal variation (+/-11%) over the peripheral areas. Over the
plateau, the precipitation rate is low (34mm/yr at 1200m.a.g.l.) with a much
larger relative seasonal variation (+/-143%). A second study that follows a
geographical division highlights the average vertical structure of
precipitation and temperature depending on the regions and their interactions
with topography. In particular, over ice-shelves, we see a strong dependence of
the distribution of precipitation on the sea-ice coverage. Finally, the
relationship between precipitation and temperature is analyzed and compared
with a simple analytical relationship. This study highlights that precipitation
is largely dependent on the advection of air masses along the topographic
slopes with an average vertical wind of 0.02m/s. This provides new diagnostics
to evaluate climate models with a three-dimensional approach of the atmospheric
structure of precipitation.