Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Functional data analysis”

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.

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

At least 1,297 records · Page 72

Operational Parameters, Considerations, and Design Decisions for Resource-Constrained Ion Trap Mass Spectrometers

Mass spectrometers are increasingly finding applications in new and unique areas, often in situations where key operational resources (i.e. power, weight and size) are limited. One such example is the Mars Organic Molecule Analyzer (MOMA). This instrument is a joint venture between NASA and the European Space Agency (ESA) to develop an ion trap mass spectrometer for chemical analysis on Mars. The constraints on such an instrument are significant as are the performance requirements. While the ideal operating parameters for an ion trap are generally well characterized, methods to maintain analytical performance with limited power and system weight need to be investigated and tested. Methods Experiments have been performed on two custom ion trap mass spectrometers developed as prototypes for the MOMA instrument. This hardware consists of quadrupole ion trap electrodes that are 70% the size of common commercial instrumentation. The trapping RF voltage is created with a custom tank circuit that can be tuned over a range of RF frequencies and is driven using laboratory supplies and amplifiers. The entire instrument is controlled with custom Lab VIEW software that allows a high degree of flexibility in the definition of the scan function defining the ion trap experiment. Ions are typically generated via an internal electron ionization source, however, a laser desorption source is also in development for analysis of larger intact molecules. Preliminary Data The main goals in this work have been to reduce the power required to generate the radio frequency trapping field used in an ion trap mass spectrometer. Generally minimizing the power will also reduce the volume and mass of the electronics to support the instrument. In order to achieve optimum performance, commercial instruments typically utilize RF frequencies in the 1 MHz range. Without much concern for power usage, they simply generate the voltage required to access the mass range of interest. In order to reduce the required RF voltage (and power), operation of the ion trap at lower RF frequencies has been investigated. Surprisingly, the performance of the instrument has only been slightly degraded at RF frequencies all the way down to 500 kHz. Mass resolution is relatively stable to this point and depending on the resonant ejection point used, the peak intensity is also quite stable. To date only masses up to m/z 200 have been fully investigated, however, additional studies are planned to verify the performance with higher mass ions. The lower frequency and voltage should reduce the pseudo potential well depth, eventually affecting the trapping efficiency of the instrument -- effect that could manifest itself in significantly limiting the mass range of trapped ions. Other methods to reduce the RF power while maintaining analytical performance are also under investigation. This includes ion ejection at lower q(sub z) values to access a given mass with a lower RF voltage. The loss of mass resolution at lower q(sub eject) points has been measured and current work is underway to leverage scan speed and the use of non-linear resonances in order to counter this trend. The overall trap performance under this range of operating conditions will be presented with a goal of identifying what trade-offs are acceptable.

Danell, Ryan M.↗

Two-Step Multi-Physics Analysis of an Annular Linear Induction Pump for Fission Power Systems

One of the key technologies associated with fission power systems (FPS) is the annular linear induction pump (ALIP). ALIPs are used to circulate liquid-metal fluid for transporting thermal energy from the nuclear reactor to the power conversion device. ALIPs designed and built to date for FPS project applications have not performed up to expectations. A unique, two-step approach was taken toward the multi-physics examination of an ALIP using ANSYS Maxwell 3D and Fluent. This multi-physics approach was developed so that engineers could investigate design variations that might improve pump performance. Of interest was to determine if simple geometric modifications could be made to the ALIP components with the goal of increasing the Lorentz forces acting on the liquid-metal fluid, which in turn would increase pumping capacity. The multi-physics model first calculates the Lorentz forces acting on the liquid metal fluid in the ALIP annulus. These forces are then used in a computational fluid dynamics simulation as (a) internal boundary conditions and (b) source functions in the momentum equations within the Navier-Stokes equations. The end result of the two-step analysis is a predicted pump pressure rise that can be compared with experimental data.

Fission Power Systems↗

Two-Step Multi-Physics Analysis of an Annular Linear Induction Pump for Fission Power Systems

