Radiative transfer corrections for accurate spectroscopic measurements of volcanic gas emissions

There is widespread use of passive remote sensing techniques to quantify trace gas column densities in volcanic plumes utilizing scattered sunlight as a light source. Examples include passive DOAS, COSPEC, and the SO₂ camera. In order to calculate trace gas concentrations or volcanic emission fluxes, knowledge about the optical path through the plume is necessary. In the past, a straight photon path through the plume has always been assumed although it was known that this is not always true. Here we present the results of model studies conducted specifically to quantify the effects of realistic radiative transfer in and around volcanic plumes on ground-based remote sensing measurements of SO₂. The results show that measurements conducted without additional information on average photon paths can be inaccurate under certain conditions, with possible errors spanning more than an order of magnitude. Both over and underestimation of the true column density can occur. Actual errors depend on parameters such as distance between instrument and plume, plume SO₂ concentration, plume aerosol load, as well as aerosol conditions in the ambient atmosphere. As an example, a measurement conducted with an SO₂ camera is discussed, the results of which can only be correctly interpreted if realistic radiative transfer is considered. Finally, a method is presented which for the first time allows the retrieval of actual average photon paths in spectroscopic (i.e. DOAS) measurements of adequate resolution. By allowing for a wavelength dependent column density during the evaluation of DOAS measurements, we show how radiative transfer effects can be corrected using information inherently available in the measured spectra, thus greatly enhancing the accuracy of DOAS measurements of volcanic emissions.     

Imaging DOAS for volcanological applications

The simultaneous quantitative determination of two-dimensional bromine monoxide (BrO) and sulphur dioxide (SO₂) distributions in volcanic gas plumes is described. Measurements at the fumarolic field on the island Vulcano (autumn 2004) and in the plume of Mt. Etna volcano (spring 2005) were carried out with an Imaging DOAS instrument. The SO₂ fluxes of several fumaroles were estimated from two-dimensional distributions of SO₂. Additionally, the first two-dimensional distributions of BrO within a volcanic plume were successfully retrieved. Slant column densities of up to 2.6 × 10¹⁴ molecules per square centimetre were detected in the plume of Mt. Etna. The investigation of the BrO/SO₂ ratio, calculated from the two-dimensional distributions of SO₂ and BrO, shows an increase from the centre to the edge of the volcanic plume. These results have significance for the involvement of ozone during BrO formation processes in volcanic emissions.     

IDOAS: A new monitoring technique to study the 2D distribution of volcanic gas emissions

Imaging Differential Optical Absorption Spectroscopy (IDOAS) is an optical remote-sensing method using scattered sunlight as light source. It combines a “pushbroom” imaging spectrometer with the DOAS technique and thus allows imaging two-dimensional trace gas distributions, e.g., in volcanic plumes. The highly sensitive and specific detection of many trace gases simultaneously (specific molecules, not just elements, e.g. SO₂, BrO, NO₂, O₃, HCHO, etc.) is possible, and the temporal and spatial variation of these gases can be measured. The IDOAS system presented here enables the taking of two-dimensional images of trace gas distributions in a volcanic plume with a spatial resolution of 100 pixels horizontally × 64 pixels vertically, each with a field of view of 0.087° in horizontal and 0.208° in vertical directions. Therefore, IDOAS provides useful information about the chemical composition and chemical variability in a volcanic plume and allows studying plume dispersal and chemical transformations. The technique was applied to map the SO₂ distribution in the plume of Mt. Etna volcano for the first time in October 2003.     

A new technique to estimate volcanic gas composition: plume measurements with a portable multi-sensor system

