Search

Daedalus MASE (Mission Assessment through Simulation Exercise) is an open-source package of scientific analysis tools aimed at research in the Lower Thermosphere-Ionosphere (LTI). It was created with the purpose to assess the performance and demonstrate closure of the mission objectives of Daedalus, a mission concept targeting to perform in-situ measurements in the LTI. However, through its successful usage as a mission-simulator toolset, Daedalus MASE has evolved to encompass numerous capabilities related to LTI science and modeling. Inputs are geophysical observables in the LTI, which can be obtained either through in-situ measurements from spacecraft and rockets, or through Global Circulation Models (GCM). These include ion, neutral and electron densities, ion and neutral composition, ion, electron and neutral temperatures, ion drifts, neutral winds, electric field, and magnetic field. In the examples presented, these geophysical observables are obtained through NCAR's Thermosphere-Ionosphere-Electrodynamics General Circulation Model. Capabilities of Daedalus MASE include: 1) Calculations of products that are derived from the above geophysical observables, such as Joule heating, energy transfer rates between species, electrical currents, electrical conductivity, ion-neutral collision frequencies between all combinations of species, as well as height-integrations of derived products. 2) Calculation and cross-comparison of collision frequencies and estimates of the effect of using different models of collision frequencies into derived products. 3) Calculation of the uncertainties of derived products based on the uncertainties of the geophysical observables, due to instrument errors or to uncertainties in measurement techniques. 4) Routines for the along-orbit interpolation within gridded datasets of GCMs. 5) Routines for the calculation of the global coverage of an in situ mission in regions of interest and for various conditions of solar and geomagnetic activity.

Daedalus MASE (Mission Assessment through Simulation Exercise) is an open-source package of scientific analysis tools aimed at research in the Lower Thermosphere-Ionosphere (LTI). It was created with the purpose to assess the performance and demonstrate closure of the mission objectives of Daedalus, a mission concept targeting to perform in-situ measurements in the LTI. However, through its successful usage as a mission-simulator toolset, Daedalus MASE has evolved to encompass numerous capabilities related to LTI science and modeling. Inputs are geophysical observables in the LTI, which can be obtained either through in-situ measurements from spacecraft and rockets, or through Global Circulation Models (GCM). These include ion, neutral and electron densities, ion and neutral composition, ion, electron and neutral temperatures, ion drifts, neutral winds, electric field, and magnetic field. In the examples presented, these geophysical observables are obtained through NCAR's Thermosphere-Ionosphere-Electrodynamics General Circulation Model. Capabilities of Daedalus MASE include: 1) Calculations of products that are derived from the above geophysical observables, such as Joule heating, energy transfer rates between species, electrical currents, electrical conductivity, ion-neutral collision frequencies between all combinations of species, as well as height-integrations of derived products. 2) Calculation and cross-comparison of collision frequencies and estimates of the effect of using different models of collision frequencies into derived products. 3) Calculation of the uncertainties of derived products based on the uncertainties of the geophysical observables, due to instrument errors or to uncertainties in measurement techniques. 4) Routines for the along-orbit interpolation within gridded datasets of GCMs. 5) Routines for the calculation of the global coverage of an in situ mission in regions of interest and for various conditions of solar and geomagnetic activity.

Daedalus MASE (Mission Assessment through Simulation Exercise) is an open-source package of scientific analysis tools aimed at research in the Lower Thermosphere-Ionosphere (LTI). It was created with the purpose to assess the performance and demonstrate closure of the mission objectives of Daedalus, a mission concept targeting to perform in-situ measurements in the LTI. However, through its successful usage as a mission-simulator toolset, Daedalus MASE has evolved to encompass numerous capabilities related to LTI science and modeling. Inputs are geophysical observables in the LTI, which can be obtained either through in-situ measurements from spacecraft and rockets, or through Global Circulation Models (GCM). These include ion, neutral and electron densities, ion and neutral composition, ion, electron and neutral temperatures, ion drifts, neutral winds, electric field, and magnetic field. In the examples presented, these geophysical observables are obtained through NCAR's Thermosphere-Ionosphere-Electrodynamics General Circulation Model. Capabilities of Daedalus MASE include: 1) Calculations of products that are derived from the above geophysical observables, such as Joule heating, energy transfer rates between species, electrical currents, electrical conductivity, ion-neutral collision frequencies between all combinations of species, as well as height-integrations of derived products. 2) Calculation and cross-comparison of collision frequencies and estimates of the effect of using different models of collision frequencies into derived products. 3) Calculation of the uncertainties of derived products based on the uncertainties of the geophysical observables, due to instrument errors or to uncertainties in measurement techniques. 4) Routines for the along-orbit interpolation within gridded datasets of GCMs. 5) Routines for the calculation of the global coverage of an in situ mission in regions of interest and for various conditions of solar and geomagnetic activity.

The lower thermosphere-ionosphere (LTI) is a key transition region between Earth’s atmosphere and space. Interactions between ions and neutrals maximize within the LTI and in particular at altitudes from 100 to 200 km, which is the least visited region of the near-Earth environment. The lack of in situ co-temporal and co-spatial measurements of all relevant parameters and their elusiveness to most remote-sensing methods means that the complex interactions between its neutral and charged constituents remain poorly characterized to this date. This lack of measurements, together with the ambiguity in the quantification of key processes in the 100–200 km altitude range affect current modeling efforts to expand atmospheric models upward to include the LTI and limit current space weather prediction capabilities. We present focused questions in the LTI that are related to the complex interactions between its neutral and charged constituents. These questions concern core physical processes that govern the energetics, dynamics, and chemistry of the LTI and need to be addressed as fundamental and long-standing questions in this critically unexplored boundary region. We also outline the range of in situ measurements that are needed to unambiguously quantify key LTI processes within this region, and present elements of an in situ concept based on past proposed mission concepts.

The lower thermosphere-ionosphere (LTI) is a key transition region between Earth’s atmosphere and space. Interactions between ions and neutrals maximize within the LTI and in particular at altitudes from 100 to 200 km, which is the least visited region of the near-Earth environment. The lack of in situ co-temporal and co-spatial measurements of all relevant parameters and their elusiveness to most remote-sensing methods means that the complex interactions between its neutral and charged constituents remain poorly characterized to this date. This lack of measurements, together with the ambiguity in the quantification of key processes in the 100–200 km altitude range affect current modeling efforts to expand atmospheric models upward to include the LTI and limit current space weather prediction capabilities. We present focused questions in the LTI that are related to the complex interactions between its neutral and charged constituents. These questions concern core physical processes that govern the energetics, dynamics, and chemistry of the LTI and need to be addressed as fundamental and long-standing questions in this critically unexplored boundary region. We also outline the range of in situ measurements that are needed to unambiguously quantify key LTI processes within this region, and present elements of an in situ concept based on past proposed mission concepts.

Daedalus MASE (Mission Assessment through Simulation Exercise) is an open-source package of scientific analysis tools aimed at research in the Lower Thermosphere-Ionosphere (LTI). It was created with the purpose to assess the performance and demonstrate closure of the mission objectives of Daedalus, a mission concept targeting to perform in-situ measurements in the LTI. However, through its successful usage as a mission-simulator toolset, Daedalus MASE has evolved to encompass numerous capabilities related to LTI science and modeling. Inputs are geophysical observables in the LTI, which can be obtained either through in-situ measurements from spacecraft and rockets, or through Global Circulation Models (GCM). These include ion, neutral and electron densities, ion and neutral composition, ion, electron and neutral temperatures, ion drifts, neutral winds, electric field, and magnetic field. In the examples presented, these geophysical observables are obtained through NCAR’s Thermosphere-Ionosphere-Electrodynamics General Circulation Model. Capabilities of Daedalus MASE include: 1) Calculations of products that are derived from the above geophysical observables, such as Joule heating, energy transfer rates between species, electrical currents, electrical conductivity, ion-neutral collision frequencies between all combinations of species, as well as height-integrations of derived products. 2) Calculation and cross-comparison of collision frequencies and estimates of the effect of using different models of collision frequencies into derived products. 3) Calculation of the uncertainties of derived products based on the uncertainties of the geophysical observables, due to instrument errors or to uncertainties in measurement techniques. 4) Routines for the along-orbit interpolation within gridded datasets of GCMs. 5) Routines for the calculation of the global coverage of an in situ mission in regions of interest and for various conditions of solar and geomagnetic activity., QC 20230627

The lower thermosphere-ionosphere (LTI) is a key transition region between Earth's atmosphere and space. Interactions between ions and neutrals maximize within the LTI and in particular at altitudes from 100 to 200 km, which is the least visited region of the near-Earth environment. The lack of in situ co-temporal and co-spatial measurements of all relevant parameters and their elusiveness to most remote-sensing methods means that the complex interactions between its neutral and charged constituents remain poorly characterized to this date. This lack of measurements, together with the ambiguity in the quantification of key processes in the 100-200 km altitude range affect current modeling efforts to expand atmospheric models upward to include the LTI and limit current space weather prediction capabilities. We present focused questions in the LTI that are related to the complex interactions between its neutral and charged constituents. These questions concern core physical processes that govern the energetics, dynamics, and chemistry of the LTI and need to be addressed as fundamental and long-standing questions in this critically unexplored boundary region. We also outline the range of in situ measurements that are needed to unambiguously quantify key LTI processes within this region, and present elements of an in situ concept based on past proposed mission concepts., QC 20230329

