SEARCH · Search NASA

Results for “Localization Error”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 991 records · Page 55

All sources NASA NTRS DOE OSTI

NASA NTRSDate not supplied

SO2 Plumes Observation with LMOL: Theory, Modeling, and Validation

LMOL, the NASA Langley Mobile Ozone Lidar, is located near NASA’s LaRC steam plant when not deployed in campaigns. The plant produces steam through the incineration of local trash, and its exhaust plume occasionally contains SO2 . SO2 is a known, regulated, pollutant that affects O3 observations in the UV, as is the case with LMOL. In this work, we show how we modified LMOL to detect the plant SO2 plumes and compute its density when the O3 background is stable; we observed densities compatible with what is expected from the typical plume exhaust. The selection of laser lines so that an O3 variation would not get confused with a SO2 detection is explained in details, and a model capable of simulating the Lidar signal with (notably) O3 and SO2 absorption for the validation of the retrieval resolution and uncertainty is presented. Finally, the comparison between the modeled and observed performances of the system is shown: the maximum altitude, resolution, and error in the modeled signal and the observed signals are the same within reasonable margins. This work demonstrates that LMOL is fully capable of working with SO2 . The next step being the addition of an additional laser channel, which is simplified by the tunable aspect of the LMOL laser, to address O3 and SO2 simultaneously, allowing the assessment of SO2 when O3 has variations. Note: This presentation is accompanied by an mp4 video of the given talk.

Guillaume Gronoff↗

NASA NTRS1998

Phase Inversion: Inferring Solar Subphotospheric Flow and Other Asphericity from the Distortion of Acoustic Waves

A method is proposed for analyzing an almost monochromatic train of waves propagating in a single direction in an inhomogeneous medium that is not otherwise changing in time. An effective phase is defined in terms of the Hilbert transform of the wave function, which is related, via the JWKB approximation, to the spatial variation of the background state against which the wave is propagating. The contaminating effect of interference between the truly monochromatic components of the train is eliminated using its propagation properties. Measurement errors, provided they are uncorrelated, are manifest as rapidly varying noise; although that noise can dominate the raw phase-processed signal, it can largely be removed by low-pass filtering. The intended purpose of the analysis is to determine the distortion of solar oscillations induced by horizontal structural variation and material flow. It should be possible to apply the method directly to sectoral modes. The horizontal phase distortion provides a measure of longitudinally averaged properties of the Sun in the vicinity of the equator, averaged also in radius down to the depth to which the modes penetrate. By combining such averages from different modes, the two-dimensional variation can be inferred by standard inversion techniques. After taking due account of horizontal refraction, it should be possible to apply the technique also to locally sectoral modes that propagate obliquely to the equator and thereby build a network of lateral averages at each radius, from which the full three-dimensional structure of the Sun can, in principle, be determined as an inverse Radon transform.

Gough, Douglas↗

NASA NTRS1997

Calibration of a Spatial-Temporal Discrimination Model from Forward, Simultaneous, and Backward Masking

We have been developing a simplified spatial-temporal discrimination model similar to our simplified spatial model in that masking is assumed to be a function of the local visible contrast energy. The overall spatial-temporal sensitivity of the model is calibrated to predict the detectability of targets on a uniform background. To calibrate the spatial-temporal integration functions that define local visible contrast energy, spatial-temporal masking data are required. Observer thresholds were measured (2IFC) for the detection of a 12 msec target stimulus in the presence of a 700 msec mask. Targets were 1, 3 or 9 c/deg sine wave gratings. Masks were either one of these gratings or two of them combined. The target was presented in 17 temporal positions with respect to the mask, including positions before, during and after the mask. Peak masking was found near mask onset and offset for 1 and 3 c/deg targets, while masking effects were more nearly uniform during the mask for the 9 c/deg target. As in the purely spatial case, the simplified model can not predict all the details of masking as a function of masking component spatial frequencies, but overall the prediction errors are small.

Ahumada, Albert J.↗

NASA NTRS1991

Pi/4-QPSK modems for satellite sound/data broadcast systems

