Investigating the Spatio-Temporal Variation of Soil Moisture and Agricultural Drought towards Supporting Water Resources Management in the Red River Basin of Vietnam

: The increasing serious droughts recently might have signiﬁcant impacts on socioeconomic development in the Red River basin (RRB). This study applied the variable inﬁltration capacity (VIC) model to investigate spatio-temporal dynamics of soil moisture in the northeast, northwest, and Red River Delta (RRD) regions of the RRB part belongs to territory of Vietnam. The soil moisture dataset simulated for 10 years (2005–2014) was utilized to establish the soil moisture anomaly percentage index (SMAPI) for assessing intensity of agricultural drought. Soil moisture appeared to co-vary with precipitation, air temperature, evapotranspiration, and various features of land cover, topography, and soil type in three regions of the RRB. SMAPI analysis revealed that more areas in the northeast experienced severe droughts compared to those in other regions, especially in the dry season and transitional months. Meanwhile, the northwest mainly suffered from mild drought and a slightly wet condition during the dry season. Different from that, the RRD mainly had moderately to very wet conditions throughout the year. The areas of both agricultural and forested lands associated with severe drought in the dry season were larger than those in the wet season. Generally, VIC-based soil moisture approach offered a feasible solution for improving soil moisture and agricultural drought monitoring capabilities at the regional scale. the VIC macroscale hydrological model was applied for simulation of the at a grid of 0.1 × 0.1 during the 10-year of for the and RRD regions of the RRB part to the territory Our hydrological model-based soil moisture approach offers an alternative solution for soil moisture observation that could address the currently limited in-situ soil moisture monitoring capabilities at the regional scale in Vietnam and other developing countries. The VIC simulation successfully captured the spatio-temporal variability of for the three regions in the RRB. This illustrated how co-varied evapotranspiration, also with of topography, type in


Introduction
The Red River basin (RRB) is one of the largest river basins in the world and the second largest river basin in Vietnam, which crosses three countries including China, Laos, and Vietnam. The RRB has a tropical monsoon climate with distinct dry and wet seasons. Droughts and floods are the major natural disasters in the RRB [1]. Droughts occur frequently in the RRB, especially during the dry season in which many areas of the RRB suffer serious water shortage that significantly affects economic activities (i.e., domestic water supply and irrigation water). For example, in the dry season of 2005, the water level in the Red River decreased to the lowest level for the last 100 years. Moreover, droughts in the RRB have also sped up the seawater intrusion in both surface water and underground water aquifers in the coastal provinces. Under the context of the increasing impact of climate change in recent years, especially during the dry season, the increased average temperature has led to greater evaporation. The rainfall variation has also become more extreme in which rainfall occurred more frequently during the rainy season with very high intensity, whereas it substantially reduced during the dry season. The frequent and serious droughts in the RRB have unfavorable impacts on agriculture, ecosystem, hydrology, and socioeconomic development. Particularly, the Red River Delta (RRD) unanswered. To address these knowledge gaps, this study employed the VIC model to explore the spatio-temporal dynamics and characteristics of soil moisture in the RRB. The novelty of this study lies in utilizing the VIC hydrological model, which not only addresses the lack of soil moisture monitoring network but also derives the soil moisture information at the regional scale for the RRB. The soil moisture dataset simulated for the 10 years (2005-2014) was then utilized to establish the soil moisture anomaly percentage index (SMAPI; [7,18]), which could be used to enhance the display of spatio-temporal variation of drought and to assess the intensity of agricultural drought. Furthermore, this study investigated the impacts of meteorological variables (precipitation and temperature) and land cover properties, topography, and soil types on the spatial and temporal variation of soil moisture and drought in different regions in the RRB. Our study is the first ever and comprehensive one in the fields of soil moisture and agricultural drought over the RRB part located in the territory of Vietnam. The key objectives of this study were to (i) explore the spatio-temporal variation and characteristics of soil moisture in the RRB; (ii) investigate the sensitivity of soil moisture to seasonal meteorological conditions, land cover properties, topography, and soil type in different regions in the RRB; and (iii) analyze the intensity of agricultural drought in relation to land uses in the RRB.

