Responses of streamflow to vegetation and climate change in southwestern Australia

Southwestern Australia has experienced recent climate change, with an increase in air temperature of 0.6°C and a reduction in mean annual precipitation of -15% since 1970. Along with the warming and drying trends, dramatic declines of streamflow have occurred across the region. However, both forest mortality and an increase in leaf area index have been observed in the southwestern forest, suggesting varied responses of vegetation to climate change. In this study, 30 catchments were analyzed using the Mann-Kendall trend test, Pettitt’s change point test and the theoretical framework of the Budyko curve to study the rainfall-runoff relationship change, and effects of climate and land cover change on streamflow. A declining trend and relatively consistent change point (2000) of streamflow were found in most catchments, with 14 catchments showing significant declines (p < 0.05, -68.1% to -35.6%) over 1970-2000 and 2001-2015. Most of the catchments have been shifting towards a more water-limited climate condition since 2000. For the period of 1970 to 2015, the dynamic of vegetation attributes (land cover/use change and growth of vegetation) dominated the decrease of streamflow in about half the study catchments. In general, a coequal role of climate and vegetation on the decline in streamflow was found in the study, suggesting the importance of vegetation management on future water management and production. Precipitation is predicted to decline in the future; therefore, some forest management intervention is required to maintain forest growth and water supply in the southwest of Australia

A Quantitative Inspection on Spatio-Temporal Variation of Remote Sensing-Based Estimates of Land Surface Evapotranspiration in South Asia

Evapotranspiration (ET) plays a key role in water resource management. It is important to understand the ET spatio-temporal pattern of South Asia for understanding and anticipating serious water resource shortages. In this study, daily ET in 2008 was estimated over South Asia by using MODerate Resolution Imaging Spectroradiometer (MODIS) products combined with field observations and Global Land Data Assimilation System (GLDAS) product through Surface Energy Balance System (SEBS) model. Monthly ET data were calculated based on daily ET and evaluated by the GLDAS ET data. Good agreements were found between two datasets for winter months (October to February) with R2 from 0.5 to 0.7. Spatio-temporal analysis of ET was conducted. Ten specific sites with different land cover types at typical climate regions were selected to analyze the ET temporal change pattern, and the result indicated that the semi-arid or arid areas in the northwest had the lowest average daily ET (around 0.3 mm) with a big fluctuation in the monsoon season, while the sites in the Indo-Gangetic Plain and in southern India has bigger daily ET (more than 3 mm) due to a large water supplement. It is suggested that the monsoon climate has a large impact on ET spatio-temporal variation in the whole region

Using Remote Sensing to Estimate Crop Water Use to Improve Irrigation Water Management

Irrigation water is scarce. Hence, accurate estimation of crop water use is necessary for proper irrigation managements and water conservation. Satellite-based remote sensing is a tool that can estimate crop water use efficiently. Several models have been developed to estimate crop water requirement or actual evapotranspiration (ETa) using remote sensing. One of them is the Mapping EvapoTranspiration at High Resolution using Internalized Calibration (METRIC) model. This model has been compared with other methods for ET estimations including weighing lysimeters, pan evaporation, Bowen Ratio Energy Balance System (BREBS), Eddy Covariance (EC), and sap flow. However, comparison of METRIC model outputs to an atmometer for ETa estimation has not yet been attempted in eastern South Dakota. The results showed a good relationship between ETa estimated by the METRIC model and estimated with atmometer (r2 = 0.87 and RMSE = 0.65 mm day-1). However, ETa values from atmometer were consistently lower than ETa values from METRIC. The verification of remotely sensed estimates of surface variables is essential for any remote-sensing study. The relationships between LAI, Ts, and ETa estimated using the remote sensing-based METRIC model and in-situ measurements were established. The results showed good agreement between the variables measured in situ and estimated by the METRIC model. LAI showed r2 = 0.76, and RMSE = 0.59 m2 m-2, Ts had r2 = 0.87 and RMSE 1.24 °C and ETa presented r2= 0.89 and RMSE = 0.71 mm day-1. Estimation of ETa using energy balance method can be challenging and time consuming. Thus, there is a need to develop a simple and fast method to estimate ETa using minimum input parameters. Two methods were used, namely 1) an energy balance method (EB method) that used input parameters of the Landsat image, weather data, a digital elevation map, and a land cover map and 2) a Kc-NDVI method that use two input parameters: the Landsat image and weather data. A strong relationship was found between the two methods with r2 of 0.97 and RMSE of 0.37 mm day-1. Hence, the Kc-NDVI method performed well for ETa estimations, indicating that Kc-NDVI method can be a robust and reliable method to estimate ETa in a short period of time. Estimation of crop evapotranspiration (ETc) using satellite remote sensing-based vegetation index such as the Normalized Difference Vegetation Index (NDVI). The NDVI was calculated using near-infrared and red wavebands. The relationship between NDVI and tabulated Kc’s was used to generate Kc maps. ETc maps were developed as an output of Kc maps multiplied by reference evapotranspiration (ETr). Daily ETc maps helped to explain the variability of crop water use during the growing season. Based on the results we can conclude that ETc maps developed from remotely sensed multispectral vegetation indices are a useful tool for quantifying crop water use at regional and field scales