One of the key technologies associated with fission power systems (FPS) is the annular linear induction pump (ALIP). ALIPs are used to circulate liquid-metal fluid for transporting thermal energy from the nuclear reactor to the power conversion device. ALIPs designed and built to date for FPS project applications have not performed up to expectations. A unique, two-step approach was taken toward the multi-physics examination of an ALIP using ANSYS Maxwell 3D and Fluent. This multi-physics approach was developed so that engineers could investigate design variations that might improve pump performance. Of interest was to determine if simple geometric modifications could be made to the ALIP components with the goal of increasing the Lorentz forces acting on the liquid-metal fluid, which in turn would increase pumping capacity. The multi-physics model first calculates the Lorentz forces acting on the liquid metal fluid in the ALIP annulus. These forces are then used in a computational fluid dynamics simulation as (a) internal boundary conditions and (b) source functions in the momentum equations within the Navier-Stokes equations. The end result of the two-step analysis is a predicted pump pressure rise that can be compared with experimental data.

Fission Power Systems↗

Tropical Simultaneous Nadir Observations for IR Sounder Evaluation and Comparison

Simultaneous Nadir Observations (SNOs) contain thousands of pairs of observations taken within 8 km and 10 minutes. SNOs have been very useful for comparisons in polar conditions. But classic SNO pairing criteria have very low yield in the tropics, making these SNOs less useful for investigating instrument performance for hotter scenes or scenes with high contrast. We introduce a modified methodology, which finds pairs of matched spectra in a wider angular range but is restricted to the tropics. We illustrate this with AIRS and CrIS data. Insight into instrument differences is gained from statistical distributions of the residual differences between the matched pairs. A sample analysis compares AIRS and CrIS brightness temperature at 900 cm-1 as function of scene brightness temperature. The proposed method may be applicable to matchups of other sensors on different spacecraft.

AIRS↗

Adaptations of Endolithic Communities to Abrupt Environmental Changes in A Hyper-Arid Desert

The adaptation mechanisms of microbial communities to natural perturbations remain relatively unexplored, particularly in extreme environments. The extremophilic communities of halite (NaCl) nodules from the hyper-arid core of the Atacama Desert are self-sustained and represent a unique opportunity to study functional adaptations and community dynamics with changing environmental conditions. We transplanted halite nodules to different sites in the desert and investigated how their taxonomic, cellular, and biochemical changes correlated with water availability, using environmental data modeling and metagenomic analyses. Salt-in strategists, mainly represented by haloarchaea, significantly increased in relative abundance at sites characterized by extreme dryness, multiple wet/dry cycles, and colder conditions. The functional analysis of metagenome-assembled genomes (MAGs) revealed site-specific enrichments in archaeal MAGs encoding for the uptake of various compatible solutes and for glycerol utilization. These findings suggest that opportunistic salt-in strategists took over the halite communities at the driest sites. They most likely benefited from metabolites newly released in the environment by the death of microorganisms least adapted to the new conditions. The observed changes were consistent with the need to maximize cellular bioenergetics when confronted with lower water availability and higher salinity, providing valuable information on microbial community adaptations and resilience to climate change.

Atacama↗

Aeroservoelastic Uncertainty Model Identification from Flight Data

Uncertainty modeling is a critical element in the estimation of robust stability margins for stability boundary prediction and robust flight control system development. There has been a serious deficiency to date in aeroservoelastic data analysis with attention to uncertainty modeling. Uncertainty can be estimated from flight data using both parametric and nonparametric identification techniques. The model validation problem addressed in this paper is to identify aeroservoelastic models with associated uncertainty structures from a limited amount of controlled excitation inputs over an extensive flight envelope. The challenge to this problem is to update analytical models from flight data estimates while also deriving non-conservative uncertainty descriptions consistent with the flight data. Multisine control surface command inputs and control system feedbacks are used as signals in a wavelet-based modal parameter estimation procedure for model updates. Transfer function estimates are incorporated in a robust minimax estimation scheme to get input-output parameters and error bounds consistent with the data and model structure. Uncertainty estimates derived from the data in this manner provide an appropriate and relevant representation for model development and robust stability analysis. This model-plus-uncertainty identification procedure is applied to aeroservoelastic flight data from the NASA Dryden Flight Research Center F-18 Systems Research Aircraft.

Brenner, Martin J.↗

Error analysis for reducing noisy wide-gap concentric cylinder rheometric data for nonlinear fluids - Theory and applications