The use of pi/4-quadrature phase-shift keying (QPSK) modems for satellite sound broadcast systems (SSBS) broadcasting to mobile or portable receivers is proposed. Three different differential detectors (including the FM-discriminator followed by integrate-sample-and-dump filter) and a novel coherent detector are discussed. The degradation caused by the frequency offset between the local oscillator (LO) and the unmodulated carrier (CR) in the baseband differential detector is studied. The performance of both coherently and differentially detected pi/4-QPSK in a Gaussian channel is also studied. It is shown that with a frequency offset of more than 3 percent of the symbol rate, the performance degradation is more than 1 dB at 0.0001. The out-of-band power of the nonlinearly amplified bandlimited pi/4-QPSK signals is reduced from -13 dB to -37 dB if a 2-dB output back-off amplifier is used instead of a hardlimiter. The performance of the pi/4-QPSK is equivalent to that of QPSK, although the pi/4-QPSK has the advantage of less spectrum restoration after nonlinear amplification. The coherent demodulator and differential decoder avoid the three-level detection and achieve the same bit-error-rate performance as DEQPSK with a simple circuit.

Liu, Chia-Liang↗

NASA NTRS2011

Wake Cycle Robustness of the Mars Science Laboratory Flight Software

The Mars Science Laboratory (MSL) is a spacecraft being developed by the Jet Propulsion Laboratory (JPL) for the purpose of in-situ exploration on the surface of Mars. The objective of MSL is to explore and quantitatively assess a local region on the Martian surface as a habitat for microbial life, past or present. This objective will be accomplished through the assessment of the biological potential of at least one target environment, the characterization of the geology and geochemistry of the landing region, an investigation of the planetary process relevant to past habitability, and a characterization of surface radiation. For this purpose, MSL incorporates a total of ten scientific instruments for which functions are to include, among others, atmospheric and descent imaging, chemical composition analysis, and radiation measurement. The Flight Software (FSW) system is responsible for all mission phases, including launch, cruise, entry-descent-landing, and surface operation of the rover. Because of the essential nature of flight software to project success, each of the software modules is undergoing extensive testing to identify and correct errors.

Whitehill, Robert↗

NASA NTRS1994

Nonvolatile random access memory

A nonvolatile magnetic random access memory can be achieved by an array of magnet-Hall effect (M-H) elements. The storage function is realized with a rectangular thin-film ferromagnetic material having an in-plane, uniaxial anisotropy and inplane bipolar remanent magnetization states. The thin-film magnetic element is magnetized by a local applied field, whose direction is used to form either a 0 or 1 state. The element remains in the 0 or 1 state until a switching field is applied to change its state. The stored information is detcted by a Hall-effect sensor which senses the fringing field from the magnetic storage element. The circuit design for addressing each cell includes transistor switches for providing a current of selected polarity to store a binary digit through a separate conductor overlying the magnetic element of the cell. To read out a stored binary digit, transistor switches are employed to provide a current through a row of Hall-effect sensors connected in series and enabling a differential voltage amplifier connected to all Hall-effect sensors of a column in series. To avoid read-out voltage errors due to shunt currents through resistive loads of the Hall-effect sensors of other cells in the same column, at least one transistor switch is provided between every pair of adjacent cells in every row which are not turned on except in the row of the selected cell.

Wu, Jiin-Chuan↗

NASA NTRS1985

Excess gamma-rays in the direction of the rho Ophiuchi cloud: An exotic object?

The COS-B X-ray data in the direction of the rho Oph dark cloud show an extended structure; at the same time, the region of highest intensity has a spatial distribution compatible with a localized source; 2CG353+16 which is designated Oph gamma. The possibility of an excess gamma ray flux over what is expected on the basis of the interaction of average density cosmic rays with an estimated cloud mass of 2 to 4 000 M is still open, pending an extended CO survey matching the gamma ray data. Estimates for this excess factor are in the range 2 to 4. While the cloud mass may be underestimated, it should be noted that an excess of the same order appears to be present in the nearby Oph-Sag area, well surveyed in CO with the Columbia dish. Possible reasons for a gamma ray excess, in view of two recent observational developments: an Einstein X-ray survey and a VLA radio survey, both covering the approx 2 deg diameter Oph gamma error box. Current interpretations link the gamma ray excess to the cloud gas, in which some active agent is present: stellar winds, or interaction with the North Polar Spur.

Montmerle, T.↗

NASA NTRSDate not supplied