Water availability and demand analysis in the Kabul River Basin, Afghanistan

Kabul River Basin (KRB), the most populated and highly heterogenic river basin of Afghanistan is the lifeline of millions of people in terms of supplying them with water for agricultural, municipal, and industrial as well as hydropower production purposes. Unfortunately, KRB is facing a multiplicity of governance, management and development relevant challenges for the last couple of decades. Detailed and reliable assessments of land use and land cover, water demand (for different sectors) as well as the available water resources are prerequisites for Integrated Water Resources Management across the basin. To achieve increased accuracy for water availability and demand analysis across the KRB, the study area was segregated into different hydrological and administrative units (provincial level, subbasin level etc.) in order to capture the heterogeneity driven by complex physiographic conditions (mainly due to huge elevation differences) and resulting in diverse cropping pattern at different reaches of the river basin. The innovative part of this study has been the concept of introducing spatial segregation of the large heterogenic river basin and using crop phenological information for evapotranspiration and land cover analysis respectively; it gave a distinct value to the output of this study. Phenologically tuned normalized difference vegetation indices (NDVI) of Aqua and Terra platforms with moderate resolution (250 m) proved to be very effective in the estimation of the land cover across the KRB with high accuracy. The phenology based segregated spatial analyses of the LULC of KRB with reference to 2003 (the base year of the study) highlighted the change in the ground coverage of main crops across the KRB e.g. wheat, barley, maize and rice. Based on the evaluation of the above results referring to the period 2003 to 2013, the rise in wheat ground coverage has been compensated by the decline in barley cultivation; maize and rice share has been almost consistent among the dominant cereals production in KRB. Upon spatial segregation, across the sub-basins (Alingar, Chak aw Logar, Ghorband aw Panjshir, Gomal, Kabul, Kunar and Shamal) Shamal, Kunar and Kabul showed highest actual evapotranspiration (ETa) throughout the study period of 2003 to 2013. The later three sub-basin host relatively large irrigated areas and production of two crops per year due to relatively favorable climatic and geographic conditions. Besides the agricultural water demand (ETa), water availability estimation through rainfall-runoff modelling by the use of the Soil and Water Assessment Tool (SWAT) has been very useful in data scarce regions like KRB. The application of the hydrological model using remote sensing products as input is the only effective choice in data scarce regions and exhibited results which are required by policy makers and investors for the strategic and sustainable planning and management of land and water resources

Changes in water and carbon in Australian vegetation in response to climate change

