Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Weight function”

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 307 records · Page 17

Validation of the Radiometric Stability of the Atmospheric Infrared Sounder

It has been widely accepted that an infrared sounder in low polar orbit is capable of producing climate quality data, if the spectral brightness temperatures have instrumental trends of less than 10 mK/yr. Achieving measurement stability at this level is not only very demanding of the design of the instrument, it is also pushes the state of art of measuring on orbit what stability is actually achieved. We discuss this using Atmospheric Infrared Sounder (AIRS) L1B data collected between 2002 and 2011. We compare the L1B brightness temperature observed in cloud filtered night tropical ocean spectra (obs) to the brightness temperature calculated based on the known surface emissivity, temperature and water vapor profiles from the ECMWF ReAnalysis (ERA) and the growth rates of CO2, N2O and Ozone. The trend in (obs-calc) is a powerful tool for the evaluation of the stability of the 2378 AIRS channels. We divided the channels into seven classes: All channels which sound in the stratosphere (at pressure levels below 150 hPa), 14 micron CO2 sounding, 4 micron CO2 P-branch sounding, 4 micron CO2 R-branch sounding, water vapor sounding, shortwave surface sounding and longwave surface sounding. The peak in the weighting function at 1050 hPa separates sounding and surface channels. The boundary between shortwave and longwave is 5 microns. Except for the stratosphere sounding channels, the remaining six groups have (obs-calc) trends of less than 20 mK/yr. The longwave surface channels have trends of 2 mK/yr, significantly less than the 8 mK/yr trend seem in the shortwave window channels. Based on the design of the instrument, trends within a group of channels should be the same. While the longwave and shortwave trends are less than the canonical 10 mK/yr, the larger trend in the shortwave channels could be an artifact of using the pre-launch determined calibration coefficients. This is currently under evaluation. The trend in (obs-calc) for the non-surface sounding channels, in particular for stratosphere sounding and upper tropospheric water channels, is dominated by artifacts created in calc, most likely due to changes in the ERA Ozone and water vapor. Based on this argument the best estimate of the trend for the channels within a channel group is given by the surface sensitive channels within the group. Based on this consideration we estimated the trend of all AIRS longwave channels as 2 mK/yr, while the shortwave channels have a trend of 8 mK/yr.

CO2↗

A Physically Based Algorithm for Non-Blackbody Correction of Cloud-Top Temperature and Application to Convection Study

Cloud-top temperature (CTT) is an important parameter for convective clouds and is usually different from the 11-micrometers brightness temperature due to non-blackbody effects. This paper presents an algorithm for estimating convective CTT by using simultaneous passive [Moderate Resolution Imaging Spectroradiometer (MODIS)] and active [CloudSat 1 Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO)] measurements of clouds to correct for the non-blackbody effect. To do this, a weighting function of the MODIS 11-micrometers band is explicitly calculated by feeding cloud hydrometer profiles from CloudSat and CALIPSO retrievals and temperature and humidity profiles based on ECMWF analyses into a radiation transfer model.Among 16 837 tropical deep convective clouds observed by CloudSat in 2008, the averaged effective emission level (EEL) of the 11-mm channel is located at optical depth; approximately 0.72, with a standard deviation of 0.3. The distance between the EEL and cloud-top height determined by CloudSat is shown to be related to a parameter called cloud-top fuzziness (CTF), defined as the vertical separation between 230 and 10 dBZ of CloudSat radar reflectivity. On the basis of these findings a relationship is then developed between the CTF and the difference between MODIS 11-micrometers brightness temperature and physical CTT, the latter being the non-blackbody correction of CTT. Correction of the non-blackbody effect of CTT is applied to analyze convective cloud-top buoyancy. With this correction, about 70% of the convective cores observed by CloudSat in the height range of 6-10 km have positive buoyancy near cloud top, meaning clouds are still growing vertically, although their final fate cannot be determined by snapshot observations.

Physical Algorithm↗

A Versatile Nonlinear Method for Predictive Modeling