Innovative Features of NASA's Celestial Mapping System to Support Exploration in the Lunar South Pole

Introduction: NASA's Celestial Mapping System (CMS) is developed to address the need for 3D tools for planetary science investigations, mission planning, in-situ operations, in a 3D-first design constructed around a unified view of a planetary globe. At present CMS provides many critical functionalities that include 1) Equipment planning and optimized placement on Lunar surface 2) Line of sight (visibility ) analysis 3) Powerful measurement tools based on 3D terrain with realistic 3D models to represent rovers, astronauts and equipment 4) Visualization of de-rived mapping products (e.g. resource maps), and 5) Data engine for hosting new observations that are not available in other contemporary lunar data tools. CMS is built on the foundation of powerful NASA WorldWind globe engines. In near future, users will be able to simultaneously deploy CMS onto multiple hardware configurations and platforms such as Windows, Linux, iOS and Android. The users will also have the flexibility to update to the latest imagery and terrain datasets as they are being acquired (in real time) before and/or during the exploration mission. CMS is also capable of consumption and analysis of data from locally hosted and external sources. It supports Open Geospatial Consorti-um (OGC) data and file standards, with current integrations of datasets from the Astrogeology Science Center of USGS which include global and local data acquired from NASA (LRO, Clementine, Lunar Orbiter) and JAXA (SELENE/Kaguya) with the capability of integrating more datasets. With development experience in both the end-user application and planetary engine side, CMS is also able to adapt to newer Lunar cartography standards as they develop and become recognized by international geospatial panels. Overcoming Polar Distortions: 3D geospatial applications traditionally suffer from significant distortion of imagery at the poles due to following reasons – 1) distortions in the source imagery 2) Incompatible tessellation algorithm on the poles 3) map projections. In the lunar context, with the focus on the South pole, this is not acceptable. The CMS team is researching ways to address polar distortion of imagery with new tessellation algorithms and by reprojecting the data using projections that are more accurate in polar scenarios. Figure1 shows the potential error introduced by different tessellation methods, represented by the red and green circles for Shoemaker crater. There is ~2 Km difference in the placement of the crater. Line of Sight Analysis and Traverse Planning: We have developed a built-in line of sight analysis (LOS) tool in CMS that analyzes the terrain profile and obstructions and provides the visibility of a given terrain for a remote observer. Figure 2 shows the viewshed analysis on the PSR in Nobile region. The PSR was created with help of HORUS generated images. The yellow pin shows the observer location outside the PSR. The yellow area shows the visible part of PSR. The obstructed area with no visibility for the observer is shown in red. This analysis was ex-tended further to set different heights for various observers and then perform the viewshed analysis. Combining the different visibility profiles can help designing improved traverses within the crater.

Geospatial Mapping↗

NASA NTRSDate not supplied

Innovative Features of NASA's Celestial Mapping System to Support Exploration in the Lunar South Pole

Geospatial Mapping↗

NASA NTRS1998

Polar Hydra Data Analysis

The science activities are: 1) Hydra is still operating successfully on orbit. 2) A large amount of analysis and discovery has occurred with the Hydra ground data processing this past year. 3) Full interdetector calibration has been implemented and documented. This intercalibration was necessitated by the incorrect installation of bias resistors in the pre-acceleration stage to the electron channeltrons. This had the effect of making the counting efficiency for electrons energy dependent as well as channeltron specific. The nature of the error had no impact on the ion detection efficiency since they have a different bias arrangement. This intercalibration is so effective, that the electron and ion moment densities are routinely produced with a level of agreement better than 20%. 4) The data processing routinely removes glint in the sensors and produces public energy time spectrograms on the web overnight. 6) Routine, but more intensive computer processing codes are operational that determine for electrons and ions, the density, the flow vector, the pressure tensor and the heat flux by numerical integration. These codes use the magnetic field to sustain the quality of their output. To gain access to this high quality magnetic field within our data stream we have monitored Russell's web page for zero levels and timing files (since his data acquisition is not telemetry synchronous) and have a local reconstruction of B for our use. We have also detected a routine anomaly in the magnetometer data stream that we have documented to Chris Russell and developed an editing algorithm to intercept these "hits" and remove them from the geophysical analysis.

Scudder, J. D.↗

NASA NTRS2019

