**Effect of Evapotranspiration on Hydrothermal Changes in Regional Scale**

Tadanobu Nakayama

Additional information is available at the end of the chapter

http://dx.doi.org/10.5772/52808

### **1. Introduction**

Evapotranspiration plays an important role not only on hydrologic cycle but also on thermal changes in various ways. In China, hydro-climate is diverse between north and south (Fig. 1). Semi-arid north is heavily irrigated and combination of increased food demand and declining water availability is creating substantial pressures in Yellow River (Brown and Halweil, 1998; Yang et al., 2004; Nakayama et al., 2006, 2010; Nakayama, 2011a, 2011b), whereas flood storage ability around lakes has decreased and impact of Three Gorges Dam (TGD) on flood occurrence in Changjiang downstream against original purpose is increasing problem in humid south (Shankman and Liang, 2003; Zhao et al., 2005; Nakayama and Watanabe, 2008b). Irrigation has a different impact on evapotranspiration changes at rotation between winter wheat and summer maize in the semi-arid region in the north (downstream of Yellow River), and double-cropping of rice in the humid south (middle of Changjiang River) in China. This mechanism changes greatly hydrologic cycle such as river discharge and groundwater, and in particular, affects extremes of flood and drought under climatic change (Nakayama, 2011a, 2011b, 2012c; Nakayama and Watanabe, 2006, 2008b; Nakayama et al., 2006, 2010).

On the other hand, urban heat island (Oke, 1987), where the urban temperature is higher than its rural surroundings, has become a serious environmental problem with the expansion of cities and industrial areas in the world (Fig. 1). Surfaces covered by concrete or asphalt can absorb a large amount of heat during the day and release it to the atmosphere at night. The evaporation of water provides an important counter to this effect, and so open parks and water surfaces are vital in urban areas for creating urban cool-island (Spronken-Smith and Oke, 1999; Chang et al., 2007). Recent researches showed that cooling effect of water-holding pavements made of new symbiotic material (consisting of porous asphalt and water-holding filler made of steel by-products based on silica compound) in addition to that of natural green

© 2013 Nakayama; licensee InTech. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. © 2013 Nakayama; licensee InTech. This is a paper distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

area on hydrothermal cycle is effective to recover sound hydrologic cycle and to create thermally-pleasing environments in eco-conscious society (Nakayama and Fujita, 2010; Nakayama and Hashimoto, 2011; Nakayama et al., 2007, 2012).

variations of LAI (leaf area index) and FPAR (fraction of photosynthetically active radiation) from satellite data (Fig. 2) (Nakayama, 2008a, 2008b, 2009, 2010, 2011a, 2011b, 2012a, 2012b, 2012c; Nakayama and Fujita, 2010; Nakayama and Hashimoto, 2011; Nakayama and Wata‐ nabe, 2004, 2006, 2008a, 2008b, 2008c; Nakayama et al., 2006, 2007, 2010, 2012). The unsaturated layer divides canopy into two layers, and soil into three layers in the vertical dimension in the SiB2 (Simple Biosphere model 2) (Sellers et al., 1996). About the saturated layer, the NICE solves three-dimensional groundwater flow for both unconfined and confined aquifers. The hillslope hydrology can be expressed by the two-layer surface runoff model including freezing/ thawing processes. The NICE connects each sub-model by considering water/heat fluxes: gradient of hydraulic potentials between the deepest unsaturated layer and the groundwater,

Effect of Evapotranspiration on Hydrothermal Changes in Regional Scale

http://dx.doi.org/10.5772/52808

129

In agricultural field, NICE is coupled with DSSAT (Decision Support Systems for Agrotechnology Transfer) (Ritchie et al., 1998), in which automatic irrigation mode supplies crop water requirement, assuming that average available water in the top layer falls below soil moisture at field capacity for cultivated fields (Nakayama et al., 2006). The model includes different functions of representative crops (wheat, maize, soybean, and rice) and simulates automatically dynamic growth processes. Potential evaporation is calculated on Priestley and Taylor equation (Priestley and Taylor, 1972), and plant growth is based on biomass formula‐ tion, which is limited by various reduction factors like light, temperature, water, and nutrient,

In urban area, NICE is coupled with UCM (Urban Canopy Model) to include the effect of hydrothermal cycle at various pavements, and with RAMS (Regional Atmospheric Modeling System) (Pielke et al., 1992) to include the hydrothermal interaction (Nakayama and Fujita, 2010; Nakayama and Hashimoto, 2011; Nakayama et al., 2012). In particular, the author expanded specific heat conductivity *cs* and heat conductivity *ks* in natural soil (Sellers et al., 1996) for engineered pavement in the following equations by including the effect of water

1

*W*

(1)

(2)

*C C W for engineered pavement*

 r

et al. (Nakayama et al., 2006; Nakayama and Watanabe, 2008b; Nakayama, 2011a).

effective precipitation, and seepage between river and groundwater.

amount on the heat characteristics in the material.

( )

1

0.5 1

(=4.18×106

); *cs* (J/m3

r

( )

( )

q

 q

*s s s ww p p s w ws*

qq

= - +× × é ù ë û = × - + × ××

1

( )

<sup>1</sup> <sup>1</sup>

<sup>=</sup> + -×

*p s ws*

= - + ××

1.5 1 1.3 <sup>1</sup> 0.4186

q

 q

0.75 0.65 0.4 <sup>1</sup>


*s s k k W for engineered pavement*

 q

> q

*Cp* (J/kg/K) is specific heat of pavement material; *Cw* (J/kg/K) is specific heat of water

heat conductivity including the effect of water; *kp* (W/m/K) is heat conductivity of pavement

/K) is specific heat conductivity including the effect of water; *ks* (W/m/K) is

*s s <sup>k</sup> for soil <sup>s</sup> <sup>W</sup>*

q

 rq

*c W C for soil*

In this way, the evapotranspiration plays an important role on hydrologic change in conti‐ nental basins where water resources are vital for human activity, and effective management of water resource is powerful for decision-making and adaptation strategy for sustainable development. This chapter represents the improvement in process-based National Integrated Catchment-based Eco-hydrology (NICE) model series (Nakayama, 2008a, 2008b, 2009, 2010, 2011a, 2011b, 2012a, 2012b, 2012c; Nakayama and Fujita, 2010; Nakayama and Hashimoto, 2011; Nakayama and Watanabe, 2004, 2006, 2008a, 2008b, 2008c; Nakayama et al., 2006, 2007, 2010, 2012) with more complex sub-systems to develop coupled human and natural systems and to analyze impact of evapotranspiration on hydrothermal changes in regional scale.

**Figure 1.** Study area in East Asia including the Changjiang and Yellow River basins in China, and the Tokyo Metropolis in Japan.

### **2. Material and methods**

### **2.1. Coupling of process-based model with complex sub-systems**

Previously, the author developed the process-based NICE model, which includes surfaceunsaturated-saturated water processes and assimilates land-surface processes describing the variations of LAI (leaf area index) and FPAR (fraction of photosynthetically active radiation) from satellite data (Fig. 2) (Nakayama, 2008a, 2008b, 2009, 2010, 2011a, 2011b, 2012a, 2012b, 2012c; Nakayama and Fujita, 2010; Nakayama and Hashimoto, 2011; Nakayama and Wata‐ nabe, 2004, 2006, 2008a, 2008b, 2008c; Nakayama et al., 2006, 2007, 2010, 2012). The unsaturated layer divides canopy into two layers, and soil into three layers in the vertical dimension in the SiB2 (Simple Biosphere model 2) (Sellers et al., 1996). About the saturated layer, the NICE solves three-dimensional groundwater flow for both unconfined and confined aquifers. The hillslope hydrology can be expressed by the two-layer surface runoff model including freezing/ thawing processes. The NICE connects each sub-model by considering water/heat fluxes: gradient of hydraulic potentials between the deepest unsaturated layer and the groundwater, effective precipitation, and seepage between river and groundwater.