A portable multi-sensor system was developed to measure volcanic plumes in order to estimate the chemical composition and temperature of volcanic gases. The multi-sensor system consists of a humidity–temperature sensor, SO₂ electrochemical sensor, CO₂ IR analyzer, pump and flow control units, pressure sensor, data logger, and batteries; the whole system is light (∼5 kg) and small enough to carry in a medium-size backpack. Volcanic plume is a mixture of atmosphere and volcanic gas; therefore volcanic gas composition and temperature can be estimated by subtracting the atmospheric gas background from the plume data. In order to obtain the contrasting data of the plume and the atmosphere, measurements were repeated in and out of the plume. The multi-sensor technique was applied to measure the plume of Tarumae, Tokachi, and Meakan volcanoes, Hokkaido, Japan. Repeated measurements at each volcano gave a consistent composition with ±10–30% errors, depending on the stability of the background atmospheric conditions. Fumarolic gas samples were also collected at the Tokachi volcano by a conventional method, and we found a good agreement (the difference <10%) between the composition estimated by the multi-sensor technique and conventional method. Those results demonstrated that concentration ratios of major volcanic gas species (i.e., H₂O, CO₂, and SO₂) and temperature can be estimated by the new technique without any complicated chemical analyses even for gases emitted from an inaccessible open vent. Estimation of a more detailed gas composition can be also achieved by the combination of alkaline filter techniques to measure Cl/F/S ratios in the plume and other sensors for H₂S and H₂.     

A CO2-rich gas trigger of explosive paroxysms at Stromboli basaltic volcano,Italy

In addition to rhythmic slug-driven Strombolian activity, Stromboli volcano occasionally produces discrete explosive paroxysms (2 per year on average for the most frequent ones) that constitute a major hazard and whose origin remains poorly elucidated. Partial extrusion of the volatile-rich feeding basalt as aphyric pumice during these events has led to consider their triggering by the fast ascent of primitive magma blobs from possibly great depth. Here I examine and discuss the alternative hypothesis that most of the paroxysms could be triggered and driven by the fast upraise of CO₂-rich gas pockets generated by bubble foam growth and collapse in the sub-volcano plumbing system. Data for the SO₂ and CO₂ crater plume emissions are used to show that Stromboli's feeding magma may originally contain as much as 2 wt.% of carbon dioxide and early coexists with an abundant CO₂-rich gas phase with high CO₂/SO₂ molar ratio (≥ 60 at 10 km depth below the vents, compared to ～ 7 in time-averaged crater emissions). Pressure-related modelling indicates that the time-averaged crater gas composition and output are well accounted for by closed system decompression of the basalt–gas mixture until the volcano–crust interface (～ 3 km depth), followed by open degassing and crystallization in the volcano conduits. However, both the low viscosity and high vesicularity of the basaltic magma permit bubble segregation and bubble foam growth at deep sill-like feeder discontinuities and at shallower physical boundaries (such as the volcano–crust interface) where the gas-rich aphyric basalt interacts with the unerupted crystal-rich and viscous magma drained back from the volcano conduits. Gas pressure build-up and bubble foam collapse at these boundaries will intermittently trigger the sudden upraise of CO₂-rich gas blobs that constitute the main driving force of the paroxysms. Deeper-sourced gas blobs, driving the most powerful explosions, will be the richest in CO₂ and have highest CO₂/SO₂ ratios. This mechanism is shown to account well for the dynamic, seismic and petrologic features of Stromboli's paroxysms and, hence, to provide a potential alternative interpretation for their genesis and their forecasting. Enhanced bubble foam leakage prior to a paroxysm, or foam emptying in several steps, should lead indeed to precursory upstream of CO₂-rich gas and increasing CO₂/SO₂ ratio in crater plume emissions. The recent detection of such signals prior to two explosions in December 2006 and March 2007 strongly supports this expectation and the model proposed in this study.     

Hydrogen emissions from Erebus volcano, Antarctica