The lower thermosphere-ionosphere (LTI) is a key transition region between Earth’s atmosphere and space. Interactions between ions and neutrals maximize within the LTI and in particular at altitudes from 100 to 200 km, which is the least visited region of the near-Earth environment. The lack of in situ co-temporal and co-spatial measurements of all relevant parameters and their elusiveness to most remote-sensing methods means that the complex interactions between its neutral and charged constituents remain poorly characterized to this date. This lack of measurements, together with the ambiguity in the quantification of key processes in the 100–200 km altitude range affect current modeling efforts to expand atmospheric models upward to include the LTI and limit current space weather prediction capabilities. We present focused questions in the LTI that are related to the complex interactions between its neutral and charged constituents. These questions concern core physical processes that govern the energetics, dynamics, and chemistry of the LTI and need to be addressed as fundamental and long-standing questions in this critically unexplored boundary region. We also outline the range of in situ measurements that are needed to unambiguously quantify key LTI processes within this region, and present elements of an in situ concept based on past proposed mission concepts.

Daedalus MASE (Mission Assessment through Simulation Exercise) is an open-source package of scientific analysis tools aimed at research in the Lower Thermosphere-Ionosphere (LTI). It was created with the purpose to assess the performance and demonstrate closure of the mission objectives of Daedalus, a mission concept targeting to perform in-situ measurements in the LTI. However, through its successful usage as a mission-simulator toolset, Daedalus MASE has evolved to encompass numerous capabilities related to LTI science and modeling. Inputs are geophysical observables in the LTI, which can be obtained either through in-situ measurements from spacecraft and rockets, or through Global Circulation Models (GCM). These include ion, neutral and electron densities, ion and neutral composition, ion, electron and neutral temperatures, ion drifts, neutral winds, electric field, and magnetic field. In the examples presented, these geophysical observables are obtained through NCAR's Thermosphere-Ionosphere-Electrodynamics General Circulation Model. Capabilities of Daedalus MASE include: 1) Calculations of products that are derived from the above geophysical observables, such as Joule heating, energy transfer rates between species, electrical currents, electrical conductivity, ion-neutral collision frequencies between all combinations of species, as well as height-integrations of derived products. 2) Calculation and cross-comparison of collision frequencies and estimates of the effect of using different models of collision frequencies into derived products. 3) Calculation of the uncertainties of derived products based on the uncertainties of the geophysical observables, due to instrument errors or to uncertainties in measurement techniques. 4) Routines for the along-orbit interpolation within gridded datasets of GCMs. 5) Routines for the calculation of the global coverage of an in situ mission in regions of interest and for various conditions of solar and geomagnetic activity.

Daedalus MASE (Mission Assessment through Simulation Exercise) is an open-source package of scientific analysis tools aimed at research in the Lower Thermosphere-Ionosphere (LTI). It was created with the purpose to assess the performance and demonstrate closure of the mission objectives of Daedalus, a mission concept targeting to perform in-situ measurements in the LTI. However, through its successful usage as a mission-simulator toolset, Daedalus MASE has evolved to encompass numerous capabilities related to LTI science and modeling. Inputs are geophysical observables in the LTI, which can be obtained either through in-situ measurements from spacecraft and rockets, or through Global Circulation Models (GCM). These include ion, neutral and electron densities, ion and neutral composition, ion, electron and neutral temperatures, ion drifts, neutral winds, electric field, and magnetic field. In the examples presented, these geophysical observables are obtained through NCAR's Thermosphere-Ionosphere-Electrodynamics General Circulation Model. Capabilities of Daedalus MASE include: 1) Calculations of products that are derived from the above geophysical observables, such as Joule heating, energy transfer rates between species, electrical currents, electrical conductivity, ion-neutral collision frequencies between all combinations of species, as well as height-integrations of derived products. 2) Calculation and cross-comparison of collision frequencies and estimates of the effect of using different models of collision frequencies into derived products. 3) Calculation of the uncertainties of derived products based on the uncertainties of the geophysical observables, due to instrument errors or to uncertainties in measurement techniques. 4) Routines for the along-orbit interpolation within gridded datasets of GCMs. 5) Routines for the calculation of the global coverage of an in situ mission in regions of interest and for various conditions of solar and geomagnetic activity.

The lower thermosphere-ionosphere (LTI) is a key transition region between Earth’s atmosphere and space. Interactions between ions and neutrals maximize within the LTI and in particular at altitudes from 100 to 200 km, which is the least visited region of the near-Earth environment. The lack of in situ co-temporal and co-spatial measurements of all relevant parameters and their elusiveness to most remote-sensing methods means that the complex interactions between its neutral and charged constituents remain poorly characterized to this date. This lack of measurements, together with the ambiguity in the quantification of key processes in the 100–200 km altitude range affect current modeling efforts to expand atmospheric models upward to include the LTI and limit current space weather prediction capabilities. We present focused questions in the LTI that are related to the complex interactions between its neutral and charged constituents. These questions concern core physical processes that govern the energetics, dynamics, and chemistry of the LTI and need to be addressed as fundamental and long-standing questions in this critically unexplored boundary region. We also outline the range of in situ measurements that are needed to unambiguously quantify key LTI processes within this region, and present elements of an in situ concept based on past proposed mission concepts.

The lower thermosphere-ionosphere (LTI) is a key transition region between Earth’s atmosphere and space. Interactions between ions and neutrals maximize within the LTI and in particular at altitudes from 100 to 200 km, which is the least visited region of the near-Earth environment. The lack of in situ co-temporal and co-spatial measurements of all relevant parameters and their elusiveness to most remote-sensing methods means that the complex interactions between its neutral and charged constituents remain poorly characterized to this date. This lack of measurements, together with the ambiguity in the quantification of key processes in the 100–200 km altitude range affect current modeling efforts to expand atmospheric models upward to include the LTI and limit current space weather prediction capabilities. We present focused questions in the LTI that are related to the complex interactions between its neutral and charged constituents. These questions concern core physical processes that govern the energetics, dynamics, and chemistry of the LTI and need to be addressed as fundamental and long-standing questions in this critically unexplored boundary region. We also outline the range of in situ measurements that are needed to unambiguously quantify key LTI processes within this region, and present elements of an in situ concept based on past proposed mission concepts.

The lower thermosphere-ionosphere (LTI) is a key transition region between Earth’s atmosphere and space. Interactions between ions and neutrals maximize within the LTI and in particular at altitudes from 100 to 200 km, which is the least visited region of the near-Earth environment. The lack of in situ co-temporal and co-spatial measurements of all relevant parameters and their elusiveness to most remote-sensing methods means that the complex interactions between its neutral and charged constituents remain poorly characterized to this date. This lack of measurements, together with the ambiguity in the quantification of key processes in the 100–200 km altitude range affect current modeling efforts to expand atmospheric models upward to include the LTI and limit current space weather prediction capabilities. We present focused questions in the LTI that are related to the complex interactions between its neutral and charged constituents. These questions concern core physical processes that govern the energetics, dynamics, and chemistry of the LTI and need to be addressed as fundamental and long-standing questions in this critically unexplored boundary region. We also outline the range of in situ measurements that are needed to unambiguously quantify key LTI processes within this region, and present elements of an in situ concept based on past proposed mission concepts.

Daedalus MASE (Mission Assessment through Simulation Exercise) is an open-source package of scientific analysis tools aimed at research in the Lower Thermosphere-Ionosphere (LTI). It was created with the purpose to assess the performance and demonstrate closure of the mission objectives of Daedalus, a mission concept targeting to perform in-situ measurements in the LTI. However, through its successful usage as a mission-simulator toolset, Daedalus MASE has evolved to encompass numerous capabilities related to LTI science and modeling. Inputs are geophysical observables in the LTI, which can be obtained either through in-situ measurements from spacecraft and rockets, or through Global Circulation Models (GCM). These include ion, neutral and electron densities, ion and neutral composition, ion, electron and neutral temperatures, ion drifts, neutral winds, electric field, and magnetic field. In the examples presented, these geophysical observables are obtained through NCAR’s Thermosphere-Ionosphere-Electrodynamics General Circulation Model. Capabilities of Daedalus MASE include: 1) Calculations of products that are derived from the above geophysical observables, such as Joule heating, energy transfer rates between species, electrical currents, electrical conductivity, ion-neutral collision frequencies between all combinations of species, as well as height-integrations of derived products. 2) Calculation and cross-comparison of collision frequencies and estimates of the effect of using different models of collision frequencies into derived products. 3) Calculation of the uncertainties of derived products based on the uncertainties of the geophysical observables, due to instrument errors or to uncertainties in measurement techniques. 4) Routines for the along-orbit interpolation within gridded datasets of GCMs. 5) Routines for the calculation of the global coverage of an in situ mission in regions of interest and for various conditions of solar and geomagnetic activit

