Search NASA⌕ Search

SEARCH · Search NASA

Results for “Ocean salinity”

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.

At least 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
NASA NTRS2018

Integration of SMAP and SMOS Observations

Soil Moisture Active Passive (SMAP) mission and the Soil Moisture and Ocean Salinity (SMOS) missions provide brightness temperature and soil moisture estimates every 2-3 days. SMAP brightness temperature observations were compared with SMOS observations at 40o incidence angle. The brightness temperatures from the two missions are not consistent. SMAP observations show a warmer TB bias (about 1.27 K: V pol and 0.62 K: H pol) as compared to SMOS. SMAP and SMOS missions use different retrieval algorithms and ancillary datasets which result in further inconsistencies between their soil moisture products. The reprocessed constant-angle SMOS brightness temperatures were used in the SMAP soil moisture retrieval algorithm to develop a consistent multi-satellite product. The integrated product has an increased global revisit frequency (1 day) and period of record that is unattainable by either one of the satellites alone. Results from the development and validation of the integrated soil moisture product will be presented.

passive microwave↗
NASA NTRS2018

Assimilation of Satellite Derived Soil Moisture Profiles into a Crop Modeling System for Robust Yield Estimates

Soil Moisture Measurement - Remote Sensing. - Microwave (MW) Remote Sensing: Physically based and quantitative in nature; Based on difference in dielectric constant; Coarse spatial resolution 25-40 km; Shallow SM estimation depth 0-5 cm (approx.); All weather capabilities (e.g. Advanced Microwave Scanning Radiometer - Earth Observing System (AMSR-E), Soil Moisture and Ocean Salinity (SMOS), Soil Moisture Active Passive (SMAP) etc.) - Thermal Infrared (TIR) Remote Sensing: Indirect SM retrieval through energy flux estimations; Relatively better spatial resolution 1-10 km; Root-zone moisture retrieval capability; Can not penetrate through clouds, hence data gaps (e.g. Surface Energy Balance Algorithm for Land (SEBAL), Atmospheric Land Exchange Inverse (ALEXI) etc.)

SERVIR↗
NASA NTRS2019

Assimilation of Satellite Soil Moisture for Improved Atmospheric Reanalyses

A newly developed, weakly coupled land and atmosphere data assimilation system for NASA’s Global Earth Observing System model is presented, and used to demonstrate the benefit of assimilating satellite soil moisture into an atmospheric reanalysis. Specifically, Advanced Scatterometer and Soil Moisture Ocean Salinity soil moisture retrievals are assimilated into a system that uses the same model, atmospheric assimilation system, and atmospheric observations as the Modern-Era Retrospective analysis for Research and Applications, version 2 (MERRA-2). The atmosphere is sensitive to soil moisture only under certain conditions. Hence, while the globally averaged model improvements were small, regionally, the soil moisture assimilation induced some substantial improvements. For example, in a large region spanning from western Europe across southern Russia, the soil moisture assimilation decreased the RMSE against independent station observations of daily maximum 2-m temperature (T(sup 2m, sub max)) by up to 0.4 K, and of 2-m specific humidity (q(exp 2m)) by up to 0.5 g/kg. Over all available stations, the mean T(sup 2m, sub max) RMSE was reduced from 2.82 to 2.79 K, while the mean q(exp 2m) RMSE was reduced from 1.25 to 1.20 g/kg. The soil moisture assimilation also reduced the mean RMSE across 29 flux tower sites from 34.2 to 32.6 Wm(exp −2) for latent heating, and from 37.7 to 36.5 Wm(exp −2) for sensible heating. For all variables evaluated, the soil moisture assimilation improved the model at monthly to seasonal, rather than daily, time scales. Based on the above experiments, it is recommended that satellite soil moisture be assimilated into future reanalyses, including the follow-on to MERRA-2.

Clara Draper↗
NASA NTRS2019

Uncertainty in Soil Moisture Retrievals: an Ensemble Approach Using SMOS L-Band Microwave Data

