The 2018 Long Rainy Season in Kenya: Hydrological Changes and Correlated Land Subsidence

The wettest 2018 long rainy season (March to May (MAM)) resulted in daily intensive rainfall events in East Africa that have seriously affected the environment and economy in many countries. Land subsidence is one of the environmental disasters that has occurred due to the long rainy season in Kenya for many years. However, it has received limited scientific attention. In this paper, we incorporate hydrological (soil moisture active–passive (SMAP) and loading models) and geodetic data (global positioning system (GPS) and interferometric synthetic aperture radar (InSAR)) to study hydrological changes and their associated subsidence potential in Kenya. Results show that widespread subsidence of more than 20 mm was associated with the MAM season in Kenya during 2018, based on SBAS InSAR measurements. The high values of land subsidence were well correlated with the areas of intense flooding during the MAM season. The widespread subsidence during the wet season has implications for the stability of the earth’s surface during the season rather than creating the possibility of potential stresses along active faults. These stresses may trigger seismicity that is expected to pose risks to urban features. The results of the current study can help governmental authorities to adopt proper urban planning that avoids or minimizes the risks of land subsidence in the areas of sinkholes.


Introduction
Following the devastating drought that struck the Horn of Africa for two years (2016-2017), the wettest 2018 long rainy season (March to May (MAM)) resulted in episodes of exceptional rainfall every day [1]. The rainfall events triggered extensive flooding with a notably high impact on society, including a significant number of casualties and the displacement of around one million people. This impact extended to environmental and economic disasters, including crop failures, livestock deaths, and water and shallow aquifer contamination across Kenya, Somalia, Ethiopia, Uganda, and South Sudan [2][3][4]. Land subsidence is one of the environmental disasters that has occurred as a result of the rainy seasons in Kenya for a long time [5]. However, few studies have been conducted into this phenomenon.
The hydroclimatic cycle (event-to-event) dominates the terrestrial hydrological regimes in East Africa, where a prolonged dry spell may be interrupted by exceptionally wet seasons [6]. East Africa experiences two rainy seasons caused by different climate circulation mechanisms: a short season, extending from October to December, and the long season that extends from March to May (MAM). The short season is typically associated with a strong Indian Ocean dipole (IOD) and the El Niño duration of the MAM season using precipitation (P) and soil moisture (SM) anomalies from across the East African countries (Kenya, Uganda, Ethiopia, Somalia, and South Sudan). Second, the vertical land deformations and flood extent are mapped using the Sentinel-1 satellite radar system. Third, the estimated magnitude of land deformations from GPS and InSAR are compared with the model calculation of load deformations.

Materials and Methods
Different hydrological, geodetic, and loading models' data are incorporated into this research as described below.

Hydrological Data (CPC, SMAP, VIC, NOAH-3.3)
The daily precipitation data of the climate prediction center (CPC) project was used to map the intensity of precipitation during the MAM season. CPC data is available globally at a spatial resolution of 0.5 • * 0.5 • based on global gauge stations [25]. Similarly, soil moisture (SM) saturation was estimated from the soil moisture active-passive (SMAP) satellites. SMAP satellites provide daily data regarding land surface soil moisture with a 2-3 days' repeat cycle. The images of 9 km spatial resolution is included in this research, which represents a daily global composite product of the surface soil moisture [26,27]. Also, two global hydrological models were used to investigate the state of the hydrological fluxes and SM in Kenya before, during, and after the MAM. These models are the variable infiltration capacity (VIC) model and NOAH-3.3 model [28,29]. The VIC model is one of four land surface models that simulate the states of water and energy budget components derived from the global land data assimilation system of version 1. The VIC model simulates various hydrological fluxes, including P, evapotranspiration (ET), and runoff (R) along with the storage of SM, canopy water, and snow water. The version implemented here has a temporal resolution of a month and a spatial resolution of 1.0 • * 1.0 • . Similarly, NOAH-3.3 simulates the same hydrological variables with different parameters. NOAH-3.3 represents the GLDAS-2.1 version with a spatial resolution of 0.25 • * 0.25 • at monthly time intervals [28][29][30]. The differences between the two versions of GLDAS rely on correcting the biases in forcing data of the agricultural meteorological modeling system (AGMET) radiation flux field and using the global precipitation climatology project (GPCP) and tropical rainfall measuring mission dataset instead of the CPC merged analysis of precipitation (CMAP) dataset [31]. To show the intensity of the precipitation, the data has been fitted to the regional empirical cumulative density functions (ECDF) [32]. ECDF is probability distribution functions of the precipitation rate and soil moisture saturation during the MAM season. It represents the cumulative distribution of rainfall and soil saturation in East Africa at monthly scale relative to the total populations (January-June 2018). It ranges from 0 to 1, where the values close to zero indicate small or no areas under flooding with low soil saturation and low rainfall. The maximum values indicate the highest saturation of the soil and high rainfall, therefore, the highest probability of flooding. The averaged time series of hydrological fluxes and SM changes across Kenya through the MAM event were obtained based on VIC and NOAH-3.3 models.