NASA Orbital Debris Engineering Model ORDEM 3.1 - Software User Guide

This National Aeronautics and Space Administration (NASA) Orbital Debris Engineering Model (ORDEM) 3.1 Software User Guide accompanies delivery of the latest upgraded version of the model, ORDEM 3.1. The user guide also provides a top-level program description and a list of capabilities. It includes descriptions of runtime error and information codes, input/output file formats, runtimes for different orbit configurations, and how to use uncertainty files. ORDEM 3.1 supersedes the previous NASA Orbital Debris Program Office (ODPO) models – ORDEM 3.0 (Stansbery, et al. 2014) and ORDEM2000 (Liou, et al. 2002). The availability of new sensor and in situ data, re-analysis of older data, and development of new analytical techniques has enabled the construction of this more comprehensive and sophisticated model. An upgraded graphical user interface (GUI) is integrated with the software. This upgraded GUI uses project-oriented organization and provides the user with graphical representations of numerous output data products. For example, these range from the conventional flux vs. average debris size (or altitude bin) for chosen analysis orbits (or views) to the more complex color-contoured, two-dimensional (2-D) directional flux diagrams in local spacecraft elevation and azimuth. The current model, ORDEM 3.1, supports spacecraft as well as telescope/radar project assessments. ORDEM 3.1 contains updated debris populations covering low Earth orbit (LEO, up to 2000 km altitude) to geosynchronous orbit (GEO, up to 40,000 km altitude) and can assess debris calculations up to year 2050, extending coverage past the previous limit of 2035 in ORDEM 3.0. Although populations differ from its predecessor, ORDEM 3.1 is functionally the same as ORDEM 3.0 and can support ORDEM 3.0 projects through backward compatibility.

Vavrin, Andrew B.↗

NASA NTRS2023

NASA Orbital Debris Engineering Model ORDEM 3.2 – Software User Guide

This National Aeronautics and Space Administration (NASA) Orbital Debris Engineering Model (ORDEM) 3.2 Software User Guide accompanies delivery of the latest upgraded version of the model, ORDEM 3.2. The user guide also provides a top-level program description and a list of capabilities. It includes descriptions of runtime error and information codes, input/output file formats, runtimes for different orbit configurations, and how to use uncertainty files. ORDEM 3.2 supersedes the previous NASA Orbital Debris Program Office (ODPO) models – ORDEM 3.0 (Stansbery, et al. 2014) and ORDEM2000 (Liou, et al. 2002). The availability of new sensor and in situ data, re-analysis of older data, and development of new analytical techniques has enabled the construction of this more comprehensive and sophisticated model. An upgraded graphical user interface (GUI) is integrated with the software. This upgraded GUI uses project-oriented organization and provides the user with graphical representations of numerous output data products. For example, these range from the conventional flux vs. average debris size (or altitude bin) for chosen analysis orbits (or views) to the more complex color-contoured, two-dimensional (2-D) directional flux diagrams in local spacecraft elevation and azimuth. The current model, ORDEM 3.2, supports spacecraft as well as telescope/radar project assessments. ORDEM 3.2 contains updated debris populations covering low Earth orbit (LEO, up to 2000 km altitude) to geosynchronous orbit (GEO, up to 40,000 km altitude) and can assess debris calculations up to year 2050, extending coverage past the previous limit of 2035 in ORDEM 3.0. Although populations differ from its predecessor, ORDEM 3.2 is functionally the same as ORDEM 3.0 and can support ORDEM 3.0 projects through backward compatibility.

Andrew Vavrin↗

NASA NTRS2009

Data Assimilation Experiments Using Quality Controlled AIRS Version 5 Temperature Soundings

