skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Analysis of the Spatial Patterns of Rainfall across the Agro-Climatic Zones of Jema Watershed in the Northwestern Highlands of Ethiopia
The association between elevation (agro-climatic zones, ACZs) and the mean annual total rainfall (MATRF) is not straightforward in different parts of the world. This study sought to estimate the amount of MATRF across four elevation zones of Jema watershed, which is situated in the northwestern highlands of Ethiopia, by employing an appropriate interpolation method. The elevation of the watershed ranges from 1895 to 3518 m a.s.l. For the sake of this study, 34 sample MATRF data were extracted from satellite and nearby gauge stations that were recorded from 1983 to 2010. These data sources were reconstructed by International Research Institute for Climate and Society at Columbia University, USA, at a scale of 10 km by 10 km. An elevation data set generated from a digital elevation model with 30-m resolution (DEM 30 m) was considered as a covariable to estimate the MATRF. To identify the optimal interpolation model, mean errors were computed using cross-validation statistics. The root-mean-square error (RMSE) analysis showed that ordinary cokriging (OCK) was the most accurate model with a predictive power of 87.3%. The root-mean-square standardized (RMSSE) analysis showed that the best precision value (0.72) occurred in OCK. Stable and Gaussian trend lines together with local polynomial types of trend removal, and an elliptical neighborhood search function could perform best to maximize the accuracy and the precision of estimating MATRF. Elevation, as a covariable, enhanced the degree of accuracy and precision of estimation. The value of the trend line function (least square) between the MATRF and elevation was very weak (R2 = 0.07), whereas the value of trend line function (least square) between the MATRF and the longitude coordinates (east–west direction) was medium (R2 = 0.34). The estimated MATRF for the entire watershed under study ranged from 1228 to 1640 mm. To conclude, elevation could contribute to the estimation of the MATRF. The value of the MATRF showed a declining pattern from the lower to higher elevation areas of the watershed.  more » « less
Award ID(s):
1639214
PAR ID:
10109676
Author(s) / Creator(s):
; ; ; ;
Date Published:
Journal Name:
Geosciences
Volume:
9
Issue:
1
ISSN:
2076-3263
Page Range / eLocation ID:
22
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    The measurement of sea ice elevation above sea level or the “freeboard” depends upon an accurate retrieval of the local sea level. The local sea level has been previously retrieved from altimetry data alone by the lowest elevation method, where the percentage of the lowest elevations over a particular segment length scale was used. Here, we provide an evaluation of the scale dependence on these local sea level retrievals using data from NASA Operation IceBridge (OIB) which took place in the Ross Sea in 2013. This is a unique dataset of laser altimeter measurements over five tracks from the Airborne Topographic Mapper (ATM), with coincidently high-spatial resolution images from the Digital Mapping System (DMS), that allows for an independent sea level validation. The local sea level is first calculated by using the mean elevation of ATM L1B data over leads identified by using the corresponding DMS imagery. The resulting local sea level reference is then used as ground truth to validate the local sea levels retrieved from ATM L2 by using nine different percentages of the lowest elevation (0.1%, 0.5%, 1%, 1.5%, 2%, 2.5%, 3%, 3.5%, and 4%) at seven different segment length scales (1, 5, 10, 15, 20, 25, and 50 km) for each of the five ATM tracks. The closeness to the 1:1 line, R2, and root mean square error (RMSE) is used to quantify the accuracy of the retrievals. It is found that all linear least square fits are statistically significant (p < 0.05) using an F test at every scale for all tested data. In general, the sea level retrievals are farther away from the 1:1 line when the segment length scale increases from 1 or 5 to 50 km. We find that the retrieval accuracy is affected more by the segment length scale than the percentage scale. Based on our results, most retrievals underestimate the local sea level; the longer the segment length (from 1 to 50 km) used, especially at small percentage scales, the larger the error tends to be. The best local sea level based on a higher R2 and smaller RMSE for all the tracks combined is retrieved by using 0.1–2% of the lowest elevations at the 1–5 km segment lengths. 
    more » « less
  2. In this study, we present a nationwide machine learning model for hourly PM2.5 estimation for the continental United States (US) using high temporal resolution Geostationary Operational Environmental Satellites (GOES-16) Aerosol Optical Depth (AOD) data, meteorological variables from the European Center for Medium Range Weather Forecasting (ECMWF) and ancillary data collected between May 2017 and December 2020. A model sensitivity analysis was conducted on predictor variables to determine the optimal model. It turns out that GOES16 AOD, variables from ECMWF, and ancillary data are effective variables in PM2.5 estimation and historical reconstruction, which achieves an average mean absolute error (MAE) of 3.0 μg/m3, and a root mean square error (RMSE) of 5.8 μg/m3. This study also found that the model performance as well as the site measured PM2.5 concentrations demonstrate strong spatial and temporal patterns. Specifically, in the temporal scale, the model performed best between 8:00 p.m. and 11:00 p.m. (UTC TIME) and had the highest coefficient of determination (R2) in Autumn and the lowest MAE and RMSE in Spring. In the spatial scale, the analysis results based on ancillary data show that the R2 scores correlate positively with the mean measured PM2.5 concentration at monitoring sites. Mean measured PM2.5 concentrations are positively correlated with population density and negatively correlated with elevation. Water, forests, and wetlands are associated with low PM2.5 concentrations, whereas developed, cultivated crops, shrubs, and grass are associated with high PM2.5 concentrations. In addition, the reconstructed PM2.5 surfaces serve as an important data source for pollution event tracking and PM2.5 analysis. For this purpose, from May 2017 to December 2020, hourly PM2.5 estimates were made for 10 km by 10 km and the PM2.5 estimates from August through November 2020 during the period of California Santa Clara Unite (SCU) Lightning Complex fires are presented. Based on the quantitative and visualization results, this study reveals that a number of large wildfires in California had a profound impact on the value and spatial-temporal distributions of PM2.5 concentrations. 
    more » « less
  3. This study was conducted within the context of providing an improved understanding of New Mexico’s food, energy, water systems (FEWS) and their behavior under variable climate and socioeconomic conditions. The goal of this paper was to characterize the relationships between production and prices of some forage crops (hay, grain sorghum, and corn) that can be used as feed supplements for beef cattle production and the potential impacts from a changing climate (precipitation, temperature) and energy inputs (crude oil production and prices). The analysis was based on 60 years of data (1958–2017) using generalized autoregressive conditional heteroscedasticity models. Hay production showed a declining trend since 2000 and in 2017, it dropped by ~33% compared to that of 2000. Crude oil production (R2 = 0.83) and beef cattle population (R2 = 0.85) were negatively correlated with hay production. A moderate declining trend in mean annual hay prices was also observed. Mean annual range conditions (R2 = 0.60) was negatively correlated with mean annual hay prices, whereas mean annual crude oil prices (R2 = 0.48) showed a positive relationship. Grain sorghum production showed a consistent declining trend since 1971 and in 2017, it dropped by ~91% compared to that of 1971. Mean annual temperature (R2 = 0.58) was negatively correlated with grain sorghum production, while beef cattle population (R2 = 0.61) and range conditions (R2 = 0.51) showed positive linear relationships. Mean annual grain sorghum prices decreased since the peak of 1974 and in 2017, they dropped by ~77% compared to those of 1974. Crude oil prices (R2 = 0.72) and beef cattle population (R2 = 0.73) were positively correlated with mean annual grain sorghum prices. Corn production in 2017 dropped by ~61% compared to the peak that occurred in 1999. Crude oil production (R2 = 0.85) and beef cattle population (R2 = 0.86) were negatively correlated with corn production. Mean annual corn prices showed a declining trend since 1974 and in 2017, they dropped by ~75% compared to those of 1974. Mean annual corn prices were positively correlated with mean annual precipitation (R2 = 0.83) and negatively correlated with crude oil production (R2 = 0.84). These finding can particularly help in developing a more holistic model that integrates FEWS components to explain their response to internal (i.e., management practices) and external (i.e., environmental) stressors. Such holistic modeling can further inform the development and adoption of more sustainable production and resource use practices. 
    more » « less
  4. null (Ed.)
    High-quality retrieval of land surface phenology (LSP) is increasingly important for understanding the effects of climate change on ecosystem function and biosphere–atmosphere interactions. We analyzed four state-of-the-art phenology methods: threshold, logistic-function, moving-average and first derivative based approaches, and retrieved LSP in the North Hemisphere for the period 1999–2017 from Copernicus Global Land Service (CGLS) SPOT-VEGETATION and PROBA-V leaf area index (LAI) 1 km V2.0 time series. We validated the LSP estimates with near-surface PhenoCam and eddy covariance FLUXNET data over 80 sites of deciduous forests. Results showed a strong correlation (R2 > 0.7) between the satellite LSP and ground-based observations from both PhenoCam and FLUXNET for the timing of the start (SoS) and R2 > 0.5 for the end of season (EoS). The threshold-based method performed the best with a root mean square error of ~9 d with PhenoCam and ~7 d with FLUXNET for the timing of SoS (30th percentile of the annual amplitude), and ~12 d and ~10 d, respectively, for the timing of EoS (40th percentile). 
    more » « less
  5. ArcticDEM provides the public with an unprecedented opportunity to access very high-spatial resolution digital elevation models (DEMs) covering the pan-Arctic surfaces. As it is generated from stereo-pairs of optical satellite imagery, ArcticDEM represents a mixture of a digital surface model (DSM) over a non-ground areas and digital terrain model (DTM) at bare grounds. Reconstructing DTM from ArcticDEM is thus needed in studies requiring bare ground elevation, such as modeling hydrological processes, tracking surface change dynamics, and estimating vegetation canopy height and associated forest attributes. Here we proposed an automated approach for estimating DTM from ArcticDEM in two steps: (1) identifying ground pixels from WorldView-2 imagery using a Gaussian mixture model (GMM) with local refinement by morphological operation, and (2) generating a continuous DTM surface using ArcticDEMs at ground locations and spatial interpolation methods (ordinary kriging (OK) and natural neighbor (NN)). We evaluated our method at three forested study sites characterized by different canopy cover and topographic conditions in Livengood, Alaska, where airborne lidar data is available for validation. Our results demonstrate that (1) the proposed ground identification method can effectively identify ground pixels with much lower root mean square errors (RMSEs) (<0.35 m) to the reference data than the comparative state-of-the-art approaches; (2) NN performs more robustly in DTM interpolation than OK; (3) the DTMs generated from NN interpolation with GMM-based ground masks decrease the RMSEs of ArcticDEM to 0.648 m, 1.677 m, and 0.521 m for Site-1, Site-2, and Site-3, respectively. This study provides a viable means of deriving high-resolution DTM from ArcticDEM that will be of great value to studies focusing on the Arctic ecosystems, forest change dynamics, and earth surface processes. 
    more » « less