The present study addresses the possible effects of soil moisture changes on the simulated daily maximum and minimum air temperatures of Australia for a duration of 13 years. Therefore, the community land model version 4.5 (CLM4.5; coupled to the RegCM4) was used to represent the soil moisture and processes associated with it. The CLM4.5 has two land-surface hydrology schemes: TOPMODEL (TOP) and Variable Infiltration Capacity (VIC) and two simulations were conducted, namely: TOP and VIC. The results showed that VIC has lower soil moisture than TOP, leading to a decrease in vegetation transpiration, evaporation, and an increase in soil evaporation relative to TOP. However, there is no considerable difference between the two simulations compared with reanalysis products. In comparison to in-situ measurements, the RegCM4 can reasonably model the climatological annual cycle of mean air temperature (TMP) and its performance varies with the study site (e.g., RegCM4 overestimates TMP by 2.76 and 5.46 °C at Yanco and Tumbarumba, respectively). In summary, the simulated maximum and minimum air temperatures are sensitive to the physical parameterization of RegCM4 rather than variations in soil moisture. Likewise, improvements to the land-surface hydrology schemes TOP/VIC are required to better model Australia's daily maximum and minimum air temperatures.

  • RegCM4 succeeds in reproducing the spatial pattern of sensible and latent heat fluxes with respect to the ERA5. However, the difference between the two simulations themselves depends on the region as well as the austral season.

  • Soil moisture changes do not impose a notable impact on the simulated surface energy balance, total cloud cover, and surface net radiation with respect to reanalysis products.

Graphical Abstract

Graphical Abstract
Graphical Abstract

The climate system comprises several interactive components such as the atmosphere, biosphere, hydrosphere, cryosphere, and geosphere. These components interact with each other on multiple time scales ranging from days, seasons, and years to millennia with complex feedback mechanisms. In particular, studying the hydrological cycle is important because of the considerable influence of climate change on the water cycle budget such as precipitation, soil moisture, surface and sub-surface runoff, and evapotranspiration (Bouraoui et al. 2004; Imbach et al. 2012; Allan et al. 2020). In return, the hydrological cycle affects the climate system through the transfer of water vapor to the atmosphere. Concerning soil moisture, the hydrological cycle can be also examined through changes in the balance between the total precipitation as an input, runoff, and total evapotranspiration as an output (Peng et al. 2019; Pereira et al. 2020). Moreover, there is a direct link between the hydrology cycle and the surface energy balance and eventually to the surface climate, because solar radiation drives the vertical transfer of water vapor from the earth's surface to the atmosphere via evaporation from bare soil and transpiration from vegetation (Siler et al. 2018).

Due to the significance of the land surface condition in modeling the regional surface climate; several researches discussed the comparison between various land surface model versions. In terms of reproducing mean air temperature and total surface precipitation, the Community Land Model version 3.5 (CLM3.5; Oleson et al. 2008) outperforms the Biosphere-Atmosphere Transfer System (BATS; Dickinson et al. 1993) as reported by Steiner et al. (2009), Wang et al. (2015), and Maurya et al. (2017). Additionally, when it comes to modeling mean air temperature and total precipitation, the Community land model version 4.5 (CLM4.5; Oleson et al. 2013) performs better than the BATS scheme (Maurya et al. 2017; Chung et al. 2018).

Soil moisture plays an important role in controlling the climate system, particularly in semi-arid and arid regions, which account for 40% of the world's area (Reynolds et al. 2007). Such importance is given to investigating the factors that control soil moisture variability (Srivastava et al. 2021a). Furthermore, soil moisture derives from physiological and biogeochemical processes such as plant transpiration and photosynthesis (Seneviratne et al. 2010; Lemoine & Budny 2022). Offline land-surface (or regional climate models; RCMs) are considered important tools to investigate the response of surface climate/terrestrial carbon fluxes to soil moisture changes. For instance, Lei et al. (2014) used the offline community land model version 4 (CLM4.0; Oleson et al. 2010) to examine the role of soil moisture changes (represented by the land-surface hydrology scheme) on the terrestrial gross primary production (GPP). They reported that soil moisture considerably affects the GPP through its influence on photosynthesis and leaf area index.

Concerning surface energy balance and surface climate, RCMs were used to address this issue. For example, Anwar et al. (2019) examined the influence of two land-surface hydrology schemes (TOP and VIC) on the surface climate over Africa using version four of RegCM (RegCM4; Giorgi et al. 2012). They found that the VIC scheme outperforms the TOP scheme compared with reanalysis products despite noted biases in the summer season. In one of the recent studies, Anwar et al. (2022) examined the influence of soil moisture changes on the surface energy balance and climate of the Amazon (AMZ) using RegCM4. They found that the VIC scheme underestimates (overestimates) latent heat (sensible heat) fluxes more than TOP. In addition, the VIC showed a high warm bias than the one observed in the TOP scheme. Moreover, no scheme performs better than the other with respect to reanalysis products. Wang et al. (2021) examined the mutual interaction between different convection and land-surface hydrology schemes in simulating the Tibetan Plateau (TP) climate using the RegCM4 model. They showed that the VIC land-surface hydrology scheme reproduces the climate aspects more than TOP. Anwar & Diallo (2021) showed that RegCM4 model bias is amplified through the influence of soil moisture on vegetation status (e.g., Leaf Area Index; LAI).

Until now, the sensitivity of the Australian's daily maximum and minimum air temperatures to varied land-surface/soil moisture variations has received less attention when using a regional climate model (e.g., RegCM4); the target of our study is to address the following points:

  • 1.

    Evaluate the potential influence of soil moisture changes on Australia's maximum and minimum air temperatures using the regional climate model (RegCM4).

  • 2.

    Compare the results of this work to those of Anwar et al. (2019, 2022) and Wang et al. (2021) to investigate how the climatic regime affects the action mechanism of the land-surface hydrological scheme.

Section 2 describes the study area, data, and methods used, while Section 3 provides the results of simulations. Section 4 provides the discussion; while Section 5 shows the final conclusion and future recommendations.

Model description

The present study used the International Center of Theoretical Physics (ICTP) regional climate model version 4.7 (RegCM-4.7.0; hereafter RegCM4; Giorgi et al. 2012). Giorgi et al. (2012) reported that the RegCM4 model performance shows a remarkable improvement in comparison with previous versions: RegCM version 3 (RegCM3; Pal et al. 2007). In addition, RegCM4 includes a tropical band mode (RegT-Band; Coppola et al. 2012) to simulate the tropical climate processes (e.g., El Niño-Southern Oscillation; ENSO). Besides, it includes an interactive chemistry module, online gas-phase chemistry scheme (Shalaby et al. 2012). The performance of the RegCM4 was extensively evaluated using different physics configurations across the globe: South America (Llopart et al. 2017), Southeast Asia (Chung et al. 2018; Zhengqi et al. 2020), and Tibetan Plateau (Gu & Wang 2020) and India (Verma & Bhatla 2021). Recently, the RegCM4 has been upgraded to work as an earth system model (RegCM4-ES; Reale et al. 2020) with implementation of a non-hydrostatic dynamical core (Coppola et al. 2021).

In this study, the following physical schemes were used: (1) Rapid Radiation Transfer Model (RRTM; Clough et al. 2005) as the radiation scheme, (2) local prognostic 1.5-order scheme of the University of Washington (UW; Grenier & Bretherton 2001) as the boundary layer scheme; and (3) convection scheme of Emanuel (1991) over land and ocean. It is worth noting that the community land model version 5 (CLM5.0; Lawrence et al. 2019) has not been yet coupled to the RegCM4. In the present study, the CLM4.5 was used as the land surface model for several reasons such as: (1) the CLM4.5 offers substantial improvements such as: a revised canopy radiation scheme, canopy scaling of leaf processes, and reducing the excessive tropical GPP (Bonan et al. 2011) and (2) in previous studies, the CLM4.5 land model was used to investigate how changes in soil moisture affected surface climate simulations of the surface climate over various geographic locations (Anwar et al. 2019, 2022).

The CLM4.5 land surface model offers two land-surface hydrology schemes: (1) TOPMODEL-based (TOP; Niu et al. 2005) and (2) Variable Infiltration Capacity (VIC; Liang et al. 1994; Srivastava et al. 2017). The reader can find more details regarding the two land-surface hydrology schemes in Anwar et al. (2019, 2021), Srivastava et al. (2018, 2020), and Kumari et al. (2021). Also, it offers different options for treating the vegetation status: static vegetation (satellite phenology – SP; Lawrence et al. 2011), dynamic vegetation cover (Thornton & Rosenbloom 2005), and dynamic vegetation fraction (DV; Levis et al. 2004). In the CLM4.5, the total sensible heat flux is calculated for bare soil and vegetation land units as:
formula
(1)
where SHF is the total sensible heat flux (W m−2), SHG is the sensible heat flux emitted from the ground, and SHV is the sensible heat flux emitted from the vegetation. Total evapotranspiration (ET; kg m−2 s−1) is the sum of water vapor flux due to vegetation evaporation (QVEGE), vegetation transpiration (QVEGT), and soil evaporation (QSOIL):
formula
(2)
The RegCM4 calculates the surface net radiation (Rn) following Peixoto & Oort (1992) as:
formula
(3)
where ‘down’ subscript refers to the energy that reaches the surface and ‘up’ the energy that leaves the surface. Rn is partitioned between SHF, latent heat flux (LE), and heat storage in soil (G) as:
formula
(4)
Note that SHF and LE are only parameterized in the RegCM4, while G is not parameterized until the present day. Also, the RegCM4 computes the maximum air temperature (TMX) at time 12:00 Universal Time Coordinated (UTC), and the minimum air temperature (TMN) at time 00:00 UTC. The mean air temperature is calculated as the average of TMX and TMN as:
formula
(5)

Study area and experimental design