The uncertainty of soil moisture (SM) retrievals from satellite brightness temperature (TB) observations depends primarily on the choice of radiative transfer model (RTM) parameters, prior SM information and TB inputs. This paper studies the sensitivity of several (quasi-)operational and experimental SM retrieval products from the Soil Moisture Ocean Salinity (SMOS) mission to these choices at 11 reference sites, located in 7 watersheds across the United States (US). Different literature-based RTM parameter sets cause large biases between retrievals. Whereas typical RTM parameter sets are calibrated for SM retrievals, it is shown that a parameter set carefully optimized for TB forward modeling can also be used for retrieving SM. It is also shown that the inclusion of dynamic prior SM estimates in a Bayesian retrieval scheme can strongly improve SM retrievals, regardless of the choice of RTM parameters, and that the use of multi-angular and multi-polarization TB does not necessarily lead to superior retrievals compared to retrievals based on TB data at a single incidence angle and polarization. The second part of this paper evaluates ensemble uncertainty metrics for SM retrievals obtained by propagating a wide range of RTM parameters through the RTM. As expected for bounded variables, the spread in the ensemble SM retrievals is smallest for wet and dry SM values and highest for intermediate SM values. After removal of the strong long-term mean bias associated with the RTM parameter values for individual ensemble members, the remaining anomaly ensemble SM spread of 0.037 cu m/cu m approximates the actual time series unbiased root-mean-square-difference of 0.042 cu m/cu m between ensemble mean retrievals and in situ reference data across the reference sites. However, the temporal variability in the anomaly ensemble spread reveals higher-order biases in the retrieval error, which should be accounted for when characterizing retrieval error.

Jan Quets↗
NASA NTRS2019

Improving Surface Soil Moisture Estimates in West Africa Through GRACE Data Assimilation

Recent annual trends of precipitation and terrestrial water storage (TWS) in West Africa have been increasing over the past decade. Despite a significant impact of soil moisture on the TWS in West Africa, there is little research on the recent spatial and temporal behaviors of surface soil moisture (SSM) along with the hydrological trends and variability in West Africa. In this study, we assimilate TWS estimates from the Gravity Recovery and Climate Experiment (GRACE) mission into the Catchment Land Surface Model (CLSM) and evaluate its impacts on SSM simulations for the years, 2002–2017. The evaluation is performed using reference datasets: the African Monsoon Multidisciplinary Analysis (AMMA) in situ soil moisture observations, three currently available microwave satellite SSM observations from the Advanced Scatterometer (ASCAT), the Soil Moisture Ocean Salinity (SMOS), and the Soil Moisture Active Passive (SMAP) satellites and also the triple collocation analysis (TCA). Overall, modeled SSM shows good agreement with the reference datasets in terms of the anomaly SSM correlations. However, both modeled and ASCAT SSM are limited in their representation of the drying rates, as observed by ground observations, SMOS and SMAP estimates. Further, GRACE data assimilation results in improved SSM simulations in the humid regions with large annual TWS variability. This study demonstrates the utility of land data assimilation to inform hydrological conditions in West Africa, where soil moisture monitoring is necessary for water resource and livestock management.

West Africa↗
NASA NTRS2016

Managing and Servicing Physical Oceanographic Data at a NASA Distributed Active Archive Center