As computational fluid dynamics techniques and tools become widely accepted for realworld practice today, it is intriguing to ask: what areas can it be utilized to its potential in the future. Some promising areas include design optimization and exploration of fluid dynamics phenomena (the concept of numerical wind tunnel), in which both have the common feature where some parameters are varied repeatedly and the computation can be costly. We are especially interested in the need for an accurate and efficient approach for handling these applications: (1) capturing complex nonlinear dynamics inherent in a system under consideration and (2) versatility (robustness) to encompass a range of parametric variations. In our previous paper, we proposed to use first-order Taylor expansion collected at numerous sampling points along a trajectory and assembled together via nonlinear weighting functions. The validity and performance of this approach was demonstrated for a number of problems with a vastly different input functions. In this study, we are especially interested in enhancing the method's accuracy; we extend it to include the second-orer Taylor expansion, which however requires a complicated evaluation of Hessian matrices for a system of equations, like in fluid dynamics. We propose a method to avoid these Hessian matrices, while maintaining the accuracy. Results based on the method are presented to confirm its validity.

nonlinear modeling↗

Triple-Pulsed Two-Micron Integrated Path Differential Absorption Lidar: A New Active Remote Sensing Capability with Path to Space

The two-micron wavelength is suitable for monitoring atmospheric water vapor and carbon dioxide, the two most dominant greenhouse gases. Recent advances in 2-micron laser technology paved the way for constructing state-of-the-art lidar transmitters for active remote sensing applications. In this paper, a new triple-pulsed 2-micron integrated path differential absorption lidar is presented. This lidar is capable of measuring either two species or single specie with two different weighting functions, simultaneously and independently. Development of this instrument is conducted at NASA Langley Research Center. Instrument scaling for projected future space missions will be discussed.

Singh, Upendra N.↗

Estimating the Distribution of Colored Dissolved Organic Matter During the Southern Ocean Gas Exchange Experiment Using Four-Dimensional Variational Data Assimilation

We diagnosed the effect of various physical processes on the distribution of mixed-layer colored dissolved organic matter (CDOM) and a sulfur hexauoride (SF6) tracer during the Southern Ocean Gas Exchange Experiment (SO GasEx). The biochemical upper ocean state estimate uses in situ and satellite biochemical and physical data in the study region, including CDOM (absorption coefcient and spectral slope), SF6, hydrography, and sea level anomaly. Modules for photobleaching of CDOM and surface transport of SF6 were coupled with an ocean circulation model for this purpose. The observed spatial and temporal variations in CDOM were captured by the state estimate without including any new biological source term for CDOM, assuming it to be negligible over the 26 days of the state estimate. Thermocline entrainment and photobleaching acted to diminish the mixed-layer CDOM with time scales of 18 and 16 days, respectively. Lateral advection of CDOM played a dominant role and increased the mixed-layer CDOM with a time scale of 12 days, whereas lateral diffusion of CDOM was negligible. A Lagrangian view on the CDOM variability was demonstrated by using the SF6 as a weighting function to integrate the CDOM elds. This and similar data assimilation methods can be used to provide reasonable estimates of optical properties, and other physical parameters over the short-term duration of a research cruise, and help in the tracking of tracer releases in large-scale oceanographic experiments, and in oceanographic process studies.

Del Castillo, C. E.↗

ATMS Radiometric Noise Characterization

The Advanced Technology Microwave Sounder (ATMS), on SNPP and NOAA-20, provides vertical atmospheric temperature and moisture profiles that are a key input to numerical weather prediction. Full characterization of radiometric noise is critical to error modeling, algorithm development and the potential development of performance requirements and product improvements. This paper presents estimates of white noise and low-frequency 1/f noise as derived from the Allan variance and scan-to-scan autocorrelations of data from the internal hot calibration. This approach allows direct comparison of on-orbit performance of the two units as well as comparison between ground calibration test data and on-orbit performance. The autocorrelation at multiple scan intervals is used as input to a theoretical model of the NET of inferred scene brightness temperature, when using weighted samples of calibration data over multiple scans. This will support selection of optimized calibration weighting functions for each channel, as well as providing further insight into the striping phenomenon. These measured autocorrelations are also compared to the autocorrelations derived from the 1/f spectra obtained both in special ground and on-orbit tests that characterize gain fluctuation. Channel cross-correlations are also presented, which can have a significant impact on effectiveness of direct assimilation performance.