Study Area
The RRB is shared by three countries including China, Vietnam, and Lao's PDR with a total watershed area about 169,020 km 2 with 81,240 km 2 (48%) in China, 1100 km 2 (0.65%) in Laos, and 86,660 km 2 (51.35%) in Vietnam [19]. The Red River with its main branch of approximately 1140 km in length, rises in a mountainous region of Southeastern China, flows into Vietnam near Lao Cai province, and flows eastward into the Tonkin Bay. In the territory of Vietnam, the RRB includes the majority of the three ecological regions in the north of Vietnam, namely northeast, northwest, and RRD, which are denoted as the region (I), (II), and (III), respectively in Figure 1. The climate conditions in the RRB are dominated by the East-Asian monsoon with the rainy (wet) season (from May to September) representing 85-90% of the annual rainfall, and the dry season (from November to March next year) representing only 10-15% of the annual rainfall. The annual average rainfall in the territory of Vietnam could largely vary from 1200 to 4800 mm. The transitional periods are April and October. This study focused on the three regions (I), (II), and (III) of the RRB, which belong to the territory of Vietnam.

VIC Soil Moisture Simulation
The VIC is a macroscale and semidistributed hydrological model based on the grid cell platform. The VIC model can solve both water and energy balance equations to estimate various parameters in the water cycle of the Earth including runoff, soil moisture, infiltration, and evapotranspiration at each grid cell with daily or subdaily interval [8,19].

VIC Soil Moisture Simulation
The VIC is a macroscale and semidistributed hydrological model based on the grid cell platform. The VIC model can solve both water and energy balance equations to estimate various parameters in the water cycle of the Earth including runoff, soil moisture, infiltration, and evapotranspiration at each grid cell with daily or subdaily interval [8,19]. In another word, these hydrometeorological parameters were derived from the VIC diagnostic output in this study. Similar to other hydrological models, the VIC model requires four basic parameters as input data including (1) topography, (2) meteorological forcing data (e.g., precipitation, maximum and minimum air temperature, and wind speed), (3) land cover land use, and (4) soil properties.
In our previous research [19], the VIC model was successfully applied to simulate the streamflow in the RRB for the period from 2005 to 2014. In detail, the 0.1 degree of resolution was used to establish the grid cell framework of the VIC model in the RRB. For data processing, the basin mask was created basing on the ESRI grid with a non-proprietary ASCII. The structure of ASCII format stores the geographic attribute of the RRB including the total number of rows and columns, which was 64 and 56, respectively. The starting point that displayed the easting and northing coordinates under the decimal degree format were 100.05 and 19.95, respectively. The unit of space was 0.1 by 0.1 decimal degree, which was the same as the geographic coordinate system above. Under the support of ASCII to the Raster tool in ArcGIS software, the raster image was generated to cover the RRB. It was a mandatory platform to integrate the other datasets into the VIC model at 0.1-degree resolution. With respect to topography, the digital elevation model (DEM) from the NASA's shuttle radar topographic mission (SRTM) (https://www2.jpl.nasa.gov/srtm/dataprod.htm (accessed on 27 January 2021)) 90 m digital elevation database version 1.4 with raster format was collected. For post-processing the resolution of this raster image, the extraction by mask tool was used to take out the elevation of each grid cell in the RRB following the basin mask, which was created in the previous step. The daily meteorological data (maximum and minimum temperature, precipitation, and wind speed) was obtained from the National Centers for Environmental Information (NCEI) with the Global Surface Summary of the Day product (https://www.ncei.noaa.gov/access/metadata/landing-page/bin/iso?id= gov.noaa.ncdc:C00516 (accessed on 27 January 2021)) to present the climate characteristic over the RRB. Our previous study [19] noted more attention to process datasets such as collect the location of the stations in the list of monitoring point network in GSOD product over the RRB and extract the meteorological dataset at each station. In order to generate the meteorological database under grid cell structure at 0.1 degree, the inverse distance weighting method was employed. The land use land cover data was taken from United States Geological Survey (USGS) with global land cover characterization (GLCC) product at 1 km and used to describe vegetation properties (e.g., feature, classification, and ecoregions interpretation) over the RRB. Meanwhile, the soil characteristic was extracted from the Harmonized World Soil Database (HWSD) of the Food and Agriculture Organization of the United Nations (FAO, http://www.fao.org/soils-portal/soil-survey/soil-maps-anddatabases/harmonizedworld-soil-database-v12/en (accessed on 27 January 2021)). This product was a combined multisource of soil information worldwide at 30 arc-seconds of a spatial scale with over 16,000 different soil mapping elements. Similar to post-processing with DEM, the GLCC and HWSD products were consolidated with the VIC model under the grid platform at 0.1 degree. Finally, each grid cell platform stored the compulsory information such as geography attribute, meteorological forcing data, land cover and land use, and soil texture to produce the input data of the VIC model.
The VIC model calibration and validation were performed for the two periods of 2005-2009 and 2010-2014, respectively using the streamflow data obtained from the hydrological monitoring network along the Red River (the main river of RRB). For calibration and validation procedures, the sensitivity parameters of the VIC model were adjusted including b_infilation, Ds, Ds_max, and Ws [10]. In our previous study [19], the reasonable parameters of the VIC model to perform the goodness fit between the estimated and e Sustainability 2021, 13, 4926 5 of 15 measured streamflows were described in detail. The key metrics for model performance evaluation showed that the VIC model could capture the hydrologic dynamic over the RRB with the Pearson's correlation coefficient reached 0.75 and the index of agreement approached 0.85.
Using the simulation outputs of VIC model, the simulated soil moisture dataset for a 10-year period of 2005-2014 was extracted for analyzing the agricultural droughts over the RRB in this study. The soil moisture dataset was processed and presented for the dry season, wet season, and transitional periods.