Geodetic Data
We used two independent geodetic data sources (GPS and SAR) to map and estimate the deformations during the 2018 flood events associated with the MAM season. The GPS data were represented by the available (DKUT) station located at longitude 36.961 • E and latitude 0.398 • S and at an altitude of 1775 m above sea level. This data was acquired from the Nevada geodetic lab [33]. The DKUT time series were organized in a processing chain as part of the global solution data, including more than 8000 stations [33]. The daily RINEX GPS files were processed using the precise point positioning method [34] implemented by the GPS inferred positioning system (GIPSY) OASIS II software. A number of conventions and models are used for correction and to obtain the final results, the estimation of white noises, tidal ocean loading, and pole tidal loading, including the satellite orbital information and calibrating the antenna phases [35,36]. The final result of the coordinates is compared to the global reference frame from the International GNSS Service 2008 (IGS08). To suppress the common-mode errors (CME) [37], which usually result from previous corrections, the singular spectrum analysis (SSA) method was used for the final results. SSA decomposes the time series into interpretable components with no prior assumptions [38]. The CME was identified as a highly varied component. The trend component was reconstructed using eigentriples one and four, and for the season component eigentriples three, five, and six were used. The rest of the 10 eigentriples represented the CME component. The analysis was applied only to the vertical component. By applying the SSA, the variance of the component was reduced from 9.6 mm in the raw data to 3.8 mm in the reconstructed components. SSA was used to model the seasonal signal from the GPS time series [38,39].
The InSAR data were represented by 74 images (Figure 2 and Table S1 in the Supplementary Information Section) acquired by a Sentinel-1 satellite from August 2016 to June 2019 with an ascending orbit of path 57 and frame 1175 (Figure 1), Level-1, single-look complex (SLC). The InSAR radar data with a temporal baseline of 100 days and spatial baseline of up to 100 m were coregistered at the pixel level using the orbital information and the 30 m digital elevation model (DEM) acquired by the shuttle radar topography mission (SRTM) [40]. The coregistered images were then used to generate a high number of interferograms, 518 pairs, that allowed monitoring of the short temporal deformations ( Figure 2). The images were aligned to the super master image captured on 26 February 2018. The interferometric phase was then obtained by (1) removing the topographic phase with the DEM of 30 m spatial resolution, (2) suppressing the noise of the phase using an adaptive filter of wavelength 300, and (3) unwrapping the remaining phase based on the statistical-cost, network-flow algorithm [41,42]. The unwrapped phases were combined to generate a coherence rate of land displacement at different temporal scales using the small-baseline subset method (SBAS) [43,44]. To estimate the contributions of the topographic and ionospheric delay, as well as to remove them, the iterative common scene method was used [45]. Finally, the rate of deformation at the acquisition dates and the mean velocity were estimated using the least squares method [41].