instrument calibration↗

Needs and Opportunities in the Development of Advanced Materials and Manufacturing Methods for Future Long-Duration Human Space Exploration

Weight, functionality, and sustainability are all critical concerns for future, long-duration, human exploration of space. Exploration missions will be mass-limited, since the amount of supplies and instruments that future astronauts will be able to take with them will be limited by launch vehicle and spacecraft mass and efficiency. Astronauts will need to have the capability to repair or replace worn-out or broke hardware and produce new components to be able to function for long periods at locations far-removed from Earth. Ultra-lightweight, multifunctional materials will be required to enable significant reductions in launch vehicle, spacecraft, and habitat mass in order to maximize payload. "Mass-less Exploration" concepts must be developed that will recycle materials and waste and convert available planetary materials into new feedstock materials and utilize in-space additive manufacturing to use these materials to produce the

composites↗

Reducing Errors in Velocity–Azimuth Display (VAD) Wind and Deformation Retrievals from Airborne Doppler Radars in Convective Environments

The present study describes methods to reduce the uncertainty of velocity–azimuth display (VAD) wind and deformation retrievals from downward-pointing, conically scanning, airborne Doppler radars. These retrievals have important applications in data assimilation and real-time data processing. Several error sources for VAD retrievals are considered here, including violations to the underlying wind field assumptions, Doppler velocity noise, data gaps, temporal variability, and the spatial weighting function of the VAD retrieval. Specific to airborne VAD retrievals, we also consider errors produced due to the radar scans occurring while the instrument platform is in motion. While VAD retrievals are typically performed using data from a single antenna revolution, other strategies for selecting data can be used to reduce retrieval errors. Four such data selection strategies for airborne VAD retrievals are evaluated here with respect to their effects on the errors. These methods are evaluated using the second hurricane nature run numerical simulation, analytic wind fields, and observed Doppler radar radial velocities. The proposed methods are shown to reduce the median absolute error of the VAD wind retrievals, especially in the vicinity of deep convection embedded in stratiform precipitation. The median absolute error due to wind field assumption violations for the along-track and for the across-track wind is reduced from 0.36 to 0.08 m s−1 and from 0.35 to 0.24 m s−1, respectively. Although the study focuses on Doppler radars, the results are equally applicable to conically scanning Doppler lidars as well.

Charles N Helms↗

Thermal Hydraulic Disaggregation of SMAP Soil Moisture Over the Continental United States

Thermal Hydraulic disaggregation of Soil Moisture (THySM) algorithm was implemented to downscale NASA’s Soil Moisture Active Passive (SMAP) Enhanced soil moisture (SM) product to 1 km over the continental United States (CONUS). This algorithm was developed by combining thermal inertia theory with a soil hydraulic-based approach that considers fine-scale SM spatial distribution driven by both heat fluxes and hydraulic conductivity in soils. Relative soil wetness values were estimated using land surface temperature and normalized difference vegetation index for the thermal inertia model and using soil properties for the hydraulic model. The relative soil wetness values at 1 km from both models were then combined by using weighting functions whereby the spatial distribution of SM was governed more by thermal fluxes during times of strong heat transport and infiltration during moisture abundant soil conditions. THySM values were evaluated using in situ SM measurements from SMAP Core Validation Sites (CVS), the USDA Soil Climate Analysis Network, and the NOAA Climate Reference Network over CONUS. THySM shows higher accuracy than the SMAP / Sentinel-1 (SPL2SMAP_S) 1 km SM product when compared to in situ measurements. The accuracy of THySM is 0.048 m3/m3 based on unbiased root mean square error (ubRMSE), outperforming SPL2SMAP_S by 0.01-0.02 m3/m3. The ubRMSE of THySM 1km SM over the SMAP grassland/rangeland-dominated CVS sites is better than 0.04 m3/m3, which meets the SMAP mission SM accuracy requirement applied at 9 and 36 km.