The lower thermosphere-ionosphere (LTI) is a key transition region between Earth’s atmosphere and space. Interactions between ions and neutrals maximize within the LTI and in particular at altitudes from 100 to 200 km, which is the least visited region of the near-Earth environment. The lack of in situ co-temporal and co-spatial measurements of all relevant parameters and their elusiveness to most remote-sensing methods means that the complex interactions between its neutral and charged constituents remain poorly characterized to this date. This lack of measurements, together with the ambiguity in the quantification of key processes in the 100 to 200 km altitude range affect current modeling efforts to expand atmospheric models upward to include the LTI and limit current space weather prediction capabilities. We present focused questions in the LTI that are related to the complex interactions between its neutral and charged constituents. These questions concern core physical processes that govern the energetics, dynamics, and chemistry of the LTI and need to be addressed as fundamental and long-standing questions in this critically unexplored boundary region. We also outline the range of in situ measurements that are needed to unambiguously quantify key LTI processes within this region, and present elements of an in situ concept based on past proposed mission concepts.

Daedalus MASE (Mission Assessment through Simulation Exercise) is an open-source package of scientific analysis tools aimed at research in the Lower Thermosphere-Ionosphere (LTI). It was created with the purpose to assess the performance and demonstrate closure of the mission objectives of Daedalus, a mission concept targeting to perform in-situ measurements in the LTI. However, through its successful usage as a mission-simulator toolset, Daedalus MASE has evolved to encompass numerous capabilities related to LTI science and modeling. Inputs are geophysical observables in the LTI, which can be obtained either through in-situ measurements from spacecraft and rockets, or through Global Circulation Models (GCM). These include ion, neutral and electron densities, ion and neutral composition, ion, electron and neutral temperatures, ion drifts, neutral winds, electric field, and magnetic field. In the examples presented, these geophysical observables are obtained through NCAR’s Thermosphere-Ionosphere-Electrodynamics General Circulation Model. Capabilities of Daedalus MASE include: (1) Calculations of products that are derived from the above geophysical observables, such as Joule heating, energy transfer rates between species, electrical currents, electrical conductivity, ion-neutral collision frequencies between all combinations of species, as well as height-integrations of derived products. (2) Calculation and cross-comparison of collision frequencies and estimates of the effect of using different models of collision frequencies into derived products. (3) Calculation of the uncertainties of derived products based on the uncertainties of the geophysical observables, due to instrument errors or to uncertainties in measurement techniques. (4) Routines for the along-orbit interpolation within gridded datasets of GCMs. (5) Routines for the calculation of the global coverage of an in situ mission in regions of interest and for various conditions of solar and geomagnetic activ

Ground-based observations of 11.072 GHz atmospheric ozone (O3) emission have been made using the Ny-Ålesund Ozone in the Mesosphere Instrument (NAOMI) at the UK Arctic Research Station (latitude 78∘55′0′′ N, longitude 11∘55′59′′ E), Spitsbergen. Seasonally averaged O3 vertical profiles in the Arctic polar mesosphere–lower thermosphere region for night-time and twilight conditions in the period 15 August 2017 to 15 March 2020 have been retrieved over the altitude range 62–98 km. NAOMI measurements are compared with corresponding, overlapping observations by the Sounding of the Atmosphere using Broadband Emission Radiometry (SABER) satellite instrument. The NAOMI and SABER version 2.0 data are binned according to the SABER instrument 60 d yaw cycles into nominal 3-month “winter” (15 December–15 March), “autumn” (15 August–15 November), and “summer” (15 April–15 July) periods. The NAOMI observations show the same year-to-year and seasonal variabilities as the SABER 9.6 µm O3 data. The winter night-time (solar zenith angle, SZA ≥ 110∘) and twilight (75∘ ≤ SZA ≤ 110∘) NAOMI and SABER 9.6 µm O3 volume mixing ratio (VMR) profiles agree to within the measurement uncertainties. However, for autumn twilight conditions the SABER 9.6 µm O3 secondary maximum VMR values are higher than NAOMI over altitudes 88–97 km by 47 % and 59 %, respectively in 2017 and 2018. Comparing the two SABER channels which measure O3 at different wavelengths and use different processing schemes, the 9.6 µm O3 autumn twilight VMR data for the three years 2017–2019 are higher than the corresponding 1.27 µm measurements with the largest difference (58 %) in the 65–95 km altitude range similar to the NAOMI observation. The SABER 9.6 µm O3 summer daytime (SZA < 75∘) mesospheric O3 VMR is also consistently higher than the 1.27 µm measurement, confirming previously reported differences between the SABER 9.6 µm channel and measurements of mesospheric O3 by other satellite instruments.

Ground-based observations of 11.072 GHz atmospheric ozone (O3) emission have been made using the Ny-Ålesund Ozone in the Mesosphere Instrument (NAOMI) at the UK Arctic Research Station (latitude 78∘55′0′′ N, longitude 11∘55′59′′ E), Spitsbergen. Seasonally averaged O3 vertical profiles in the Arctic polar mesosphere–lower thermosphere region for night-time and twilight conditions in the period 15 August 2017 to 15 March 2020 have been retrieved over the altitude range 62–98 km. NAOMI measurements are compared with corresponding, overlapping observations by the Sounding of the Atmosphere using Broadband Emission Radiometry (SABER) satellite instrument. The NAOMI and SABER version 2.0 data are binned according to the SABER instrument 60 d yaw cycles into nominal 3-month “winter” (15 December–15 March), “autumn” (15 August–15 November), and “summer” (15 April–15 July) periods. The NAOMI observations show the same year-to-year and seasonal variabilities as the SABER 9.6 µm O3 data. The winter night-time (solar zenith angle, SZA ≥ 110∘) and twilight (75∘ ≤ SZA ≤ 110∘) NAOMI and SABER 9.6 µm O3 volume mixing ratio (VMR) profiles agree to within the measurement uncertainties. However, for autumn twilight conditions the SABER 9.6 µm O3 secondary maximum VMR values are higher than NAOMI over altitudes 88–97 km by 47 % and 59 %, respectively in 2017 and 2018. Comparing the two SABER channels which measure O3 at different wavelengths and use different processing schemes, the 9.6 µm O3 autumn twilight VMR data for the three years 2017–2019 are higher than the corresponding 1.27 µm measurements with the largest difference (58 %) in the 65–95 km altitude range similar to the NAOMI observation. The SABER 9.6 µm O3 summer daytime (SZA < 75∘) mesospheric O3 VMR is also consistently higher than the 1.27 µm measurement, confirming previously reported differences between the SABER 9.6 µm channel and measurements of mesospheric O3 by other satellite instruments.

Ground-based observations of 11.072 GHz atmospheric ozone (O3) emission have been made using the Ny-Ålesund Ozone in the Mesosphere Instrument (NAOMI) at the UK Arctic Research Station (latitude 78∘55′0′′ N, longitude 11∘55′59′′ E), Spitsbergen. Seasonally averaged O3 vertical profiles in the Arctic polar mesosphere–lower thermosphere region for night-time and twilight conditions in the period 15 August 2017 to 15 March 2020 have been retrieved over the altitude range 62–98 km. NAOMI measurements are compared with corresponding, overlapping observations by the Sounding of the Atmosphere using Broadband Emission Radiometry (SABER) satellite instrument. The NAOMI and SABER version 2.0 data are binned according to the SABER instrument 60 d yaw cycles into nominal 3-month “winter” (15 December–15 March), “autumn” (15 August–15 November), and “summer” (15 April–15 July) periods. The NAOMI observations show the same year-to-year and seasonal variabilities as the SABER 9.6 µm O3 data. The winter night-time (solar zenith angle, SZA ≥ 110∘) and twilight (75∘ ≤ SZA ≤ 110∘) NAOMI and SABER 9.6 µm O3 volume mixing ratio (VMR) profiles agree to within the measurement uncertainties. However, for autumn twilight conditions the SABER 9.6 µm O3 secondary maximum VMR values are higher than NAOMI over altitudes 88–97 km by 47 % and 59 %, respectively in 2017 and 2018. Comparing the two SABER channels which measure O3 at different wavelengths and use different processing schemes, the 9.6 µm O3 autumn twilight VMR data for the three years 2017–2019 are higher than the corresponding 1.27 µm measurements with the largest difference (58 %) in the 65–95 km altitude range similar to the NAOMI observation. The SABER 9.6 µm O3 summer daytime (SZA < 75∘) mesospheric O3 VMR is also consistently higher than the 1.27 µm measurement, confirming previously reported differences between the SABER 9.6 µm channel and measurements of mesospheric O3 by other satellite instruments.

