Research Plan
Research Plan
Research Plan
INTRODUCTION
1.1. Background
Accurate measurement of rainfall on global and regional scale are critical for understanding the
climate and hydrological cycle, simulating land surface hydrological processes, water resources
management and security (Stage et al., 2017; Tan et al., 2016) and it is also an imperative input
for hydrological models (Habib et al., 2014).
Variations in rainfall often influenced by non-linear relations between numerous factors like local
variations of topography, the orientation of mountains and aspects (Haile et al., 2009) which clues
directly to flexible discharge output (Biemans, 2012) and fluctuations in water storage (Narjary
and Kamra, 2013). Finally, this leads to incorrect simulation results and sometimes even in wrong
conclusions (Vrgut et al., 2005). Precipitation data can be collected directly through ground-based
observations using rain gauges and weather radars (Zhang et al., 2014), and indirectly through
satellites (Lopez et al., 2017). Still, ground rainfall stations in many parts of the world and most of
parts of Ethiopia are very sparse and unevenly distributed (Ayehu et al., 2017).
Thus, rainfall estimation is unavoidably subject to error due to different factors depending on the
type of measurement (Alemohammed et al., 2015). To overcome some of the above-mentioned
restrictions of ground-based precipitation measurements, space-based (satellite) estimations of
precipitations provide a hopeful alternative source (Bigiarini et al., 2017). Satellite-based
precipitation estimates may offer plentiful information with spatio-temporal (near-real-time) high-
resolution over wide-ranging regions wherever conventional rainfall data are unusual or absent
(Kidd and Huffman, 2011; Kidd et al., 2009).
Since combined precipitation products mostly have more improved quality than individual data
set (Xie and Xiong, 2011), different emerged precipitation products (reanalysis-gauge and
satellite-gauge) are gradually being developed (Duan et al., 2017).It is principally worth nothing
that the lately released Climate Hazards Group Infra-Red Precipitation with Station data
(CHIRPS_2) provides rainfall at the finest spatial resolution of 0.05o and 0.25o (Funk et al., 2015).
At different times many scholars have been carried out a study to assess the performance of satellite
precipitation estimates for hydrological modeling in different regions of the world including
Ethiopia and indicated that high resolution satellite precipitation products have potential use for
hydrological modeling based on their product types.
Hence, this study was intended to assess the suitability of CHIRPS_2 and Tropical Rainfall
Measuring Mission (TRMM) Multisatellite Precipitation Analysis,TMPA_3B42v7 satellite
precipitation products besides in-situ precipitation products for stream flow simulation in Dabus
watershed, Upper Blue Nile basin, Ethiopia in a more comprehensive way in-terms of evaluated
aspects, temporal and spatial scales.
Subsequently measurement of all parameters that affect watershed’s stream flow simulation is
difficult, choosing a suitable model of simple structure, minimum input data requests and
reasonable precision is essential (Sharifi et al., 2004).
One of the hydrological models that meet the listed criteria around the globe is Hydrological
Engineering Center-Hydrologic Modeling System (HEC-HMS). Due to the above reasons and its
confirmed competence in the Upper Blue Nile basin (Fekadu et al., 2017; Onyutha, 2016; Gebre,
2015), this study used HEC-HMS hydrological model to evaluate the highly resolution satellite
precipitation products for stream flow simulations.
Spatial distribution and extent of rainfall are imperative for establishing a concrete flood
prediction, warning services and drought monitoring. But, in many regions of the world the number
of ground measuring stations is very scarce and randomly distributed which make countries to
have poor (humble) water resources assessment and flood prediction (Shrestha et al, 2011; Yong
et al., 2010).
In order to challenge the above doubts, a large number of academics have been exploring that
satellite rainfall assessments have potential use for hydrological modeling, but model performance
and uncertainty depends on way of the model developed, physical condition of application,
external factors (data), satellite product type, watershed size, hydro climatic regimes ( Liu, 2015;
Gebregiorgis and Hossain, 2013) and some behaviors of models are still difficult to totally
understand and can hide surprising properties (Berthet et al., 2010a).
Subsequently, all the above evaluations did not include Dabus watershed of Ethiopia which is
rugged in topography and have more sparsely in-situ rain gauge networks. Additionally, most of
the studies were limited to satellite precipitation products of early version and only event-based
analysis for examples, in (Tesfaye et al., 2017). Moreover, in Dabus watershed
the influences of water withdrawal, land cover and climate change have not been quantified yet
due to scarce rainfall data for rainfall-runoff modeling.
As a result, this study concentrates on evaluating satellite rainfall products of CHIRPS_2 and
TMPA_3B42v7 estimate by using in-situ rain gauge products based on their hydrologic simulation
result as well as identification of sensitive parameters of model out puts when calibrated and
validated with input specific satellite precipitation estimates. The HEC-HMS hydrologic model
was used to simulate the stream flow for both rain gauge and satellite rainfall products.
1.3. Objectives
The main objective of this research is to evaluate CHIRPS_2 and TMPA_3B42v7 satellite rainfall
products for hydrologic simulation in Dabus watershed, Abbay basin, Ethiopia.
To evaluate the ability of satellite rainfall products to characterize rainfall patterns and
capture the magnitude of rainfall over Dabus watershed.
To assess the skill of satellite rainfall products as an input into a hydrological model for
stream flow simulation.
To evaluate the performance of HEC-HMS model using satellite and ground based
precipitation products.
1.4. Research Questions
Even if, in-situ gauging stations exist in the study area the distribution is very inadequate, random
and less reliable. Therefore, to overcome such type of problems this study tried to evaluate high
spatial and temporal resolution of satellite rainfall products for hydrological simulation in order to
add confidence on reality of satellite products in place of inaccessible (data scarcity) area. The
study also may be helpful for hydrologists, climate scientists and researchers in real life (design)
applications mainly, for the study area.
This study is limited on evaluating CHIRPS_2 and TMPA_3B42v7 satellite precipitation products
by using in-situ rain gauge precipitation data’s with HEC-HMS hydrologic model for Dabus
watershed, Abbay basin, Ethiopia. Primarily, simulation of observed rainfall data was takes place;
secondly, simulation of satellite precipitation data’s was performed.
Finally, comparison of in-situ rainfall data’s with satellite rainfall products and comparison
between the two satellite precipitation products were achieved. This study doesn’t take in to
account stream flow changes caused due to model inputs other than rainfall.
This thesis work is organized in five chapters, references and appendices. The first chapter consist
of introduction, background, statement of the problem, objectives of the study, research questions,
significance and scope of the study. Chapter two review the work done earlier on related topics.
Chapter three describes about location and background of the study area, methodology, data source
and methods for data collection, analysis systems and conceptual framework.
Chapter four presents the outcomes and discussions based on the objectives of the study. Chapter
five contain the conclusions and recommendations of the study depending on the results, whereas
the list of references includes the references of different works cited in the manuscript. Finally, list
of the appendices in the form of tables and figures are presented as a supporting document to this
study.
2. LITERATURE REVIEW
2.1. Rainfall
Rainfall is a main component of the climate system and plays a critical role in the Earth’s
hydrological cycle and energy balance. Conventionally, the rain gauge is the main source of
precipitation data, which has been the most accurate and reliable approach for rainfall data
measurement (Ayehu et al., 2017). However, ground rainfall stations in many parts of the world
and most parts of Ethiopia especially, in relatively remote and mountainous areas are very spars
and unevenly distributed.
Correspondingly, rainfall data with very fine resolution and high quality are important for
hydrological computations (Zhou et al., 2017), in turn hydrological model output is highly affected
by precipitation input uncertainty ( Kamali et al., 2017; Haile et al., 2010). In spite of its
significance role in our lives and global climate, the measurement of precipitation is very difficult
due to its high spatial and temporal variability’s (Rodriguez et al., 2015b).
In practice, the density and configuration of rain gauge network is determined based on the
availability of funds, accessibility of site and the purpose of the network. As a result, most of the
time in-situ rain gauges instruments in developing countries are universally installed in towns that
located nearby main roads to provide accessibility.
In contrast, satellite remote sensing can make available the spatial precipitation data over large
areas in a temporally continuous way for various hydrological studies since 1970s’ (Liu et al.,
2014). However, due to the lack of quantitative confidence in the accuracy, these satellite rainfall
products are not often used for different operations (Burkea et al., 2017).
Over the past decades, countless efforts have been made to generate gridded precipitation products,
thus leading to the increasing accessibility of rainfall dataset at different spatial and temporal
resolution over the global or quasi-global scale (Schneider et al., 2015) and can be broadly
categorized into four classes (Duan et al., 2016).
1). Gauge-only products that build only on observations from in-situ rain gauge stations using
diverse interpolation methods which are often accessible at coarser spatial resolution (> 0.5º) 2).
Rainfall products from numerical weather forecasts or atmospheric models that uses a blend of
satellite and in-situ observations of several atmospheric properties as inputs (Ebert et al., 2007).
3). Satellite remote sensing-only products. 4). Satellite-gauge products that combine two
individual gauge-only and satellite-only data’s together through different bias correction or
blending technique and which products are often existing at the spatial resolution of 0.25º or finer.
Satellite is an object which has been purposely placed into an orbit to do specific work. Satellite
system that provide weather (precipitation) related observations are separated in to geostationary
GEO satellites and low earth orbiting LEO satellites (Aoki, 2017). Imagery from geostationary
satellite systems, when compared to polar-orbiting satellite information, normally, results in a
reduction of the sampling errors at all temporal scales.
Consequently, at the present time, estimates of global, tropical and subtropical precipitation at
short timescales (from few hours to perhaps a few weeks) are most accurately obtained from
geostationary satellites (Tang et al., 2016). Algorithms to estimate precipitation data from satellite
observations are based either on thermal infrared (TIR) bands or on passive microwave (PMW)
sensors (Duan et al., 2016), which have their own specific strengths and limitations.
Thermal infrared sensors on geostationary satellite produce high spatial coverage with continuous
temporal coverage, covering nearly the complete globe every one hour or less. Thermal infrared
information based on the cloud-top brightness is frequently valuable for the distinction between
raining and non-raining time (Yang et al., 2017). Several studies displayed that TIR-based
precipitation evaluations may have a large uncertainty when some types of cloud or warm clouds
are present (Huang et al., 2017).
On the other hand, the passive microwave sensors on polar-orbiting satellite recognize the rainfalls
particles by scattering due to the large ice particles present in the clouds. PMW is better at
approximating the precipitation quantity owing to the increased direct physical connection
between the sensor signal and precipitation, but runs at a much lower temporal frequency and
coarser spatial resolution (Bergsted et al., 2017).
To deal with these technical restrictions, the more recently satellite-based rainfall products, such
as CHIRPS_2 and TMPA_3B42v7 combine multiple data sources as TIR-PMW based rainfall data
sets.
2.2.1. CHIRPS_2
The Climate Hazards Group Infra-Red Precipitation with Station data (CHIRPS version 2) is a
global daily, pentadal (5), monthly precipitation product and a third generation rainfall procedures
clearly designed for monitoring agricultural drought and global environmental change over land
(Funk et al., 2015), which was developed by the U.S Geological Survey and the Climate Hazards
Group at the University of California, Santa Barbara in 2015.
CHIRPS product which belongs to the “satellite-gauge” rainfall source category combines
remotely sensed precipitation of geo-synchronous and polar orbiting satellites, from five different
most recommended satellite precipitation products, with more than 2000 station records to
calibrate global Cold Cloud Duration (CCD) precipitation estimates (Funk et al., 2015). The
product structures at a spatial resolution of 0.05º and 0.25ºfrom 50º S to 50º N across all longitudes
with more than 30-year final monthly rainfall records (1981-present).
The station statistics from the near-real time World Meteorological Organization’s Global
Telecommunication System (GTS) are constantly used to update roughly every two days, and
validate the remote sensors using extra evidence (physiographic and remotely sensed earth energy
emissions that agree to the location and intensity of rainfall) through moving window geo
statistical regression (Funk et al., 2015).The main used datasets for the creation of CHIRPS_2
product includes: -
1). The monthly Climate Hazard Precipitation climatology (CHPclim) that is molded with in-situ
rain gauge stations.
CHPclim is composed from the Agro-met Group of the Food and Agriculture Organization of the
United Nations (FAO), the Global Historical Climate Network (GHCN) and the Cold Cloud
Duration (CCD) information based on thermal infrared data archived from Climate Prediction
Center (CPC) and National Oceanographic and Atmospheric Administration (NOAA) National
Climate Data Center (NCDC).
2). The Version_7 TRMM 3B42 data. 3). The Version_2 atmospheric model precipitation data
field from the NOAA Climate Forecast System (CFS) and 4). The rain gauge stations data from
several sources. Primarily, the CCD records are calibrated with TRMM-3B42 to generate the five
daily CCD-based precipitation estimates which are further reformed to the fractions of the long-
term mean precipitation estimates.
The fractions are then multiplied with CHPclim data to eliminate the systematic bias and the
derived product is called CHIRP product. Finally, the CHIRP rainfall product is merged with rain
gauge stations information using a modified Inverse Distance Weighting (IDW) algorithm to
produce the CHIRPS_2.
All the process indicated above are accomplished at the five-daily (pentadal) time scales. The daily
CCD measurements and daily CFS measurements are finally used to disaggregate the 5-daily
products to daily rainfall estimates using a simple redistribution method which data is available at:
http://chg.geog.ucsb.edu/data/chirps/ with resolution of 0.05ºand 0.25º.
2.2.2. TMPA_3B42v7
The Tropical Rainfall Measuring Mission (TRMM) Multisatellite Precipitation Analysis (TMPA)
is a joint US and Japan satellite mission, which was launched in 1997 to monitor tropical and
subtropical precipitation data and to estimate its related latent heating covering the latitude band
50º N – 50º S.
It also intended to provide a best estimate of quasi-global precipitation data from a wide variety of
modern satellite-born rainfall related sensors. Its rainfall estimates are provided at relatively high
spatial resolution (0.25º×0.25º) and 3-hourly time stapes, in both real and post-real time to meet a
wide range of research needs (Huffman et al., 2010).
This satellite rainfall estimates combines IR data’s from geosynchronous earth orbit (GEO)
satellites with four passive microwave (PMW) sensors, namely TRMM Microwave Imager (TMI),
Special Sensor Microwave/Imager (SSM/I), Advanced Microwave Sounding Unit-B (AMSUB),
and the Advanced Microwave Scanning Radiometer – EOS (AMSR-E).
The TMPA_3B42v7 satellite precipitation estimates include the version (6 and 7) real-time
products 3B42RT (3B42RTv6 and 3B42RTv7), the 3-hourly research data 3B42 (3B42v6 and
3B42v7), and the daily accumulated precipitation products 3B42_daily (v6 and v7).
The 3B42 algorithm is executed in four steps:1). PMW rainfall estimates are calibrated and
combined; 2). TIR precipitation estimates are generated using the calibrated PMW data; 3). Both
TIR and PMW data are then combined and 4). Rescaled on a monthly basis using two sources of
rain gauge records: such as,
i). The Global Precipitation Climatology Center (GPCC) monthly rain gauge analysis and ii). The
Climate Assessment and Monitoring System (CAMS) monthly rain gauge analysis developed by
CPC (Bigiarini et al., 2017).
The latest version of TMPA_3B42v7 rainfall product was released in June 2012, and recent studies
indicated that 3B42v7 estimates improve upon 3B42v6 (Bigiarini et al., 2017; Chen et al., 2013).
Nowadays, after more than 17 years of data collection, the instruments turned off on 8, April, 2015,
after the spacecraft exhausted its fuel reserves and replaced by Global Precipitation Measurement
GPM mission. For this thesis work TMPA_3B42v7 daily data’s was used and can be easily
accessed from (https://pmm.nasa.gov/).
Several studies have been carried out to evaluate a single or multiple satellite precipitation products
on stream flow simulation capability with and without using hydrological models in scale from the
global or quasi-global (Massari et al., 2017), continental scale (Maidment et al.,2017), Country-
wide (Maidment et al., 2017), regional scale (Haile et al., 2016), basin scales (Duan et al., 2016)
and watershed scales (Liu et al., 2014; Stage et al., 2014).
However, an insufficient studies have been conducted in Ethiopia, for instances, (Bhatti et al.,
2016) evaluated CMORPH satellite rainfall products for Upper Gilgel Abay catchment. Habib et
al., (2014) assessed effect of bias correction on CMORPH satellite product. Haile et al., (2012)
evaluated CMORPH satellite rainfall products for Gilgel Abay watershed.
Gebremichaiel and Romilly, (2011) evaluated three satellite rainfall products (CMORPH,
PERSIAN and TRMM) for six Ethiopian river basins. The condition is primarily remarking that
because of latest availability of the (CHIRPS_2 and TMPA_3B42v7) high resolution satellite
rainfall products, very few investigations have been done to evaluate and compare the products.
Table 2.1 below show earlier studies about satellite precipitation products, the method they used
and their main gaps.
Table 2.1: Earlier Studies around the World and Ethiopia
In conclusion, based on the studies carried out in different regions of the world as well as in
Ethiopia, the evaluated satellite precipitation estimates showed that satellite precipitation products
have potential use for hydrologic modeling, nevertheless the performance of the
hydrological model for stream flow simulation were depends on the type of satellite rainfall
product, basin area, hydro-climatic regions, length of simulation periods and topography of the
study area.
Hydrological models are picture of the real world system and can be classified as physical,
conceptual and empirical based on the degree of complexity and physical completeness in the
establishment of the structure. Models are further classified as lumped, distributed and semi-
distributed depending on the degree of reorganization when describing the terrain in the basin
(Cunderlik, 2003).
Lumped models:- In this models the hydrologic parameters do not vary spatially within the
watershed and its response is evaluated only at the outlet, without explicitly account for the
response of individual sub-basins (Cunderlic, 2003).
Distributed models are generally require large amount of data in each grid cell. If governing
physical process are modeled in detail and suitably applied, distributed models can provide the
highest degree of accuracy. In semi-distributed model, the parameters of the models are allowed
to vary partially in space by dividing the basin in to a number of smaller sub-basins.
In deterministic models only one outcome is permit from a simulation with one set of inputs and
parameter values, whereas stochastic models allow for some randomness or uncertainty in the
possible results due to uncertainty in input variables, boundary conditions or model parameters
(Beven, 2000).
Conceptual models are constructed on the limited picture of the physical process acting to produce
the hydrological outputs, while physically based models are built more solidly on understanding
of the related process (Ward and Robinson, 2000).
According to Beven, (2000) many rainfall-run of models are carried out purely for research
purposes as a means of improving knowledge about hydrological systems. In addition to this,
models are developed as tools for simulation and prediction aiming ultimately to allow decision
makers to improve decision making about hydrological problems.
It is designed to simulate the broad hydrologic process of dendritic watershed system by including
many traditional hydrologic analysis procedures such as event infiltration, unit hydrographs, and
hydrologic routing.
It also includes processes necessary for continuous and gridded runoff simulation including
evapotranspiration, snowmelt and soil moisture accounting. The program features have an
integrated work environment with data entry utilities, computation engines and outcome reporting
tools.
Simulation results are stored in the data storage system HEC-DSS and can be used in union with
others software for studies of water accessibility, urban drainage, flow forecasting, future
urbanization influence, reservoir spillway design, flood damage reduction, flood plain regulation
and system operations (USACE, 2016).
The Geospatial Hydrologic Modeling System extension (HEC-GeoHMS) tools has been
developed as a geospatial hydrology toolkit for engineers and hydrologists with inadequate GIS
understanding (Alazzy et al., 2017). Examining digital terrain data, HEC-GeoHMS transforms the
drainage paths and sub-basin boundaries into a hydrologic data structure that describes the
drainage network.
The model allows user to visualize spatial information, to document watershed features, to perform
spatial investigation, and to delineate sub-basins. The currently free program used for this thesis
was version 4 (HEC-HMS 4.2.1) and focused principally on new computation features which has
an extensive array of capabilities for conducting hydrologic simulations.
HEC-HMS hydrological model is selected on the basis of applicability and restrictions of each
method, availability of data, well recognized, stable, widely acceptable, a lot of scholar
recommendations, simple, powerful tool and its versatility. The model can care to save time and
money in attaining the runoff data rather than measurement of runoff in the watershed.
The HEC-HMS hydrologic model was found precise in spatially and temporally forecasting
watershed response in event based and continuous simulation along with simulating various
scenarios in flood forecasting and early warnings. It is freely accessible for download from the
HEC-HMS website (www.hec.usace.army.mil) and is supported by the US Army Corps of
Engineers. Previous studies used method and range of parameters are presented in below table 2.2
and table 2.3 respectively.
Table 2.2: Method available in the sub basin and reach hydrologic elements
Table 2.3: Calibration parameters value of HEC-HMS model used in the past studies
Every simulation system has its own restriction due to the selections prepared in the design and
development of the software (USACE, 2016). Some of the limitation of HEC-HMS models are; -
1. Simplified model formulation - simplifying the model formulation allows to the program to
complete simulation very quickly whereas generating accurate and precise results.
2. Simplified flow representation – simplifying the flow representation aids in keeping the compute
process efficient and reduces repetition of capability in the HEC software suite.
There are two main different calibration methods for simulating stream flow from satellite
precipitation inputs:
1). Calibrating the model using rain gauge rainfall measurements and validating the model with
satellite rainfall products (Elgamal et al., 2017).
2). Calibrating and validating the model with the satellite rainfall products (Bajracharya et al.,
2017; Stisen and Sandholt, 2010).
However, compared to calibrating hydrologic model using in-situ rain gauge precipitation,
calibrating the model with consistent satellite precipitation data increases the performance of
satellite rainfall stream flow simulation but decrease the performances of satellite evapo-
transpiration simulations (Bitew and Gebremichael, 2011a).