remote Droughts and Floods in the La Plata Basin in Soil Moisture Data and GRACE

: The mapping and forecasting of droughts and floods is an important potential field of application of global soil moisture and water storage products from satellites and models. Especially when extremes in near-surface soil moisture propagate into extremes in total water storage, agricultural production and water supply can be severely impacted. This study relates soil moisture from the WaterGAP Global Hydrology Model (WGHM) and the satellite sensors Advanced Microwave Scanning Radiometer—Earth Observing System (AMSR-E) and Advanced Scatterometer (ASCAT) to total water storage variations from the satellite gravity mission GRACE. A particular focus is on destructive hydrological extreme events, as listed in the International Disaster Database EM-DAT. Data sets are analyzed via correlation, time shift, and principal component analyses. The study area is the La Plata Basin in South America. The results indicate that most of the soil moisture anomalies are linked to periods of El Niño and La Niña and associated natural disasters. For the La Plata drought of 2008/2009 and the El Niño flooding of 2009/2010, soil moisture serves as an indicator for the later deficit or surplus in total water storage. These hydrological anomalies were strongest in the southern, central, and eastern parts of the basin, but more than one hundred thousand people were also affected in the northwestern part.


Introduction
Mapping and understanding soil moisture dynamics under extreme hydrologic conditions is important for agricultural production, disaster management, and water supply [1].In contrast to conventional in situ measurements, remote sensing techniques have the advantage of mapping variations in soil moisture on a global scale with a frequent revisit time [2].Recent satellite sensors providing information on global surface soil moisture include the Advanced Microwave Scanning Radiometer-Earth Observing System (AMSR-E) onboard the NASA mission AQUA (2002)(2003)(2004)(2005)(2006)(2007)(2008)(2009)(2010)(2011), the Advanced Scatterometer (ASCAT) onboard the EUMETSAT mission METOP (2006-present), and the Microwave Imaging Radiometer using Aperture Synthesis (MIRAS) onboard the ESA mission SMOS (2010-present).Various studies have analyzed how well the soil moisture products of these sensors map hydrological extremes.The use of soil moisture data for agricultural drought monitoring is a major contributor to the investigation of study sites in, e.g., the USA, Africa, Canada, and China [1,[3][4][5][6].In the context of floods, research has focused on the satellite-based monitoring of soil wetness [7].The relationship between soil moisture and runoff [8] and the use of satellite soil moisture data for runoff modelling [9,10] has also been investigated.The growing need for further research on the use of large-scale soil moisture data, specifically for drought monitoring and socioeconomic modelling and forecasting (e.g., addressing issues of drought impact and food security; [11]), has been recently emphasized in a review on the state of the art in large-scale soil moisture monitoring [12].
During severe weather events, extremes in soil moisture can propagate into extremes in total continental water storage (TWS).The satellite gravity mission GRACE (Gravity Recovery And Climate Experiment) allows for mapping changes in TWS during droughts and floods (e.g., [13][14][15][16]).Soil moisture and TWS are closely linked, as soil moisture is an essential component of TWS in addition to surface water, snow, ice, and groundwater.Furthermore, both parameters strongly depend on the difference between inflow (precipitation) and outflow (percolation, drainage, and evapotranspiration) of the soil and its water-holding capacity (as a function of, e.g., soil physical properties).It was shown that the dynamics of GRACE TWS and those of near-surface satellite soil moisture are strongly correlated in many regions worldwide, with correlation coefficients above 0.5 if all data sets are harmoniously processed [17].In high-agreement regions, soil moisture either dominates the TWS signal or simultaneously changes with surface or ground water.
This study aims to better understand the interaction between soil moisture and TWS, specifically under extreme hydrologic conditions.As soil moisture is closely linked to precipitation and evaporation, it is expected to provide information on the initial phase of hydrologic extremes.Variations of TWS are analyzed to identify possible complementary information on water cycle dynamics during hydrological extremes by considering the sum of all water storage components.Consequently, GRACE TWS is not disaggregated into individual TWS components such as soil moisture, groundwater, and surface water (as for example investigated by [18,19]).Instead the integral information contained in TWS is used.
The study area extends over the La Plata Basin in South America (Section 2).Earlier studies have shown that GRACE TWS can be used to analyze large-scale hydrologic events such as droughts and floods in this basin ( [20][21][22][23]) and the adjacent Amazon Basin ( [23][24][25][26][27][28]). Furthermore satellite soil moisture data from AMSR-E have been used to record flooding along the Paraná River in the southeastern part of the La Plata Basin [29].
In this study, the dynamics of soil moisture and TWS are analyzed over the four-year period September 2007-August 2011.The data sets comprise GRACE TWS and three different soil moisture data sets: root zone soil moisture from the WaterGAP Global Hydrology Model (WGHM) and surface soil moisture from the active microwave sensor ASCAT and the passive microwave sensor AMSR-E (Sections 3.1-3.3).An additional data source-the International Disaster Database EM-DAT-is included.It provides information on natural disasters that occurred in the La Plata Basin (Section 3.5).The definition of a natural disaster is based on the number of people affected or killed by an event; moreover, it applies if a state declares an emergency or calls for international assistance.Therefore, the definition is completely independent of any hydrological and meteorological parameter.The EM-DAT list of natural disasters is integrated in this study to focus on hydrologic extreme events that had a significant impact on society.
The analysis is structured as follows.First, surface soil moisture from ASCAT and AMSR-E are scaled with respect to the root zone soil moisture from WGHM (Section 4.1) and the different soil moisture data sets are compared by correlation and time-shift analysis (Sections 4.3 and 4.4).Consistencies and inconsistencies among the different data sets are discussed on the basis of six sub-regions, representing major ecosystems in the La Plata Basin (Section 5.1).
Second, a comparative analysis of soil moisture and GRACE TWS is performed.To this end, the different data sets on soil moisture and TWS are first harmoniously processed to receive comparable data representations (Section 4.2).Relationships of monthly soil moisture, TWS, and precipitation anomalies with destructive extreme events, listed in the International Disaster Database EM-DAT, are analyzed (Section 5.2).Furthermore, the relationship with El Niño and La Niña patterns is discussed on the basis of the Oceanic Niño Index (Section 3.4).A Principal Component Analysis (PCA; Section 4.3) is performed to identify and compare dominant spatial patterns of the various data sets.Finally, interrelations between the dynamics of soil moisture and of TWS are analyzed on the basis of correlation patterns (Section 5.3).

Study Area: La Plata Basin, South America
The study area is the La Plata Basin in South America.With its three main contributing rivers (the Paraguay, Paraná, and Uruguay Rivers), it is the world's fifth largest river basin and covers parts of Argentina, Brazil, Bolivia, Paraguay, and Uruguay [20].The region was selected because the topographic complexity is low despite its large size of approximately 3,100,000 km 2 [21].This is favorable, especially for signals from active microwave systems that are degraded by mountainous terrain.Moreover, previous studies showed that GRACE can monitor extreme weather events in this river basin and AMSR-E data proved to be sensible for flooding events along the Paraná River (References are provided in Section 1).The main land cover classes of the La Plata Basin are shown in Figure 1.The mean annual cycle of precipitation and river runoff varies significantly within the La Plata Basin [30].The following description of the sub-regions outlined in Figure 1 is summarized from [31] (unless otherwise cited), who describe the climatology of the La Plata Basin in great detail.
 Southern cropland (Pampas; Figure 1, sub-region 1) and southern grassland (Campos; Figure 1, sub-region 6): The Pampas and the Campos are two major sub-regions of the major South American grassland, which stretches from the south of Brazil to Argentina, comprising more than 700,000 km 2 .The climate changes from west to east from humid to dry-subhumid [32].The plain contains the most fertile soils of the La Plata Basin.The Campos has continental climate with frequent frost in winter and daily maximum temperature above 36 °C in summer.Rainfall is distributed throughout the year but changes strongly between years [33].The annual mean precipitation in the North (1600 mm) is larger than in the South (1200 mm).The Campos is mainly used for grassland-based livestock production [34].This is different for the Pampas.This region is intensively used by agriculture.The climate is similar to the Chaco with the main precipitation falling from October to March. Western deciduous broadleaf forest (Chaco; Figure 1, sub-region 2): The Chaco is a subtropical area and located west of the Paraguay River and east of the Andes [35].It is the warmest region in the La Plata Basin.In summer, the daily maximum temperature often exceeds 45°.The rainy season usually takes place between September and April.The semi-arid plains of the Chaco are dominated by dry woodland, which (approaching to the east) transform into more open savannas (Figure 1). Wetland Pantanal (Figure 1, sub-region 3): The Pantanal is one of the world's largest wetlands with enormous biodiversity, covering about 100,000 km 2 of the Paraguay Basin [36].It behaves as a regulator for the entire La Plata Basin, by slowing the streamflow of the Paraguay River before its junction with the Paraná River [29].Annual totals of precipitation exceed 1300 mm.In December, the precipitation maximum takes place, reaching about 200 mm.In July and August, the precipitation is close to zero [37]. Northern savanna (Cerrado; Figure 1, sub-region 4): The Cerrado is a rather moist savanna region with an average annual precipitation of around 1500 mm.The climate is seasonal with a very strong dry period, spanning from April to September.The temperatures are mild, ranging from 22° to 27°.Large parts of the Cerrado have been transformed into pasture or are used for cash cropping [38]. Eastern evergreen broadleaf forest (Mata Atlantica; Figure 1, sub-region 5): Originally the Atlantic forest was one of the world's largest rainforests, covering approximately 150 million ha.Nowadays, only smaller fragments exist due to intensive deforestation and agricultural use [39].The region has the highest average annual precipitation (above 2000 mm) in the La Plata Basin [30] and receives precipitation during the whole year.The climate is tropical to subtropical.The month of maximum precipitation varies from north (boreal summer) to south (boreal winter).Therefore, the annual cycle of precipitation for the whole region shows three peaks [40]: one at the beginning of spring, one in the middle of summer, and one in winter (with respect to boreal seasons).Due to the height of the region, frost occurs frequently during winter and precipitation may fall in some sub-region as snow (e.g., high places of Santa Cantarina State).

Data
Table 1 gives an overview of the different data products used in this study and highlights differences in spatial resolution, units, and representation.For the analysis, all data sets were processed for the overlapping time span (September 2007-August 2011) with a temporal resolution of one month and a spatial resolution of 1°.If soil moisture or precipitation data sets were compared with GRACE, a harmonization procedure was applied (see Section 4.2).

Satellite Soil Moisture
For this study, two satellite soil moisture data sets with a large overlap in mission operation were selected.The first data set of the active microwave scatterometer ASCAT is provided by EUMETSAT and contains additional information on quality control.We apply the same thresholds as proposed by [17], excluding data points with topographic complexity, inundation and wetland fractions greater than or equal to 50%, and soil moisture error greater than or equal to 10%.Aiming at a better understanding of data quality and a large spatial coverage, these criteria are less stringent as those used in data assimilation studies.
The second data set was developed by Vrije Universiteit Amsterdam and NASA's Goddard Space Flight Center and originates from the passive microwave sensor AMSR-E.As in many other studies [45][46][47], only descending tracks were used herein.For quality control, data points were excluded if any of the flags (which are delivered with the data set) indicating high vegetation in the X-band, high vegetation in the C-band, or frozen ground were set.
ASCAT and AMSR-E, which operate in the C-band at 5.3 and 6.8 GHz, respectively, provide information for the uppermost portion (<5 cm) of the top soil layer.ASCAT data are provided as percentages of dryness, where 0% corresponds to the driest soil surface conditions (lowest backscatter coefficients) and 100% corresponds to wettest soil surface conditions (highest backscatter coefficients).
AMSR-E data give information on the volumetric soil moisture content [m 3 H2O/m 3 soil].

WaterGAP Global Hydrology Model (WGHM)
The global hydrology and water use model WaterGAP has been developed for the assessment of global water resources and water balances of large river basins [48].In addition to water use models for different sectors, WaterGAP comprises the Global Hydrology Model (WGHM), which calculates daily water storage and water flow for all continental areas except Antarctica with a spatial resolution of 0.5°.WGHM accounts for all main continental water storage components including groundwater, soil moisture, snow, canopy storage, and surface waters in lakes, rivers, wetlands, and reservoirs.The water mass of all components adds up to the total continental water storage, which was analyzed by [49] and compared to TWS measurements from GRACE in various studies (e.g., [50][51][52][53]).In this study, the latest version of WGHM in WaterGAP2.2was applied [44].For model calibration against observed river discharge time series, 1319 gauging stations worldwide are used in WGHM.The precipitation data set used as model input in this study was the Monitoring Product Version 4.0 taken from the Global Precipitation Climatology Center (GPCC) [54].Other climate forcing data such as temperature, cloudiness, and number of rainy days per month were taken from operational forecast or analysis data of the European Centre for Medium-Range Weather Forecasts (ECMWF).Soil moisture in WGHM is calculated for one soil layer of which the thickness varies with the land cover depending on the vegetation rooting depth.The maximum available soil water capacity is calculated by multiplying the rooting depth by the total available water capacity in the uppermost meter of the soil [48].

GRACE
In this study, we used GRACE RL05 Level-2 data provided by the GFZ German Research Centre for Geosciences [41].Original GRACE observations of the time-variable gravity include the effects of all mass changes in the atmosphere, on the Earth's surface, and in its interior.During pre-processing, the gravity signal is reduced for mass changes within the atmosphere and the ocean as well as for solid Earth tides [41].The remaining signal predominantly provides information on changes in continental water storage.However, for some regions, post-glacial rebound [55] or other non-hydrological signals may need to be considered [56,57].It is assumed that in the La Plata Basin, the GRACE signal mainly maps changes in soil moisture, surface water, and groundwater storage.
GRACE data are usually provided in spherical harmonic coefficients with a temporal resolution of one month.In this study, spherical harmonic coefficients were processed up to degree and order 70, corresponding to a spatial scale of approximately 300 km.GRACE data show longitudinal stripes due to correlations among certain spherical harmonic coefficients [58].Further errors result from the orbit configuration of the twin satellites [59] and sub-monthly mass changes that are not captured by the background models in the course of the abovementioned reduction process (i.e., aliasing of high-frequency mass variations) [60].To remove correlated errors in the coefficients, the spherical harmonics of the GRACE data were filtered with a least-squares polynomial filter [58].Furthermore, an isotropic Gauss filter [59] with 300-km half-wavelength radius was applied.Subsequently, the spherical harmonic coefficients were converted back into geographical grids of 1° × 1°.
The main remaining errors in the GRACE data arise from spatial leakage and signal loss, as the filtering attenuates both the noise and parts of the signal and enhances leakage of mass variations from the surrounding regions into the study area.In this study, we assume that all data sets are similarly affected by spatial leakage due to identical processing (see Section 4.2).

El Niño/La Niña Index
The Climate Prediction Center of the United States National Oceanic and Atmospheric Administration (NOAA) provides the Oceanic Niño Index (ONI) to identify El Niño and La Niña events in the tropical Pacific [61].The ONI is calculated from the three-month mean sea surface temperature anomaly for the Niño 3.4 region (5°N-5°S, 120°-170°W).If five consecutive overlapping three-month periods exceed a certain threshold, they are classified as El Niño (warm, positive anomaly) or La Niña (cool, negative anomaly) events [62].The threshold is set to 1.5 for strong events, 1.0-1.4 for moderate events, and 0.5-0.9 for weak events [63].

EM-DAT, the International Disaster Database
EM-DAT (the OFDA/CRED International Disaster Database-www.emdat.be-Universitécatholique de Louvain, Brussels, Belgium) is a global database on natural and technological disasters, and it is maintained by the Centre for Research on the Epidemiology of Disasters (CRED).From 1900 onward, more than 17,000 disasters that meet at least one of the following criteria have been registered [64]: -ten or more people killed; -one hundred or more people affected; -declaration of a state of emergency; and -call for international assistance.The main EM-DAT disaster types considered in this study were hydrometeorological in nature and classified as storms, floods, extreme temperature, and droughts.We neglected wet mass movements, landslides, avalanches, and land subsidence as their impacts are limited to local scales.We also neglected wildfires.A further criterion for the selection of events was the size of the affected area (events were not considered, e.g., if only one city was affected).Figure 2 shows a rough spatial distribution of the selected disasters within the La Plata Basin.Most of the disasters were registered in densely populated regions such as the east coast of South Brazil and the surroundings of the metropolitan area Gran Asunción in Paraguay [65].In sparsely populated regions such as the Chaco, the Pantanal, and the Campos, only a few events were listed.This is due to the fact that extreme weather events are only classified as natural disasters within EM-DAT if they had a high impact on society.

Scaling
ASCAT and AMSR-E only map the near-surface soil moisture, whereas WGHM provides the soil moisture content for the entire root zone.All three data sets cannot be compared directly in terms of magnitude.Additionally, ASCAT data are provided as percentages of dryness.To harmonize the units and approximate the root zone soil moisture for all data sets, we scaled ASCAT and AMSR-E data with respect to the soil moisture content provided by WGHM.A pixel-wise scaling algorithm was applied because the depth of the soil column may vary among the pixels.The time series of each pixel was scaled with respect to the minimum and maximum soil moisture values of the same pixel, as provided in the reference soil moisture data (WGHM herein).
The simple scaling has the advantage that it does not influence the correlation between data sets, as shown by [47].For example, the correlation between scaled ASCAT soil moisture and GRACE TWS yields the same correlation coefficients as the correlation between unscaled ASCAT soil moisture and GRACE TWS.However, the simple scaling has the disadvantage that it does not account for the heterogeneous vertical distribution of soil moisture.Also the maximum and minimum values of root zone soil moisture, taken from the WGHM for the scaling, may contain errors.Due to this limitation, soil moisture values are not discussed in absolute terms in this study.
Alternative algorithms exist, which derive root zone soil moisture from satellite surface soil moisture.For example, several studies suggest the use of an exponential filter to calculate the so called soil water index (SWI) [66,67].The algorithm comprises a parameter (called time characteristic length), which depends on the depth of the root zone and the soil properties.Such an algorithm is not used here, as it requires reference data (which is optimally in situ data) for calibration.Furthermore, in the case of filtered data, correlation results are not any more independent from the reference data.For example, the correlation between the SWI from ASCAT and GRACE TWS does not yield the same correlation coefficients as the correlation between unscaled ASCAT soil moisture and GRACE TWS.

Harmonization
The data sets from satellite gravimetry, soil moisture remote sensing, and hydrological models differ in temporal and spatial resolution, representation (grid format vs. spherical harmonics), and processing (e.g., filtering).To avoid these differences that strongly impact the analysis, all data sets were harmonized in a first step, as thoroughly described in [17].Due to the particular processing of GRACE data (see Section 3.3), all soil moisture data sets were converted into spherical harmonics of degree and order 70 and an isotropic Gauss filter with a 300-km half-wavelength radius was applied.Subsequently, the spherical harmonic coefficients were converted back into geographical grids of 1° × 1°.Regions affected by snow, ice, and post-glacial rebound or very low variation in soil moisture must be masked out when comparing soil moisture data with GRACE [17].In the case of the La Plata Basin, no masking is necessary.
The main impact of the harmonization procedure is the spatial smoothing of the data, which causes damping of the signal, loss of detail, and reduction of noise.This leads to an increase in the correlation coefficients between the data sets.Therefore, data are harmonized in this study only if a comparison with GRACE data is sought (e.g., comparing GRACE with satellite soil moisture products).Otherwise (e.g., for the comparison of different remote sensing soil moisture products), only the spatial resolution is harmonized by computing the simple average of all data points falling into a 1° × 1° grid cell.In any case, the temporal resolution of all data sets is aggregated monthly, which is similar to the standard GRACE resolution.

Correlation and Principal Component Analysis
Two well-known methods were used for the comparison of data sets.The first one is the Pearson product-moment correlation coefficient [68], which is defined as the covariance of two variables divided by the product of their standard deviations.As in other studies, which focus on the correlation of various soil moisture data sets [69][70][71], Fisher's z-transformation is used to calculate confidence limits for the Pearson product-moment coefficients.
The second method is the principal component analysis (PCA).With the help of PCA, a spatio-temporal signal (e.g., a time series of maps) is split into a series of empirical orthogonal functions (EOFs) that describe the spatial base pattern of the signal.The corresponding coefficients are called principle components (PCs) and describe the temporal behavior of the signal.The first mode (first pair of EOF and PC in the series) describes the most dominant part of the signal.Higher modes have decreasing importance [72].PCA has already been used in arid environments to compare satellite soil moisture products and GRACE [73].

Sample Cross-Covariance Function
To identify possible time shifts in the signals of different data pairs, the time lag maximizing the cross correlation between two signals was calculated.Therefore, the sample cross-covariance function described by [74] was used.It provides an estimate of the covariance (c) between two time series xt ([x1,x2,… xT]) and yt ([y1,y2,… yT]) for a defined time lag k: where ̅ is the mean of the first time series and is the mean of the second time series.The formula for estimating the cross correlation (r) is: where sx is the variance of the time series xt ( 0 ) and sy is the variance of the time series yt ( 0 ).In this study, xt is the signal of one data set, yt is the signal of another data set, T = 48 months, and n = 6 months.

Comparison of Soil Moisture Data Sets
Before linking the signatures of soil moisture to extreme weather conditions, the different global soil moisture data sets of the hydrological model WGHM and the two satellite products ASCAT and AMSR-E were compared.Figures 3 and 4 show that ASCAT and WGHM agree well in most parts of the La Plata Basin with correlation coefficients above 0.7 for the monthly time series (90% confidence interval ranges from 0.6 to 1).Exclusively in the evergreen broadleaf forest in the eastern part of the study area, the soil moisture signal from WGHM lags several months behind the ASCAT signal (lag −3-−6) and therefore has to be shifted forward in time (Figure 4, row 2, column 3).However, even if the time shift is implemented, the correlation does not improve much (Figure 4, row 3, column 3).For soil moisture from AMSR-E and WGHM, the correlation coefficients are lower than those for ASCAT and WGHM.Especially in the western part of the basin, they are negative or close to zero since the seasonal signals are out of phase (Figure 3, column 1, rows 3-4 and Figure 4, row 1, column 2).The poor correlation values agree with the results of [51], which reports large errors for the AMSR-E signal in the La Plata Basin.In the eastern evergreen broadleaf forest, neither ASCAT nor AMSR-E agree well with soil moisture from WGHM.Although WGHM model deficits cannot be excluded, this indicates that in this region, the microwave signal received by ASCAT and AMSR-E is severely degraded by the topographic complexity [69] and/or the very high foliage density, which shield the soil moisture signal [75].
Figure 3 shows that in the southeastern part of the basin (rows 6 and 7), the seasonal pattern of precipitation and soil moisture for the relatively short study period considered here is dominated by the anomaly of the El Niño event in the boreal winter of 2009/2010.The absence of a marked seasonality in this region in the absence of extreme events is also described by [30] for streamflow data.The strong anomaly in the Campos region (column 2, row 6) is consistent with a severe flood that Uruguay suffered from in November 2009, affecting 22,000 people.In the eastern evergreen broadleaf forest (column 2, row 5) the anomaly has approximately the same amplitude as the seasonal signal and shows high fluctuation.This region, including the state of São Paulo, Brazil, is densely populated [65].Both facts explain why the region was struck by five major disasters within only four years (compare Figure 2).

Hydrological Extremes
The comparison of the anomalies of different hydrological parameters in the La Plata Basin (Figure 5) shows the close connection between soil moisture and precipitation.Furthermore, the anomalies are in agreement with time variations of water levels at different locations in the La Plata Basin, as presented by [23].The comparatively smooth GRACE time series is an expression of the integrative nature of TWS that, to some extent, balances out short-time variations in individual storage components.Furthermore, anomalies of GRACE TWS tend to be shifted by a few months with respect to soil moisture and precipitation, as it takes time until precipitation extremes accumulate within surface and ground water storage.A drastic decrease in soil moisture can therefore be an indicator for an upcoming deficit or surplus in total water storage.One example is the La Plata drought in 2009, wherein soil moisture was lowest between December 2008 and January 2009, and total water storage reached its minimum in May 2009.The minimum soil moisture in December 2008 agrees with the fact that the drought was declared and registered in January 2009 as a natural disaster within the International Disaster Database EM-DAT.Furthermore, the local minima for soil moisture in May 2009 and the decrease in total water storage until May 2009 illustrate that the drought lasted for several months until approximately November, 2009.The clear signal of the La Plata drought of 2008/2009 in GRACE TWS was also observed by [20].The flood period in the boreal winter of 2009/2010 (as supported by [22]) and the flood period in the beginning of 2011 are also visible in the GRACE data.Apart from these events, it was not possible to link anomalies of GRACE TWS directly to the listed natural disasters.This is due to the coarse spatial resolution of the GRACE sensor and its low sensitivity to small mass changes.Also short-time events such as storms might not affect all water storage components contained in TWS.
In contrast, almost all droughts and floods could be identified as single events or as accumulated events (e.g., winter 2009/2010 and March and April 2011) in at least one of the soil moisture data sets from ASCAT, AMSR-E, or WGHM.The only exception is the flood in Brazil in November 2010, which affected the northeastern fringe area of the basin, including parts of Minas Gerais, São Paulo, and Rio de Janeiro.Furthermore, some meteorological events (storms and cyclones) could be linked to anomalies in the soil moisture data, e.g., the storm in Paraguay in August 2008.During this month, precipitation data from GPCC and the two soil moisture data sets from AMSR-E and ASCAT show local maxima.
Figure 6 emphasizes that specific extreme events can be identified more clearly when focusing on smaller regions.For example, the flood in Bolivia from December 2007 to February 2008 is clearly visible in all soil moisture signatures when averaging over the western sub-region in the Chaco only.At the same time, it became clear that the minima of the AMSR-E signal in June 2010, as seen in the basin average in Figure 5, is not associated with the drought in Bolivia.Thus, a differentiation into hydrological extreme events of basin-wide character or of regional extent is only possible with the given data sets.
Figure 5 does not show the expected below average precipitation in the La Plata Basin during the moderate La Niña period of 2007/2008.However, atypical precipitation patterns during this La Niña period were also observed in southeastern Australia [76].From the summer of 2008 onward, there is a strong relationship between variations in soil moisture and the ONI (Figure 5).The first phase of the La Plata drought in 2009 fell into the weak La Niña period during the boreal winter of 2008/2009, which is also supported by [20].The dominant spatial patterns of the various data sets and the major differences among them were further analyzed by PCA.The first, second, and third modes explain approximately 60%, 25%, and 5% of the TWS signal from GRACE and the soil moisture signals from ASCAT and WGHM, respectively.For AMSR-E, the first three modes account for 50%, 25%, and 10% of the signal.The first mode mainly represents the seasonal part of the signal (Figure 7, row 1).In the PCs, no time shift between the seasonal cycles of the different soil moisture products is visible.For GRACE, the phase of the seasonal signal is approximately two months later.The EOFs of the first mode show a characteristic opposing pattern in the northern part of the La Plata Basin versus in the southeastern part (Figure 8, row 1).The inverse annual signal for AMSR-E visible in Figure 3 (column 1, rows 3-4) in the western and northern parts of the basin is reflected by the inverse signs of the EOFs of AMSR-E compared with the EOFs of the other soil moisture products in this region.The nonconformity between AMSR-E and other soil moisture data sets is also reflected in the PCs and EOFs of the second mode.The PCs of the second and third modes show secondary minima for soil moisture from ASCAT and WGHM during the La Plata drought in 2009 (Figure 7, rows 2-3).This indicates that the southern, central (Figure 8, row 2, columns 2-3), and eastern parts (Figure 8, row 3, columns 2-3) of the basin were mostly affected.The described spatial pattern is in agreement with precipitation maps for the La Plata Basin in winter and spring 2009 from the Argentinean National Meteorological Service (SMN), as presented by [22].
Similarly, the flood period associated with El Niño in the boreal winter of 2009/2010 can be identified in the PCs of the third and second modes of the soil moisture data.As the anomaly in the third mode (Figure 7, row 3) appears before the anomaly in the second mode (Figure 7, row 2), the eastern part of the basin (Figure 8, row 3, columns 2-3) followed by the southeastern part (Figure 8, row 2, columns 2-3) were affected.This sequence of spatial patterns is also supported by the precipitation anomalies of SMN [22].Furthermore they are reflected in the signatures of precipitation and soil moisture in Figure 3 (column 2, rows 5-7).
For GRACE, the two periods of extremes are also clearly visible in the PCs of the second and third modes.For the La Plata drought in 2009, the second mode is shifted by a few months and the third mode by several months with respect to the signatures of soil moisture.The EOF of the second mode (Figure 8, row 2, column 4) shows that the northeastern and central parts of the basin were first affected.Then, based on the EOF of the third mode (Figure 8, row 3, column 4), the southern part was affected.The temporal development of the drought provided by the PCs is similar to the results of [17]; the minimum values of water storage in the entire basin occurred during the first half of 2009 followed by strong but regionally confined low water storage in the southern part in July and August.
For the first half of 2009, the soil moisture data show that the largest anomalies were in the southern, central and eastern parts of the basin.However, within the International Disaster Database EM-DAT, the largest impact of the drought was registered in the central and northwestern parts of the basin, affecting more than 227,000 people in Bolivia and Paraguay.It is possible that the International Disaster Database EM-DAT lacks information on the impact of the drought; for example, the database does not contain any information on the drought situation in Uruguay.A situation report of the United Nations Development Program in Uruguay (dated at the beginning of the drought in January 2009) clarifies that Uruguay was affected by the drought (e.g., through food shortages and lack of water for livestock) [77].However, the report also states that there was a lack of integral information on the impact of the drought and the economic consequences.This example shows that databases like EM-DAT might be incomplete due to the difficulty in collecting and receiving quantitative data on an international level.Moreover, it should be noted that the disaster classification given by EM-DAT explicitly considers the vulnerability of the region, which is defined by demographic factors such as population density (which is relatively low in Uruguay [65]), environmental factors (e.g., state of resource degradation and depletion), social factors (e.g., traditional knowledge systems), and economic factors (e.g., economic status of individuals, communities, and nations) [78].In contrast, the soil moisture and TWS products only reflect the physical factors in terms of a subset of hydrometeorological hazards.

Soil Moisture and GRACE TWS
Figure 9 highlights two regions where the dynamics of soil moisture from ASCAT and TWS from GRACE are interrelated.In both regions, soil moisture is not the dominant component of TWS, but it is closely linked to the most dominant water storage components (surface water and groundwater herein).The first region with high agreement is the northern wetland Pantanal (correlation coefficients are around 0.8, 90% confidence interval ranges from 0.7 to 0.9).In the Pantanal, the shallow groundwater table may rise and saturate the surface soil or fall and leave the surface soil dry [79].Therefore, near-surface soil moisture dynamics monitored by ASCAT can be expected to be a good indicator of TWS variations.The second region with a good correlation between ASCAT and GRACE (correlation coefficients range from 0.6 to 0.7; 90% confidence interval is 0.5 to 0.8) ranges from the junction of the Paraguay and Paraná Rivers, which is surrounded by wetlands (Figure 1), to the southeastern coast of the basin.Again, in the wetland region, the dynamics of groundwater and soil moisture tend to be closely linked.Approaching the coast, the connection to surface water is stronger.
The correlation of ASCAT time series to precipitation is high for most parts of the study area (correlation coefficients are above 0.7 as shown in Figure 9 and the 90% confidence interval ranges from 0.6 to 0.9).This is reasonable as near-surface soil moisture variations seen by ASCAT are directly triggered by the precipitation dynamics.Nevertheless, the correlation tends to decrease toward the southeastern downstream parts of the basin, whereas the correlation of ASCAT with WGHM simulation-based surface water storage simultaneously increases.In this region the groundwater recharge is high (>100 mm) [23].This implies that the relevance of local rainfall-induced soil moisture variations decreases in favor of soil moisture variations caused by allochthonous processes at the regional scale, i.e., surface water inflow from upstream areas of the river basin.
Figure 9. Correlation after data harmonization for soil moisture from ASCAT with respect to precipitation from GPCC (column 1), total water storage from GRACE (column 2), groundwater from WGHM (column 3), and surface water storage from WGHM (column 4).

Conclusions
In this paper, we investigated how total water storage from GRACE and soil moisture from two remote sensing satellites sensors (ASCAT and AMSR-E) and one hydrological model (WGHM) map extreme weather events in the La Plata Basin of South America.Information from the International Disaster Database EM-DAT was included in the analysis to place special emphasis on extreme hydrometeorological events that were highly destructive.Furthermore, the link between variations in soil moisture and TWS was analyzed.Therefore, the distinct representations and formats of all data sets were harmonized by conversion in spherical harmonics and filtering whenever any data set was compared with GRACE.For comparison, we correlated different data pairs, analyzed the time shift between seasonal cycles, and performed PCA to split the signal in spatial base functions and temporal components.
First, the different soil moisture data sets were compared.The results indicate that soil moisture products from ASCAT and WGHM are strongly correlated (correlation coefficients above 0.7) in the majority of the sub-regions of the La Plata Basin.For AMSR-E, the correlation coefficients take values above 0.6.However, in the eastern deciduous broadleaf forest, the western evergreen broadleaf forest, and the wetland Pantanal, signatures of AMSR-E differ strongly from those of both WGHM and ASCAT.In the eastern forested region, neither ASCAT nor AMSR-E correlate well with soil moisture from WGHM, possibly due to signal disturbance by the topographic complexity and dense vegetation in this region.
Second, extreme meteorological and hydrological conditions, specifically those that have been classified as natural disasters within the International Disaster Database EM-DAT, were brought into focus.Thus, two independent data sources were linked.On the one hand, signatures of hydrological variables were analyzed to gather information on hydrometeorological extremes.On the other hand, natural disasters that affected the La Plata Basin during the analyzed time period were identified to include information on the number of people affected or killed by extreme events.The analysis of both data sets showed that most of the strong variations in soil moisture match temporally with natural disasters (droughts, floods, or storms).For two hydrologic extreme periods that led to several natural disasters (the La Plata drought in 2009 and the flood period in the winter of 2009/2010), it could be observed that soil moisture first changes drastically followed by TWS from GRACE with a delay of a few months.In these two cases, soil moisture serves as an indicator for the later deficit or surplus in TWS.
Out of four El Niño and La Niña events during the analyzed time period, three events match with extreme soil moisture and TWS conditions.These events include the weak La Niña event in 2008/2009, the moderate El Niño event in 2009/2010, and the strong La Niña event in 2010/2011.The most meteorologically and hydrologically affected areas were the southern, central, and eastern parts of the La Plata Basin.However, based on the International Disaster Database, the main impacts were registered in the central and northwestern parts of the basin.This reveals that several limitations have to be considered when linking extremes in soil moisture or total water storage to natural disasters.On the one hand, not all extreme hydrologic events have large impacts on society.This might be due to the fact that the destructiveness of a disaster not only depends on the magnitude of the hydrometeorological extreme indicated by GRACE, ASCAT or AMSRE but also on the vulnerability of the region considered in EM-DAT, which largely depends on other factors such as the preparedness of the population and the population density.On the other hand, the collection of quantitative information on the impacts of disasters is challenging on an international level.Therefore, limited data availability in the EM-DAT database for the southern parts of the study area cannot be excluded.
Analyzing the hydrological signals as basin averages and with a temporal resolution of only one month makes the identification of regional or short-term events difficult.Nevertheless, it could be shown that the remote sensing products used here allow for differentiating hydrological extreme events of basinwide character or of local to regional extents.
Finally, the discussion focused on possible interdependencies between the dynamics in soil moisture and the dynamics in TWS.A close link between soil moisture and total water storage was identified in regions where groundwater or inflow of surface water from upstream areas influences the dynamics of soil moisture.

Figure 1 .
Figure 1.Land cover map of the La Plata Basin in South America (MODIS (Moderate-resolution Imaging Spectroradiometer) Land Cover Type Climate Modeling Grid product (MCD12C1) with IGBP (International Geosphere-Biosphere Programme) land cover classification, https://lpdaac.usgs.gov)with its three main contributing rivers (Paraguay, Paraná, and Uruguay).The sub-regions selected for this study are shown with contours and include: 1, southern cropland (Pampas); 2, western deciduous broadleaf forest (Chaco); 3, wetland Pantanal; 4, northern savanna (Cerrado); 5, eastern evergreen broadleaf forest (Mata Atlantica); and 6, southern grassland (Campos).The black continuous line shows the outline of the La Plata Basin and the black dashed lines indicate the boundaries of the countries.

Figure 2 .
Figure 2. Location of natural disasters in the La Plata Basin from October 2007 to July 2011 based on the International Disaster Database EM-DAT.

Figure 3 .
Figure 3. Monthly mean (column 1) and monthly anomaly (column 2) for soil moisture from ASCAT, AMSR-E and WGHM (ASCAT and AMSR-E are scaled with respect to WGHM), and for precipitation from GPCC for the entire La Plata Basin (row 1) and the different sub-regions delineated in Figure1[including the northern savanna (row 2), the wetland Pantanal (row 3), the western deciduous forest (row 4), the eastern evergreen broadleaf forest (row 5), the southern grassland (row 6), and the southern cropland (row 7)].

Figure 4 .
Figure 4. Correlation (row 1), time shift in months (row 2; negative values indicate that the seasonality of the first soil moisture product lags behind the second product), and correlation after the implementation of time shift (row 3) for soil moisture from AMSRE and ASCAT (column 1), AMSRE and WGHM (column 2), and ASCAT and WGHM (column 3).

Figure 5 .
Figure 5. Disasters as registered by the International Disaster Database EM-DAT for the La Plata Basin (top), number of people affected and killed for each disaster (second from top), monthly anomalies as basin averages after data harmonization for TWS from GRACE, precipitation from GPCC (third from top), and soil moisture from ASCAT, AMSR-E, and WGHM (ASCAT and AMSR-E are scaled with respect to WGHM; fourth from top), and El Niño and La Niña anomalies as described by NOAA through the Oceanic Niño Index (ONI; bottom).
Furthermore, the moderate El Niño period in the boreal winter of 2009/2008 coincides with the very high soil moisture values in December 2009 and several floods affecting Argentina, Brazil, Bolivia, and Uruguay.The anomaly is also visible in TWS from GRACE and initiated by high precipitation during November and December 2009.The strong La Niña period during the boreal winter of 2010/2011 is shown as local minima during November 2010 for all parameters (except for AMSR-E).

Figure 6 .
Figure 6.Disasters as registered by the International Disaster Database EM-DAT for the deciduous broadleaf forest (Chaco) in the western part of the La Plata Basin (see Figure 1) (top) and monthly anomalies of soil moisture from ASCAT, AMSR-E, and WGHM (ASCAT and AMSR-E are scaled with respect to WGHM; bottom).

Figure 7 .
Figure 7. Principal components for the first, second, and third modes for soil moisture from AMSR-E, ASCAT, WGHM (AMSR-E and ASCAT data are scaled with respect to WGHM), and GRACE after data harmonization.For the second and third modes, the shaded areas indicate the La Plata drought in 2009 (yellow) and the flood period associated with El Niño in the boreal winter of 2009/2010 (blue).

Figure 8 .
Figure 8. Empirical orthonormal functions for the first, second, and third modes for soil moisture from AMSR-E (column 1), ASCAT (column 2), and WGHM (column 3) (AMSR-E and ASCAT data are scaled with respect to WGHM) and total water storage from GRACE (column 4) after data harmonization.

Table 1 .
Specifications of the main data sets of soil moisture and total water storage.