Australia continent's climate is characterized by a large spatial variability owing to its geographical nature. The majority of its land is classified as arid to semi-arid. An important future of Australia is its relatively flat landscape. Along the coastal strip, the climate regime shows a considerable gradient from humid-subtropical at its north to temperate at its south (Peel et al. 2007). An important feature of Northern Australia is its tropical savanna climate with a monsoonal wet season. Specht (1970) reported that land unit/vegetation of the Australian continent can be categorized as closed-forest, open-forest, woodland, shrub-land, scrub, health, and herbland. The RegCM4 model domain was customized with 25°S latitude, 132°E longitude, 60 km horizontal grid spacing and 80 grid points in the zonal direction, and 60 grid points in the meridional direction and 18 vertical levels.

Figure 1 shows the surface model elevation including location of six stations: Adelaide River, Alice Springs, Howard Springs, Robson Creek, Yanco, and Tumbarumba (to evaluate the RegCM4 performance in simulating the mean air temperature; Section 3.5), while Supplementary Figure S1 shows the seven states of Australia.
Figure 1

The surface model elevation of the Australian domain (in meters) including the six stations from the Ozflux network.

Figure 1

The surface model elevation of the Australian domain (in meters) including the six stations from the Ozflux network.

Close modal

Note that, elevation above 1,200 m will not appear (see Figure 1) because of using quite coarse resolution (60 km). Two 13-year simulations were conducted and integrated from 01 January 1998 till 31 December 2010. This period was chosen based on availability of the computational power at the time of conducting the RegCM4 simulations. The first two years of each simulation were not considered in the analysis to allow equilibration of the soil moisture (Anwar & Diallo 2022). Therefore, the actual period of analysis starts at 01 January 2000 till 31 December 2010. The two simulations were driven by the National Centre for Environmental Prediction/National Centre for Atmospheric Reanalysis version 2 of 1.8 × 1.8 degrees (NCEP/NCAR-R2; Kanamitsu et al. 2002) and the National Oceanic and Atmospheric Administration (NOAA) Optimum Interpolation Sea Surface Temperature of 1 × 1 degrees (OISST) (Reynolds et al. 2002). It is worth mentioning that the utilized version of the RegCM4 does not support ERA5 as a lateral boundary condition. Also, sensitivity of RegCM4 to atmospheric forcing (and its horizontal resolution) is beyond the scope of the present study.

The two simulations considered the satellite phenology mode (SP; Lawrence et al. 2011). In this mode, the vegetation parameters are prescribed from the Moderate Resolution Imaging Spectroradiometer dataset (MODIS; Lawrence & Chase 2007). Besides, soil moisture changes and climate variability do not affect the vegetation status (Anwar 2019). The first simulation was referred to as TOP (control simulation), while the other simulation was designated as VIC (test simulation). The difference between each simulation and the observational dataset was referred to as TOP-OBS for the TOP simulation, VIC-OBS for the VIC simulation, and the difference between the two simulations themselves (Diff). OBS is the reanalysis product (with which RegCM4 simulations were compared).

Observational data

To assess the RegCM4 model performance, various datasets were utilized. For instance, the ERA5 reanalysis product was used to evaluate the simulated surface net radiation, total cloud cover (TCC), sensible, and latent heat fluxes. ERA5 was used in the present study because of its higher temporal and spatial resolutions (0.25° horizontal grid spacing; Hersbach et al. 2020; Tarek et al. 2020). Also, it was used to evaluate the RegCM4 performance over the AMZ (Anwar et al. 2022). ERA5-land product was used to evaluate the simulated sensible heat flux as an example and it is available at 0.1° grid spacing (Muñoz-Sabater et al. 2021). Another reason to use the ERA5-land product is to check whether the simulated sensible heat flux is sensitive to the spatial resolution of the ERA5/ERA5-land reanalysis product. To account for the uncertainty associated with the observational dataset, the simulated daily maximum and minimum air temperatures were compared with two reanalysis products. The first one is the NOAA-CIRES-DOE 20th Century Reanalysis V3 (20CRv3; Compo et al. 2011; for short Century).

Century is available at 1° × 1° global grid and it covers the period of January 1836 to December 2015. In addition, it reasonably produces good estimates of the atmospheric variables throughout the 20th Century. Moreover, it has a good performance on long time scales (e.g., climate scale; Slivinski et al. 2019). The second product is the Australian Gridded Climate Data (AGCD; Jones et al. 2009). AGCD is the Bureau of Meteorology's official dataset, and it offers a high-quality set of historical and ongoing real-time climate analyses for Australia for the meteorological variables: rainfall, temperature (maximum and minimum), as well as vapor pressure at daily and monthly timescales. Also, it is available in 0.05° horizontal grid spacing and it covers the period of 1910–2020.

It is worth mentioning that the accuracy of AGCD is limited to regions with a low density of stations. To evaluate the RegCM4 model performance in simulating the monthly and climatological seasonal cycle of the mean air temperature (TMP) over the period 2000–2010, station data of the Ozflux network were used (Srivastava et al. 2021b). To serve the purpose of the present study, the time period of 2000–2010 was chosen for the ERA5, ERA5-land, Century, and AGCD. Moreover, the aforementioned reanalysis products were bi-linearly interpolated on the RegCM4-CLM4.5 curvilinear grid (Anwar & Diallo 2022).

In this section, two steps were carried out to analyze the role of land-surface hydrology schemes in controlling the maximum, minimum, and mean air temperatures of Australia following Anwar et al. (2022). Firstly, a diagnostic analysis was applied to the following variables: infiltration rate (QINFL), soil moisture of depth of 10 cm (SM10), vegetation transpiration (QVEGT), vegetation evaporation (QVEGE), soil evaporation (QSOIL), and solar radiation absorbed by vegetation and ground (SABV and SABG), sensible heat emitted from the ground and vegetation (SHG and SHV), as well as relative humidity (RH). Such analysis was performed by calculating the significant difference between the two simulations using the Student's t-test with a significance level of 0.05 (α = 5%).

These variables were retrieved from the CLM4.5 output results using the clm4.5_1dto2d program of the RegCM4. The second step considers a quantitative evaluation of the surface net radiation (Rn), TCC, surface energy balance (sensible and latent heat fluxes; SHF and LE), and maximum and minimum air temperatures (TMX and TMN) with respect to ERA5/ERA5-land, Century, and AGCD reanalysis products. The significant bias/difference was also calculated using a Student's t-test always with a significance level of 0.05 (α = 5%).

Influence on hydrology and energy variables

Figure 2 shows the simulated infiltration rate, soil moisture of depth 10 cm, vegetation transpiration, vegetation evaporation, soil evaporation, and RH for the two simulations as well as the significant difference between them. From Figure 2, it can be observed that both simulations exhibit low values of the infiltration rate (Figure 2(a) and 2(b)). Also, VIC has a lower QINFL than TOP by 0.5–2 mm day−1 approaching 3–4 mm day−1 over coastal regions. This difference is quite small if it is compared with the results reported by Anwar et al. (2019, 2022). This minor difference (Figure 2(c)) can be attributed to the following reasons: (1) Australia's continent is dominated by evergreen shrubs, short, tall grass vegetation species, and semi-desert nature, (2) rough calibration of the one-degree VIC surface dataset with respect to the Global Runoff Data Centre (Huang & Liang 2006), and (3) using high values of the maximum precipitation efficiency over land (epmax_lnd = 0.999) and rainfall speed (omtrain = 90 m s−1) in comparison with the values adopted in Africa and Amazon (Anwar et al. 2019, 2022).
Figure 2

Averaged infiltration rate (QINFL; in mm day−1; a–c); soil moisture at depth 10 cm (SM10; in mm; d–f); vegetation transpiration (QVEGT; in mm day−1; g–i); vegetation evaporation (QVEGE; in mm day−1; j–l); soil evaporation (QSOIL; fifth row in mm day−1; m–o); 2-m relative humidity (RH; sixth row in %; p–r) over the period 2000–2010. For each row, TOP simulation is on the left, VIC simulation is in the middle and the difference between the two simulations on the right (VIC minus TOP; Diff). Significant changes are indicated in black dots using the Student's t-test with α = 0.05.

Figure 2

Averaged infiltration rate (QINFL; in mm day−1; a–c); soil moisture at depth 10 cm (SM10; in mm; d–f); vegetation transpiration (QVEGT; in mm day−1; g–i); vegetation evaporation (QVEGE; in mm day−1; j–l); soil evaporation (QSOIL; fifth row in mm day−1; m–o); 2-m relative humidity (RH; sixth row in %; p–r) over the period 2000–2010. For each row, TOP simulation is on the left, VIC simulation is in the middle and the difference between the two simulations on the right (VIC minus TOP; Diff). Significant changes are indicated in black dots using the Student's t-test with α = 0.05.

Close modal

Therefore, the two simulations show a quite similar pattern of SM10 (Figure 2(d) and 2(e)), and the difference between them ranges from 2 to 4 mm approaching 6 mm over Western Australia and coastal regions (Figure 2(f)). The noted behavior of the QINFL and SM10 is reflected in partitioning the total ET budget. For instance, the two simulations show a quite similar pattern of the QVEGT (Figure 2(g) and 2(h)) and the difference between them is around 0.2–0.4 mm day−1 (Figure 2(i)). For the QVEGE, the two simulations match each other (Figure 2(j) and 2(k)) and there is no obvious difference between them (Figure 2(l)). The situation for the QSOIL is quite different (Figure 2(m) and 2(n)) because the difference between them is ranging between 0.4 and 0.8 mm day−1 approaching 1 mm day−1 over all the continent except for the Northern part of West Australia where the difference ranges from −0.2 to −0.4 mm day−1 (Figure 2(o)).