Materials and Methods
Different hydrological, geodetic, and loading models' data are incorporated into this research as described below. The daily precipitation data of the climate prediction center (CPC) project was used to map the movement of the mass center of a given geographic location compared to the global mass center based on water storage from the land surface discharge model (LSDM) [49]. The total mass represents the masses of the snow, soil moisture, surface water, rivers, lakes, and wetlands. The data were resampled at a 0.125° * 0.125° degree, and available from 1976 to present. Here we used the weighted loading area from the model for central Kenya in one year extending from December 2017 to December 2018. The data represent the load deformation resulting from changes in flood plains, soil moisture, and groundwater during one year of dry and wet seasons in central Kenya.  Table S1.  Table S1.
To map the spatial extent of the MAM flooding in Kenya, the same Sentinel-1 scenes were first processed for the radiometric calibration, in which the backscattering coefficients were produced from the pixel values [46]. Speckle filtering using a boxcar filter of window size 3 * 3 was applied to the absolute phase of the scattered coefficients to minimize the noise levels in the images [47]. To separate the flooded areas from the non-flooded areas, the filtered, backscattered histogram was analyzed.
The lower values represented the flooded areas, and the higher values were the pixels showing no impact from the flood. Different threshold values were tested to produce an effective value, using the values of known water bodies. The value that produced the 90 th percentile was adopted to discriminate the flooded and the non-flooded areas. The newly created images were geometrically corrected using the range-doppler terrain correction [48].

Loading Model (HYDL)
To investigate the flood-loading deformation associated with the MAM season in Kenya, we used the hydrological loading model. The daily global vertical deformation model was developed based on weather models that included the runoff and the total water storage available at a 0.5 • * 0.5 • regular grid at daily temporal resolution. The input data for the model were outputs from the European Center for Medium-Range Weather Forecasts (ECMWF) [49]. The model represented the movement of the mass center of a given geographic location compared to the global mass center based on water storage from the land surface discharge model (LSDM) [49]. The total mass represents the masses of the snow, soil moisture, surface water, rivers, lakes, and wetlands. The data were resampled at a 0.125 • * 0.125 • Remote Sens. 2020, 12, 1390 6 of 16 degree, and available from 1976 to present. Here we used the weighted loading area from the model for central Kenya in one year extending from December 2017 to December 2018. The data represent the load deformation resulting from changes in flood plains, soil moisture, and groundwater during one year of dry and wet seasons in central Kenya.

Hydrological Responses to the Long Rainy Season MAM 2018
Using the ECDF method, the distribution of the precipitation and soil moisture changes were mapped across the East African countries from January to June 2018. During January to the end of February, the region experienced low rainfall and dry soil (represented by bright colors in Figure 4) in Kenya, Somalia, Ethiopia, South Sudan, Uganda, Rwanda, and Burundi. Only in Tanzania was the precipitation higher, and the soil was saturated during January (represented by the dark color in Figure 4). Precipitation was below average for the two consecutive rainy seasons of 2016 and 2017, and the prolonged drought affected the whole Horn of Africa, resulting in hunger and humanitarian crises [50]. The 2016-2017 prolonged drought was associated with the La Niña event [51]. Following the dry years, precipitation increased to be higher than normal in February in Rwanda, Burundi, and Ethiopia. However, it was not enough to fully saturate the soil. Starting from March and continuing to the end of May, most of the region experienced a highly wet season where ECDF ranged from 0.75 to 1. Similarly, the saturation of the soil ranged from 0.75 to 1, except in South Sudan where the saturation did not reach high values during April. The rainfall and saturation of the soil reached the optimum level during April, when the wet soil could not absorb any more water, producing overland flow and flooding in Kenya, Somalia, Ethiopia, Eretria, Djibouti, Uganda, Rwanda, and Burundi. By the end of May, the situation was back to normal (saturation near zero) across Kenya and Somalia. However, in South Sudan and Ethiopia, the precipitation was still high. The MAM season of 2018 is considered one of the wettest seasons on record in Kenya. This is generally attributed to the combined effect of the active MJO events in phases 2-4 and a tropical cyclone event in the southern Indian Ocean [1].
The hydrological variations before and after the long rainy season in the water budget in the study area including the input (P) and outputs (evapotranspiration, surface runoff, and soil moisture) based on VIC and NOAH-3.3 models are presented in Figure 3. When comparing the precipitation, it can be seen that the peak of the wet season in 2016, as estimated by the NOAH-3.3 model, is higher than that of 2018. The MAM season of 2016 represented the wet season that interrupted the 2015 drought. However, this trend is significantly overturned using the estimation of P based on the VIC model (Figure 3), in which the peak of P in the 2018 MAM season is much higher than that of 2016. Both the NOAH-3.3 and VIC models show a dry season in 2017. The maximum peak of P through the 2018 MAM season occurred during April:~243 mm and~205 mm as estimated from NOAH-3.3 and VIC models, respectively. For ET, the peaks of the two rainy seasons in 2016 and 2017 were smaller than that of the 2018 season, which is consistent with the region's dominant hydroclimatic changes. However, the estimated peaks of the ET generated by the NOAH-3.3 and VIC models are highly comparable in the 2016 season: 70-74 mm/month (Figure 3). Comparable values are also observed for the 2018 season in the two models (i.e., 80-105 mm/month) (Figure 3). For R, which includes the total runoff (surface flow and subsurface flow), both models are consistent with the drought and flood changes. However, the NOAH-3.3 runoff is higher than that shown by the VIC model. The peaks of the R during the MAM 2018 were 46 mm/month for the NOAH-3.3 and 25 mm/month for the VIC model. The subsurface flow constitutes more than 50% of the total runoff estimated by the NOAH-3.3. NOAH-3.3 may overestimate R and underestimate the ET in the region, whereas the VIC model may underestimate the R or overestimate the ET. In general, according to the models, the disaggregation of P into ET and R in land surface models varies with climate regions, and the estimates are sensitive to the ET and R parameters [52].
relatively high values of 55 mm (Figure 4). The drastic deficit in the root zone soil moisture was estimated using the NOAH-3.3 model only. The hydrological fluxes are aligned with the recorded flooding that occurred through the 2018 MAM season. Many consequences were recorded in association with the flooding of this season. One of these consequences was the land deformation that we monitored and measured in the study area located in central Kenya.