Ground-based observations of 11.072 GHz atmospheric ozone (O3) emission have been made using the Ny-Ålesund Ozone in the Mesosphere Instrument (NAOMI) at the UK Arctic Research Station (latitude 78∘55′0′′ N, longitude 11∘55′59′′ E), Spitsbergen. Seasonally averaged O3 vertical profiles in the Arctic polar mesosphere–lower thermosphere region for night-time and twilight conditions in the period 15 August 2017 to 15 March 2020 have been retrieved over the altitude range 62–98 km. NAOMI measurements are compared with corresponding, overlapping observations by the Sounding of the Atmosphere using Broadband Emission Radiometry (SABER) satellite instrument. The NAOMI and SABER version 2.0 data are binned according to the SABER instrument 60 d yaw cycles into nominal 3-month “winter” (15 December–15 March), “autumn” (15 August–15 November), and “summer” (15 April–15 July) periods. The NAOMI observations show the same year-to-year and seasonal variabilities as the SABER 9.6 µm O3 data. The winter night-time (solar zenith angle, SZA ≥ 110∘) and twilight (75∘ ≤ SZA ≤ 110∘) NAOMI and SABER 9.6 µm O3 volume mixing ratio (VMR) profiles agree to within the measurement uncertainties. However, for autumn twilight conditions the SABER 9.6 µm O3 secondary maximum VMR values are higher than NAOMI over altitudes 88–97 km by 47 % and 59 %, respectively in 2017 and 2018. Comparing the two SABER channels which measure O3 at different wavelengths and use different processing schemes, the 9.6 µm O3 autumn twilight VMR data for the three years 2017–2019 are higher than the corresponding 1.27 µm measurements with the largest difference (58 %) in the 65–95 km altitude range similar to the NAOMI observation. The SABER 9.6 µm O3 summer daytime (SZA < 75∘) mesospheric O3 VMR is also consistently higher than the 1.27 µm measurement, confirming previously reported differences between the SABER 9.6 µm channel and measurements of mesospheric O3 by other satellite instruments.

The lower-thermosphere-ionosphere (LTI) system consists of the upper atmosphere and the lower part of the ionosphere and as such comprises a complex system coupled to both the atmosphere below and space above. The atmospheric part of the LTI is dominated by laws of continuum fluid dynamics and chemistry, while the ionosphere is a plasma system controlled by electromagnetic forces driven by the magnetosphere, the solar wind, as well as the wind dynamo. The LTI is hence a domain controlled by many different physical processes. However, systematic in situ measurements within this region are severely lacking, although the LTI is located only 80 to 200 km above the surface of our planet. This paper reviews the current state of the art in measuring the LTI, either in situ or by several different remote-sensing methods. We begin by outlining the open questions within the LTI requiring high-quality in situ measurements, before reviewing directly observable parameters and their most important derivatives. The motivation for this review has arisen from the recent retention of the Daedalus mission as one among three competing mission candidates within the European Space Agency (ESA) Earth Explorer 10 Programme. However, this paper intends to cover the LTI parameters such that it can be used as a background scientific reference for any mission targeting in situ observations of the LTI.

The lower-thermosphere-ionosphere (LTI) system consists of the upper atmosphere and the lower part of the ionosphere and as such comprises a complex system coupled to both the atmosphere below and space above. The atmospheric part of the LTI is dominated by laws of continuum fluid dynamics and chemistry, while the ionosphere is a plasma system controlled by electromagnetic forces driven by the magnetosphere, the solar wind, as well as the wind dynamo. The LTI is hence a domain controlled by many different physical processes. However, systematic in situ measurements within this region are severely lacking, although the LTI is located only 80 to 200 km above the surface of our planet. This paper reviews the current state of the art in measuring the LTI, either in situ or by several different remote-sensing methods. We begin by outlining the open questions within the LTI requiring high-quality in situ measurements, before reviewing directly observable parameters and their most important derivatives. The motivation for this review has arisen from the recent retention of the Daedalus mission as one among three competing mission candidates within the European Space Agency (ESA) Earth Explorer 10 Programme. However, this paper intends to cover the LTI parameters such that it can be used as a background scientific reference for any mission targeting in situ observations of the LTI.

The lower-thermosphere-ionosphere (LTI) system consists of the upper atmosphere and the lower part of the ionosphere and as such comprises a complex system coupled to both the atmosphere below and space above. The atmospheric part of the LTI is dominated by laws of continuum fluid dynamics and chemistry, while the ionosphere is a plasma system controlled by electromagnetic forces driven by the magnetosphere, the solar wind, as well as the wind dynamo. The LTI is hence a domain controlled by many different physical processes. However, systematic in situ measurements within this region are severely lacking, although the LTI is located only 80 to 200 km above the surface of our planet. This paper reviews the current state of the art in measuring the LTI, either in situ or by several different remote-sensing methods. We begin by outlining the open questions within the LTI requiring high-quality in situ measurements, before reviewing directly observable parameters and their most important derivatives. The motivation for this review has arisen from the recent retention of the Daedalus mission as one among three competing mission candidates within the European Space Agency (ESA) Earth Explorer 10 Programme. However, this paper intends to cover the LTI parameters such that it can be used as a background scientific reference for any mission targeting in situ observations of the LTI.

The lower-thermosphere-ionosphere (LTI) system consists of the upper atmosphere and the lower part of the ionosphere and as such comprises a complex system coupled to both the atmosphere below and space above. The atmospheric part of the LTI is dominated by laws of continuum fluid dynamics and chemistry, while the ionosphere is a plasma system controlled by electromagnetic forces driven by the magnetosphere, the solar wind, as well as the wind dynamo. The LTI is hence a domain controlled by many different physical processes. However, systematic in situ measurements within this region are severely lacking, although the LTI is located only 80 to 200 km above the surface of our planet. This paper reviews the current state of the art in measuring the LTI, either in situ or by several different remote-sensing methods. We begin by outlining the open questions within the LTI requiring high-quality in situ measurements, before reviewing directly observable parameters and their most important derivatives. The motivation for this review has arisen from the recent retention of the Daedalus mission as one among three competing mission candidates within the European Space Agency (ESA) Earth Explorer 10 Programme. However, this paper intends to cover the LTI parameters such that it can be used as a background scientific reference for any mission targeting in situ observations of the LTI.

The lower-thermosphere-ionosphere (LTI) system consists of the upper atmosphere and the lower part of the ionosphere and as such comprises a complex system coupled to both the atmosphere below and space above. The atmospheric part of the LTI is dominated by laws of continuum fluid dynamics and chemistry, while the ionosphere is a plasma system controlled by electromagnetic forces driven by the magnetosphere, the solar wind, as well as the wind dynamo. The LTI is hence a domain controlled by many different physical processes. However, systematic in situ measurements within this region are severely lacking, although the LTI is located only 80 to 200 km above the surface of our planet. This paper reviews the current state of the art in measuring the LTI, either in situ or by several different remote-sensing methods. We begin by outlining the open questions within the LTI requiring high-quality in situ measurements, before reviewing directly observable parameters and their most important derivatives. The motivation for this review has arisen from the recent retention of the Daedalus mission as one among three competing mission candidates within the European Space Agency (ESA) Earth Explorer 10 Programme. However, this paper intends to cover the LTI parameters such that it can be used as a background scientific reference for any mission targeting in situ observations of the LTI.

The lower-Thermosphere-ionosphere (LTI) system consists of the upper atmosphere and the lower part of the ionosphere and as such comprises a complex system coupled to both the atmosphere below and space above. The atmospheric part of the LTI is dominated by laws of continuum fluid dynamics and chemistry, while the ionosphere is a plasma system controlled by electromagnetic forces driven by the magnetosphere, the solar wind, as well as the wind dynamo. The LTI is hence a domain controlled by many different physical processes. However, systematic in situ measurements within this region are severely lacking, although the LTI is located only 80 to 200 km above the surface of our planet. This paper reviews the current state of the art in measuring the LTI, either in situ or by several different remote-sensing methods. We begin by outlining the open questions within the LTI requiring high-quality in situ measurements, before reviewing directly observable parameters and their most important derivatives. The motivation for this review has arisen from the recent retention of the Daedalus mission as one among three competing mission candidates within the European Space Agency (ESA) Earth Explorer 10 Programme. However, this paper intends to cover the LTI parameters such that it can be used as a background scientific reference for any mission targeting in situ observations of the LTI..

Recently, a model for medium energy (30–1000 keV) radiation belt‐driven electron precipitation (ApEEP) has been put forward for use in decadal to century‐long climate model runs as part of the Climate Modelling Intercomparison Project, phase 6 (CMIP6). The ApEEP model is based on directly observed precipitation data spanning 2002‐2012 from the constellation of low Earth orbiting Polar Operational Environmental Satellites (POES). Here we test the ApEEP model's ability using its magnetic local time variant, ApEEP_MLT, to accurately represent electron precipitation fluxes from the radiation belts during a large geomagnetic storm that occurred outside of the span of the development dataset. In a study of narrow band sub‐ionospheric VLF transmitter data collected during March 2015, continuous phase observations have been analyzed throughout the entire St. Patrick's Day geomagnetic storm period for the first time. Using phase data from the UK transmitter, call‐sign GVT (22.1 kHz), received in Reykjavik, Iceland, electron precipitation fluxes from L=2.8‐5.4 are calculated around magnetic local noon (12 MLT), and magnetic midnight (00 MLT). VLF‐inferred >30 keV fluxes are similar to the equivalent directly‐observed POES fluxes. The ApEEP_MLT >30 keV fluxes for L<5.5 describe the overall St Patrick's Day geomagnetic storm‐driven flux enhancement well, although they are a factor of 1.7 (1.3) lower than POES (VLF‐inferred) fluxes during the recovery phase.. Such close agreement in >30 keV flux levels during a large geomagnetic storm, using three different techniques, indicates this flux forcing are appropriate for decadal climate simulations for which the ApEEP model was created.