Regarding RH, the two simulations show low values between 20 and 30% in the inland regions of Northern Territory, Queensland, South Australia, and New South Wales with high values of 70–90% in the coastal borders of Western Australia, Northern Territory, and Queensland (Figure 2(p) and 2(q)). Likewise, the significant difference between the two simulations is ∼5% (Figure 2(r)). This quite small difference can be attributed to the small difference between the two simulations regarding QVEGT. The influence of soil moisture changes was also examined for the energy fluxes: SABG, SABV, SHG, and SHV. Figure 3 shows the simulated SABG, SABV, SHG, and SHV for the two simulations as well as the significant difference between them. From Figure 3(a)–3(c), it can be observed that the VIC has SABG than TOP by 10–30 W m−2 because VIC has lower SM10 than TOP as discussed earlier. For SABV (Figure 3(d)–3(f)), it can be shown that there is almost no obvious difference between the two simulations because vegetation cover and fraction are prescribed from the MODIS (Section 2.2). Furthermore, the VIC has a lower SHG than TOP by 10–20 W m−2 (Figure 3(g)–3(i)). Similar to the SABV, there is no notable difference between the two simulations regarding the SHV (Figure 3(j)–3(l)).
Figure 3

Averaged solar radiation absorbed by ground (SABG; in W m−2; a–c); solar radiation absorbed by ground (SABV; in W m−2; d–f); sensible heat flux from ground (SHG; in W m−2; g–i); and sensible heat flux from vegetation (SHV; in W m−2; j–l) over the period 2000–2010. For each row, TOP simulation is on the left, VIC simulation is in the middle, and the difference between the two simulations on the right (VIC minus TOP; Diff). Significant changes are indicated in black dots using the Student's t-test with α = 0.05.

Figure 3

Averaged solar radiation absorbed by ground (SABG; in W m−2; a–c); solar radiation absorbed by ground (SABV; in W m−2; d–f); sensible heat flux from ground (SHG; in W m−2; g–i); and sensible heat flux from vegetation (SHV; in W m−2; j–l) over the period 2000–2010. For each row, TOP simulation is on the left, VIC simulation is in the middle, and the difference between the two simulations on the right (VIC minus TOP; Diff). Significant changes are indicated in black dots using the Student's t-test with α = 0.05.

Close modal

Influence on surface energy balance

In this section, the potential influence of soil moisture changes in simulating the LE and SHF with respect to the ERA5 for the austral seasons: autumn (March–April–May; MAM), winter (June–July–August; JJA), spring (September–October–November; SON), and summer (December–January–February; DJF) is investigated. Figure 4 shows the seasonal climatology of the simulated LE (in W m−2) in comparison with the ERA5. From Figure 4, it can be observed that RegCM4 can reproduce the spatial pattern of the LE in all austral seasons. However, the RegCM4 bias, as well as the difference between the two simulations, vary with the season and region. For instance, the two simulations show a positive bias of 10–40 W m−2 overall in the Australia continent except for Queensland which exhibits a negative bias of 10–45 W m−2 in the austral MAM season (Figure 4(a)–4(e)).
Figure 4

Averaged latent heat flux (LE) over the period 2000–2010 (in W m−2) for each austral season (MAM a–f; JJA g–l; SON m–r; and DJF s–x). For each row, TOP is the first panel on the left, VIC is the second panel from the left, and ERA5 reanalysis product is the third from left, then the fourth (TOP minus ERA5) and fifth (VIC minus ERA5) panels are the model bias for each simulation against the ERA5 data and the difference between the two simulations (Diff). Significant changes are indicated in black dots using the Student's t-test with α = 0.05.

Figure 4

Averaged latent heat flux (LE) over the period 2000–2010 (in W m−2) for each austral season (MAM a–f; JJA g–l; SON m–r; and DJF s–x). For each row, TOP is the first panel on the left, VIC is the second panel from the left, and ERA5 reanalysis product is the third from left, then the fourth (TOP minus ERA5) and fifth (VIC minus ERA5) panels are the model bias for each simulation against the ERA5 data and the difference between the two simulations (Diff). Significant changes are indicated in black dots using the Student's t-test with α = 0.05.

Close modal

Furthermore, it can be observed that the VIC has a higher LE than the TOP by 5–15 W m−2 over the Northern Territory and it is small elsewhere. In the austral JJA and SON seasons, both simulations exhibit a negative bias of 10–25 W m−2 over all the continent except for Western Australia which has a positive bias of 20–30 W m−2 (Figure 4(g)–4(k) and 4(m)–4(q)). Also, the VIC differs from the TOP by 5–15 W m−2 over the same region in the austral JJA season (Figure 4(l)) and by 5–10 W m−2 in the lower region of Western Australia in the austral SON season (Figure 4(r)).

In addition, both simulations exhibit a positive bias of 15–50 W m−2 over the majority regions and by ∼60 W m−2 over Western Australia (Figure 4(s)–4(w)). Quantitatively, VIC has a higher LE than TOP by 5–25 W m−2 over the Continent except for the Western region (Figure 4(x)). It can be noticed that the influence of soil moisture changes depends on the region of interest as well as the austral season. Also, the noted positive bias of the LE (in the austral MAM and DJF seasons) can be attributed to the fact that the VIC has higher soil evaporation than the TOP (Section 3.1).

Figure 5 shows the seasonal climatology of the simulated SHF (in W m−2) in comparison with the ERA5. Similar to the LE, the RegCM4 is able to reproduce the spatial pattern of the SHF in comparison with the ERA5 in all austral seasons (Figure 5(a)–5(c), 5(g)–5(i), 5(m)–5(o), and 5(s)–5(u)). In the austral MAM season, the two simulations exhibit a negative bias of 15–20 W m−2 in the regions of Western Australia and a positive bias of 15–40 W m−2 elsewhere (Figure 5(d) and 5(e)). Quantitatively, the difference between the two simulations themselves is around ±10 W m−2 (Figure 5(f)). While in the austral JJA season, the RegCM4 model shows a positive bias of 15–30 W m−2 mainly over Northern Territory and Queensland (Figure 5(j) and 5(k)) and a negative bias of 5–10 W m−2 between the two simulations themselves (Figure 5(l)). Furthermore, the two simulations exhibit a positive bias of 15–40 W m−2 over Northern Territory, Queensland, New South Wales, and Victoria as well as the coast of Western Australia in the austral SON season (Figure 5(p) and 5(q)). Also, there is no significant difference between the two simulations themselves (Figure 5(r)). Lastly, in the austral DJF season, the two simulations show a negative bias of 20–60 W m−2 over Western Australia, Northern Territory, and Queensland (Figure 5(v) and 5(w)).
Figure 5

Averaged sensible heat flux (SHF) over the period 2000–2010 (in W m−2) for each austral season (MAM a–f; JJA g–l; SON m–r; and DJF s–x). For each row, TOP is the first panel on the left, VIC is the second panel from the left, and ERA5 reanalysis product is the third from left, then the fourth (TOP minus ERA5) and fifth (VIC minus ERA5) panels are the model bias for each simulation against the ERA5 data and the difference between the two simulations (Diff). Significant changes are indicated in black dots using the Student's t-test with α = 0.05.

Figure 5

Averaged sensible heat flux (SHF) over the period 2000–2010 (in W m−2) for each austral season (MAM a–f; JJA g–l; SON m–r; and DJF s–x). For each row, TOP is the first panel on the left, VIC is the second panel from the left, and ERA5 reanalysis product is the third from left, then the fourth (TOP minus ERA5) and fifth (VIC minus ERA5) panels are the model bias for each simulation against the ERA5 data and the difference between the two simulations (Diff). Significant changes are indicated in black dots using the Student's t-test with α = 0.05.

Close modal

This negative bias can be attributed to the fact that the RegCM4 shows the same bias order of magnitude with opposite signs in simulating the LE. In addition, VIC has a lower SHF than TOP by 10–15 W m−2 (Figure 5(x)). From Figures 4 and 5, it can be observed that there is no obvious difference between the two simulations (except for some regions), suggesting that soil moisture changes do not determine the surface energy budget. However, there are some differences between the two simulations themselves; which vary with the region under investigation as well as the austral season. In comparison with the ERA5-land, the RegCM4 shows a bias similar to the one observed using ERA5 (Supplementary Figure S2). Such a comparison suggests that the spatial resolution of the observed data does not considerably affect the model bias. Therefore, observed data with quite coarse resolution can be used to evaluate the model performance.

Influence on TCC and surface net radiation

In the present section, the possible influence of soil moisture changes (on the surface net radiation and TCC) is examined. In general, the RegCM4 can reproduce the spatial pattern of the TCC (in %) in comparison with the ERA5 in all austral seasons (Figure 6(a)–6(c), 6(g)–6(i), 6(m)–6(o), and 6(s)–6(u)). However, the RegCM4 bias varies with the region as well as the austral season. For instance, in the austral MAM season, both simulations show a negative bias of 30–40% in the coastal regions of Northern Territory, Queensland, and Victoria; while in the inner regions, the negative bias ranges from 10 to 20% (Figure 6(d) and 6(e)). In addition, there is no significant difference between the two simulations themselves (Figure 6(f)). In the austral JJA and SON seasons, the RegCM4 shows a negative bias of 30–40% over the coastal regions of Western Australia, Northern Territory, Queensland, and Victoria (Figure 6(j), 6(k), 6(p), and 6(q)). Qualitatively, the VIC shows a higher TCC than the TOP by ∼5%.
Figure 6

Averaged total cloud cover (TCC) over the period 2000–2010 (in %) for each austral season (MAM a–f; JJA g–l; SON m–r; and DJF s–x). For each row, TOP is the first panel on the left, VIC is the second panel from the left, and ERA5 reanalysis product is the third from left, then the fourth (TOP minus ERA5) and fifth (VIC minus ERA5) panels are the model bias for each simulation against the ERA5 data and the difference between the two simulations (Diff). Significant changes are indicated in black dots using the Student's t-test with α = 0.05.

Figure 6

Averaged total cloud cover (TCC) over the period 2000–2010 (in %) for each austral season (MAM a–f; JJA g–l; SON m–r; and DJF s–x). For each row, TOP is the first panel on the left, VIC is the second panel from the left, and ERA5 reanalysis product is the third from left, then the fourth (TOP minus ERA5) and fifth (VIC minus ERA5) panels are the model bias for each simulation against the ERA5 data and the difference between the two simulations (Diff). Significant changes are indicated in black dots using the Student's t-test with α = 0.05.