The continuous measurement of molecular hydrogen (H₂) emissions from passively degassing volcanoes has recently been made possible using a new generation of low-cost electrochemical sensors. We have used such sensors to measure H₂, along with SO₂, H₂O and CO₂, in the gas and aerosol plume emitted from the phonolite lava lake at Erebus volcano, Antarctica. The measurements were made at the crater rim between December 2010 and January 2011. Combined with measurements of the long-term SO₂ emission rate for Erebus, they indicate a characteristic H₂ flux of 0.03?kg s^–1 (2.8?Mg? day^–1). The observed H₂ content in the plume is consistent with previous estimates of redox conditions in the lava lake inferred from mineral compositions and the observed CO₂/CO ratio in the gas plume (～0.9 log units below the quartz–fayalite–magnetite buffer). These measurements suggest that H₂ does not combust at the surface of the lake, and that H₂ is kinetically inert in the gas/aerosol plume, retaining the signature of the high-temperature chemical equilibrium reached in the lava lake. We also observe a cyclical variation in the H₂/SO₂ ratio with a period of ～10?min. These cycles correspond to oscillatory patterns of surface motion of the lava lake that have been interpreted as signs of a pulsatory magma supply at the top of the magmatic conduit.     

Fumarole/plume and diffuse CO2 emission from Sierra Negra caldera, Galapagos archipelago

Measurements of visible and diffuse gas emission were conducted in 2006 at the summit of Sierra Negra volcano, Galapagos, with the aim to better characterize degassing after the 2005 eruption. A total SO₂ emission of 11?±?2?t day^?1 was derived from miniature differential optical absorption spectrometer (mini-DOAS) ground-based measurements of the plume emanating from the Mini Azufral fumarolic area, the most important site of visible degassing at Sierra Negra volcano. Using a portable multigas system, the H₂S/SO₂, CO₂/SO₂, and H₂O/SO₂ molar ratios in the Mina Azufral plume emissions were found to be 0.41, 52.2, and 867.9, respectively. The corresponding H₂O, CO₂, and H₂S emission rates were 562, 394, and 3?t day^?1, respectively. The total output of diffuse CO₂ emissions from the summit of Sierra Negra volcano was 990?±?85?t day^?1, with 605?t day^?1 being released by a deep source. The diffuse-to-plume CO₂ emission ratio was about 1.5. Mina Azufral fumaroles released gasses containing 73.6?mol% of H₂O; the main noncondensable components amounted to 97.4?mol% CO₂, 1.5?mol% SO₂, 0.6?mol% H₂S, and 0.35?mol%?N₂. The higher H₂S/SO₂ ratio values found in 2006 as compared to those reported before the 2005 eruption reveal a significant hydrothermal contribution to the fumarolic emissions. ³He/⁴He ratios measured at Mina Azufral fumarolic discharges showed values of 17.88?±?0.25?R _A , indicating a mid-ocean ridge basalts (MORB) and a Galapagos plume contribution of 53 and 47?%, respectively.     

Impact of volcanic fluoride and SO2 emissions from moderated activity volcanoes on the surrounding vegetation

Studies in the regions of the volcanoes Etna (Italy) and Masaya (Nicaragua) show that the continuous emissions of gaseous pollutants (HF and SO₂) from moderated activity volcanoes causes a chronic pollution in the surrounding vegetation with certain economical and ecological consequences.Reciprocally the measure of the pollutants in the plants growing in volcanic regions may be a simple and fast method to investigate some characteristics of the volcanic plume: for example, intensity of the emissions of gas, direction and extent of the plume.     

Interactions between gravity currents and convective dissolution