Recently, a model for medium energy (30–1000 keV) radiation belt‐driven electron precipitation (ApEEP) has been put forward for use in decadal to century‐long climate model runs as part of the Climate Modelling Intercomparison Project, phase 6 (CMIP6). The ApEEP model is based on directly observed precipitation data spanning 2002‐2012 from the constellation of low Earth orbiting Polar Operational Environmental Satellites (POES). Here we test the ApEEP model's ability using its magnetic local time variant, ApEEP_MLT, to accurately represent electron precipitation fluxes from the radiation belts during a large geomagnetic storm that occurred outside of the span of the development dataset. In a study of narrow band sub‐ionospheric VLF transmitter data collected during March 2015, continuous phase observations have been analyzed throughout the entire St. Patrick's Day geomagnetic storm period for the first time. Using phase data from the UK transmitter, call‐sign GVT (22.1 kHz), received in Reykjavik, Iceland, electron precipitation fluxes from L=2.8‐5.4 are calculated around magnetic local noon (12 MLT), and magnetic midnight (00 MLT). VLF‐inferred >30 keV fluxes are similar to the equivalent directly‐observed POES fluxes. The ApEEP_MLT >30 keV fluxes for L<5.5 describe the overall St Patrick's Day geomagnetic storm‐driven flux enhancement well, although they are a factor of 1.7 (1.3) lower than POES (VLF‐inferred) fluxes during the recovery phase.. Such close agreement in >30 keV flux levels during a large geomagnetic storm, using three different techniques, indicates this flux forcing are appropriate for decadal climate simulations for which the ApEEP model was created.

The Ku-band microwave frequencies (10.70–14.25 GHz) overlap emissions from ozone (O3) at 11.072 GHz and hydroxyl radical (OH) at 13.441 GHz. These important chemical species in the polar middle atmosphere respond strongly to high-latitude geomagnetic activity associated with space weather. Atmospheric model calculations predict that energetic electron precipitation (EEP) driven by magnetospheric substorms produces large changes in polar mesospheric O3 and OH. The EEP typically peaks at geomagnetic latitudes of ∼65∘ and evolves rapidly with time longitudinally and over the geomagnetic latitude range 60–80∘. Previous atmospheric modelling studies have shown that during substorms OH abundance can increase by more than an order of magnitude at 64–84 km and mesospheric O3 losses can exceed 50 %. In this work, an atmospheric simulation and retrieval study has been performed to determine the requirements for passive microwave radiometers capable of measuring diurnal variations in O3 and OH profiles from high-latitude Northern Hemisphere and Antarctic locations to verify model predictions. We show that, for a 11.072 GHz radiometer making 6 h spectral measurements with 10 kHz frequency resolution and root-mean-square baseline noise of 1 mK, O3 could be profiled over 8×10−4–0.22 hPa (∼98–58 km) with 10–17 km height resolution and ∼1 ppmv uncertainty. For the equivalent 13.441 GHz measurements with vertical sensor polarisation, OH could be profiled over 3×10−3–0.29 hPa (∼90–56 km) with 10–17 km height resolution and ∼3 ppbv uncertainty. The proposed observations would be highly applicable to studies of EEP, atmospheric dynamics, planetary-scale circulation, chemical transport, and the representation of these processes in polar and global climate models. Such observations would provide a relatively low-cost alternative to increasingly sparse satellite measurements of the polar middle atmosphere, extending long-term data records and also providing “ground truth” calibration data

The Ku-band microwave frequencies (10.70–14.25 GHz) overlap emissions from ozone (O3) at 11.072 GHz and hydroxyl radical (OH) at 13.441 GHz. These important chemical species in the polar middle atmosphere respond strongly to high-latitude geomagnetic activity associated with space weather. Atmospheric model calculations predict that energetic electron precipitation (EEP) driven by magnetospheric substorms produces large changes in polar mesospheric O3 and OH. The EEP typically peaks at geomagnetic latitudes of ∼65∘ and evolves rapidly with time longitudinally and over the geomagnetic latitude range 60–80∘. Previous atmospheric modelling studies have shown that during substorms OH abundance can increase by more than an order of magnitude at 64–84 km and mesospheric O3 losses can exceed 50 %. In this work, an atmospheric simulation and retrieval study has been performed to determine the requirements for passive microwave radiometers capable of measuring diurnal variations in O3 and OH profiles from high-latitude Northern Hemisphere and Antarctic locations to verify model predictions. We show that, for a 11.072 GHz radiometer making 6 h spectral measurements with 10 kHz frequency resolution and root-mean-square baseline noise of 1 mK, O3 could be profiled over 8×10−4–0.22 hPa (∼98–58 km) with 10–17 km height resolution and ∼1 ppmv uncertainty. For the equivalent 13.441 GHz measurements with vertical sensor polarisation, OH could be profiled over 3×10−3–0.29 hPa (∼90–56 km) with 10–17 km height resolution and ∼3 ppbv uncertainty. The proposed observations would be highly applicable to studies of EEP, atmospheric dynamics, planetary-scale circulation, chemical transport, and the representation of these processes in polar and global climate models. Such observations would provide a relatively low-cost alternative to increasingly sparse satellite measurements of the polar middle atmosphere, extending long-term data records and also providing “ground truth” calibration data

While observed mesospheric polar nitric acid enhancements have been attributed to energetic particle precipitation through ion cluster chemistry in the past, this phenomenon is not reproduced in current whole-atmosphere chemistry-climate models. We investigate such nitric acid enhancements resulting from energetic electron precipitation events using a recently developed variant of the Whole Atmosphere Community Climate Model (WACCM) that includes a sophisticated ion chemistry tailored for the D-layer of the ionosphere (50-90 km), namely WACCM-D. Using the specified-dynamics mode, i.e., nudging dynamics in the troposphere and stratosphere to meteorological re-analyses, we perform a one-year long simulation (July 2009-June 2010) and contrast WACCM-D with the standard WACCM. Both WACCM and WACCM-D simulations are performed with and without forcing from medium-to-high energy electron precipitation, allowing a better representation of the energetic electrons penetrating into the mesosphere. We demonstrate the effects of the strong particle precipitation events which occurred during April and May 2010 on nitric acid and on key ion cluster species, as well as other relevant species of the nitrogen family. The one-year-long simulation allows the event-related changes in neutral and ionic species to be placed in the context of their annual cycle. We especially highlight the role played by medium-to-high energy electrons in triggering ion cluster chemistry and ion-ion recombinations in the mesosphere and lower thermosphere during the precipitation event, leading to enhanced production of nitric acid and raising its abundance by two orders of magnitude from 10-4 to a few 10-2 ppb.

While observed mesospheric polar nitric acid enhancements have been attributed to energetic particle precipitation through ion cluster chemistry in the past, this phenomenon is not reproduced in current whole-atmosphere chemistry-climate models. We investigate such nitric acid enhancements resulting from energetic electron precipitation events using a recently developed variant of the Whole Atmosphere Community Climate Model (WACCM) that includes a sophisticated ion chemistry tailored for the D-layer of the ionosphere (50-90 km), namely WACCM-D. Using the specified-dynamics mode, i.e., nudging dynamics in the troposphere and stratosphere to meteorological re-analyses, we perform a one-year long simulation (July 2009-June 2010) and contrast WACCM-D with the standard WACCM. Both WACCM and WACCM-D simulations are performed with and without forcing from medium-to-high energy electron precipitation, allowing a better representation of the energetic electrons penetrating into the mesosphere. We demonstrate the effects of the strong particle precipitation events which occurred during April and May 2010 on nitric acid and on key ion cluster species, as well as other relevant species of the nitrogen family. The one-year-long simulation allows the event-related changes in neutral and ionic species to be placed in the context of their annual cycle. We especially highlight the role played by medium-to-high energy electrons in triggering ion cluster chemistry and ion-ion recombinations in the mesosphere and lower thermosphere during the precipitation event, leading to enhanced production of nitric acid and raising its abundance by two orders of magnitude from 10-4 to a few 10-2 ppb.

This paper describes the recommended solar forcing dataset for CMIP6 and highlights changes with respect to CMIP5. The solar forcing is provided for radiative properties, namely total solar irradiance (TSI), solar spectral irradiance (SSI), and the F10.7 index as well as particle forcing, including geomagnetic indices Ap and Kp, and ionization rates to account for effects of solar protons, electrons, and galactic cosmic rays. This is the first time that a recommendation for solar-driven particle forcing has been provided for a CMIP exercise. The solar forcing datasets are provided at daily and monthly resolution separately for the CMIP6 preindustrial control, historical (1850–2014), and future (2015–2300) simulations. For the preindustrial control simulation, both constant and time-varying solar forcing components are provided, with the latter including variability on 11-year and shorter timescales but no long-term changes. For the future, we provide a realistic scenario of what solar behavior could be, as well as an additional extreme Maunder-minimum-like sensitivity scenario. This paper describes the forcing datasets and also provides detailed recommendations as to their implementation in current climate models. For the historical simulations, the TSI and SSI time series are defined as the average of two solar irradiance models that are adapted to CMIP6 needs: an empirical one (NRLTSI2–NRLSSI2) and a semi-empirical one (SATIRE). A new and lower TSI value is recommended: the contemporary solar-cycle average is now 1361.0 W m−2. The slight negative trend in TSI over the three most recent solar cycles in the CMIP6 dataset leads to only a small global radiative forcing of −0.04 W m−2. In the 200–400 nm wavelength range, which is important for ozone photochemistry, the CMIP6 solar forcing dataset shows a larger solar-cycle variability contribution to TSI than in CMIP5 (50 % compared to 35 %). We compare the climatic effects of the CMIP6 solar forcing dataset to its CMI

