Estimation of Surface Turbulent Fluxes From Land Surface Moisture and Temperature Via a Variational Data Assimilation Framework
Explore the source record for details and available documents.
SEARCH · Engineering Papers
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.
The aim of this article is to estimate surface latent heat fluxes in the vicinity of a rapidly deepening cyclone before and during its period of most rapid intensification. This is done with a bulk parameterization scheme and remotely sensed input data. A method for estimating the difference in specific humidity between the surface and a 10-m height is investigated using the Special Sensor Microwave/Imager-observed integrated water vapor field and a sea surface temperature analysis. The surface latent heat flux fields generally have estimated errors below 40 percent south of 40 deg N and outside the region of high integrated water vapor values associated with frontal bands. The method of estimating surface latent heat fluxes for the case study was found to be usable in most regions of the northwest Atlantic Ocean except for those locations directly adjacent to coastlines in instances of offshore flow and in the vicinity of surface fronts.
Thermal-based remote sensing of surface energy fluxes has traditionally relied on high spatial resolution satellite data with revisit frequencies on the order of weeks. In this study, we evaluate a biophysics-based analytical surface energy balance model for predicting latent energy (LE) and sensible heat (H) fluxes using proximal sensing observations. The Surface Temperature Initiated Closure (STIC1.2) model has been extensively validated across a wide range of spatial and temporal scales using various satellite-derived thermal datasets. Here we extend this validation by applying STIC at sub-hourly temporal resolution over multiple growing seasons for four distinct agricultural systems. We further develop and evaluate novel STIC variants that incorporate machine learning (ML) techniques to eliminate the need for specific surface energy balance observations, specifically net radiation and soil heat flux, thereby enhancing model applicability in data-sparse settings. The integration of an ML component to estimate surface available energy is shown to have strong predictive performance for both LE (R2 = 0.81-0.94) and H (R2 = 0.46-0.72) across all agricultural systems examined here, demonstrating the potential of hybrid biophysical – machine learning approaches for surface energy balance modeling with minimal data requirements. This study concludes with a novel application of explainable machine learning (exML) to diagnose sources of model error. This exML framework attributes residual prediction errors to both model input variables and environmental drivers not explicitly included in the simulation experiments. This approach provides a new pathway for improving model design and integrating previously overlooked yet influential variables into future model iterations.
The upper Colorado River basin is the primary source of water for 40 million people. With declining snowpack in the basin, forecasting hydrological budgets in the Southwest United States is more important than ever. However, due in part, to a lack of reliable observations of precipitation in complex terrain, hydrological models struggle to assess and forecast snowpack snow water equivalent (SWE) in the upper Colorado River basin (UCRB). Therefore, the need for more reliable SWE forecasts in the UCRB motivated the U.S. Department of Energy Atmospheric Radiation Measurement Facility’s Surface Atmospheric Integrated Field Laboratory (SAIL) that occurred from June 2021 to June 2023. During SAIL, the X-band precipitation radar from Colorado State University conducted volume scans sampling the precipitation properties over the UCRB. The ARM facility developed a gridded Surface Quantitative Precipitation Estimates (SQUIRE) product from the radar observations. To do this, various daily SWE estimates from radar using the radar reflectivity factor Z e and specific differential phase K dp were compared against ground-based precipitation gauges. SWE in precipitation calculated from Wolfe and Snider’s S–Z e estimator was in best agreement with the rain gauges for the days when SWE < 12 mm. For days with SWE > 12 mm, the WSR-88D Intermountain West relationship had the best agreement with the precipitation gauges. Airborne snow depth observations show that SQUIRE captures regions of orographic enhancement in the mountains to the west and northwest of the SAIL study area, indicating that the scientific community should focus on understanding and ultimately simulating orographic atmospheric precipitation processes to improve UCRB snowpack SWE assessment and forecasting.
Thermal infrared-based remote sensing of surface energy fluxes has traditionally relied on high spatial resolution satellite data with revisit frequencies on the order of weeks. In this study, we evaluate a biophysics-based analytical surface energy balance model for predicting latent energy ( LE ) and sensible heat ( H ) fluxes using proximal sensing observations. The Surface Temperature Initiated Closure (STIC1.2) model has been extensively validated across a wide range of spatial and temporal scales using various satellite-derived thermal infrared data sets. Here we extend this validation by applying STIC at sub-hourly temporal resolution over multiple growing seasons for four distinct agricultural systems. We further develop and evaluate novel STIC variants that incorporate machine learning (ML) techniques to eliminate the need for surface energy balance observations, specifically net radiation and soil heat flux, thereby enhancing model applicability in data-sparse settings. The integration of a ML component to estimate surface available energy is shown to have strong predictive performance for both LE (R 2 = 0.81–0.94) and H (R 2 = 0.46–0.72) across all agricultural systems examined here, demonstrating the potential of hybrid biophysical-machine learning approaches for surface energy balance modeling with minimal data requirements. This study concludes with a novel application of explainable machine learning (exML) to diagnose sources of model error. This exML framework attributes residual prediction errors to both model input variables and environmental drivers not explicitly included in the simulation experiments. This approach provides a new pathway for improving model design and integrating previously overlooked yet influential variables into future model iterations.
Soil moisture is an important component of the water and energy balances of the Earth's surface. Furthermore, it has been identified as a parameter of significant potential for improving the accuracy of large-scale land surface-atmosphere interaction models. However, accurate estimates of surface soil moisture are often difficult to make, especially at large spatial scales. Soil moisture is a highly variable land surface parameter, and while point measurements are usually accurate, they are representative only of the immediate site which was sampled. Simple averaging of point values to obtain spatial means often leads to substantial errors. Since remotely sensed observations are already a spatially averaged or areally integrated value, they are ideally suited for measuring land surface parameters, and as such, are a logical input to regional or larger scale land process models. A nine-year database of surface soil moisture is being developed for the Central United States from satellite microwave observations. This region forms much of the GCIP study area, and contains most of the Mississippi, Rio Grande, and Red River drainages. Daytime and nighttime microwave brightness temperatures were observed at a frequency of 6.6 GHz, by the Scanning Multichannel Microwave Radiometer (SMMR), onboard the Nimbus 7 satellite. The life of the SMMR instrument spanned from Nov. 1978 to Aug. 1987. At 6.6 GHz, the instrument provided a spatial resolution of approximately 150 km, and an orbital frequency over any pixel-sized area of about 2 daytime and 2 nighttime passes per week. Ground measurements of surface soil moisture from various locations throughout the study area are used to calibrate the microwave observations. Because ground measurements are usually only single point values, and since the time of satellite coverage does not always coincide with the ground measurements, the soil moisture data were used to calibrate a regional water balance for the top 1, 5, and 10 cm surface layers in order to interpolate daily surface moisture values. Such a climate-based approach is often more appropriate for estimating large-area spatially averaged soil moisture because meteorological data are generally more spatially representative than isolated point measurements of soil moisture. Vegetation radiative transfer characteristics, such as the canopy transmissivity, were estimated from vegetation indices such as the Normalized Difference Vegetation Index (NDVI) and the 37 GHz Microwave Polarization Difference Index (MPDI). Passive microwave remote sensing presents the greatest potential for providing regular spatially representative estimates of surface soil moisture at global scales. Real time estimates should improve weather and climate modelling efforts, while the development of historical data sets will provide necessary information for simulation and validation of long-term climate and global change studies.
Through analyses of the model simulated data-base, we developed a technique to estimate surface soil moisture under HYDROS radar sensor (L-band multi-polarizations and 40deg incidence) configuration. This technique includes two steps. First, it decomposes the total backscattering signals into two components - the surface scattering components (the bare surface backscattering signals attenuated by the overlaying vegetation layer) and the sum of the direct volume scattering components and surface-volume interaction components at different polarizations. From the model simulated data-base, our decomposition technique works quit well in estimation of the surface scattering components with RMSEs of 0.12,0.25, and 0.55 dB for VV, HH, and VH polarizations, respectively. Then, we use the decomposed surface backscattering signals to estimate the soil moisture and the combined surface roughness and vegetation attenuation correction factors with all three polarizations.
The Langley cosmic ray transport code and the Langley nucleon transport code (BRYNTRN) are used to quantify the transport and attenuation of galactic cosmic rays (GCR) and solar proton flares through the Martian atmosphere. Surface doses are estimated using both a low density and a high density carbon dioxide model of the atmosphere which, in the vertical direction, provides a total of 16 g/sq cm and 22 g/sq cm of protection, respectively. At the Mars surface during the solar minimum cycle, a blood-forming organ (BFO) dose equivalent of 10.5 to 12 rem/yr due to galactic cosmic ray transport and attenuation is calculated. Estimates of the BFO dose equivalents which would have been incurred from the three large solar flare events of August 1972, November 1960, and February 1956 are also calculated at the surface. Results indicate surface BFO dose equivalents of approximately 2 to 5, 5 to 7, and 8 to 10 rem per event, respectively. Doses are also estimated at altitudes up to 12 km above the Martian surface where the atmosphere will provide less total protection.
Regional to continental scale magnetic anomaly maps are becoming increasingly available from airborne, shipborne, and terrestrial surveys. Satellite data are commonly considered to fill the coverage gaps in regional compilations of these near-surface surveys. For the near-surface Antarctic magnetic anomaly map being produced by the Antarctic Digital Magnetic Anomaly Project (ADMAP), we show that near-surface magnetic anomaly estimation is greatly enhanced by the joint inversion of the near-surface data with the satellite observations relative to the conventional technique such as minimum curvature. Orsted observations are especially advantageous relative to the Magsat data that have order-of-magnitude greater measurement errors, albeit at much lower orbital altitudes. CHAMP is observing the geomagnetic field with the same measurement accuracy as the Orsted mission, but at the lower orbital altitudes covered by Magsat. Hence, additional significant improvement in predicting near-surface magnetic anomalies can result as these CHAMP data are available. Our analysis also suggests that considerable new insights on the magnetic properties of the lithosphere may be revealed by a further order-of-magnitude improvement in the accuracy of the magnetometer measurements at minimum orbital altitude.
The overall goal of this project was to increase our understanding of processes which determine the temporally varying distributions of surface chlorophyll pigment concentration and surface temperature in the California Current System (CCS) on the time-scale of 'events', i.e., several days to several weeks. We also proposed to investigate seasonal and regional differences in these events. Additionally, we proposed to evaluate methods of estimating surface velocities and horizontal transport of pigment and heat from sequences of AVHRR and CZCS images. The four specific objectives stated in the original proposal were to: (1) test surface current estimates made from sequences of both SST and color images using variations of the statistical method of Emery et al. (1986) and estimate the uncertainties in these satellite-derived surface currents; (2) characterize the spatial and temporal relationships of chlorophyll and temperature in rapidly evolving features for which adequate imagery exist and evaluate the contribution of these events to monthly and seasonal averages; (3) use the methods tested in (1) to determine the nature of the velocity fields in the CCS; and (4) compare the currents, temperature, and currents in different seasons and in different geographic regions.
Space-borne microwave sensors provide critical rain information used in several global multi-satellite rain products, which in turn are used for a variety of important studies, including landslide forecasting, flash flood warning, data assimilation, climate studies, and validation of model forecasts of precipitation. This study employs four years (2003-2006) of satellite data to assess the relative performance and skill of SSM/I (F13, F14 and F15), AMSU-B (N15, N16 and N17), AMSR-E (Aqua) and the TRMM Microwave Imager (TMI) in estimating surface rainfall based on direct instantaneous comparisons with ground-based rain estimates from Tropical Rainfall Measuring Mission (TRMM) Ground Validation (GV) sites at Kwajalein, Republic of the Marshall Islands (KWAJ) and Melbourne, Florida (MELB). The relative performance of each of these satellite estimates is examined via comparisons with space- and time-coincident GV radar-based rain rate estimates. Because underlying surface terrain is known to affect the relative performance of the satellite algorithms, the data for MELB was further stratified into ocean, land and coast categories using a 0.25 terrain mask. Of all the satellite estimates compared in this study, TMI and AMSR-E exhibited considerably higher correlations and skills in estimating/observing surface precipitation. While SSM/I and AMSU-B exhibited lower correlations and skills for each of the different terrain categories, the SSM/I absolute biases trended slightly lower than AMSRE over ocean, where the observations from both emission and scattering channels were used in the retrievals. AMSU-B exhibited the least skill relative to GV in all of the relevant statistical categories, and an anomalous spike was observed in the probability distribution functions near 1.0 mm/hr. This statistical artifact appears to be related to attempts by algorithm developers to include some lighter rain rates, not easily detectable by its scatter-only frequencies. AMSU-B, however, agreed well with GV when the matching data was analyzed on monthly scales. These results signal developers of global rainfall products, such as the TRMM Multi-Satellite Precipitation Analysis (TMPA), and the Climate Data Center s Morphing (CMORPH) technique, that care must be taken when incorporating data from these input satellite estimates in order to provide the highest quality estimates in their products. 3
Space-borne microwave sensors provide critical rain information used in several global multi-satellite rain products, which in turn are used for a variety of important studies, including landslide forecasting, flash flood warning, data assimilation, climate studies, and validation of model forecast of precipitation. This study employs four years (2003-2006) of satellite data to assess the relative performance and skill of SSM/I (F13, F14 and F15), AMSU-B (N15, N16 and N17), AMSR-E (AQUA) and the TRMM Microwave Imager (TMI) in estimating surface rainfall based on direct instantaneous comparison with ground-based rain estimates from Tropical Rainfall Measuring Mission (TRMM) Ground Validation (GV) sites at Kwajalein, Republic of the Marshall Islands (KWAJ) and Melbourne, Florida (MELB). The relative performance of each of these satellites is examined via comparisons with GV radar-based rain rate estimates. Because underlying surface terrain is known to affect the relative performance of the satellite algorithms, the data for MELB was further stratified into ocean, land and coast categories using a 0.25 terrain mask. Of all the satellite estimates compared in this study, TMI and AMSR-E exhibited considerably higher correlations and skills in estimating/observing surface precipitation. While SSM/I and AMSU-B exhibited lower correlations and skills for each of the different terrain categories, the SSM/I absolute biases trended slightly lower than AMSRE over ocean, where the observations from both emission and scattering channels were used in the retrievals. AMSU-B exhibited the least skill relative to GV in all of the relevant statistical categories, and an anomalous spike was observed in the probability distribution functions near 1.0 mm hr-1. This statistical artifact appears to be related to attempts by algorithm developers to include some lighter rain rates, not easily detectable by its scatter-only frequencies. AMSU-B, however, agreed well with GV when the matching data was analyzed on monthly scales. These results signal developers of global rainfall products, such as the TRMM Multi-Satellite Precipitation Analysis (TMPA), and the Climate Data Center s Morphing (CMORPH) technique, that care must be taken when incorporating data from these input satellite estimates in order to provide the highest quality estimates in their products.
Space-borne microwave sensors provide critical rain information used in several global multi-satellite rain products, which in turn are used for a variety of important studies, including landslide forecasting, flash flood warning, data assimilation, climate studies, and validation of model forecasts of precipitation. This study employs four years (2003-2006) of satellite data to assess the relative performance and skill of SSM/I (F13, F14 and F15), AMSU-B (N15, N16 and N17), AMSR-E (Aqua) and the TRMM Microwave Imager (TMI) in estimating surface rainfall based on direct instantaneous comparisons with ground-based rain estimates from Tropical Rainfall Measuring Mission (TRMM) Ground Validation (GV) sites at Kwajalein, Republic of the Marshall Islands (KWAJ) and Melbourne, Florida (MELB). The relative performance of each of these satellite estimates is examined via comparisons with space- and time-coincident GV radar-based rain rate estimates. Because underlying surface terrain is known to affect the relative performance of the satellite algorithms, the data for MELB was further stratified into ocean, land and coast categories using a 0.25deg terrain mask. Of all the satellite estimates compared in this study, TMI and AMSR-E exhibited considerably higher correlations and skills in estimating/observing surface precipitation. While SSM/I and AMSU-B exhibited lower correlations and skills for each of the different terrain categories, the SSM/I absolute biases trended slightly lower than AMSR-E over ocean, where the observations from both emission and scattering channels were used in the retrievals. AMSU-B exhibited the least skill relative to GV in all of the relevant statistical categories, and an anomalous spike was observed in the probability distribution functions near 1.0 mm/hr. This statistical artifact appears to be related to attempts by algorithm developers to include some lighter rain rates, not easily detectable by its scatter-only frequencies. AMSU-B, however, agreed well with GV when the matching data was analyzed on monthly scales. These results signal developers of global rainfall products, such as the TRMM Multi-Satellite Precipitation Analysis (TMPA), and the Climate Data Center s Morphing (CMORPH) technique, that care must be taken when incorporating data from these input satellite estimates in order to provide the highest quality estimates in their products.
The Hurricane Imaging Radiometer (HIRAD) is an experimental C-band passive microwave radiometer designed to map the horizontal structure of surface wind speed fields in hurricanes. New data processing and customized retrieval approaches were developed after the 2015 Tropical Cyclone Intensity (TCI) experiment, which featured flights over Hurricanes Patricia, Joaquin, Marty, and the remnants of Tropical Storm Erika. These new approaches produced maps of surface wind speed that looked more realistic than those from previous campaigns. Dropsondes from the High Definition Sounding System (HDSS) that was flown with HIRAD on a WB-57 high altitude aircraft in TCI were used to assess the quality of the HIRAD wind speed retrievals. The root mean square difference between HIRAD-retrieved surface wind speeds and dropsonde-estimated surface wind speeds was 6.0 meters per second. The largest differences between HIRAD and dropsonde winds were from data points where storm motion during dropsonde descent compromised the validity of the comparisons. Accounting for this and for uncertainty in the dropsonde measurements themselves, we estimate the root mean square error for the HIRAD retrievals as around 4.7 meters per second. Prior to the 2015 TCI experiment, HIRAD had previously flown on the WB-57 for missions across Hurricanes Gonzalo (2014), Earl (2010), and Karl (2010). Configuration of the instrument was not identical to the 2015 flights, but the methods devised after the 2015 flights may be applied to that previous data in an attempt to improve retrievals from those cases.
Two components of fault slip are uniquely determined from two line-of-sight (LOS) radar interferograms by assuming that the fault-normal component of displacement is zero. We use this approach with ascending and descending interferograms from the ERS satellites to estimate surface slip along the Hector Mine earthquake rupture. The LOS displacement is determined by visually counting fringes to within 1 km of the outboard ruptures. These LOS estimates and uncertainties are then transformed into strike- and dip-slip estimates and uncertainties; the transformation is singular for a N-S oriented fault and optimal for an E-W oriented fault. In contrast to our previous strike-slip estimates, which were based only on a descending interferogram, we now find good agreement with the geological measurements, except at the ends of the rupture. The ascending interferogram reveals significant west-sidedown dip-slip (approximately 1.0 m) which reduces the strike-slip estimates by 1 to 2 m, especially along the northern half of the rupture. A spike in the strike-slip displacement of 6 m is observed in central part of the rupture. This large offset is confirmed by subpixel cross correlation of features in the before and after amplitude images. In addition to strike slip and dip slip, we identify uplift and subsidence along the fault, related to the restraining and releasing bends in the fault trace, respectively. Our main conclusion is that at least two look directions are required for accurate estimates of surface slip even along a pure strike-slip fault. Models and results based only on a single look direction could have major errors. Our new estimates of strike slip and dip slip along the rupture provide a boundary condition for dislocation modeling. A simple model, which has uniform slip to a depth of 12 km, shows good agreement with the observed ascending and descending interferograms.
Aerosol optical depth (AOD), an indirect estimate of particle matter using satellite observations, has shown great promise in improving estimates of PM 2.5 air quality surface. Currently, few studies have been conducted to explore the optimal way to apply AOD data to improve the model accuracy of PM 2.5 surface estimation in a real-time air quality system. We believe that two major aspects may be worthy of consideration in that area: 1) the approach to integrate satellite measurements with ground measurements in the pollution estimation, and 2) identification of an optimal temporal scale to calculate the correlation of AOD and ground measurements. This paper is focused on the second aspect on the identifying the optimal temporal scale to correlate AOD with PM2.5. Five following different temporal scales were chosen to evaluate their impact on the model performance: 1) within the last 3 days, 2) within the last 10 days, 3) within the last 30 days, 4) within the last 90 days, and 5) the time period with the highest correlation in a year. The model performance is evaluated for its accuracy, bias, and errors based on the following selected statistics: the Mean Bias, the Normalized Mean Bias, the Root Mean Square Error, Normalized Mean Error, and the Index of Agreement. This research shows that the model with the temporal scale of within the last 30 days displays the best model performance in this study area using 2004 and 2005 data sets.
Column abundance estimates of sodium, and analogously, potassium, in Mercury's exosphere are strongly correlated to the surface reflection model used to calibrate the spectral data and the surface reflection model incorporated into the atmospheric radiative transfer solution. Depending on the surface reflection model parameters used, there can be differences in calibration factors of up to +/- 30% and differences in estimated column abundance of up to +/- 35%. Although the surface reflectance may not be used in the calibration of spacecraft measurements, the interaction between the reflected surface light and the atmospheric brightness remains important.
We provide the first geometrically accurate (i.e., 3-D) temperature maps of the entire windward surface of the Space Shuttle during hypersonic reentry. To accomplish this task we began with estimated surface temperatures derived from CFD models at integral high Mach numbers and used them, the Shuttle's surface properties and reasonable estimates of the sensor-to-target geometry to predict the emitted spectral radiance from the surface (in units of W sr-1 m-2 nm-1). These data were converted to sensor counts using properties of the sensor (e.g. aperture, spectral band, and various efficiencies), the expected background, and the atmosphere transmission to inform the optimal settings for the near-infrared and midwave IR cameras on the Cast Glance aircraft. Once these data were collected, calibrated, edited, registered and co-added we formed both 2-D maps of the scene in the above units and 3-D maps of the bottom surface in temperature that could be compared with not only the initial inputs but also thermocouple data from the Shuttle itself. The 3-D temperature mapping process was based on the initial radiance modeling process. Here temperatures were guessed for each node in a well-resolved 3-D framework, a radiance model was produced and compared to the processed imagery, and corrections to the temperature were estimated until the iterative process converged. This process did very well in characterizing the temperature structure of the large asymmetric boundary layer transition the covered much of the starboard bottom surface of STS-119 Discovery. Both internally estimated accuracies and differences with CFD models and thermocouple measurements are at most a few percent. The technique did less well characterizing the temperature structure of the turbulent wedge behind the trip due to limitations in understanding the true sensor resolution. (Note: Those less inclined to read the entire paper are encouraged to read an Executive Summary provided at the end.)