Geological storage of carbon dioxide (CO₂) is a promising technology for reducing atmospheric emissions. The large discrepancy in the time- and length-scales between up-dip migration of buoyant supercritical CO₂ and the sinking fingers of dissolved CO₂ poses a challenge for numerical simulations aimed at describing the fate of the plume. Hence, several investigators have suggested methods to simplify the problem, but to date there has been no reference solution with which these simplified models can be compared. We investigate the full problem of Darcy-based two-phase flow with gravity-current propagation and miscible convective mixing, using high-resolution numerical simulations. We build on recent developments of the Automatic Differentiation - General Purpose Research Simulator (AD-GPRS) at Stanford. The results show a CO₂ plume that travels for 5000 years reaching a final distance of 14 km up-dip from the injection site. It takes another 2000 years before the CO₂ is completely trapped as residual (40%) and dissolved (60%) CO₂. Dissolution causes a significant reduction of the plume speed. While fingers of dissolved CO₂ appear under the propagating gravity current, the resident brine does not become fully saturated with CO₂ anywhere under the plume. The overall mass transfer of CO₂ into the brine under the plume remains practically constant for several thousands of years. These results can be used as a benchmark for verification, or improvements, of simplified (reduced-dimensionality, upscaled) models. Our results indicate that simplified models need to account for: (i) reduced dissolution due to interaction with the plume, and (ii) gradual reduction of the local dissolution rate after the fingers begin to interact with the bottom of the aquifer.     

Migration of seismic activity during the 1998 volcanic unrest at Iwate volcano,northeastern Japan,with reference to P and S wave velocity anomaly and crustal deformation

We describe the seismicity at Iwate volcano, northeastern Japan, during the volcanic unrest of 1998 with reference to a three-dimensional P and S wave velocity model from tomographic analysis. The abnormal seismic activity beneath Iwate volcano started under the caldera in February, 1998 and migrated westward in the period February to August, 1998. Previous geodetic modeling [Sato and Hamaguchi, Chikyu Monthly 21 (1999) 312–317] suggested the growth of a dike in the time of the seismic activity. Comparing the seismicity and dike extension with the tomographic images of the P and S wave velocity structure, we find that the trace of the growing dike coincides with the region of the high V_p and high V_p/V_s ratio beneath the volcano. The seismic and geodetic data are consistent with an intrusion of magma or other fluid under the caldera in 1998. Another pressure source causing the predominant crustal deformation at Iwate volcano was detected from geodetic data, which was located in the region with high V_p/V_s ratio under the western end of the volcano through the period from February to August. It is suggested that the activation of the point pressure source probably associated with the inflation of a hot fluid reservoir relate to a geothermal region adjacent to the western edge of the volcano.     

参量阵浅地层剖面仪在海底羽状流探测中的应用——以ATLAS P70在马克兰海域调查为例

海底渗漏的羽状流是沉积层赋存天然气水合物的重要证据之一,基于非线性水声学原理的参量阵浅地层剖面仪作为海洋探测的重要设备,对于获得羽状流在水体中的物性特征和渗漏点的浅地层信息有着重大意义.本文根据ATLAS P70浅地层剖面仪在马克兰海域调查中得到的浅地层剖面数据,结合多道地震数据、多波束数据以及地质样品等资料,刻画了研究区内羽状流形态特征,分析了羽状流区海底地层流体运移的通道以及近海底微地形地貌特征.通过研究发现在羽状流区伴随泥火山喷发,自生碳酸盐岩发育,剥蚀海底松散沉积物形成大小不一的麻坑,滋生生物群落等特征.反映在浅剖初始高频(20 kHz)数据界面上羽状流表现为柱状浑浊反射异常,形态呈火焰状,高度由80 m到1500 m不等;对应在次级低频(4 kHz)信号界面可以清晰显示流体渗漏的浅地层结构特征,从中不仅可以识别出流体的运移通道,如泥火山和管状通道等,而且揭示了流体逸散的残留地貌,如麻坑构造和海底滑坡等.本文依托参量阵浅地层剖面数据,对巴基斯坦马克兰海域羽状流有了较全面的认识,为天然气水合物的研究垫定了基础.     

The impact of moderate-scale explosive eruptions on stratospheric gas injections