This work discusses the propagation of errors for the recovery of the shear rate from wide-gap concentric cylinder viscometric measurements of non-Newtonian fluids. A least-square regression of stress on angular velocity data to a system of arbitrary functions is used to propagate the errors for the series solution to the viscometric flow developed by Krieger and Elrod (1953) and Pawlowski (1953) ('power-law' approximation) and for the first term of the series developed by Krieger (1968). A numerical experiment shows that, for measurements affected by significant errors, the first term of the Krieger-Elrod-Pawlowski series ('infinite radius' approximation) and the power-law approximation may recover the shear rate with equal accuracy as the full Krieger-Elrod-Pawlowski solution. An experiment on a clay slurry indicates that the clay has a larger yield stress at rest than during shearing, and that, for the range of shear rates investigated, a four-parameter constitutive equation approximates reasonably well its rheology. The error analysis presented is useful for studying the rheology of fluids such as particle suspensions, slurries, foams, and magma.

Borgia, Andrea↗

Optimal spatial filtering and transfer function for SAR ocean wave spectra

The Seasat Synthetic Aperture Radar (SAR) has proved to be an instrument of great utility in the sensing of ocean conditions on a global scale. An analysis of oceanographic and atmospheric aspects of Seasat data has shown that the features observed in the imagery are linked to ocean phenomena such as storm sources and their resulting swell systems. However, there remains one central problem which has not been satisfactorily solved to date. This problem is related to the accurate measurement of wind-generated ocean wave spectra. Investigations addressing this problem are currently being conducted. The problem has two parts, including the accurate measurement of the image spectra and the inference of actual surface wave spectra from these measurements. A description is presented of the progress made towards solving the first part of the problem, taking into account a digital rather than optical computation of the image transforms.

Goldfinger, A. D.↗

Analysis and retracking of continental ice sheet radar altimeter waveforms

The Seasat-1 radar altimeter data set acquired over both the Antarctic and Greenland continental ice sheets is analyzed to obtain corrected ranges to the ice surface. The radar altimeter functional response over the continental ice sheets is considerably more complex than over the oceans. Causal factors identified in this complicated response include sloping surfaces, undulating ice surfaces with characteristic wavelengths on the same spatial scale as the altimeter beam-limited footprint, off-track reflections, and dynamic lag of the altimeter tracking circuit. Retracking methods using the altimeter return pulse waveforms give range corrections that are typically several meters. The entire set of Seasat-1 altimetry over the continental ice sheets is being retracked by fitting a multi-parameter function to each waveform. Many waveforms have double ramps indicating near-normal reflections from two distinct portions of the ice surface within the altimeter beam. Two independent range measurements differing by less than 25 m are obtained from retracking the double-ramp waveforms.

Martin, T. V.↗

Waste water processing technology for Space Station Freedom - Comparative test data analysis

Comparative tests were conducted to choose the optimum technology for waste water processing on SSF. A thermoelectric integrated membrane evaporation (TIMES) subsystem and a vapor compression distillation subsystem (VCD) were built and tested to compare urine processing capability. Water quality, performance, and specific energy were compared for conceptual designs intended to function as part of the water recovery and management system of SSF. The VCD is considered the most mature and efficient technology and was selected to replace the TIMES as the baseline urine processor for SSF.

Miernik, Janie H.↗

Ascent Heating Thermal Analysis on Spacecraft Adaptor Fairings

When the Crew Exploration Vehicle (CEV) is launched, the spacecraft adaptor (SA) fairings that cover the CEV service module (SM) are exposed to aero heating. Thermal analysis is performed to compute the fairing temperatures and to investigate whether the temperatures are within the material limits for nominal ascent aeroheating case. The ascent heating is analyzed by using computational fluid dynamics (CFD) and engineering codes at Marshall Space Flight Center. The aeroheating environment data used for this work is known as Thermal Environment 3 (TE3) heating data. One of the major concerns is with the SA fairings covering the CEV SM and the SM/crew launch vehicle (CLV) flange interface. The TE3 heating rate is a function of time, wall temperature, and the spatial locations. The implementation of the TE3 heating rate as boundary conditions in the thermal analysis becomes challenging. The ascent heating thermal analysis on SA fairings and SM/CLV flange interface are performed using two commercial software packages: Cullimore & Ring (C&R) Thermal Desktop (TD) 5.1 and MSC Patran 2007r1 b. TD is the pre-and post-processor for SINDA, which is a finite-difference-based solver. In TD, the geometry is built and meshed, the boundary conditions are defined, and then SINDA is used to compute temperatures. MSC Pthermal is a finite-element- based thermal solver. MSC Patran is the pre- and post-processor for Pthermal. Regarding the boundary conditions, the convection, contact resistance, and heat load can be imposed in different ways in both programs. These two software packages are used to build the thermal model for the same analysis to validate each other and show the differences in the modeling details.

