An Intercomparison Study of Algorithms for Downscaling SMAP Radiometer Soil Moisture Retrievals
Explore the source record for details and available documents.
SEARCH · Search NASA
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
The mountain regions of the northeastern United States are a critical socioeconomic resource for Vermont, New York State, New Hampshire, Maine, and southern Quebec. While global climate models (GCMs) are important tools for climate change risk assessment at regional scales, even the increased spatial resolution of statistically downscaled GCMs (commonly approximately 1/ 8 deg) is not sufficient for hydrologic, ecologic, and land-use modeling of small watersheds within the mountainous Northeast. To address this limitation, an ensemble of topographically downscaled, high-resolution (30"), daily 2-m maximum air temperature; 2-m minimum air temperature; and precipitation simulations are developed for the mountainous Northeast by applying an additional level of downscaling to intermediately downscaled (1/ 8 deg) data using high-resolution topography and station observations. First, observed relationships between 2-m air temperature and elevation and between precipitation and elevation are derived. Then, these relationships are combined with spatial interpolation to enhance the resolution of intermediately downscaled GCM simulations. The resulting topographically downscaled dataset is analyzed for its ability to reproduce station observations. Topographic downscaling adds value to intermediately downscaled maximum and minimum 2-m air temperature at high-elevation stations, as well as moderately improves domain-averaged maximum and minimum 2-m air temperature. Topographic downscaling also improves mean precipitation but not daily probability distributions of precipitation. Overall, the utility of topographic downscaling is dependent on the initial bias of the intermediately downscaled product and the magnitude of the elevation adjustment. As the initial bias or elevation adjustment increases, more value is added to the topographically downscaled product.
The Soil Moisture Active Passive (SMAP) mission is dedicated toward global soil moisture mapping. Typically, an L-band microwave radiometer has spatial resolution on the order of 36-40 km, which is too coarse for many specific hydro-meteorological and agricultural applications. With the failure of the SMAP active radar within three months of becoming operational, an intermediate (9-km) and finer (3-km) scale soil moisture product solely from the SMAP mission is no longer possible. Therefore, the focus of this study is a disaggregation of the 36-km resolution SMAP passive-only surface soil moisture (SSM) using the Soil Evaporative Efficiency (SEE) approach to spatial scales of 3-km and 9-km. The SEE was computed using thermal-infrared (TIR) estimation of surface evaporation over Continental U.S. (CONUS). The disaggregation results were compared with the 3 months of SMAP-Active (SMAP-A) and Active/Passive (AP) products, while comparisons with SMAP-Enhanced (SMAP-E), SMAP-Passive (SMAP-P), as well as with more than 180 Soil Climate Analysis Network (SCAN) stations across CONUS were performed for a 19 month period. At the 9-km spatial scale, the TIR-Downscaled data correlated strongly with the SMAP-E SSM both spatially (r = 0.90) and temporally (r = 0.87). In comparison with SCAN observations, overall correlations of 0.49 and 0.47; bias of −0.022 and −0.019 and unbiased RMSD of 0.105 and 0.100 were found for SMAP-E and TIR-Downscaled SSM across the Continental U.S., respectively. At 3-km scale, TIR-Downscaled and SMAP-A had a mean temporal correlation of only 0.27. In terms of gain statistics, the highest percentage of SCAN sites with positive gains (>55%) was observed with the TIR-Downscaled SSM at 9-km. Overall, the TIR-based downscaled SSM showed strong correspondence with SMAP-E; compared to SCAN, and overall both SMAP-E and TIR-Downscaled performed similarly, however, gain statistics show that TIR-Downscaled SSM slightly outperformed SMAP-E.
The Water Cycle Observation Mission (WCOM) will build upon previous L and C band passive microwave soil moisture satellite missions. WCOM will consist of a passive microwave synthetic aperture radiometer operating at L, S, and C bands. The WCOM requirements for passive soil moisture are to estimate soil moisture in the top 5 cm of soil layer with an error less than 0.04 cu m/cu m, at 15 km resolution and with a 3-day revisit. A new set of algorithms for these multi-frequency platforms will need to be developed for estimating the data products at the desired resolution. To accomplish this, a brightness temperature (TB) downscaling methodology is developed that uses passive S-band TB (30 km) to downscale L-band TB (50 km) and to estimate soil moisture at a 30 km resolution, based on the linear relationships between the passive signals of L-band and S-band. To test this downscaling method, analysis was performed using PALS data from the Soil Moisture Experiments in 2002 (SMEX02). For this study, 4 km L-band observations were downscaled to 800 m. The root mean square errors (RMSE) between the downscaled TBL at 800m with the observed TBL at 800m are 2.63 K and 1.60 K for H and V polarizations respectively. The results also showed that it was possible to use these disaggregated TB to estimate soil moisture to meet the mission requirement of 0.04 cu m/cu m. These results showed that we can obtain higher resolution soil moisture from L band passive TB with a high accuracy (<0.04 cu m/cu m) by using S-band information.
Regional quantification of energy and water balance fluxes depends inevitably on the estimation of surface and rootzone soil moisture. The simulation of soil moisture depends on the soil retention characteristics, which are difficult to estimate at a regional scale. Thus, the present study proposes a new method to estimate high-resolution Soil Hydraulic Parameters (SHPs) which in turn help to provide high-resolution (spatial and temporal) rootzone soil moisture (RZSM) products. The study is divided into three phases—(I) involves the estimation of finer surface soil moisture (1 km) from the coarse resolution satellite soil moisture. The algorithm utilizes MODIS 1 km Land Surface Temperature (LST) and 1 km Normalized difference vegetation Index (NDVI) for downscaling 25 km C-band derived soil moisture from AMSR-2 to 1 km surface soil moisture product. At one of the test sites, soil moisture is continuously monitored at 5, 20, and 50 cm depth, while at 44 test sites data were collected randomly for validation. The temporal and spatial correlation for the downscaled product was 70% and 83%, respectively. (II) In the second phase, downscaled soil moisture product is utilized to inversely estimate the SHPs for the van Genuchten model (1980) at 1 km resolution. The numerical experiments were conducted to understand the impact of homogeneous SHPs as compared to the three-layered parameterization of the soil profile. It was seen that the SHPs estimated using the downscaled soil moisture (I-d experiment) performed with similar efficiency as compared to SHPs estimated from the in-situ soil moisture data (I-b experiment) in simulating the soil moisture. The normalized root mean square error (nRMSE) for the two treatments was 0.37 and 0.34, respectively. It was also noted that nRMSE for the treatment with the utilization of default SHPs (I-a) and AMSR-2 soil moisture (I-c) were found to be 0.50 and 0.43, respectively. (III) Finally, the derived SHPs were used to simulate both surface soil moisture and RZSM. The final product, RZSM which is the daily 1 km product also showed a nearly 80% correlation at the test site. The estimated SHPs are seen to improve the mean NSE from 0.10 (I-a experiment) to 0.50 (I-d experiment) for the surface soil moisture simulation. The mean nRMSE for the same was found to improve from 0.50 to 0.31.
With its unique and trusted earth observations, NASA is a critical source in informing decisions that will help achieve the U.S. goal of Net-Zero Greenhouse Gas (GHG) Emissions by 2050. NASA’s Prediction of Worldwide Energy Resource (POWER) project facilitates the use of NASA Earth Science data holdings within the energy, agricultural, and building heating/cooling design industries. POWER packages solar and meteorological data from several NASA projects in a user friendly GIS-enabled web services system (https://power.larc.nasa.gov). As part of the development of new data products to support the energy and building heating/cooling design communities, we estimate the changes in energy required to heat and cool buildings in the future climate at 14 different NASA site locations spread throughout the continental United States, as projected by CMIP6 climate models under different emissions scenarios. Bias-corrected downscaled time series of meteorological variables are taken from NASA Earth Exchange (NEX) Global Daily Downscaled Projections (GDDP-CMIP6) downscaled climate model data. The spread between the different model projections is accounted for by analyzing both the ensemble average of 22 CMIP6 models and 6 representative models with different climate sensitivities and different interannual variability. Changes in energy use are estimated in two ways. First, changes in the total annual heating and cooling degree days (HDD and CDD, respectively) are calculated relative to the current climate. This is done at all 14 sites. Second, the downscaled time series are used as inputs into RETScreen(R), a clean energy management decision tool, to give an estimate of heating/cooling energy use for a typical office building. We use this estimation method with model data at Langley Research Center. In the next 50 years, the annual total of HDD (CDD) is projected to decrease by 8-38% (increase by 5-28%) at all sites, with the increase in CDD typically a larger magnitude the decrease in HDD. For a typical small office building at Langley Research Center, the amount of energy needed to cool increases by 33-54% and the amount of energy to heat decreases by 29-40%. POWER is working to develop long term climate data services based on these results to include in future data products to provide to users.
A new computational analysis tool, downscaling test, is introduced and applied for studying the convergence rates of truncation and discretization errors of nite-volume discretization schemes on general irregular (e.g., unstructured) grids. The study shows that the design-order convergence of discretization errors can be achieved even when truncation errors exhibit a lower-order convergence or, in some cases, do not converge at all. The downscaling test is a general, efficient, accurate, and practical tool, enabling straightforward extension of verification and validation to general unstructured grid formulations. It also allows separate analysis of the interior, boundaries, and singularities that could be useful even in structured-grid settings. There are several new findings arising from the use of the downscaling test analysis. It is shown that the discretization accuracy of a common node-centered nite-volume scheme, known to be second-order accurate for inviscid equations on triangular grids, degenerates to first order for mixed grids. Alternative node-centered schemes are presented and demonstrated to provide second and third order accuracies on general mixed grids. The local accuracy deterioration at intersections of tangency and in flow/outflow boundaries is demonstrated using the DS tests tailored to examining the local behavior of the boundary conditions. The discretization-error order reduction within inviscid stagnation regions is demonstrated. The accuracy deterioration is local, affecting mainly the velocity components, but applies to any order scheme.
Global soil moisture products that are completely independent of any type of ancillary data and solely rely on satellite observations are presented. Additionally, we further develop an existing downscaling technique that enhances the spatial resolution of such products to approximately 11 km. These products are based on internal modules of the Land Parameter Retrieval Model (LPRM), an algorithm that uses the radiative transfer equation to link soil moisture, vegetation optical depth and land surface temperature to observed brightness temperatures. The soil moisture product that is independent of any type of ancillary data uses the internally calculated dielectric constant as a soil moisture proxy. This data product is not influenced by errors associated with coarse-scale global soil property maps or by any other type of forcing (e.g. re-analysis) data and is therefore solely based on satellite microwave observations. The second step builds upon recent developments to increase the spatial resolution of the LPRM retrievals using a smoothing filter downscaling method. With this method we can attain a spatial resolution that can be more useful at the scale of local and regional hydrological studies as well. The steps presented in this paper were applied to observations from the Advanced Microwave Scanning Radiometer-Earth Observing System (AMSR-E). The newly derived data sets were validated using ground-based observations from the International Soil Moisture Network (ISMN). The internally calculated dielectric constant product results in significantly more days with valid retrievals than the original soil moisture data products, in particular over arid regions. The dielectric constant product resulted in similar correlations with in situ data as the original soil moisture data product. Together, these findings demonstrate the usefulness of this new dielectric constant product for the hydrological modeling community and climate studies. A case study on the Australian Fitzroy catchment demonstrated that the downscaled data product has a more detailed spatial description of soil moisture, especially during wet and dry conditions with more pronounced dry and wet regions within the catchment. The increased resolution data products could therefore improve runoff predictions and this study demonstrated the potential added value of a transitioning from a spatial resolution of 56 km toward a higher resolution of 11 km. The hydrological implications of these newly developed data records are not only linked to AMSR-E satellite data, but also to the next generation Soil Moisture Active and Passive (SMAP) mission where a 9 km spatial resolution is the target resolution for satellite soil moisture products. The new data products will not replace the current LPRM products, but will be added to the existing array of data products and will become publicly available through our data portals.
NASA's Global Modeling and Assimilation Office at Goddard Space Flight Center is undertaking a series of very computationally intensive Nature Runs and a downscaled reanalysis. The nature runs use the GEOS-5 as an Atmospheric General Circulation Model (AGCM) while the reanalysis uses the GEOS-5 in Data Assimilation mode. This paper will present computational challenges from three runs, two of which are AGCM and one is downscaled reanalysis using the full DAS. The nature runs will be completed at two surface grid resolutions, 7 and 3 kilometers and 72 vertical levels. The 7 km run spanned 2 years (2005-2006) and produced 4 PB of data while the 3 km run will span one year and generate 4 BP of data. The downscaled reanalysis (MERRA-II Modern-Era Reanalysis for Research and Applications) will cover 15 years and generate 1 PB of data. Our efforts to address the big data challenges of climate science, we are moving toward a notion of Climate Analytics-as-a-Service (CAaaS), a specialization of the concept of business process-as-a-service that is an evolving extension of IaaS, PaaS, and SaaS enabled by cloud computing. In this presentation, we will describe two projects that demonstrate this shift. MERRA Analytic Services (MERRA/AS) is an example of cloud-enabled CAaaS. MERRA/AS enables MapReduce analytics over MERRA reanalysis data collection by bringing together the high-performance computing, scalable data management, and a domain-specific climate data services API. NASA's High-Performance Science Cloud (HPSC) is an example of the type of compute-storage fabric required to support CAaaS. The HPSC comprises a high speed Infinib and network, high performance file systems and object storage, and a virtual system environments specific for data intensive, science applications. These technologies are providing a new tier in the data and analytic services stack that helps connect earthbound, enterprise-level data and computational resources to new customers and new mobility-driven applications and modes of work. In our experience, CAaaS lowers the barriers and risk to organizational change, fosters innovation and experimentation, and provides the agility required to meet our customers' increasing and changing needs
South and Southeast Asia is subject to significant hydrometeorological extremes, including drought. Under rising temperatures, growing populations, and an apparent weakening of the South Asian monsoon in recent decades, concerns regarding drought and its potential impacts on water and food security are on the rise. Reliable sub-seasonal to seasonal (S2S) hydrological forecasts could, in principle, help governments and international organizations to better assess risk and act in the face of an oncoming drought. Here, we leverage recent improvements in S2S meteorological forecasts and the growing power of Earth Observations to provide more accurate monitoring of hydrological states for forecast initialization. Information from both sources is merged in a South and Southeast Asia sub-seasonal to seasonal hydrological forecasting system (SAHFS-S2S), developed collaboratively with the NASA SERVIR program and end-users across the region. This system applies the Noah-MultiParameterization (NoahMP) Land Surface Model (LSM) in the NASA Land Information System (LIS), driven by downscaled meteorological fields from the Global Data Assimilation System (GDAS) and Climate Hazards InfraRed Precipitation products (CHIRP and CHIRPS) to optimize initial conditions. The NASA Goddard Earth Observing System Model - sub-seasonal to seasonal (GEOS-S2S) forecasts, downscaled using the National Center for Atmospheric Research (NCAR) General Analog Regression Downscaling (GARD) tool and quantile mapping, are then applied to drive 5-km resolution hydrological forecasts to a 9-month forecast time horizon. Results show that the skillful predictions of root zone soil moisture can be made one to two months in advance for forecasts initialized in rainy seasons and up to 8 months when initialized in dry seasons. The memory of accurate initial conditions can positively contribute to forecast skills throughout the entire 9-month prediction period in areas with limited precipitation. This SAHFS-S2S has been operationalized at the International Centre for Integrated Mountain Development (ICIMOD) to support drought monitoring and warning needs in the region.
Streamflow is one of the key variables in the hydrological cycle. Simulation and forecasting of streamflow are challenging tasks for hydrologists, especially in sparsely gauged areas. Coarse spatial resolution remote sensing soil moisture products (equal to or larger than 9 km) are often assimilated into hydrological models to improve streamflow simulation in large catchments. This study uses the Ensemble Kalman Filter (EnKF) technique to assimilate SMAP soil moisture products at the coarse spatial resolution of 9 km (SMAP 9 km), and downscaled SMAP soil moisture product at the higher spatial resolution of 1 km (SMAP 1 km), into the Soil and Water Assessment Tool (SWAT) to investigate the usefulness of different spatial and temporal resolutions of remotely sensed soil moisture products in streamflow simulation and forecasting. The experiment was set up for eight catchments across the tropical climate of Vietnam, with varying catchment areas from 267 to 6430 km^2 during the period 2017–2019. We comprehensively evaluated the EnKF-based SWAT model in simulating streamflow at low, average, and high flow. Our results indicated that high-spatial resolution of downscaled SMAP 1 km is more beneficial in the data assimilation framework in aiding the accuracy of streamflow simulation, as compared to that of SMAP 9 km, especially for the small catchments. Our analysis on the impact of observation resolution also indicates that the improvement in the streamflow simulation with data assimilation is more significant at catchments where downscaled SMAP 1 km has fewer missing observations. This study is helpful for adding more understanding of performances of soil moisture data assimilation based hydrological modelling over the tropical climate region, and exhibits the potential use of remote sensing data assimilation in hydrology.
The soil moisture (SM) data retrieved from the Soil Moisture Active and Passive (SMAP) satellite are available at a 9 km grid spacing since April 2015. This product can provide valuable information for research and applications in hydrology and other related fields. However, the resolution may be too coarse for applications at catchment or field scale. In this study, an established downscaling methodology, which had a major modification regarding its application on the SMAP 33 km domain, was implemented to develop a 1 km soil moisture product based on the SMAP 9 km data. The algorithm proposed here is based on the thermal inertia principle and developed by modeling the relationship between surface temperature difference and SM for different Normalized Difference Vegetation Index (NDVI) classes. The model functions were established and tuned using data from the NASA’s Land Information System (LIS) North America Land Data Assimilation System (NLDAS) and remotely sensed VISible/InfRared (VIS/IR) reflectance data from Long Term Data Record (LTDR) AVHRR (Advanced Very High Resolution Radiometer) for the growing season months of April-September 1981–2018. These were then implemented using the MODIS (Moderate Resolution Imaging Spectroradiometer) data over the Continental United States (CONUS) domain. Validation activities were carried out using in situ measurements distributed through the International Soil Moisture Network (ISMN). The validation results computed using the 1 km SM data showed that the R2, unbiased RMSE (root mean square error) and bias were improved relative to the 9 km SMAP product by 0.045, 0.018m3/m3 and 0.001m3/m3, respectively. The 1 km SM also exhibited a strong time-series autocorrelation. Further accuracy assessment analyses indicated that precipitation might contribute to the uncertainties in both the 9 km SMAP and 1 km downscaled SMAP SM products.
NASA Earth Exchange (NEX) combines state-of-the-art supercomputing, Earth system modeling, and NASA remote sensing data feeds to deliver a work environment for exploring and analyzing petabyte-scale datasets covering large regions, continents, or the globe. As an accessible platform, NEX can accelerate fundamental research, develop new applications, and reduce overall project costs by providing the research community with data, software, and high-end computing power. Two research thrusts for the NEX community are developing and distributing the NEX-GDDP-CMIP6 downscaled climate dataset and the community-driven wildfire research. NEX-GDDP-CMIP6 dataset is comprised of global downscaled climate scenarios derived from the General Circulation Model (GCM) runs conducted under the Coupled Model Intercomparison Project Phase 6 (CMIP6) and across two of the four “Tier 1” greenhouse gas emissions scenarios known as Shared Socioeconomic Pathways (SSPs). The wildfire research thrust leverages geostationary and low earth orbit remote sensing platforms and advanced modeling (WRF). We present a high-level overview of the NEX community, discuss relevant datasets and show several visualization techniques of data used in the wildfire research activity.
The method of elevation classes, in which the ice surface model is run at multiple elevations within each grid cell, has proven to be a useful way for a low-resolution atmosphere inside a general circulation model (GCM) to produce high-resolution downscaled surface mass balance fields for use in one-way studies coupling atmospheres and ice flow models. Past uses of elevation classes have failed to conserve mass and energy because the transformation used to regrid to the atmosphere was inconsistent with the transformation used to downscale to the ice model. This would cause problems for two-way coupling. A strategy that resolves this conservation issue has been designed and is presented here. The approach identifies three grids between which data must be regridded and five transformations between those grids required by a typical coupled atmosphere-ice flow model. This paper develops a theoretical framework for the problem and shows how each of these transformations may be achieved in a consistent, conservative manner. These transformations are implemented in Glint2, a library used to couple atmosphere models with ice models. Source code and documentation are available for download. Confounding real-world issues are discussed, including the use of projections for ice modeling, how to handle dynamically changing ice geometry, and modifications required for finite element ice models.
This work details use of the North American Multi-Model Ensemble (NMME) experimental forecasts as drivers for Decision Support Systems (DSSs) in the NASA / USAID initiative, SERVIR (a Spanish acronym meaning "to serve"). SERVIR integrates satellite observations, ground-based data and forecast models to monitor and forecast environmental changes and to improve response to natural disasters. Through the use of DSSs whose "front ends" are physically based models, the SERVIR activity provides a natural testbed to determine the extent to which NMME monthly to seasonal projections enable scientists, educators, project managers and policy implementers in developing countries to better use probabilistic outlooks of seasonal hydrologic anomalies in assessing agricultural / food security impacts, water availability, and risk to societal infrastructure. The multi-model NMME framework provides a "best practices" approach to probabilistic forecasting. The NMME forecasts are generated at resolution more coarse than that required to support DSS models; downscaling in both space and time is necessary. The methodology adopted here applied model output statistics where we use NMME ensemble monthly projections of sea-surface temperature (SST) and precipitation from 30 years of hindcasts with observations of precipitation and temperature for target regions. Since raw model forecasts are well-known to have structural biases, a cross-validated multivariate regression methodology (CCA) is used to link the model projected states as predictors to the predictands of the target region. The target regions include a number of basins in East and South Africa as well as the Ganges / Baramaputra / Meghna basin complex. The MOS approach used address spatial downscaling. Temporal disaggregation of monthly seasonal forecasts is achieved through use of a tercile bootstrapping approach. We interpret the results of these studies, the levels of skill by several metrics, and key uncertainties.
An analysis of observed trends in African annual-average near-surface temperatures over the last five decades reveals drastic increases, particularly over parts of the subtropics and central tropical Africa. Over these regions, temperatures have been rising at more than twice the global rate of temperature increase. An ensemble of high-resolution downscalings, obtained using a single regional climate model forced with the sea-surface temperatures and sea-ice fields of an ensemble of global circulation model (GCM) simulations, is shown to realistically represent the relatively strong temperature increases observed in subtropical southern and northern Africa. The amplitudes of warming are generally underestimated, however. Further warming is projected to occur during the 21st century, with plausible increases of 4-6 C over the subtropics and 3-5 C over the tropics by the end of the century relative to present-day climate under the A2 (a low mitigation) scenario of the Special Report on Emission Scenarios. High impact climate events such as heat-wave days and high fire-danger days are consistently projected to increase drastically in their frequency of occurrence. General decreases in soil-moisture availability are projected, even for regions where increases in rainfall are plausible, due to enhanced levels of evaporation. The regional downscalings presented here, and recent GCM projections obtained for Africa, indicate that African annual-averaged temperatures may plausibly rise at about 1.5 times the global rate of temperature increase in the subtropics, and at a somewhat lower rate in the tropics. These projected increases although drastic, may be conservative given the model underestimations of observed temperature trends. The relatively strong rate of warming over Africa, in combination with the associated increases in extreme temperature events, may be key factors to consider when interpreting the suitability of global mitigation targets in terms of African climate change and climate change adaptation in Africa.