Soil Moisture Anomaly Percentage Index (SMAPI)
Many studies have used precipitation-based drought indices to measure the drought severity in which the two common ones have been the standard precipitation index (SPI; [20]) and the Palmer drought severity index (PDSI; [14]). However, according to [21], droughts occurring at the basin scale were influenced by the decreased precipitation and other factors governing water balance state. Particularly, for the agricultural sector, soil moisture has been considered as the onset of droughts [22] and more vital than precipitation as different stages of crop growth are largely impacted by soil moisture. The four indices that are typically used for studying agricultural drought include the soil moisture index (SMI), the crop drought index (CDI), the crop moisture index (CMI), and the soil moisture anomaly percentage index (SMAPI). The SMI reflects the relative dryness of soil, which is calculated as the percentage of actual soil moisture attributed to the field capacity, while the CDI is estimated basing on the ecological and physiological responses of crops to drought using satellite-based measurements over a short period of time. Both the SMI and CDI are more appropriate for analyzing drought in small areas with specific soil types and crops. The CMI is determined as the sum of evapotranspiration deficit and soil demand, being a byproduct of PDSI with the assumption that land use, land cover, and vegetation are soil properties that are evenly distributed in the whole climatic zone, which may not reflect the real conditions in the studied region. Whereas, the SMAPI has been developed for characterization of agricultural drought at the regional scale [7], and thus is a more suitable metric for the purpose of this study. The SMAPI is defined as the following: where, SM and SM* express the current value of the simulated daily soil moisture and the mean soil moisture of the study period, respectively, using the standard VIC mode diagnostic of ET. The classification of SMAPI values was made following the approach by [7,12] and the codes were assigned to represent different drought/wet levels for more visual analysis of spatial SMAPI maps as presented in Table 1.