The NASA Earth Science Data Information Systems Project funds and operates 12 Distributed Active Archive Center(s) (DAAC) throughout the United States. Of these 12 centers, the Physical Oceanography DAAC (PO.DAAC) is committed to providing long term archival, distribution and stewardship for NASA physical oceanographic data, primarily derived from space-born satellite systems, but also including a growing set of recent and future in situ observations from the SPURS-1 and SPURS-2 campaigns. Notable NASA missions supported include: Seasat, TOPEX/Poseidon, NSCAT, QuikSCAT, ISS-RapidScat, Jason-1, Jason-2/OSTM, GRACE, Aquarius, GHRSST, and MODIS. The following interagency and international missions are also supported by PO.DAAC: AVHRR, Coriolis, DMSP, MetOp-A, MetOp-B, Oceansat-2. The PO.DAAC currently holds 525 datasets in public distribution, spanning the following observational parameters: sea surface temperature, sea surface salinity, ocean color, ocean surface currents, ocean surface wind speed, ocean surface wind direction, sea surface height, significant wave height, ocean water mass/thickness, and sea ice age. A hundred of these datasets are available in near-real-time. Datasets are distributed through a variety of open-source access protocols including FTP, OPeNDAP, and THREDDS. FTP will soon be phased out in favor of a recently introduced HTTPS PO.DAAC Drive interface that supports WebDAV and interoperable machine-to-machine communication. OPeNDAP supports remote data/metadata query, subset, and download. THREDDS provides the features of OPeNDAP with the additional feature of temporal aggregation. PO.DAAC also offers proprietary tools and services to further enhance the data discovery, visualization and analysis experience, including but not limited to: State of the Ocean, Web Services (data/metadata discovery and extraction), HiTIDE Level-2 subsetter, Live Access Server (LAS), Webification (w10nsci), and Rich Site Summary (RSS) Datacasting. To assist with provenance of datasets, PO.DAAC has implemented DOIs for the data it distributes so that they can be properly cited. There is a user forum and helpdesk that contains data recipes and via which users can get guidance. In summary, this presentation aims to provide a general overview of PO.DAAC’s web portal and data holdings along with a set of illustrative examples leading prospective data users into the practical utility of its tools and services.

Moroni, David F.↗
NASA NTRS2019

Integrated SMAP and SMOS Soil Moisture Observations

Soil Moisture Active Passive (SMAP) mission and the Soil Moisture and Ocean Salinity (SMOS) missions provide brightness temperature and soil moisture estimates every 2-3 days. SMAP brightness temperature observations were compared with SMOS observations at 40o incidence angle. The brightness temperatures from the two missions are close to each other but SMAP observations show a warmer TB bias (about 0.64 K: V pol and 1.14 K: H pol) as compared to SMOS. SMAP and SMOS missions use different retrieval algorithms and ancillary datasets which result in further inconsistencies between their soil moisture products. The reprocessed constant-angle SMOS brightness temperatures (SMOS-SMAP) were used in the SMAP soil moisture retrieval algorithm to develop a consistent multi-satellite product. The integrated product has an increased global revisit frequency (1 day) and period of record that is unattainable by either one of the satellites alone. The SMOS-SMAP soil moisture retrievals compared with in situ observations show a retrieval accuracy of less than 0.04 m3/m3. Results from the development and validation of the integrated soil moisture product will be presented.

SMOS↗
NASA NTRS2019

Assimilation of Satellite Soil Moisture for Improved Atmospheric Reanalyses

Satellite soil moisture observations are assimilated into NASA's Global Earth Observing System using a weakly coupled land/atmosphere data assimilation system. Specifically, Advanced Scatterometer and Soil Moisture Ocean Salinity soil moisture retrievals are assimilated with an Ensemble Kalman Filter into a system that uses the same model, atmospheric assimilation system (3DVar), and atmospheric observations as the Modern-Era Retrospective Analysis for Research and Applications, version 2 (MERRA-2). The first set of experiments investigate the impact of the soil moisture assimilation on the model states and fluxes within the cycling DA system. Comparison to independent ground-based observations from global networks of soil moisture, latent and sensible heating, and 2m temperature and specific humidity showed small, but positive, improvements in the global mean statistics. Additionally, in some regions, the soil moisture assimilation induced more substantial improvements. For all variables evaluated, the soil moisture assimilation improved the model at monthly to seasonal, rather than daily, time scales. Based on the above experiments, it is recommended that satellite soil moisture be assimilated into future reanalyses, including the follow-on to MERRA-2.

Draper, Clara↗
NASA NTRS2019

Global Assimilation of L-Band Brightness Temperature Observations from SMAP and SMOS into the Catchment Land Surface Model and Contribution to the Skill of Soil Moisture Estimates