Close modal

In the DJF austral season, the RegCM4 model bias is restricted to the coastal regions of the Northern Territory and Victoria. Furthermore, the negative bias noted in the MAM, JJA, and SON austral seasons is quite shrunk (Figure 6(v) and 6(w)). Moreover, there is no notable difference between the two simulations themselves (Figure 6(x)). The quite small difference of the two simulations either with respect to the ERA5 or between themselves can be attributed to the small difference between the two simulations regarding the RH. Furthermore, the negative bias noted in all austral seasons, particularly those over the coastal regions, can be attributed to the convection parameterization of the RegCM4.

Figure 7 shows the simulated seasonal climatology of the surface net radiation (Rn; in W m−2) in comparison with the ERA5. It can be noticed that the RegCM4 captures well the spatial pattern of the Rn with respect to the ERA5 in all austral seasons, particularly the high values located in the region of 130–150°E and 15–35°S (Figure 7(a)–7(c), 7(g)–7(i), 7(m)–7(o), and 7(s)–7(u)). However, the RegCM4 model bias varies with the region under investigation as well as the season. In the MAM austral season, both simulations show a positive bias range from 20 to 80 W m−2 particularly over the lower boundary of Western Australia, Northern Territory, New South Wales, and South Australia (Figure 7(d) and 7(e)). Quantitatively, the VIC has a higher Rn than the TOP by 5 W m−2 (Figure 7(f)). In the austral JJA season, the situation is quite different because the two simulations show a positive bias of 20–40 W m−2 over all the continent except for Queensland and New South Wales where the Rn bias reaches 60 W m−2 (Figure 7(j) and 7(k)). In addition, there is no considerable difference between the two simulations (Figure 7(l)).
Figure 7

Averaged surface net radiation (Rn) over the period 2000–2010 (in W m−2) for each austral season (MAM a–f; JJA g–l; SON m–r; and DJF s–x). For each row, TOP is the first panel on the left, VIC is the second panel from the left, and ERA5 reanalysis product is the third from left, then the fourth (TOP minus ERA5) and fifth (VIC minus ERA5) panels are the model bias for each simulation against the ERA5 data and the difference between the two simulations (Diff). Significant changes are indicated in black dots using the Student's t-test with α = 0.05.

Figure 7

Averaged surface net radiation (Rn) over the period 2000–2010 (in W m−2) for each austral season (MAM a–f; JJA g–l; SON m–r; and DJF s–x). For each row, TOP is the first panel on the left, VIC is the second panel from the left, and ERA5 reanalysis product is the third from left, then the fourth (TOP minus ERA5) and fifth (VIC minus ERA5) panels are the model bias for each simulation against the ERA5 data and the difference between the two simulations (Diff). Significant changes are indicated in black dots using the Student's t-test with α = 0.05.

Close modal

In the austral SON season, the RegCM4 model bias ranges from 20 to 40 W m−2 over the majority of the continent except for the coastal region of Western Australia where the bias ranges from 60 to 80 W m−2 approaching 80 to 100 W m−2 over Victoria (Figure 7(p) and 7(q)). Furthermore, there is no significant difference between the two simulations themselves (Figure 7(r)). Finally, in the austral DJF season, the RegCM4 model bias ranges from +20 to 60 W m−2 approaching 80 W m−2 particularly over Victoria, New South Wales, South Australia, and Western Australia (Figure 7(v) and 7(w)). It is worth mentioning that the noted bias of the simulated LE, SHF, Rn, and TCC is not only because of the uncertainty associated with the RegCM4 physical parameterization, but also with the uncertainty associated with the ERA5 because ERA5 uses the TESSEL land surface model (Anwar et al. 2022). It can be noted that, the VIC has a higher Rn than the TOP by ∼10 W m−2 over Western Australia, Queensland, and Northern Territory (Figure 7(r)).

Furthermore, it can be observed that the two simulations show a high Rn because of the low TCC over all the continent suggesting that Rn is mainly constrained by the convection scheme rather than the soil moisture changes. Also, there is no considerable difference between the two simulations because the vegetation status is prescribed in both simulations (i.e., vegetation parameters retrieved from the MODIS; Section 2.2); therefore, there is no change in simulated total albedo and eventually Rn. It is important to mention that Equation (4) has been examined to check whether the surface net radiation equals the sum of the sensible, latent heat fluxes, and the heat storage in soil. After some mathematical calculations, it was found that the surface net radiation is larger than the sum of the sensible and latent heat fluxes by ∼150 W m−2 suggesting that the remaining part can be attributed to heat storage in soil. Also, shortwave/longwave parameterization is an important source of uncertainty in the simulated Rn.

Influence on the maximum and minimum air temperatures

The possible influence of soil moisture changes on the maximum and minimum air temperatures is examined by comparing the two simulations with respect to the Century/AGCD reanalysis dataset as well as between the two simulations themselves. Figure 8 shows the seasonal climatology of the maximum air temperature (TMX) simulated by the RegCM4 model with respect to the Century reanalysis dataset. In general, the RegCM4 can reproduce the spatial pattern of the TMX with respect to the Century in all austral seasons (Figure 8(a)–8(c), 8(g)–8(i), 8(m)–8(o), and 8(s)–8(u)). Also, the RegCM4 model bias varies with the season as well as the region of interest. For instance, in the austral MAM season, the two simulations show a cold bias of 1–2 °C over the coast of Northern Territory and Western Australia and a warm bias of 2–4 °C over Western Australia, South Australia, and Northern Territory (Figure 8(d), 8(e), 8(j), and 8(k)). Between the two simulations themselves, there is no significant difference (Figure 8(f)). In the austral JJA season, the RegCM4 has a warm bias of 1–2 °C approaching 4 °C in some regions (Figure 8(j) and 8(k)).
Figure 8

Averaged 2-m maximum air temperature (TMX) over the period 2000–2010 (in °C) for each austral season (MAM a–f; JJA g–l; SON m–r; and DJF s–x). For each row, TOP is the first panel on the left, VIC is the second panel from the left, and Century reanalysis product is the third from left, then the fourth (TOP minus Century) and fifth (VIC minus Century) panels are the model bias for each simulation against the Century data and the difference between the two simulations (Diff). Significant changes are indicated in black dots using the Student's t-test with α = 0.05.

Figure 8

Averaged 2-m maximum air temperature (TMX) over the period 2000–2010 (in °C) for each austral season (MAM a–f; JJA g–l; SON m–r; and DJF s–x). For each row, TOP is the first panel on the left, VIC is the second panel from the left, and Century reanalysis product is the third from left, then the fourth (TOP minus Century) and fifth (VIC minus Century) panels are the model bias for each simulation against the Century data and the difference between the two simulations (Diff). Significant changes are indicated in black dots using the Student's t-test with α = 0.05.

Close modal

Furthermore, the difference between the two simulations is less than 1 °C (Figure 8(l)). Furthermore, the two simulations exhibit a warm bias of 1–2 °C over all the continent and 3 °C over Western Australia and the Northern Territory (Figure 8(p) and 8(q)). Similar to the austral MAM and JJA seasons, there is no clear difference between the two simulations themselves (Figure 8(r)). It can be noticed that the RegCM4 behavior in the austral DJF is quite different. For instance, the two simulations show a cold bias of 1–2 °C approaching 4–5 °C over the borders of Western Australia and Northern Territory (Figure 8(v) and 8(w)). This cold bias can be attributed to the fact that the RegCM4 overestimates (underestimates) the latent (sensible) heat flux in this austral season (see Figures 4 and 5) because of the soil moisture surplus in this region (Sharmila & Hendon 2020). Moreover, the TMX cold/warm bias can be attributed to biases of the TCC, Rn, and SHF because low TCC leads to high Rn, SHF and therefore high TMX.

To account for the uncertainty associated with the observed dataset, the TOP simulation was taken as an example to be compared with Century and AGCD reanalysis products (Supplementary Figure S3). From Supplementary Figure S3, it can be noted that TOP can reproduce the spatial pattern of the TMX with respect to Century and AGCD (Supplementary Figure S3(a)–S3(c), S3(g)–S3(i), S3(m)–S3(o), and S3(s)–S3(u)). Furthermore, it can be noted that Century is colder than AGCD by 1–3 °C in all austral seasons except for the DJF, where Century is warmer than AGCD by 1–4 °C particularly over Western Australia (Supplementary Figure S3(f), S3(l), S3(r), and S3(x)). In the austral MAM, it can be observed that TOP shows a lesser bias (by 1–2 °C) when it is compared with AGCD than Century (Supplementary Figure S3(d) and S3(e)). Also, the TOP shows a lower bias (by 2 °C) when it is compared with AGCD relative to Century in the austral JJA and SON (Supplementary Figure S3(j), S3(k), S3(p), and S3(q)). Lastly in the austral DJF, the TOP shows quite similar bias for both Century and AGCD except for Western Australia where the bias with AGCD is lower than Century by 1–2 °C (Supplementary Figure S3(v) and S3(w)).

Figure 9 shows the simulated seasonal climatology of the daily minimum air temperature (TMN) with respect to the Century product. It can be observed that the RegCM4 can reproduce well the spatial pattern of the TMN in all austral seasons (Figure 9(a)–9(c), 9(g)–9(i), 9(m)–9(o), and 9(s)–9(u)). However, the RegCM4 bias varies with the region as well as the season. For instance, the two simulations show a warm bias of 1–4 °C mostly over Western Australia, South Australia, New South Wales, and Victoria with a cold bias of 1–2 °C over the Northern part of Western Australia in the austral MAM season (Figure 9(d) and 9(e)). Quantitatively, the difference between the two simulations themselves is less than 1 °C (Figure 9(f)). In the austral JJA season, the situation is quite different because the warm bias (1–3 °C) is restricted to the upper region of Western Australia and the Northern Territory, while over the rest of the continent the bias is ±1 °C (Figure 9(j) and 9(k)). In addition, there is no considerable difference between the two simulations themselves (Figure 9(l)). Moreover, the RegCM4 shows a warm bias of 1–5 °C over all the continent in the austral SON season (Figure 9(p) and 9(q)); while there is no difference between the two simulations (Figure 9(r)). Finally, in the austral DJF season, the two simulations show a cold bias of 1–5 °C over Western Australia with a warm bias of 2–4 °C over South Australia, New South Wales, and Victoria (Figure 9(v) and 9(w)). Like the other seasons, there is no notable difference between the two simulations (Figure 9(x)). The noted warm/cold bias of the TMN can be attributed to biases of TCC and Rn in all austral seasons, since the low TCC contributes to the increase of Rn and eventually increase of TMN.
Figure 9