area on hydrothermal cycle is effective to recover sound hydrologic cycle and to create thermally-pleasing environments in eco-conscious society (Nakayama and Fujita, 2010;

In this way, the evapotranspiration plays an important role on hydrologic change in conti‐ nental basins where water resources are vital for human activity, and effective management of water resource is powerful for decision-making and adaptation strategy for sustainable development. This chapter represents the improvement in process-based National Integrated Catchment-based Eco-hydrology (NICE) model series (Nakayama, 2008a, 2008b, 2009, 2010, 2011a, 2011b, 2012a, 2012b, 2012c; Nakayama and Fujita, 2010; Nakayama and Hashimoto, 2011; Nakayama and Watanabe, 2004, 2006, 2008a, 2008b, 2008c; Nakayama et al., 2006, 2007, 2010, 2012) with more complex sub-systems to develop coupled human and natural systems and to analyze impact of evapotranspiration on hydrothermal changes in regional scale.

> Impact of Megalopolis

> > Urban Area

Yellow River Japan

**Figure 1.** Study area in East Asia including the Changjiang and Yellow River basins in China, and the Tokyo Metropolis

Previously, the author developed the process-based NICE model, which includes surfaceunsaturated-saturated water processes and assimilates land-surface processes describing the

Tokyo Metropolis

Regional Area

Nakayama and Hashimoto, 2011; Nakayama et al., 2007, 2012).

North China Plain

China

Changjiang River

**2.1. Coupling of process-based model with complex sub-systems**

Drought

128 Evapotranspiration - An Overview

Flood

**2. Material and methods**

in Japan.

In agricultural field, NICE is coupled with DSSAT (Decision Support Systems for Agrotechnology Transfer) (Ritchie et al., 1998), in which automatic irrigation mode supplies crop water requirement, assuming that average available water in the top layer falls below soil moisture at field capacity for cultivated fields (Nakayama et al., 2006). The model includes different functions of representative crops (wheat, maize, soybean, and rice) and simulates automatically dynamic growth processes. Potential evaporation is calculated on Priestley and Taylor equation (Priestley and Taylor, 1972), and plant growth is based on biomass formula‐ tion, which is limited by various reduction factors like light, temperature, water, and nutrient, et al. (Nakayama et al., 2006; Nakayama and Watanabe, 2008b; Nakayama, 2011a).

In urban area, NICE is coupled with UCM (Urban Canopy Model) to include the effect of hydrothermal cycle at various pavements, and with RAMS (Regional Atmospheric Modeling System) (Pielke et al., 1992) to include the hydrothermal interaction (Nakayama and Fujita, 2010; Nakayama and Hashimoto, 2011; Nakayama et al., 2012). In particular, the author expanded specific heat conductivity *cs* and heat conductivity *ks* in natural soil (Sellers et al., 1996) for engineered pavement in the following equations by including the effect of water amount on the heat characteristics in the material.

$$\begin{aligned} \mathcal{C}\_s &= \left[ 0.5 \left( 1 - \theta\_s \right) + \theta\_s \cdot \mathbf{W}\_1 \right] \mathbf{C}\_w \cdot \boldsymbol{\rho}\_w \quad \text{for soil} \\ &= \mathbf{C}\_p \cdot \boldsymbol{\rho}\_p \left( 1 - \theta\_s \right) + \mathbf{C}\_w \cdot \boldsymbol{\rho}\_w \cdot \theta\_s \cdot \mathbf{W}\_1 \quad \text{for engineered parameter} \end{aligned} \tag{1}$$

$$\begin{aligned} k\_s &= 0.4186 \frac{1.5 \left(1 - \theta\_s\right) + 1.3 \theta\_s \cdot W\_1}{0.75 + 0.65 \theta\_s - 0.4 \theta\_s \cdot W\_1} \quad \text{for soil} \\ &= k\_p \left(1 - \theta\_s\right) + k\_w \cdot \theta\_s \cdot W\_1 \quad \text{for engineered parameter} \end{aligned} \tag{2}$$