This paper describes the recommended solar forcing dataset for CMIP6 and highlights changes with respect to CMIP5. The solar forcing is provided for radiative properties, namely total solar irradiance (TSI), solar spectral irradiance (SSI), and the F10.7 index as well as particle forcing, including geomagnetic indices Ap and Kp, and ionization rates to account for effects of solar protons, electrons, and galactic cosmic rays. This is the first time that a recommendation for solar-driven particle forcing has been provided for a CMIP exercise. The solar forcing datasets are provided at daily and monthly resolution separately for the CMIP6 preindustrial control, historical (1850–2014), and future (2015–2300) simulations. For the preindustrial control simulation, both constant and time-varying solar forcing components are provided, with the latter including variability on 11-year and shorter timescales but no long-term changes. For the future, we provide a realistic scenario of what solar behavior could be, as well as an additional extreme Maunder-minimum-like sensitivity scenario. This paper describes the forcing datasets and also provides detailed recommendations as to their implementation in current climate models. For the historical simulations, the TSI and SSI time series are defined as the average of two solar irradiance models that are adapted to CMIP6 needs: an empirical one (NRLTSI2–NRLSSI2) and a semi-empirical one (SATIRE). A new and lower TSI value is recommended: the contemporary solar-cycle average is now 1361.0 W m−2. The slight negative trend in TSI over the three most recent solar cycles in the CMIP6 dataset leads to only a small global radiative forcing of −0.04 W m−2. In the 200–400 nm wavelength range, which is important for ozone photochemistry, the CMIP6 solar forcing dataset shows a larger solar-cycle variability contribution to TSI than in CMIP5 (50 % compared to 35 %). We compare the climatic effects of the CMIP6 solar forcing dataset to its CMI

This paper describes the recommended solar forcing dataset for CMIP6 and highlights changes with respect to CMIP5. The solar forcing is provided for radiative properties, namely total solar irradiance (TSI), solar spectral irradiance (SSI), and the F10.7 index as well as particle forcing, including geomagnetic indices Ap and Kp, and ionization rates to account for effects of solar protons, electrons, and galactic cosmic rays. This is the first time that a recommendation for solar-driven particle forcing has been provided for a CMIP exercise. The solar forcing datasets are provided at daily and monthly resolution separately for the CMIP6 preindustrial control, historical (1850–2014), and future (2015–2300) simulations. For the preindustrial control simulation, both constant and time-varying solar forcing components are provided, with the latter including variability on 11-year and shorter timescales but no long-term changes. For the future, we provide a realistic scenario of what solar behavior could be, as well as an additional extreme Maunder-minimum-like sensitivity scenario. This paper describes the forcing datasets and also provides detailed recommendations as to their implementation in current climate models. For the historical simulations, the TSI and SSI time series are defined as the average of two solar irradiance models that are adapted to CMIP6 needs: an empirical one (NRLTSI2–NRLSSI2) and a semi-empirical one (SATIRE). A new and lower TSI value is recommended: the contemporary solar-cycle average is now 1361.0 W m−2. The slight negative trend in TSI over the three most recent solar cycles in the CMIP6 dataset leads to only a small global radiative forcing of −0.04 W m−2. In the 200–400 nm wavelength range, which is important for ozone photochemistry, the CMIP6 solar forcing dataset shows a larger solar-cycle variability contribution to TSI than in CMIP5 (50 % compared to 35 %). We compare the climatic effects of the CMIP6 solar forcing dataset to its CMI

This paper describes the recommended solar forcing dataset for CMIP6 and highlights changes with respect to CMIP5. The solar forcing is provided for radiative properties, namely total solar irradiance (TSI), solar spectral irradiance (SSI), and the F10.7 index as well as particle forcing, including geomagnetic indices Ap and Kp, and ionization rates to account for effects of solar protons, electrons, and galactic cosmic rays. This is the first time that a recommendation for solar-driven particle forcing has been provided for a CMIP exercise. The solar forcing datasets are provided at daily and monthly resolution separately for the CMIP6 preindustrial control, historical (1850–2014), and future (2015–2300) simulations. For the preindustrial control simulation, both constant and time-varying solar forcing components are provided, with the latter including variability on 11-year and shorter timescales but no long-term changes. For the future, we provide a realistic scenario of what solar behavior could be, as well as an additional extreme Maunder-minimum-like sensitivity scenario. This paper describes the forcing datasets and also provides detailed recommendations as to their implementation in current climate models. For the historical simulations, the TSI and SSI time series are defined as the average of two solar irradiance models that are adapted to CMIP6 needs: an empirical one (NRLTSI2–NRLSSI2) and a semi-empirical one (SATIRE). A new and lower TSI value is recommended: the contemporary solar-cycle average is now 1361.0 W m−2. The slight negative trend in TSI over the three most recent solar cycles in the CMIP6 dataset leads to only a small global radiative forcing of −0.04 W m−2. In the 200–400 nm wavelength range, which is important for ozone photochemistry, the CMIP6 solar forcing dataset shows a larger solar-cycle variability contribution to TSI than in CMIP5 (50 % compared to 35 %). We compare the climatic effects of the CMIP6 solar forcing dataset to its CMI

This paper describes the recommended solar forcing dataset for CMIP6 and highlights changes with respect to CMIP5. The solar forcing is provided for radiative properties, namely total solar irradiance (TSI), solar spectral irradiance (SSI), and the F10.7 index as well as particle forcing, including geomagnetic indices Ap and Kp, and ionization rates to account for effects of solar protons, electrons, and galactic cosmic rays. This is the first time that a recommendation for solar-driven particle forcing has been provided for a CMIP exercise. The solar forcing datasets are provided at daily and monthly resolution separately for the CMIP6 preindustrial control, historical (1850–2014), and future (2015–2300) simulations. For the preindustrial control simulation, both constant and time-varying solar forcing components are provided, with the latter including variability on 11-year and shorter timescales but no long-term changes. For the future, we provide a realistic scenario of what solar behavior could be, as well as an additional extreme Maunder-minimum-like sensitivity scenario. This paper describes the forcing datasets and also provides detailed recommendations as to their implementation in current climate models. For the historical simulations, the TSI and SSI time series are defined as the average of two solar irradiance models that are adapted to CMIP6 needs: an empirical one (NRLTSI2–NRLSSI2) and a semi-empirical one (SATIRE). A new and lower TSI value is recommended: the contemporary solar-cycle average is now 1361.0 W m−2. The slight negative trend in TSI over the three most recent solar cycles in the CMIP6 dataset leads to only a small global radiative forcing of −0.04 W m−2. In the 200–400 nm wavelength range, which is important for ozone photochemistry, the CMIP6 solar forcing dataset shows a larger solar-cycle variability contribution to TSI than in CMIP5 (50 % compared to 35 %). We compare the climatic effects of the CMIP6 solar forcing dataset to its CMI

This paper describes the recommended solar forcing dataset for CMIP6 and highlights changes with respect to CMIP5. The solar forcing is provided for radiative properties, namely total solar irradiance (TSI), solar spectral irradiance (SSI), and the F10.7 index as well as particle forcing, including geomagnetic indices Ap and Kp, and ionization rates to account for effects of solar protons, electrons, and galactic cosmic rays. This is the first time that a recommendation for solar-driven particle forcing has been provided for a CMIP exercise. The solar forcing datasets are provided at daily and monthly resolution separately for the CMIP6 preindustrial control, historical (1850–2014), and future (2015–2300) simulations. For the preindustrial control simulation, both constant and time-varying solar forcing components are provided, with the latter including variability on 11-year and shorter timescales but no long-term changes. For the future, we provide a realistic scenario of what solar behavior could be, as well as an additional extreme Maunder-minimum-like sensitivity scenario. This paper describes the forcing datasets and also provides detailed recommendations as to their implementation in current climate models. For the historical simulations, the TSI and SSI time series are defined as the average of two solar irradiance models that are adapted to CMIP6 needs: an empirical one (NRLTSI2–NRLSSI2) and a semi-empirical one (SATIRE). A new and lower TSI value is recommended: the contemporary solar-cycle average is now 1361.0 W m−2. The slight negative trend in TSI over the three most recent solar cycles in the CMIP6 dataset leads to only a small global radiative forcing of −0.04 W m−2. In the 200–400 nm wavelength range, which is important for ozone photochemistry, the CMIP6 solar forcing dataset shows a larger solar-cycle variability contribution to TSI than in CMIP5 (50 % compared to 35 %). We compare the climatic effects of the CMIP6 solar forcing dataset to its CMI