Spatio-Temporal Variation of Soil Moisture in Relation to Meteorological Variables
The spatio-temporal variation of soil moisture, meteorological variables (precipitation and air temperature), and evapotranspiration during the dry and wet seasons, and transitional months for the period of 2005-2014 in the RRB are shown in Figure 2. In order to gain an improved understanding of the impacts of meteorological variables (precipitation and air temperature) and evapotranspiration on the spatio-temporal variation of soil moisture, the seasonal covariability of soil moisture and the above variables was analyzed. Our results revealed that there were moderate to good agreements of spatial correlations between soil moisture vs. precipitation, air temperature, and evapotranspiration during the dry and wet seasons, and transitional months for both the northeast and northwest regions, except those for some areas of the RRD (Figure 3). The mean air temperature in the wet season was much higher than those in the dry season and transitional months, subsequently leading to the larger evapotranspiration in the wet season than those in the dry season and transitional months. Both precipitation and air temperature had significant impacts on the soil moisture changes in the northeast and northwest regions, which was reported in other studies [6,9,13,17,[23][24][25][26][27]. However, the impact of temperature on soil moisture appeared to be larger than that of precipitation in the two regions, which was also reported in the previous study [28]. The relatively weak correlation results in the RRD region suggested that not only precipitation, air temperature, and/or evapotranspiration, the other unique meteorological conditions in the RRD might also have significant impacts on the seasonal variation of soil moisture in this region. During the early spring period (Feb-Mar), the typical weather condition of the RRD is characterized by the frequent drizzle with extremely humid air. In addition, the prevailing easterly and southeasterly winds for most of the time in a year could bring a huge amount of In general, there was a strongly spatial consistency among the soil moisture, precipitation, mean temperature, and evapotranspiration maps. As seen in Figure 2, there was a clear seasonal regularity of soil moisture, especially for the northeast and northwest regions of the RRB with the highest values during the wet season, followed by those during the transitional period and dry season. The soil moisture in the RRD showed a less pronounced seasonal trend, although the soil moisture during the wet season was slightly higher than those during the transitional months and dry season. The spatial maps of precipitation, mean temperature, and evapotranspiration indicated an obviously increasing pattern for these variables during the wet season for all of the regions in the RRB. During the wet season, in the upper part of the northeast region and in the lower part of northwest region, the lower precipitation amount might cause the decreased recharge in soil layers, and subsequently contribute to the lower soil moisture in these areas in comparison to those in the other areas of the RRB. Moreover, the significantly increasing evapotranspiration induced by the increasing air temperature during the wet season has also contributed to the lower soil moisture in these areas in comparison to those in the remaining areas of the RRB. Clearly, on average, the spatial distribution analysis revealed that the RRD always had the maximum soil moisture throughout the year.
In order to gain an improved understanding of the impacts of meteorological variables (precipitation and air temperature) and evapotranspiration on the spatio-temporal variation of soil moisture, the seasonal covariability of soil moisture and the above variables was analyzed. Our results revealed that there were moderate to good agreements of spatial correlations between soil moisture vs. precipitation, air temperature, and evapotranspiration during the dry and wet seasons, and transitional months for both the northeast and northwest regions, except those for some areas of the RRD (Figure 3). The mean air temperature in the wet season was much higher than those in the dry season and transitional months, subsequently leading to the larger evapotranspiration in the wet season than those in the dry season and transitional months. Both precipitation and air temperature had significant impacts on the soil moisture changes in the northeast and northwest regions, which was reported in other studies [6,9,13,17,[23][24][25][26][27]. However, the impact of temperature on soil moisture appeared to be larger than that of precipitation in the two regions, which was also reported in the previous study [28]. The relatively weak correlation results in the RRD region suggested that not only precipitation, air temperature, and/or evapotranspiration, the other unique meteorological conditions in the RRD might also have significant impacts on the seasonal variation of soil moisture in this region. During the early spring period (Feb-Mar), the typical weather condition of the RRD is characterized by the frequent drizzle with extremely humid air. In addition, the prevailing easterly and southeasterly winds for most of the time in a year could bring a huge amount of moist air from the east sea westward to the coastal provinces in the RRD. These factors would contribute significantly in maintaining a stable humid condition, and simultaneously high soil moisture in the RRD throughout the year. Besides the influences of meteorological variables, the spatio-temporal variation of soil moisture in the RRB could be also largely impacted by the other factors (i.e., land cover, seasonal variation of vegetation, topography, and soil types), which would be further explored in the following section.

Relationship between Soil Moisture vs. Land Cover, Topography, and Soil Type
In addition to the impacts of meteorological variables, according to [4], soil moisture variation is also impacted by other factors such as (i) land cover type and vegetation, which affect the runoff, infiltration, and evapotranspiration through root water uptake [4,[29][30][31]; (ii) topography, which influences the runoff and infiltration [32,33]; and (iii) soil type, which impacts the water-holding capacity [34,35]. Vegetation type, density, and uniformity are key characteristics that contribute to soil moisture changes at different spatial and temporal scales. The partial and full canopy coverages were associated with the high and low variability of soil moisture [4,36]. Topography could contribute significantly to the spatial variation of soil moisture at different scales [36][37][38]. Variations in key topographical features such as relative elevation, slope, curvature, and aspect have impacts on the variation of soil moisture near the surface. The variation of soil moisture at the hill-slope and small catchment scales could be influenced by water-routing processes, heterogeneity in soil characteristics, and vegetation [4]. The previous studies [36][37][38][39][40][41] reported that soil heterogeneity, which was determined by the variations in soil texture, porosity, structure, and organic matter content could influence the distribution of soil moisture through governing soil water transmission and retention [42]. Due to the variation in soil particle or pore sizes, substantial variation of soil moisture could occur over very small areas, and subsequently affect the transient downward and upward soil water fluxes [4]. In order to further explore the combined effects of these factors on the spatio-temporal variation of soil moisture, the maps of land cover and vegetation, topography, and soil type in the RRB are shown in Figure 4. moist air from the east sea westward to the coastal provinces in the RRD. These factors would contribute significantly in maintaining a stable humid condition, and simultaneously high soil moisture in the RRD throughout the year. Besides the influences of meteorological variables, the spatio-temporal variation of soil moisture in the RRB could be also largely impacted by the other factors (i.e., land cover, seasonal variation of vegetation, topography, and soil types), which would be further explored in the following section.

Relationship between Soil Moisture vs. Land Cover, Topography, and Soil Type
In addition to the impacts of meteorological variables, according to [4], soil moisture variation is also impacted by other factors such as (i) land cover type and vegetation, which affect the runoff, infiltration, and evapotranspiration through root water uptake [4,[29][30][31]; (ii) topography, which influences the runoff and infiltration [32,33]; and (iii) soil With respect to land cover and topography, the northeast region has the most dominant land cover and vegetation type is the mixed forest, followed by the cropland/grassland mosaic (Figure 4a) and the elevation decreasing from north to south (Figure 4b). The uppermost areas near the border with China are mainly covered by the cropland/grassland mosaic and have the mountainous and steep landform with an average elevation of 1000-1200 m. The cropland/grassland mosaic could favor soil evaporation and has the weak effect of slowing down the runoff [17,43] while the mountainous landform could facilitate the surface runoff and reduce the infiltration, which resulted in the low soil moisture in these areas for all three seasons (Figure 2). A similar condition could be observed for the lowermost areas of the northeast region with the mainland cover type is the cropland/grassland mosaic. Meanwhile, most of the remaining areas of the northeast region are covered by the mixed forest. The previous studies [17,[44][45][46][47] have shown that forests Sustainability 2021, 13, 4926 9 of 15 could slow down the runoff and thus increase the infiltration, and subsequently increase the soil moisture. Thus, it is inferred that the mixed forest coupled with the relatively steep landform could be the key factors affecting the distribution of soil moisture in these areas. In the northwest region, there are different land cover types. The upper and lower western parts of the northwest region are mainly covered by the evergreen broadleaf forest and deciduous broadleaf forest, respectively. Due to the high shading effect of the year-round full canopy of the evergreen broadleaf forest with dense vegetation, the soil evaporation could be reduced significantly, which could lead to the soil moisture in the upper western part being relatively stable and high [48], especially during the wet season. The high evapotranspiration in the upper western part of the northwest region was probably attributed to the increased transpiration. Whereas the lower western part of the northwest region is mainly covered by the savanna and deciduous broadleaf forest with the leafless period occurring in the dry season, which could induce the increased soil evaporation, and then the decreased soil moisture in these areas. Meanwhile, the remaining areas (the eastern and southern parts) of the northwest region are mainly covered by the mixed forest, which could reduce the soil evaporation and maintain relatively stable soil moisture in these areas. In the RRD region, the major land cover type is the irrigated cropland and pasture. The studies by [17,49,50] have shown that the regulation of surface runoff of cropland is poorer than that of forests. Furthermore, this region has a plain landform, which does not have a strong effect on the formation of surface runoff likes those in the other regions in the RRB. As a result, instead of occurrence of the surface runoff, a part of precipitation could be infiltrated into the soil making the soil moisture increased. On the other hand, as the largest agricultural production area in the RRB, the intensive irrigation activity was also expected to be one of the key factors contributing to the high soil moisture in the RRD. type, which impacts the water-holding capacity [34,35]. Vegetation type, density, and uniformity are key characteristics that contribute to soil moisture changes at different spatial and temporal scales. The partial and full canopy coverages were associated with the high and low variability of soil moisture [4,36]. Topography could contribute significantly to the spatial variation of soil moisture at different scales [36][37][38]. Variations in key topographical features such as relative elevation, slope, curvature, and aspect have impacts on the variation of soil moisture near the surface. The variation of soil moisture at the hillslope and small catchment scales could be influenced by water-routing processes, heterogeneity in soil characteristics, and vegetation [4]. The previous studies [36][37][38][39][40][41] reported that soil heterogeneity, which was determined by the variations in soil texture, porosity, structure, and organic matter content could influence the distribution of soil moisture through governing soil water transmission and retention [42]. Due to the variation in soil particle or pore sizes, substantial variation of soil moisture could occur over very small areas, and subsequently affect the transient downward and upward soil water fluxes [4].
In order to further explore the combined effects of these factors on the spatio-temporal variation of soil moisture, the maps of land cover and vegetation, topography, and soil type in the RRB are shown in Figure 4. With respect to land cover and topography, the northeast region has the most dominant land cover and vegetation type is the mixed forest, followed by the cropland/grassland mosaic (Figure 4a) and the elevation decreasing from north to south (Figure 4b). The uppermost areas near the border with China are mainly covered by the cropland/grassland mosaic and have the mountainous and steep landform with an average elevation of 1000-1200 m. The cropland/grassland mosaic could favor soil evaporation and has the weak effect of slowing down the runoff [17,43] while the mountainous landform could facilitate the surface runoff and reduce the infiltration, which resulted in the low soil moisture in these areas for all three seasons (Figure 2). A similar condition could be observed for the lowermost areas of the northeast region with the mainland cover type is the cropland/grassland mosaic. Meanwhile, most of the remaining areas of the northeast region are covered by the mixed forest. The previous studies [17,[44][45][46][47] have shown that forests could slow down the runoff and thus increase the infiltration, and subsequently increase the soil moisture. Thus, it is inferred that the mixed forest coupled with the relatively steep landform could be the key factors affecting the distribution of soil moisture in these areas. In the northwest region, there are different land cover types. The upper and lower western parts of the northwest region are mainly covered by the evergreen broadleaf forest and deciduous broadleaf forest, respectively. Due to the high shading effect of the year-round full canopy of the evergreen broadleaf forest with dense vegetation, the soil evaporation could be reduced significantly, which could lead to the soil moisture in Regarding soil types, the most dominant soil is the sandy clay loam, followed by the loam and loamy sand for both the northeast and northwest regions. Meanwhile, the RRD has a mix of clay (light), loam, and loamy sand (Figure 4c). Storage and redistribution of soil water are influenced by soil texture and structure. Soils that have higher clay contents would have a larger pore space and a greater water holding capacity. It was reported that smaller particle soils (i.e., clay) have larger surface areas than those of larger particle soils (i.e., sand), which allow them to hold more water [41]. Therefore, the presence of more clay-rich soil with the higher water-holding capacity in the RRD could be a key factor contributing to the more stable and higher soil moisture in this region compared to those in the northeast and northwest regions that have soil types with the lower water-holding capacity. On the other hand, there are dense river network systems (Figure 1) in the RRD, which could increase the groundwater recharge and maintain stable soil moisture throughout the year in this region. However, as this study has not taken into account the groundwater recharge mechanism for the VIC model design, and therefore the relationship between groundwater recharge and variation of soil moisture in the RRD should be further investigated in the future studies.

Spatio-Temporal Variation of SMAPI and Implications for Agroforestry
The seasonally averaged SMAPI maps for the RRB during the period 2005-2014 are displayed in Figure 5. It is seen that the SMAPI maps were able to capture the spatiotemporal variability of drought conditions in different regions of the RRB. It was generally found that the northeast region had different drought intensities (mild and severe) and near normal condition during the dry season and transitional months, respectively. Among the regions in the RRB, the northeast was also found to be the region with more areas experienced severe droughts compared to those in other regions, especially in the dry season and transitional months. During the wet season, more areas in the northeast region were classified as having the slightly wet condition. The northwest region mainly suffered from mild drought and slightly wet condition during the dry season. Whereas, most of the areas in the northwest region had slightly/moderately wet conditions during the wet season and transitional months. Different from the northeast and northwest regions, the majority of the RRD consistently had moderately to very wet conditions for the whole year in which wet conditions during the dry season and transitional months were found more pronounced than that during the wet season. As drought and wet conditions have important implications for agroforestry activities, the relationships between land cover types and drought/wet intensities (code D1-D4 and W1-W4 categories for drought and wet conditions, respectively- Table 1) in the RRB were further analyzed. Based on the land cover and SMAPI maps, the percentage of land cover areas corresponding to different drought and wet intensities for the dry and wet seasons was statistically calculated and presented in Table 2. Clearly, land cover types could be separated into two major groups including: (1) those (i.e., urban and built-up land, water bodies, and herbaceous wetland) that were less affected by drought and wet conditions and (2) those (the remaining land cover types) that were variably affected by drought and wet conditions. It is generally seen that there were insignificant differences in the percentage of areas that suffered from the mild drought or slightly/moderately wet conditions between the dry and wet seasons for the group (1). In contrast, there were large variations in the percentage of areas impacted by different drought/wet intensities between the dry and wet seasons for the land cover types under the group (2).
In terms of the major agricultural and pasture lands, the irrigated cropland and pasture (15.78% of the RRB's total area, mainly located in the RRD), and cropland/grassland mosaic (15.49% of total area, mainly located in the lower part of the northeast region) showed the increased areas associated with the slightly-moderately wet conditions (W1-W2) in the wet season in comparing to those in the dry season. Meanwhile, the areas that suffered from the mild drought (D1) and severe drought (D3) for these two types of agricultural and pasture land in the dry season (4.31% of total area) were 1.28 times larger than those in the wet season (3.38% of total area). This implies that proper management of irrigation water should be considered for the dry season in the RRB.
With respect to the major forested lands, the deciduous broadleaf forest (6.6% of the As drought and wet conditions have important implications for agroforestry activities, the relationships between land cover types and drought/wet intensities (code D1-D4 and W1-W4 categories for drought and wet conditions, respectively- Table 1) in the RRB were further analyzed. Based on the land cover and SMAPI maps, the percentage of land cover areas corresponding to different drought and wet intensities for the dry and wet seasons was statistically calculated and presented in Table 2. Clearly, land cover types could be separated into two major groups including: (1) those (i.e., urban and built-up land, water bodies, and herbaceous wetland) that were less affected by drought and wet conditions and (2) those (the remaining land cover types) that were variably affected by drought and wet conditions. It is generally seen that there were insignificant differences in the percentage of areas that suffered from the mild drought or slightly/moderately wet conditions between the dry and wet seasons for the group (1). In contrast, there were large variations in the percentage of areas impacted by different drought/wet intensities between the dry and wet seasons for the land cover types under the group (2). In terms of the major agricultural and pasture lands, the irrigated cropland and pasture (15.78% of the RRB's total area, mainly located in the RRD), and cropland/grassland mosaic (15.49% of total area, mainly located in the lower part of the northeast region) showed the increased areas associated with the slightly-moderately wet conditions (W1-W2) in the wet season in comparing to those in the dry season. Meanwhile, the areas that suffered from the mild drought (D1) and severe drought (D3) for these two types of agricultural and pasture land in the dry season (4.31% of total area) were 1.28 times larger than those in the wet season (3.38% of total area). This implies that proper management of irrigation water should be considered for the dry season in the RRB.
With respect to the major forested lands, the deciduous broadleaf forest (6.6% of the RRB's total area, mainly located in the lower western part of the northwest region), the evergreen broadleaf forest (11.76% of total area, mainly located in the upper western part of the northwest region), and the mixed forest (25.75% of total area, mainly located in both the northeast and northwest regions), all showed the increasing patterns of both areas associated with the wet conditions in the wet season. Regarding the drought conditions, while the areas having the severe drought (D3) were not very different from the dry season (4.46% of total area) and the wet season (4.60% of total area), the areas suffered to the mild drought (D1) in the dry season were about 4.2, 1.8, and 1.9 times larger than those in the wet season for the deciduous broadleaf, evergreen broadleaf, and mixed forests, respectively. These results were consistent with the seasonal covariation of soil moisture and forest types in the RRB as discussed in the previous section. These findings imply that to effectively improve the water conservation capacity and to formulate the drought preparedness and response plan for the RRB region, the forest protection and management policies should pay adequate attention towards controlling soil moisture stabilization in the basin.

Conclusions
In this study, the VIC macroscale hydrological model was applied for simulation of the daily soil moisture at a grid resolution of 0.1 • × 0.1 • during the 10-year period of 2005-2014 for the northeast, northwest, and RRD regions of the RRB part belongs to the territory of Vietnam. Our hydrological model-based soil moisture approach offers an alternative solution for soil moisture observation that could address the currently limited in-situ soil moisture monitoring capabilities at the regional scale in Vietnam and other developing countries. The VIC simulation successfully captured the spatio-temporal variability of soil moisture for the three regions in the RRB. This study illustrated how soil moisture co-varied with not only precipitation, air temperature, and evapotranspiration, but also with various features of land cover type and vegetation, topography, and soil type in three regions of the RRB during the dry and wet seasons, and transitional months. However, it should be noted that the irrigation scheme was not considered in the VIC model set up in this study. In order to have a better understanding on the extent that irrigation may influence the water budget and soil moisture in the RRB, particularly the RRD with intensively irrigated land, future studies should take irrigation scheme into account for VIC model design.
Analysis of the spatio-temporal variation of soil moisture-based SMAPI revealed that more areas in the northeast region experienced severe drought compared to those in the northwest and the RRD regions, especially in the dry season and transitional months. Meanwhile, the northwest region mainly had mild drought and slightly wet condition during the dry season. Different from the northeast and northwest regions, the majority of the RRD appeared to have moderately to very wet conditions throughout the year. The areas of both the agricultural and forested lands associated with severe drought in the dry season were larger than those in the wet season, which provide the important implications for the water resource as well as forest protection and management policies in responding to droughts in the RRB.