The Soil Moisture Active Passive (SMAP) and Soil Moisture Ocean Salinity (SMOS) missions provide global observations of L-band (1.4 GHz) passive microwave brightness temperature (Tb) observations at a resolution of ~40 km every 2-3 days. These observations have been available since 2015 from SMAP and since 2010 from SMOS.Soil moisture estimates from the separate assimilation of SMAP and SMOS Tb observations into land surface models were previously shown to improve over model-only estimates, thereby demonstrating the value of assimilating L-band Tb observations for soil moisture estimation. The assimilation experiments documented in the literature do not, however, establish whether the joint assimilation of SMAP and SMOS Tbs improves the skill of the resulting soil moisture estimates beyond what can be achieved with either set of observations alone. Moreover, because the published SMAP and SMOS results used different assimilation systems and configurations and applied different evaluation data and procedures, it is unclear whether the assimilation of SMAP-only or SMOS-only Tbs results in better soil moisture skill. In this presentation, we compare the results of three separate assimilation experiments for the period from April 2015 to present. The experiments utilize the SMAP Level-4 Soil Moisture (L4_SM) algorithm, which assimilates L-band Tb observations into the NASA Catchment land surface model using a spatially distributed ensemble Kalman filter. Specifically, the three experiments presented here assimilate (i) SMAP Tbs only, as in the L4_SM product, (ii) SMOS Tbs only, after interpolation of 40° incidence angle, and (iii) both SMAP and SMOS Tbs. In all other respects, the configuration of three experiments is identical, as is the validation of the output. Preliminary results suggest that the joint assimilation of SMAP and SMOS Tbs yields the most skillful soil moisture estimates when compared to independent in situ measurements. Moreover, the skill of the SMAP-only assimilation exceeds that of the SMOS-only assimilation. The presentation provides an in-depth evaluation of the skill of the assimilation estimates vs. independent in situ and satellite measurements and in terms of statistics of the observation-minus-forecast residuals and analysis increments.

Reichle, Rolf↗
NASA NTRS2020

Improved Groundwater Table and L-Band Brightness Temperature Estimates for Northern Hemisphere Peatlands Using New Model Physics and SMOS Observations in a Global Data Assimilation Framework

There is an urgent need to include northern peatland hydrology in global Earth system models to better understand land-atmosphere interactions and sensitivities of peatland functions to climate change, and, ultimately, to improve climate change predictions. In this study, we introduced for the first time peatland-specific model physics into an assimilation scheme for L-band brightness temperature (Tb) data from the Soil Moisture Ocean Salinity (SMOS) mission to improve groundwater table estimates. We conducted two sets of model-only and data assimilation experiments using the Catchment Land Surface Model (CLSM), applying (over peatlands only) in one of them a peatland-specific adaptation (PEATCLSM). The evaluation against in-situ measurements of peatland groundwater table depth indicates the superiority of PEATCLSM model physics and additionally improved performance after assimilating SMOS Tb observations. The better performance of PEATCLSM over nearly all Northern Hemisphere peatlands is further supported by the better agreement between SMOS Tb observations and Tb estimates from the model-only and data assimilation runs. Within the data assimilation scheme, PEATCLSM reduces Tb observation-minus-forecast residuals and leads to reduced data assimilation updates of water storage components and, thus, reduced water budget imbalances in the assimilation system.

SMOS↗
NASA NTRS2021

Global Soil Water Estimates as Landslide Predictor: The Effectiveness of SMOS, SMAP, and GRACE Observations, Land Surface Simulations, and Data Assimilation

This global feasibility study assesses the potential of coarse-scale, gridded soil water estimates for the probabilistic modeling of hydrologically-triggered landslides, using Soil Moisture Ocean Salinity (SMOS), Soil Moisture Active Passive (SMAP) and Gravity Recovery and Climate Experiment (GRACE) remote sensing data, Catchment Land Surface Model (CLSM) simulations and six data products based on the assimilation of SMOS, SMAP, and/or GRACE observations into CLSM. SMOS or SMAP observations (~40-km resolution) are only available for less than 20% of the globally reported landslide events, because they are intermittent and uncertain in regions with complex terrain. GRACE terrestrial water storage estimates include 75% of the reported landslides but have coarse spatial and temporal resolutions (monthly, ~300-km). CLSM soil water simulations have the added advantage of complete spatial and temporal coverage, and are found to be able to distinguish between “stable slope” (no landslide) conditions and landslide-inducing conditions in a probabilistic way. Assimilating SMOS and/or GRACE data increases the landslide probability estimates based on soil water percentiles for the reported landslides, relative to model-only estimates at 36-km resolution for the period 2011-2016, unless the CLSM model-only soil water content is already high (≥ 50th percentile). The SMAP Level 4 data assimilation product (at 9-km resolution, period 2015-2019) more generally updates the soil water conditions towards higher landslide probabilities for the reported landslides, but is similar to model-only estimates for the majority of landslides where SMAP data cannot easily be converted to soil moisture owing to complex terrain.

