the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
A calibration procedure which accounts for nonlinearity in singlemonochromator Brewer ozone spectrophotometer measurements
Zahra Vaziri Zanjani
Omid Moeini
Tom McElroy
David Barton
Vladimir Savastiouk
It is now known that singlemonochromator Brewer spectrophotometer ozone and sulfur dioxide measurements suffer from nonlinearity at large ozone slant column amounts due to the presence of instrumental stray light caused by scattering within the optics of the instrument. Because of the large gradient in the ozone absorption spectrum in the nearultraviolet, the atmospheric spectra measured by the instrument possess a very large gradient in intensity in the 300 to 325 nm wavelength region. This results in a significant sensitivity to stray light when there is more than 1000 Dobson units (DU) of ozone in the light path. As the light path (air mass) through ozone increases, the straylight effect on the measurements also increases. The measurements can be of the order of 10 %, low for an ozone column of 600 DU and an air mass factor of 3 (1800 DU slant column amount), which is an example of conditions that produce large slant column amounts.
Primary calibrations for the Brewer instrument are carried out at Mauna Loa Observatory in Hawaii and Izana Observatory in Tenerife. They are done using the Langley plot method to extrapolate a set of measurements made under a constant ozone vertical column to an extraterrestrial calibration constant. Since the effects of a small nonlinearity at moderate ozone paths may still be important, a better calibration procedure should account for the nonlinearity of the instrument response. Studies involving the scanning of a laser source have been used to characterize the straylight response of the Brewer (Fioletov et al., 2000), but until recently these data have not been used to elucidate the relationship between the straylight response and the ozone measurement nonlinearity.
In a study done by Karppinen et al. (2015), a method for correcting stray light has been presented that uses an additive correction, which is determined via instrument slit characterization and a radiative transfer model simulation and is then applied to the single Brewer data (Karppinen et al., 2015). The European Brewer Network is also applying straylight corrections, which includes an iterative process that results in correcting the single Brewer data to agree with double Brewer data (Rimmer et al., 2018; Redondas et al., 2018). The first model requires measurements of the slit function and the latter method relies on a calibrated instrument, such as a double Brewer, to characterize the instrument and to determine a correction for stray light.
This paper presents a simple and practical method to correct for the effects of stray light, which includes a mathematical model of the instrument response and a nonlinear retrieval approach that calculates the best values for the model parameters. The model can then be used in reverse to provide more accurate ozone values up to a defined maximum ozone slant path. The parameterization used was validated using an instrument physical model simulation. This model can be applied independently to any Brewer instrument and correct for the effects of stray light.
The Brewer spectrophotometer instrument is a diffractiongrating polychromatic spectrophotometer which produces monochromatic light at a set of six exit slits. It measures total column amounts of ozone, sulfur dioxide and aerosol optical thickness (Silva and Kirchhoff, 2004) by direct sun measurements at five wave bands centred at the approximate wavelengths of 306.3, 310.1, 313.5, 316.8 and 320.0 nm (Kerr et al., 1984). The Brewer spectrophotometer was designed as a replacement for the Dobson instrument in the 1970s and became commercially available in the 1980s (Redondas et al., 2014).
There are two types of Brewers in use today: single Brewers (Mk II, Mk IV, Mk V) and double Brewers (Mk III). The double Brewer (DBr) spectrophotometer is the combination of two single Brewer (SBr) optical frames wherein the exit slits of the first monochromator are the entrance slits of the second one. The first monochromator disperses the light, while the second one recombines it, giving the instrument a much greater capability for straylight rejection (Gröbner et al., 1998).
As in all optical instruments, stray light is also problematic in Brewer spectrophotometers. The source of stray light in a monochromator is light scattering from the various optical surfaces and walls of the instrument. Light scattered from the atmosphere is also a source of stray light, which is called the skyscattered radiation. Skyscattered radiation can be eliminated by taking direct sun measurements very close to the sun and subtracting them from the focused sun measurements (Josefsson, 1992). Instrumental stray light becomes much more significant when the source of energy or detector sensitivity changes rapidly as a function of wavelength. In the Brewer, the source of energy is the sun or moon. When measuring ozone, as the solar zenith angle increases the importance of instrumental stray light increases due to the increasing gradient in intensity towards longer wavelengths. One way to quantify the effect of stray light is to observe how the measurements deviate from linearity according to Beer's law at large ozone slant paths (Slavin, 1963). In a Brewer spectrophotometer, instrumental stray light arises mainly from the holographic diffraction grating and the collimating and focusing mirrors (Silva and Kirchhoff, 2004). Comparisons between measurements of SBr and DBr have been made by Bais et al. (1996), which show a much lower effect of stray light in measurements from the DBr. At wavelengths below 300 nm the SBr shows a 10 % underestimation in absolute irradiances when there is more than 1000 DU of ozone slant column present (Bais et al., 1996). At wavelengths above 300 nm, the presence of stray light may be problematic as well. This paper addresses this issue.
According to the WOUDC website, over 200 Brewers, single and double, are being used in more than 40 countries and 100 stations to measure ozone column amounts (Savastiouk, 2006; WOUDC, 2016). Many single Brewers are being replaced by double Brewers, and therefore, new measurements may show a false increase in ozone column amounts due to the effect of lower stray light in the DBr, particularly for measurements made at large solar zenith angles. In this paper a method is described that accounts for the effect of stray light in the data from the SBr to provide more accurate, reprocessed SBr data.
This paper presents a practical method to improve Brewer measurements by compensating for the effect of stray light. This model expands the Langley method to include a term representing the contribution of the nonlinearity in retrieved ozone resulting from the effects of stray light. The simulation of the Brewer performance using a physical model (Moeini et al., 2018) provided validation of the parameterization of the analysis model.
2.1 Physical model
In the physical model, the slit function of the instrument measured is used to characterize the stray light. For this purpose, a HeCd laser at 325 nm was used as a source and all wavelengths starting at 290 nm for the SBr and DBr were measured. The photon count rate measured by the Brewer is the integral of the spectral intensities on all wavelengths weighted by the slit function. Figure 1 shows the measured slit function of a SBr no. 009 and DBr no. 119 reported by Moeini et al. (2018). The dots show the measurements made by the HeCd laser and the solid lines show the best fit to the measurements. Under ideal conditions without the effect of stray light, the fit to the slit function would be a trapezoid with its wings extending to zero. However, in reality, as seen in the figure, the fit is trapezoidal from full width half maximum to the peak, but the wings form a Lorentzian function that extends to a horizontal line, which is not at zero. The difference between the wings of the SBr and DBr slit functions clearly shows the presence of stray light. The horizontal fit line is approximately 10^{−4} for the SBr and 10^{−6} for the DBr which shows that the effect of stray light is much more pronounced in the SBr measurements than the DBr (Moeini et al., 2018).
2.2 Mathematical model
To compensate for the effects of stray light on the measured ozone column amount, a new technique is described in this paper. A sensitive method for measuring the effect of instrumental stray light is using the deviation of the measurements from linearity according to Beer's law at large ozone slant paths (Slavin, 1963).
Beer's law states that the attenuation of light by a material increases exponentially with an increase in path length in a uniformly absorbing medium. Equation (1) shows this relation, where I_{0} is the intensity of light before entering the layer of material, I is the intensity of light after going through the layer of material in question, with an absorption coefficient of α, a column amount of X and a slant path length or air mass of μ. In the Brewer instrument, I_{0} is the intensity of the extraterrestrial light before it enters the atmosphere. This absolute intensity measured by the groundbased instrument relies on the knowledge of the extraterrestrial source of the light, which is difficult to determine due to scattering and absorption by clouds and aerosols. Therefore, the absorbances at two different wavelengths, one short and one long, are calculated to construct a measurement function which is independent of absolute intensity. This method is called differential absorption spectroscopy and is used for measurements made by the Dobson spectrophotometer. The Beer's law, which is linear with respect to air mass, can be calculated as shown in Eq. (2) where $\mathrm{\Delta}\mathit{\alpha}={\mathit{\alpha}}_{\mathrm{S}}{\mathit{\alpha}}_{\mathrm{L}}$ is the difference in absorption coefficients at I_{S} and I_{L}, the intensities of the short and long wavelengths, and I_{S0} and I_{L0}, the intensities of the short and long wavelength at zero air mass (no medium to absorb). The logarithm of the ratio of the two intensities is denoted as F and is called the absorption function shown in Eq. (3).
In the Brewer instrument, five wavelengths are used instead of just the two. These wavelengths are 306.3, 310.1, 313.5, 316.8, and 320.1 nm. A weighting is applied to the logarithm of the counts at each wavelength and a sum is used to calculate the absorption function. This method eliminates the absolute dependence on intensity and suppresses anything that is linear with respect to wavelength. This weighting will also produce measurements that are more sensitive to ozone. The weighting is shown in Eq. (4), where the absorption function, F, is calculated by the product of the weighting vector, w, by L, which is a 5×1 matrix composed of the logarithm of the counts measured at each wavelength, log(C_{λ}), multiplied by 10^{4}. A study done by Savastiouk and McElroy (2005) shows the calculations for deriving the weighting vector (Savastiouk and McElroy, 2005):
Under realistic conditions, Eq. (3) deviates from its linear form at large air mass values in the presence of instrumental stray light. To account for this nonlinearity, a new model for the absorption function measured by the instrument is defined. The form of the model, incorporating a correction of the absorption function which is approximately cubic in ozone, was determined empirically by testing different corrections.
This new model accounts for the nonlinearity in measuring ozone column amount. It also accounts for filter changes in the instrument. This instrument model is presented in Eq. (5), where F_{m} is the model absorption function, α is the absorption coefficient of ozone, μ is ozone air mass, X is ozone column amount, F_{0} is the absorption function at zero air mass, γ is the nonlinearity factor, b_{j} is the filter change factor and ND_{j} is the filter vector, which is zero for filter numbers (j) not used and has a value of 1 for filter numbers used. Different forms of the model have been experimented with and the model described in Eq. (5) was found to have good agreement with the observations.
The components of the model (Eq. 5) to be determined are ${v}_{k}=(x,\mathit{\gamma},{b}_{j}$, F_{0}), where k indicates the number of components of v that are to be retrieved. The Langley method is used to determine the k components of vector v by finding suitable values for v that minimize the square error between the model and the observations. This method is described in the following. In the first step, initial values for v_{k} are estimated. To predict an initial value for F_{0}, the conventional Langley plot is used, where the absorption function vs. air mass is plotted. As an example, the absorption function vs. air mass of the single Brewer no. 009 is plotted in Fig. 2, where the dots show the instrument measurements. The measured values tend to deviate from the linear model as the air mass increases. The plot is quite linear at air mass values smaller than 2; therefore, this part of the data is used to apply the leastsquares method using the linear model (Eq. 3) to find the slope, α.x and the intercept at zero air mass, F_{0}. The initial values of γ and b_{j} should be very small; therefore an initial value close to zero is used.
The nonlinear Langley method uses the leastsquares method to determine v_{k}. In this process the square error, SE, between the absorption function measured by the instrument, F_{i}, and the modelled absorption function, F_{mi}, will be minimized for all the N observations. The index i denotes the observation number from 1 to N. Therefore, the derivative of the difference between the observations and the modelled absorption function (Eq. 4) with respect to all the components of vector v_{k} should be zero as shown in Eq. (6).
Because the model is nonlinear in the ozone, several iterations of the solution are made to arrive at an accurate result. The maximum number of iterations can be set to 50 or more but the answer reaches a useful convergence between 5 and 10 iterations. Equations (6) to (13) show the Langley method loop and the estimation of v_{k}, where F_{i}−F_{mi} is replaced by ΔF_{i} in Eq. (7), M_{ik} is the derivative of the model absorption function (Eq. 4) with respect to the k components of v_{k} at each measurement point, i. M_{in} denotes the derivative of the model with respect to the nth component of v_{k} at each measurement point, i. The components of the Jacobian, M_{ik}, have been provided in the Appendix.
To improve accuracy, a weighting of the observations is needed. From Fig. 2, it can be seen that the measurements are denser at smaller air mass than at larger air mass; therefore the model will put more emphasis on these points and estimate the components of v_{k} accordingly. To reduce this problem the data are multiplied by $\frac{\mathrm{1}}{\mathit{\mu}}$ and the uncertainty, $\frac{\mathrm{1}}{{\mathit{\sigma}}^{\mathrm{2}}}$, associated with the counts as shown in Eq. (14) before applying the model. j is the number of wavelengths, which in this case is 1 to 5, Counts_{j} is the photon count for each of the five wavelengths and w_{j} is the weighting for each wavelength. With this weighting the difference between the model and measurements is exaggerated at large air mass and, therefore, carries more weight.
2.3 Deadtime correction
The dead time of an instrument is important to account for. In a Brewer instrument the detector photomultiplier pulses have a finite width of ∼30 ns. If two or more photons arrive at the detector within 30 ns they will be counted as one count. This will cause an error in counting the photons. This error associated with the dead time increases with count rate (Kerr, 2010). The dead time varies from instrument to instrument and can also change with time. For lowozone slant columns and highintensity solar radiation, a small change of 10 ns in the dead time can cause an error of up to 5 % in the measured total ozone column amount (Fountoulakis et al., 2016). The deadtime correction is modelled in Eq. (15) where N is the measured photon count, N_{0} is the corrected photon count and τ is the dead time (Kerr, 2010).
To calculate dead time, measurements are done with the Brewer through two exit slits simultaneously and separately through each exit slit. For example, measurements are done with exit slit 2 and 4 resulting in N_{2} and N_{4} and measurements are done with both exit slits open resulting in N_{2+4}. If Eq. (15) is written for each case, three equations with four unknowns result, which are N_{02}, N_{04}, N_{02+4} and τ. A fourth equation, ${N}_{\mathrm{02}+\mathrm{4}}={N}_{\mathrm{02}}+{N}_{\mathrm{04}}$, will help find the dead time (Kerr, 2010). This procedure is not done on a regular basis and the dead time of an instrument changes somewhat over time; therefore a correction prediction for dead time must be performed during the Langley process. To account for the deadtime prediction, when calculating the derivative of the square error, SE, the derivative of F_{i} with respect to dead time must also be considered. In this case, v_{k} is τ and thus Eq. (8) becomes Eq. (16), which shows how this derivative is included in the Langley method.
The Langley method was applied to the data collected by two Brewer instruments no. 009 (single Brewer) and no. 119 (double Brewer) stationed at Mauna Loa at −155.5^{∘} longitude and 19.5^{∘} latitude, and as an example the results for 14 July and 25 October 2010 are presented. The nonlinear mathematical model (Eq. 5) was applied to the data and model parameters where retrieved. The nonlinearity term, filter change term and deadtime term were applied in reverse to the data to calculate corrected values for the absorption function and ozone values in Dobson units.
The observed, modelled and corrected absorption function (F) for single Brewer no. 009 is plotted in Fig. 3. In the observed and modelled F, you can clearly see the filter change causing a jump in the intensity of light received by the detector. This sudden fluctuation in intensity is modelled using the nonlinear model and is corrected to give a more linear F with respect to ozone air mass. The same has been plotted for the double Brewer no. 119 in Fig. 4, where the observed modelled and corrected observations are very close to linear with respect to air mass, and thus the correction is less than 1 %.
Using the corrected absorption function, the corrected ozone values were calculated. The measured and corrected ozone column amounts in Dobson units vs. ozone air mass for both the single (no. 009) and double (no. 119) Brewers for 14 July and 25 October 2010 are illustrated in Figs. 5 and 7 respectively. There is some nonlinearity observed for the single Brewer at low air mass which is corrected to an amount slightly higher than the double Brewercorrected amounts. The ratios of single to double Brewer measurements before and after applying corrections for 14 July and 25 October 2010 are plotted in Figs. 6 and 8 respectively. The corrected ratio in red is closer to linearity than the uncorrected ratio in black. As a comparison the obtained parameters from the model and the standard parameters used to retrieve ozone for 14 July 2010 can be seen in Table 1.
There is a difference of up to 10 % in the ozone column amount measured by the single and double Brewers (single lower than double) as observed at large ozone slant paths and illustrated in Fig. 5 for 14 July 2010 and Fig. 7 for 25 October 2010. The mathematical model presented in this paper accounts for this nonlinearity in the data collected by single Brewers. Applying corrections to the historical data of single Brewers is essential, as it will eliminate errors due to the effect of stray light in the ozone measurements.
The double Brewer more accurately rejects straylight effects (Gröbner et al., 1998), and as the single Brewers in ozone measurement stations around the world are being replaced by the more advanced double Brewer, a slight increase in the ozone amounts may be observed, particularly in those months in which observations must necessarily be made at large solar zenith angles. This may lead to a false assumption that the total ozone column amounts are showing an increasing trend when this may not be the case.
The departure from linearity is more prominent for large air mass (μ) values as is seen in Figs. 3, 5 and 7. This nonlinearity at large air masses becomes more prominent in highlatitude regions such as the Arctic, where the sun is at large solar zenith angles most of the year. It is also difficult to take reference Brewer instruments to the Arctic for calibration purposes. To provide an accurate onsite method to correct the ozone data for nonlinearity, the process described in this paper is recommended.
The future steps for this research will be to apply this method to historical databases and produce an applicable software that would output corrected daily ozone values. Data processing using the new methodology leads to a more accurate absolute calibration on single Brewers and to a way to properly transfer calibration constants between single and double Brewers.
Data available at

