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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
NASA NTRS2014

Transition Marshall Space Flight Center Wind Profiler Splicing Algorithm to Launch Services Program Upper Winds Tool

NASAs LSP customers and the future SLS program rely on observations of upper-level winds for steering, loads, and trajectory calculations for the launch vehicles flight. On the day of launch, the 45th Weather Squadron (45 WS) Launch Weather Officers (LWOs) monitor the upper-level winds and provide forecasts to the launch team via the AMU-developed LSP Upper Winds tool for launches at Kennedy Space Center (KSC) and Cape Canaveral Air Force Station. This tool displays wind speed and direction profiles from rawinsondes released during launch operations, the 45th Space Wing 915-MHz Doppler Radar Wind Profilers (DRWPs) and KSC 50-MHz DRWP, and output from numerical weather prediction models.The goal of this task was to splice the wind speed and direction profiles from the 45th Space Wing (45 SW) 915-MHz Doppler radar Wind Profilers (DRWPs) and KSC 50-MHz DRWP at altitudes where the wind profiles overlap to create a smooth profile. In the first version of the LSP Upper Winds tool, the top of the 915-MHz DRWP wind profile and the bottom of the 50-MHz DRWP were not spliced, sometimes creating a discontinuity in the profile. The Marshall Space Flight Center (MSFC) Natural Environments Branch (NE) created algorithms to splice the wind profiles from the two sensors to generate an archive of vertically complete wind profiles for the SLS program. The AMU worked with MSFC NE personnel to implement these algorithms in the LSP Upper Winds tool to provide a continuous spliced wind profile.The AMU transitioned the MSFC NE algorithms to interpolate and fill data gaps in the data, implement a Gaussian weighting function to produce 50-m altitude intervals in each sensor, and splice the data together from both DRWPs. They did so by porting the MSFC NE code written with MATLAB software into Microsoft Excel Visual Basic for Applications (VBA). After testing the new algorithms in stand-alone VBA modules, the AMU replaced the existing VBA code in the LSP Upper Winds tool with the new algorithms. They then tested the code in the LSP Upper Winds tool with archived data. The tool will be delivered to the 45 WS after the 50-MHz DRWP upgrade is complete and the tool is tested with real-time data. The 50-MHz DRWP upgrade is expected to be finished in October 2014.

Space launch↗
NASA NTRS2014

Technology Development for a Hyperspectral Microwave Atmospheric Sounder (HyMAS)

The Hyperspectral Microwave Atmospheric Sounder (HyMAS) is being developed at Lincoln Laboratories and accommodated by the Goddard Space Flight Center for a flight opportunity on a NASA research aircraft. The term hyperspectral microwave is used to indicate an all-weather sounding instrument that performs equivalent to hyperspectral infrared sounders in clear air with vertical resolution of approximately 1 km. Deploying the HyMAS equipped scanhead with the existing Conical Scanning Microwave Imaging Radiometer (CoSMIR) shortens the path to a flight demonstration. Hyperspectral microwave is achieved through the use of independent RF antennas that sample the volume of the Earths atmosphere through various levels of frequencies, thereby producing a set of dense, spaced vertical weighting functions.

sounding↗
NASA NTRS2015

Hyperspectral Microwave Atmospheric Sounder (HyMAS) - New Capability in the CoSMIR-CoSSIR Scanhead

Lincoln Laboratory and NASA's Goddard Space Flight Center have teamed to re-use an existing instrument platform, the CoSMIR/CoSSIR system for atmospheric sounding, to develop a new capability in hyperspectral filtering, data collection, and display. The volume of the scanhead accomodated an intermediate frequency processor(IFP), that provides the filtering and digitization of the raw data and the interoperable remote component (IRC) adapted to CoSMIR, CoSSIR, and HyMAS that stores and archives the data with time tagged calibration and navigation data. The first element of the work is the demonstration of a hyperspectral microwave receiver subsystem that was recently shown using a comprehensive simulation study to yield performance that substantially exceeds current state-of-the-art. Hyperspectral microwave sounders with approximately 100 channels offer temperature and humidity sounding improvements similar to those obtained when infrared sensors became hyperspectral, but with the relative insensitivity to clouds that characterizes microwave sensors. Hyperspectral microwave operation is achieved using independent RF antenna/receiver arrays that sample the same area/volume of the Earth's surface/atmosphere at slightly different frequencies and therefore synthesize a set of dense, finely spaced vertical weighting functions. The second, enabling element of the proposal is the development of a compact 52-channel Intermediate Frequency processor module. A principal challenge in the development of a hyperspectral microwave system is the size of the IF filter bank required for channelization. Large bandwidths are simultaneously processed, thus complicating the use of digital back-ends with associated high complexities, costs, and power requirements. Our approach involves passive filters implemented using low-temperature co-fired ceramic (LTCC) technology to achieve an ultra-compact module that can be easily integrated with existing radio frequency front-end technology. This IF processor is universally applicable to other microwave sensing missions requiring compact IF spectrometry. The data include 52 operational channels with low IF module volume (less than 100 cubic centimeters) and mass (less than 300 grams) and linearity better than 0.3 percent over a 330,000 dynamic range.

multispectral↗
NASA NTRS2015

Hyperspectral Microwave Atmospheric Sounder (HyMas) - New Capability in the CoSMIR-CoSSIR Scanhead

Lincoln Laboratory and NASA's Goddard Space Flight Center have teamed to re-use an existing instrument platform, the CoSMIRCoSSIR system for atmospheric sounding, to develop a new capability in hyperspectral filtering, data collection, and display. The volume of the scanhead accomodated an intermediate frequency processor(IFP), that provides the filtering and digitization of the raw data and the interoperable remote component (IRC) adapted to CoSMIR, CoSSIR, and HyMAS that stores and archives the data with time tagged calibration and navigation data.The first element of the work is the demonstration of a hyperspectral microwave receiver subsystem that was recently shown using a comprehensive simulation study to yield performance that substantially exceeds current state-of-the-art. Hyperspectral microwave sounders with 100 channels offer temperature and humidity sounding improvements similar to those obtained when infrared sensors became hyperspectral, but with the relative insensitivity to clouds that characterizes microwave sensors. Hyperspectral microwave operation is achieved using independent RF antennareceiver arrays that sample the same areavolume of the Earths surfaceatmosphere at slightly different frequencies and therefore synthesize a set of dense, finely spaced vertical weighting functions. The second, enabling element of the proposal is the development of a compact 52-channel Intermediate Frequency processor module. A principal challenge in the development of a hyperspectral microwave system is the size of the IF filter bank required for channelization. Large bandwidths are simultaneously processed, thus complicating the use of digital back-ends with associated high complexities, costs, and power requirements. Our approach involves passive filters implemented using low-temperature co-fired ceramic (LTCC) technology to achieve an ultra-compact module that can be easily integrated with existing RF front-end technology. This IF processor is universally applicable to other microwave sensing missions requiring compact IF spectrometry.The data include 52 operational channels with low IF module volume (100cm3) and mass (300g) and linearity better than 0.3 over a 330K dynamic range.

sounding↗
NASA NTRS2012

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↗
NASA NTRS2014

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↗
NASA NTRS2015

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↗
NASA NTRS2015

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.↗
NASA NTRS2017

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.↗
NASA NTRS2018

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↗
NASA NTRS2018

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↗
NASA NTRS2020

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↗
NASA NTRS2022

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↗
NASA NTRS2022

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↗
NASA NTRSDate not supplied

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↗
NASA NTRS2023

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↗
NASA NTRS2023

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↗