| Summary: | Long-term continuous time series of <inline-formula><math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><semantics><msub><mi>SO</mi><mn>2</mn></msub></semantics></math></inline-formula> emissions are considered critical elements of both volcano monitoring and basic research into processes within magmatic systems. One highly successful framework for computing these fluxes involves reconstructing a representative time-averaged <inline-formula><math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><semantics><msub><mi>SO</mi><mn>2</mn></msub></semantics></math></inline-formula> plume from which to estimate the <inline-formula><math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><semantics><msub><mi>SO</mi><mn>2</mn></msub></semantics></math></inline-formula> source flux. Previous methods within this framework have used ancillary wind datasets from reanalysis or numerical weather prediction (NWP) to construct the mean plume and then again as a constrained parameter in the fitting. Additionally, traditional <inline-formula><math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><semantics><msub><mi>SO</mi><mn>2</mn></msub></semantics></math></inline-formula> datasets from ultraviolet (UV) sensors lack altitude information, which must be assumed, to correctly calibrate the <inline-formula><math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><semantics><msub><mi>SO</mi><mn>2</mn></msub></semantics></math></inline-formula> data and to capture the appropriate NWP wind level which can be a significant source of error. We have made novel modifications to this framework which do not rely on prior knowledge of the winds and therefore do not inherit errors associated with NWP winds. To perform the plume rotation, we modify a rudimentary computer vision algorithm designed for object detection in medical imaging to detect plume-like objects in gridded <inline-formula><math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><semantics><msub><mi>SO</mi><mn>2</mn></msub></semantics></math></inline-formula> data. We then fit a solution to the general time-averaged dispersion of <inline-formula><math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><semantics><msub><mi>SO</mi><mn>2</mn></msub></semantics></math></inline-formula> from a point source. We demonstrate these techniques using <inline-formula><math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><semantics><msub><mi>SO</mi><mn>2</mn></msub></semantics></math></inline-formula> data generated by a newly developed probabilistic layer height and column loading algorithm designed for the Cross-track Infrared Sounder (CrIS), a hyperspectral infrared sensor aboard the Joint Polar Satellite System’s Suomi-NPP and NOAA-20 satellites. This <inline-formula><math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><semantics><msub><mi>SO</mi><mn>2</mn></msub></semantics></math></inline-formula> data source is best suited to flux estimates at high-latitude volcanoes and at low-latitude, but high-altitude volcanoes. Of particular importance, IR <inline-formula><math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><semantics><msub><mi>SO</mi><mn>2</mn></msub></semantics></math></inline-formula> data can fill an important data gap in the UV-based record: estimating <inline-formula><math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><semantics><msub><mi>SO</mi><mn>2</mn></msub></semantics></math></inline-formula> emissions from high-latitude volcanoes through the polar winters when there is insufficient solar backscatter for UV sensors to be used.
|
|---|