Wang, Xiao Yen↗

Simultaneous measurement of stratospheric O3, H2O, CH4, and N2O profiles from infrared limb thermal emissions

Thermal emission measurements of the earth's stratospheric limb were made with a cryogenically cooled high-resolution Michelson interferometer on a balloon flight launched from Palestine, TX, on Nov. 6, 1984. Infrared spectra for complete limb sequences were obtained over portions of the 700-1940/cm range with an unapodized spectral resolution of 0.03/cm for tangent heights varying from 13 to 39 km. The observed data from 1125 to 1425/cm have been analyzed for simultaneous measurement of O3, H2O, CH4, and N2O profiles. The analysis employs line-by-line and layer-by-layer radiative-transfer calculations, including curvature and refraction effects. The optimum use of geometric and spectral effects is made to obtain sharply peaked weighting functions. Contributions from stratospheric aerosol are included by measuring the light extinction within the window regions of the observed spectra. The retrieved constituent profiles are compared with measurements made with a variety of techniques by other groups. The comparison shows good agreement with the published data for all gases, indicating the capability of retrieving trace gas profiles from high-resolution thermal emission limb measurements.

Abbas, M. M.↗

File Specification for GEOS-5 FP (Forward Processing)

The GEOS-5 FP Atmospheric Data Assimilation System (GEOS-5 ADAS) uses an analysis developed jointly with NOAA's National Centers for Environmental Prediction (NCEP), which allows the Global Modeling and Assimilation Office (GMAO) to take advantage of the developments at NCEP and the Joint Center for Satellite Data Assimilation (JCSDA). The GEOS-5 AGCM uses the finite-volume dynamics (Lin, 2004) integrated with various physics packages (e.g, Bacmeister et al., 2006), under the Earth System Modeling Framework (ESMF) including the Catchment Land Surface Model (CLSM) (e.g., Koster et al., 2000). The GSI analysis is a three-dimensional variational (3DVar) analysis applied in grid-point space to facilitate the implementation of anisotropic, inhomogeneous covariances (e.g., Wu et al., 2002; Derber et al., 2003). The GSI implementation for GEOS-5 FP incorporates a set of recursive filters that produce approximately Gaussian smoothing kernels and isotropic correlation functions. The GEOS-5 ADAS is documented in Rienecker et al. (2008). More recent updates to the model are presented in Molod et al. (2011). The GEOS-5 system actively assimilates roughly 2 × 10(exp 6) observations for each analysis, including about 7.5 × 10(exp 5) AIRS radiance data. The input stream is roughly twice this volume, but because of the large volume, the data are thinned commensurate with the analysis grid to reduce the computational burden. Data are also rejected from the analysis through quality control procedures designed to detect, for example, the presence of cloud. To minimize the spurious periodic perturbations of the analysis, GEOS-5 FP uses the Incremental Analysis Update (IAU) technique developed by Bloom et al. (1996). More details of this procedure are given in Appendix A. The assimilation is performed at a horizontal resolution of 0.3125-degree longitude by 0.25- degree latitude and at 72 levels, extending to 0.01 hPa. All products are generated at the native resolution of the horizontal grid. The majority of data products are time-averaged, but four instantaneous products are also available. Hourly data intervals are used for two-dimensional products, while 3-hourly intervals are used for three-dimensional products. These may be on the model's native 72-layer vertical grid or at 42 pressure surfaces extending to 0.1 hPa. This document describes the gridded output files produced by the GMAO near real-time operational FP, using the most recent version of the GEOS-5 assimilation system. Additional details about variables listed in this file specification can be found in a separate document, the GEOS-5 File Specification Variable Definition Glossary. Documentation about the current access methods for products described in this document can be found on the GMAO products page: http://gmao.gsfc.nasa.gov/products/.

GSI↗

Local Elastic Constants for Epoxy-Nanotube Composites from Molecular Dynamics Simulation