Land Deformation Related to MAM Season
The land motions associated with the exceptional rainy MAM season demonstrate varying magnitudes of subsidence, shown by both GPS and InSAR measurements (Figures 5 and 6). The precipitation rates between 2017 and 2018 were plotted against the estimates of the load deformations, as determined from the hydrological loading model (HYDL) and the nonlinear trend of the GPS station, filtered by the singular spectral analysis (SSA) method. For the first three months (December 2017 to February 2018), the precipitation level showed an average of 23 mm/month. Throughout this temporal interval, the surface displacement was marginal, with an average of 1.0 mm uplift as determined from both the HYDL and the GPS station. Starting from March 2018, the precipitation increased markedly, with a magnitude of ~180 mm, and it continued to increase, reaching its highest magnitude of ~250 mm in April. This significant increase in rainfall greatly developed terrestrial water storage, resulting in mass load-induced land subsidence. Both the HYDL model and the GPS measurements reflected sharp land subsidence of 5 mm through March and April, during the highest hydrologic mass loading. This magnitude of land subsidence remained almost consistent through May 2018. Following the end of the rainy season, the land had started an isostatic rebound by the end of May 2018, which recovered 50% of the subsiding land by August 2018. In the meantime, the GPS station showed stable movement by the end of April 2018 through May. However, the subsidence continued gradually from June until August 2018, reaching 9.0 mm ( Figure 5). The GPS measurements indicated a surface rebound in September 2019. The HYDL model is a product of all summations of the mass loading of all storage components that respond in time to changes in rainfall rates. The GPS displacement represents the direct movement of the earth as a result of increasing mass loading as the rain gradually increases. The InSAR results demonstrated further the land motion estimated by both the HYDL model and the GPS, but in different magnitudes. Regarding the change in SM, the dryness of the soil was persistent during the drought period of 2015. The 2015 drought was followed by the relatively long wet season (October to December) and continued with the MAM 2016. From the summer of 2016 to the MAM 2018, the prolonged drought continued, as suggested by the estimates of both the surface soil moisture and the root zone moisture contents from the two models ( Figure 3). The wetness of the soil reached its peak from April to May of 2018 with a maximum of 87 mm and 51 mm as estimated by the VIC and the NOAH-3.3 models, respectively. Similarly, the moisture of the root zone reached its peak in the same periods with relatively high values of 55 mm (Figure 3). The drastic deficit in the root zone soil moisture was estimated using the NOAH-3.3 model only. The hydrological fluxes are aligned with the recorded flooding that occurred through the 2018 MAM season. Many consequences were recorded in association with the flooding of this season. One of these consequences was the land deformation that we monitored and measured in the study area located in central Kenya.

