Quantifying the Storm Time Thermospheric Neutral Density Variations Using Model and Observations

Accurate determination of thermospheric neutral density holds crucial importance for satellite drag calculations. The problem is twofold and involves the correct estimation of the quiet time climatology and storm time variations. In this work, neutral density estimations from two empirical and three physicsâ based models of the ionosphereâ thermosphere are compared with the neutral densities along the Challenging Microâ Satellite Payload satellite track for six geomagnetic storms. Storm time variations are extracted from neutral density by (1) subtracting the mean difference between model and observation (bias), (2) setting climatological variations to zero, and (3) multiplying model data with the quiet time ratio between the model and observation. Several metrics are employed to evaluate the model performances. We find that the removal of bias or climatology reveals actual performance of the model in simulating the storm time variations. When bias is removed, depending on event and model, storm time errors in neutral density can decrease by an amount of 113% or can increase by an amount of 12% with respect to error in models with quiet time bias. It is shown that using only average and maximum values of neutral density to determine the model performances can be misleading since a model can estimate the averages fairly well but may not capture the maximum value or vice versa. Since each of the metrics used for determining model performances provides different aspects of the error, among these, we suggest employing mean absolute error, prediction efficiency, and normalized root mean square error together as a standard set of metrics for the neutral density.Plain Language SummaryThermospheric neutral density is the largest source of uncertainty in atmospheric drag calculations. Consequently, mission and maneuver planning, satellite lifetime predictions, collision avoidance, and orbit determination depend on the accurate estimation of the thermospheric neutral density. Thermospheric neutral density varies in different timescales. In short timescales, the largest variations occur due to the geomagnetic storms. Several empirical and physicsâ based models of the ionosphereâ thermosphere system are used for estimating the variations in the neutral density. However, the storm time responses from the models are clouded by the climatology (background variations), upon which the effect of geomagnetic storms is superimposed. In this work, we show that it is critical to use reference levels for the neutral density to extract the true performance of the models for the evaluation of the storm time performances. We demonstrate that mean absolute error, prediction efficiency, and normalized root mean square error should be considered together for the performance evaluations, since each of them provides different aspects of the error.Key PointsUsing the average and maximum values of neutral densities to determine the model performances can be misleadingRemoving the quiet time trend from the neutral density reveals the actual performance of the model in simulating the storm time variationsMean absolute error, prediction efficiency, and normalized root mean square error should be considered together for the evaluationsPeer Reviewedhttps://deepblue.lib.umich.edu/bitstream/2027.42/148396/1/swe20816_am.pdfhttps://deepblue.lib.umich.edu/bitstream/2027.42/148396/2/swe20816-sup-0001-2018SW002033-SI.pdfhttps://deepblue.lib.umich.edu/bitstream/2027.42/148396/3/swe20816.pd

icebearcanada/cavsiopy: v1.1.1

patch release: 'attitude_3d_ground_quiver' has been enhanced to display a line connecting the subsatellite point with the ground target on the ground map. name changes for several functions in auxiliary directory

icebearcanada/cavsiopy: v1.2.2

slew_example.py: added to examples in attitude_analysis.py: functions in utils.py and complement_missing_sofa.py are now embedded in attitude_analysis.py use_rotation_matrices.py: utils.py and missing_complement_sofa.py imports removed attitude_analysis.spacecraft_distance_from_a_point: fixed a minor bug, which caused the distance array to return empty. requirements.txt: version numbers were added for dependencies. .readthedocs.yaml: version number and installation instructions corrected