Averaged 2-m minimum air temperature (TMN) over the period 2000–2010 (in °C) for each austral season (MAM a–f; JJA g–l; SON m–r; and DJF s–x). For each row, TOP is the first panel on the left, VIC is the second panel from the left, and Century reanalysis product is the third from left, then the fourth (TOP minus Century) and fifth (VIC minus Century) panels are the model bias for each simulation against the Century data and the difference between the two simulations (Diff). Significant changes are indicated in black dots using the Student's t-test with α = 0.05.

Figure 9

Averaged 2-m minimum air temperature (TMN) over the period 2000–2010 (in °C) for each austral season (MAM a–f; JJA g–l; SON m–r; and DJF s–x). For each row, TOP is the first panel on the left, VIC is the second panel from the left, and Century reanalysis product is the third from left, then the fourth (TOP minus Century) and fifth (VIC minus Century) panels are the model bias for each simulation against the Century data and the difference between the two simulations (Diff). Significant changes are indicated in black dots using the Student's t-test with α = 0.05.

Close modal

Similar to the TMX, the simulated TMN of the TOP was compared with the Century and AGCD reanalysis products (Supplementary Figure S4). From Supplementary Figure S4, it can be noted that Century is warmer than AGCD by 1–6 °C, particularly in the austral JJA and SON (Supplementary Figure S4(f), S4(l), S4(r), and S4(x)). In the austral MAM, the TOP shows a bias of 1–4 °C in comparison with Century particularly over South Australia and New South Wales (Supplementary Figure S4(d)), while TOP shows a high warm bias 1–6 °C over the same regions when it is compared with AGCD (Supplementary Figure S4(e)). Compared with Century, the TOP shows a bias up to 1–2 °C and 2–4 °C in the austral JJA and SON, respectively (Supplementary Figure S4(j) and S4(p)), while in comparison with the AGCD, the TOP shows a high warm bias (3–8 °C) over the Australian continent (Supplementary Figure S4(k) and S4(q)). In the austral DJF, the TOP shows a warm bias of 1–3 °C over a majority of the Australian continent and a cold bias of 1–3 °C over Western Australia in comparison with the Century (Supplementary Figure S4(v)). Moreover, the TOP shows a warm bias of 1–4 °C over a majority of Australia's continent and up to 6 °C over Queensland and New South Wales (Supplementary Figure S4(w)).

The ability of the RegCM4 model to capture the spatial pattern of the LE, SHF, Rn, TCC, and TMP was quantitatively evaluated by plotting a spatial map of the Pearson Correlation Coefficient (PCC) with the ERA5/Century reanalysis products (Supplementary Figure S5). From Supplementary Figure S5, it can be observed that both simulations show a gradient of PCC in simulating the LE: 0.4–0.6 over the interior region of Western Australia, 0.6–0.8 over the upper region of Western Australia, South Australia, and New South Wales as well as 0.8–1 over Queensland and Northern Territory (Supplementary Figure S5(a) and S5(f)). Regarding SHF, the RegCM4 shows also a gradient of PCC ranging from 0.4 to 0.6 over the upper half of Western Australia, Northern Territory, and Queensland; 0.6–1 over lower half of Western Australia, South Australia, and New South Wales (Supplementary Figure S5(b) and S5(g)). Moreover, both simulations show a PCC of 0.8–1 overall the continent regarding the Rn and TMP (Supplementary Figure S5(c), S5(h), S5(e), and S5(j)). For the TCC, both simulations show a gradient of PCC similar to the one observed in the simulated LE (Supplementary Figure S5(d) and S5(i)).

Evaluating the RegCM4 performance with in-situ observation

As reported in Section 3.4, there is no considerable difference between the two simulations regarding TMX and TMN. Therefore, the TOP simulation was taken as an example to evaluate the simulated mean air temperature (TMP) with respect to the in-situ observation. Evaluation of the TOP simulation (hereafter RegCM4) was done in two routes: (1) quantitative evaluation of the RegCM4 on a monthly basis using two statistical metrics (mean bias: MB and Pearson Correlation Coefficient: PCC) and (2) plotting the climatological annual cycle of the simulated TMP versus in-situ observation.

Figure 10 shows the climatological annual cycle of the TMP in comparison to the in-situ observation. At Adelaide River, the RegCM4 can reproduce the climatological seasonal cycle of the TMP in comparison with the observation. However, the RegCM4 underestimates the TMP from July until November with an MB of −1.026 °C and PCC = 0.768 (Figure 10(a)). Over Alice Springs, the RegCM4 reproduces well the seasonal cycle and overestimates the TMP with respect to the observation. Quantitatively, this is confirmed as the RegCM4 shows an MB of 3.9 °C and PCC of 0.972 (Figure 10(b)). At Howard Springs, the RegCM4 shows performance similar to the one observed over the Adelaide River with an MB of −1.16 °C and PCC of 0.768 (Figure 10(c)). In addition, the RegCM4 overestimates the TMP by 4.09 °C when it is compared with the observation over Robson Creek (Figure 10(d)).
Figure 10

Climatological seasonal cycle of the 2-m mean air temperature (in °C) as simulated by the RegCM4 model (in blue) with respect to station observation retrieved from the Ozflux network (in red) over the period 2000–2010. (a) Adelaide River, (b) Alice Springs, (c) Howard Springs, (d) Robson Creek, (e) Yanco, and (f) Tumbarumba. Please refer to the online version of this paper to see this figure in colour: https://dx.doi.org/10.2166/wcc.2023.512.

Figure 10

Climatological seasonal cycle of the 2-m mean air temperature (in °C) as simulated by the RegCM4 model (in blue) with respect to station observation retrieved from the Ozflux network (in red) over the period 2000–2010. (a) Adelaide River, (b) Alice Springs, (c) Howard Springs, (d) Robson Creek, (e) Yanco, and (f) Tumbarumba. Please refer to the online version of this paper to see this figure in colour: https://dx.doi.org/10.2166/wcc.2023.512.

Close modal

Moreover, the RegCM4 shows a smoother annual cycle than the observation (indicated by PCC of 0.847). Over Yanco and Tumbarumba, the RegCM4 overestimates the TMP by 2.76 and 5.46 °C, respectively. However, the RegCM4 shows a smoother seasonal cycle than the one noted in the observation and it has PCC of 0.9 over both locations (Figure 10(e) and 10(f)). It can be noticed that the RegCM4 performance varies with both location and time. Furthermore, the noted bias can be attributed to the horizontal resolution adopted in this study (i.e., 60 km) as well as the RegCM4 physical parameterization. This can be applied to any location considered in this study. Another reason is the mismatch in representing the land unit/vegetation between the RegCM4 and in-situ observation.

To explore the RegCM4 performance for capturing the monthly variability with respect to the in-situ observation, the Taylor diagram (Taylor 2001) was used to quantify such performance. Taylor's diagram is based on three statistical metrics: standard deviation ratio (STD), mean bias (MB), and spatial correlation coefficient (CORR). From Figure 11, it can be noted that all RegCM4 shows a high CORR for all stations ranging from 0.85 to 0.98. Also, the bias percentage is between 1 and 5%. However, the STD ratio varies considerably among all stations. For instance, Yanco and Tumbarumba show the least ratio compared with the other stations (around 0.35). The rest of the stations show an STD ranging from 0.63 to 0.82. Such performance of the RegCM4 can be attributed to the physical parameterization and the horizontal grid spacing (i.e., 60 km).
Figure 11

Taylor diagram for RegCM4 performance in capturing the monthly variability of TMP for the stations: Adelaide River, Alice Springs, Howard Springs, Robson Creek, Yanco, and Tumbarumba.

Figure 11

Taylor diagram for RegCM4 performance in capturing the monthly variability of TMP for the stations: Adelaide River, Alice Springs, Howard Springs, Robson Creek, Yanco, and Tumbarumba.

Close modal

Soil moisture plays an important role in controlling the surface energy balance and surface climate. Such a role has been documented in Llopart et al. (2017). In addition, the difference between land-surface hydrology schemes of the CLM4.5 land surface model was examined across the globe: Africa (Anwar et al. 2019), Amazon (Anwar et al. 2022), and Tibetan Plateau (Wang et al. 2021). Anwar et al. (2019, 2022) reported that the difference between the two land-surface hydrology schemes is considerable regarding the hydrological variables, surface energy balance, and surface climate. Furthermore, this difference varies with the region of study as well as the climate regime.

Anwar et al. (2019) found that the VIC scheme outperforms the TOP scheme in simulating the surface energy balance and climate. Also, Anwar et al. (2022) reported that there is no land-surface hydrology scheme that performs better than the other over the Amazon. Also, the role of RCMs in simulating the Australian surface climate is documented in Di Virgilio et al. (2019). Our main target is to examine the role of soil moisture changes in simulating the TMX and TMN of Australia using a regional climate model (RegCM4). To achieve this goal, two 13-year simulations were conducted and driven with the NCEP/NCAR reanalysis version 2. After that, the model output was evaluated with respect to various reanalysis products and in-situ observations. Also, the results of the present study were compared with those conducted over Africa and Amazon (Anwar et al. 2019, 2022).

