Issue 
J. Space Weather Space Clim.
Volume 9, 2019



Article Number  A43  
Number of page(s)  19  
DOI  https://doi.org/10.1051/swsc/2019041  
Published online  10 December 2019 
Research Article
The GOESR EUVS model for EUV irradiance variability
^{1}
Laboratory for Atmospheric and Space Physics, University of Colorado, Boulder, CO 80303, USA
^{2}
Cooperative Institute for Research in Environmental Sciences, University of Colorado, Boulder, CO 80309, USA
^{3}
National Centers for Environmental Information, National Oceanic and Atmospheric Administration, Boulder, CO 80305, USA
^{4}
Space Weather Prediction Center, National Oceanic and Atmospheric Administration, Boulder, CO 80305, USA
^{*} Corresponding author: thiemann@lasp.colorado.edu
Received:
3
May
2019
Accepted:
20
November
2019
The Geostationary Operational Environmental Satellite R (GOESR) series of four satellites are the next generation NOAA GOES satellites. Once on orbit and commissioned, they are renamed GOES 16–19, making critical terrestrial and space weather measurements through 2035. GOES 16 and 17 are currently on orbit, having been launched in 2016 and 2018, respectively. The GOESR satellites include the Extreme Ultraviolet (EUV) and Xray Irradiance Sensors (EXIS) instrument suite, which measures calibrated solar irradiance in eight lines or bands between 25 nm and 285 nm with the Extreme Ultraviolet Sensors (EUVS) instrument. EXIS also includes the XRay Sensor (XRS) instrument, which measures solar soft Xray irradiance at the legacy GOES bands. The EUVS Measurements are used as inputs to the EUVS Model, a solar spectral irradiance model for space weather operations that predicts irradiance in twentytwo 5 nm wide intervals from 5 nm to 115 nm, and one 10 nm wide interval from 117 to 127 nm at 30 s cadence. Once fully operational, NOAA will distribute the EUVS Model irradiance with 1 min latency as a primary space weather data product, ushering in a new era of rapid dissemination and measurement continuity of EUV irradiance spectra. This paper describes the EUVS Model algorithms, data sources, calibration methods and associated uncertainties. Typical model (relative) uncertainties are less than ~5% for variability at timescales longer than 6 h, and are ~25% for solar flare induced variability. The absolute uncertainties, originating from the instruments used to calibrate the EUVS Model, are ~10%. Examples of model results are presented at both subdaily and multiyear timescales to demonstrate the model’s capabilities and limitations. Example solar flare irradiances are also modeled.
Key words: spectral irradiance / space weather / modeling / space environment / solar activity
© E.M.B. Thiemann et al., Published by EDP Sciences 2019
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1 Introduction
Solar Extreme Ultraviolet (EUV, 10–121 nm) irradiance is the primary energy input into the Earth’s upper atmosphere at low to mid latitudes and at all latitudes during geomagnetically quiet periods. By ionizing gases, EUV radiation creates the ionosphere and heats the thermosphere. The solar EUV irradiance varies significantly due to the evolving 11year solar cycle, 27day solar rotation period and transient solar flares; resulting in corresponding changes in the density, temperature and composition of the thermosphere and ionosphere. Additionally, both the solar spectrum and atmospheric gas absorption crosssections are highly structured at EUV wavelengths, influencing the altitudes at which solar EUV radiation is absorbed. EUVinduced changes in the thermosphere can extend to Low Earth Orbit (160–2000 km), directly impacting satellite drag, where a hotter, denser thermosphere exerts a greater drag force on the satellites orbiting within it (Jachhia, 1959). With regard to the ionosphere, EUVinduced changes modulate its index of refraction, impacting transionospheric communication and navigation signals (Davies, 1990). The response of both the ionosphere and the thermosphere to EUV changes are relatively fast, with the ionosphere responding near instantaneously and the thermosphere responding in 2–4 h (e.g. Mendillo et al., 1974; Qian et al., 2010). As such, near real time solar EUV spectral irradiance information is needed for accurate space weather forecasting. To meet this demand, the National Oceanic and Atmospheric Administration (NOAA) has included an operational (continuous, high time cadence, low latency) EUV irradiance data product as part of its Geostationary Operational Environmental Satellite R (GOESR) series program. The GOESR satellites are scheduled to make observations from 2016 to 2035, providing nearly two decades of realtime continuous solar EUV irradiance measurements, and changing the paradigm for the availability and dissemination of spectral EUV irradiance data.
The NOAA GOESR series program is the latest iteration of the GOES satellite constellation, and consists of four satellites, each carrying a suite of identical instruments designed to monitor terrestrial and space weather continuously. Once commissioned, the GOESR series satellites will become GOES 16–19, with the first in the series, GOES16, having launched on 19 November 2016, and the last in the series, GOES19, anticipated to launch in the early 2020s and to be operational until at least 2035. The GOES satellites have monitored solar soft Xray irradiance in two bands since their inception in 1975 with the XRay Sensor (XRS) instruments, and GOES 13–15 have measured EUV irradiance in several bands with the Extreme Ultraviolet Sensor (EUVS) instruments (Viereck et al., 2007). For the GOESR series satellites, new versions of the XRS and EUVS instruments have been built as part of the EUV and Xray Irradiance Sensors (EXIS) instrument suite. The EXIS instrument suites were built at the Laboratory for Atmospheric and Space Physics (LASP) at the University of Colorado at Boulder. Each EXIS suite consists of the EUVS (Eparvier et al., 2009), which measures solar emissions at eight lines or bands between 25 nm and 285 nm, and the XRS (Chamberlin et al., 2009), which measures soft Xray irradiance in the legacy 0.1–0.8 nm and 0.05–0.4 nm bands.
The solar atmosphere is comprised of four distinct regions, which are (from nearest to furthest from the surface): the photosphere, chromosphere, transition region and corona. They can be equivalently categorized according to temperature with the photosphere being the coolest (~5700 K) and the corona being the hottest (>1 MK). Solar radiation emitted from the different regions of the solar atmosphere tends to vary differently as a result of the different processes driving the dynamics within them. For example, bright plasma tends to be concentrated regionally in active region magnetic loops in the corona, whereas in the chromosphere, bright plasma is distributed more uniformly across a magnetic network (Antia et al., 2003).
The emission lines observed by EUVS were selected to span a broad range of emission formation temperatures in the solar atmosphere in order to capture a broad range of irradiance variability. Specifically, EUVSA measures the He II 25.6 nm, Fe XV 28.4 nm and the He II 30.4 nm lines; EUVSB measures the C III 117.5 nm, H I 121.6 nm (Lymanα), C II 133.5 nm and the Si IV/O IV (blended) 140.5 nm lines; and EUVSC measures the Mg II emission line coretowing ratio (Heath & Schlesinger, 1986, hereafter, the Mg II Index) around 280 nm. These eight fullycalibrated EUVS Measurements (hereafter, the EUVS Measurements) are used as inputs to the EUVS Level 1B solar irradiance model (hereafter, the EUVS Model), which predicts solar EUV spectral irradiance continuously at 30 s cadence and 1 min latency in twentytwo 5 nm wide intervals from 5–115 nm and a single 10 nm interval from 117 nm–127 nm.
Solar EUV and Far Ultraviolet (FUV, 122–200 nm) irradiance is completely absorbed in the Earth’s upper atmosphere, requiring it to be measured by sophisticated spacebased instrumentation that is prone to degradation. This measurement difficulty has resulted in extended periods of time when direct spectral irradiance measurements are unavailable, necessitating models of solar EUV and FUV irradiance to bridge the observational gaps. Additionally, because calibrations often differ from instrument to instrument, models of solar irradiance provide the capability to estimate irradiance over a long timespan, while avoiding discontinuities that occur at the edges of concatenated datasets due to calibration differences.
Hinteregger (1981) proposed the first widely adopted EUV spectral irradiance variability model, which was based on the Atmospheric ExplorerE EUV Spectrophotometer (AEE/EUVS) measurements. In his paper, Hinteregger reported two models. The first model used the Fe XVI 33.5 nm and H Lymanβ 102.6 nm emissions to model coronal and chromospheric variability, respectively, with a set of regression coefficients that related these two emissions to the EUV irradiance spectrum. Fe XVI 33.5 nm and H Lymanβ 102.6 nm were only available from AEE/EUVS during this epoch and therefore this model was of limited use, but the success of the method laid the foundation for models that followed, including the GOES EUVS Model described here. The second model proposed by Hinteregger (1981) used regression coefficients relating the AEE/EUVS Measurements and the 10.7 cm solar radio flux (F10.7) daily and 81day mean values to estimate the spectrum. Later, Torr & Torr (1985) repartitioned the Hinteregger (1981) model spectra into 37 spectral intervals based on the absorption crosssections of major constituent gases in the Earth’s atmosphere.This repartitioning used twenty 5 nm intervals from 5 nm to 105 nm and then 17 narrow intervals to model the irradiance of individual emissions lines near regions where atmospheric absorption crosssections are highly structured. A significant advance occurred when Richards et al. (1994) introduced the EUV flux model for Aeronomic Calculations (EUVAC), which was essentially a recalibration of the Torr & Torr (1985) model coefficients to be consistent with more accurate rocket measurements. Later work by Tobiska & Eparvier (1998) and Tobiska et al. (2000) added spaceborne model inputs in addition to F10.7; and incorporated new suborbital rocket measurements and data from the San Marco 5 satellite into the model calibration dataset. Regular subdaily EUV measurements made by the Solar EUV Experiment (SEE) onboard the Thermosphere Ionosphere Mesosphere Energetics and Dynamics (TIMED) satellite were incorporated into the Flare Irradiance Spectral Model (FISM) developed by Chamberlin et al. (2007). Chamberlin et al. (2008) incorporated the capability for estimating solar flare irradiance based on the GOES/XRS irradiance into FISM. FISM was recently updated by Thiemann et al. (2017a) to use measurements from the EUV Monitor (EUVM) onboard the Mars Atmosphere and EvolutioN (MAVEN) probe as inputs and included spectral irradiance data from the EUV Variability Experiment (EVE) onboard the Solar Dynamics Observatory (SDO) in the model calibration dataset. The Solar Spectral Proxy Irradiance from GOES (SSPRING, Suess et al., 2016) model uses irradiance measurements from GOES15 as inputs. The EUVS Model presented here is an iteration of these preceding empirical models for solar spectral irradiance, and the first to be implemented in an operational environment.
The long term dataset of EUV spectral irradiance observed by TIMED/SEE provided an opportunity to rigorously test the idea that solar EUV variability can be decomposed into groups of emissions which vary similarly depending on the layer of the solar atmosphere from which they originate. A number of studies investigated these assumptions statistically (Kretzschmar et al., 2006; Chamberlin et al., 2007; Lilensten et al., 2007; Amblard et al., 2008; Dudok de Wit et al., 2008; Dudok de Wit et al., 2009; Cessateur et al., 2011). Notably, Amblard et al. (2008) used statistical methods to deconstruct the solar spectrum into elementary components and found that a minimum of three elementary spectra are needed to reconstruct the EUV spectral irradiance, representative of the inactive Sun, cool chromosphere and hot corona. Chamberlin et al. (2007) used the longterm TIMED/SEE dataset to determine model error for the integrated 0.1–193 nm band is reduced by ~30% when spaceborne modelinputs are used in lieu of F10.7. Later, Cessateur et al. (2011) used statistical methods to analyze the TIMED/SEE dataset and determined that empirical model error could be reduced by a factor of two using four spectral passbands measured from space rather than the commonly used F10.7 and Mg II indices.
This paper presents the EUVS Model, its coefficients and the methods used to derive them, as well as the model uncertainty. In Section 2, the model is described and the model equations are presented. Section 3 presents the historical data used to compute the model coefficients, the crosscalibration between historical and EXIS measurements, and the associated systematic uncertainty. Sections 4 and 5 present the long and shortterm model components and their uncertainty; and Section 6 discusses the overall results.
The primary purpose of this paper is to show how the EUVS Model coefficients and uncertainties are derived, rather than serve as a reference for their exact values. The coefficients and uncertainties presented here are for GOES16 at the time of this writing and may change with updates to the GOES16/EXIS calibration. Additionally, the coefficients and uncertainties for the GOES 17–19 EUVS Models will possibly differ from those presented here. As such, current values for the EUVS Model coefficients and uncertainties will be made available via the web address listed in the Acknowledgements section.
2 Model description
The EUVS Model estimates irradiance in the nth wavelength interval, E _{ n } (t), by solving the equation,
(1)where the summations are over the eight EUVS Measurements, E _{ n,0} is the offset for the nth wavelength interval, and P _{ i }(t) and Q _{ i }(t) are the long and shortterm components of the ith EUVS Measurement, and j _{ i,n } and k _{ i,n } are the long and shortterm regression coefficients, respectively. Note, the EXIS XRS measurements could also be included, but are not used in practice because their inclusion increases model uncertainty due to the variability at soft Xray wavelengths typically being much larger than that at EUV wavelengths. The long and shortterm components are defined by
(3)where X _{ i } (t) is the ith EUVS Measurement with reference offset, X _{ i,0}, and (lagging) moving average, , over time, T. The fractional units of P _{ i }(t) and Q _{ i }(t) reduces the sensitivity to subtle differences between spectral resolution and preflight calibrations of the historical model training data and the EUVS Measurements.
The shortterm component is for modeling rapid solar flare induced variability whereas the longterm component is for modeling more gradual nonflaring variability. The separation of flaring and nonflaring variability is done because different emission lines may contribute to the irradiance in a given spectral band during flares as a result of the significant temperature difference between flare and nonflare plasma. For example, the EUVS Measurement at 25.6 nm is primarily due to He II emissions during nonflaring times but dominated by Fe XIV emission during flares. As such, T is taken to be 6 h in order to be significantly longer than most solar flares. Specifically, T is ~4 times larger than the expected Xclass flare duration as determined by Veronig et al. (2002), who analyzed nearly 50,000 flares and found 90% of Xclass flares last less than 98 min. Modeled irradiance for long duration flares, which can last for several hours, will be subject to increased error if the flare irradiance contributes measurably to the 6hour average.
The offsets and coefficients in equations (1)–(3) are determined from historical data as follows: X _{ i,0} is chosen as the minimum value for the historical data representing the ith EUVS Measurement. E _{ n,0} and j _{ i,n } are found by multiple linear regression fitting of historical data representing P _{ i }(t) and the longterm average (nominally 6 h) of E _{ n }(t), ; and k _{ i,n } is found by multiple linear regression fitting of historical data representing Q _{ i }(t) and (). These methods are described in further detail in Sections 4 and 5, and the historical data are described in Section 3.
The EUVS Model wavelength intervals are 5 nm wide from 5 nm to 115 nm (i.e. 5–10 nm, 10–15 nm, …, 110–115 nm), with one 10 nm wide interval from 117 nm to 127 nm. This is the official operational wavelength scheme produced in real time by NOAA and is the focus of this paper. The cadence of X _{ i }(t) (and, hence, the EUVS Model) is 30 s, where X _{ i }(t) is a 30second centered running average of the ith EUVS Measurement. The EUVS Measurements all nominally have a 1second cadence except for EUVSC, which has a 3second cadence. is a running average of the preceding 6 h of X _{ i }(t). The EUVS Model latency is 1 min (i.e. EUVS Model irradiances are publicly available 1 min after the corresponding EUVS Measurements are made).
The operational nature of the EUVS Model requires the continuous production of irradiances, even during times when some EUVS Measurements are unavailable. This requires the derivation of independent sets of coefficients for the most likely contingencies. For brevity and clarity, we only focus on the nominal model, which assumes all EUVS Measurements are available. The same methods described here are used to derive coefficients and model uncertainty for contingency cases, which will be reported on the web link given in the Acknowledgements section. The model uncertainty will generally be higher if the number of inputs is reduced, in particular, if there is no substitute measurement with a similar formation temperature available as discussed in Section 1. It is important to note that NOAA plans to maintain redundant satellites on orbit, reducing the likelihood that the EUVS Model will be operational at reduced performance for an extended period of time.
3 Model training data and crosscalibration
3.1 Historical datasets
The historical training datasets used to calculate the EUVS Model coefficients, and the time and spectral ranges over which they span are shown in Figure 1, where Figure 1a shows the spectral irradiance datasets used to create the simulated spectral intervals (E _{ n }) and Figure 1b shows the datasets used for the simulated EUVS Measurements (X _{ i }). Two separate time intervals of data are used to generate model coefficients, corresponding to periods when either TIMED/SEE or SDO/EVE are well calibrated. These intervals are distinguished in Figure 1 with red and blue shading. Note, the TIMED/SEE dataset is well calibrated through 2011, but the TIMED/SEE (redshaded) interval stops at January 1, 2008 to prevent any sampling bias associated with including a disproportionate amount of data near solar minimum in the training dataset. Also shown in Figure 1b, with gray shading, is the time period and measurements used to crosscalibrate the historical training datasets with GOES16. The corresponding F10.7 flux is superimposed in the background of both panels for context. Additionally, Figure 1 indicates cases where a dataset is only used for either the short or longterm regression coefficients. For example, in the TIMEDEGS block in Figure 1a the designation “j _{ i,n }” implies this dataset is only used to find the longterm regression coefficient, j _{ i,n }. Table 1 provides further information on the historical datasets, including data levels, versions, sampling and accuracy for the datasets shown in Figure 1. The Purpose column in Table 1 indicates whether the dataset is used to find the shortterm or longterm model regression coefficients, or for crosscalibrating the EXIS measurements with the historical datasets.
Fig. 1 The time periods and wavelength ranges spanned by (a) spectral irradiance and (b) EXIS measurement historical datasets used to determine EUVS Model coefficients. Wavelengths or bands spanned are on the vertical axis and time is on the horizontal axis. Table 1 provides details on the datasets used. 
Fig. 2 Scatterplots and bestfits for the EUVSA line irradiances from FISM (a proxy for historical EVE data) and GOES16/EXIS on the horizontal and vertical axes, respectively. FISM is used inlieu of EVE/MEGSA because EVE/MEGSA and GOES16/EXIS did not make measurements contemporaneously. Daily average data from 1 February 2017 through 10 April 2018 are used. 
Fig. 3 Scatterplots and bestfits of the GOES16 EUVSB line irradiances and SORCE/SOLSTICE data. Daily average data from 1 February 2017 through 10 April 2018 are used. 
Fig. 4 Scatterplot and bestfit of the GOES16 EUVSC Mg II Index and Bremen Composite Mg II Index. Daily average data from 1 February 2017 through 10 April 2018 are used. The Bremen Composite Mg II Index source for this time period is from the GOME2B instrument onboard the MetOpB satellite. 
Fig. 5 Sample modelmeasurement comparisons for intervals derived from SDO/EVE data. 
Fig. 6 Sample modelmeasurement comparisons for intervals derived from TIMED/EGS data. 
Fig. 7 Histograms of all M class flares observed by MEGS A (dotted), flares used to to compute model coefficients in the MEGSA (dashed) and MEGSB (solid) wavelength ranges. 
Fig. 8 Scatterplots of measured and modeled peak flare irradiances for the flares used to calibrate the EUVS Model shortterm component. 
Fig. 9 A comparison of measured (black) and modeled (red) shortterm variability for a 24 h period beginning on 5 July 2012 12:00 UT. MeasurementModel differences are shown in gray. Five intervals are shown that were selected because they typically show relatively large enhancements during flares; the interval wavelength range is given on each panel. 
Historical training datasets used to determine the EUVS Model coefficients. The Purpose column indicates if the data are used to compute the longterm and/or shortterm model coefficients, or to crosscalibrate the historical data with the GOES16 EUVS Measurements.
3.2 EXIS crosscalibration with historical datasets
The historical data used to represent the EXIS measurements (shown in Fig. 1b) are recalibrated to best represent the GOES16/EXIS measurements prior to finding the GOES model coefficients as follows: for EUVSA and EUVSB, the historical spectral irradiance data are interpolated to match the pixel scales of the GOES16/EXIS detectors and then integrated over the same pixel masks used to compute the GOES16/EXIS Level 1b line irradiances. These values are then crosscalibrated against the GOES16/EXIS Level 1b line irradiances measured for the first 9 months of its mission using a first order Total Least Squares fit, which is implemented using the method of Van Huffel (1989). For EUVSA, measurements from the EVE/MEGSA channel do not overlap in time with GOES16/EXIS. In lieu of direct measurements for crosscalibrating EVE/MEGSA with GOES16, a version of the FISM model that is calibrated to EVE (Chamberlin et al., 2018) is used as an intermediary. For EUVSC, the Bremen Composite Mg II Index is fit against the GOES16 Mg II Index native scale. The Bremen Composite Mg II Index is sourced from the GOME2B instrument onboard the MetOpB satellite for the crosscalibration time period.
Figures 2–4 show the crosscalibration data and fits for EUVSA, B and C, respectively. The fit coefficients are given in Table 2 along with the crosscalibration uncertainty, σ _{CC,i }. Note that σ _{CC,i } only captures the uncertainty during the period of overlap, and uncorrected calibration drifts occurring in the training dataset, or future uncorrected drifts in the EUVS Measurements could increase the systematic error. These fits are used to simulate the EUVS Measurements using the historical training data. For all channels except EUVSA, σ _{CC,i } is the uncertainty of the linear fit. For EUVSA lines used in the model intervals from 5 nm to 35 nm, σ _{CC,i } is the quadrature sum of the crosscalibration uncertainty and the FISM model uncertainty, where the latter is reported in Thiemann et al. (2017a). For EUVSA lines used in the model intervals from 35 nm to 127 nm, an additional error term quantifying differences in the TIMED/SEE and SDO/EVE calibrations is included in the quadrature sum (in addition to the two aforementioned terms); this additional term is the uncertainty of the fit of the EUVSA line bandpass measured by SDO/EVE to that measured by TIMED/SEE. Therefore, two values are reported for σ _{CC,i } for the EUVSA lines in Table 2; the value in parenthesis corresponds with model intervals from 35 nm to 127 nm whereas the other corresponds with model intervals from 5 nm to 35 nm. There are also two values reported for the EUVSB 121.6 nm line; the value in parenthesis is for the GOES15/EUVS E channel crosscalibration, whereas the other value corresponds with SORCE/SOLSTICE.
Crosscalibration fit coefficients and uncertainty for the eight EUVS Measurements. The first two columns indicate the EUVSchannel and linecenter wavelength. The slope and offset correspond with the line equation coefficients to estimate the EUVS Measurements from the historical data shown in Figure 1, and the corresponding uncertainty is shown in the rightmost column.
3.3 Model uncertainties
The available historical training datasets are partitioned roughly in half, with onehalf used to fit the model coefficients and the other half used to compute the uncertainties. The uncertainty, σ, of a linear fit, say, , between two linearly related, Nvalued parameters y _{ l } and x _{ i } with fit coefficients, m and b is given by (e.g. Taylor, 1997),
Equation (4) is used in this study to compute the absolute uncertainties, with corresponding with the lth model output and y _{ l } corresponding with the lth observation. In order to account for the crosscalibration uncertainty (σ _{cc,i }, listed in Table 2), gaussian noise with a standard deviation of σ _{cc,i } is added to the ith model input prior to being used to compute .
It is useful to compare uncertainties in percent units. The percent uncertainty can be approximated as:
An additional useful metric is the percent uncertainty relative to the typical variability as represented by the standard deviation,
(6)where s is the standard deviation of the observations. It is important to keep in mind that σ _{%} and σ _{ s } are relative terms and, therefore, vary depending on the context in which they are computed. σ _{%} will be larger (smaller) for lower (higher) irradiance values. Similarly, σ _{ s } will be larger (smaller) when computed over a weak (strong) solar cycle. s and are computed over the time periods spanned by the historical data as shown in Figure 1.
In this study, the Bootstrap method (Efron, 1979) is used to estimate the model uncertainty. Onethousand synthetic datasets are generated by resampling (with replacement) the available pairs of values, and values for σ are computed from each synthetic dataset. It is found that the standard deviation of σ is typically ~3% (i.e. a reported 10% uncertainty has a spread of ~0.3%) for the longterm model component and ~15% for the short term model component. σ values and their standard deviations as determined by Bootstrapping are reported in Sections 4 and 5.
4 6hour average (longterm) model component
4.1 Longterm model formulation
The first two terms on the righthand side of equation (1) comprise the longterm model component. The historical spectral irradiance data are interpolated to 1nm sampling for determining the model coefficients, E _{ n,0} and j _{ i,n }. The regression coefficients are found using the Interactive Data Language (IDL) multiple linear regression program, regress.pro, via an iterative scheme: regressors with negative coefficients or that contribute to less than 5% of the fit (as a result of having relatively small coefficients) are set to zero, and the fits are recomputed until all coefficients are positive and all regressors contribute to at least 5% of the fit. It was found that for particularly noisy intervals, multiple linear regression can result in model error that is larger than that from single linear regression using a single EUVS Measurement. To account for this, a multiple linear regression model using the combined eight EUVS Measurements and eight single linear regression models using the eight EUVS Measurements individually are computed. The corresponding model errors, defined as the quadrature sum of the mean and standard deviation of the modelmeasurement difference, are compared, and the model with the smallest error is selected for each 1nm interval.
The coefficients for the 1nm intervals corresponding with a given EUVS Model interval are then summed over the wavelength range defined by the interval. An example is given to demonstrate how the fit coefficients determined from the 1 nm wide intervals are combined to correspond with a 5 nm wide interval: the coefficients for the 5–10 nm interval, E _{5–10,0} and j _{i,5–10} are found from the five sets of 1nm sampling coefficients as follows:
The factor of rescales the coefficients for spectral irradiance from 1nm to 5nm intervals, with units of irradiance per nm. For the 117–127 nm interval, the coefficients are scaled by .
4.2 Longterm model component results
In order to more accurately characterize the model uncertainty, the historical datasets are partitioned approximately in half, with the half corresponding with lower solar activity being used to compute the model coefficients, while the half corresponding with higher solar activity is used for modelmeasurement comparison to quantify the model uncertainty. The longterm component model uncertainties, σ _{LTM}, for each EUVS Model interval are reported in Table 3 in absolute units, units of percent uncertainty and percent variability. Additionally, Table 3 reports the standard deviation of the uncertainties, s _{ σ }, in absolute units and the longterm component instrument percent uncertainties, σ _{LTI}, the latter of which correspond with the signalweighted average instrument accuracies of the historical spectralirradiance measurements that are used to train the EUVS Model.
Longterm model uncertainty in absolute units along with its standard deviation. Uncertainties are also reported in percent units, both relative to the mean and observed variability (as measured by the standard deviation).
σ _{LTM} is the random uncertainty (i.e. the classical precision) of the EUVS Model longterm component, whereas σ _{LTI} is the Model’s systematic uncertainty (i.e. the classical accuracy). In other words, σ _{LTI} characterizes the magnitude of a constant bias in the model estimates, while σ _{LTM} characterizes the error in the 6hour variability estimated by the model. It is important to note the model absolute accuracy strictly depends on the absolute accuracy of the historical training data, and for a number of wavelength intervals, the random uncertainty (σ _{LTM}) is smaller than the uncertainty of the systematic bias (σ _{LTI}).
The sensitivity of the model parameters to the model training period used is investigated by exchanging the model training and validation periods. When the model is trained using the dataset half corresponding with higher solar activity and validated with the less active half, σ _{LTM} is smaller for all intervals except three (the 50–55 nm, 55–60 nm and 75–80 nm intervals). The average difference in σ _{LTM} is 15%. Additionally, the model coefficients found from the two different training periods tend to be different, and the sets of optimal model inputs for a model interval also tend to differ. This suggests that model error can be reduced by using training data from periods that have solar activity levels comparable to that of the period for which the model is being applied.
Sample model measurement comparisons for six intervals are shown in Figures 5 and 6 for model intervals derived from SDO/EVE and TIMED/EGS data, respectively. The modeled and measured irradiances are shown in the lefthand columns with red and black curves, respectively. The righthand columns show time series of the model uncertainty in percent units. Here, the percent uncertainty is defined as σ _{LT} divided by the daily modeled irradiance multiplied by 100%. The training and validation time periods can be inferred from Figures 5 and 6: times when data are shown correspond with the validation dataset times; times when data are absent between 30 April 2010 and 26 May 2014 (8 February 2002 and 1 January 2008) correspond with the training dataset times for the SDO/EVE (TIMED/EGS) data interval.
5 Minuteaverage (shortterm) model component
5.1 Shortterm model formulation
The shortterm model component is the second term on the right hand side of equation (1). The model coefficients are found using multiple linear regression at EUVS Model resolution rather than 1 nm sampling to improve the SignaltoNoise Ratio (SNR) of short time scale variability observations. Because the largest variations in solar EUV variability at short timescales are due to solar flares, the regression fitting is done using the peak irradiance enhancements occurring during solar flares. Typical fits result in a small constant term, which is collected in E _{ n,0} in equation (1) along with the constant term from the longterm component. The SDO/MEGS spectral irradiances are used for the 6–105 nm wavelength range, except for the 35–40 nm interval where the TIMED/EGS spectral irradiances are used. TIMED/EGS spectral irradiances are also used for the 105–115 nm wavelength range, and the SORCE/SOLSTICE irradiance is used for the 117–127 nm wavelength interval, each of which will be discussed independently in the following paragraphs.
The EVE/MEGS dataset is the most extensive dataset measured to date of EUV variability at short timescales, having measured EUV variability at 0.1 Hz for thousands of flares. MEGS A made measurements near continuously, while MEGS B has a lower duty cycle for degradation mitigation, resulting in MEGS B having made fewer flare observations. For every Mclass or larger flare observed by MEGS A, the preflare background and flare peak irradiance values are identified for each 1 nm interval in the MEGS wavelength range. This is done by manually identifying the time range corresponding with the background and peak irradiances using the bright flaring 13.29 Fe XXIII line and then using an automated program to average the background and identify the peak values. This same process was used by Thiemann et al. (2018), and the reader is referred to that paper for further details.
A total of 249 (50) flares are identified in the MEGS A (B) wavelength range. Of these flares, the 100 flares with the largest enhancement at 13.5 nm are used to compute model coefficients in the MEGS A wavelength range, while 30 flares with the largest enhancement at 97.5 nm are used for finding model coefficients in the MEGS B wavelength range. The MEGS B sample size is limited to 30 because including more flares increases the model error across many intervals as a result of having a low SNR. However, the fit uncertainties approach an asymptotic value after a sample size of ~10, so 30 flares is a large enough sample size to accurately characterize the model error. 13.5 nm and 97.5 nm are used for bright flare identification because they both show significant enhancement during flares, and hence, have high SNR. The 13.5 nm interval is dominated by a coronal Fe XXIII line during flares, and hence has relatively small opacity during flares. Although, it is important to note that even hot coronal lines, typically considered to be optically thin, limb darken during flares (Thiemann et al., 2018). On the other hand, the 97.5 nm interval in the MEGS B wavelength range is dominated by C III and Lymanγ, the latter of which strongly limb darkens due to the high abundance of H in the Sun’s atmosphere. This inherently biases the sample set of flares in the MEGS B wavelength range to those having originated away from the limb. As an alternative, the 49.5 nm interval, which is dominated by the hot coronal and, hence, less opaque Si XII emission, was also considered for flare identification. However, it was found that the SNR of the training set greatly decreased for many intervals in the MEGS B range as a result of their being dominated by more optically thick transition region emissions. Therefore, it was decided to use the 97.5 nm interval in order to have higher model accuracy for flares located nearer to disk center, as they tend to be more geoeffective (Qian et al., 2010) at a cost of higher model uncertainty for flares located nearer to the limb. The full set of flares, which includes both diskcenter and limb flares, is used to compute the model uncertainty. As such, the reported uncertainties are consistent with the average flare, and likely overestimated (underestimated) for diskcenter (limb) flares. Figure 7 shows a histogram of the flare disklocation for all 249 flares considered as well as those used in the MEGS A and MEGS B wavelength ranges. Note, in both cases, flares at 80° or further from disk center were discarded to avoid using flares that may be partially occulted by the limb.
The MEGS dataset only contains the EUVSA lines, and there are no corresponding high timecadence measurements made for the EUVSB and EUVSC lines except for the 121.6 nm Lymanα line, which is measured by the GOES15/EUVSE instrument. As a result, only simulated EXIS measurements for the EUVSA lines and Lymanα are available for generating the model coefficients for the intervals ranging from 5 nm to 105 nm (except the 35–40 nm interval).
The model intervals ranging from 105 nm to 115 nm are derived from the TIMEDSEE dataset. The 35–40 nm interval is also derived from the TIMEDSEE dataset because the MEGSB flare measurements had very low SNR over this wavelength range. Flares of magnitude M5 or larger occurring between 30 October 2002 and 18 October 2008 during a SEE observation are used to derive the model coefficients (where the flare timing is taken to be that reported in the GOES events list returned by the pr_gev.pro SolarSoft routine). This results in a training dataset of 36 flares.
The 117–127 nm interval is derived from the SORCE/SOLSTICE dataset using all X1 or larger flares that occurred during SOLSTICE observations between mid2003 and 2011, resulting in 19 flares. Because SOLSTICE does not measure shortward of 117 nm, only simulated EUVSB lines are used to generate the model coefficients for this wavelength interval.
5.2 Shortterm model component results
Table 4 shows the corresponding shortterm model uncertainty (σ _{STM}) values in absolute units and their standard deviations (as determined by Bootstrapping). Values for for σ _{STM} are also given in percent units relative to the mean and standard deviation of the flare peak intensities of the training dataset. The relative magnitude of the mean flare peak values to the peak daily average values are given in the rightmost column in percent units for context. As is the case with the longterm model component, σ _{STM} is the random uncertainty, and the shortterm instrument uncertainty (σ _{STI}, not shown) is the systematic uncertainty of the shortterm model component. The values for σ _{STI} are similar to σ _{LTI}, but somewhat smaller from 35 nm to 105 nm due to the smaller absolute uncertainty of SDO/EVE relative to TIMED/SEE over these wavelengths.
Shortterm model uncertainty in absolute units along with its standard deviation. Uncertainties are also reported in percent units, both relative to the mean and observed flare peak intensity of the training dataset. The rightmost column shows the mean flare peak intensity relative to the daily average in percent units for context.
Figure 8 shows scatterplots of the measured and modeled peak flare enhancements from the training data for all 23 EUVS model intervals. From Figure 8 and Table 4, it is evident that the EUVS model predicted peak flare enhancements are highly correlated with measurements for most intervals.
In Figures 9 and 10, the shortterm EUVS Model predictions are compared with SDO/EVE measurements for two sample days, beginning at 12:00 UT on 5 July 2012 and 4 August 2011, respectively. Each figure shows measured and modeled shortterm irradiances for a 24hour period using black and red curves, respectively, for five wavelength intervals chosen because they tend to show a relatively large flare enhancement. The curves represent the shortterm solar variability, i.e. the rightmost term in equation (1). Measurementmodel differences are shown in each panel with gray curves. The SDO/EVE and simulated EXIS measurement data have the 6hour moving average removed and are plotted at 10 s cadence with a 30second (threesample) moving average applied. The missing SDO/EVE data in Panels e of both figures is because SDO/EVE reduces its duty cycle at this wavelength range for degradation mitigation.
Considering first the nonflaring shortterm variability that is most prevalent in Figure 9, the EUVS Model captures this variability well in the 5–10 nm and the 30–35 nm wavelength intervals, but to a lesser degree in the 10–15 nm and 15–20 nm intervals. The EUVS Model is expected to model variability in the 30–35 nm interval well because it is dominated by the 30.4 nm He II emission and this emission is used as a direct input to the EUVS Model. The emissions forming from 5 nm to 20 nm are formed primarily at hot temperatures in the solar corona. Therefore, their corresponding emission measures can be invisible to the EXIS measured emissions, which are formed at relatively cooler temperatures. The shortterm variability in these intervals is modeled predominantly using the 25.6 nm and 28.4 nm EXIS Measurements. The reason for the relatively large magnitude of the 25.6 nm measurement coefficient for these wavelength regions dominated by coronal forming emissions is because the 25.6 nm EXIS measurement is dominated by an Fe XXIV emission at 25.53 nm during solar flares. As such, the flare enhancements are well modeled in the 10–20 nm intervals, when the Fe XXIV emission dominates the EXIS measurement, but the model performance decreases during nonflaring periods when the cooler forming He II emission dominates.
Figure 8 quantifies how well the EUVS Model predicts the peak flare emission. From Figure 8, there appears to be some slight nonlinearity between the model predictions and measurements. For example, in the 5–10 nm panel, the slope is steeper below 10 μW/m^{2} than it is above 10 μW/m^{2}, resulting in the model underpredicting enhancements for (less frequent) larger flares.
Figures 9 and 10 provide insight into how well the EUVS Model captures the time evolution of solar flare emissions. Flare emissions at EUV wavelengths can pass through a number of phases as they evolve in time (Woods et al., 2011), including an initial impulsive phase followed by a gradual phase, which can be delayed in time depending of the emission formation temperature (Thiemann et al., 2017b); later flare phases typically associated with eruptive flares include a dimming phase and a late phase. Examples of all these phases except the late phase can be seen in Figures 9 and 10. Three impulsive flares are evident in Figure 9, occurring near 6 × 10^{4}, 8 × 10^{4} and 13 × 10^{4} s. These are all wellmodeled in the 30–35 nm interval, again due to it being driven by the 30.4 nm EXIS measurement. The impulsive phase for these flares is modeled poorly in the 15–20 nm interval, where it is underpredicted for the first two flares and overpredicted for the third flare. This is likely because this interval is driven predominantly by the Fe XXIV 25.53 nm emission during flares, which typically lacks an impulsive phase. The model does not adequately capture the gradual phase delay in the 10–20 nm intervals for the first two flares in Figure 9. In the 10–15 nm interval, the model overpredicts the gradual phase delay, whereas in the 15–20 nm interval the gradual phase delay is slightly under predicted, for both flares. Coronal dimming occurs after both flares in Figure 10 and is responsible for the emission decrease in panels c and d following one or both flares. The apparent dimming in panels a and b is a result of subtracting off the lagging 6hour average from the data and the relatively long duration of the two flares. The dimming is underpredicted in the 15–20 nm range for both flares because the Fe XXIV emission typically does not dim. For the second flare, dimming appears in the model but not in the data because this interval is driven by the 30.4 nm EXIS measurement, which showed dimming for this flare.
6 Discussion
All coefficients needed to compute EUVS Model irradiances from the EUVS Measurements are reported in Appendix. The EUVS model coefficients determined here from irradiance measurements made during Solar Cycles 23 and 24 are expected to apply to future solar cycles because the relationship between EUV emissions and the thermal structure of the solar atmosphere are expected to be invariant between solar cycles. Because the less (solar) active half of the available data were used to determine the model coefficients and the more active half of the data were used to characterize its uncertainty, the model performance has been validated for periods of higher activity than that occurring during the period used to compute the model coefficients. However, the available historical datasets span relatively moderate levels of solar activity when compared to previous solar cycles (e.g. Solar Cycles 20–22) and the accuracy of the model may decrease during periods when solar activity exceeds that of the historical datasets. The model performance during future stronger solar cycles can be improved by recomputing the model coefficients using the more active half of the historical data.
The EUVS Model uncertainties for the longterm component are comparable to those from the most recently published FISM update by Thiemann et al. (2017a). A key difference between the FISM and EUVS models is that FISM uses a single measurement as an input for each model interval and decomposes the daily average measurement and spectral irradiance predictions into solar cycle and solar rotation components, whereas the EUVS Model uses multiple measurements as inputs for each model interval and does not decompose them according to longerterm solar variability. This indicates that, at least at 5 nm sampling, daily average irradiance can be accurately predicted without decomposing solar variability into solar cycle and solar rotation components if the model inputs span a broad range of temperatures in the Sun’s atmosphere using model coefficients derived from multiple linear regression analysis. This has important implications for realtime solar spectral irradiance estimation, which does not have the future knowledge required when decomposing solar irradiance variability into solar cycle and solar rotation components.
On the other hand, the EUVS Model uncertainties for the shortterm component are a significant improvement over those from FISM. For example, Thiemann et al. (2017a) reported that the FISMM flare (i.e. shortterm) model uncertainty is 35% and 70% in the 13.5 nm and 30.5 nm (1nm wide) intervals, respectively. These values should be compared to the 10–15 nm and 30–35 nm interval uncertainties of 26.8% and 12.7%, respectively, reported in Table 4. Although some of the difference in uncertainties is due to the EUVS Model having larger intervals than FISM, by comparing the model measurement comparisons in Figure 8 with those for FISMM in Figure 4 of Thiemann et al. (2017a), it is clear that a significant part of the uncertainty differences are due to FISM using soft Xray irradiance variability, which varies by factors of 10–100, to estimate EUV irradiance variability, which varies by tens of percent. This is most pronounced in the intervals that contain the He II 30.4 nm emission line. The uncertainties of the EUVS Model and FISM for the intervals containing the Fe XXIII 13.3 nm emission line are similar, and likely due to Fe XXIII tending to evolve similarly to the 0.1–0.8 nm band (Thiemann et al., 2017b). In order to reduce model uncertainty, future empirical flare irradiance models should use EUV irradiance measurements as model inputs when possible. However, it is important to note that Thiemann et al. (2018) showed that the GOES XRS peak emission measure is linearly proportional to peak EUV irradiance for hotforming emission lines (T > 9 MK) during flares, unlike the XRS irradiances, which have a highly nonlinear relationship with EUV emissions as discussed above and reported in Thiemann et al. (2017a). As such, the XRS emission measure is another alternative for reducing flare irradiance model uncertainty.
The EUVS Model can be improved in a number of ways: an intermediate nonflaring timescale could be added to accurately capture variability that occurs at timescales of a few hours that is not the result of flares as seen, for example, in Figure 9. This could be a third term in equation (1) at an intermediate time scale, say 2 h. The current EUVS Model does not explicitly attempt to model the various EUV emission flare phases. Figures 9 and 10 suggest that the k _{ i } coefficients are weighted to capture some of the flare phases accurately. For example, model intervals that tend to show an impulsive flare phase are predominantly driven by the impulsive He II 30.4 nm model input, and those intervals that do not tend to be impulsive tend to be predominantly driven by the Fe XXIV 25.3 nm (nominally the He II EXIS 25.6 nm line measurement). Further, from Figure 9, the dimming phase appears to be accurately modeled for the two flares shown although it is important to note that a false dimming phase tends to appear in the 10–15 nm interval for long duration flares. The gradual phase delay is not properly modeled using the current EUVS Model algorithm. Thiemann et al. (2017b) showed that by applying a simple differential equation to emission line measurements from hotter EUV lines, the delay and broadening apparent in cooler EUV lines can be predicted if the time difference between the line peaks can be determined independently. Finally, this study did not consider the degree to which the EUVS Model can reproduce the EUV Late Phase. The Fe XV 28.4 nm EUVS Measurement does tend to show the EUV Late Phase and, as such, could be used to drive model intervals that also show the EUV Late Phase.
The same methods used to derive the coefficients and uncertainties for the relatively coarse wavelength intervals presented here can be applied to other more common wavelength schemes used in upper atmospheric research such as those introduced in Richards et al. (1994) and Solomon & Qian (2005). The broadband intervals of Richards et al. (1994) correspond directly with the broadband intervals presented here, and the remaining work involves finding coefficients for the 17 emission line intervals. The Solomon & Qian (2005) intervals do not match those presented here and also include wavelengths below 5 nm. As such, producing coefficients for the Solomon & Qian (2005) intervals would require resampling the source datasets and expanding the source data to include wavelengths below 5 nm.
The results presented in this paper are for the EXIS instrument on the GOES16 satellite. The reader should follow the link in the Acknowledgements section for possible updates, as well as values for other flight models as they become available. The EXIS dependent model uncertainties, σ _{CC,i}, are expected to be similar between EXIS models on GOES 17–19. As such, it is expected the EUVS Model uncertainties presented here for GOES16/EXIS will hold for the entire GOESR series EXIS instruments. The EUVS Model coefficients are expected to be similar between GOES 16–19 because the EXIS derived inputs, P _{ i }(t) and Q _{ i }(t), are independent of absolute calibration. However, P _{ i }(t) and Q _{ i }(t) for GOES 17–9 should be fit to those from GOES16 to confirm this, and k _{ i } and j _{ i } should be scaled accordingly if the corresponding slopes differ from unity.
Acknowledgments
The editor thanks two anonymous referees for their assistance in evaluating this paper.
This work was funded by the NOAA GOES Program contract to the University of Colorado, NNG07HW00C. The data used in this study are publicly available for download through the following locations: SDO/EVE data are available from http://lasp.colorado.edu/home/eve/data/dataaccess/. TIMED/SEE data are available from http://lasp.colorado.edu/home/see/data/. SORCE/SOLSTICE data are available from http://lasp.colorado.edu/home/sorce/data/. Mg II c/w data are available from http://www.iup.unibremen.de/gome/gomemgii.html.
Once the EXIS data have been approved for public release, EUVS Model spectra will be publicly available through the National Centers for Environmental Information website, https://www.ngdc.noaa.gov/stp/spaceweather.html. Additionally, future updates to the model coefficients and uncertainties for GOESR, as well as future releases of those parameters for GOESS through GOESU will be available through https://data.ngdc.noaa.gov/platforms/solarspaceobservingsatellites/goes/goesXX/l1b/exisl1bsfeu/model_coefficients/, where XX = 16–19.
Appendix
EUVS model coefficients
EUVS Measurement reference values.
Longterm model coefficients.
Shortterm model coefficients.
References
 Amblard PO, Moussaoui S, Dudok de Wit T, Aboudarham J, Kretzschmar M, Lilensten J, Auchère F. 2008. The EUV Sun as the superposition of elementary Suns. A&A 487(2): L13–L16. https://doi.org/10.1051/00046361:200809588. [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Antia H, Bhatnagar A, Ulmschneider P. 2003. Lectures on solar physics, vol 619, Springer Science & Business Media, Berlin, Germany, 619 p. https://doi.org/10.1007/3540369635. [Google Scholar]
 Cessateur G, Dudok de Wit T, Kretzschmar M, Lilensten J, Hochedez JF, Snow M. 2011. Monitoring the solar UV irradiance spectrum from the observation of a few passbands. A&A 528: A68. https://doi.org/10.1051/00046361/201015903. [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Chamberlin P, Woods T, Didkovsky L, Eparvier F, Jones A, et al. 2018. Solar ultraviolet irradiance observations of the solar flares during the intense September 2017 storm period. Space Weather 16(10): 1470–1487. https://doi.org/10.1029/2018SW001866. [CrossRef] [Google Scholar]
 Chamberlin PC, Woods TN, Eparvier FG. 2007. Flare irradiance spectral model (FISM): daily component algorithms and results. Space Weather 5(7). https://doi.org/10.1029/2007SW000316. [Google Scholar]
 Chamberlin PC, Woods TN, Eparvier FG. 2008. Flare irradiance spectral model (FISM): flare component algorithms and results. Space Weather 6(5). https://doi.org/10.1029/2007SW000372. [Google Scholar]
 Chamberlin PC, Woods TN, Eparvier FG, Jones AR. 2009. Next generation xray sensor (XRS) for the NOAA GOESR satellite series. In: Proc. SPIE 7438, Solar Physics and Space Weather Instrumentation III, 743802 (26 August 2009), San Diego, California, United States. https://doi.org/10.1117/12.826807. [Google Scholar]
 Davies K. 1990. Ionospheric radio. IEE electromagnetic waves series, vol. 31, Peter Peregrinus Ltd., London, UK. https://doi.org/10.1049/PBEW031E. [CrossRef] [Google Scholar]
 Dudok de Wit T, Kretzschmar M, Lilensten J, Woods T. 2009. Finding the best proxies for the solar UV irradiance. Geophys Res Lett 36(10). https://doi.org/10.1029/2009GL037825. [CrossRef] [Google Scholar]
 Dudok de Wit T, Kretzschmar M, Aboudarham J, Amblard PO, Auchère F, Lilensten J. 2008. Which solar EUV indices are best for reconstructing the solar EUV irradiance? Adv Space Res 42(5): 903–911. https://doi.org/10.1016/j.asr.2007.04.019. [CrossRef] [Google Scholar]
 Efron B. 1979. Bootstrap methods: another look at the jackknife annals of statistics. Ann Stat 7(1): 1–26. https://doi.org/10.1214/aos/1176344552. [Google Scholar]
 Eparvier FG, Crotser D, Jones AR, McClintock WE, Snow M, Woods TN. 2009. The extreme ultraviolet sensor (EUVS) for GOESR. In: Proc. SPIE 7438, Solar Physics and Space Weather Instrumentation III, 743804 (23 September 2009), San Diego, California, United States. https://doi.org/10.1117/12.826445. [Google Scholar]
 Heath DF, Schlesinger BM. 1986. The Mg 280nm doublet as a monitor of changes in solar ultraviolet irradiance. J Geophys Res: Atmos 91(D8): 8672–8682. https://doi.org/10.1029/JD091iD08p08672. [Google Scholar]
 Hinteregger H. 1981. Representations of solar EUV fluxes for aeronomical applications. Adv Space Res 1(12): 39–52. https://doi.org/10.1016/02731177(81)904166. [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Jachhia L. 1959. Two atmospheric effects in the orbital acceleration of artificial satellites. Nature 183(4660): 526–527. https://doi.org/10.1038/183526a0. [CrossRef] [Google Scholar]
 Kretzschmar M, Lilensten J, Aboudarham J. 2006. Retrieving the solar EUV spectral irradiance from the observation of 6 lines. Adv Space Res 37(2): 341–346. https://doi.org/10.1016/j.asr.2005.02.029. [NASA ADS] [CrossRef] [Google Scholar]
 Lilensten J, Dudok de Wit T, Amblard PO, Aboudarham J, Auchère F, Kretzschmar M. 2007. Recommendation for a set of solar EUV lines to be monitored for aeronomy applications. Ann Geophys 25(6): 1299–1310. https://doi.org/10.5194/angeo2512992007. [CrossRef] [Google Scholar]
 McClintock WE, Rottman GJ, Woods TN. 2005. SolarStellar Irradiance Comparison Experiment II (SOLSTICE II): Instrument concept and design. Sol Phys 230(1–2): 225–258. https://doi.org/10.1007/s112070057432x. [Google Scholar]
 Mendillo M, Klobuchar J, Fritz R, Da Rosa A, Kersley L, et al. 1974. Behavior of the ionospheric F region during the great solar flare of August 7, 1972. J Geophys Res 79(4): 665–672. https://doi.org/10.1029/JA079i004p00665. [CrossRef] [Google Scholar]
 Qian L, Burns AG, Chamberlin PC, Solomon SC. 2010. Flare location on the solar disk: modeling the thermosphere and ionosphere response. J Geophys Res: Space Phys 115(A9). https://doi.org/10.1029/2009JA015225. [Google Scholar]
 Richards P, Fennelly J, Torr D. 1994. EUVAC: a solar EUV flux model for aeronomic calculations. J Geophys Res: Space Phys 99(A5): 8981–8992. https://doi.org/10.1029/94JA00518. [NASA ADS] [CrossRef] [Google Scholar]
 Solomon SC, Qian L. 2005. Solar extremeultraviolet irradiance for general circulation models. J Geophys Res: Space Phys 110(A10). https://doi.org/10.1029/2005JA011160. [Google Scholar]
 Suess K, Snow M, Viereck R, Machol J. 2016. Solar Spectral Proxy Irradiance from GOES (SSPRING): a model for solar EUV irradiance. J Space Weather Space Clim 6: A10. https://doi.org/10.1051/swsc/2016003. [CrossRef] [Google Scholar]
 Taylor J. 1997. Introduction to error analysis, the study of uncertainties in physical measurements, University Science Books, Sausalito, CA, USA. [Google Scholar]
 Thiemann E, Chamberlin P, Eparvier F, Epp L. 2018. Centertolimb variability of hot coronal EUV emissions during solar flares. Sol Phys 293(2): 19. https://doi.org/10.1007/s1120701812442. [NASA ADS] [CrossRef] [Google Scholar]
 Thiemann EM, Chamberlin PC, Eparvier FG, Templeman B, Woods TN, Bougher SW, Jakosky BM. 2017a. The MAVEN EUVM model of solar spectral irradiance variability at Mars: algorithms and results. J Geophys Res: Space Phys 122(3): 2748–2767. https://doi.org/10.1002/2016JA023512. [Google Scholar]
 Thiemann EM, Eparvier FG, Woods TN. 2017b. A time dependent relation between EUV solar flare lightcurves from lines with differing formation temperatures. J Space Weather Space Clim 7: A36. https://doi.org/10.1051/swsc/2017037. [CrossRef] [EDP Sciences] [Google Scholar]
 Tobiska WK, Eparvier F. 1998. EUV97: improvements to EUV irradiance modeling in the soft Xrays and FUV. Sol Phys 177(1–2): 147–159. https://doi.org/10.1023/A:10049314. [NASA ADS] [CrossRef] [Google Scholar]
 Tobiska WK, Woods T, Eparvier F, Viereck R, Floyd L, Bouwer D, Rottman G, White O. 2000. The SOLAR2000 empirical solar irradiance model and forecast tool. J Atmos SolarTerr Phys 62(14): 1233–1250. https://doi.org/10.1016/S13646826(00)000705. [NASA ADS] [CrossRef] [Google Scholar]
 Torr MR, Torr D. 1985. Ionization frequencies for solar cycle 21: revised. J Geophys Res: Space Phys 90(A7): 6675–6678. https://doi.org/10.1029/JA090iA07p06675. [Google Scholar]
 Van Huffel S. 1989. The extended classical total least squares algorithm. J Comput Appl Math 25(1): 111–119. https://doi.org/10.1016/03770427(89)900800. [CrossRef] [Google Scholar]
 Veronig A, Temmer M, Hanslmeier A, Otruba W, Messerotti M. 2002. Temporal aspects and frequency distributions of solar soft Xray flares. A&A 382(3): 1070–1080. https://doi.org/10.1051/00046361:20011694. [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Viereck R, Hanser F, Wise J, Guha S, Jones A, McMullin D, Plunket S, Strickland D, Evans S. 2007. Solar extreme ultraviolet irradiance observations from GOES: design characteristics and initial performance. In: Proc. SPIE 6689, Solar Physics and Space Weather Instrumentation II, 66890K (20 September 2007), San Diego, California, United States. https://doi.org/10.1117/12.734886. [Google Scholar]
 Woods T, Eparvier F, Hock R, Jones A, Woodraska D, et al. 2010. Extreme Ultraviolet Variability Experiment (EVE) on the Solar Dynamics Observatory (SDO): overview of science objectives, instrument design, data products, and model developments. In: The solar dynamics observatory, Pesnell WD, Thompson BJ, Chamberlin PC (Eds.), Springer Science and Business Media, New York, NY, USA, pp. 115–143. https://doi.org/10.1007/s1120701198413. [CrossRef] [Google Scholar]
 Woods TN, Eparvier FG, Bailey SM, Chamberlin PC, Lean J, Rottman GJ, Solomon SC, Tobiska WK, Woodraska DL. 2005. Solar EUV Experiment (SEE): mission overview and first results. J Geophys Res: Space Phys 110(A1): A01312. https://doi.org/10.1029/2004JA010765. [CrossRef] [Google Scholar]
 Woods TN, Hock R, Eparvier F, Jones AR, Chamberlin PC, et al. 2011. New solar extremeultraviolet irradiance observations during flares. Astrophys J 739(2): 59. https://doi.org/10.1088/0004637X/739/2/59. [NASA ADS] [CrossRef] [Google Scholar]
 Woods TN, Tobiska WK, Rottman GJ, Worden JR. 2000. Improved solar Lyman α irradiance modeling from 1947 through 1999 based on UARS observations. J Geophys Res: Space Phys 105(A12): 27195–27215. https://doi.org/10.1029/2000JA000051. [Google Scholar]