Land Deformation Related to MAM Season
The land motions associated with the exceptional rainy MAM season demonstrate varying magnitudes of subsidence, shown by both GPS and InSAR measurements (Figures 5 and 6). The precipitation rates between 2017 and 2018 were plotted against the estimates of the load deformations, as determined from the hydrological loading model (HYDL) and the nonlinear trend of the GPS station, filtered by the singular spectral analysis (SSA) method. For the first three months (December 2017 to February 2018), the precipitation level showed an average of 23 mm/month. Throughout this temporal interval, the surface displacement was marginal, with an average of 1.0 mm uplift as determined from both the HYDL and the GPS station. Starting from March 2018, the precipitation increased markedly, with a magnitude of~180 mm, and it continued to increase, reaching its highest magnitude of~250 mm in April. This significant increase in rainfall greatly developed terrestrial water storage, resulting in mass load-induced land subsidence. Both the HYDL model and the GPS measurements reflected sharp land subsidence of 5 mm through March and April, during the highest hydrologic mass loading. This magnitude of land subsidence remained almost consistent through May 2018. Following the end of the rainy season, the land had started an isostatic rebound by the end of May 2018, which recovered 50% of the subsiding land by August 2018. In the meantime, the  Regarding the change in SM, the dryness of the soil was persistent during the drought period of 2015. The 2015 drought was followed by the relatively long wet season (October to December) and continued with the MAM 2016. From the summer of 2016 to the MAM 2018, the prolonged drought continued, as suggested by the estimates of both the surface soil moisture and the root zone moisture contents from the two models (Figure 4). The wetness of the soil reached its peak from April to May throughout the studied frame, indicating no significant land deformation ( Figure 6). However, some relatively high values can be noticed in scattered spots throughout the analyzed frame. Their scattered pattern may or may not indicate potential significance. In 2018, the InSAR measurements presented high values of land subsidence ranging from 20.0 mm to 40.0 mm/year. Interestingly, the subsiding area is wide and has a NE-SW trend. In 2019, it seemed that most of the subsided areas had recovered and showed a smaller magnitude of land deformation, comparable to that seen in 2017. The area of the highest subsidence is occupied mainly by the Kapiti-Athi and Thika plains.   relatively high values can be noticed in scattered spots throughout the analyzed frame. Their scattered pattern may or may not indicate potential significance. In 2018, the InSAR measurements presented high values of land subsidence ranging from 20.0 mm to 40.0 mm/year. Interestingly, the subsiding area is wide and has a NE-SW trend. In 2019, it seemed that most of the subsided areas had recovered and showed a smaller magnitude of land deformation, comparable to that seen in 2017. The area of the highest subsidence is occupied mainly by the Kapiti-Athi and Thika plains.   The InSAR results in Figure 6 show that the subsidence was widespread in 2018 across the analyzed frame, in contrast with low deformation rates recorded in 2017 and 2019. Through the temporal interval of 2016-2017, the values of the InSAR measurements were consistently low throughout the studied frame, indicating no significant land deformation ( Figure 6). However, some relatively high values can be noticed in scattered spots throughout the analyzed frame. Their scattered pattern may or may not indicate potential significance. In 2018, the InSAR measurements presented high values of land subsidence ranging from 20.0 mm to 40.0 mm/year. Interestingly, the subsiding area is wide and has a NE-SW trend. In 2019, it seemed that most of the subsided areas had recovered and showed a smaller magnitude of land deformation, comparable to that seen in 2017. The area of the highest subsidence is occupied mainly by the Kapiti-Athi and Thika plains.