This study presents a new ion-neutral chemical model coupled into the Whole Atmosphere Community Climate Model (WACCM). The ionospheric D region (altitudes ~ 50–90 km) chemistry is based on the Sodankylä Ion and Neutral Chemistry (SIC) model, a 1-dimensional model containing 306 ion-neutral and ionrecombination reactions of neutral species, positive and negative ions, and electrons. The SIC mechanism was reduced using the Simulation Error Minimization Connectivity Method (SEM-CM) to produce a reaction scheme of 181 ion-molecule reactions. This scheme describes the concentration profiles at altitudes between 20 km and 120 km of a set of major neutral species (HNO3, O3, H2O2, NO, NO2, HO2, OH, N2O5) and ions (O2+, O4+, NO+, NO+(H2O), O2+(H2O), H+(H2O), H+(H2O)2, H+(H2O)3, H+(H2O)4, O3−, NO2−, O−, O2, OH−, O2−(H2O), O2−(H2O)2, O4−, CO3−, CO3−(H2O), CO4−, HCO3−, NO2−, NO3−, NO3−(H2O), NO3(H2O)2, NO3−(HNO3), NO3−(HNO3)2, Cl−, ClO−), which agree with the full SIC mechanism within a 5 % tolerance. Four 3D model simulations were then performed, using the impact of the January 2005 Solar Proton Event (SPE) on D region HOx and NOx chemistry as a test case of four different model versions: the standard WACCM (no negative ions and a very limited set of positive ions); WACCM-SIC (standard WACCM with the full SIC chemistry of positive and negative ions); WACCM-D (standard WACCM with a heuristic reduction of the SIC chemistry, recently used to examine HNO3 formation following an SPE); and WACCM-rSIC (standard WACCM with a reduction of SIC chemistry using the SEM-CM Method). Standard WACCM misses the HNO3 enhancement during the SPE, while the full and reduced model versions predict significant NOx, HOx and HNO3 enhancements in the mesosphere during solar proton events. The SEM-CM reduction also identifies the important ion-molecule reactions that affect the partitioning of odd nitrogen (NOx), odd hydrogen (HOx), and O3 in the stratosphere and mesosphere.

This study presents a new ion-neutral chemical model coupled into the Whole Atmosphere Community Climate Model (WACCM). The ionospheric D region (altitudes ~ 50–90 km) chemistry is based on the Sodankylä Ion and Neutral Chemistry (SIC) model, a 1-dimensional model containing 306 ion-neutral and ionrecombination reactions of neutral species, positive and negative ions, and electrons. The SIC mechanism was reduced using the Simulation Error Minimization Connectivity Method (SEM-CM) to produce a reaction scheme of 181 ion-molecule reactions. This scheme describes the concentration profiles at altitudes between 20 km and 120 km of a set of major neutral species (HNO3, O3, H2O2, NO, NO2, HO2, OH, N2O5) and ions (O2+, O4+, NO+, NO+(H2O), O2+(H2O), H+(H2O), H+(H2O)2, H+(H2O)3, H+(H2O)4, O3−, NO2−, O−, O2, OH−, O2−(H2O), O2−(H2O)2, O4−, CO3−, CO3−(H2O), CO4−, HCO3−, NO2−, NO3−, NO3−(H2O), NO3(H2O)2, NO3−(HNO3), NO3−(HNO3)2, Cl−, ClO−), which agree with the full SIC mechanism within a 5 % tolerance. Four 3D model simulations were then performed, using the impact of the January 2005 Solar Proton Event (SPE) on D region HOx and NOx chemistry as a test case of four different model versions: the standard WACCM (no negative ions and a very limited set of positive ions); WACCM-SIC (standard WACCM with the full SIC chemistry of positive and negative ions); WACCM-D (standard WACCM with a heuristic reduction of the SIC chemistry, recently used to examine HNO3 formation following an SPE); and WACCM-rSIC (standard WACCM with a reduction of SIC chemistry using the SEM-CM Method). Standard WACCM misses the HNO3 enhancement during the SPE, while the full and reduced model versions predict significant NOx, HOx and HNO3 enhancements in the mesosphere during solar proton events. The SEM-CM reduction also identifies the important ion-molecule reactions that affect the partitioning of odd nitrogen (NOx), odd hydrogen (HOx), and O3 in the stratosphere and mesosphere.

Observations have shown that mesospheric hydroxyl (OH) is affected by energetic electron precipitation (EEP) at magnetic latitudes connected to the outer radiation belt. It is not clear, however, if the current satellite-based electron flux observations can be used to accurately describe EEP in atmospheric models. We use the Sodankylä Ion and Neutral Chemistry (SIC) model to reproduce the changes in OH and ozone observed by the Microwave Limb Sounder (MLS/Aura) during four strong EEP events. The daily mean electron energy-flux spectrum, needed for ionization rate calculations, is determined by combining the Medium Energy Proton and Electron Detector fluxes and spectral form from the instrument for the detection of particles high-energy electron detector on board the DEMETER satellite. We show that in general SIC is able to reproduce the observed day-to-day variability of OH and ozone. In the lower mesosphere, the model tends to underestimate the OH concentration, possibly because of uncertainties in the electron spectra for energies >300 keV. The model predicts OH increases at 60–80 km, reaching several hundred percent at 70–80 km during peak EEP forcing. Increases in OH are followed by ozone depletion, up to several tens of percent. The magnitude of modeled changes is similar to those observed by MLS and comparable to effects of individual solar proton events. Our results suggest that the combined satellite observations of electrons can be used to model the EEP effects above 70 km during geomagnetic storms, without a need for significant adjustments. However, for EEP energies >300 keV impacting altitudes <70 km, correction factors may be required.

Observations have shown that mesospheric hydroxyl (OH) is affected by energetic electron precipitation (EEP) at magnetic latitudes connected to the outer radiation belt. It is not clear, however, if the current satellite-based electron flux observations can be used to accurately describe EEP in atmospheric models. We use the Sodankylä Ion and Neutral Chemistry (SIC) model to reproduce the changes in OH and ozone observed by the Microwave Limb Sounder (MLS/Aura) during four strong EEP events. The daily mean electron energy-flux spectrum, needed for ionization rate calculations, is determined by combining the Medium Energy Proton and Electron Detector fluxes and spectral form from the instrument for the detection of particles high-energy electron detector on board the DEMETER satellite. We show that in general SIC is able to reproduce the observed day-to-day variability of OH and ozone. In the lower mesosphere, the model tends to underestimate the OH concentration, possibly because of uncertainties in the electron spectra for energies >300 keV. The model predicts OH increases at 60–80 km, reaching several hundred percent at 70–80 km during peak EEP forcing. Increases in OH are followed by ozone depletion, up to several tens of percent. The magnitude of modeled changes is similar to those observed by MLS and comparable to effects of individual solar proton events. Our results suggest that the combined satellite observations of electrons can be used to model the EEP effects above 70 km during geomagnetic storms, without a need for significant adjustments. However, for EEP energies >300 keV impacting altitudes <70 km, correction factors may be required.

Energetic particle precipitation leads to enhancement of odd hydrogen (HOx) below 80 km altitude through water cluster ion chemistry. Using measurements from the Microwave Limb Sounder (MLS/Aura) and Medium Energy Proton and Electron Detector (MEPED/POES) between 2004–2009, we study variations of nighttime OH caused by radiation belt electrons at geomagnetic latitudes 55–65°. For those months with daily mean 100–300 keV electron count rate exceeding 150 counts/s in the outer radiation belt, we find a strong correlation (r ≥ 0.6) between OH mixing ratios at 70–78 km (0.046–0.015 hPa) and precipitating electrons. Correlations r ≥ 0.35, corresponding to random chance probability p ≤ 5%, are observed down to52 km (0.681 hPa), while no clear correlation is observed at altitudes below. This suggests that the fluxes of ≥3 MeV electrons were not high enough to cause observable changes in OH mixing ratios. At 75 km, in about 34% of the 65 months analyzed we find a correlation r ≥ 0.35. Although similar results are obtained for both hemispheres in general, in some cases the differences in atmospheric conditions make the OH response more difficult to detect in the South. Considering the latitude extent of electron forcing, we find clear effects on OH at magnetic latitudes 55–72°, while the lower latitudes are influenced much less. Because the time period 2004–2009 analyzed here coincided with an extended solar minimum, and the year 2009 was anomalously quiet, it is reasonable to assume that our results provide a lower-limit estimation of the importance of energetic electron precipitation at the latitudes considered.