Volcanic gases such as SO ₂, H ₂S, HCl and COS emitted during explosive eruptions significantly affect atmospheric chemistry and therefore the Earth's climate. We have evaluated the dependence of volcanic gas emission into the atmosphere on altitude, latitude, and tectonic setting of volcanoes and on the season in which eruptions occurred. These parameters markedly influence final stratospheric gas loading. The latitudes and altitudes of 360 active volcanoes were compared to the height of the tropopause to calculate the potential quantity of volcanic gases injected into the stratosphere. We calculated a possible stratospheric gas loading based on different volcanic plume heights (6, 10, and 15 km) generated by moderate-scale explosive eruptions to show the importance of the actual plume height and volcano location. At a plume height of 15 km for moderate-scale explosive eruptions, a volcano at sea level can cause stratospheric gas loading because the maximum distance to the tropopause is 15–16 km in the equatorial region (0–30°). Eruptions in the tropics have to be more powerful to inject gas into the stratosphere than eruptions at high latitudes because the tropopause rises from ca. 9–11 km at the poles to 15–16 km in the equatorial region (0–30°N and S). The equatorial region is important for stratospheric gas injection because it is the area with the highest frequency of eruptions. Gas injected into the stratosphere in equatorial areas may spread globally into both hemispheres.     

Automated,high time-resolution measurements of SO<Subscript>2</Subscript> flux at Soufrière Hills Volcano,Montserrat

We report here the first results from an automated, telemetered UV scanning spectrometer system for monitoring SO ₂ emission rates at Soufrière Hills Volcano, Montserrat. Two spectrometers receive light by way of a motor-driven stepping prism and telescope in order to make vertical scans of the volcanic plume. Spectral data from these spectrometers, situated 2,800 m apart and 4,500 m from the volcano, are relayed back to the observatory every 4–5 s via radio modems. A full scan of the plume is accomplished every 1–6 min by the (time-synchronised) spectrometers and a SO ₂ emission rate is calculated using the SO ₂ slant concentrations, scan angles and plume speeds estimated from the wind speed from a telemetered weather station near to the volcano. The plume's position and dimensions are calculated using the angular data from the two spectrometers. The plume height varies significantly diurnally and seasonally and is important in order to minimise the error on SO ₂ emission rates. The new scanning system (Scanspec) provides SO ₂ emission rates from 08:00 to 16:00 h local time every day. Preliminary results highlight a number of features of the SO ₂ time series and plume dynamics and give our first indications of the errors and limits of detection of this system. SO ₂ emission rates vary widely on all time scales (minutes, days, months). This new system has already provided the first real and consistent indication that SO ₂ emission rates vary on a minutes to hours basis, which can be correlated with volcanic activity (for example, rockfall and pyroclastic flow activity). It is anticipated that this system at Soufrière Hills will yield information on shallow processes occurring on short time scales (periods of minutes to hours) as well as deep processes relating to magma supply rates, which will be associated with longer wavelength SO ₂ signals of weeks to months.     

New clues on the contribution of Earth’s volcanism to the global mercury cycle

Active volcanoes are thought to be important contributors to the atmospheric mercury (Hg) budget, and this chemical element is one of the most harmful atmospheric pollutants, owing to its high toxicity and long residence time in ecosystems. There is, however, considerable uncertainty over the magnitude of the global volcanic Hg flux, since the existing data on volcanogenic Hg emissions are sparse and often ambiguous. In an attempt to extend the currently limited dataset on volcanogenic Hg emissions, we summarize the results of Hg flux measurements at seven active open-conduit volcanoes; Stromboli, Asama, Miyakejima, Montserrat, Ambrym, Yasur, and Nyiragongo.. Data from the dome-building Soufriere Hills volcano are also reported. Using our determined mercury to SO₂ mass ratios in tandem with the simultaneously-determined SO₂ emission rates, we estimate that the 7 volcanoes have Hg emission rates ranging from 0.2 to 18 t yr^-1 (corresponding to a total Hg flux of ~41 t·yr^-1). Based on our dataset and previous work, we propose that a Hg/SO₂ plume ratio ~10^-5 is best-representative of gas emissions from quiescent degassing volcanoes. Using this ratio, we infer a global volcanic Hg flux from persistent degassing of ~95 t·yr^-1 .     