Discussion
Land subsidence can be a naturally driven (e.g., tectonic movements [54] and permafrost processes [55]), or a human-driven process (e.g., aquifer compaction [56][57][58][59] and hydrocarbon extraction [60,61]). Sinkholes are areas that subside naturally due to different mechanisms across many geographical locations [62,63]. In Kenya, the relationship between land subsidence and sinkholes has been generally considered in the context of geologic and geomorphic studies. The study area where land subsidence took place in response to the MAM rainy season comprises four main regions: the Kapiti, Athi, Embakasi, and Thika plains [64] (Figure 7). This area is comprised of a low-lying plain with low slopes toward the west and northwest [65]. In response to the volcanic activities associated with the African Rift Valley, the low-lying plains have been filled by fissure lava eruptions [64,66]. These plains were formed by a variety of volcanic eruptions during different temporal stages [64,66]. In response to high temperature and annual rainfall ranging from~700 to 1000 mm, the volcanic rocks have been eroded, forming clay soils in different colors and at different thicknesses [67,68]. Sandy clay loams also occur as pockets and bands in the Athi river area [67]. Tectonically, the faulting in the rift zone and its margins comprises three major periods: major faults, grid faults, and cross faults [68]. The major faults have N-S trends and preceded the volcanic activities in the region. The grid faults developed in the Middle Pleistocene era and marked the end of the main volcanic events in the area. The cross faults and minor faults have affected all other structural features in the Nairobi area [68]. Below fault scarps or at fault-line intersections, sinkholes exist [68]. Under heavy rain, soils that fill the sinkhole subside as the bottom materials are flushed downward [68]. Ngecu and Nyambok [12] first monitored the ground subsidence at the Nakuru area in the central Rift Valley northwest of the study area. This ground subsidence occurs along parallel fault zones in several areas in Nakuru township during or after heavy rain events [5,12]. They attributed the ground subsidence to the oversaturation of the common volcanoclastic sediments, resulting in a reduction of the shear strengths of these sediments. The mass hydrologic loads of the saturated sediments induced subterranean erosion along faults, resulting in eventual falling of these unconsolidated sediments into subsurface channels, forming the sinkholes [12] ( Figure S1). Soil subsidence and formation of sinkholes may also occur in areas where no faults exist [68] when subsidence occurs within soft and incohesive soils after rain. One of the main characteristics of the Rift Valley zone and its environs is the occurrence of daily tremors and seismic activities [5,68]. The relationship between land subsidence and seismic activities was studied in the land subsidence of 1985 in the Nakuru area. The study indicated that no lateral displacements were detected along the subsidence zones of the 1985 incident [5]. The interpretation of the occurrence of land subsidence as a result of subterranean erosion during and after heavy rains has been outweighed in the Nakuru area [5,68]. Considering the magnitude and distribution of land subsidence, we may speculate that the subsidence during the rainy seasons could trigger future seismicity rather than increasing the spatial extents and magnitude of flooding in the study area. The link between heavy rainfall and the occurrence of earthquakes was primarily established by Huang et al. [69] in southern California along the San Andreas fault. They observed a correlation between the intense dry/wet cycle (and therefore the surface load) and the occurrences of earthquakes with magnitudes greater than six. Furthermore, land subsidence can potentially intensify flooding by changing flood plain boundaries and base flood drainage [22,70] and by lowering flood control structures [71]. The blue dots in Figure 8 represent the flooded areas between March and June 2018. The maximum spatial distribution of the points was mapped in the second week of April. Before March or after May 2018, the scattering points are limited and represent only the normal water bodies in the study area. From the comparison between the HDYL model and the GPS measurements and the comparison between the InSAR measurements and flood mapping extents using radar data, it is possible to conclude that there is a potential relationship between heavy rainfall and ground subsidence in the study area. The occurrence of heavy rainfall events and flooding over volcanoclastic regions covered mainly by clay soils would intensify the sinkhole processes in which subterranean channel erosion actively flushes bottom materials downwards. Accordingly, land subsidence can take place. The effects of the spatial locations, the spatial pattern of the active faults, and the seismic activities of these faults on the land subsidence in the study area require more detailed investigation. Although this paper may have helped provide an empirical example for the climate-driven deformations in Kenya, many limitations can be outlined. For instance, the absence of dense GPS stations prohibited proper validation of InSAR measurements of land deformation. In addition, the unavailability of in situ rainfall and flood gauge stations to substantiate the various hydrological parameters extracted from earth observation satellites and hydrologic models' data is another limitation. In addition, given the severity of the 2018 rainy season and the widespread viscoelastic sediments and intense faulting along the East African rift system, it is likely expected that similar deformations may have occurred in other regions. However, it is beyond our capability to address these phenomenal earth processes across the vast areas covered by the heavy rains of the MAM season. To investigate the relationship between the flooding of the MAM season and land subsidence in the studied area (Figure 7), the flooded areas were estimated using the scattering coefficient technique, as mentioned in Section 2.2. The spatial correlation implied a strong relationship between the areas of high magnitudes of land subsidence ( Figure 6) and those areas that experienced intense flooding during the MAM (Figure 8).