https://doi.org/10.5281/zenodo.2531911 (Vaziri Zanjani, 2019a),

https://doi.org/10.5281/zenodo.2531905 (Vaziri Zanjani, 2019b),

https://doi.org/10.5281/zenodo.2531955 (Vaziri Zanjani, 2019c) and

https://doi.org/10.5281/zenodo.2531951 (Vaziri Zanjani, 2019d).
ZVZ contributed by developing the mathematical model and MATLAB software to analyze and correct for the nonlinearity caused by stray light present in the ozone measurements by the Brewer instrument. OM developed the physical model of the Brewer instrument. TME, DB and VS assisted in writing the paper and provided valuable discussions about the Brewer retrievals and measurements.
The authors declare that they have no conflict of
interest.
Edited by: Andreas
Hofzumahaus
Reviewed by: three anonymous referees
Bais, A. F., Zerefos, C. S., and McElroy, C. T.: Solar UVB measurements with the double and singlemonochromator Brewer ozone spectrophotometers, Geophys. Res. Lett., 23, 833–836, https://doi.org/10.1029/96GL00842, 1996.
Fioletov, V. E., Kerr, J. B., Wardle, D. I., and Wu, E.: Correction of stray light for the Brewer single monochromator, Proc. Quadrennial Ozone Symposium, 21, 369–370, 2000.
Fountoulakis, I., Redondas, A., Bais, A. F., RodriguezFranco, J. J., Fragkos, K., and Cede, A.: Dead time effect on the Brewer measurements: correction and estimated uncertainties, Atmos. Meas. Tech., 9, 1799–1816, https://doi.org/10.5194/amt917992016, 2016.
Gröbner, J., Wardle, D. I., McElroy, C. T., and Kerr, J. B.: Investigation of the wavelength accuracy of brewer spectrophotometers, Appl. Optics, 37, 8352–8360, 1998.
Josefsson, W. A. P.: Focused Sun Observations Using a Brewer Ozone Spectrophotometer, J. Geophys. Res., 97, 15813–15817, 1992.
Karppinen, T., Redondas, A., García, R. D., Lakkala, K., McElroy, C. T., and Kyrö, E. Compensating for the effects of stray light in singlemonochromator brewer spectrophotometer ozone retrieval, Atmos. Ocean, 53, 66–73, https://doi.org/10.1080/07055900.2013.871499, 2015.
Kerr, J. B.: The Brewer Spectrophotometer, in: UV Radiation in Global Climate Change, edited by: Gao, W., Schmoldt, D. L., and Slusser, J. R., Tsinghua University Press, Beijing and Springer, 160–191, 2010.
Kerr, J. B., McElroy, C. T., Wardle, D. I., Olafson, R. A., and Evans, W. F. J. The automated Brewer spectrophotometer, in: Quadrennial Ozone Symposium, edited by: Zerefos, C. S. and Ghazi, A., Halkidiki, Greece, 396–401, 1984.
Moeini, O., Vaziri, Z., McElroy, C. T., Tarasick, D. W., Evans, R. D., Petropavlovskikh, I., and Feng, K.H.: The Effect of Instrumental Stray Light on Brewer and Dobson Total Ozone Measurements, Atmos. Meas. Tech. Discuss., https://doi.org/10.5194/amt20182, in review, 2018.
Redondas, A., Evans, R., Stuebi, R., Köhler, U., and Weber, M.: Evaluation of the use of five laboratorydetermined ozone absorption cross sections in Brewer and Dobson retrieval algorithms, Atmos. Chem. Phys., 14, 1635–1648, https://doi.org/10.5194/acp1416352014, 2014.
Redondas, A., Carreño, V., LeónLuis, S. F., HernándezCruz, B., LópezSolano, J., RodriguezFranco, J. J., Vilaplana, J. M., Gröbner, J., Rimmer, J., Bais, A. F., Savastiouk, V., Moreta, J. R., Boulkelia, L., Jepsen, N., Wilson, K. M., Shirotov, V., and Karppinen, T.: EUBREWNET RBCCE Huelva 2015 Ozone Brewer Intercomparison, Atmos. Chem. Phys., 18, 9441–9455, https://doi.org/10.5194/acp1894412018, 2018.
Rimmer, J. S., Redondas, A., and Karppinen, T.: EuBrewNet – A European Brewer network (COST Action ES1207), an overview, Atmos. Chem. Phys., 18, 10347–10353, https://doi.org/10.5194/acp18103472018, 2018.
Savastiouk, V.: Improvements To the DirectSun Ozone Observations Taken With the Brewer Spectrophotometer, York University, available at: https://www.esrl.noaa.gov/gmd/grad/neubrew/docs/publications/VladimirSavastiouk_PhD_thesis.pdf (last access: 2018), 2006.
Savastiouk, V. and McElroy, C. T.: Brewer spectrophotometer total ozone measurements made during the 1998 Middle Atmosphere Nitrogen Trend Assessment (MANTRA) Campaign, Atmos. Ocean, 43, 315–324, https://doi.org/10.3137/ao.430403, 2005.
Silva, A. A. and Kirchhoff, V. W. J. H.: Aerosol optical thickness from Brewer spectrophotometers and an investigation into the straylight effect, Appl. Optics, 43, 2484–2489, 2004.
Slavin, W.: Stray Light in Ultraviolet, Visible and NearInfrared Spectrophotometry, Anal. Chem., 35, 561–566, 1963.
Vaziri Zanjani, Z.: Brewer 009_25 Oct 2010 [Data set], Zenodo, https://doi.org/10.5281/zenodo.2531911, 2019a.
Vaziri Zanjani, Z.: Brewer 009_14 Jul 2010 [Data set], Zenodo, https://doi.org/10.5281/zenodo.2531905, 2019b.
Vaziri Zanjani, Z.: Brewer 119_14 Jul 2010_AM, https://doi.org/10.5281/zenodo.2531955, 2019c.
Vaziri Zanjani, Z.: Brewer 119_ 25 Oct 2010_ AM [Data set], Zenodo, https://doi.org/10.5281/zenodo.2531951, 2019d.
WOUDC: instrument list is available at: http://www.woudc.org/data/instruments/ (last access: 1 May 2017), 2016.