A method from molecular dynamics simulation is developed for determining local elastic constants of an epoxy/nanotube composite. The local values of C11, C33, K12, and K13 elastic constants are calculated for an epoxy/nanotube composite as a function of radial distance from the nanotube. While the results possess a significant amount of statistical uncertainty resulting from both the numerical analysis and the molecular fluctuations during the simulation, the following observations can be made. If the size of the region around the nanotube is increased from shells of 1 to 6 in thickness, then the scatter in the data reduces enough to observe trends. All the elastic constants determined are at a minimum 20 from the center of the nanotube. The C11, C33, and K12 follow similar trends as a function of radial distance from the nanotube. The K13 decreases greater distances from the nanotube and becomes negative which may be a symptom of the statistical averaging.

Frankland, S. J. V.↗

Time and Measurement Days

Questions in data analysis involving the concepts of time and measurement are often pushed into the background or reserved for a philosophical discussion. Some examples are: a) Is causality a consequence of the laws of physics, or can the arrow of time be reversed? b) Can we determine the arrow of time of an event? c) Do we need the continuum hypothesis for the underlying function in any measurement process? d) Can we say anything about the analyticity of the underlying process of an event? e) Would it be valid to model a non-analytical process as function of time? f) What are the implications of all these questions for classical Fourier techniques? However, in the age of big data gathered either from space missions supplying ultra-precise long time series, or e.g. LIGO data from the ground, the moment to bring these questions to the foreground seems arrived. The limitations of our understanding of some fundamental processes is emphasized by the lack of solution for problems open for more than 2 decades, such as the non-detection of solar g-modes, or the modal identification of main sequence stellar pulsators like delta Scuti stars. Flicker noise or 1/f noise, for example, attributed in solar-like stars to granulation, is analyzed mostly only to apply noise reduction techniques, neither considering the classical problem of 1/f noise that was introduced a 100 years ago, nor taking into account ergodic or non-ergodic solutions that make inapplicable spectral analysis techniques in practice. This topic was discussed by Nicholas W. Watkins during the ITISE meeting held in Granada in 2016. There he presented preliminary results of his research on Mandelbrot's related work. We reproduce here his quotation of Mandelbrot (1999) "There is a sharp contrast between a highly anomalous ("non-white") noise that proceeds in ordinary clock time and a noise whose principal anomaly is that it is restricted to fractal time", suggesting a connection with the above proposed topics that could be phrased as the following additional questions:a) Is self-organized criticality (SOC) frequent in astrophysical phenomena? b) Could all fractals in nature be considered stochastic? c) Could we establish mathematical/physical relationships between chaotic and fractal behaviors in time series? d) Could the differences between fractals and chaos in terms of analyticity be used to understand the residuals of the fitting of stellar light curves? In this meeting we would like to approximate these problems from a holistic and multidisciplinary perspective, taking into account not only technical issues but also the deeper implications. In particular the concept of connectivity (introduced in Pascual-Granado et al. A&A, 2015) could be used to implement, within the framework of ARMA processes, an "arrow of time" (see attached document), and so studying the possible implications in the concept of time as envisaged by Watkins.

data analysis↗

Combining Data with Physical Knowledge for Uncertainty Quantification in Certification and Reliability Analysis

Unifying empirical data with predictive models can enable engineering cost-savings through certification by analysis and reliability-based design. Both concepts require rigorous uncertainty quantification (UQ) and robust understanding and treatment of relevant physics. Combining sampling-based UQ algorithms with high-fidelity simulations creates a computational bottleneck that is often alleviated through the use of machine learning (ML). ML can be used to create computationally efficient surrogates for simulations of complex or high-dimensional physical interactions (e.g., multi-phase interactions associated with melt pools in laser powder bed fusion or spatially-dependent material properties in functionally graded materials). However, negative side effects of ML may include a lack of interpretability and negative correlation between event rarity and simulation accuracy due to a lack of training data. As such, it is important to infuse ML algorithms with physics-based guardrails to provide confidence in their predictions. This talk will provide a brief review of recent NASA research at this intersection of physics-based simulation, ML, and UQ with a focus on certification and reliability analysis.

uncertainty quantification↗

Linking Improvements in Sulfur Dioxide Emissions to Decreasing Sulfate Wet Deposition by Combining Satellite and Surface Observations with Trajectory Analysis

