Global hydrological models that include human activities are powerful tools for assessing water availability and use at global and continental scales. Such models are typically applied at a spatial resolution of 30 arcmin (approximately 50 km). In recent years, some
5 arcmin (9 km) applications have been reported but with numerous technical challenges, including the validation of calculations for more than
Global hydrological models have successfully incorporated major natural hydrological components and a number of important human activities, including agricultural water use (Döll and Siebert, 2002; Wada et al., 2013), industrial and domestic water use (Alcamo et al., 2003a; Flörke et al., 2013), reservoir operation (Hanasaki et al., 2006; Shin et al., 2019), groundwater pumping (Wada et al., 2010), desalination (Hanasaki et al., 2016), and aqueduct water transfer (Hanasaki et al., 2018b). Several global hydrological models are capable of simulating global water availability and use, including most of these human activities, for the purpose of water management. Such models include WaterGAP2 (Alcamo et al., 2003a; Müller Schmied et al., 2021), LPJmL (Rost et al., 2008; Schaphoff et al., 2018), H08 (Hanasaki et al., 2008a, 2018b), WBMplus (Wisser et al., 2010), PCR-GLOBWB (Wada et al., 2011; Sutanudjaja et al., 2018), HiGW-MAT (Pokhrel et al., 2012; Shin et al., 2019); VIC5 (Droppers et al., 2020), and the Community Water Model (CWatM; Burek et al., 2020). These models have been widely employed for water scarcity assessment (Alcamo et al., 2003b; Hanasaki et al., 2018a), climate change impact assessment (Schewe et al., 2014), water footprinting (Hanasaki et al., 2010; Dalin et al., 2017), estimation of the Earth's carrying capacity (Gerten et al., 2020), and many other purposes.
Among various challenges facing global hydrological models, increasing spatial resolution is perhaps the most straightforward and most difficult. For the past decade, most global hydrological simulations have been conducted at a spatial resolution of 30 arcmin (0.5
Recently, some models have been applied at spatial resolutions considerably finer than 30 arcmin. To the authors' knowledge, the first example was reported by Eisner and Flörke (2015). They developed the WaterGAP3 model, based on WaterGAP2 (Alcamo et al., 2003a), and applied it globally at 5 arcmin resolution. Wada et al. (2016) applied the PCR-GLOBWB model globally at 6 arcmin spatial resolution. They estimated sector-wise water usage, along with the supply of surface and groundwater, and then calculated a conventional water stress index as the ratio of usage to supply. They found that water stress is unrealistically concentrated in grid cells containing the downtown areas of the largest cities (e.g., Paris and New York) and therefore recommended assessment across larger spatial domains, such as subbasins and counties. Sutanujaja et al. (2018) developed PCR-GLOBWB2 as an upgrade of PCR-GLOBWB, which has two standard spatial resolutions, i.e., 30 and 5 arcmin. They reported that the latter resolution outperformed the former, mainly due to better representation of snow processes by its more detailed air temperature distribution in mountainous areas. Burek et al. (2020) developed a new global hydrological model called CWatM, and it is based on PCR-GLOBWB and LISFLOOD (De Roo et al., 2000). Although no detailed results were reported, similar to PCR-GLOBWB2, boundary conditions for 5 arcmin resolution have been determined for this model.
Although several successful applications have been described, global hydrological modeling at scales finer than 30 arcmin is still in its infancy. Such modeling efforts face three major challenges. The first challenge is validation. To demonstrate that the models are capable of operating at a certain spatial resolution, reliable observations must be collected in advance. However, aside from river discharge, reliable global observations of hydrological variables are challenging to obtain, especially in the form of fine-scale data for water use and management. The second challenge is model refinement. To express site-specific hydrological processes and water use practices, parameter regionalization (Beck et al., 2016; Yoshida et al., 2021) and model and data localization processes are required, and these processes can easily become unmanageable unless conducted systematically and strategically. The third challenge is the inter-grid-cell connection. Because water supply and drainage systems for urban and irrigation water distribution typically cover a greater extent than the 5 arcmin grid cell, interconnections must be considered.
To investigate the applicability of a global hydrological model at a fine spatial resolution, we applied the latest version of the H08 global hydrological model (Hanasaki et al., 2018b) to Kyushu island, Japan, at 1 arcmin spatial resolution (approximately 2 km at the Equator). The study area provides sufficient quantity and quality of data for simulation and validation with reasonable efforts. The site includes major rivers, dams, irrigation projects, active industries, and densely populated cities, providing good example locations for the application of the hyper-resolution global hydrological model. Similar efforts have been reported by Mateo et al. (2014) and Masood et al. (2015), who applied H08 at 5 arcmin spatial resolution to the Chao Phraya River in Thailand and the Ganges–Brahmaputra–Meghna rivers of the Indian subcontinent, respectively. The present study uses the latest H08 model at a 5 times higher resolution than previous works, with a focus on water management and use. Throughout this study, we assess water resources in terms of whether water is available when and where it is needed. Other assessments, such as flood risk, are outside the scope of this research. The present study addresses the following three questions. How well can the default global H08 reproduce the hydrology and water management of Kyushu island? Does the localization of the data and model improve the performance? What is the benefit of conducting hyper-resolution hydrological modeling from the viewpoint of water resources management and assessment?
Kyushu is the third largest island in Japan, with an area of 36 750 km
In this study, we targeted nine rivers with catchment areas exceeding 1000 km
Study domain. Large and small squares represent manned and automated weather stations, respectively. Numbers in circles indicate gauging stations along major rivers. See Table 1 for a list of rivers and gauging stations. Colored areas are the catchment areas of the gauging stations. Numbers in triangles indicate dams located on the major rivers. See Table 2 for a list of the dams.
Major rivers of Kyushu island. The gauging stations are the official reference points of rivers, where a long-term streamflow record is obtained for the planning of water management and development, with two exceptions. The official reference point for the Chikugo River is Arase (1443 km
Major dams along the rivers included in this study. We excluded the Ryumon Dam due to a lack of historical observation data. We excluded the Matsubara Dam because it is located just below the Shimouke Dam, complicating its use for validation.
In this study, we applied the global hydrological model H08 (Hanasaki et al., 2018b) to Kyushu island in Japan at a spatial resolution of 1 arcmin. We ran the H08 in the following two ways: using the standard global model and using a modified model. H08 requires two types of boundary conditions, namely meteorological information at daily intervals and static geographic information. Boundary conditions were also prepared in two ways. The first set of conditions was obtained using global data. We attempted to find data at 1 arcmin or finer resolution, but most data were available only at 5 arcmin or coarser resolution, which we linearly interpolated into 1 arcmin data. The second method for obtaining boundary conditions employs local data from Japanese governmental agencies and other sources, which are compiled into 1 arcmin resolution. The combination of the default model with global boundary conditions is the current typical implementation of global hyper-resolution simulations, while the localized model with local boundary conditions represents an ideal implementation when data accessibility issues have been resolved (i.e., assuming hydrometeorological data are available at high resolution for any location worldwide).
H08 was originally developed to simulate the global hydrological cycle with major human activities for advanced water resource assessment at the global scale (Hanasaki et al., 2008a, 2018b). The model is grid based, and its standard spatial resolution is 30 arcmin, but it was designed for application to any spatial domain and resolution. The model consists of six submodels, namely land surface hydrology, river routing, crop growth, water abstraction, reservoir operation, and environmental flow. These components enable simulation of water availability and use at daily intervals in a fully consistent manner and division of the water used by three sectors (agricultural, industrial, and municipal) into five surface water sources (river, aqueduct, local reservoir, seawater desalination, and unspecified surface water, which is defined below) and two groundwater sources (renewable and non-renewable groundwater).
The crop growth and water abstraction submodels are detailed below, as they play particularly important roles in this study. The crop growth submodel estimates the daily growth of crops using the heat unit approach and biomass accumulation based on photosynthesis. The formulation of this submodel is identical to the widely used Soil and Water Assessment Tool (SWAT; Arnold et al., 2012). The crop growth submodel has two functions in H08. First, it is a component of H08 that outputs crop-related variables fully consistent with global water cycle simulation, such as crop yield. Second, it is a standalone model for estimating the crop calendar in the study area. Crop calendars, specifically planting and harvesting dates, are essential for determining the timing of irrigation, but local information is difficult to obtain at any scale, as these dates are rarely monitored and recorded. The principal purposes of the standalone simulation were to iterate the crop yield estimation 365 times by shifting the planting date from 1 January to 31 December in all grid cells and, thereby, to identify the date that results in the maximum crop yield.
The water abstraction submodel also plays an essential role in H08. The fundamental concept used for modeling of each water source is explained below. Water abstraction from rivers indicates the removal of a portion of the streamflow to meet the water requirement within a grid cell. Streamflow is estimated by accumulating runoff through the river network. Water abstraction from an aqueduct indicates transfer of streamflow beyond a grid cell through the artificial aqueduct network. This network is defined either explicitly (i.e., manual connections are made between source and destination grid cells) or implicitly (i.e., the origin and destination of aqueducts is inferred by assuming gravitational water transfer via an algorithm; see Hanasaki et al., 2018b, for details). Water abstraction from a local reservoir indicates abstraction of water stored in relatively small reservoirs or irrigation ponds within a grid cell that does not affect the streamflow of the river network. Note that large reservoirs are incorporated into the river network and directly regulate streamflow. Water abstraction from seawater desalination indicates abstraction of water processed in desalination plants, which is possible only in designated regions. Water abstraction from unspecified surface water indicates the volume of the daily water requirement that was from the sources listed above; it indicates that the water requirement was not met (i.e., water deficit) or that errors occurred in simulation (e.g., underestimation of runoff, overestimation of demand, or misrepresentation of water infrastructure). Water abstraction from renewable groundwater indicates water abstraction from the groundwater reservoir of the cell, which can be recharged. Abstraction from non-renewable groundwater indicates deep groundwater pumping after the depletion of the groundwater reservoir. For simplicity, we excluded local reservoirs and seawater desalination from this study. Thus, the water requirement is divided into river, aqueduct, unspecified surface water, renewable groundwater, and non-renewable groundwater.
To reflect basin characteristics more accurately, we localized some submodels of the original H08 model. For the land surface submodel, we tuned the hydrological parameters. Numerous methods have been proposed for the tuning of parameters in hydrological models, and we adopted a simple Monte-Carlo-based method described previously (Mateo et al., 2014; Masood et al., 2015). This process identifies the combination of four sensitive hydrological parameters that minimize the Nash–Sutcliffe efficiency (NSE; Nash and Sutcliffe, 1970) of the daily streamflow simulation among three possible values for each parameter (total
The reservoir operation submodel incorporated into H08 employs simple reservoir operation rules to determine daily water release and storage for individual reservoirs from limited globally available information (Hanasaki et al., 2006). The validity of the original algorithm as a generic global model has been demonstrated in numerous previous works (Masaki et al., 2017; Yassin et al., 2019; Shin et al., 2019; Turner et al., 2020). The original reservoir operation submodel (see Hanasaki et al., 2006) releases water at a constant rate throughout the year at a rate equivalent to the mean annual inflow. This process results in increased storage during the wet season and reduced storage in the dry season. When no other control is applied for subannual operation, storage can reach the storage capacity of the reservoir (and then begin to overflow) or become completely depleted (so that release becomes identical to inflow), which seldom occurs in reality (Masaki et al., 2017).
To enable more realistic behavior, we added the release curve and the upper and lower storage curves to the model (Fig. 2). The release curve assigns different release rates for the wet and dry seasons. The upper and lower storage curves determine the maximum and minimum storage in a year. The storage curves are subdivided into four phases, namely maintaining high storage, releasing, maintaining low storage, and recharging. In the releasing phase, the maximum and minimum storage values are gradually decreased to avoid a sudden increase in water release. In the recharging phase, the maximum and minimum storage values are increased to allow fast and effective recovery. The release curve is also subdivided into four phases, namely the dry season, dry to wet transition, the wet season, and wet to dry transition. To represent the various shapes of release and storage curves, we constructed a generic curve, as shown in Fig. 2. To represent these curves, we specified four points (R1–4; in time and release coordinates) on the release curve, and a total of eight points (U1–4 and L1–4; in time and storage coordinates) on the upper and lower curves. R1–4 were determined from the long-term mean monthly release levels of individual reservoirs. U1–4 and L1–4 were determined from the actual operation curves compiled by reservoir operators. Storage curves of reservoirs managed by central and local governments and public agencies are disclosed due to legislation in Japan. The parameter values used for the six dams are listed in Tables S3–5.
Schematic diagrams of
In reality, to abstract streamflow for water use, water is withdrawn at weirs and transferred directly to the point of use (e.g., irrigated cropland) or to water treatment facilities and then supplied to urban areas through pipelines. Because the identification of weirs and manual delineation of water transfer routes is infeasible for large study domains, H08 can infer the route of inter-grid-cell gravity-based water transfer (Hanasaki et al., 2018b). The original algorithm was designed to work at a resolution of 30 arcmin, and therefore, we modified it for better performance at a fine spatial resolution.
We modified two components of this algorithm. First, the original method allows transfer of only one grid cell ahead but was modified to allow transfers a maximum of five grid cells ahead (approximately 10 km). Next, interbasin transfer was not possible in the original model but was allowed in this study if the catchment area of the destination is smaller than two grid cells. This modification adds flexibility for water transfer, especially in grid cells near coastlines, where basins are generally small (Fig. 3). Along with these implicit aqueducts, the origins and destinations of known major water transfer facilities (designated explicit aqueducts' see Sect. 3.2.) were prepared and implemented into the river network (Fig. S1; Table S6). These explicit aqueducts transfer water irrelevant to elevation.
Implicit aqueducts.
To run the H08 model, the following seven meteorological variables must be supplied at daily intervals: precipitation, air temperature, humidity, wind speed, surface pressure, precipitation, and downward longwave and shortwave radiation. We prepared the data using two methods, including one based on global data representing the current practice in high-resolution global modeling and one using local data.
For the global data, we extracted and linearly interpolated the WFDEI global meteorological dataset (Weedon et al., 2014) for the study domain. The dataset contains all seven variables at daily intervals from 1979–2016. WFDEI is a hybrid product of monthly observation-based global gridded meteorological data, for example, CRU TS3.1 (Climatic Research Unit gridded Time Series v3.1; Harris et al., 2014) and the Global Precipitation Climatology Centre (GPCC; Schneider et al., 2014) and the ERA-Interim daily global reanalysis data (Dee et al., 2011).
For local data, we used published meteorological records. Kyushu island contains 44 observatories (district and local meteorological observatories) and 121 automatic weather stations (Automated Meteorological Data Acquisition System – AMeDAS;
For H08, the boundary conditions listed in Table 3 must be provided. We first prepared all of these parameters by extracting and interpolating their values from the default global simulation or using the functions of H08, which were described in detail by Hanasaki et al. (2018b). Then, we prepared another set of data, utilizing local information whenever possible, as detailed below.
We first developed a flow direction map of Kyushu island. We used the digital elevation model SRTM
Geographic data used in this study.
We adopted the data of Kohyama et al. (2003) for the distributions of crop land, crop type, and irrigated area, which are crucial factors in this simulation. They developed these data by combining National Land Numerical Information (NLNI) provided by the Ministry of Land Infrastructure and Transportation (MLIT) with the census of agriculture and forestry (CAF;
The default H08 model subdivides the grid cells into four land use types, namely double-crop irrigated cropland (for warm regions only), single-crop irrigated cropland, rainfed cropland, and non-cropland (Fig. 4a). The primary and secondary crop types are assumed to be cultivated across the entire grid cell (Hanasaki et al., 2008a). In this study, we first subdivided cropland into two categories, paddy and upland, as irrigation water use differs sharply between these two farming methods. Then, we subdivided the grid cells into five land-use types with specific crop rotations, namely double-crop irrigated paddy (rice and barley), single-crop irrigated paddy (rice), single-crop irrigated upland (summer crop), rainfed upland (summer crop), and non-cropland (Fig. 4b). The crop parameters of the crop growth submodel were taken from SWAT (Arnold et al., 2012). We used the generic values for summer crop.
Grid cell categorization, with
To estimate the crop calendars of three crops, we used the standalone crop growth submodel (see Sect. 3.1.1). In previous global applications of H08, the estimated calendar of major crops agrees well with the reported calendar, which supports the validity of this approach. However, in application to Kyushu, we noted that special treatment is needed for rice. The H08 model estimates the optimal planting date for rice as being early May for most areas of Kyushu. However, the actual cropping practice indicates planting in early June. Therefore, we shifted the planting date of rice by 30 d. For other crops, we used the estimated crop calendar as provided.
Irrigation efficiency is defined as the rate of water consumption (i.e., evapotranspiration) to total water withdrawal. In general, irrigation efficiency is determined by irrigation equipment (i.e., high efficiency for drip irrigation and low efficiency for flood irrigation). Döll and Siebert (2002) adopted this approach and proposed a value of 0.35 for East Asia. Through comparison between simulation results and local observations (reported in Sect. 4.1.2), we obtained a value of 0.45 for surface irrigation in this study, which is by far the predominant irrigation type used on Kyushu island. Last, for the fraction of irrigation obtained from groundwater, we were unable to find reliable statistics that cover the entire island. Therefore, we did not prepare local data and used only the global dataset (Siebert et al., 2010).
Industrial water withdrawal is provided as a census (Ministry of Economy Trade and Industry – METI, 2018). The finest geographical data are provided at the prefecture level and are further divided into 24 industrial classes. As the prefecture level is too crude (with seven in the study domain), we downscaled the data to the municipality level (numbering 233) through weighting by the value of exports, which is available for each municipality and industrial class. Finally, gridded industrial water withdrawal data were developed by assuming that the spatial distribution of water usage is proportional to population (Statistic Bureau of Japan, 2017). All data were prepared for the year 2015 to maintain consistency between municipality-based and grid-based data.
Domestic water withdrawal data were obtained from the Waterworks Statistics (Japan Water Works Association – JWWA, 2018). As the finest resolution available was the prefecture level, gridded data were developed in a manner similar to that of industrial water. The fraction of groundwater withdrawal for industrial and domestic use was obtained from the yearbook (Ministry of Land Infrastructure Transportation and Tourism – MLIT, 2021). These data are available only at the regional scale (Kyushu island is subdivided into two regions) for the industrial and domestic sectors (data for the agriculture sector are not available).
To validate the major outputs, we collected the observational data listed in Table 4. We used the water information system (
Daily irrigation water withdrawal is not published in Japan. We requested daily irrigation water withdrawal data from the Kyushu Regional Agricultural Administration Office (KRAAO) of the Ministry of Agriculture, Forestry, and Fisheries (MAFF). KRAAO kindly provided the daily irrigation water withdrawal data for 18 weirs in seven national irrigation projects on Kyushu, for 11 years, from 2007 to 2018. The data were provided under following conditions: first, the location of the weirs shall not be disclosed. Second, the raw data shall not be published in any document. Third, the redistribution of data shall not be permitted for any reason. Due to these conditions, we hereafter refer to the irrigation projects as A–G.
Validation data. MLIT and MAFF represent the Ministry of Land Infrastructure, and Transportation and the Ministry of Agriculture, Forestry, and Fisheries, respectively.
H08 was set up at 1 arcmin spatial resolution and applied to Kyushu island (Fig. 1). Simulation was conducted for 20 years from 1995 to 2014. Due to the availability of irrigation water withdrawal data, we primarily analyzed the results of the final 8 years (2007 to 2014). Calibration of hydrological parameters was performed using simulation results for 2001–2006. We conducted the following two simulations: a global simulation (hereafter GLB), which involved running the H08 model at a high spatial resolution using the default model and interpolated standard global data and the localized simulation (LOC), which uses localized H08 model with localized boundary conditions (Table 5).
We conducted LOC variant simulations to investigate the effects of the aqueduct settings. Due to substantial differences in water availability (i.e., streamflow) and water demand (affected by irrigated area and population) among grid cells, inter-grid water transfer plays a crucial role in balancing these terms. We conducted simulations with different settings for aqueducts and examined the effect on the simulation results. LOC-I1 employs only the original implicit aqueduct algorithm. Aqueducts transfer water to only one downstream grid cell (approximately 2 km). Because traditional weirs and canal systems can transfer water across several to dozens of kilometers, LOC-I5 extends the maximum reach of implicit aqueducts to five grid cells (9 km). Water-insufficient grid cells tend to occur in coastal zones, where the catchment area is rather limited but heavily populated. LOC-I5O allows water transfer outside of the basin, as illustrated in Sect. 3.1.4. Finally, LOC-I5OE includes explicit aqueducts (i.e., major weirs and municipal water systems; listed in Table S6).
Simulations. Key differences between the GLB and LOC simulations.
Simulated daily river discharge at Senoshita on the Chikugo River is shown in Fig. 5 (results for other stations are shown in Fig. S2). Monthly river discharge obtained from LOC simulations reproduced historical observations well throughout the simulation period. Aspects of the timing and the magnitude of the largest peaks in July and some minor ones, as well as the magnitude of discharge in dry months, were similar to the observations. Daily river discharge also agreed with the observations for both wet and dry years (i.e., the years with minimum and maximum annual rainfall in the validation period), but the simulated peaks tended to be underestimated, especially with respect to minor peaks during the dry period (e.g., from September to May 2012). Monthly river discharge from GLB simulations also agreed with the observations, but the daily results show considerable discrepancies. The magnitude of peaks does not match the observed levels, and the duration of high discharge is apparently longer than the observed period (i.e., the shape of curves is flattened) after a prolonged precipitation event in the GLB simulation (not shown).
River discharge at Senoshita on the Chikugo River. The unit is cubic meters per second (m
The results for other basins showed a similar general tendency (Fig. S2). The results from LOC for both monthly and daily simulations reproduced the observations. GLB performed moderately well for monthly simulations but poorly for daily simulations. The NSE values (Nash and Sutcliffe, 1970) of all seven rivers are shown in Table 6. The NSE of LOC consistently exceeded 0.6, while all values for GLB were below 0.4.
Nash–Sutcliffe efficiency of daily river discharge simulations.
We speculate that the difference in performance between the two simulations is due to the quality of precipitation data. Although the hydrological parameters of LOC were selected from 81 possible combinations, their effect is limited because only one parameter set was applied across the entire study domain (Sect. 3.1.2). The daily variation in precipitation for GLB is based on the ERA-Interim global reanalysis product with a resolution of
approximately
Figure 6 shows the estimated irrigation water demand at irrigation project A (results for the other six sites are shown in Fig. S3). The simulated monthly water requirement of LOC was overestimated by a factor of 2, but the duration of the irrigation and its temporal variations agreed well with observed water withdrawal. GLB had a similar performance to LOC in terms of magnitude, but as discussed above, the irrigation period begins and ends earlier than the observed period by 1 month.
The daily water requirement results for 2007 and 2010 highlight the differences between the two simulations. First, the daily variations in the water requirement are much closer to reality for LOC than for GLB. This difference is due to the meteorological data and hydrological parameters used for each simulation. For example, the smaller
Irrigation water requirement at irrigation project A.
Table 7 shows the monthly NSE for the irrigation water requirements of seven irrigation projects. In general, LOC performed well, except for projects A (Fig. 6) and F. Excluding those two projects, LOC exceeded NSE values of 0.6 and, thus, substantially outperformed GLB, mainly due to the timing of irrigation. The low NSE values in A and F are due to the overestimation of the water requirement. The reasons for such an overestimation have not been fully explored and may be attributed to various possible factors, including the difference between the land use data period (1995) and simulation period (2007–2014) and region-specific irrigation practices. To investigate the influence of settings, a systematic sensitivity simulation for irrigation requirement was conducted (Appendix D).
Monthly NSE of the irrigation requirement.
Figure 7 shows the results for reservoir operation at the Shimouke Dam on the Chikugo River (results for the other five dams are shown in Fig. S4). First, the monthly inflow of LOC agrees well with the observations but shows a tendency toward underestimation of the peaks, which is most apparent in 2009, 2010, 2011, 2012, and 2013. The results of GLB also agree fairly well with observations, but the tendency toward underestimation is more prevalent, including during the dry period.
The monthly release estimates of LOC agree moderately well with observations in terms of the arrival of peaks and seasonal variations, while the estimates of GLB do not. In general, GLB tends to overestimate release in winter (October to March; apparent in 2010–2012), leading to frequent depletion (particularly in the winters of 2007, 2009, and 2013). This overestimation is due to the original algorithm of the reservoir operation submodel of H08, which was designed to release water constantly throughout the year at the rate of mean annual inflow (see Sect. 3.1.3). As numerous previous applications have demonstrated, the algorithm works well for the largest reservoirs in the world, in which total storage capacity is greater than or comparable to mean annual inflow (Hanasaki et al., 2006; Masaki et al., 2017; Shin et al., 2019). However, in this reservoir, the storage capacity is only 14 % (5 %–9 % for the remaining reservoirs assessed here) of the mean annual inflow at the dam and thus is easily depleted during dry periods and overflows in wet periods. The algorithm uses several functions to avoid frequent depletion and overflow events (see Hanasaki et al., 2006), but they are not applicable to this case because the storage capacity is too small to buffer strong temporal variations in inflow.
LOC and GLB contrast sharply in the performance of their monthly storage simulations. The LOC simulation matched the observed amplitude and temporal variations in storage, while GLB results differed markedly; for example, they included extended periods of depletion (e.g., the winters of 2008, 2010, and 2012). As discussed in the previous paragraphs, the hydrographs for inflow and release generally match well. Such a run-of-the river operation is often used for relatively small reservoirs located in rivers with high streamflow. To add the maximum possible flow control to the run-of-the-river operation, incorporation of the upper and lower storage curves (Fig. 2) is beneficial. The seasonal variations in release (Fig. 2) also played an important role in stabilizing storage (i.e., not following the upper or lower curve) and release (i.e., not falling below the intended release level) by reducing the gap between inflow and release.
Reservoir operation at the Shimouke Dam on the Chikugo River.
Table 8 shows NSE values for six reservoirs. In general, the performance of LOC in terms of release and inflow is good (i.e., NSE
NSE values for monthly storage, release, and inflow.
The results presented thus far indicate that the simulation using localized data and model parameters (LOC) exhibited satisfactory performance for the estimation of streamflow in major river basins, irrigation water requirements of major projects, and the operation of major dams. Here, we assess whether water availability and use are balanced. The background for this assessment is that sufficient water is generally available on Kyushu island, except for a few locations and events, and therefore, we expect neither a large spatial extent nor a long duration of water scarcity in the study domain. We quantified water balance by calculating the water sustainability index (WSI), as defined in Eq. (1).
The distribution of mean (2007–2014) WSI is shown in Fig. 8. The results of LOC-I1 (localized model with implicit aqueducts transferring water only 1 grid cell ahead within the same basin) are shown in Fig. 8a. Under this condition, a substantial number of grid cells indicate low to severe water scarcity, including those very close to major rivers (e.g., the lower Chikugo River east of Saga city). Next, the results of LOC-I5 (same as LOC-I1 but with a transfer of up to five grid cells ahead) are shown in Fig. 8b. Although a substantial number of water insufficient grid cells remained, water scarcity was attenuated across the entire domain because water could be transferred a maximum of five grid cells ahead, thereby connecting grid cells with large water requirements to those with large water availability (e.g., northeast of Fukuoka). Then, the results of LOC-I5O (same as LOC-I5 but allowing water transfer outside of basins with small catchment areas) are shown in Fig. 8c. Because water demand is concentrated along the seashore, interbasin water transfer reduced water scarcity in coastal areas dramatically (e.g., southwest of Kumamoto). Finally, the results of LOC-I5OE (same as LOC-I5O but including explicit aqueducts) are shown in Fig. 8d. Under this condition, water scarcity was further attenuated (e.g., the lower reach of the Chikugo River near Saga city). Places with exceptionally low WSI are discussed in the next subsection.
The numbers of grid cells with WSI values below 0.99, 0.80, and 0.50 are shown in Table 9. As the total number of grid cells in the study domain is 13 105, approximately 6 %–8 % of grid cells are categorized as experiencing at least low water scarcity (WSI
Water sustainability index.
Numbers of grid cells below the WSI thresholds.
We investigated locations with medium water scarcity (WSI
The Saga plain contains vast areas of irrigated cropland mainly supplied from the Kasegawa River (catchment area 368 km
Shiroishi plain is located to the southwest of Saga city, along Ariake Sea. A considerable portion of the plain is reclaimed land that is distant from major rivers. Historically, the area has been irrigated using groundwater pumping. After completion of the Kasegawa Dam in 2012, which is located below the Kitagawa Dam, a new surface water irrigation project was initiated to divert water from the Kasegawa River to this area, but that project was
not included in the present study because it was operational only during the last 3 years of the simulation period (2012–2014). Because no detailed data on the spatial distribution of irrigation water sources were available, we the used global data at
The Kitano plain is located southeast of Fukuoka (and northeast of Saga) on the right bank of the middle reach of the Chikugo River. The project in that area is the Ryochiku Heiya Irrigation Project, which is supplied mainly by two tributaries of the Chikugo River, namely the Koishiwaragawa River (catchment area of 85.9 km
Usa city is located along the downstream portion of the Yakkan River (catchment area, 430 km
In this study, we validated the simulated streamflow, irrigation water withdrawal, and reservoir operation for the downstream gauging station of the nine largest rivers, seven national irrigation projects, and six reservoirs on Kyushu island. The results were reasonably good, indicating that the overall hydrological processes and human water use and management were properly simulated at the scales used for validation. Nonetheless, the water supply and demand balance were not fully closed for all grid cells, indicating that further efforts are needed.
The discussion in the previous subsection indicates that cases of possible water scarcity occur in river basins, including reservoirs smaller than the threshold size (i.e.,
In this study, we fixed the surface water and groundwater usage rate to a uniform value across a vast study domain, as the spatial distribution of data was unavailable. As discussed for Shiroishi plain, some regions use groundwater intensively. Another apparent exception to the uniform value is Kumamoto city, where 100 % of the municipal water is supplied by groundwater (Kumamoto City, 2015). Aside from usage, the recharge and flow of groundwater may be additional sources of uncertainty. For example, Mount Aso, which is located at the center of Kyushu island, is an active volcano. The soil in the surrounding area is highly permeable and forms abundant groundwater. The lack of hydrogeological detail in this model is another source of uncertainty.
In this study, we have introduced the implicit and explicit aqueducts which allow water withdrawal from surrounding grid cells where water resources are abundant. This is particularly important for urban areas where municipal and industrial water is highly concentrated. Although improved compared to earlier approaches (e.g., Wada et al., 2016) which assumed water withdrawal and use occurred in a single grid cell, there is still room for progress. Kyushu island contains seven cities, with populations greater than 400 000, and many other settlements. The urban areas of such cities cover dozens of grid cells, which are served by municipal water supply and sewage systems. For urban areas, water is withdrawn at a limited number of inlets to treatment plants and released at outlets from sewage treatment plants. In contrast, H08 allows for water abstraction and drainage from every grid cell, as it was originally developed for a spatial resolution of 30 arcmin, at which the effects of urban water systems become negligible. Thus, water scarcity may be underestimated in major cities in this simulation, which allows for usage of water from one grid cell upstream, or an intensive cascade of drainage among neighboring cells, but completely ignores the presence of municipal water systems. This problem will be largely solved by including waterwork facilities (i.e., water treatment plants, water and sewer pipe networks, and wastewater treatment plants) into the model. This will make the simulation more realistic in terms of representation of the urban water scarcity problem.
While this study has mainly focused on the processes of water use and management, natural hydrological processes also need to be further improved. Snow falls only in limited mountainous areas of Kyushu island, and therefore, the performance of snow process estimation was not evaluated in this study. Groundwater is expressed in highly conceptual way in H08 (treated as a hypothetical tank), which hampers the validation of groundwater simulation. As for the lateral flows of subsurface water, Ji et al. (2017) reported that they matter to the distribution of soil moisture and evapotranspiration at the spatial resolution of 1000 m and finer. The inclusion of the process is crucial, particularly in the regions where groundwater is the major water source. As such, the dominant hydrological processes differ region by region. Further application and investigation of the model under various environment is indispensable to fully realizing globally applicable hyper-resolution modeling.
This study presents a method of validation, localization, and utilization for water resource assessment based on hyper-resolution global hydrological data. The study domain is 37 500 km
River discharge is the easiest factor to validate, but the process has two prerequisites. First, a reliable flow direction map, which reflects the actual shapes of basins and catchment areas of major river gauging points, is needed. Second, a standardized river discharge publication system is required, and must be well designed for (semi-) automatic data retrieval. The water information system used in this study is a good example of such a publication platform but is optimized for manual data retrieval through mouse operation. Also it is only available in Japanese, hindering data retrieval by non-Japanese-speaking users. Such detailed data provisioning systems are commonly available for developed countries but not for other regions. Although efforts have been made to collect historical river discharge records for a considerable number of river basins (van Dijk et al., 2013; Beck et al., 2013), further database development is needed for global hydrological research.
The validation of reservoir operations is more challenging than validation of river discharge. First, reservoir operation records are far less frequently published. This difference is due in part to the fact that reservoirs are operated by multiple agencies. For example, of the 29 reservoirs on Kyushu (i.e.,
Elevation and river topography data are available globally at 3 arcsec or 90 m resolution (Yamazaki et al., 2017, 2019). A central data issue is the upscaling of such fine-scale data to coarser (e.g., 3 to 30 arcsec) resolution, which is not a trivial problem (Yamazaki et al., 2009). In addition, proper transfer of gauging station and infrastructure (e.g., reservoirs, weirs, and canal networks) data is essential. Collection and georeferencing of data for aqueducts and weirs is a challenging issue. Although previous attempts to estimate the water transfer of large cities (McDonald et al., 2014) and mega-projects (Shumilova et al., 2018) have been reported, far more extensive efforts are needed to achieve global coverage, including small projects (neither of the previously cited works covered any of the aqueducts included in this study).
Most global agricultural data are available at 5 arcmin resolution (Ramankutty et al., 2008; Monfreda et al., 2008; Siebert et al., 2005, 2010). As noted above, 5 arcmin resolution is too coarse for use with 1 arcmin or finer resolution. Because land use and crop type significantly affect water demand, their areal information should be obtained from a national agricultural census. Digitizing and georeferencing census results is crucially important (Kohyama et al., 2003). Some factors, including the irrigation water withdrawal, fraction of groundwater use, irrigation efficiency, and crop calendar, are hard to obtain. Systematic collection of data on such everyday agricultural practices is crucial.
Perhaps the most important data type in terms of simulation performance is meteorological data. The overall performance of this simulation depends on good meteorological data, especially high-quality precipitation data. Most global meteorological data currently available have spatial resolution of 30 arcmin (Weedon et al., 2014; Cucchi et al., 2020). Recent efforts have made some of meteorological data available at 15 arcmin resolution (Beck et al., 2017; Hersbach et al., 2020; Yatagai et al., 2012), but a gap in the spatial resolution remains. A community effort with regard to data standardization and sharing is essential for future progress.
In this study, the model was localized in four aspects, namely hydrological parameters, irrigation parameters, reservoir operation settings, and implicit canal estimation configuration. First, hydrological parameter optimization can be more systematically applied. Yoshida et al. (2021) demonstrated the semi-automatic parameter optimization of H08 based on 777 river gauges globally, which is a promising technique for hyper-resolution global application. For irrigation and reservoir parameter settings, the applicability of the methodology devised in this study to other regions remains unknown. Further modification or improvement will likely be needed for application to other parts of Japan and the world, as hydrometeorological conditions differ significantly among regions. Identification of key formulas and parameters to represent various water use practices will continue. Rigorous databasing of large-scale aqueducts and further algorithm validation would make the implicit canal inference process applicable to a broad range of regions. Last, the validity of implicit canal estimation was demonstrated only indirectly in this study. The results indicate great improvement of the balance between water supply and demand in the study domain, with realistic results for Kyushu island, but direct validation against a more detailed database of weirs and canals is required for generalizing the origins and destinations of water diversion networks in Japan.
High-resolution water assessment will provide novel insights into issues related to water resources. For example, the water scarcity distribution shown in Fig. 8 is a reminder of how sensitive water availability is to climatic and geographic conditions. Water scarcity may occur around tributaries of major rivers (e.g., the Kitano plain), although such areas have been largely overlooked in low-resolution modeling efforts. The present study also highlights the importance of aqueducts and reservoirs, which play crucial roles in balancing water demand and supply. Ignoring water infrastructure may result in misleading water resource assessment results.
In this study, we applied the global hydrological model H08 to Kyushu island, at a spatial resolution of 1 arcmin, to address three research questions. First, in terms of the reproducibility of the default H08 model and simulation settings, the detailed validation results indicate highly negative values, suggesting that a simple increase in the spatial resolution of global simulation does not lead to reliable results at the local scale.
The results highlight the importance of rigorous validation of the most important variables and locations. Second, the process of model localization was demonstrated as promising and effective. Localization of the data and model improve the performance to a satisfactory level (e.g., daily streamflow NSE
This study contributes to hydrological science in three ways. First, this is the first study to provide an example application of a global water resources model to 1 arcmin spatial resolution. Although the study domain was rather limited (
For practical climate change adaptation planning and water resource development, local human water use and management are non-negligible problems that can be addressed with sophisticated hyper-resolution global modeling. The detailed information obtained from this study can be utilized to support the recent trend toward identifying environmental impacts of corporate activities in a meaningful manner (e.g., CDP water security;
We converted sunshine duration (h;
We estimated downward longwave radiation (W m
We have decomposed the differences between GLB and LOC simulations into the following three factors: (A) calibration, (B) local meteorological observation, and (C) high spatial resolution of data. The combination of factors is shown in Table C1.
YYY and NNN are identical to the LOC and the GLB simulations, respectively. YYY and NYY use the high-resolution daily local meteorological observation (AMeDAS). YNN and NNN use the global meteorological data (WFDEI; Weedon et al. 2014). YYN and NYN use the low-resolution AMeDAS, which was prepared as follows: we aggregated the AMeDAS data for
The results are shown in Fig. C1. In general, the calibration increases the NSE for many cases if other conditions are identical (21 out of 27), but the effect (the change in NSE) is limited. Exceptions are seen, for example, in Stations 1, 6, and 8, where NYY (without calibration) performs slightly better than YYY (with calibration). This is considered that due to possible overfitting during the calibration period. Next, high spatial resolution increases the NSE for most of the cases (16 out of 18) but again the effect is limited. Exceptions are seen in Stations 2 and 5, where NYN (low resolution) outperforms NYY (high resolution). Finally, local meteorological data significantly increases the NSE for all the cases (18 out of 18 by the median of 0.43).
The results of sensitivity test clearly indicate that the usage of local meteorological observation is the key factor in the simulation performance. As mentioned earlier, the timing and magnitude of submonthly weather events in the global meteorological data (WFDEI) are based on reanalysis which are not always successful in reproducing the timing and the location where weather events occurred.
Sensitivity simulation. The simulation names and the combinations of factors are given.
The results of sensitivity analysis on river discharge. See Table C1 for simulation runs. YYY and NNN correspond to the GLB and LOC simulations, respectively.
The reproducibility of reanalysis is being improved generation by generation. By using the latest global meteorological data based on the latest generation of reanalysis, the GLB simulation is expected to perform better (e.g., the WFDE5, Cucchi et al., 2020, which uses the ERA5 reanalysis; Hersbach et al., 2020).
We have decomposed the differences between GLB and LOC simulations into the following three factors: adoption of (A) local crop calendar, (B) local cropping practices, and (C) local hydrometeorological conditions. The combination of factors is shown in Table D1.
The sensitivity simulations SI1, SI2, SI3 were identical to LOC, except for the aforementioned three factors. SI1 replaced the adjusted crop calendar for rice with the originally simulated one (Sect. 3.2.3). SI2 is same as SI1 but additionally replaced the localized settings of the land use (Sect. 3.2.3, Fig. 4) and the irrigation efficiency parameter (Sect. 3.2.3) with those of GLB. Here, crop calendar, land use, and irrigation efficiency are all essential in estimating irrigation water requirement but virtually solely available from governmental reports and, hence, hard to obtain globally. We singled out the crop calendar (i.e., separated SI1 and SI2) because it would be available in the future globally since some concrete methods have been proposed by utilizing the latest satellite remote sensing techniques (Kotsuki and Tanaka, 2015). SI3 is same as SI2 but additionally replaced the local meteorological data with the global ones and the calibrated hydrological parameters with the default ones (Sect. 3.1.2 and 3.2.1).
The results are shown in Fig. D1. In general, the NSE scores for SI1, SI2, and SI3 are deteriorated in this order and approaching those for GLB. The apparent exception is the influence of local cropping practices (i.e., SI2). The irrigation efficiency of GLB is set at 0.35 uniformly, while LOC is at 0.45–0.8. The adoption of low efficiency inflates irrigation water requirement which drastically decreases the NSE score. The effect is canceled out when the global hydrometeorological conditions are applied, which suppresses the total irrigation water requirement (Sect. 4.1.2).
In contrast to the sensitivity test for discharge shown in Appendix C, the influence of adopting local meteorological condition looks marginal. This is, however, largely due to the considerable drop in the NSE score caused by other factors. This sensitivity simulation highlights the importance of reliable local cropping practices. They can alter the shape, the phase, and the magnitude of simulation time series and are, hence, highly influential in the results and performances.
Sensitivity simulation on irrigation requirement.
The results of sensitivity analysis. See Table D1 for simulation runs. The NSE for irrigation project F (SI2) is
We have decomposed the differences between GLB and LOC simulations into the following three factors: adoption of (A) upper and lower storage rule curves (Fig. 2a), (B) release curve (Fig. 2b), and (C) local hydrometeorological conditions. The combination of factors is shown in Table E1.
The sensitivity simulations SD1, SD2, and SD3 were identical to LOC except for the three factors. SD1 removed the setting of upper and lower storage curves (Sect. 3.1.3, Fig. 2). This setting allows the storage to fluctuate, at any time, between being completely depleted and completely full. SD2 is same as SD1 but, additionally, removed the release curve. This setting fixes the release throughout a year at the mean annual inflow (see Hanasaki et al., 2006, for the complete formulations). Although these curves are essential in the actual reservoir operation, it is hard to obtain globally. The advances in satellite global reservoir monitoring would offer the opportunity to estimate the upper and lower storage curves for major reservoirs (Busker et al., 2019). Hence, we considered that the chance of earning the storage curves (SD1) is more likely than the release one (SD2). SD3 is same as SD2 but additionally replaced the local meteorological data with the global data and the calibrated hydrological parameters with the default parameters (Sect. 3.1.2 and 3.2.1).
Sensitivity simulation on reservoir operation.
Sensitivity analysis for reservoir storage. See Table E1 for the simulation runs.
The results for storage are shown in Fig. E1. In general, the NSE scores for SD1, SD2, and SD3 are deteriorated in this order and approaching that of GLB. Note that SD3 and GLB do not completely agree because of other settings between LOC and GLB (e.g., water withdrawal and diversion). The removal of storage curves has the largest impact on the performance of storage simulation, which can be easily expected. The storage capacity of the dams is relatively small compared to their inflow, which typically causes strong storage variations. Once the constraints due to storage rule curves are lost, the simulated storage showed considerable variations by reflecting errors in inflow and release.
The results for release are shown in Fig. E2. Again, the NSE scores for SD1, SD2, and SD3 are deteriorated in this order and approaching that of GLB. The impacts due to the adoption of the storage curves, the release curve, and the local hydrometeorological data are not easily discernible, except in the case for the Midorikawa Dam and the Tsuruda Dam. For these two reservoirs, the adoption of local hydrometeorological data dominantly influenced the simulation performance.
Overall, specifying the key factors in the performance of reservoir operation simulation is harder than the discharge (Appendix C) and irrigation requirement (Appendix D). This is likely due to the overall low reproducibility of the reality and, in particular, the storage variation. Further fundamental improvement is needed in reservoir modeling to reproduce the operations of the past.
Sensitivity analysis for reservoir release. See Table E1 for the simulation runs.
Code and data are available at
The supplement related to this article is available online at:
NH designed the experiments, and NH, HM, and MF carried them out. NH, HM, and MF developed the model code and performed the simulations. YH, SS, SK, and TO validated the model. NH prepared the paper, with contributions from all co-authors.
The contact author has declared that neither they nor their co-authors have any competing interests.
Publisher’s note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The authors are grateful to Noemi Vergopolan, the two anonymous referees, and Luca Brocca, for providing comments on this work. This research has been supported by KAKENHI (grant nos. 16H06291, 21H05178, and 21H05002) of the Japan Society for the Promotion of Science. The authors are responsible for any errors.
This research has been supported by the Japan Society for the Promotion of Science (KAKENHI; grant nos. 16H06291, 21H05178, and 21H05002).
This paper was edited by Xing Yuan and reviewed by Noemi Vergopolan and two anonymous referees.