Australia has experienced pronounced climate change since 1950, especially in forested areas where a reducing trend in annual precipitation has occurred. However, the interaction between forests and water at multiple scales, in different geographical locations, under different management regimes and in different forest types with diverse species is not fully understood. Therefore, some interactions between forests and hydrological variables, and in particular whether the changes are mediated by management or climate, remain controversial. This thesis investigates the responses of Australia’s terrestrial ecosystems to both historical and projected climate change using remote sensing data and ecohydrological models. The thesis is structured in seven chapters, and contains five research chapters. Vegetation dynamics and sensitivity to precipitation change on the Australian continent for the past long drought period (2002-2010) are explored in Chapter 2 using multi-resource vegetation indices (VIs; normalized difference vegetation index (NDVI) and leaf area index (LAI)) and gridded climate data. During drought, precipitation and VIs declined across 90% and 80% of the whole continent, respectively, compared to the baseline period of 2000-2001. The most dramatic declines in VIs occurred in open shrublands near the centre of Australia and in southwestern Australia coinciding with significant reductions in precipitation and soil moisture. Overall, a strong relationship between water (precipitation and soil moisture) and VIs was detected in places where the decline in precipitation was severe. For five major vegetation types, cropland showed the highest sensitivity to water change, followed by grassland and woody savanna. Open shrublands showed moderate sensitivity to water change, while evergreen broadleaf forests only showed a slight sensitivity to soil moisture change. Although there was no consistent significant relationship between precipitation and VIs of evergreen broadleaf forests, forests in southeastern Australia, where precipitation had declined since 1997, appear to have become more sensitive to precipitation change than in southwestern Australia. The attribution of impacts from climate change and vegetation on streamflow change at the catchment scale for southwestern Australia are described in Chapter 3. This region is characterized by intensive warming and drying since 1970. Along with these significant climate changes, dramatic declines in streamflow have occurred across the region. Here, 79 catchments were analyzed using the Mann-Kendall trend test, Pettitt’s change point test, and the theoretical framework of the Budyko curve to study changes in the rainfall-runoff relationship, and effects of climate and vegetation change on streamflow. A declining trend and relatively consistent change point (2000) of streamflow were found in most catchments, with over 40 catchments showing significant declines (p < 0.05, -20% to -80%) between the two periods of 1982-2000 and 2001-2011. Most of the catchments have been shifting towards a more water-limited climate condition since 2000. Although streamflow is strongly related to precipitation for the period of 1982 to 2011, change of vegetation (land cover/use change and growth of vegetation) dominated the decrease in streamflow in about two-thirds of catchments. The contributions of precipitation, temperature and vegetation to streamflow change for each catchment varied with different catchment characters and climate conditions. In Chapter 4, the magnitude and trend of water use efficiency (WUE) of forest ecosystems in Australia, and their response to drought from 1982 to 2014, were analyzed using a modified version of the Community Atmosphere Biosphere Land Exchange (CABLE) land surface model in the BIOS2 modelling environment. Instead of solely relying on the ratio of gross primary productivity (GPP) to evapotranspiration (ET) as WUE (GPP/ET), the ratio of net primary productivity (NPP) to Transpiration (ETr) (NPP/ETr) was also adopted to more comprehensively understand the response of vegetation to drought. For the study period, national average annual forest WUE was 1.39 ± 0.80 g C kg−1 H2O for GPP/ET and 1.48 ± 0.28 g C kg−1 H2O for NPP/ETr. The WUE increased in the entire study area during this period (with a rate of 0.003 g C kg−1 H2O yr-1 for GPP/ET; p < 0.005 and a rate of 0.0035 g C kg−1 H2O yr-1 for NPP/ETr; p < 0.01), whereas different trends were detected in different biomes. A significantly increasing trend of annual WUE was only found in woodland areas due to higher magnitudes of increases in GPP and NPP than ET and ETr. The exception was in eucalyptus open forest area where ET and ETr decreased more than reductions in GPP and NPP. The response of WUE to drought was further analyzed using 1-48 month scales standardised precipitation-evapotranspiration index (SPEI). More severe (SPEI < -1) and frequent droughts (over ca. 8 years) occurred in the north than in the southwest and southeast of Australia since 1982. The response of WUE to drought varied significantly regionally and across forest types. The response of WUE to drought varied significantly regionally and across forest types, due to the different responses of carbon sequestration and water consumption to drought. The cumulative lagged effect of drought on monthly WUE derived from NPP/ETr was consistent and relatively short and stable between biomes (< 4 months), but notably varied for WUE based on GPP/ET, with a long time lag (mean of 16 months). As Chapters 2-4 confirmed that climate change has been playing an important role in the water yield and vegetation dynamics in Australia, the response of water yield and carbon sequestration to projected future climate change scenarios were integrated using the Water Supply Stress Index and Carbon model (WaSSI-C) ecohydrology model in Chapter 5. This model was calibrated with the latest water and carbon observations from the OzFlux network. The performance of the WaSSI-C model was assessed with measures of Q from 222 Hydrologic Reference Stations (HRSs) in Australia. Across the 222 HRSs, the WaSSI-C model generally captured the spatial variability of mean annual and monthly Q as evaluated by the Correlation Coefficient (R2 = 0.1-1.0), Nash-Sutcliffe Efficiency (NSE = -0.4-0.97), and normalized Root Mean Squared Error by Q (RMSE/Q = 0.01-2.2). Then 19 Global Climate Models (GCMs) from the Coupled Model Intercomparison Project phase 5 (CMIP5), across all Representative Concentration Pathways (RCPs) (RCP2.6, RCP4.5, RCP6.0 and RCP8.5), were used to investigate the potential impacts of climate change on water and carbon fluxes. Compared with the baseline period of 1995-2015 across the 222 HRSs, the temperature was projected to rise by an average of 0.56 to 2.49 ˚C by 2080, while annual precipitation was projected to vary significantly. All RCPs demonstrated a similar spatial pattern of change of projected Q and GPP by 2080, however, the magnitude varied widely among the 19 GCMs. Overall, future climate change may result in a significant reduction in Q but may be accompanied by an increase in ecosystem productivity. Mean annual Q was projected to decrease by 5 - 211 mm yr-1 (34% - 99%) by 2080, with over 90% of the watersheds declining. On the contrary, GPP was projected to increase by 17 - 255 g C m-2 yr-1 (2% - 17%) by 2080 in comparison with 1995-2015 in southeastern Australia. A significant limitation of WaSSI-C model is that it only runs serially. High resolution simulations at the continental scale are therefore not only computationally expensive but also present a run-time memory burden. In Chapter 6, using distributed (Message Passing Interface, MPI) and shared (Open Multi-Processing, OpenMP) memory parallelism techniques, the model was parallelized (and renamed as dWaSSI-C), and this approach was very effective in reducing the computing run-time and memory use. By using the parallelized model, several experiments were carried out to simulate water and carbon fluxes over the Australian continent to test the sensitivity of the model to input data-sets of different resolutions, as well as the sensitivity of the model to its WUE parameter for different vegetation types. These simulations were completed within minutes using dWaSSI-C, and this would not have been possible with the serial version. Results show that the model is able to simulate the seasonal cycle of GPP reasonably well when compared to observations at 4 eddy flux sites in Australia. The sensitivity analysis showed that simulated GPP was more sensitive to WUE during the Australian summer as compared to winter, and woody savannas and grasslands showed higher sensitivity than evergreen broadleaf forests and shrublands. With the parallelized dWaSSI-C model, it will now be much easier and faster to conduct continental scale analyses of the impacts of climate change and land cover change on water and carbon. Overall, vegetation and water of Australian ecosystems have become very sensitive to climate change after a considerable decline in streamflow. Australian ecosystems, especially in temperate Australia, are projected to experience warmer and drier climate conditions with increasing drought risk. However, the prediction from different models varied significantly due to the uncertainty of each climate model. The impacts of different forest management scenarios should be studied to find the best land use pattern under the changing climate. Forest management methods, such as thinning and reforestation, may be conducted to mitigate the impacts of drought on water yield and carbon sequestration in the future