The AIRS Science Team Version 5 retrieval algorithm has been finalized and is now operational at the Goddard DAAC in the processing (and reprocessing) of all AIRS data. The AIRS Science Team Version 5 retrieval algorithm contains a number of significant improvements over Version 4. Two very significant improvements are described briefly below. 1) The AIRS Science Team Radiative Transfer Algorithm (RTA) has now been upgraded to accurately account for effects of non-local thermodynamic equilibrium on the AIRS observations. This allows for use of AIRS observations in the entire 4.3 micron CO2 absorption band in the retrieval algorithm during both day and night. Following theoretical considerations, tropospheric temperature profile information is obtained almost exclusively from clear column radiances in the 4.3 micron CO2 band in the AIRS Version 5 temperature profile retrieval step. These clear column radiances are a derived product that are indicative of radiances AIRS channels would have seen if the field of view were completely clear. Clear column radiances for all channels are determined using tropospheric sounding 15 micron CO2 observations. This approach allows for the generation of accurate values of clear column radiances and T(p) under most cloud conditions. 2) Another very significant improvement in Version 5 is the ability to generate accurate case-by-case, level-by-level error estimates for the atmospheric temperature profile, as well as for channel-by-channel clear column radiances. These error estimates are used for quality control of the retrieved products. Based on error estimate thresholds, each temperature profiles is assigned a characteristic pressure, pg, down to which the profile is characterized as good for use for data assimilation purposes. We have conducted forecast impact experiments assimilating AIRS quality controlled temperature profiles using the NASA GEOS-5 data assimilation system, consisting of the NCEP GSI analysis coupled with the NASA FVGCM, at a spatial resolution of 0.5 deg by 0.5 deg. Assimilation of Quality Controlled AIRS temperature profiles down to pg resulted in significantly improved forecast skill compared to that obtained from experiments when all data used operationally by NCEP, except for AIRS data, is assimilated. These forecasts were also significantly better than to those obtained when AIRS radiances (rather than temperature profiles) are assimilated, which is the way AIRS data is used operationally by NCEP and ECMWF.

Susskind, Joel↗

NASA NTRS2013

Smart-Divert Powered Descent Guidance to Avoid the Backshell Landing Dispersion Ellipse

A smart-divert capability has been added into the Powered Descent Guidance (PDG) software originally developed for Mars pinpoint and precision landing. The smart-divert algorithm accounts for the landing dispersions of the entry backshell, which separates from the lander vehicle at the end of the parachute descent phase and prior to powered descent. The smart-divert PDG algorithm utilizes the onboard fuel and vehicle thrust vectoring to mitigate landing error in an intelligent way: ensuring that the lander touches down with minimum- fuel usage at the minimum distance from the desired landing location that also avoids impact by the descending backshell. The smart-divert PDG software implements a computationally efficient, convex formulation of the powered-descent guidance problem to provide pinpoint or precision-landing guidance solutions that are fuel-optimal and satisfy physical thrust bound and pointing constraints, as well as position and speed constraints. The initial smart-divert implementation enforced a lateral-divert corridor parallel to the ground velocity vector; this was based on guidance requirements for MSL (Mars Science Laboratory) landings. This initial method was overly conservative since the divert corridor was infinite in the down-range direction despite the backshell landing inside a calculable dispersion ellipse. Basing the divert constraint instead on a local tangent to the backshell dispersion ellipse in the direction of the desired landing site provides a far less conservative constraint. The resulting enhanced smart-divert PDG algorithm avoids impact with the descending backshell and has reduced conservatism.

Carson, John M.↗

NASA NTRS2014

The Seasonal Cycle of Water Vapour on Mars from Assimilation of Thermal Emission Spectrometer Data

We present for the first time an assimilation of Thermal Emission Spectrometer (TES) water vapour column data into a Mars global climate model (MGCM). We discuss the seasonal cycle of water vapour, the processes responsible for the observed water vapour distribution, and the cross-hemispheric water transport. The assimilation scheme is shown to be robust in producing consistent reanalyses, and the global water vapour column error is reduced to around 2-4 pr micron depending on season. Wave activity is shown to play an important role in the water vapour distribution, with topographically steered flows around the Hellas and Argyre basins acting to increase transport in these regions in all seasons. At high northern latitudes, zonal wavenumber 1 and 2 stationary waves during northern summer are responsible for spreading the sublimed water vapour away from the pole. Transport by the zonal wavenumber 2 waves occurs primarily to the west of Tharsis and Arabia Terra and, combined with the effects of western boundary currents, this leads to peak water vapour column abundances here as observed by numerous spacecraft. A net transport of water to the northern hemisphere over the course of one Mars year is calculated, primarily because of the large northwards flux of water vapour which occurs during the local dust storm around L(sub S) = 240-260deg. Finally, outlying frost deposits that surround the north polar cap are shown to be important in creating the peak water vapour column abundances observed during northern summer.