The gas content and buoyancy of strombolian ash plumes

Plinian plumes erupt with a bulk density greater than that of air, and depend upon air entrainment during their gas-thrust phase to become buoyant; if entrainment is insufficient, the column collapses into a potentially deadly pyroclastic flow. This study shows that strombolian ash plumes can be erupted in an initially buoyant state due to their extremely high initial gas content, and in such cases are thus impervious to column collapse. The high gas content is a consequence of decoupled gas rise in the conduit, in which particles are ultimately incidental. The relations between conduit gas flow, eruption style and plume density are explored here for strombolian scenarios and contrasted with conventional wisdom derived from plinian eruptions. Considering the inherent relation between gas content and initial plume density together with detailed measurements of plume velocities can help unravel ambiguities surrounding conduit processes, eruption styles and hazards at poorly understood volcanoes. Analysis of plume dynamics at Santiaguito volcano, Guatemala adds further support for a model involving decoupled gas rise in the conduit.     

Degassing patterns of Tungurahua volcano (Ecuador) during the 1999–2006 eruptive period,inferred from remote spectroscopic measurements of SO2 emissions

This paper presents the results of 7 years (Aug. 1999–Oct. 2006) of SO₂ gas measurements during the ongoing eruption of Tungurahua volcano, Ecuador. From 2004 onwards, the operation of scanning spectrometers has furnished high temporal resolution measurements of SO₂ flux, enabling this dataset to be correlated with other datasets, including seismicity. The emission rate of SO₂ during this period ranges from less than 100 to 35,000 tonnes/day (t d^− 1) with a mean daily emission rate of 1458 t d^− 1 and a standard deviation of ± 2026 t d^− 1. Average daily emissions during inferred explosive phases are about 1.75 times greater than during passive degassing intervals. The total amount of sulfur emitted since 1999 is estimated as at least 1.91 Mt, mostly injected into the troposphere and carried westwards from the volcano. Our observations suggest that the rate of passive degassing at Tungurahua requires SO₂ exsolution of an andesitic magma volume that is two orders of magnitude larger than expected for the amount of erupted magma. Two possible, and not mutually exclusive, mechanisms are considered here to explain this excess degassing: gas flow through a permeable stagnant-magma-filled conduit and gas escape from convective magma overturning in the conduit. We have found that real-time gas monitoring contributes significantly to better eruption forecasting at Tungurahua, because it has provided improved understanding of underlying physical mechanisms of magma ascent and eruption.     

Estimation of SO<Subscript>2</Subscript> abundance in the eruption plume of Mt. Etna using two MIVIS thermal infrared channels: a case study from the Sicily-1997 Campaign

In this paper, an algorithm is developed based on the split-window technique, to estimate the SO₂ abundance in the plume of Mt. Etna volcano using the multispectral infrared and visible imaging spectrometer (MIVIS). The MIVIS data were remotely sensed in the thermal infrared (TIR) during the Sicily-1997 Campaign. In this study, the MODTRAN 3.5 code has been used to simulate the radiance at the sensor; the radiative transfer model was input along with the data of radio-sounding performed simultaneously with the MIVIS flight using a mobile radio-theodolite. From the SO₂ map, derived from the MIVIS image, the SO₂ flux along the axis of the plume was computed knowing the wind speed at the plume altitude. The SO₂ flux is variable along the plume axis. The average SO₂ flux (about 45 kg s^-1 on 12 June and about 30 kg s^-1 on 16 June) emitted from the vents is compared with the correlation spectrometer (COSPEC) measurements carried out by other teams (from the ground and from a light aircraft flying under the plume) during the MIVIS flight. Finally, by means of this algorithm it should be easier, with respect to the previously described procedure to monitor the SO₂ flux of a specific volcano such as Mt. Etna.     

Diffuse CO2 efflux from Iwojima volcano,Izu-Ogasawara arc,Japan