Anne Felsberg↗
NASA NTRS2020

Toward P-Band Passive Microwave Sensing of Soil Moisture

Currently, near-surface soil moisture at a global scale is being provided using National Aeronautics and Space Administration’s (NASA’s) Soil Moisture Active Passive (SMAP) and European Space Agency’s (ESA’s) Soil Moisture and Ocean Salinity (SMOS) satellites, both of which utilize L-band (1.4 GHz; 21 cm wavelength ) passive microwave remote sensing techniques. However, a fundamental limitation of this technology is that the water content can only be measured for approximately the top 5-cm layer of soil moisture, and only over low-to-moderate vegetation covered areas in order to meet the 0.04 m3/m3 target accuracy, limiting its applicability. Consequently, a longer wavelength radiometer is being explored as a potential solution for measuring soil moisture in a deeper surface layer of soil and under denser vegetation. It is expected that P-band ( wavelength of 40 cm and frequency of 750 MHz) could potentially provide soil moisture information for the top 10-cm layer of soil, being one-tenth to one-quarter of the wavelength. In addition, P-band is expected to have higher soil moisture retrieval accuracy due to its reduced sensitivity to vegetation water content and surface roughness. To demonstrate the potential of P-band passive microwave soil moisture remote sensing, a short-term airborne field experiment was conducted over a center pivot irrigated farm at Cressy in Tasmania, Australia, in January 2017. First results showing a comparison of airborne P-band brightness temperature observations against airborne L-band brightness temperature observations and ground soil moisture measurements are presented. The P-band brightness temperature was found to have a similar but stronger response to soil moisture compared to L-band.

Nan Ye↗
NASA NTRS2021

Validation of Soil Moisture Data Products from the NASA SMAP Mission

The National Aeronautics and Space Administration (NASA)Soil Moisture Active Passive(SMAP) mission has been validatingits soil moisture (SM) products since the start of data production onMarch 31, 2015. Prior to launch, the mission defined a set of criteria for core validation sites (CVS) that enable the testing of the key mission SM accuracy requirement(unbiased root-mean-square error <0.04 m3/m3). Thevalidation approach also includes other (“sparse network”) in situSM measurements, satellite SM products, model-based SM products, and field experiments. Over the past six years, the SMAP SM products have been analyzed with respect to thesereference data,and the analysis approaches themselves have been scrutinizedin an effort to best understand the products’ performance. Validation of themost recent SMAP Level 2 and 3 SMretrievalproducts (R17000) shows that the L-band (1.4 GHz) radiometer-based SM record continues to meet mission requirements. The products aregenerallyconsistentwith SM retrievals from the European Space Agency (ESA)Soil Moisture Ocean Salinity mission, althoughthere aredifferencesin some regions. The high-resolution (3-km) SM retrieval product,generated by combining Copernicus Sentinel-1 data with SMAP observations,performswithin expectations. Currently, however,there is limited availability of3-kmCVSdatato support extensive validation at this spatial scale. The most recent (version 5)SMAP Level4 SMdata assimilation productprovidingsurface and root-zone SM with complete spatio-temporal coverageat 9-km resolution also meets performance requirements. The SMAP SM validation program will continue throughout the mission life; futureplans include expanding ittoforestedand high-latituderegions

SMAP↗
NASA NTRS2021

Study of a Strong L-Band RFI Source

