Drop Size Distribution Variability in Central Argentina during RELAMPAGO-CACTI

The Remote sensing of Electrification, Lightning, And Meso-scale/micro-scale Processes with Adaptive Ground Observations (RELAMPAGO) and the Cloud, Aerosol, and Complex Terrain Interactions Experiment Proposal (CACTI) field campaigns provided an unprecedented thirteendisdrometer dataset in Central Argentina during the Intensive (IOP, 15 November to 15 December 2018) and Extended (EOP, 15 October 2018 to 30 April 2019) Observational Periods. The drop size distribution (DSD) parameters and their variability were analyzed across the region of interest, which was divided into three subregions characterized by the differing proximity to the Sierras de Córdoba (SDC), in order to assess the impact of complex terrain on the DSD parameters. A rigorous quality control of the data was first performed. The frequency distributions of DSD-derived parameters were analyzed, including the normalized intercept parameter (logNw), the mean volume diameter (D0), the mean mass diameter (Dm), the shape parameter (μ), the liquid water content (LWC), and the rain rate (R). The region closest to the SDC presented higher values of logNw, lower D0, and higher μ, while the opposite occurred in the farthest region, i.e., the concentration of small drops decreased while the concentration of bigger drops increased with the distance to the east of the SDC. Furthermore, the region closest to the SDC showed a bimodal distribution of D0: the lower values of D0 were associated with higher values of logNw and were found more frequently during the afternoon, while the higher D0 were associated with lower logNw and occurred more frequently during the night. The data were analyzed in comparison to the statistical analysis of Dolan et al. 2018 and sorted according to the classification proposed in the cited study. The logNw-D0 and LWC-D0 two-dimensional distributions allowed further discussion around the applicability of other mid-latitude and global precipitation classification schemes (startiform/convection) in the region of interest. Finally, three precipitation case studies were analyzed with supporting polarimetric radar data in order to relate the DSD characteristics to the precipitation type and the microphysical processes involved in each case.