SMAP↗

Validations of Satellite Ozone Profiles in Austral Spring Using Ozonesonde Measurements in the Jang Bogo Station, Antarctica

Using ozonesonde measurements from 2015 to 2018 at the Jang Bogo station located in the southeastern Antarctic region, we evaluate ozone profiles retrieved from the three satellite measurements that are widely used: Ozone Monitoring Instrument (OMI), Microwave Limb Sounder (MLS), and Ozone Mapping Profiler Suite (OMPS) data. For the fair validation, ozonesonde profiles are smoothed using the weighting function of each satellite retrieval algorithm (i.e., convolution process). Compared with limb-viewing MLS and OMPS ozone profiles, the OMI ozone profiles are relatively less qualified: coarser vertical resolution and larger inter-annual variation. Nevertheless, our validation reveals that the quality of all three satellite ozone profiles looks comparable; In general, difference from ozonesonde profile is ∼1 ppm absolutely, and −20 to 30% relatively at maximum. This quantitative range well corresponds to previous work, meaning that our new validation confirms the reliability of satellite ozone profiles in the southeastern Antarctic region where the measurement data for the validation were not enough. Another interesting feature is the role of a priori ozone profile; Nadir-viewing OMI satellite can have qualified ozone profiles by a proper assumption of a priori ozone profile. Since the performance of limb-viewing ozone profiles is better, however, the careful usage of nadir-viewing ozone profile is still required. We think that the simultaneous usage of multiple satellite ozone profiles can contribute to better understanding of Antarctic ozone characteristics.

Antarctica↗

Implementation of A New Microwave Scattering Database and A Forward Model for Active Microwave Sensors in CRTM

Radiative transfer models are extensively used for the assimilation of satellite observations into NWP models as well as retrieving geophysical products from satellite measurements. CRTM is a community model developed by NOAA JCSDA and widely used for different purposes requiring RT calculations. CRTM requires bulk scattering lookup tables in order to perform all-sky RT calculations. However, the current CRTM lookup tables for microwave frequencies were generated based on the Mie theory by assuming spherical frozen particles. The scattering lookup tables generated using the DDA technique has shown to largely improve the RT scattering calculations in the MW region. This presentation targets (i) the implementation and validation of a DDA database that was originally developed for the ARTS RT model into CRTM, and (ii) developing the CRTM active sensor module that takes advantage of the backscattering coefficients computed using the DDA method. The DDA database only provides single scattering properties of different habits, while CRTM requires bulk scattering properties. The CRTM cloud coefficients were previously generated based on the effective radius for representing the size of the particles. However, effective radius is neither measurable nor provided by the NWP models, thus need to be estimated from other geophysical variables such as water content. Therefore, in addition to calculating the CRTM bulk scattering properties from the DDA single scattering database, the CRTM was also largely modified to use cloud water content (kg.m-3), instead of effective radius, for performing the interpolation over size/mass of the particles. CRTM already requires water content as input, thus no extra variables are required for performing scattering calculations using the new ARTS DDA database. The CRTM scattering modules search for effective radius in cloud coefficient files and will use the cloud water content if the effective radius dimension is not found in the cloud coefficient files. Figure 1 shows the CRTM simulated brightness temperatures computed using different cloud coefficients versus ATMS observed values over Hurricane Irma on September 7, 2017 at 18:00 UTC. We used all the cloud water content values included in ERA5 with default CRTM/DDA habits for water, rain, snow, ice, hail, and graupel. ERA5 does not provide separate water content values for ice, hail, and graupel, thus the ice water content values were divided between ice, hail, and graupel clouds similar to what was explained in the previous section. In channels with a frequency lower than 90 GHz, emission from water and rain clouds can compensate for cloud scattering so that cloud contaminated Tbs are larger than corresponding clear sky Tbs. The DDA simulations for channels 1-7 largely perform better than the Mie simulations. The DDA simulations show a mix of small negative and positive simulated minus observed values, while the Mie results show large negative biases. The weighting functions for some of the ATMS temperature sounding channels (channels 9-15) peak mostly above the clouds, therefore the measured Tbs become less sensitive to clouds so that the results of both Mie and DDA become very similar. The Mie lookup tables generate excessive scattering for channel 16, but not enough scattering for the water vapor channels. In the specific case of Hurricane Maria, the DDA lookup tables do not generate enough scattering for channel 16, but the DDA results are much more consistent with observations for water vapor channels than for channel 16. It should be noted that the results may vary if we use other habits to represent snow, hail, and graupel in the DDA simulations. Although these results clearly show the advantage of the DDA database over the Mie dataset, different error sources such as error in the observations, displacement of clouds in the ERA5 reanalysis, and also lack of convective clouds or in general errors in the input atmospheric and cloud profiles contribute to the differences between the simulated and observed values. Aside from the improvements in the simulations, a major advantage of the new dataset is a large number of habits that can be used to tune the data assimilation systems to perform well in different weather conditions.