CN006 is the name given by the European Space Agency (ESA) in their reports on radio frequency interference (RFI) to an instance of RFI in the L-band spectral window at 1.413 GHz protected for passive use only. The source of the interference is located east of the Chinese city of Hangzhou in an area with radar installations. Its effect was initially indistinguishable from that of many such RFI sources observed by the ESA's SMOS and NASA's Soil Moisture Active Passive (SMAP) radiometers over China. However, in July 2020, the level of radiation increased dramatically reaching levels exceeding 1700000 K as reported by Soil Moisture and Ocean Salinity (SMOS) becoming the strongest source observed by SMOS. This article describes an analysis of this very strong source (a radar). It provides information to provide insight into an example of a well-defined source of RFI and illustrates the powerful capability of the SMAP radiometer receiver and RFI processing incorporated in it for identifying and understanding interference.

Radio Frequency Interface↗
NASA NTRS2021

The Analysis of Using Satellite Soil Moisture Observations for Flood Detection, Evaluating Over the Thailand’s Great Flood of 2011

A flood monitoring and warning system provides critical information that can protect property and save lives. A basin-scale flood monitoring system requires an effective observation platform that offers extensive ground coverage of flood conditions, low latency, and high spatiotemporal resolution. While satellite imagery offers substantial spatial flood extent in detail due to its high spatial resolution, the coarse temporal resolution and cloud obstruction limit its near real-time application. Daily soil moisture data derived from satellite sensors at a scale of a few km can be used to monitor extreme wet surface conditions arising in flood occurrences. This study analyses the flood detection capabilities of several sources of soil moisture information, including the Soil Moisture and Ocean Salinity mission (SMOS), the Advanced Microwave Scanning Radiometer on EOS (AMSR-E), the Advanced SCATterometer (ASCAT) on MetOp, the Global Land Data Assimilation System (GLDAS), and the WaterGAP Global Hydrology Model (WGHM). In addition to soil moisture, the analysis includes measurements of surface reflectance from the Moderate Resolution Imaging Spectroradiometer (MODIS), precipitation measurements from the Tropical Rainfall Measuring Mission (TRMM), and terrestrial water storage estimates from the Gravity Recovery And Climate Experiment (GRACE) as proxies for flood inundations. The analysis was conducted over the Chao Phraya River Basin (CPB) in Thailand, where the Great Flood of 2011 led to one of the most significant economic losses in the country's history. Satellite-derived soil moisture exhibits a stronger correlation with the flood inundations than the precipitation, model-derived soil moisture, and terrestrial water storage data. SMOS soil moisture observation agrees best with the MODIS-derived flood extent/occurrence, both in terms of spatial distribution and timing, and providing approximated flood lead-time of one week or longer. A neural network constructed from SMOS and MODIS data is used to predict flood intensity/occurrence (given soil moisture input) with a predicted time window from eight days to thirty-two days. The short-term prediction (e.g., eight days) achieves the highest accuracy with an averaged recovery rate of approximately 60% (correlation coefficient). This study's results suggest a potential application of satellite soil moisture data in assisting flood monitoring and warning systems.

Satellite soil moisture↗
NASA NTRS2021

Land Surface Modeling Over the Dry Chaco: The Impact of Model Structures, and Soil, Vegetation and Land Cover Parameters

In this study, we tested the impact of a revised set of soil, vegetation and land cover parameters on the performance of three different state-of-the-art land surface models (LSMs) within the NASA Land Information System (LIS). The impact of this revision was tested over the South American Dry Chaco, an ecoregion characterized by deforestation and forest degradation since the 1980s. Most large-scale LSMs may lack the ability to correctly represent the ongoing deforestation processes in this region, because most LSMs use climatological vegetation indices and static land cover information. The default LIS parameters were revised with improved soil parameters, satellite-based interannually varying vegetation indices (leaf area index and green vegetation fraction) instead of climatological vegetation indices, and yearly land cover information instead of static land cover. A relative comparison in terms of water budget components and “efficiency space” for various baseline and revised experiments showed that large regional and long-term differences in the simulated water budget partitioning relate to different LSM structures, whereas smaller local differences resulted from updated soil, vegetation and land cover parameters. Furthermore, the different LSM structures redistributed water differently in response to these parameter updates. A time-series comparison of the simulations to independent satellite-based estimates of evapotranspiration and brightness temperature (Tb) showed that no LSM setup significantly outperformed another for the entire region and that not all LSM simulations improved with updated parameter values. However, the revised soil parameters generally reduced the bias between simulated surface soil moisture and pixel-scale in situ observations and the bias between simulated Tb and regional Soil Moisture Ocean Salinity (SMOS) observations. Our results suggest that the different hydrological responses of various LSMs to vegetation changes may need further attention to gain benefits from vegetation data assimilation.

