Development of a Multimodel Ensemble Prediction System to Specify Ionospheric Weather
Reviews the development of a multimode space weather system. The system assimilates a range of ionospheric, electrodynamic and thermospheric data models.
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.
Reviews the development of a multimode space weather system. The system assimilates a range of ionospheric, electrodynamic and thermospheric data models.
A three-dimensional (3-D) Global Assimilative Ionospheric Model (GAIM) is currently being developed by a joint University of Southern California and Jet Propulsion Laboratory (JPL) team. To estimate the electron density on a global grid, GAIM uses a first-principles ionospheric physics model and the Kalman filter as one of its possible estimation techniques.
Increasingly accurate and detailed global 3-D specification of the Earth's space plasma environment is required to further understand its intricate organization and behavior. For a long time space physics and aeronomy research has been data starved due to the great variety of natural time scales involved in the plasma phenomenology. We have started developing a new approach to the global ionospheric specification called Real-Time Assimilative Mapping (RTAM). The IRI-RTAM will use data from the Global Ionospheric Radio Observatory (GIRO) to smoothly transform International Reference Ionosphere's (IRI's) background empirical maps of the ionospheric characteristics to match the observations. Such empirical assimilative modeling will provide a high-resolution, global picture of the ionospheric response to various short-term events observed during periods of storm activity or the impact of gravity waves coupling the ionosphere to the lower atmosphere, including timelines of the vertical restructuring of the plasma distribution. It will also contribute to the challenging task of providing a rapid insight into the temporal and spatial space weather development using the real-time GIRO data streams. The new assimilation technique "updates" the IRI electron density distribution while preserving the overall integrity of IRI s typical ionospheric feature representations. The technique adjusts the coefficients of the spherical/diurnal expansions used by the CCIR and URSI-88 model to obtain the global sub-peak electron density distribution. The set of global corrected coefficients can be generated as frequently as every 15 min and easily disseminated using a single real-time RTAM server operated by GIRO.
High-latitude electric potential and ionospheric conductance patterns are presented and discussed for the Northern Hemisphere during the SUNDIAL period of September 23-26, 1986 using the Assimilative Mapping of Ionospheric Electrodynamics (AMIE) technique of Richmond and Kamide (1988). Data sources used for the model are satellite and ground magnetometers, electron precipitation instruments, incoherent scatter radars, and ionospheric coherent radars. The period was characterized by many substorms, and a wide variety of instantaneous patterns of electrodynamic parameters on a hemispheric scale are derived throughout this period, of which only a few examples are displayed. The entire set of electric potentials and conductances are being made available through the NCAR CEDAR Data Base for further analysis and utilization in simulation models.
The International Space Station (ISS) orbits near the F-peak of the ionosphere (approximately 400 km altitude). Generally, satellites orbiting at this altitude would have a floating potential (FP) of approximately -1 V due to the electron temperature (Te). However, the ISS has 8 large negatively grounded 160 V solar array wings (SAW) that collect a significant electron current from the ionosphere. This current drives the ISS FP much more negative during insolation and is highly dependent on the electron density (Ne). Also, due to the size of the ISS, magnetic inductance caused by the geomagnetic field produces a delta potential up to 40 V across the truss, possibly producing positive potentials. During Extravehicular Activity (EVA) the negative FP can lead to an arcing hazard when it exceeds -45.5 V, and the positive FP can produce a DC current high enough to stimulate the astronaut's muscles and also cause a hazard. Data collected from the Floating Potential Monitoring Unit (FPMU) have shown that the probability of either of these hazards occurring during times with quiet to moderately disturbed geomagnetic activity is low enough to no longer be considered a risk. However, a study of the ionosphere Ne during severe geomagnetic storm activity has shown that the Ne can be enhanced by a factor of 6 in the ISS orbit. As a result, the ISS Safety Review Panel (SRP) requires that ionospheric conditions be monitored using the FPMU in conjunction with the ISS Plasma Interaction Model (PIM) to determine if a severe geomagnetic storm could result in a plasma environment that could produce a hazard. A 'Real-Time' plasma hazard assessment process was developed to support ISS Program real-time decision making providing constraint relief information for EVAs planning and operations. This process incorporates 'real time' ionospheric conditions, ISS solar arrays' orientation, ISS flight attitude, and where the EVA will be performed on the ISS. This assessment requires real time data that is presently provided by the FPMU including ISS floating potential, along with ionospheric Ne and Te, in order to determine the present environment. Once the present environment conditions are known to be either above, below, or near the current IRI values, the IRI is used to forecast what the environment could become in the event of a severe geomagnetic storm. If the FPMU should fail, the Space Environments team needs another source of data which is utilized to support a short-term forecast for EVAs. The IRI Real-Time Assimilative Mapping (IRTAM) model is an ionospheric model that uses real time measurements from approximately 70 digisondes to produce foF2 and hmF2 global maps in 15 minute cadence. The Boeing Space Environments team has used the IRI coefficients produced in IRTAM to calculate the Ne along the ISS orbital track. The results of the IRTAM model have been compared to FPMU measurements and show excellent agreement (figure 1). IRTAM has been identified as a potential FPMU back-up system will be used as a backup for the FPMU to support the ISS Program following completion of an FPMU/IRTAM validation campaign.
We introduce a new framework called Machine Learning (ML) based Auroral Ionospheric electrodynamics Model (ML-AIM). ML-AIM solves a current continuity equation by utilizing the ML model of Field Aligned Currents of Kunduri et al. (2020, https://doi.org/10.1029/2020JA027908), the FAC-derived auroral conductance model of Robinson et al. (2020, https://doi.org/10.1029/2020JA028008), and the solar irradiance conductance model of Moen and Brekke (1993, https://doi.org/10.1029/92gl02109). The ML-AIM inputs are 60-min time histories of solar wind plasma, interplanetary magnetic fields (IMF), and geomagnetic indices, and its outputs are ionospheric electric potential, electric fields, Pedersen/Hall currents, and Joule Heating. We conduct two ML-AIM simulations for a weak geomagnetic activity interval on 14 May 2013 and a geomagnetic storm on 7–8 September 2017. ML-AIM produces physically accurate ionospheric potential patterns such as the two-cell convection pattern and the enhancement of electric potentials during active times. The cross polar cap potentials (ΦPC) from ML-AIM, the Weimer (2005, https://doi.org/10.1029/2004ja010884) model, and the Super Dual Auroral Radar Network (SuperDARN) data-assimilated potentials, are compared to the ones from 3204 polar crossings of the Defense Meteorological Satellite Program F17 satellite, showing better performance of ML-AIM than others. ML-AIM is unique and innovative because it predicts ionospheric responses to the time-varying solar wind and geomagnetic conditions, while the other traditional empirical models like Weimer (2005, https://doi.org/10.1029/2004ja010884) designed to provide a quasi-static ionospheric condition under quasi-steady solar wind/IMF conditions. Plans are underway to improve ML-AIM performance by including a fully ML network of models of aurora precipitation and ionospheric conductance, targeting its characterization of geomagnetically active times.
The three-cornered hat (3CH) method estimates the uncertainties of three different co-located model or observational data sets (Anthes and Rieckh, 2018; Sjoberg et al., 2021). Rieckh et al. (2021) used the 3CH method to compare the random error statistics of different global forecast and reanalysis models, as well as radio occultation (RO) and radiosonde observations. That study showed that the MERRA-2 reanalysis, while having smaller errors in the stratosphere than its predecessor MERRA, had larger errors in the troposphere than many of the other data sets analyzed. The MERRA-2 errors were particularly large in the tropics. In a collaborative effort between UCAR’s COSMIC (Constellation Observing System for Meteorology, Ionosphere and Meteorology) program and NASA’s Global Modeling and Assimilation Office (GMAO), we carried out further 3CH error diagnostics to help isolate the causes of these larger errors and help guide the development of an improved reanalysis. This presentation summarizes random error statistics associated with MERRA-2, ECMWF’s ERA5 reanalysis, and COSMIC-2 (C2) RO observations. We compute 3CH error variance estimates of refractivity, as well as temperature and specific humidity using UCAR’s COSMIC Data Analysis and Archive Center (CDAAC) improved 1D-variational (1D-Var) retrieval (wetPf2) over 15 latitude bands from 45S to 45N. The 1D-Var retrievals of specific humidity and temperature for C2 use NCEP’s Global Forecast System (GFS) as the background. Anthes et al. (2021) showed that it gives accurate estimates of temperature and specific humidity in the tropics and subtropics, even in the challenging environment of intense Hurricane Dorian (2019). This presentation confirms the previous results that MERRA-2 has significantly larger errors in the tropics and subtropics than either C2 or ERA5. Its errors are larger between 30S and 30N compared to 30-45 N-S latitudes, and are also larger over land compared to oceans. Most of the MERRA-2 refractivity errors come from specific humidity, except over land below 3 km where temperature errors are large. These results suggest that moist convection and atmospheric boundary layer physics in MERRA-2 may be responsible for a significant part of the higher uncertainties. These results are being used to guide GMAO in developing an improved next-generation reanalysis, as shown in a companion presentation submitted to this conference (El Akkraoui et al., 2021), which extends this study and describes improvements to MERRA-2 leading to the next GMAO reanalysis.
Explore the source record for details and available documents.
No abstract available
The Cyclone Global Navigation Satellite System (CYGNSS) mission was designed to study hurricane intensification by measuring wind speeds in tropical cyclones. However, the delay-Doppler maps (DDM) that are produced can be used to estimate the sea surface height (SSH) at the specular reflection point on the ocean surface. Proofof-concept studies that DDMs are suitable to solve for SSH have been recently reported (Clarizia et al., 2016; Mashburn et al., 2018), based on data acquired by the demonstration satellite experiment Tech Demo Sat - 1 (TDS-1) carrying a GNSS-R receiver similar to the ones onboard CYGNSS. Although the precision of each 1sec averaged SSH is considerably lower than that of the existing satellite altimeters, by virtue of the dense coverage and frequent revisit time exhibited by the constellation of 8 microsats, the error may be smoothed down considerably by optimal interpolation (Li et al., 2016). Hence the CYGNSS dataset presents a potential opportunity to sample the tropical oceans, and investigate the sensitivity of the SSH measurements to mesoscale eddies. Our objective is to analyze the CYGNSS data and test a suite of retrieval algorithms, including that of Mashburn et al., 2018 to obtain SSH, to understand the error sources and possible corrections when available. Because of the limited antenna gain of the CYGNSS instrument, the power SNR for general wind conditions is expected to be too low to accurately track the reflection point contribution, resulting in altimetry errors of the order of several meters. Other significant error sources are the CYGNSS satellites orbit knowledge, the ionosphere, the troposphere, the mean sea surface and the tides; some of these errors can be reduced by using models. The corrected SSH data will be assimilated into the ROMS high-resolution ocean model to investigate their ability to reproduce mesoscale eddies in the tropical oceans.
Topics covered include: Measurements of Ultra-Stable Oscillator (USO) Allan Deviations in Space; Gaseous Nitrogen Orifice Mass Flow Calculator; Validation of Proposed Metrics for Two-Body Abrasion Scratch Test Analysis Standards; Rover Low Gain Antenna Qualification for Deep Space Thermal Environments; Automated, Ultra-Sterile Solid Sample Handling and Analysis on a Chip; Measuring and Estimating Normalized Contrast in Infrared Flash Thermography; Spectrally and Radiometrically Stable, Wideband, Onboard Calibration Source; High-Reliability Waveguide Vacuum/Pressure Window; Methods of Fabricating Scintillators With Radioisotopes for Beta Battery Applications; Magnetic Shield for Adiabatic Demagnetization Refrigerators (ADR); CMOS-Compatible SOI MESFETS for Radiation-Hardened DC-to-DC Converters; Silicon Heat Pipe Array; Adaptive Phase Delay Generator; High-Temperature, Lightweight, Self-Healing Ceramic Composites for Aircraft Engine Applications; Treatment to Control Adhesion of Silicone-Based Elastomers; High-Temperature Adhesives for Thermally Stable Aero-Assist Technologies; Rockballer Sample Acquisition Tool; Rock Gripper for Sampling, Mobility, Anchoring, and Manipulation; Advanced Magnetic Materials Methods and Numerical Models for Fluidization in Microgravity and Hypogravity; Data Transfer for Multiple Sensor Networks Over a Broad Temperature Range; Using Combustion Synthesis to Reinforce Berms and Other Regolith Structures; Visible-Infrared Hyperspectral Image Projector; Three-Axis Attitude Estimation With a High-Bandwidth Angular Rate Sensor Change_Detection.m; AGATE: Adversarial Game Analysis for Tactical Evaluation; Ionospheric Simulation System for Satellite Observations and Global Assimilative; Modeling Experiments (ISOGAME); An Extensible, User- Modifiable Framework for Planning Activities; Mission Operations Center (MOC) - Precipitation Processing System (PPS) Interface Software System (MPISS); Automated 3D Damaged Cavity Model Builder for Lower Surface Acreage Tile on Orbiter; Mixed Linear/Square-Root Encoded Single-Slope Ramp Provides Low-Noise ADC with High Linearity for Focal Plane Arrays; RUSHMAPS: Real-Time Uploadable Spherical Harmonic Moment Analysis for Particle Spectrometers; Powered Descent Guidance with General Thrust-Pointing Constraints; X-Ray Detection and Processing Models for Spacecraft Navigation and Timing; and Extreme Ionizing-Radiation-Resistant Bacterium
We present data from the Floating Potential Measurement Unit (FPMU), that is deployed on the starboard (S1) truss of the International Space Station. The FPMU is a suite of instruments capable of redundant measurements of various plasma parameters. The instrument suite consists of: a Floating Potential Probe, a Wide-sweeping spherical Langmuir probe, a Narrow-sweeping cylindrical Langmuir Probe, and a Plasma Impedance Probe. This paper gives a brief overview of the instrumentation and the received data quality, and then presents the algorithm used to reduce I-V curves to plasma parameters. Several hours of data is presented from August 5th, 2006 and March 3rd, 2007. The FPMU derived plasma density and temperatures are compared with the International Reference Ionosphere (IRI) and USU-Global Assimilation of Ionospheric Measurement (USU-GAIM) models. Our results show that the derived in-situ density matches the USU-GAIM model better than the IRI, and the derived in-situ temperatures are comparable to the average temperatures given by the IRI.
A principal obstacle to fielding automated planning systems is the difficulty of modeling. Physical systems are modeled conventionally based on specification documents and the modeler's understanding of the system. Thus, the model is developed in a way that is disconnected from the system's actual behavior and is vulnerable to manual error. Another obstacle to fielding planners is testing and validation. For a space mission, generated plans must be validated often by translating them into command sequences that are run in a simulation testbed. Testing in this way is complex and onerous because of the large number of possible plans and states of the spacecraft. Though, if used as a source of domain knowledge, the simulator can ease validation. This paper poses a challenge: to ground planning models in the system physics represented by simulation. A proposed, interactive model development environment illustrates the integration of planning and simulation to meet the challenge. This integration reveals research paths for automated model construction and validation.
We analyze retrospectively/prospectively the transient variations of six different physical parameters in the atmosphere/ionosphere during the M7.8 and M7.3 earthquakes in Nepal, namely: 1) outgoing longwave radiation (OLR) at the top of the atmosphere (TOA); 2) GPS/TEC; 3) the very-low-frequency (VLF/LF) signals at the receiving stations in Bishkek (Kyrgyzstan) and Varanasi (India); 4) Radon observations; 5) Atmospheric chemical potential from assimilation models; and; 6) Air Temperature from NOAA ground stations. We found that in mid-March 2015, there was a rapid increase in the radiation from the atmosphere observed by satellites. This anomaly was located close to the future M7.8 epicenter and reached a maximum on April 21–22. The GPS/TEC data analysis indicated an increase and variation in electron density, reaching a maximum value during April 22–24. A strong negative TEC anomaly in the crest of EIA (Equatorial Ionospheric Anomaly) occurred on April 21, and a strong positive anomaly was recorded on April 24, 2015. The behavior of VLF-LF waves along NWC-Bishkek and JJY-Varanasi paths has shown abnormal behavior during April 21–23, several days before the first, stronger earthquake. Our continuous satellite OLR analysis revealed this new strong anomaly on May 3, which was why we anticipated another major event in the area. On May 12, 2015, an M7.3 earthquake occurred. Our results show coherence between the appearance of these pre-earthquake transient’s effects in the atmosphere and ionosphere (with a short time-lag, from hours up to a few days) and the occurrence of the 2015 M7.8 and M7.3 events. The spatial characteristics of the pre-earthquake anomalies were associated with a large area but inside the preparation region estimated by Dobrovolsky-Bowman. The preearthquake nature of the signals in the atmosphere and ionosphere was revealed by simultaneous analysis of satellite, GPS/TEC, and VLF/LF and suggest that they follow a general temporal-spatial evolution pattern that has been seen in other large earthquakes worldwide
This study is a preliminary analysis of the effectiveness (in terms of altimeter calibration accuracy) of various ionosphere models and the Global Positioning System (GPS) to calibrate single frequency altimeter height measurements for ionospheric path delay. In particular, the research focused on ingesting GPS Total Electron Content (TEC) data into the physical Parameterized Real-Time Ionospheric Specification Model (PRISM), which estimates the composition of the ionosphere using independent empirical and physical models and has the capability of adjusting to additional ionospheric measurements. Two types of GPS data were used to adjust the PRISM model: GPS receiver station data mapped from line-of-sight observations to the vertical at the point of interest and a grid map (generated at the Jet Propulsion Laboratory) of GPS derived TEC in a sun-fixed longitude frame. The adjusted PRISM TEC values, as well as predictions by the International Reference Ionosphere (IRI-90), a climatological (monthly mean) model of the ionosphere, were compared to TOPEX dual-frequency TEC measurements (considered as truth) for a number of TOPEX sub-satellite tracks. For a 13.6 GHz altimeter, a Total Electron Content (TEC) of 1 TECU 10(exp 16) electrons/sq m corresponds to approximately 0.218 centimeters of range delay. A maximum expected TEC (at solar maximum or during solar storms) of 10(exp 18) electrons/sq m will create 22 centimeters of range delay. Compared with the TOPEX data, the PRISM predictions were generally accurate within the TECU when the sub-satellite track of interest passed within 300 to 400 km of the GPS TEC data or when the track passed through a night-time ionosphere. If neither was the case, in particular if the track passed through a local noon ionosphere, the PRISM values differed by more than 10 TECU and by as much as 40 TECU. The IRI-90 model, with no current ability to unseat GPS data, predicted TEC to a slightly higher error of 12 TECU. The performance of PRISM is very promising for predicting TEC and will prove useful for calibrating single frequency altimeter height measurements for ionospheric path delay. When adjusted to the GPS line-of-sight data the PRISM URSI empirical model predicted TEC over a day's period to within a global error of 8.60 TECU rms during a nighttime ionosphere and 9.74 TECU rms during the day. When adjusted to the GPS derived TEC grid, the PRISM parametrized model predicted TEC to within an error of 8.47 TECU rms for a nighttime ionosphere and 12.83 TECU rms during the day. However, the grid cannot be considered globally due to the lack of sufficient numbers of GPS stations and large latitude gaps in GPS data. It is the opinion of the authors that using the PRISM model and adjusting to the global sun-fixed TEC grid regenerated with a localized weighted interpolation offers the best possibility of meeting the 10 TECU global rms (or 2 cm at 13.6 GHz) ionosphere range correction accuracy requirement of TOPEX/Poseidon and should be the subject of further study. However, it is clear that the anticipated requirement of 34 TECU global rms for TOPEX/Poseidon Follow-On (corresponding to the TOPEX/Poseidon performance) can not be met with any realizable combination of existing models and data assimilation schemes.
This study is a preliminary analysis of the accuracy of various ionosphere models to correct single frequency altimeter height measurements for Ionospheric path delay. In particular, research focused on adjusting empirical and parameterized ionosphere models in the parameterized real-time ionospheric specification model (PRISM) 1.2 using total electron content (TEC) data from the global positioning system (GPS). The types of GPS data used to adjust PRISM included GPS line-of-sight (LOS) TEC data mapped to the vertical, and a grid of GPS derived TEC data in a sun-fixed longitude frame. The adjusted PRISM TEC values, as well as predictions by IRI-90, a climatotogical model, were compared to TOPEX/Poseidon (T/P) TEC measurements from the dual-frequency altimeter for a number of T/P tracks. When adjusted with GPS LOS data, the PRISM empirical model predicted TEC over 24 1 h data sets for a given local time to with in a global error of 8.60 TECU rms during a midnight centered ionosphere and 9.74 TECU rms during a noon centered ionosphere. Using GPS derived sun-fixed TEC data, the PRISM parameterized model predicted TEC within an error of 8.47 TECU rms centered at midnight and 12.83 TECU rms centered at noon. From these best results, it is clear that the proposed requirement of 3-4 TECU global rms for TOPEX/Poseidon Follow-On will be very difficult to meet, even with a substantial increase in the number of GPS ground stations, with any realizable combination of the aforementioned models or data assimilation schemes.
Detailed meteorological analyses based on observations extending through the middle atmosphere (~15 to 100 km altitude) can provide key information to whole atmosphere modelling systems regarding the physical mechanisms linking day-to-day changes in ionospheric electron density to meteorological variability near the Earth’s surface. However, the extent to which independent middle atmosphere analyses differ in their representation of wave-induced coupling to the ionosphere is unclear. To begin to address this issue, we present the first intercomparison among four such analyses, JAGUAR-DAS, MERRA-2, NAVGEM-HA, and WACCMX+DART, focusing on the Northern Hemisphere (NH) 2009-2010 winter, which includes a major sudden stratospheric warming (SSW). This intercomparison examines the altitude, latitude, and time dependences of zonal mean zonal winds and temperatures among these four analyses over the 1 December 2009 – 31 March 2010 period, as well as latitude and altitude dependences of monthly mean amplitudes of the diurnal and semidiurnal migrating solar tides, the eastward propagating diurnal zonal wave number 3 nonmigrating tide, and traveling planetary waves associated with the quasi-5-day and quasi-2-day Rossby modes. Our results show generally good agreement among the four analyses up to the stratopause (~50 km altitude). Large discrepancies begin to emerge in the mesosphere and lower thermosphere owing to (1) differences in the types of satellite data assimilated by each system and (2) differences in the details of the global atmospheric models used by each analysis system. The results of this intercomparison provide initial estimates of uncertainty in analyses commonly used to constrain middle atmospheric meteorological variability in whole atmosphere model simulations.
Space weather nowcasting and forecasting models require assimilation of near‐real time (NRT) space environment data to improve the precision and accuracy of operational products. Typically, these models begin with a climatological model to provide "most probable distributions" of environmental parameters as a function of time and space. The process of NRT data assimilation gently pulls the climate model closer toward the observed state (e.g. via Kalman smoothing) for nowcasting, and forecasting is achieved through a set of iterative physics‐based forward‐prediction calculations. The issue of required space weather observatories to meet the spatial and temporal requirements of these models is a complex one, and we do not address that with this poster. Instead, we present some examples of how tethered satellites can be used to address the shortfalls in our ability to measure critical environmental parameters necessary to drive these space weather models. Examples include very long baseline electric field measurements, magnetized ionospheric conductivity measurements, and the ability to separate temporal from spatial irregularities in environmental parameters. Tethered satellite functional requirements will be presented for each space weather parameter considered in this study.