*Cp* (J/kg/K) is specific heat of pavement material; *Cw* (J/kg/K) is specific heat of water (=4.18×106 ); *cs* (J/m3 /K) is specific heat conductivity including the effect of water; *ks* (W/m/K) is heat conductivity including the effect of water; *kp* (W/m/K) is heat conductivity of pavement material; *kw* (W/m/K) is heat conductivity of water (=0.59); *ρ<sup>p</sup>* (kg/m3 ) is specific gravity of pavement material; *ρ<sup>w</sup>* (kg/m3 ) is density of water (=1,000), *Wi* is the soil moisture fraction of the *i*-th layer (=*θ<sup>i</sup>* /*θs*); *θs* (m3 /m3 ) is volumetric soil moisture in the *i*-th layer; *θs* (m3 /m3 ) is the value of *θ* at saturation, respectively.

ma and Watanabe, 2006, 2008b; Nakayama et al., 2006, 2010). For a multi-scaled model in the Tokyo area, hourly observation data from AMeDAS (Automated Meteorological Data Acquisition System) data (Japan Meteorological Agency, 2005-2006) were assimilated with the model. At the lateral boundaries of regional area, some meteorological data are input to the model from the ECMWF (European Centre for Medium-Range Weather Forecasts) with a resolution of 1°x 1° and from the MSM (Meso Scale Model) with a resolution of 10km x 10km (Japan Meteorological Agency, 2005-2006). Mean elevation was calculated by using a global digital elevation model (DEM; GTOPO30) (U.S. Geological Survey, 1996). Digital land cover data were categorized such as forests, grasses, bushes, shrubs, paddy fields, and cultivated fields. About 50 vegetation and soil parameters were calculated on the basis of vegetation class and soil maps (Chinese Academy of Sciences, 1988; Digital National Land Information GIS data of Japan, 2002). The geological structures were divided into four types on the basis of hydraulic conductivity, the specific storage of porous material, and specific yield by scanning and digitizing the geological material (Geological Atlas of China, 2002; Nakayama et al., 2007) and core-sampling data. Artificial augmentation of waterworks and sewerage systems, and anthropogenic sensible, latent, and sewage heats generated by buildings and factories were input into the model (Nakayama and Fujita, 2010; Nakayama and Hashimoto, 2011;

Effect of Evapotranspiration on Hydrothermal Changes in Regional Scale

http://dx.doi.org/10.5772/52808

131

At the upstream boundaries, a reflecting condition on the hydraulic head was used assuming that there is no inflow from the mountains in the opposite direction (Nakayama and Watanabe, 2004). Time-series of tidal level was input as a variable head at the sea boundary (Nakayama, 2011a; Nakayama et al., 2006, 2007, 2010). Vertical geological structures were divided into 10– 20 layers by using sample database. The hydraulic head values parallel to the ground level were input as the initial conditions. In river grids decided by digital river network, inflows or outflows from the riverbeds were simulated at each time step depending on the difference in the hydraulic heads of groundwater and river. The hydraulic head values parallel to the observed ground level were input as initial conditions for the groundwater sub-model. In river grids decided by digital river network from topographic maps, inflows or outflows from the riverbeds were simulated at each time step depending on the difference in the hydraulic heads

The simulation area covered 3,000 km by 2,000 km with a grid spacing of 10 km, covering the entire Changjiang and Yellow River Basins (Fig. 1). The vertical layer was discretized in thickness with depth, with each layer increased in thickness by a factor of 1.1 (Nakayama, 2011b; Nakayama and Watanabe, 2008b; Nakayama et al., 2006). In the Tokyo area, the simulation was conducted with multi-scaled levels in horizontally regional area (260 km wide by 260 km long with a grid spacing of 2 km covering Kanto region) nesting with one way to urban area (36 km wide by 26 km long with 200 m grid covering the Kawasaki City) (Nakayama and Hashimoto, 2011; Nakayama et al., 2007, 2012) (Fig. 1). These areas were discretized with a grid spacing of 200 m – 5 km in the horizontal direction. The NICE simulation was conducted on a NEC SX-8 supercomputer. The first 6 months were used as a warm-up period until equilibrium water levels were reached, and parameters were estimated by a comparison of simulated steady-state value with that published in previous literatures. A time step of the

Nakayama et al., 2007, 2012).

of groundwater and river.

Downward short- and long-wave radiation, precipitation, atmospheric pressure, air temper‐ ature, air humidity, and wind speed simulated by the atmospheric model are input into the UCM, whereas momentum, sensible, latent, and long-wave flux simulated by the UCM are input into the atmospheric model at each time step. This procedure means that the feedback process about water and heat transfers between atmospheric region and land surface are implicitly included in the simulation process.

**Figure 2.** National Integrated Catchment-based Eco-hydrology (NICE) model.

#### **2.2. Model input data and running the simulation**

Six-hour re-analyzed data were input into the model after interpolation of ECMWF (European Centre for Medium-Range Weather Forecasts) in inverse proportion to the distance backcalculated in each grid. Because the ECMWF precipitation had the least reliability and underestimated observed peak values, rain gauge daily precipitation collected at meteorolog‐ ical stations were used to correct the ECMWF value (Nakayama, 2011a, 2011b, 2012c; Nakaya‐ ma and Watanabe, 2006, 2008b; Nakayama et al., 2006, 2010). For a multi-scaled model in the Tokyo area, hourly observation data from AMeDAS (Automated Meteorological Data Acquisition System) data (Japan Meteorological Agency, 2005-2006) were assimilated with the model. At the lateral boundaries of regional area, some meteorological data are input to the model from the ECMWF (European Centre for Medium-Range Weather Forecasts) with a resolution of 1°x 1° and from the MSM (Meso Scale Model) with a resolution of 10km x 10km (Japan Meteorological Agency, 2005-2006). Mean elevation was calculated by using a global digital elevation model (DEM; GTOPO30) (U.S. Geological Survey, 1996). Digital land cover data were categorized such as forests, grasses, bushes, shrubs, paddy fields, and cultivated fields. About 50 vegetation and soil parameters were calculated on the basis of vegetation class and soil maps (Chinese Academy of Sciences, 1988; Digital National Land Information GIS data of Japan, 2002). The geological structures were divided into four types on the basis of hydraulic conductivity, the specific storage of porous material, and specific yield by scanning and digitizing the geological material (Geological Atlas of China, 2002; Nakayama et al., 2007) and core-sampling data. Artificial augmentation of waterworks and sewerage systems, and anthropogenic sensible, latent, and sewage heats generated by buildings and factories were input into the model (Nakayama and Fujita, 2010; Nakayama and Hashimoto, 2011; Nakayama et al., 2007, 2012).

material; *kw* (W/m/K) is heat conductivity of water (=0.59); *ρ<sup>p</sup>* (kg/m3

/m3

**Figure 2.** National Integrated Catchment-based Eco-hydrology (NICE) model.

Six-hour re-analyzed data were input into the model after interpolation of ECMWF (European Centre for Medium-Range Weather Forecasts) in inverse proportion to the distance backcalculated in each grid. Because the ECMWF precipitation had the least reliability and underestimated observed peak values, rain gauge daily precipitation collected at meteorolog‐ ical stations were used to correct the ECMWF value (Nakayama, 2011a, 2011b, 2012c; Nakaya‐

**2.2. Model input data and running the simulation**

) is density of water (=1,000), *Wi*

Downward short- and long-wave radiation, precipitation, atmospheric pressure, air temper‐ ature, air humidity, and wind speed simulated by the atmospheric model are input into the UCM, whereas momentum, sensible, latent, and long-wave flux simulated by the UCM are input into the atmospheric model at each time step. This procedure means that the feedback process about water and heat transfers between atmospheric region and land surface are

) is volumetric soil moisture in the *i*-th layer; *θs* (m3

pavement material; *ρ<sup>w</sup>* (kg/m3

/*θs*); *θs* (m3

implicitly included in the simulation process.

value of *θ* at saturation, respectively.

the *i*-th layer (=*θ<sup>i</sup>*

130 Evapotranspiration - An Overview

) is specific gravity of

/m3

) is the

is the soil moisture fraction of

At the upstream boundaries, a reflecting condition on the hydraulic head was used assuming that there is no inflow from the mountains in the opposite direction (Nakayama and Watanabe, 2004). Time-series of tidal level was input as a variable head at the sea boundary (Nakayama, 2011a; Nakayama et al., 2006, 2007, 2010). Vertical geological structures were divided into 10– 20 layers by using sample database. The hydraulic head values parallel to the ground level were input as the initial conditions. In river grids decided by digital river network, inflows or outflows from the riverbeds were simulated at each time step depending on the difference in the hydraulic heads of groundwater and river. The hydraulic head values parallel to the observed ground level were input as initial conditions for the groundwater sub-model. In river grids decided by digital river network from topographic maps, inflows or outflows from the riverbeds were simulated at each time step depending on the difference in the hydraulic heads of groundwater and river.

The simulation area covered 3,000 km by 2,000 km with a grid spacing of 10 km, covering the entire Changjiang and Yellow River Basins (Fig. 1). The vertical layer was discretized in thickness with depth, with each layer increased in thickness by a factor of 1.1 (Nakayama, 2011b; Nakayama and Watanabe, 2008b; Nakayama et al., 2006). In the Tokyo area, the simulation was conducted with multi-scaled levels in horizontally regional area (260 km wide by 260 km long with a grid spacing of 2 km covering Kanto region) nesting with one way to urban area (36 km wide by 26 km long with 200 m grid covering the Kawasaki City) (Nakayama and Hashimoto, 2011; Nakayama et al., 2007, 2012) (Fig. 1). These areas were discretized with a grid spacing of 200 m – 5 km in the horizontal direction. The NICE simulation was conducted on a NEC SX-8 supercomputer. The first 6 months were used as a warm-up period until equilibrium water levels were reached, and parameters were estimated by a comparison of simulated steady-state value with that published in previous literatures. A time step of the simulation was changed from Δt = 1.5 sec to 1 h depending on spatial scale and the sub-model. Simulations were validated against various hydrothermal observed variables such as river discharge, soil moisture, groundwater level, air temperature, surface temperature, and heatflux budget, et al.

### **3. Result and discussion**

#### **3.1. Effect of irrigation on hydrologic change**

After the verification procedure (Nakayama, 2011b; Nakayama and Watanabe, 2008b), the model simulated effect of irrigation on evapotranspiration at rotation between winter wheat and summer maize in the downstream of Yellow River, and double-cropping of rice in the middle of Changjiang River (Fig. 3). Because more water is withdrawn during winter-wheat period due to small rainfall in the semi-arid north, the irrigation in this period affects greatly the increase in evapotranspiration (Fig. 3a), which is supplied by the limited water resources of river discharge and groundwater there. In particular, most of the irrigation is withdrawn from aquifer in the North China Plain (NCP) because surface water is seriously limited there (Nakayama, 2011a, 2011b; Nakayama et al., 2006). In the south, the irrigation is usually from the river to fill the paddy fiddle as ponding water depth (Nakayama, 2012c; Nakayama and Watanabe, 2008b), which increases evapotranspiration more in the drier season (Fig. 3b). This implies that energy supply is abundant relative to the water supply and the hydrological process is more sensitive to precipitation in the north, whereas the water supply is abundant relative to the energy supply and sun duration has a more significant impact in the south (Cong et al., 2010).

The model also simulated groundwater level in both Changjiang and Yellow River basins (Fig. 3c). The level decreases rapidly around the source area and the Qinghai Tibet Plateau, indicating that there are many sources of spring water in this region. The value is very low in the downstream because of the low elevation and overexploitation, in particular, in the NCP (Nakayama, 2011a, 2011b; Nakayama et al., 2006). This result indicates that hydrologic cycle including groundwater level is highly related not only to the topography but also to the irrigation water use. The NICE is effective to provide better evaluation of hydrological trends in longer period including 'evaporation paradox' (Roderick and Farquhar, 2002; Cong et al., 2010) together with observation networks because the model does not need the crop coefficient (depending on a growing stage and a kind of crop) for the calculation of actual evaporation and simulates it directly without detailed site-specific information or empirical relation to calculate effective precipitation (Nakayama, 2011a; Nakayama et al., 2006).

The mean TINDVI (Time-Integrated Normalized Difference Vegetation Index) gradients during 1982-1999 in various field crops (wheat, maize, and rice) at 4 stations were compared with trends of crop yields in the previous research (Tao et al., 2006) (Fig. 4a). The correlation of both values is relatively good (*r*<sup>2</sup> = 0.986) and the TINDVI gradient has a linear relation to the yield trend. The spatial pattern of the mean TINDVI gradient in agricultural fields shows a generally increasing tendency, particularly in the Yellow downstream and the NCP (Fig.

Fig. 3. Impact of irrigation and ponding water depth on hydrological cycle; (a) evapotranspiration change at rotation between winter wheat and summer maize in the downstream of Yellow River, (b) evapotranspiration change at doublecropping of rice in the middle of Changjiang River, and (c) simulated results of annual-averaged groundwater level (a.s.l.) in the Changjiang and Yellow River basins. In Fig. 3a and 3b, right axis shows a period of each crop (WH; wheat,

**Figure 3.** Impact of irrigation and ponding water depth on hydrological cycle; (a) evapotranspiration change at rota‐ tion between winter wheat and summer maize in the downstream of Yellow River, (b) evapotranspiration change at double-cropping of rice in the middle of Changjiang River, and (c) simulated results of annual-averaged groundwater level (a.s.l.) in the Changjiang and Yellow River basins. In Fig. 3a and 3b, right axis shows a period of each crop (WH;

MZ; maize, and RI; rice, respectively).

wheat, MZ; maize, and RI; rice, respectively).

0 120 240 360 480 600 720

0 120 240 360 480 600 720

Date

0 120 240 360 480 600 720

Date

WH(=1) or MZ(=2)

WH(=1) or MZ(=2)

(a) (b)

WH(=1) or MZ(=2)

Evap. (mm/day)

RI-RI (Irrigated)

RI-RI (Irrigated)

RI-RI (Non-irrigated)

Effect of Evapotranspiration on Hydrothermal Changes in Regional Scale

Evap. (mm/day)

Water depth (mm)

(c)

0 120 240 360 480 600 720

0 120 240 360 480 600 720

Date

0 120 240 360 480 600 720

Date

G.W.L.(m)

RI(=3)

http://dx.doi.org/10.5772/52808

RI(=3)

133

RI(=3)

Date

Date

Irrigation (mm/day)

Evap. (mm/day)

WH-MZ (Irrigated)

WH-MZ (Irrigated)

WH-MZ (Non-irrigated)

Evap. (mm/day)

Effect of Evapotranspiration on Hydrothermal Changes in Regional Scale http://dx.doi.org/10.5772/52808 133

simulation was changed from Δt = 1.5 sec to 1 h depending on spatial scale and the sub-model. Simulations were validated against various hydrothermal observed variables such as river discharge, soil moisture, groundwater level, air temperature, surface temperature, and heat-

After the verification procedure (Nakayama, 2011b; Nakayama and Watanabe, 2008b), the model simulated effect of irrigation on evapotranspiration at rotation between winter wheat and summer maize in the downstream of Yellow River, and double-cropping of rice in the middle of Changjiang River (Fig. 3). Because more water is withdrawn during winter-wheat period due to small rainfall in the semi-arid north, the irrigation in this period affects greatly the increase in evapotranspiration (Fig. 3a), which is supplied by the limited water resources of river discharge and groundwater there. In particular, most of the irrigation is withdrawn from aquifer in the North China Plain (NCP) because surface water is seriously limited there (Nakayama, 2011a, 2011b; Nakayama et al., 2006). In the south, the irrigation is usually from the river to fill the paddy fiddle as ponding water depth (Nakayama, 2012c; Nakayama and Watanabe, 2008b), which increases evapotranspiration more in the drier season (Fig. 3b). This implies that energy supply is abundant relative to the water supply and the hydrological process is more sensitive to precipitation in the north, whereas the water supply is abundant relative to the energy supply and sun duration has a more significant impact in the south (Cong

The model also simulated groundwater level in both Changjiang and Yellow River basins (Fig. 3c). The level decreases rapidly around the source area and the Qinghai Tibet Plateau, indicating that there are many sources of spring water in this region. The value is very low in the downstream because of the low elevation and overexploitation, in particular, in the NCP (Nakayama, 2011a, 2011b; Nakayama et al., 2006). This result indicates that hydrologic cycle including groundwater level is highly related not only to the topography but also to the irrigation water use. The NICE is effective to provide better evaluation of hydrological trends in longer period including 'evaporation paradox' (Roderick and Farquhar, 2002; Cong et al., 2010) together with observation networks because the model does not need the crop coefficient (depending on a growing stage and a kind of crop) for the calculation of actual evaporation and simulates it directly without detailed site-specific information or empirical relation to

The mean TINDVI (Time-Integrated Normalized Difference Vegetation Index) gradients during 1982-1999 in various field crops (wheat, maize, and rice) at 4 stations were compared with trends of crop yields in the previous research (Tao et al., 2006) (Fig. 4a). The correlation of both values is relatively good (*r*<sup>2</sup> = 0.986) and the TINDVI gradient has a linear relation to the yield trend. The spatial pattern of the mean TINDVI gradient in agricultural fields shows a generally increasing tendency, particularly in the Yellow downstream and the NCP (Fig.

calculate effective precipitation (Nakayama, 2011a; Nakayama et al., 2006).

flux budget, et al.

132 Evapotranspiration - An Overview

et al., 2010).

**3. Result and discussion**

**3.1. Effect of irrigation on hydrologic change**

Fig. 3. Impact of irrigation and ponding water depth on hydrological cycle; (a) evapotranspiration change at rotation between winter wheat and summer maize in the downstream of Yellow River, (b) evapotranspiration change at doublecropping of rice in the middle of Changjiang River, and (c) simulated results of annual-averaged groundwater level **Figure 3.** Impact of irrigation and ponding water depth on hydrological cycle; (a) evapotranspiration change at rota‐ tion between winter wheat and summer maize in the downstream of Yellow River, (b) evapotranspiration change at double-cropping of rice in the middle of Changjiang River, and (c) simulated results of annual-averaged groundwater level (a.s.l.) in the Changjiang and Yellow River basins. In Fig. 3a and 3b, right axis shows a period of each crop (WH; wheat, MZ; maize, and RI; rice, respectively).

(a.s.l.) in the Changjiang and Yellow River basins. In Fig. 3a and 3b, right axis shows a period of each crop (WH; wheat,

MZ; maize, and RI; rice, respectively).

4b), which is closely related to increasing tendency for winter wheat production in the downstream (U.S. Department of Agriculture, 1994), increase in irrigation water use (Yang et al., 2004) and chemical fertilizer, changes in crop varieties, improvements in technology such as agricultural machines, and other agronomic changes. On the other hand, it shows a generally decreasing tendency, particularly in the mid-lower reaches and around the lakes in the Changjiang River. This is caused mainly by an increase in lake reclamation, levee construction, and the resultant relative decrease in rice productivity in the lower reaches (Shankman and Liang, 2003; Zhao et al., 2005; Nakayama and Watanabe, 2008b). The decrease in the TINDVI gradient near the Bohai Sea, the East China Sea, and the Taihang Mountains was due to several effects including groundwater degradation, seawater intrusion, and rapid urbanization in the areas surrounding bigger cities (Brown and Halweil, 1998). Generally, these results suggest that the increase in irrigation water use is one of the reasons for the increase in crop production (Yang et al., 2004).

of vaporization by using the water-holding pavement, and this trend continues during 5 days after the fulfillment of water-holding pavement. The rapid increase of latent heat (438 W/m2

(Yamagata et al., 2008). The simulated maximum sensible and latent heat fluxes in the water-

previous research. The model could simulate reasonably the general trend that the latent heat in the infiltration pavement was smaller than that in the water-holding pavement because the

(a)

Date

Aug 3 4 5 6 7 8

**Figure 5.** Simulated heat-flux budget in the symbiotic urban pavements after water irrigation at 3 August 2007; (a) infiltration pavement, and (b) water-holding pavement. Dotted line, net radiation (NR); solid line, sensible heat flux

The author predicted the hydrothermal changes in the symbiotic urban scenarios (Fig. 6) (Nakayama and Hashimoto, 2011; Nakayama et al., 2012). NICE correctly predicted the much lower surface temperature of the water-holding block than those of the other pavements for several days after rainfall in comparison with simplified AUSSSM (Nakayama and Fujita, 2010), which was caused mainly by the rapid increase of evaporation after the rainfall (Fig. 6b). The predicted surface temperature on the scenario of water-holding pavement shows drastic decrease in the entire Kawasaki City (Fig. 6c). In particular, the business district beside the sea, where the urban heat island is predominant mainly due to the paved surface and the

Date

Aug 3 4 5 6 7 8

(b)

and 337 W/m2

holding pavement just after water irrigation were 130 W/m2 and 345 W/m2

) was reported in the field observation at sprinkling roads

Effect of Evapotranspiration on Hydrothermal Changes in Regional Scale

NR SH LH GH

NR SH LH GH

relative to sensible heat (127 W/m2

relatively good approximations of 161 W/m2

infiltration was more predominant than the evaporation.


0

200

Heat-budget (W/m

2

)

400

600


0

200

Heat-budget (W/m

2

)

400

600

Infiltration, Cal.

Water-holding, Cal.

(SH); bold line, latent heat flux (LH); dash-dotted line, ground transfer heat flux (GH), respectively.

greater anthropogenic heat sources, has an effective cooling on this scenario.

)

135

, which had

in observed values and this

http://dx.doi.org/10.5772/52808

**Figure 4.** TINDVI over the agricultural fields by NOAA/AVHRR satellite images in the agricultural fields; (a) comparison of TINDVI with trend in observed yields at 4 stations (Pt.I; Changsha, Pt.II; Hefei, Pt.III; Zhengzhou, Pt.IV; Tianshui), and (b) spatial patterns in mean TINDVI gradient (per year) during 1982-1999.

#### **3.2. Effective management of water resources to reduce urban heat island**

The simulated heat-flux budget was compared at the infiltration and water-holding pavements after the verification with the analyzed value (Fig. 5). The water-holding pavement takes about 10.0 °C lower of the surface temperature than the infiltration pavement. The observed daily cycles of the road surface and air temperature at 1.5 m height were compared with those simulated by NICE after water irrigation in the same way as the previous research (Nakayama and Fujita, 2010). The model generally captured the observed water amount, and the associated diurnal cycle of the road surface and air temperature, and could estimate the temperature decrease trend with high accuracy. This cooling temperature is closely related to the promotion of vaporization by using the water-holding pavement, and this trend continues during 5 days after the fulfillment of water-holding pavement. The rapid increase of latent heat (438 W/m2 ) relative to sensible heat (127 W/m2 ) was reported in the field observation at sprinkling roads (Yamagata et al., 2008). The simulated maximum sensible and latent heat fluxes in the waterholding pavement just after water irrigation were 130 W/m2 and 345 W/m2 , which had relatively good approximations of 161 W/m2 and 337 W/m2 in observed values and this previous research. The model could simulate reasonably the general trend that the latent heat in the infiltration pavement was smaller than that in the water-holding pavement because the infiltration was more predominant than the evaporation.

4b), which is closely related to increasing tendency for winter wheat production in the downstream (U.S. Department of Agriculture, 1994), increase in irrigation water use (Yang et al., 2004) and chemical fertilizer, changes in crop varieties, improvements in technology such as agricultural machines, and other agronomic changes. On the other hand, it shows a generally decreasing tendency, particularly in the mid-lower reaches and around the lakes in the Changjiang River. This is caused mainly by an increase in lake reclamation, levee construction, and the resultant relative decrease in rice productivity in the lower reaches (Shankman and Liang, 2003; Zhao et al., 2005; Nakayama and Watanabe, 2008b). The decrease in the TINDVI gradient near the Bohai Sea, the East China Sea, and the Taihang Mountains was due to several effects including groundwater degradation, seawater intrusion, and rapid urbanization in the areas surrounding bigger cities (Brown and Halweil, 1998). Generally, these results suggest that the increase in irrigation water use is one of the reasons for the increase in crop production

> DTINDVI (1/yr) -1.00 -0.75 -0.50 -0.25 0.00 0.25 0.50 0.75 1.00

**Figure 4.** TINDVI over the agricultural fields by NOAA/AVHRR satellite images in the agricultural fields; (a) comparison of TINDVI with trend in observed yields at 4 stations (Pt.I; Changsha, Pt.II; Hefei, Pt.III; Zhengzhou, Pt.IV; Tianshui), and

The simulated heat-flux budget was compared at the infiltration and water-holding pavements after the verification with the analyzed value (Fig. 5). The water-holding pavement takes about 10.0 °C lower of the surface temperature than the infiltration pavement. The observed daily cycles of the road surface and air temperature at 1.5 m height were compared with those simulated by NICE after water irrigation in the same way as the previous research (Nakayama and Fujita, 2010). The model generally captured the observed water amount, and the associated diurnal cycle of the road surface and air temperature, and could estimate the temperature decrease trend with high accuracy. This cooling temperature is closely related to the promotion

**Pt.IV**

**Pt.II**

**Pt.I**

**Pt.III**

(a) (b)

(Yang et al., 2004).

134 Evapotranspiration - An Overview


D

TINDVI (1/yr)


Yield (kg/ha/yr)

(b) spatial patterns in mean TINDVI gradient (per year) during 1982-1999.

**3.2. Effective management of water resources to reduce urban heat island**

Y=0.092X, R2

=0.99

Wheat

**Pt.IV**

Rice

**Pt.II**

Rice

**Pt.I**

D

Maize

**Pt.III**

Wheat

**Figure 5.** Simulated heat-flux budget in the symbiotic urban pavements after water irrigation at 3 August 2007; (a) infiltration pavement, and (b) water-holding pavement. Dotted line, net radiation (NR); solid line, sensible heat flux (SH); bold line, latent heat flux (LH); dash-dotted line, ground transfer heat flux (GH), respectively.

The author predicted the hydrothermal changes in the symbiotic urban scenarios (Fig. 6) (Nakayama and Hashimoto, 2011; Nakayama et al., 2012). NICE correctly predicted the much lower surface temperature of the water-holding block than those of the other pavements for several days after rainfall in comparison with simplified AUSSSM (Nakayama and Fujita, 2010), which was caused mainly by the rapid increase of evaporation after the rainfall (Fig. 6b). The predicted surface temperature on the scenario of water-holding pavement shows drastic decrease in the entire Kawasaki City (Fig. 6c). In particular, the business district beside the sea, where the urban heat island is predominant mainly due to the paved surface and the greater anthropogenic heat sources, has an effective cooling on this scenario.

The predicted groundwater level change in August is interesting in contrast to the simulated temperature (Fig. 6d). In the simulation on the scenario of water-holding pavement, all the necessary water to fill the water-holding pavement was automatically simulated by consid‐ ering the difference between precipitation and evaporation, and withdrawn from the under‐ neath groundwater in the NICE, which means that the pavement was always saturated. The groundwater level would decrease drastically at the commercial and industrial areas beside the sea and at the inland residential area, in exchange for a drastic cooling in the corresponding and the surrounding areas. The increase in recharge rate and the consequent increase in groundwater level on the scenario of a natural zone and green area are also effective to promote the groundwater resources in the urban area covered by impermeable pavement (Nakayama et al., 2007, 2012). The predicted temperature and groundwater level are greatly affected not only at the business district in the Kawasaki area but also at the Tokyo metropolitan area in the northern side of this study area. This result is very important from the political point of view, which indicates that we have to estimate more precisely the arrangement of symbiotic urban scenario in the study area together with the neighbouring administrations.

**3.3. Discussion**

estimation of uncertainty (Milly et al., 2008).

the achievement in sustainable development in urban regions.

ment, 2004; Nakayama and Hashimoto, 2011; Nakayama et al., 2007, 2012).

Simulated results about the impact of irrigation on evapotranspiration change showed a clear difference between the Changjiang and Yellow Rivers. Because surface water and groundwater are administered separately by different authorities in China (Nakayama, 2011a), water management becomes further complicated if we consider it in both the Changjiang and Yellow Rivers. Any change in water accounting procedures may need to be negotiated through agreements brokered at relatively high levels of government, because surface water and groundwater are physically closely related to each other. The future development of irrigated and unirrigated fields and the associated crop production would affect greatly hydro-climate change and usable irrigation water from river and aquifer, and vice versa (Nakayama, 2011b). This research presented the lateral subsurface flow also has an important effect on the hydrologic cycle even in the continental scale, which extends traditional 'dynamic equilibrium' with atmospheric forcing (Maxwell and Kollet, 2008). From this point of view, nonstationarity models of relevant environmental variables have to be further developed to incorporate water infrastructure and water users including agricultural and energy sectors with a careful

Effect of Evapotranspiration on Hydrothermal Changes in Regional Scale

http://dx.doi.org/10.5772/52808

137

This study also showed that the water-holding pavement scenario is effective for the urban area of floating subways, stations, and buildings in order to use moderately groundwater and to ameliorate the severe heat island through promotion of evaporation, whereas alternative land cover scenario of green area is also effective for the urban area of groundwater depression and seawater intrusion in order to promote the infiltration and to cool temperature. Therefore, the effective management of water resources including groundwater as a heat sink, particu‐ larly during the summer, would be attractive for both recovering a sound hydrologic cycle and tackling the urban heat island phenomenon (Ministry of Environment, 2004; Nakayama et al., 2007, 2012), so-called, 'Win-Win' approach. We are convinced that integrated manage‐ ment of both surface water and groundwater by using NICE in a political scenario for the effective selection and use of ecosystem service sites (Millennium Ecosystem Assessment, 2005) would play an important role in the creation of thermally-pleasing environments and

Recently, environmental pollution is becoming intertwined with various aspects (Fig. 7). Land subsidence is still a serious environmental problem around the foot of mountain near the Tokyo Metropolis due to delayed regulatory enforcement and continuation of extensive groundwater extraction. Water contaminant also shows some relation to this heterogeneity in addition to the thermal environment through complex chemical reaction, which implies further importance to evaluate these problems synthetically. Present result indicates effective management of water resource including evapotranspiration is also powerful for mitigation of heat island, which implies a possibility of achieving win-win solution about hydrothermal pollutions in eco-conscious society. The procedure to construct integrated assessment system would be also valuable for adaptation to climate change and urbanization in global scale, proposal of sustainability index, and providing eco-conscious society (Ministry of Environ‐

**Figure 6.** Prediction of hydrothermal changes at 5 August 2006; (a) present surface temperature, (b) evaporation si‐ mulated by NICE and simplified AUSSSM above water-holding pavement during 8-19 August 2006, and (c-d) predic‐ tions of surface temperature and groundwater level changes after water-holding pavement.

### **3.3. Discussion**

The predicted groundwater level change in August is interesting in contrast to the simulated temperature (Fig. 6d). In the simulation on the scenario of water-holding pavement, all the necessary water to fill the water-holding pavement was automatically simulated by consid‐ ering the difference between precipitation and evaporation, and withdrawn from the under‐ neath groundwater in the NICE, which means that the pavement was always saturated. The groundwater level would decrease drastically at the commercial and industrial areas beside the sea and at the inland residential area, in exchange for a drastic cooling in the corresponding and the surrounding areas. The increase in recharge rate and the consequent increase in groundwater level on the scenario of a natural zone and green area are also effective to promote the groundwater resources in the urban area covered by impermeable pavement (Nakayama et al., 2007, 2012). The predicted temperature and groundwater level are greatly affected not only at the business district in the Kawasaki area but also at the Tokyo metropolitan area in the northern side of this study area. This result is very important from the political point of view, which indicates that we have to estimate more precisely the arrangement of symbiotic

urban scenario in the study area together with the neighbouring administrations.

0.0 0.2 0.4 0.6 0.8 1.0

**Figure 6.** Prediction of hydrothermal changes at 5 August 2006; (a) present surface temperature, (b) evaporation si‐ mulated by NICE and simplified AUSSSM above water-holding pavement during 8-19 August 2006, and (c-d) predic‐

(c) (d)

Date

**[m] 0.0 -0.2 -0.4 -0.6 -0.8 -1.0 -1.2 -1.4**

Aug 8 9 10 11 12 13 14 15 16 17 18 19

NICE-URBAN:without water effect NICE-URBAN:include water effect

(b)

Water-holding AUSSSM

Evaporation (kg)

**[ ]**

(a)

136 Evapotranspiration - An Overview

**[ ]**

tions of surface temperature and groundwater level changes after water-holding pavement.

Simulated results about the impact of irrigation on evapotranspiration change showed a clear difference between the Changjiang and Yellow Rivers. Because surface water and groundwater are administered separately by different authorities in China (Nakayama, 2011a), water management becomes further complicated if we consider it in both the Changjiang and Yellow Rivers. Any change in water accounting procedures may need to be negotiated through agreements brokered at relatively high levels of government, because surface water and groundwater are physically closely related to each other. The future development of irrigated and unirrigated fields and the associated crop production would affect greatly hydro-climate change and usable irrigation water from river and aquifer, and vice versa (Nakayama, 2011b). This research presented the lateral subsurface flow also has an important effect on the hydrologic cycle even in the continental scale, which extends traditional 'dynamic equilibrium' with atmospheric forcing (Maxwell and Kollet, 2008). From this point of view, nonstationarity models of relevant environmental variables have to be further developed to incorporate water infrastructure and water users including agricultural and energy sectors with a careful estimation of uncertainty (Milly et al., 2008).

This study also showed that the water-holding pavement scenario is effective for the urban area of floating subways, stations, and buildings in order to use moderately groundwater and to ameliorate the severe heat island through promotion of evaporation, whereas alternative land cover scenario of green area is also effective for the urban area of groundwater depression and seawater intrusion in order to promote the infiltration and to cool temperature. Therefore, the effective management of water resources including groundwater as a heat sink, particu‐ larly during the summer, would be attractive for both recovering a sound hydrologic cycle and tackling the urban heat island phenomenon (Ministry of Environment, 2004; Nakayama et al., 2007, 2012), so-called, 'Win-Win' approach. We are convinced that integrated manage‐ ment of both surface water and groundwater by using NICE in a political scenario for the effective selection and use of ecosystem service sites (Millennium Ecosystem Assessment, 2005) would play an important role in the creation of thermally-pleasing environments and the achievement in sustainable development in urban regions.

Recently, environmental pollution is becoming intertwined with various aspects (Fig. 7). Land subsidence is still a serious environmental problem around the foot of mountain near the Tokyo Metropolis due to delayed regulatory enforcement and continuation of extensive groundwater extraction. Water contaminant also shows some relation to this heterogeneity in addition to the thermal environment through complex chemical reaction, which implies further importance to evaluate these problems synthetically. Present result indicates effective management of water resource including evapotranspiration is also powerful for mitigation of heat island, which implies a possibility of achieving win-win solution about hydrothermal pollutions in eco-conscious society. The procedure to construct integrated assessment system would be also valuable for adaptation to climate change and urbanization in global scale, proposal of sustainability index, and providing eco-conscious society (Ministry of Environ‐ ment, 2004; Nakayama and Hashimoto, 2011; Nakayama et al., 2007, 2012).

**Acknowledgements**

acknowledged.

**Author details**

UK

**References**

Tadanobu Nakayama1,2

The author thanks Dr. M. Watanabe, Keio University, Japan, Dr. Y. Yang, Shijiazhuang Institute of Agricultural Modernization of the Chinese Academy of Sciences (CAS), China, and Dr. K. Xu, National Institute for Environmental Studies (NIES), Japan, for valuable comments about the study areas. Some of the simulations in this study were run on an NEC SX–8 supercomputer at the Center for Global Environmental Research (CGER), NIES. The support of the Environ‐ mental Technology Development Fund from the Japanese Ministry of Environment is also

Effect of Evapotranspiration on Hydrothermal Changes in Regional Scale

http://dx.doi.org/10.5772/52808

139

1 National Institute for Environmental Studies (NIES), Onogawa, Tsukuba, Ibaraki, Japan

2 Centre for Ecology & Hydrology (CEH), Crowmarsh Gifford, Wallingford, Oxfordshire,

[1] Brown, L. R, & Halweil, B. (1998). China's water shortage could shake world food

[2] Chang, C. R, Li, M. H, & Chang, S. D. (2007). A preliminary study on the local coolisland intensity of Taipei city parks. *Landscape Urban Plan.*, 0169-2046, 80, 386-395. [3] Chinese Academy of Sciences(1988). *Administrative division coding system of the People's*

[4] Cong, Z, Zhao, J, Yang, D, & Ni, G. (2010). Understanding the hydrological trends of river basins in China. *J. Hydrol.*, doi:j.jhydrol.2010.05.013, 0022-1694, 388, 350-356. [5] Digital National Land Information GIS Data of Japan(2002). *Database of groundwater in Japan.* Ministry of Land Infrastructure and Transport of Japan, http://

[7] Japan Meteorological Agency (JMA)(2005). *AMeDAS (Automated Meteorological Data Acquisition System) Annual Reports 20052006*Japan Meteorological Business Support

[8] Japan Meteorological Agency (JMA)(2005). *MSM (Meso Scale Model) Objective Analysis*

[6] Geological Atlas of China(2002). Geological Publisher, Beijing (in Chinese)

*Data 20052006*Japan Meteorological Business Support Center (CD-ROM)

security. *World Watch*, 0896-0615, 11, 10-18.

*Republic of China*, Beijing

www.nla.go.jp/ksj/

Center (CD-ROM)

**Figure 7.** Intertwined environmental pollution in urban and surrounding areas.

### **4. Conclusion**

This study coupled National Integrated Catchment-based Eco-hydrology (NICE) model series with complex sub-models to develop coupled human and natural systems and to analyze impact of evapotranspiration on hydrothermal changes in regional scale. The model includes different functions of representative crops (wheat, maize, soybean, and rice) and simulates automatically dynamic growth processes and biomass formulation. The simulated result showed impact of irrigation on eco-hydrological processes. The spatial pattern of Time-Integrated NDVI (TINDVI) gradient in agricultural fields indicated heterogeneous character‐ istics of crop yield, which implied the increase in irrigation water use is one of the reasons for the increase in crop production. NICE also reproduced reasonably the observed hydrothermal characteristics including the water and heat budgets in various pavements, evaluated the role of a new surface material (water-holding pavement) in promoting evaporation and cooling temperature to counter the urban heat island phenomenon, and predicted the hydrothermal changes under alternative land cover scenarios. These results suggest strongly the needs of trans-boundary and multi-disciplinary solutions of water management for sustainable development under sound socio-economic conditions contributory to national and global securities.

### **Acknowledgements**

The author thanks Dr. M. Watanabe, Keio University, Japan, Dr. Y. Yang, Shijiazhuang Institute of Agricultural Modernization of the Chinese Academy of Sciences (CAS), China, and Dr. K. Xu, National Institute for Environmental Studies (NIES), Japan, for valuable comments about the study areas. Some of the simulations in this study were run on an NEC SX–8 supercomputer at the Center for Global Environmental Research (CGER), NIES. The support of the Environ‐ mental Technology Development Fund from the Japanese Ministry of Environment is also acknowledged.

### **Author details**

Tadanobu Nakayama1,2

1 National Institute for Environmental Studies (NIES), Onogawa, Tsukuba, Ibaraki, Japan

2 Centre for Ecology & Hydrology (CEH), Crowmarsh Gifford, Wallingford, Oxfordshire, UK

### **References**

**Simulated G.W.L.** 

**Figure 7.** Intertwined environmental pollution in urban and surrounding areas.

**Tokyo Sta.**

138 Evapotranspiration - An Overview

**Anchors (130x200kg) are buried to prevent from floating**

**4. Conclusion**

securities.

**relative to ground surface Groundwater nitrogen conc. interpolated by obs. data**

This study coupled National Integrated Catchment-based Eco-hydrology (NICE) model series with complex sub-models to develop coupled human and natural systems and to analyze impact of evapotranspiration on hydrothermal changes in regional scale. The model includes different functions of representative crops (wheat, maize, soybean, and rice) and simulates automatically dynamic growth processes and biomass formulation. The simulated result showed impact of irrigation on eco-hydrological processes. The spatial pattern of Time-Integrated NDVI (TINDVI) gradient in agricultural fields indicated heterogeneous character‐ istics of crop yield, which implied the increase in irrigation water use is one of the reasons for the increase in crop production. NICE also reproduced reasonably the observed hydrothermal characteristics including the water and heat budgets in various pavements, evaluated the role of a new surface material (water-holding pavement) in promoting evaporation and cooling temperature to counter the urban heat island phenomenon, and predicted the hydrothermal changes under alternative land cover scenarios. These results suggest strongly the needs of trans-boundary and multi-disciplinary solutions of water management for sustainable development under sound socio-economic conditions contributory to national and global

**<Imbalance of hydrologic cycle> <Water contamination>**

**Thermal environment map (Tokyo Metro. Gov.)**

**<Heat island>**


[9] Maxwell, R. M, & Kollet, S. J. (2008). Interdependence of groundwater dynamics and land-energy feedbacks under climate change. *Nat. Geosci.*, doi:ngeo315, 1752-0894, 1, 665-669.

N. Gallo & Marco H. Ferrari (Eds.), Nova Science Publishers, Inc.,

Effect of Evapotranspiration on Hydrothermal Changes in Regional Scale

http://dx.doi.org/10.5772/52808

141

[23] Nakayama, T. (2010). Simulation of hydrologic and geomorphic changes affecting a

[24] Nakayama, T, & Fujita, T. (2010). Cooling effect of water-holding pavements made of new materials on water and heat budgets in urban areas. *Landscape Urban Plan.*,

[25] Nakayama, T, Sun, Y, & Geng, Y. (2010). Simulation of water resource and its relation to urban activity in Dalian City, Northern China. *Global Planet. Change*, doi:j.gloplacha.

[26] Nakayama, T. (2011a). Simulation of complicated and diverse water system accompa‐ nied by human intervention in the North China Plain. *Hydrol. Process.*, doi:hyp.8009,

[27] Nakayama, T. (2011b). Simulation of the effect of irrigation on the hydrologic cycle in the highly cultivated Yellow River Basin. *Agr. Forest Meteorol.*, doi:j.agrformet.

[28] Nakayama, T, & Hashimoto, S. (2011). Analysis of the ability of water resources to reduce the urban heat island in the Tokyo megalopolis. *Environ. Pollut.*, doi:j.envpol.

[29] Nakayama, T. (2012a). Visualization of missing role of hydrothermal interactions in Japanese megalopolis for win-win solution. *Water Sci. Technol.*, doi:wst.2012.205,

[30] Nakayama, T. and regime shift of mire ecosystem in northern Japan. *Hydrol. Process.*,

[31] Nakayama, T. (2012c). Impact of anthropogenic activity on eco-hydrological process in continental scales. *Proc. Environ. Sci.*, doi:j.proenv.2012.01.008, 1878-0296, 13, 87-94.

[32] Nakayama, T, Hashimoto, S, & Hamano, H. (2012). Multi-scaled analysis of hydro‐ thermal dynamics in Japanese megalopolis by using integrated approach. *Hydrol.*

[33] Oke, T. R. (1987). *Boundary layer climates*, Methuen Press, 978-0-41504-319-9London.

[34] Pielke, R. A, Cotton, W. R, Walko, R. L, Tremback, C. J, Lyons, W. A, Grasso, L. D, Nicholls, M. E, Moran, M. D, Wesley, D. A, Lee, T. J, & Copeland, J. H. (1992). A omprehensive meteorological modeling system-RAMS. *Meteorol. Atmos. Phys.*,

[35] Priestley, C. H. B, & Taylor, R. J. (1972). On the assessment of surface heat flux and evaporation using large-scale parameters. *Mon. Weather Rev.*, 0027-0644, 100, 81-92.

shrinking mire. *River Res. Appl.*, doi:rra.1253, 1535-1459, 26, 305-321.

doi:j.landurbplan.2010.02.003, 0169-2046, 96, 57-67.

978-1-60456-643-7New York, 1-89.

2010.06.001, 0921-8181, 73, 172-185.

2010.11.006, 0168-1923, 151, 314-327.

2010.11.016, 0269-7491, 159, 2164-2173.

doi:hyp.9347, 0885-6087, 26, 2455-2469.

*Process.*, doi:hyp.9290, 0885-6087, 26, 2431-2444.

0885-6087, 0885-6087.

0273-1223, 66, 409-414.

0177-7971, 49, 69-91.


N. Gallo & Marco H. Ferrari (Eds.), Nova Science Publishers, Inc., 978-1-60456-643-7New York, 1-89.

[23] Nakayama, T. (2010). Simulation of hydrologic and geomorphic changes affecting a shrinking mire. *River Res. Appl.*, doi:rra.1253, 1535-1459, 26, 305-321.

[9] Maxwell, R. M, & Kollet, S. J. (2008). Interdependence of groundwater dynamics and land-energy feedbacks under climate change. *Nat. Geosci.*, doi:ngeo315, 1752-0894, 1,

[10] Millennium Ecosystem Assessment(2005). *Strengthening capacity to manage ecosystems sustainability for human well-being*, http://www.millenniumassessment.org/en/

[11] Milly, P. C. D, Betancourt, J, Falkenmark, M, Hirsch, R. M, Kundzewicz, Z. W, Letten‐ maier, D. P, & Stouffer, R. J. (2008). Stationarity is dead: Whither water management ?.

[12] Ministry of Environment(2004). *Report on heat-island measures by controlling anthropo‐ genic exhaustion heat in the urban area*, http://www.env.go.jp/air/report/hin Japanese) [13] Nakayama, T, & Watanabe, M. (2004). Simulation of drying phenomena associated with vegetation change caused by invasion of alder (Alnus japonica) in Kushiro Mire. *Water*

[14] Nakayama, T, & Watanabe, M. (2006). Simulation of spring snowmelt runoff by considering micro-topography and phase changes in soil layer. *Hydrol. Earth Syst. Sci.*

[15] Nakayama, T, Yang, Y, Watanabe, M, & Zhang, X. (2006). Simulation of groundwater dynamics in North China Plain by coupled hydrology and agricultural models. *Hydrol.*

[16] Nakayama, T, Watanabe, M, Tanji, K, & Morioka, T. (2007). Effect of underground urban structures on eutrophic coastal environment. *Sci. Total Environ.*, doi:j.scitotenv.

[17] Nakayama, T. (2008a). Factors controlling vegetation succession in Kushiro Mire. *Ecol.*

[18] Nakayama, T. (2008b). Shrinkage of shrub forest and recovery of mire ecosystem by river restoration in northern Japan. *Forest Ecol. Manag.*, doi:j.foreco.2008.07.017,

[19] Nakayama, T, & Watanabe, M. (2008a). Missing role of groundwater in water and nutrient cycles in the shallow eutrophic Lake Kasumigaura, Japan. *Hydrol. Process.*,

[20] Nakayama, T, & Watanabe, M. (2008b). Role of flood storage ability of lakes in the Changjiang River catchment. *Global Planet. Change*, doi:j.gloplacha.2008.04.002,

[21] Nakayama, T, & Watanabe, M. (2008c). Modelling the hydrologic cycle in a shallow

[22] Nakayama, T. (2009). Simulation of Ecosystem Degradation and its Application for Effective Policy-Making in Regional Scale, In: *River Pollution Research Progress*, Mattia

eutrophic lake. *Verh. Internat. Verein. Limnol.*, 0368-0770, 30, 345-348.

*Science*, doi:science.1151915, 0036-8075, 319, 573-574.

*Resour. Res.*, W08402, doi:WR003174, 0043-1397, 40

*Process.*, doi:hyp.6142, 0885-6087, 20, 3441-3466.

*Model.*, doi:j.ecolmodel.2008.02.017, 0304-3800, 215, 225-236.

*Discuss.*, 1027-5606, 3, 2101-2144.

2006.11.033, 0048-9697, 373, 270-288.

doi:hyp.6684, 0885-6087, 22, 1150-1172.

0378-1127, 256, 1927-1938.

0921-8181, 63, 9-22.

665-669.

140 Evapotranspiration - An Overview

index.aspx


[36] Ritchie, J. T, Singh, U, Godwin, D. C, & Bowen, W. T. (1998). Cereal growth, develop‐ ment and yield, In: *Understanding Options for Agricultural Production*, Tsuji, G.Y., Hoogenboom, G. & Thornton, P.K. (Eds.), Kluwer, 0-79234-833-8Britain, 79-98.

**Chapter 8**

**A Parametric Model for Potential**

Aristoteles Tegos, Andreas Efstratiadis and

Additional information is available at the end of the chapter

tion, while in cold hydroclimates this percentage is lower that 30%.

Demetris Koutsoyiannis

http://dx.doi.org/10.5772/52927

**1. Introduction**

**Evapotranspiration Estimation Based on a Simplified**

The estimation of evapotranspiration is a typical task in several disciplines of geosciences and engineering, including hydrology, meteorology, climatology, ecology and agricultural engineering. According to the specific area of interest, we distinguish between three differ‐ ent aspects of evapotranspiration, i.e. actual, potential and reference crop. In particular, the actual evapotranspiration refers to the quantity of water that is actually removed from a sur‐ face due to the combined processes of evaporation and transpiration. At the global scale, it is estimated that about 60% of precipitation reaching the Earth's terrain returns to the atmos‐ phere by means of water losses due to evapotranspiration [1]. Obviously this percentage is not constant but exhibits significant variability both in space and time. In particular, in semiarid basins the annual percentage of evapotranspiration losses may be 70-80% of precipita‐

The actual evapotranspiration is one of the most difficult processes to measure in the field, except for experimental, small-scale areas, for which lysimeter observations can be assumed representative. For the spatial scale of practical applications (e.g., river basin), there is a growing interest on remote sensing technologies, which attempt to estimate actual evapo‐ transpiration at regional scales, by combining ground measurements with satellite-derived data [2]. Otherwise, the only reliable method is based on the calculation of the water balance of the basin. Apparently, this approach is valid only when actual evapotranspiration is the single unknown component of the water balance. In addition, it can be implemented only for large time scales (annual and over-annual), for which storage regulation effects can be

> © 2013 Tegos et al.; licensee InTech. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/3.0), which permits unrestricted use,

© 2013 Tegos et al.; licensee InTech. This is a paper distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

distribution, and reproduction in any medium, provided the original work is properly cited.

**Formulation of the Penman-Monteith Equation**


**Chapter 8**