Cite this article as: Thiemann EMB, Eparvier FG, Woodraska D, Chamberlin PC, Machol J, 2019. The GOESR EUVS model for EUV irradiance variability. J. Space Weather Space Clim. 9, A43.
All Tables
Historical training datasets used to determine the EUVS Model coefficients. The Purpose column indicates if the data are used to compute the longterm and/or shortterm model coefficients, or to crosscalibrate the historical data with the GOES16 EUVS Measurements.
Crosscalibration fit coefficients and uncertainty for the eight EUVS Measurements. The first two columns indicate the EUVSchannel and linecenter wavelength. The slope and offset correspond with the line equation coefficients to estimate the EUVS Measurements from the historical data shown in Figure 1, and the corresponding uncertainty is shown in the rightmost column.
Longterm model uncertainty in absolute units along with its standard deviation. Uncertainties are also reported in percent units, both relative to the mean and observed variability (as measured by the standard deviation).
Shortterm model uncertainty in absolute units along with its standard deviation. Uncertainties are also reported in percent units, both relative to the mean and observed flare peak intensity of the training dataset. The rightmost column shows the mean flare peak intensity relative to the daily average in percent units for context.
All Figures
Fig. 1 The time periods and wavelength ranges spanned by (a) spectral irradiance and (b) EXIS measurement historical datasets used to determine EUVS Model coefficients. Wavelengths or bands spanned are on the vertical axis and time is on the horizontal axis. Table 1 provides details on the datasets used. 