Sulfur dioxide (SO2), a criteria pollutant, and sulfate (SO42) deposition are major environmental concerns in the eastern U.S. and both have been on the decline for two decades. In this study, we use satellite column SO2 data from the Ozone Monitoring Instrument (OMI), and SO42 wet deposition data from the NADP (National Atmospheric Deposition Program) to investigate the temporal and spatial relationship between trends in SO2 emissions and the downward sulfate wet deposition over the eastern U.S. from 2005 to 2015. To establish the relationship between SO2 emission sources and receptor sites, we conducted a Potential Source Contribution Function (PSCF) analysis using HYSPLIT back trajectories for five selected Air Quality System (AQS) sites - (Hackney, OH, Akron, OH, South Fayette, PA, Wilmington, DE, and Beltsville, MD) - in close proximity to NADP sites with large downward SO42 trends since 2005. Back trajectories were run for three summers (JJA) and three winters (DJF) and used to generate seasonal climatology PSCFs for each site. The OMI SO2 and interpolated NADP sulfate deposition trends were normalized and overlapped with the PSCF, to identify the areas that had the highest contribution to the observed drop. The results suggest that emission reductions along the Ohio River Valley have led to decreases in sulfate deposition in eastern OH and western PA (Hackney, Akron and South Fayette). Farther to the east, emission reductions in southeast PA resulted in improvements in sulfate deposition at Wilmington, DE, while for Beltsville, reductions in both the Ohio River Valley and nearby favorably impacted sulfate deposition. For Beltsville, sources closer than 300km from the site contribute roughly 56% observed deposition trends in winter, and 82% in summer, reflecting seasonal changes in transport pattern as well as faster oxidation and washout of sulfur in summer. This suggests that emissions and wet deposition are linked through not only the location of sources relative to the observing sites, but also to photochemistry and the weather patterns characteristic to the region, as evidenced by a west to east shift in the contribution between winter and summer. The method developed here is applicable to other regions with significant trends such as China and India, and can be used to estimate the potential benefits of emission reduction in those areas.

Acid depostion↗

Linking Improvements in Sulfur Dioxide Emissions to Decreasing Sulfate Wet Deposition By Combining Satellite and Surface Observations With Trajectory Analysis

Sulfur dioxide (SO2), a criteria pollutant, and sulfate (SO42−) deposition are major environmental concerns in the eastern U.S. and both have been on the decline for two decades. In this study, we use satellite column SO2 data from the Ozone Monitoring Instrument (OMI), and SO42− wet deposition data from the NADP (National Atmospheric Deposition Program) to investigate the temporal and spatial relationship between trends in SO2 emissions and the downward sulfate wet deposition over the eastern U.S. from 2005 to 2015. To establish the relationship between SO2 emission sources and receptor sites, we conducted a Potential Source Contribution Function (PSCF) analysis using HYSPLIT back trajectories for five selected Air Quality System (AQS) sites - (Hackney, OH, Akron, OH, South Fayette, PA, Wilmington, DE, and Beltsville, MD) - in close proximity to NADP sites with large downward SO42− trends since 2005. Back trajectories were run for three summers (JJA) and three winters (DJF) and used to generate seasonal climatology PSCFs for each site. The OMI SO2 and interpolated NADP sulfate deposition trends were normalized and overlapped with the PSCF, to identify the areas that had the highest contribution to the observed drop. The results suggest that emission reductions along the Ohio River Valley have led to decreases in sulfate deposition in eastern OH and western PA (Hackney, Akron and South Fayette). Farther to the east, emission reductions in southeast PA resulted in improvements in sulfate deposition at Wilmington, DE, while for Beltsville, reductions in both the Ohio River Valley and nearby favorably impacted sulfate deposition. For Beltsville, sources closer than 300 km from the site contribute roughly 56% observed deposition trends in winter, and 82% in summer, reflecting seasonal changes in transport pattern as well as faster oxidation and washout of sulfur in summer. This suggests that emissions and wet deposition are linked through not only the location of sources relative to the observing sites, but also to photochemistry and the weather patterns characteristic to the region, as evidenced by a west to east shift in the contribution between winter and summer. The method developed here is applicable to other regions with significant trends such as China and India, and can be used to estimate the potential benefits of emission reduction in those areas.

sulfur dioxide emissions↗