As reported in Section 2.2, vegetation parameters (e.g., LAI) are retrieved from the MODIS. This means that the seasonal cycle of LAI is repeated every year and it is not affected by climate or soil moisture variability. Therefore, it can induce a source of uncertainty in the simulated surface energy balance, TMX and TMN. Moreover, the role of vegetation fraction changes (i.e., by enabling the dynamic vegetation; DV) is not considered in the present study. In comparison with the results reported by Anwar et al. (2019, 2022), the difference between the two land-surface hydrology schemes over Australia (as represented by QINFL and SM10) is quite small. Such a difference can be attributed for the reasons reported in Section 3.1. Therefore, a small difference is noted between the two simulations regarding the SHF and LE with respect to the ERA5 (Section 3.2). Also, there is no considerable difference between the two simulations regarding the Rn and TCC with respect to the ERA5. In addition, warm/cold biases can be attributed to the RegCM4 physical parameterization rather than soil moisture changes.

The regional climate model (RegCM4) was used to examine the role of the land-surface hydrology schemes (of the CLM4.5 land surface model) in controlling the daily maximum and minimum air temperatures of Australia with respect to reanalysis products and in-situ observations. The results are summarized in the following points:

  • 1.

    VIC has lower QINFL and SM10 than TOP. Therefore, VIC shows lower QVEGT and QVEGE with higher QSOIL than TOP. In addition, VIC has lower SABG and SHG than TOP.

  • 2.

    RegCM4 succeeds in reproducing the spatial pattern of SHF and LE with respect to the ERA5. However, the difference between the two simulations themselves depends on the region as well as the austral season.

  • 3.

    Soil moisture changes do not impose a notable impact on the simulated SHF and LE, Rn and TCC, TMX and TMN in comparison with the ERA5, Century and AGCD.

  • 4.

    RegCM4 can reproduce the climatological seasonal cycle of the mean air temperature (TMP) with respect to the in-situ observation.

  • 5.

    The simulated surface net radiation is larger than the sum of SHF and LE (for either TOP or VIC simulation).

Considering the recommendation of Anwar et al. (2022), the role of land-surface hydrology schemes need to be further examined across the globe. Our future study will include additional sensitivity studies for a long period (e.g., 30 years); revising the land-surface hydrology schemes TOP/VIC, addressing the role of different bias-correction methods to possibly reduce the surface energy fluxes and eventually the daily maximum and minimum air temperatures of Australia and finally considering the role of vegetation cover changes following Mehboob et al. (2020).