Conclusions
Using SMAP satellite and hydrological models, we mapped the surplus precipitation and soil moisture that accompanied the long MAM rainy season of 2018 across the Horn of Africa. In Kenya, for the first time, we have shown the solid-earth responses to the MAM season using loading models, GPS, and interferometric analysis of the synthetic aperture radar (InSAR) data. We found that the MAM rainy season resulted in widespread land subsidence up to > 20 mm in Kenya during 2018 based on the SBAS InSAR with ~500 interferograms. The rate of land subsidence measured from InSAR is higher than the estimated rates from both the nearby GPS site and the hydrological loading model. The areas showing a high magnitude of ground subsidence were spatially correlated with the intensively flooded areas during the MAM season of 2018. Recovery had occurred by the end of 2018 and the beginning of 2019. Considering the magnitude and distribution of the land subsidence, we The blue dots in Figure 8 represent the flooded areas between March and June 2018. The maximum spatial distribution of the points was mapped in the second week of April. Before March or after May 2018, the scattering points are limited and represent only the normal water bodies in the study area. From the comparison between the HDYL model and the GPS measurements and the comparison between the InSAR measurements and flood mapping extents using radar data, it is possible to conclude that there is a potential relationship between heavy rainfall and ground subsidence in the study area. The occurrence of heavy rainfall events and flooding over volcanoclastic regions covered mainly by clay soils would intensify the sinkhole processes in which subterranean channel erosion actively flushes bottom materials downwards. Accordingly, land subsidence can take place. The effects of the spatial locations, the spatial pattern of the active faults, and the seismic activities of these faults on the land subsidence in the study area require more detailed investigation. Although this paper may have helped provide an empirical example for the climate-driven deformations in Kenya, many limitations can be outlined. For instance, the absence of dense GPS stations prohibited proper validation of InSAR measurements of land deformation. In addition, the unavailability of in situ rainfall and flood gauge stations to substantiate the various hydrological parameters extracted from earth observation satellites and hydrologic models' data is another limitation. In addition, given the severity of the 2018 rainy season and the widespread viscoelastic sediments and intense faulting along the East African rift system, it is likely expected that similar deformations may have occurred in other regions. However, it is beyond our capability to address these phenomenal earth processes across the vast areas covered by the heavy rains of the MAM season.

Conclusions
Using SMAP satellite and hydrological models, we mapped the surplus precipitation and soil moisture that accompanied the long MAM rainy season of 2018 across the Horn of Africa. In Kenya, for the first time, we have shown the solid-earth responses to the MAM season using loading models, GPS, and interferometric analysis of the synthetic aperture radar (InSAR) data. We found that the MAM rainy season resulted in widespread land subsidence up to > 20 mm in Kenya during 2018 based on the SBAS InSAR with~500 interferograms. The rate of land subsidence measured from InSAR is higher than the estimated rates from both the nearby GPS site and the hydrological loading model. The areas showing a high magnitude of ground subsidence were spatially correlated with the intensively flooded areas during the MAM season of 2018. Recovery had occurred by the end of 2018 and the beginning of 2019. Considering the magnitude and distribution of the land subsidence, we may speculate that the subsidence during the rainy seasons could trigger future seismicity. From economic and societal perspectives, the solid-earth responses to the surface loading of the wet and dry cycles have a severe environmental and humanitarian impact. Governments are advised to adopt a strategic plan when constructing roads and building across draining systems that may not function adequately to drain rainwater.