Michiel Maertens↗
NASA NTRS2021

Study of a Strong RFI Source at L-Band Using SMAP Radiometer Data

This paper presents an analysis of Radio Frequency Interference (RFI) in the 1.400-1.427 GHz frequency band. The study considers the sudden and strong increase of interference from a particular emitter in China that has been observed in July 2020 by radiometers from both ESA's SMOS (Soil Moisture Ocean Salinity) and NASA's SMAP (Soil Moisture Active Passive) missions. It provides an example of the characterization of a source of RFI and illustrates the capabilities of the SMAP radiometer receiver and RFI processing incorporated in it to identify and understand interference.

Radio Frequency Interface↗
NASA NTRSDate not supplied

Adding GPU Support to the Markov Chain Monte Carlo Code Catmip

In geophysics, we are confronted with many under-determined inverse problems. For example, all of our observations of earthquakes are made at the Earth’s surface. So, when we try to infer how slip during an earthquake evolves in space and time, we find that there are many potential slip histories that are consistent with our limited observations and our understanding of earthquake physics. One way to approach these problems is with Bayesian analysis which allows us to infer the ensemble of all potential slip models that satisfy the observations and our prior knowledge of earthquake physics. In Bayesian analysis, our prior knowledge is known as the prior probability density function or prior PDF, the fit to the data is known as the data likelihood, and the target PDF that satisfies both the prior PDF and data likelihood is known as the posterior PDF. However, simulating the posterior PDF typically requires using Markov Chain Monte Carlo (MCMC) to draw tens of billions of random realizations of earthquake slip models, which may not be computationally feasible. To make this and similar geophysical inversions computationally tractable, we developed the Cascading Adaptive Transitional Metropolis In Parallel (CATMIP) algorithm. CATMIP is an efficient parallel Markov Chain Monte Carlo (MCMC) sampler that is used for model fitting and uncertainty quantification in geophysics. Example use cases are earthquake rupture modeling, determining mineral composition on Mars, reconstructing the history of ocean salinity, and historical earthquake relocation. CATMIP employs many parallel instances of the Metropolis algorithm for sampling in a transitioning framework. Transitioning is a process in which a set of random samples at equilibrium with a known probability density function (PDF) are used as seeds for the Markov chains to sample successive target PDFs that incrementally move the distribution from the starting seeds to the final desired PDF that describes the relative plausibility of potential values for the model parameters. The algorithm is implemented as a Master-Worker model employing MPI for communication. The worker processes are loosely coupled with global parameters periodically optimized by the master process. This provides a very high amount of parallelism with little communication between updates. During the presentation we will discuss the history of the algorithm and elaborate the earthquake rupture modeling use case for the CATMIP package. Our first step toward GPU optimization was to optimize the code for the CPU. CPU profiling revealed that most of the compute time is spent in calls to level 2 BLAS routines and calls to GSL random number generators. We revised the algorithm to employ level 3 BLAS routines instead. In our presentation we will describe how this was accomplished. Adding GPU support to CATMIP consisted mostly of replacing the calls to GSL with calls to GPU vendor-provided library routines. A small number of loops were directly implemented in CUDA. In the presentation will provide implementation details. Finally, we will discuss methods for profiling and opportunities for further optimizing GPU execution. By creating a code with the flexibility to run on either a CPU or GPU architecture, CATMIP can be used on systems ranging from large CPU-based HPC environments to single servers with GPU acceleration and everything in between.

HECC↗