Iwojima volcano, located on the southernmost part of the Izu-Ogasawara arc, is characterized by the extrusion of trachyte or trachy andesite lavas and pyroclastic rocks of Holocene and surface thermal manifestations. Small phreatic explosions have been recorded frequently during the last 100 years with the most recent in 1999 and 2001. In order to elucidate the behavior of volcanic volatiles and to assess the potential activity of this volcano, diffuse CO₂ efflux, CO₂ content and δ¹³C–CO₂ in soil gas, and soil temperature at 30 cm depth were measured at 272 sites in March 2000, 112 sites in December 2000 and 40 sites in December 2001. We found that high CO₂ efflux values, of more than 100 g m⁻² day⁻¹, occurred at several locations on Motoyama volcano corresponding with high soil temperatures (more than 60 °C at 30 cm depth) region and with areas where CO₂ with magmatic δ¹³C was observed. Here, the magmatic δ¹³C determined for fumarolic CO₂ data ranged from −2‰ to +3‰, which is clearly higher than magmatic gas values (−8‰ to −2‰) typically found in island arc settings around the world. However, this can be explained in terms of carbon-isotope fractionation between calcite and CO₂ under subsurface temperature and pressure conditions at Iwojima. A total efflux of CO₂ for Iwojima volcano is estimated to be 760 t day⁻¹, with a magmatic contribution of about 450 t day⁻¹. This value is rather high compared with other volcanoes in island arc settings. Since Iwojima has no visible plume, almost all volcanic CO₂ is released as diffuse efflux through the volcanic edifice.     

Comparison of SO<Subscript>2</Subscript> column retrievals from BRD and DOAS algorithms

Atmospheric SO₂ has a significant impact on the urban environment and global climate. Band Residual Difference Algorithm (BRD) and Differential Optical Absorption Spectroscopy (DOAS) were used respectively by NASA and ESA science team to derive SO₂ columns from satellite observations, but there are few studies on the comparison and validation of BRD and DOAS SO₂ retrievals under the same observation conditions. In this study, the radiative transfer model SCIATRAN was firstly used to validate the accuracies of BRD and DOAS SO₂ retrievals, and analyse the uncertainty of SO₂ retrieval caused by band selection, O3 absorption, aerosol, surface reflectance, solar and viewing zenith angle. Finally, BRD and DOAS algorithms were applied to the same radiances from satellite observations, and comparisons of BRD and DOAS SO₂ retrievals were conducted over volcanic eruption and North China. Results show that, for the case with low SO₂ columns, BRD SO₂ retrievals have higher retrieval accuracy than DOAS, but typical seasonal variation with high SO₂ column in winter and low in summer can be more clearly discernible in DOAS SO₂ retrievals than BRD from satellite observations. For the case with high SO₂ columns, the differences between BRD (310.8–314.4 nm) and DOAS (315–327 nm) retrievals are large, and the value and accuracy of BRD (310.8–314.4 nm) SO₂ retrievals are lower than those of DOAS (315–327 nm) retrievals. Compared with the SO₂ inputs in forward model, both BRD (310.8–314.4 nm) and DOAS (315–327 nm) SO₂ retrievals are underestimated for the case with high SO₂ columns. The selection of wavelength range can significantly affect the accuracy of SO₂ retrieval. The error of BRD SO₂ retrieval from 310.8–314.4 nm is lower than other bands in the ultraviolet spectral region (306–327 nm). The increase of wavelength in the ultraviolet spectral region 306–330 nm can reduce the underestimation of DOAS SO₂ retrievals in the case of high SO₂ column, but slight overestimation of SO₂ retrieval is found from the 315–327 nm range in the case of low SO₂ column. The values of BRD and DOAS SO₂ retrieval decrease with atmospheric O₃ column and aerosol optical depth increasing, but increase with surface reflectance increasing. Large solar zenith angle and viewing zenith angle can introduce more errors to the BRD and DOAS SO₂ retrievals. This study is important for the improvement of retrieval algorithm and the application of SO₂ products from satellite observations.