In the text 
Fig. 2 Scatterplots and bestfits for the EUVSA line irradiances from FISM (a proxy for historical EVE data) and GOES16/EXIS on the horizontal and vertical axes, respectively. FISM is used inlieu of EVE/MEGSA because EVE/MEGSA and GOES16/EXIS did not make measurements contemporaneously. Daily average data from 1 February 2017 through 10 April 2018 are used. 

In the text 
Fig. 3 Scatterplots and bestfits of the GOES16 EUVSB line irradiances and SORCE/SOLSTICE data. Daily average data from 1 February 2017 through 10 April 2018 are used. 

In the text 
Fig. 4 Scatterplot and bestfit of the GOES16 EUVSC Mg II Index and Bremen Composite Mg II Index. Daily average data from 1 February 2017 through 10 April 2018 are used. The Bremen Composite Mg II Index source for this time period is from the GOME2B instrument onboard the MetOpB satellite. 

In the text 
Fig. 5 Sample modelmeasurement comparisons for intervals derived from SDO/EVE data. 

In the text 
Fig. 6 Sample modelmeasurement comparisons for intervals derived from TIMED/EGS data. 

In the text 
Fig. 7 Histograms of all M class flares observed by MEGS A (dotted), flares used to to compute model coefficients in the MEGSA (dashed) and MEGSB (solid) wavelength ranges. 

In the text 
Fig. 8 Scatterplots of measured and modeled peak flare irradiances for the flares used to calibrate the EUVS Model shortterm component. 

In the text 
Fig. 9 A comparison of measured (black) and modeled (red) shortterm variability for a 24 h period beginning on 5 July 2012 12:00 UT. MeasurementModel differences are shown in gray. Five intervals are shown that were selected because they typically show relatively large enhancements during flares; the interval wavelength range is given on each panel. 

In the text 
Fig. 10 Same as Figure 9, but for a period beginning on 4 August 2011 12:00 UT. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.