Isaac Moradi↗

Optimal Estimation Inversion of Ionospheric Electron Density from GNSS-POD Limb Measurements: Part I-Algorithm and Morphology

GNSS-LEO radio links from Precise Orbital Determination (POD) and Radio Occultation (RO) antennas have been used increasingly in characterizing the global 3D distribution and variability of ionospheric electron density (N e ). In this study, we developed an optimal estimation (OE) method to retrieve N e profiles from the slant total electron content (hTEC) measurements acquired by the GNSS-POD links at negative elevation angles (ε < 0°). Although both OE and onion-peeling (OP) methods use the Abel weighting function in the N e inversion, they are significantly different in terms of performance in the lower ionosphere. The new OE results can overcome the large N e oscillations, sometimes negative values, seen in the OP retrievals in the E-region ionosphere. In the companion paper in this Special Issue, the HmF2 and NmF2 from the OE retrieval are validated against ground-based ionosondes and radar observations, showing generally good agreements in NmF2 from all sites. Nighttime hmF2 measurements tend to agree better than the daytime when the ionosonde heights tend to be slightly lower. The OE algorithm has been applied to all GNSS-POD data acquired from the COSMIC-1 (2006–2019), COSMIC-2 (2019–present), and Spire (2019–present) constellations, showing a consistent ionospheric N e morphology. The unprecedented spatiotemporal sampling of the ionosphere from these constellations now allows a detailed analysis of the frequency–wavenumber spectra for the N e variability at different heights. In the lower ionosphere (~150 km), we found significant spectral power in DE1, DW6, DW4, SW5, and SE4 wave components, in addition to well-known DW1, SW2, and DE3 waves. In the upper ionosphere (~450 km), additional wave components are still present, including DE4, DW4, DW6, SE4, and SW4. The co-existence of eastward- and westward-propagating wave4 components implies the presence of a stationary wave4 (SPW4), as suggested by other earlier studies. Further improvements to the OE method are proposed, including a tomographic inversion technique that leverages the asymmetric sampling about the tangent point associated with GNSS-LEO links.

GNSS constellation↗

Optimal Estimation Inversion of Ionospheric Electron Density from GNSS-POD Limb Measurements: Part II-Validation and Comparison Using NmF2 and hmF2