Introduction
The characterization of drop size distributions (DSDs) is fundamental for both the remote sensing of precipitation and its representation in numerical models. It is directly related to the radar reflectivity (Z) and the rain rate (R), which makes it an essential parameter for quantitative precipitation estimates [1][2][3][4]. The parameterizations of microphysical processes in numerical models also need to assume particle size distributions, generally represented by a gamma distribution with one to three prognostic parameters [5][6][7], among others.
The DSD varies with the different microphysics and precipitation regimes across the world [8]. The classification of convective/stratiform precipitation from DSD parameters has been the subject of many studies, and different classification schemes have been proposed, for instance setting thresholds on variables such as the standard deviation of R [9] or establishing separations in the two-parameter space of logN w -D 0 , related to the number concentration and mean size of drops, respectively [10,11] (BR09, TH15). Recently, Ref. [12] (D18) applied principal component analysis (PCA) to a global DSD dataset, revealing the main modes of variability of the DSD in high, mid-, and low latitudes. Based on the PCA analysis, the dataset was divided into six different groups with similar DSD characteristics related to the dominant microphysical processes in each one of them.
Previous works have analyzed DSDs from different regions and climatic regimes and identified the characteristics of continental, maritime, and tropical environments [9][10][11]13,14]. Distinctive features of orographic DSDs, such as the elevated concentration of small drops compared to the DSDs found in the adjacent plains, were found especially for low to moderate R [15][16][17][18]. In addition, several studies have analyzed the DSD in severe weather situations [16] and mesoscale convective systems in various regions of the world [19][20][21].
Deep moist convection over Southeastern South America (SESA) (blue box in Figure 1a) has awoken interest among the scientific community on a global scale, since it produces extreme storms of great vertical development [22], which are usually associated with high-impact weather [23][24][25] like giant hail [26], lightning [27,28], intense wind gusts, and, occasionally, tornadoes [29,30]. Mulholland et al. (2018) [31] identified a large number of discrete and multicellular convective events close to Sierras de Córdoba (SDC) (Figure 1b) using ground radar scanning. Most of the extreme convection observed in the area initiates on the high central terrain and the northern tip of the SDC [32]. Due to the environmental conditions, eastward propagation is preferentially observed. The most active months for deep convective storms near the SDC are between November and February or during the transition from austral spring into summer. Discrete convective modes tend to be favored earlier in the austral spring season, whereas multicellular convective modes are skewed toward later in the spring season and into the summer season. The field campaigns Remote sensing of Electrification, Lightning, And Meso-scale/ micro-scale Processes with Adaptive Ground Observations (RELAMPAGO) [33] and the Cloud, Aerosol and Complex Terrain Interactions Experiment Proposal (CACTI) [34] took place between 15 October 2018 and 30 April 2019, leeward of the Andes Mountains in Central Argentina (http://relampago-cacti.org (accessed on 1 October 2019)). RELAMPAGO consisted of intensive operations between 1 November and 16 December 2018 and an extended period of hydro-meteorological observations between May 2018 and May 2019. CACTI studied the atmospheric boundary layer, aerosols, and clouds in a fixed location in Yacanto de Calamuchita, next to the highest peaks of the SDC, between October 2018 and April 2019. Both experiments have built a unique and unprecedented dataset in the region. In the context of the RELAMPAGO-CACTI field campaigns, the present study compiled a dataset including thirteen disdrometers in different locations of Central Argentina, covering over seven months of data.
The main objective of this study was to present a robust characterization of the DSD and its variability in Central Argentina, making use of the exceptional amount of data that was gathered during RELAMPAGO-CACTI field campaigns, and makes a contribution to the understanding of physical precipitation processes during the warm season in the mid-latitudes.
Section 2 includes a description of the dataset and the methodologies for the quality control and the estimation of the DSD parameters. In Section 3.1, the distribution of the different DSD parameters is analyzed as a function of the distance to the SDC, and the data are sorted according to the D18 classification. Section 3.2 presents three case studies, which, with supporting polarimetric radar data, relate the DSD characteristics described in the previous section with the precipitation regime and microphysics involved. Finally, Section 4 presents the summary and conclusions of the study.

Materials and Methods
The DSD sample compiles thirteen OTT Parsivel disdrometer datasets from Central Argentina. Twelve disdrometers were installed during RELAMPAGO-CACTI, with the support of the National Science Foundation (NSF), NASA Global Precipitation Measurement-Ground Validation (GPM-GV), and the Department of Energy of the United States (DOE), and one belongs to the Sistema Nacional de Radares Meteorológicos (SINARAME). The RE-LAMPAGO data were downloaded from the Earth Observing Laboratory (EOL-NCAR) website (https://www.eol.ucar.edu (accessed on 1 October 2019)); the CACTI data were downloaded from the Atmospheric Radiation Measurement program (ARM-DOE) website (https://adc.arm.gov/discovery (accessed on 1 October 2019)); the GPM disdrometer data were provided by the Global Precipitation Measurement Ground Validation Office (https://trmm-fc.gsfc.nasa.gov/ftp/pub/tokay/relampago/ (accessed on 1 October 2019)); and the SINARAME data were available under a cooperation agreement between the Secretaría de Infraestructura y Política Hídrica de la Nación and the Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires. The thirteen disdrometers were located between 31.3 • S and 33.2 • S and between 60.9 • W and 64.7 • W, as shown in Figure 1b, which also shows the three proposed regions in which the study area was subdivided in order to analyze the impact of complex terrain and the proximity to the SDC on the DSD parameters. Table 1 details information about the sites, project, location, height above mean sea level (AMSL), time frame, sampling time, time between samples, and the length of each dataset. During the field experiment, five C-band radars were available: CSAPR (DOE-CACTI), CSU (NSF-RELAMPAGO), COW (NSF-RELAMPAGO), RMA1 (SINARAME), and PAR (Instituto Nacional de Tecnología Agropecuaria (INTA)). Their locations and 120 km range are shown in Figure 1. However, only CSU and PAR were used in the present study, as a visual aid in case studies to assist with event description. For further details about the radar capabilities, see Nesbitt et al. 2021 [33].
The OTT Parsivel disdrometer [35][36][37] is an optical sensor that measures the attenuation generated by the precipitation particles that go through a laser beam during a sampling time, which is usually one minute. The amplitude of the attenuation is related to the particle size, and the length of the signal is related to its fall speed; therefore, the raw data provided by the disdrometer are a two-dimensional histogram of drop sizes D i (i = 1, . . . , 32) and fall speeds V j (j = 1, . . . , 32). Then, the DSD can be computed summing over j.
Before calculating the DSD, a quality control was performed on the data. First, it was necessary to filter out possible solid precipitation particles, non-meteorological signals such as spiders, spider webs, or insects, and errors related to border effects, such as drops that were captured only partially and were therefore classified as smaller drops with anomalously high fall speeds, or the opposite case where two or more drops fall through the sampling area simultaneously and were therefore classified as a single bigger drop with an anomalously low fall speed [16,38]. This was achieved following the methodology used in many previous studies [15,[39][40][41][42], which consisted of removing the drops whose measured fall speed V j differed by more than 50% from the expected terminal velocity for a drop of its size, given by the following expression, which is a polynomial fit to the experimental results obtained by Gunn and Kinzer (1949) [43]: Equation (1) was developed considering sea level pressure conditions. Foote and Du Toit (1969) [44] showed that, between 1013 and 900 mb, terminal raindrop fall speeds do not change significantly due to air density effects, and all disdrometers are located approximately within that pressure range. That is why air density effects were not considered. Figure 2 shows the normalized frequency of each diameter (D i ) and fall speed (V j ) for the whole-sample thirteen-disdrometer dataset before quality control. The gray curve represents the expected terminal velocity given by Equation (1), and the shaded area indicates the ±50% thresholds. Although the maximum frequencies fell mostly inside the shaded area, the unfiltered sample presented above-zero frequencies for almost all sizes and speeds in the spectrum (zero values were masked), which proved the importance of this step. The largest portion of the drops removed presented diameters between 0.3 mm and 1.6 mm and fall speeds lower than half the terminal velocity expected for raindrops of those sizes, and they were non-precipitation particles, small drizzle particles that may be affected by turbulence, or particles that go through the sampling area almost horizontally due to strong winds.
Following D18 and TH15, only 1 min DSD data with at least 100 total drops and R > 0.05 mm h −1 during at least 3 consecutive minutes (except for FDC, which has 1 observation every 10 min) were analyzed. This prevented a large number of very small samples from skewing the analysis towards an underestimation of the DSD parameters [45,46]. This step removed 56% of the data points, but only 6% of the total rainfall. The 1 min DSD was then computed from the raw data according to the following expression: where S e f f (D i ) is the effective sampling area (m 2 ), t is the sampling time (60 s in our dataset), ∆D i is the bin width (mm), n ij is the number of drops of size D i and fall speed V j , and v D i (m s −1 ) is the terminal velocity given by Equation (1). The units of N(D) are mm −1 m −3 , so that multiplying by ∆D i (mm) yields the number of drops of size between accounts for the border effects previously mentioned, where L = 180 mm and W = 30 mm were the length and width of the laser beam, respectively. v D i was used instead of the measured fall speed V j because the Parsivel disdrometer underestimated the speed of drops, particularly for the smaller diameters, around 1.1 mm [37], and this could have a significant effect on the estimated DSD and its parameters.
The liquid water content (LWC) can be obtained by integrating N(D) as follows: where ρ w is the water density (g mm −3 ). R is also computed integrating N(D): The four-parameter modified gamma distribution is often used to represent size distributions [47] of cloud drops and precipitation particles, both liquid and solid. The exponential distribution [48], the power law [49], and the normalized gamma with three free parameters [50] are special cases of the former. The last one is often the best fit for rain DSDs [51,52] and is a function of the normalized intercept parameter N w (mm −1 m −3 ), the mean mass diameter D m (mm), and the shape parameter µ following: where f (µ) contains the gamma function: D m is the first moment of the mass spectrum m(D) = πρ w 6×10 3 N(D)D 3 and can be written as the ratio between the fourth and third moment of N(D): and N w is the intercept parameter of an exponential function (µ = 0) with the same LWC: The normalized gamma distribution can also be expressed in terms of the mean volume diameter D 0 , which is the diameter that accumulates half LWC: and is related to D m and N w through the following equations: The shape parameter µ was obtained by minimizing the mean-squared error between the measured DSD and the normalized gamma function defined by the estimated N w and D m , with µ varying between −10 and 20. The gamma fit and associated parameters were calculated using the library PyDSD (https://github.com/josephhardinee/PyDSD (accessed on 1 October 2019)) [53].
The normalized gamma fit was used here in order to compare our results with the large number of previous studies that also take this approach, although it may not fully capture the true spectrum of DSDs [54]. The gamma fitted DSDs were further restricted to only those with µ values in the range from −4 to 15, to ensure the gamma fit was a reasonable assumption, following D18. This condition removed 39% of observations (after the filtering described earlier) and 7% of the total rainfall, and the proportion increased with the proximity to the SDC. The removed samples presented much lower values of D m and R than the remaining samples, which might indicate drizzle or very small drops going through the sampling area repeatedly because of turbulent motion [55]. These conditions generated a database of 72,395 raining minutes representative of a Southern Hemisphere warm season in mid latitudes.

Characterization of the DSD during RELAMPAGO-CACTI
In order to evaluate how the proximity to the complex terrain in the SDC affects the DSD parameters, the study area was subdivided into three regions, as shown in Figure 1b: Region A is the closest to the SDC and includes the disdrometers APU, FDC, CCT, 014, and 015, all located above 400 m AMSL; Region B includes disdrometers 001 and 015, both between 200 m and 400 m AMSL; and Region C includes disdrometers 003, 004, 005, 007, 008, and 010, all of them below 200 m AMSL and well within the plains, approximately 350 km from the SDC. It is worth mentioning that Region B had only 8616 data points, while Regions A and B had 27,668 and 36,111, respectively, so this should be taken into consideration in the following analysis. Figure 3 shows the violin plots of different DSD parameters for each region. The dashed lines indicate the median, and the dotted lines indicate the first and third quartiles (P25 and P75, respectively). The data were normalized so that each violin had the same area, which made it easier to compare samples of different sizes. Additionally, Table 2 lists the mean value, mode, standard deviation, minimum, maximum, and fifth (P05) and ninety-fifth percentiles (P95) for each DSD parameter at each region.  From Region A to Region C, the median, P25, and P75 of D m adopted larger values (Figure 3a), as well as the mean and P95 ( Table 2). The same happened with D 0 (Figure 3b), which, in addition, presented a bimodal distribution in Region A. The maximum diameter D max was the diameter of the biggest drop in each individual sample, so it was the only parameter that was directly measured by the disdrometer, which classified the drops in size bins of varying width. That is the reason why the violin plots in Figure 3c exhibit irregular borders. In Regions B and C, P25 and the median were higher than in Region A, and in Region C, P75 was higher than in Regions A and B. Moreover, Table 2 shows that the mean of D max grew from A to C. For the three parameters described, Welch's t-test (t-test for unequal variances and/or unequal sample sizes) indicated that the mean in Region B was larger than in Region A and that the mean in Region C was larger than the mean in Regions A and B, with a 99% confidence level. This means that the probability of finding DSDs with bigger drops increased eastward from SDC.
It should be noted that Welch's t-test assumes that the means of different samples from a population are normally distributed, and usually, the sample means are normally distributed for sample sizes above 30, which is why the test is reliable despite the fact that the DSD parameters were not normally distributed.
The median and interquartile range of logN w decreased from A to C (Figure 3d), and the mean and P05 and P95 behaved similarly ( Table 2). Welch's t-test indicated that the mean in A was larger than the mean in B and C and that the mean in B was larger than in C, with a 99% confidence level. Therefore, the parameters analyzed pointed out that the DSDs with higher number concentrations and smaller drops were more frequent near the SDC, and with the growing distance to the east, the number concentrations decreased while the concentration of bigger drops increased. However, Region A presented the broadest range in all four parameters analyzed, which means that although the DSDs with a higher concentration of small drops were more frequent, it was also possible to find DSDs whose characteristics resembled those of Region C.
Regarding µ, there were no noticeable differences between Regions A and B. However, the difference between these two regions and Region C, whose median and interquartile range (Figure 3e), as well as the mean and P95 (Table 2) adopted smaller values, was evident. Welch's t-test indicated that the mean in Region A was larger than in Region C and that the mean in Region B was also larger than the mean in Region C, with a 99% confidence level. Smaller µ values were related to broader distributions with higher concentrations of bigger drops, while larger µ values were related to narrower distributions. Hence, the spatial variability of µ was consistent with the behavior of the parameters analyzed previously.
The median and interquartile range of LWC adopted larger values from A to C (Figure 3f). The mean, mode, P05, P95, and maximum value exhibited the same behavior ( Table 2). Figure 3g shows the violin plots of R. The median and interquartile range adopted larger values from A to C, and the dispersion also increased. The mean, mode, and P95 showed the same behavior (Table 2). Both for LWC and R, Welch's t-test indicated that the mean in Region A was smaller than in Regions B and C and that the mean in Region B was smaller than in Region C, with a 99% confidence level. That is, the probability of finding larger values of LWC and R increased with the distance to the east of the SDC, which was consistent with the increasing concentration of larger drops.
D18 proposed a classification of the DSD based on the principal component analysis (PCA) of a global dataset. The input variables for the PCA analysis were: D m , logN w , σ m (square root of the second moment of the mass spectrum-Equation (8)), logR, logLWC, and logN t (total number of drops). First, the data were normalized by subtracting the mean and dividing by the standard deviation, then the covariance matrix was constructed with the input variables, and the orthogonal empirical functions (EOFs) were computed. By setting thresholds (±1.5) on the first two PCs, they defined six groups with independent characteristics. The physical processes in the six groups were explained through supporting radar observations. Group 1 was related to convective precipitation, Group 2 to stratiform precipitation, Group 3 warm rain showers, Group 4 heavy stratiform precipitation, Group 5 (characteristic of the low-latitudes) robust warm rain processes, and Group 6 (more frequent in the mid-latitudes) ice-based convection, although it may also appear in heavy stratiform rain. Readers are encouraged to read D18 for further details.
The classification shown in Figure 4 was obtained by projecting the data of the present study onto the EOFs derived from the D18 world dataset and applying similar thresholds on the first two PCs. A slightly higher threshold was selected here (±2.3) in order to avoid overlapping between modes. The scatterplots in Figure 4 correspond to the entire dataset (i.e., Regions A, B, and C). The dotted and solid lines show the convective (above)/stratiform (below) separations proposed by BR09 and TH15, respectively. The latter was proposed specifically for oceanic regions of tropical latitudes, so it was not applicable here, and it is shown for reference. Note that the groups resulting from the PCA analysis in D18 were consistent with the BR09 and TH15 separations. Table 3 shows the percentage of observations in each group for the whole dataset and for each region. Sixty-four percent of the data were below the imposed thresholds, falling into the "Ambiguous" category, which means they resembled multiple EOFs. Group 4 had the largest number of observations (18%), and it was followed by Group 6 (9%) associated with mixed-phase storms. The proportion of observations in these two groups increased from Region A to Region C. Six percent of the data were classified as Group 2 with the proportion decreasing eastward of the SDC. Finally, Group 1 was the least frequent, with 4% of the data falling in the category and the proportion increasing eastward of the SDC.   Figure 5a-c shows the 2D histograms of logN w -D 0 for each region. Region A presented two local maxima: one centered approximately at D 0 = 0.9 mm and logN w = 3.4 and another one at D 0 = 1.3 mm and logN w = 3.1. The maximum in Region B was centered at D 0 = 1.3 and logN w = 3.1, and the maximum in Region C was centered at D 0 = 1.4 mm and logN w = 3.0. Therefore, from Region A to Region C, the maximum frequencies adopted smaller values of logN w and larger D 0 , which was consistent with our findings so far. Figure 5 additionally shows that the frequencies in the area of the diagram associated with Group 6 increased from A to C, as well as the frequencies above the dotted line, which were associated with convective DSDs. In the histograms of LWC-D 0 (Figure 5d-f), the shift in the position of the maximum was less apparent. However, there was an increase in the frequency of DSDs with LWC above 0.1 g m −3 and D 0 above 2 mm, which was the area of the diagram associated with Group 6. Region A presented again a bimodal distribution, with a local maximum centered approximately at D 0 = 0.9 mm and LWC = 0.03 g m −3 and another local maximum centered approximately at D 0 = 1.3 mm and LWC = 0.05 g m −3 .
The histograms of D as a function of time (UTC) ( Figure 6) show the diurnal cycle of the DSD for Regions A, B, and C separately. The three regions presented a maximum of observations at night (LT is UTC-4hs), which is characteristic of the study region [56][57][58]. However, the maximum frequencies (above 0.6%) in Region A were centered between 0400 and 0600 UTC; in Region B, the maximum frequencies were centered between 0400 and 0800 UTC; and in Region C, the maximum frequencies were centered between 0600 and 0900 UTC. This might be related to a zonal propagation of the systems. Region A presented a second local maximum in the afternoon, centered approximately at 1800 UTC (i.e., 1400 LT), which was not evident in the other regions. Figure 7 shows the 2D histograms of logN w -D 0 and LWC-D 0 only for Region A and every three-hour interval. The night maximum was associated with lower logN w (maximum centered approximately around three) and larger D 0 (maximum centered approximately around 1.6 mm) (Figure 7b,c), while the afternoon maximum (1500 to 1800 and 1800 to 2100 UTC) was associated with higher logN w (maximum centered approximately at 3.5) and smaller D 0 values (maximum centered approximately at 1 mm) (Figure 9f,g). LWC did not vary significantly, although from 0300 to 0600 UTC (Figure 9b) and from 1800 to 2100 UTC (Figure 9g), there was a significant increase of DSDs with moderate to large D 0 and large LWC, as well as higher frequencies above the dotted line in the logN w -D 0 diagrams, which were associated with convection. From 0900 to 1500 UTC and from 2100 to 0300 UTC, the histograms showed a transition between the described regimes.

Case Studies
In this section, we analyzed three case studies with supporting radar data from radars CSU and PAR, relating the DSD characteristics with the type of precipitation and the inferred microphysical processes behind this. The first two cases are examples of the two regimes observed in Region A, and the third case is representative of Region C. The radar data were processed using the PyART library (https://arm-doe.github.io/pyart/ (accessed on 1 October 2019)) [59] and used only as a visual aid to assist with event description.

Case 1: 26 January 2019
This case study was measured by the CCT disdrometer in Region A, which had coincident observations from the CSU radar (see Figure 1b). The sample consisted of 201 consecutive minutes.
At 0200 UTC, a line of convective cells passed over the CCT disdrometer, as shown by the CSU radar (Figure 8g), with maximum reflectivities of around 50 dBZ. The peaks of R, D 0 , and logN w at the beginning of the time series in Figure 8a-c mark the passing of the squall line. The low values of µ were consistent with the wide spectrum of sizes (shading in Figure 8a) and high concentration of big drops. The color of the markers indicated the corresponding group following the classification proposed by D18, here obtained after projecting the data onto the D18 world EOFs as described in the previous section. During this initial convective phase, DSDs were classified as Group 1 (convection) and Group 6 (intense ice-based convection). This was consistent with the presence of big drops, which were probably formed by the melting of ice particles. All data points fell in the convective side of BR09 classification (Figure 8e), with moderate to high logN w and D 0 and high LWC (Figure 8f). Following the initial convective phase, a transition phase followed during which µ presented a large variability, logN w oscillated around three, and D 0 and R decreased. Most DSDs in this phase did not belong to any of the groups proposed by D18 and were classified as "Ambiguous".
Finally, at 0500 UTC, a stratiform system passed over the disdrometer. The PPI in Figure 8h shows a widespread system, with Z H around 30 dBZover the disdrometer and the bright band signature following approximately the 120 km range, with Z H over 40 dBZ. As shown in Figure 8a-d, R became stable, and so did N(D); D 0 increased slightly; logN w decreased slightly; and µ oscillated between approximately two and five, somewhat higher values than those during the convective phase, which was consistent with the narrower size distributions (shading in Figure 8a). Almost all data points during this final stratiform phase were classified as Group 4 (intense stratiform), which was consistent with the moderate values of D 0 and logN w and the strong bright band signature. The BR09 separation line also classified this phase as stratiform (Figure 8e), and LWC was significantly lower than in the previous convective phase (Figure 8f), which was consistent with the lower values of logN w and D 0 .
Both D18 and BR09 classifications were successful in this case study, with the initial convective phase well classified as convective and the final stratiform phase well classified as stratiform in both cases. Clearly, TH15 classification did not apply in this case and would not be applicable in the following cases either, since, as was mentioned before, it was specifically developed for the oceanic regions of tropical latitudes. It is included however in all the logN w -D 0 diagrams for visual reference.

Case 2: 17 January 2019
Case 2 was a case of weak, shallow, and unorganized convection in Region A, measured by CCT with coincident observations from the CSU radar. The sample consisted of 337 consecutive minutes. Figure 9g,h shows a number of isolated convective cells, with reflectivities no higher than 30-35 dBZ, as observed by the CSU radar. The time series in Figure 9a-d show great variability, consistent with the successive passing of convective cells over the disdrometer. The size spectrum of drops (shading in Figure 9a) was much narrower than in the previous case; D 0 was significantly smaller, oscillating between 1 and 2 mm, and logN w was significantly higher, oscillating between three and four.
Most data points were classified as ambiguous according to the D18 classification applied to the present database. Only a few were classified as Group 2 (stratiform) and Group 4 (intense stratiform). In the logN w -D 0 space, all data points fell in the stratiform area of the BR09 separation ( Figure 9e). Neither classification worked in this case. However, the characteristics of the event were compatible to those of Group 3 in the D18 classification, associated with weak echoes and shallow convection. It is possible that the method failed in this case because of the underrepresentation of subtropical latitudes and mountainous regions in the D18 world dataset, as well as the fact that all locations in that dataset belonged to the Northern Hemisphere. It is also possible that there were additional modes of variability of the DSD that were better explained by higher order EOFs and therefore were not described by the proposed classification (D18).
The cases analyzed so far were examples of the two regimes identified in Region A: Case 1, which took place at night, was characterized by lower N w and larger D 0 , while Case 2, which took place during the afternoon, presented significantly higher values of N w and much smaller D 0 .

Case 3: 4 March 2019
Case 3 was similar to Case 1 in that both consisted of an initial convective phase and a final stratiform phase. It was measured by disdrometer 005, which was within the 240 km range of the PAR radar (Figure 1b). The sample consisted of 449 consecutive minutes.
The initial convective phase in Case 3 extended from the beginning of the time series up to approximately 0500 UTC ( Figure 10), with elevated R, D 0 around 3 mm, logN w around 3, and low µ values, consistent with the wide spectrum of sizes (shading in Figure 10a). The PPI taken by the PAR radar at 0330 UTC (Figure 10g) showed an intense squall line, with maximum reflectivities around 50 dBZ. Most data points of the initial phase were classified as Group 6 (intense ice-based convection) in the D18 classification and some as Group 1 (convection), and all of them fell in the convective area of the BR09 separation (Figure 10e), with the highest LWC of the time series (Figure 10f). After the initial convective phase, a transition period followed during which R decreased significantly, as well as D 0 . LogN w remained approximately constant, and µ increased significantly. Most of those data points were classified as ambiguous and some as Group 4 (stratiform), although there was no bright band signature observed yet.
The stratiform phase took place later, as observed by the PAR radar at 0830 UTC. The PPI in Figure 10h shows a widely spread system with uniform horizontal characteristics and moderate reflectivities of around 30 dBZ above the disdrometer. Figure 10a-d shows a stabilization of R and all other parameters. D 0 oscillated between 2 and 2.5 mm, exceeding 3 mm in some extreme cases. LogN w oscillated around 2.5, and µ oscillated around two, adopting slightly smaller values than during the stratiform phase in Case 1, which was consistent with the much wider size spectrums (shading in Figure 10a). Data points during this final phase were classified as Group 4 (intense stratiform) and Group 6 (intense icebased convection). Although Group 6 was mostly related to convective processes, it might also be associated with very intense stratiform precipitation (D18), which was clearly the case. Most data points in this final phase fell in the stratiform area of the BR09 separation, but some of the points classified as Group 6 crossed the separation, falling in the convective zone. That is why, in this case, the D18 classification worked better than BR09, since the latter would classify as convective some DSDs that were stratiform without doubt.

Discussion
The variability of the different DSD parameters was analyzed as a function of the distance to the SDC. Region A, which was the closest to the SDC, presented the highest number concentrations, smallest mean and maximum diameters, higher values of µ, and low R and LWC, while Region C, which was the farthest to the east of the SDC, exhibited much lower values of N w , significantly larger D 0 , D m , and D max , i.e., higher concentrations of bigger drops also denoted by lower values of µ, and, consequently, the largest values of LWC and R. Region B, located in the middle between A and C, presented intermediate characteristics.
This suggests that the number of drops per unit volume decreased with the distance to the east of the SDC, while the mean size of drops increased, leading to an increase in LWC and R.
Then, the data were projected onto the D18 world EOFs and classified accordingly into Groups 1 to 6. Sixty-four percent of the data were below the imposed thresholds, falling into the "Ambiguous" category. The most frequent Groups were 4 (18%) and 6 (9%), related to intense stratiform rain associated with aggregation and riming processes and intense ice-based convection, respectively. Six percent of the data were classified as Group 2 (stratiform rain associated mainly with vapor deposition processes), and four percent were classified as Group 1 (convection). No data points were classified as Group 3 or Group 5. The latter was expected, since this category was almost non-existent in the mid-latitude dataset in D18 and is found mainly in the low latitudes. The lack of points classified as Group 3, however, might be due to the underrepresentation in the D18 world dataset of mountainous regions, subtropical latitudes, and the Southern Hemisphere. It might also be related to the lower values of N w and the larger D 0 of the RELAMPAGO-CACTI dataset compared to the mid-latitude dataset in D18, which might also explain the fact that most non-ambiguous data points were classified as Group 4 and Group 6 and very few as Group 2 and Group 1.
Two distinct modes of variability were identified in Region A: at night, the DSD presented lower N w and larger D 0 , while in the afternoon, N w increased and D 0 decreased. The elevated number concentrations and smaller mean diameters were compatible with orographic DSDs described in previous literature [17,18], and since this mode of variability was more frequent in the afternoon, it might be related to diurnal heating on the hillsides as a forcing mechanism of convection.
Finally, three case studies were presented. The first two were representative of Region A, one of the night regime (Case 1) and the other one of the afternoon regime (Case 2). The third one was representative of Region C. The aim was to relate the DSD characteristics with the type of precipitation and microphysical processes involved, with the support of radar observations, in the context of the analysis proposed in D18. In Cases 1 and 3, the D18 classification worked fine, with most convective times being classified as Group 6 (intense ice-based convection) and some as Group 1 (convection), most stratiform data points being classified as Group 4 (intense stratiform), and some as Group 6, especially in Case 3 (Region C), where drops were bigger and N w was lower, which might be related to large melted ice particles. In Case 2, however, most points fell into the "Ambiguous" category, and a few were classified as Group 2 and Group 4. Case 2 was an example of weak shallow convection in the afternoon in Region A, with N w higher than in the other cases and smaller D 0 . Radar observations suggested Case 2 was compatible with Group 3 in the D18 classification, although here, logN w was significantly lower, which was why most data points fell in the "Ambiguous" category rather than being classified as Group 3. Again, the misclassification in this case might be explained by the differences between the D18 world dataset and the RELAMPAGO-CACTI dataset, as well as the underrepresentation of the study region in the former.

Conclusions
We presented a robust characterization of the DSD in Central Argentina based on the large amount of data collected during RELAMPAGO-CACTI, dividing the study region into three subregions in order to assess the impact of the proximity to the complex terrain of the Sierras de Córdoba (SDC).
The variability of the different DSD parameters was analyzed as a function of the distance to the SDC. The distribution of the different parameters suggested that the number of drops per unit volume decreased with the distance to the east of the SDC, while the mean size of drops increased, leading to an increase in LWC and R.
Furthermore, the disdrometers located closer to the SDC presented a bimodal distribution of D 0 : the larger D 0 were associated with lower logN w and occurred more frequently during the night, while the smaller D 0 were associated with higher values of logN w , compatible with orographic DSDs, and were found more frequently during the afternoon.
The data were classified following the methodology proposed in Dolan et al. (2018) (D18). The most frequent categories were those related to intense stratiform and intense convection, both associated with ice-based processes; fewer data points were classified as convective and stratiform alone, and no data points were classified as weak convection or convection associated with warm rain processes. The latter was expected, since that mode of variability of the DSD was found mainly at low latitudes. However, the lack of points classified as weak convection (small D 0 and high logN w ), as well as the high frequency of intense stratiform (low logN 0 and moderate D 0 ) and intense convection (low logN w and large D 0 ) might be explained by the larger D 0 and lower logN w of the RELAMPAGO-CACTI dataset compared to the global dataset analyzed in the cited article.
Finally, three case studies were analyzed. The first two corresponded to the two regimes observed in the region closest to the SDC: Case 1 corresponded to the night regime and Case 2 to the afternoon regime. Case 3 was representative of the easternmost region. The events were described with the aid of coincident radar observations. Case 1 and Case 3 consisted of an initial convective phase followed by a stratiform period, which was reflected on the temporal evolution of the DSD parameters, as well as the radar images, and the classification proposed in D18 worked well in both cases. However, Bringi et al.'s (2009) (BR09) separation misclassified some of the stratiform points in Case 3 mainly because of the larger D 0 , characteristic of the study region and associated with the presence of big melted ice particles. Neither the D18 classification, nor the BR09 separation were able to correctly classify the data in Case 2, which consisted of weak shallow convection. According to the D18 classification, most data points were ambiguous, and the BR09 separation classified most data points as stratiform, when the DSD parameters and radar images indicated that it was a convective event.