In order to make best use of the opportunities provided by space missions such as the Radiation Belt Storm Probes, we determine the response of complementary subionospheric radiowave propagation measurements (VLF), riometer absorption measurements, cosmic noise absorption, and GPS-produced total electron content (vTEC) to different energetic electron precipitation (EEP). We model the relative sensitivity and responses of these instruments to idealized monoenergetic beams of precipitating electrons, and more realistic EEP spectra chosen to represent radiation belts and substorm precipitation. In the monoenergetic beam case, we find riometers are more sensitive to the same EEP event occurring during the day than during the night, while subionospheric VLF shows the opposite relationship, and the change in vTEC is independent. In general, the subionospheric VLF measurements are much more sensitive than the other two techniques for EEP over 200 keV, responding to flux magnitudes two-three orders of magnitude smaller than detectable by a riometer. Detectable TEC changes only occur for extreme monoenergetic fluxes. For the radiation belt EEP case, clearly detectable subionospheric VLF responses are produced by daytime fluxes that are ∼10 times lower than required for riometers, while nighttime fluxes can be 10,000 times lower. Riometers are likely to respond only to radiation belt fluxes during the largest EEP events and vTEC is unlikely to be significantly disturbed by radiation belt EEP. For the substorm EEP case both the riometer absorption and the subionospheric VLF technique respond significantly, as does the change in vTEC, which is likely to be detectable at ∼3–4 total electron content units.

Using ground-based subionospheric radio wave propagation data from two very low frequency (VLF) receiver sites, riometer absorption data, and THEMIS satellite observations, we examine in detail energetic electron precipitation (EEP) characteristics associated with two substorm precipitation events that occurred on 28 May 2010. In an advance on the analysis undertaken by Clilverd et al. (2008), we use phase observations of VLF radio wave signals to describe substorm-driven EEP characteristics more accurately than before. Using a >30 keV electron precipitation flux of 5.6 X 107 el. cm-2 sr-1 s-1 and a spectral gradient consistent with that observed by THEMIS, it was possible to accurately reproduce the peak observed riometer absorption at Macquarie Island (L = 5.4) and the associated NWC radio wave phase change observed at Casey, Antarctica, during the second, larger substorm. The flux levels were near to 80% of the peak fluxes observed in a similar substorm as studied by Clilverd et al. (2008). During the initial stages of the second substorm, a latitude region of 5 < L < 9 was affected by electron precipitation. Both substorms showed expansion of the precipitation region to 4 < L < 12 more than 30 min after the injection. While both substorms occurred at similar local times, with electron precipitation injections into approximately the same geographical region, the second expanded in an eastward longitude more slowly, suggesting the involvement of lower-energy electron precipitation. Each substorm region expanded westward at a rate slower than that exhibited eastward. This study shows that it is possible to successfully combine these multi-instrument observations to investigate the characteristics of substorms.

Energetic particle precipitation leads to enhancement of odd hydrogen (HOx) below 80 km altitude through water cluster ion chemistry. Using measurements from the Microwave Limb Sounder (MLS/Aura) and Medium Energy Proton and Electron Detector (MEPED/POES) between 2004–2009, we study variations of nighttime OH caused by radiation belt electrons at geomagnetic latitudes 55–65°. For those months with daily mean 100–300 keV electron count rate exceeding 150 counts/s in the outer radiation belt, we find a strong correlation (r ≥ 0.6) between OH mixing ratios at 70–78 km (0.046–0.015 hPa) and precipitating electrons. Correlations r ≥ 0.35, corresponding to random chance probability p ≤ 5%, are observed down to52 km (0.681 hPa), while no clear correlation is observed at altitudes below. This suggests that the fluxes of ≥3 MeV electrons were not high enough to cause observable changes in OH mixing ratios. At 75 km, in about 34% of the 65 months analyzed we find a correlation r ≥ 0.35. Although similar results are obtained for both hemispheres in general, in some cases the differences in atmospheric conditions make the OH response more difficult to detect in the South. Considering the latitude extent of electron forcing, we find clear effects on OH at magnetic latitudes 55–72°, while the lower latitudes are influenced much less. Because the time period 2004–2009 analyzed here coincided with an extended solar minimum, and the year 2009 was anomalously quiet, it is reasonable to assume that our results provide a lower-limit estimation of the importance of energetic electron precipitation at the latitudes considered.

In order to make best use of the opportunities provided by space missions such as the Radiation Belt Storm Probes, we determine the response of complementary subionospheric radiowave propagation measurements (VLF), riometer absorption measurements, cosmic noise absorption, and GPS-produced total electron content (vTEC) to different energetic electron precipitation (EEP). We model the relative sensitivity and responses of these instruments to idealized monoenergetic beams of precipitating electrons, and more realistic EEP spectra chosen to represent radiation belts and substorm precipitation. In the monoenergetic beam case, we find riometers are more sensitive to the same EEP event occurring during the day than during the night, while subionospheric VLF shows the opposite relationship, and the change in vTEC is independent. In general, the subionospheric VLF measurements are much more sensitive than the other two techniques for EEP over 200 keV, responding to flux magnitudes two-three orders of magnitude smaller than detectable by a riometer. Detectable TEC changes only occur for extreme monoenergetic fluxes. For the radiation belt EEP case, clearly detectable subionospheric VLF responses are produced by daytime fluxes that are ∼10 times lower than required for riometers, while nighttime fluxes can be 10,000 times lower. Riometers are likely to respond only to radiation belt fluxes during the largest EEP events and vTEC is unlikely to be significantly disturbed by radiation belt EEP. For the substorm EEP case both the riometer absorption and the subionospheric VLF technique respond significantly, as does the change in vTEC, which is likely to be detectable at ∼3–4 total electron content units.

Using ground-based subionospheric radio wave propagation data from two very low frequency (VLF) receiver sites, riometer absorption data, and THEMIS satellite observations, we examine in detail energetic electron precipitation (EEP) characteristics associated with two substorm precipitation events that occurred on 28 May 2010. In an advance on the analysis undertaken by Clilverd et al. (2008), we use phase observations of VLF radio wave signals to describe substorm-driven EEP characteristics more accurately than before. Using a >30 keV electron precipitation flux of 5.6 X 107 el. cm-2 sr-1 s-1 and a spectral gradient consistent with that observed by THEMIS, it was possible to accurately reproduce the peak observed riometer absorption at Macquarie Island (L = 5.4) and the associated NWC radio wave phase change observed at Casey, Antarctica, during the second, larger substorm. The flux levels were near to 80% of the peak fluxes observed in a similar substorm as studied by Clilverd et al. (2008). During the initial stages of the second substorm, a latitude region of 5 < L < 9 was affected by electron precipitation. Both substorms showed expansion of the precipitation region to 4 < L < 12 more than 30 min after the injection. While both substorms occurred at similar local times, with electron precipitation injections into approximately the same geographical region, the second expanded in an eastward longitude more slowly, suggesting the involvement of lower-energy electron precipitation. Each substorm region expanded westward at a rate slower than that exhibited eastward. This study shows that it is possible to successfully combine these multi-instrument observations to investigate the characteristics of substorms.

Utilizing observations from the Medium Energy Proton and Electron Detector (MEPED) on board the Polar Orbiting Environmental Satellite (POES) and the Microwave Limb Sounder (MLS) on board the Aura satellite, we demonstrate that there is a strong link between 100-300 keV loss cone electron count rates observed in the outer radiation belt and nighttime OH concentrations in the middle mesosphere at 71-78 km altitude. In theory, this can be expected because the ionization caused by energetic electron precipitation (EEP) leads to odd hydrogen (HOx) production through ionic reactions. However, this is the first time that OH production due to EEP has been observed. We consider daily mean data from 2 months, March 2005 and April 2006, which were selected because of (1) relatively high count rates of radiation belt electrons observed and (2) the absence of solar proton events that could mask the EEP effects. The results show that at 55-65 degrees magnetic latitude (equivalent to McIlwain L shells 3.0-5.6) increases in electron count rates by 2 orders of magnitude are accompanied by increases in nighttime OH concentration of 100%. There is a high correlation between MEPED and MLS data such that 56-87% of the OH variation can be explained by changes in EEP. Because the relation between MEPED count rate observations and the flux of electrons actually entering the atmosphere is not trivial, we discuss the possibility of using OH observations to obtain an estimate of EEP forcing that could be used in atmospheric modeling.

Utilizing observations from the Medium Energy Proton and Electron Detector (MEPED) on board the Polar Orbiting Environmental Satellite (POES) and the Microwave Limb Sounder (MLS) on board the Aura satellite, we demonstrate that there is a strong link between 100-300 keV loss cone electron count rates observed in the outer radiation belt and nighttime OH concentrations in the middle mesosphere at 71-78 km altitude. In theory, this can be expected because the ionization caused by energetic electron precipitation (EEP) leads to odd hydrogen (HOx) production through ionic reactions. However, this is the first time that OH production due to EEP has been observed. We consider daily mean data from 2 months, March 2005 and April 2006, which were selected because of (1) relatively high count rates of radiation belt electrons observed and (2) the absence of solar proton events that could mask the EEP effects. The results show that at 55-65 degrees magnetic latitude (equivalent to McIlwain L shells 3.0-5.6) increases in electron count rates by 2 orders of magnitude are accompanied by increases in nighttime OH concentration of 100%. There is a high correlation between MEPED and MLS data such that 56-87% of the OH variation can be explained by changes in EEP. Because the relation between MEPED count rate observations and the flux of electrons actually entering the atmosphere is not trivial, we discuss the possibility of using OH observations to obtain an estimate of EEP forcing that could be used in atmospheric modeling.