A growing number of SmallSat/CubeSat constellations with high-rate (50–100 Hz) global navigation satellite system radio occultations (GNSS-RO) as well as low-rate (1 Hz) precise orbit determination (GNSS-POD) limb-viewing capabilities provide unprecedented spatial and temporal sampling rates for ionospheric studies. In the F-region electron density (N 𝑒 ) retrieval process, instead of the conventional onion-peeling (OP) inversion, an optimal estimation (OE) inversion technique was recently developed using total electron content measurements acquired by GNSS-POD link. The new technique is applied to data acquired from the COSMIC-1, COSMIC-2, and Spire constellations. Although both OE and OP techniques use the Abel weighting function in N 𝑒 inversion, OE significantly differs in its performance, especially in the lower F- and E-regions. In this work, we evaluate and compare newly derived data sets using F2 peak properties with other space-based and ground-based observations. We determine the F2 peak N 𝑒 (NmF2) and its altitude (hmF2), and compare them with the OP-retrieved values. Good agreement is observed between the two techniques for both NmF2 and hmF2. In addition, we also utilize autoscaled F2 peak measurements from a number of worldwide Digisonde stations (∼30). The diurnal sensitivity and latitudinal variability of the F2 peak between the two techniques are carefully studied at these locations. Good agreement is observed between OE-retrieved NmF2 and Digisonde-measured NmF2. However, significant differences appear between OE-retrieved hmF2 and Digisonde-measured hmF2. During the daytime, Digisonde-measured hmF2 remains ∼25–45 km below the OE-retrieved hmF2, especially at mid and high latitudes. We also incorporate F-region N 𝑒 measurements from two incoherent scatter radar observations at high latitudes, located in the North American (Millstone Hill) and European (EISCAT at Tromso) sectors. The radar measurements show good agreement with OE-retrieved values. Although there are several possible sources of error in the ionogram-derived N 𝑒 profiles, our further analysis on F1 and F2 layers indicates that the low Digisonde hmF2 is caused by the autoscaled method, which tends to detect a height systematically below the F2 peak when the F1 layer is present.

GNSS radio occultation↗

Hybrid Temperature-Reactivity PID Controllers for Nuclear Thermal Propulsion Startup

In this work, a novel temperature-reactivity controller is introduced to rotate control drums based on a chamber temperature signal. It consists of two components: (1) a reactivity-driven proportional controller converting the temperature demand into a reactivity signal, and (2) a temperature-driven proportional-derivative controller which relies on a derivative term to anticipate sudden changes in demand and limit temperature overshoot. A weighting function is proposed to continuously transition from one component to the other.

33 ADVANCED PROPULSION SYSTEMS↗

Hybrid Temperature-Reactivity PID Controllers for Nuclear Thermal Propulsion Startup

In this work, a novel temperature-reactivity controller is introduced to rotate control drums based on a chamber temperature signal. It consists of two components: (1) a reactivity-driven proportional controller converting the temperature demand into a reactivity signal, and (2) a temperature-driven proportional-derivative controller which relies on a derivative term to anticipate sudden changes in demand and limit temperature overshoot. A weighting function is proposed to continuously transition from one component to the other.

33 ADVANCED PROPULSION SYSTEMS↗

Efficient generalized boundary detection

Fast, efficient, and robust compression-based methods for detecting boundaries in arbitrary datasets, including sequences (1D datasets), are desired. The methods, each employing three simple algorithms, approximate the information distance between two adjacent sliding windows within a dataset. One of the algorithms calculates an initial ordered list of subsequences; while a second algorithm updates the ordered list of subsequences by dropping a first entry and appending a last entry rather than calculating completely new ordered lists with each iteration. Large values in the distance metric are indicative of boundary locations. A smoothed z-score or a wavelet-based algorithm may then be used to locate peaks in the distance metric, thereby identifying boundary locations. An adaptive version of the method employs a collection of window sizes and corresponding weighting functions, making it more amenable to real datasets with unknown, complex, and changing structures.

Ting, Christina↗

Adaptive spectra-to-exposure conversion using ridge regularized polynomial response models

Real-time gamma spectra-to-exposure conversion in aerial and ground monitoring commonly relies on calibration-derived, detector- or system-specific conversion coefficients that are assumed to generalize across operational environments. In practice, deployment specific differences in spectral composition and transport conditions can introduce systematic bias relative to reference instruments, motivating methods that adapt coefficients using minimal field supervision while explicitly limiting overfitting. In this work, we present a conservative coefficient adaptation framework that updates a baseline polynomial energy-weighting function using ridge-regularized regression, with leave-one-out cross-validation (LOOCV) used to select the regularization strength. The findings support ridge-constrained minimal-supervision adaptation as a practical mechanism to suppress site-specific bias without destabilizing a calibration-derived baseline.

61 RADIATION PROTECTION AND DOSIMETRY↗