The Egyptian Meteorological Authority (EMA) is acknowledged for providing the computational power to conduct the RegCM4 model simulations. The climate group in ESP team is acknowledged for providing the RegCM4 model code and input data (http://clima-dods.ictp.it/RegCM4) to run the model. ECMWF is acknowledged for providing the ERA5/ERA5-land reanalysis product from their website at https://cds.climate.copernicus.eu.

20th Century Reanalysis V3 data was provided by the NOAA/OAR/ESRL PSL, Boulder, Colorado, USA, from their Web site at https://www.psl.noaa.gov/data/gridded/data.20thC_ReanV3.html.

The mean air temperature station data are retrieved from the Ozflux network. AGCD product is available at http://climate-cms.wikis.unsw.edu.au/AGCD. NCL plotting software is available at https://www.ncl.ucar.edu.

S.A.A.R. designed the simulation, analyzed the results, and wrote the manuscript. A.S. and B.Z. participated in analyzing the results, writing the manuscript. Furthermore, A.S. provided the observed mean air temperature data. The final version of the manuscript was revised by all authors.

The RegCM4 model code is available at https://github.com/ictp-esp/RegCM4.

We give our consent for the publication of identifiable details, which can include a photograph(s) and/or videos and/or case history and/or details within the text (‘Material’) to be published in the Journal of Water and Climate Change. We confirm that we have seen and been given the opportunity to read both the Material and the Article to be published by the Journal of Water and Climate Change.

All relevant data are included in the paper or its Supplementary Information.

The authors declare there is no conflict.

Allan
R. P.
,
Barlow
M.
,
Byrne
M. P.
,
Cherchi
A.
,
Douville
H.
,
Fowler
H. J.
,
Gan
T. Y.
,
Pendergrass
A. G.
,
Rosenfeld
D.
,
Swann
A. L. S.
,
Wilcox
L. J.
&
Zolina
O.
2020
Advances in understanding large-scale responses of the water cycle to climate change
.
Ann. N. Y. Acad. Sci.
1472
(
1
),
49
75
.
Anwar
S. A.
2019
Understanding the contribution of the vegetation-runoff system for simulating the African climate using the RegCM4 model
.
Theor. Appl. Climatol.
138
,
1219
1230
.
https://doi.org/10.1007/s00704-019-02885-x
.
Anwar
S. A.
&
Diallo
I.
2022
Modelling the Tropical African Climate using a state-of-the-art coupled regional climate-vegetation model
.
Clim. Dyn.
58
,
97
113
.
https://doi.org/10.1007/s00382-021-05892-9
.
Anwar
S. A.
,
Zakey
A. S.
,
Robaa
S. M.
&
Wahab
M. M.
2019
The influence of two land-surface hydrology schemes on the regional climate of Africa using the RegCM4 model
.
Theor. Appl. Climatol.
136
,
1535
.
https://doi.org/10.1007/s00704-018-2556-8
.
Anwar
S. A.
,
Mamadou
O.
,
Diallo
I.
&
Sylla
M. B.
2021
On the influence of vegetation cover changes and vegetation-runoff systems on the simulated summer potential evapotranspiration of tropical Africa using RegCM4
.
Earth Syst. Environ.
5
,
883
897
.
https://doi.org/10.1007/s41748-021-00252-3
.
Anwar
S. A.
,
Reboita
M. S.
&
Llopart
M.
2022
On the sensitivity of the Amazon surface climate to two land-surface hydrology schemes using a high-resolution regional climate model (RegCM4)
.
Int. J. Climatol.
42
(
4
),
2311
2327
.
https://doi.org/10.1002/joc.7367
.
Bonan
G. B.
,
Lawrence
P. J.
,
Oleson
K. W.
,
Levis
S.
,
Jung
M.
,
Reichstein
M.
,
Lawrence
D. M.
&
Swenson
S. C.
2011
Improving canopy processes in the Community Land Model (CLM4) using global flux fields empirically inferred from FLUXNET data
.
J. Geophys. Res.
116
,
G02014
.
https://doi.org/10.1029/2010JG001593
.
Bouraoui
F.
,
Grizzetti
B.
,
Granlund
K.
,
Rekolainen
S.
&
Bidoglio
G.
2004
Impact of climate change on the water cycle and nutrient losses in a Finnish catchment
.
Clim. Change
66
(
1
),
109
126
.
Chung
J. X.
,
Juneng
L.
,
Tangang
F.
&
Jamaluddina
A. F.
2018
Performances of BATS and CLM land-surface schemes in RegCM4 in simulating precipitation over CORDEX Southeast Asia domain
.
Int. J. Climatol.
38
,
794
810
.
https://doi.org/10.1002/joc.5211
.
Clough
S. A.
,
Shephard
M. W.
,
Mlawer
E. J.
&
Delamere
J. S.
2005
Atmospheric radiative transfer modeling: a summary of the AER codes, short communication
.
J. Quant. Spectrosc. Radiat. Transfer
91
,
233
244
.
Compo
G. P.
,
Whitaker
J. S.
,
Sardeshmukh
P. D.
,
Matsui
N.
,
Allan
R. J.
,
Yin
X.
,
Gleason
B. E.
,
Vose
R. S.
,
Rutledge
G.
,
Bessemoulin
P.
,
Brönnimann
S.
,
Brunet
M.
,
Crouthamel
R. I.
,
Grant
A. N.
,
Groisman
P. Y.
,
Jones
P. D.
,
Kruk
M.
,
Kruger
A. C.
,
Marshall
G. J.
,
Maugeri
M.
,
Mok
H. Y.
,
Nordli
Ø.
,
Ross
T. F.
,
Trigo
R. M.
,
Wang
X. L.
,
Woodruff
S. D.
&
Worley
S. J.
2011
The twentieth century reanalysis project
.
Q. J. R. Meteorol. Soc.
137
,
1
28
.
http://dx.doi.org/10.1002/qj.776
.
Coppola
E.
,
Giorgi
F.
,
Mariotti
L.
&
Bi
X.
2012
RegT-Band: a tropical band version of RegCM4
.
Clim. Res.
52
,
115
133
.
http://dx.doi.org/10.3354/cr01078
.
Coppola
E.
,
Stocchi
P.
,
Pichelli
E.
,
Torres Alavez
J. A.
,
Glazer
R.
,
Giuliani
G.
,
Di Sante
F.
,
Nogherotto
R.
&
Giorgi
F.
2021
Non-hydrostatic RegCM4 (RegCM4-NH): model description and case studies over multiple domains
.
Geosci. Modell. Dev.
14
,
7705
7723
.
https://doi.org/10.5194/gmd-14-7705-2021
.
Dickinson
R. E.
,
Henderson-Sellers
A.
&
Kennedy
P. J.
1993
Biosphere-Atmosphere Transfer Scheme (BATS) Version 1e as Coupled to the NCAR Community Climate Model (No. NCAR/TN-387 + STR)
.
University Corporation for Atmospheric Research
.
https://doi.org/10.5065/D67W6959
.
Di Virgilio
G.
,
Evans
J. P.
,
Di Luca
A.
,
Olson
R.
,
Argüeso
D.
,
Kala
J.
,
Andrys
J.
,
Hoffmann
P.
,
Katzfey
J. J.
&
Rockel
B.
2019
Evaluating reanalysis-driven CORDEX regional climate models over Australia: model performance and errors
.
Clim. Dyn.
53
,
2985
3005
.
https://doi.org/10.1007/s00382-019-04672-w
.
Giorgi
F.
,
Coppola
E.
,
Solmon
F.
,
Mariotti
L.
,
Sylla
M. B.
,
Bi
X.
,
Elguindi
N.
,
Diro
G. T.
,
Nair
V.
,
Giuliani
G.
,
Turuncoglu
U. U.
,
Cozzini
S.
,
Güttler
I.
,
O'Brien
T. A.
,
Tawfik
A. B.
,
Shalaby
A.
,
Zakey
A. S.
,
Steiner
A. L.
,
Stordal
F.
,
Sloan
L. C.
&
Brankovic
C.
2012
RegCM4: model description and preliminary tests over multiple CORDEX domains
.
Clim. Res.
52
,
7
29
.
Gu
H.
&
Wang
X.
2020
Performance of the RegCM4.6 for high-resolution climate and extreme simulations over Tibetan Plateau
.
Atmosphere
11
,
1104
.
https://doi.org/doi:10.3390/atmos11101104
.
Hersbach
H.
,
Bell
B.
,
Berrisford
P.
,
Hirahara
S.
,
Hor_anyi
A.
,
Muñoz-Sabater
J.
,
Nicolas
J.
,
Peubey
C.
,
Radu
R.
,
Schepers
D.
,
Simmons
A.
,
Soci
C.
,
Abdalla
S.
,
Abellan
X.
,
Balsamo
G.
,
Bechtold
P.
,
Biavati
G.
,
Bidlot
J.
,
Bonavita
M.
,
De Chiara
G.
,
Dahlgren
P.
,
Dee
D.
,
Diamantakis
M.
,
Dragani
R.
,
Flemming
J.
,
Forbes
R.
,
Fuentes
M.
,
Geer
A.
,
Haimberger
L.
,
Healy
S.
,
Hogan
R.J.
,
H_olm
E.
,
Janiskov_a
M.
,
Keeley
S.
,
Laloyaux
P.
,
Lopez
P.
,
Lupu
C.
,
Radnoti
G.
,
de Rosnay
P.
,
Rozum
I.
,
Vamborg
F.
,
Villaume
S.
&
Thépaut
J.-N.
2020
The ERA5 global reanalysis
.
Q. J. R. Meteorol. Soc.
146
,
1999
2049
.
Imbach
P.
,
Molina
L.
,
Locatelli
B.
,
Roupsard
O.
,
Mahe
G.
,
Neilson
R.
,
Corrales
L.
,
Scholze
M.
&
Ciais
P.
2012
Modeling potential equilibrium states of vegetation and terrestrial water cycle of Mesoamerica under climate change scenarios
.
J. Hydrometeorol.
13
(
2
),
665
680
.
Jones
D. A.
,
Wang
W.
&
Fawcett
R.
2009
High-quality spatial climate data-sets for Australia
.
Aust. Meteorol. Oceanogr. J.
58
,
233
248
.
Kanamitsu
M.
,
Ebisuzaki
W.
,
Woollen
J.
,
Yang
S. K.
,
Hnilo
J. J.
,
Fiorino
M.
&
Potter
G. L.
2002
NCEP-DOE AMIP-II reanalysis (R-2)
.
Bull. Am. Meteorol. Soc.
83
,
1631
1643
.
Kumari
N.
,
Srivastava
A.
,
Sahoo
B.
,
Raghuwanshi
N. S.
&
Bretreger
D.
2021
Identification of suitable hydrological models for streamflow assessment in the Kangsabati River Basin, India, by using different model selection scores
.
Nat. Resour. Res.
30
,
4187
4205
.
https://doi.org/10.1007/s11053-021-09919-0
.
Lawrence
P. J.
&
Chase
T. N.
2007
Representing a MODIS consistent land surface in the community land model (CLM 3.0)
.
J. Geophys. Res.
112
,
G01023
.
https://doi.org/10.1029/2006JG000168
.
Lawrence
D. M.
,
Fisher
R. A.
,
Koven
C. D.
,
Oleson
K. W.
,
Swenson
S. C.
,
Bonan
G.
,
Collier
N.
,
Ghimire
B.
,
van Kampenhout
L.
,
Kennedy
D.
,
Kluzek
E.
,
Lawrence
P. J.
,
Li
F.
,
Li
H.
,
Lombardozzi
D.
,
Riley
W. J.
,
Sacks
W. J.
,
Shi
M.
,
Vertenstein
M.
,
Wieder
W. R.
,
Xu
C.
,
Ali
A. A.
,
Badger
A. M.
,
Bisht
G.
,
van den Broeke
M.
,
Brunke
M. A.
,
Burns
S. P.
,
Buzan
J.
,
Clark
M.
,
Craig
A.
,
Dahlin
K.
,
Drewniak
B.
,
Fisher
J. B.
,
Flanner
M.
,
Fox
A. M.
,
Gentine
P.
,
Hoffman
F.
,
Keppel-Aleks
G.
,
Knox
R.
,
Kumar
S.
,
Lenaerts
J.
,
Leung
L. R.
,
Lipscomb
W. H.
,
Lu
Y.
,
Pandey
A.
,
Pelletier
J. D.
,
Perket
J.
,
Randerson
J. T.
,
Ricciuto
D. M.
,
Sanderson
B. M.
,
Slater
A.
,
Subin
Z. M.
,
Tang
J.
,
Thomas
R. Q.
,
Val
M. M.
&
Zeng
X.
2011
Parameterization improvements and functional and structural advances in version 4 of the community land model
.
J. Adv. Model. Earth Syst.
3
,
27
.
Lawrence
D. M.
,
Oleson
K. W.
,
Flanner
M. G.
,
Thornton
P. E.
,
Swenson
S. C.
,
Lawrence
P. J.
,
Zeng
X.
,
Yang
Z-L.
,
Levis
S.
,
Sakaguchi
K.
,
Bonan
G. B.
&
Slater
A. G.
2019
The community land model version 5: description of new features, benchmarking, and impact of forcing uncertainty
.
J. Adv. Model. Earth Syst.
https://doi.org/10.1029/2018MS001583
.
Lei
H.
,
Huang
M.
,
Leung
L. R.
,
Yang
D.
,
Shi
X.
,
Mao
J.
,
Hayes
D. J.
,
Schwalm
C. R.
,
Wei
Y.
&
Liu
S.
2014
Sensitivity of global terrestrial gross primary production to hydrologic states simulated by the Community Land Model using two runoff parameterizations
.
J. Adv. Model. Earth Syst.
6
,
658
679
.
https://doi.org/10.1002/2013MS000252
.
Levis
S.
,
Bonan
G. B.
,
Vertenstein
M.
&
Oleson
K. W.
2004
The Community Land Model's Dynamic Vegetation Model (CLM–DGVM): Technical Description and User's Guide. NCAR Technical Note TN-459 + IA
.
Liang
X.
,
Lettenmaier
D. P.
,
Wood
E. F.
&
Burges
S. J.
1994
A simple hydrologically based model of land surface water and energy fluxes for general circulation models
.
J. Geophys. Res.
99
,
14415
14428
.
Llopart
M.
,
Rosmeri
P.
,
da Rocha
M. R.
&
Santiago
C.
2017
Sensitivity of simulated South America climate to the land surface schemes in RegCM4
.
Clim. Dyn.
49
,
3975
3987
.
https://doi.org/10.1007/s00382-017-3557-5
.
Maurya
R. K. S.
,
Sinha
P.
,
Mohanty
M. R.
&
Mohanty
U. C.
2017
Coupling of community land model with RegCM4 for Indian summer monsoon simulation
.
Pure Appl. Geophys.
174
,
4251
4270
.
https://doi.org/10.1007/s00024-017-1641-8
.
Mehboob
M. S.
,
Kim
Y.
,
Lee
J.
,
Um
M. J.
,
Erfanian
A.
&
Wang
G.
2020
Projection of vegetation impacts on future droughts over West Africa using a coupled RegCM4-CLM-CN-DV
.
Clim. Change
163
(
2
).
https://doi.org/10.1007/s10584-020-02879-z
.
Muñoz-Sabater
J.
,
Dutra
E.
,
Agustí-Panareda
A.
,
Albergel
C.
,
Arduini
G.
,
Balsamo
G.
,
Boussetta
S.
,
Choulga
M.
,
Harrigan
S.
,
Hersbach
H.
,
Martens
B.
,
Miralles
D. G.
,
Piles
M.
,
Rodríguez-Fernández
N. J.
,
Zsoter
E.
,
Buontempo
C.
&
Thépaut
J. N.
2021
ERA5-Land: a state-of-the-art global reanalysis dataset for land applications
.
Earth Syst. Sci. Data
13
,
4349
4383
.
https://doi.org/10.5194/essd-13-4349-2021
.
Niu
G. Y.
,
Yang
Z. L.
,
Dickinson
R. E.
&
Gulden
L. E.
2005
A simple TOPMODEL-based runoff parameterization (SIMTOP) for use in global climate models
.
J. Geophys. Res.
110
,
D21106
.
https://doi.org/10.1029/2005JD006111
.
Oleson
K. W.
,
Niu
G.
,
Yang
Z. L.
,
Lawrence
D. M.
,
Thornton
P. E.
,
Lawrence
P. J.
,
Stöckli
R.
,
Dickinson
R. E.
,
Bonan
G. B.
,
Levis
S.
,
Dai
A.
&
Qian
T.
2008
Improvements to the community land model and their impact on the hydrologic cycle
.
J. Geophys. Res.
113
,
G01021
.
https://doi.org/10.1029/2007JD000563
.
Oleson
K. W.
,
Lawrence
D. M.
,
Bonan
G. B.
,
Drewniak
B.
,
Huang
M.
,
Koven
C. D.
,
Levis
S.
,
Li
F.
,
Riley
W. J.
,
Subin
Z. M.
,
Swenson
S.
,
Thornton
P. E.
,
Bozbiyik
A.
,
Fisher
R.
,
Heald
C. L.
,
Kluzek
E.
,
Lamarque
J.-F.
,
Lawrence
P. J.
,
Leung
L. R.
,
Lipscomb
W.
,
Muszala
S. P.
,
Ricciuto
D. M.
,
Sacks
W. J.
,
Sun
Y.
,
Tang
J.
&
Yang
Z.-L.
2010
Technical Description of Version 4.0 of the Community Land Model (CLM). NCAR Technical Note NCAR/TN-503 + STR. National Center for Atmospheric Research, Boulder
.
Oleson
K. W.
,
Lawrence
D. M.
,
Bonan
G. B.
,
Drewniak
B.
,
Huang
M.
,
Koven
C. D.
,
Levis
S.
,
Li
F.
,
Riley
W. J.
,
Subin
Z. M.
,
Swenson
S.
,
Thornton
P. E.
,
Bozbiyik
A.
,
Fisher
R.
,
Heald
C. L.
,
Kluzek
E.
,
Lamarque
J.-F.
,
Lawrence
P. J.
,
Leung
L. R.
,
Lipscomb
W.
,
Muszala
S. P.
,
Ricciuto
D. M.
,
Sacks
W. J.
,
Sun
Y.
,
Tang
J.
&
Yang
Z.-L.
2013
Technical Description of Version 4.5 of the Community Land Model (CLM). NCAR Technical Note NCAR/TN-503 + STR. National Center for Atmospheric Research, Boulder
.
Pal
J. S.
,
Giorgi
F.
,
Bi
X.
,
Elguindi
N.
,
Solmon
F.
,
Gao
X.
,
Rauscher
S. A.
,
Francisco
R.
,
Zakey
A.
,
Winter
J.
,
Ashfaq
M.
,
Syed
F. S.
,
Bell
J. L.
,
Diffenbaugh
N. S.
,
Karmacharya
J.
,
Konaré
A.
,
Martinez
D.
,
Da Rocha
R. P.
,
Sloan
L. C.
&
Steiner
A. L.
2007
Regional climate modeling for the developing world: the ICTP RegCM3 and RegCNET
.
Bull. Am. Meteorol. Soc.
88
,
1395
1409
.
https://doi.org/10.1175/BAMS-88-9-1395
.
Peel
M. C.
,
Finlayson
B. L.
&
McMahon
T. A.
2007
Updated world map of the Köppen-Geiger climate classification
.
Hydrol. Earth Syst. Sci.
11
,
1633
1644
.
https://doi.org/10.5194/hess-11-1633-2007
.
Peixoto
J. P.
&
Oort
A. H.
1992
Physics of Climate
.
AIP Press
,
New York
.
Reale
M.
,
Giorgi
F.
,
Solidoro
C.
,
Di Biagio
V.
,
Di Sante
F.
,
Mariotti
L.
,
Farneti
R.
&
Sannino
G.
2020
The Regional Earth System Model RegCM4-ES: evaluation of the Mediterranean climate and marine biogeochemistry
.
J. Adv. Model. Earth. Syst.
12
,
e2019MS001812
.
https://doi.org/10.1029/2019MS001812
.
Reynolds
R. W.
,
Rayner
N. A.
,
Smith
T. M.
,
Stokes
D. C.
&
Wang
W.
2002
An improved in situ and satellite SST analysis for climate
.
J. Clim.
15
,
1609
1625
.
Reynolds
J. F.
,
Smith
D. M. S.
,
Lambin
E. F.
,
Turner
B. L.
,
Mortimore
M.
,
Batterbury
S. P. J.
,
Downing
T. E.
,
Dowlatabadi
H.
,
Fernandez
R. J.
,
Herrick
J. E.
,
Huber-Sannwald
E.
,
Jiang
H.
,
Leemans
R.
,
Lynam
T.
,
Maestre
T.
,
Ayarza
M.
&
Walker
B.
2007
Global desertification: building a science for dryland development
.
Science
316
(
5826
),
847
851
.
https://doi.org/10.1126/science.1131634
.
Seneviratne
S. I.
,
Corti
T.
,
Davin
E. L.
,
Hirschi
M.
,
Jaeger
E. B.
,
Lehner
I.
,
Orlowsky
B.
&
Teuling
A. J.
2010
Investigating soil moisture–climate interactions in a changing climate: a review
.
Earth Sci. Rev.
99
(
3–4
),
125
161
.
https://doi.org/10.1016/j.earscirev.2010.02.004
.
Shalaby
A.
,
Zakey
A. S.
,
Tawfik
A. B.
,
Solmon
F.
,
Giorgi
F.
,
Stordal
F.
,
Sillman
S.
,
Zaveri
R. A.
&
Steiner
A. L.
2012
Implementation and evaluation of online gas-phase chemistry within a regional climate model (RegCM4-CHEM4)
.
Geosci. Model Dev.
5
,
741
760
.
Sharmila
S.
&
Hendon
H. H.
2020
Mechanisms of multiyear variations of Northern Australia wet-season rainfall
.
Sci. Rep.
10
,
508
.
https://doi.org/10.1038/s41598-020-61482-5
.
Slivinski
L. C.
,
Compo
G. P.
,
Whitaker
J. S.
,
Sardeshmukh
P. D.
,
Giese
B. S.
,
McColl
C.
,
Allan
R.
,
Yin
X.
,
Vose
R.
,
Titchner
H.
,
Kennedy
J.
,
Spencer
L. J.
,
Ashcroft
L.
,
Brönnimann
S.
,
Brunet
M.
,
Camuffo
D.
,
Cornes
R.
,
Cram
T. A.
,
Crouthamel
R.
,
Domínguez-Castro
F.
,
Freeman
J. E.
,
Gergis
J.
,
Hawkins
E.
,
Jones
P. D.
,
Jourdain
S.
,
Kaplan
A.
,
Kubota
H.
,
Le Blancq
F.
,
Lee
T.
,
Lorrey
A.
,
Luterbacher
J.
,
Maugeri
M.
,
Mock
C. J.
,
Moore
G. W. K.
,
Przybylak
R.
,
Pudmenzky
C.
,
Reason
C.
,
Slonosky
V. C.
,
Smith
C.
,
Tinz
B.
,
Trewin
B.
,
Valente
M. A.
,
Wang
X. L.
,
Wilkinson
C.
,
Wood
K.
&
Wyszyn'ski
P.
2019
Towards a more reliable historical reanalysis: improvements for version 3 of the Twentieth Century Reanalysis system
.
Q. J. R. Meteorol. Soc
.
https://doi.org/10.1002/qj.3598
.
Specht
R. L.
,
1970
Vegetation
. In:
Australian Environment
(
Leeper
G. W.
, ed.).
Melbourne University Press
,
Melbourne
, pp.
44
67
.
Srivastava
A.
,
Sahoo
B.
,
Raghuwanshi
N. S.
&
Singh
R.
2017
Evaluation of variable-infiltration capacity model and MODIS-terra satellite-derived grid-scale evapotranspiration estimates in a river basin with tropical monsoon-type climatology
.
J. Irrig. Drain. Eng.
143
(
8
),
04017028
.
https://doi.org/10.1061/(ASCE)IR.1943-4774.0001199
.
Srivastava
A.
,
Sahoo
B.
,
Raghuwanshi
N. S.
&
Chatterjee
C.
2018
Modelling the dynamics of evapotranspiration using variable infiltration capacity model and regionally calibrated Hargreaves approach
.
Irrig. Sci.
36
(
4
),
289
300
.
https://doi.org/10.1007/s00271-018-0583-y
.
Srivastava
A.
,
Kumari
N.
&
Maza
M.
2020
Hydrological response to agricultural land use heterogeneity using variable infiltration capacity model
.
Water Resour. Manage.
34
,
3779
3794
.
https://doi.org/10.1007/s11269-020-02630-4
.
Srivastava
A.
,
Rodriguez
J. F.
,
Saco
P. M.
,
Kumari
N.
&
Yetemen
O.
2021a
Global analysis of atmospheric transmissivity using cloud cover, aridity and flux network datasets
.
Remote Sens.
13
,
1716
.
https://doi.org/10.3390/rs13091716
.
Srivastava
A.
,
Saco
P. M.
,
Rodriguez
J. F.
,
Kumari
N.
,
Chun
K. P.
&
Yetemen
O.
2021b
The role of landscape morphology on soil moisture variability in semi-arid ecosystems
.
Hydrol. Processes
35
(
1
),
e13990
.
Steiner
A. L.
,
Pal
J. S.
,
Rauscher
S. A.
,
Bell
J. L.
,
Diffenbaugh
N. S.
,
Boone
A.
,
Sloan
L. C.
&
Giorgi
F.
2009
Land surface coupling in regional climate simulations of the West African monsoon
.
Clim. Dyn.
33
(
6
),
869
892
.
https://doi.org/10.1007/s00382-009-0543-6
.
Tarek
M.
,
Brissette
F. P.
&
Arsenault
R.
2020
Evaluation of the ERA5 reanalysis as a potential reference dataset for hydrological modelling over North America
.
Hydrol. Earth Syst. Sci.
24
,
2527
2544
.
https://doi.org/10.5194/hess-24-2527-2020
.
Wang
X.
,
Yang
M.
&
Pang
G.
2015
Influences of two land-surface schemes on RegCM4 precipitation simulations over the Tibetan Plateau
.
Adv. Meteorol.
2015
,
12
.
Article ID 106891, http://dx.doi.org/10.1155/2015/106891
.
Wang
X.
,
Chen
D.
,
Pang
G.
,
Anwar
S. A.
,
Ou
T.
&
Yang
M.
2021
Effects of cumulus parameterization and land-surface hydrology schemes on Tibetan Plateau climate simulation during the wet season: insights from the RegCM4 model
.
Clim. Dyn.
57
,
1853
1879
.
https://doi.org/10.1007/s00382-021-05781-1
.
Zhengqi
W.
,
Xuejie
G.
,
Zhenyu
H.
,
Jia
W.
,
Ying
X.
&
Liew
J.
2020
Assessing the sensitivity of RegCM4 to cumulus and ocean surface schemes over the Southeast Asia domain of the coordinated regional climate downscaling experiment
.
Atmos. Oceanic Sci. Lett.
13
(
1
),
71
79
.
https://doi.org/10.1080/16742834.2020.1697615
.
This is an Open Access article distributed under the terms of the Creative Commons Attribution Licence (CC BY 4.0), which permits copying, adaptation and redistribution, provided the original work is properly cited (http://creativecommons.org/licenses/by/4.0/).