Thermal Emission Spectrometer↗

NASA NTRS1994

Eddy current correction in volume-localized MR spectroscopy

The quality of volume-localized magnetic resonance spectroscopy is affected by eddy currents caused by gradient switching. Eddy currents can be reduced with improved gradient systems; however, it has been suggested that the distortion due to eddy currents can be compensated for during postprocessing with a single-frequency reference signal. The authors propose modifying current techniques for acquiring the single-frequency reference signal by using relaxation weighting to reduce interference from components that cannot be eliminated by digital filtering alone. Additional sequences with T1 or T2 weighting for reference signal acquisition are shown to have the same eddy current characteristics as the original signal without relaxation weighting. The authors also studied a new eddy current correction method that does not require a single-frequency reference signal. This method uses two free induction decays (FIDs) collected from the same volume with two sequences with opposite gradients. Phase errors caused by eddy currents are opposite in these two FIDs and can be canceled completely by combining the FIDs. These methods were tested in a phantom. Eddy current distortions were corrected, allowing quantitative measurement of structures such as the -CH = CH- component, which is otherwise undetectable.

NASA Discipline Musculoskeletal↗

NASA NTRS2016

Monitoring Orbital Precession of EO-1 Hyperion With Three Atmospheric Correction Models in the Libya-4 PICS

Spaceborne spectrometers require spectral-temporal stability characterization to aid in validation of derived data products. Earth Observation 1 (EO-1) began orbital precession in 2011 after exhausting onboard fuel resources. In the Libya-4 pseudo-invariant calibration site (PICS), this resulted in a progressive shift from a mean local equatorial crossing time of approximately10:00 A.M. in 2011 to approximately 8:30 A.M. in late 2015. Here, we studied precession impacts to Hyperion surface reflectance products using three atmospheric correction approaches from 2004 to 2015. Combined difference estimates of surface reflectance were less than 5 percent in the visible near infrared (VNIR) and less than 10 percent for most of the shortwave infrared (SWIR). Combined coefficient of variation estimates in the VNIR ranged from 0.025 to 0.095, and in the SWIR it ranged from 0.025 to 0.06, excluding bands near atmospheric absorption features. Reflectances produced with different atmospheric models were correlated (R squared) in VNIR from 0.25 to 0.94 and in SWIR from 0.12 to 0.88 (p value (calculated probability) less than 0.01). The uncertainties in all the models increased with a terrain slope up to 15 degrees and selecting dune flats could reduce errors. We conclude that these data remain a valuable resource over this period for sensor intercalibration despite orbital decay.

Neigh, Christopher S. R.↗

NASA NTRS2010

Direct Retrieval of Sulfur Dioxide Amount and Altitude from Spaceborne Hyperspectral UV Measurements: Theory and Application

We describe the physical processes by which a vertically localized absorber perturbs the top-of-atmosphere solar backscattered ultraviolet (UV) radiance. The distinct spectral responses to perturbations of an absorber in its column amount and layer altitude provide the basis for a practical satellite retrieval technique, the Extended Iterative Spectral Fitting (EISF) algorithm, for the simultaneous retrieval of these quantities of a SO2 plume. In addition, the EISF retrieval provides an improved UV aerosol index for quantifying the spectral contrast of apparent scene reflectance at the bottom of atmosphere bounded by the surface and/or cloud; hence it can be used for detection of the presence or absence of UV absorbing aerosols. We study the performance and characterize the uncertainties of the EISF algorithm using synthetic backscattered UV radiances, retrievals from which can be compared with those used in the simulation. Our findings indicate that the presence of aerosols (both absorbing and nonabsorbing) does not cause large errors in EISF retrievals under most observing conditions when they are located below the SO2 plume. The EISF retrievals assuming a homogeneous field of view can provide accurate column amounts for inhomogeneous scenes, but they always underestimate the plume altitudes. The EISF algorithm reduces systematic errors present in existing linear retrieval algorithms that use prescribed SO2 plume heights. Applying the EISF algorithm to Ozone Monitoring Instrument satellite observations of the recent Kasatochi volcanic eruption, we demonstrate the successful retrieval of effective plume altitude of volcanic SO2, and we also show the improvement in accuracy in the corresponding SO2 columns.

Yang, Kau↗