Engineering topics

Wang, Y.

Publications and source records attributed to Wang, Y..

At least 37 records · Page 2

NASA NTRS2014

Observation of Mountain Lee Waves with MODIS NIR Column Water Vapor

Mountain lee waves have been previously observed in data from the Moderate Resolution Imaging Spectroradiometer (MODIS) "water vapor" 6.7 micrometers channel which has a typical peak sensitivity at 550 hPa in the free troposphere. This paper reports the first observation of mountain waves generated by the Appalachian Mountains in the MODIS total column water vapor (CWV) product derived from near-infrared (NIR) (0.94 micrometers) measurements, which indicate perturbations very close to the surface. The CWV waves are usually observed during spring and late fall or some summer days with low to moderate CWV (below is approx. 2 cm). The observed lee waves display wavelengths from3-4 to 15kmwith an amplitude of variation often comparable to is approx. 50-70% of the total CWV. Since the bulk of atmospheric water vapor is confined to the boundary layer, this indicates that the impact of thesewaves extends deep into the boundary layer, and these may be the lowest level signatures of mountain lee waves presently detected by remote sensing over the land.

satellite remote sensing↗

NASA NTRS2013

High Resolution Aerosol Data from MODIS Satellite for Urban Air Quality Studies

The Moderate Resolution Imaging Spectroradiometer (MODIS) provides daily global coverage, but the 10 km resolution of its aerosol optical depth (AOD) product is not suitable for studying spatial variability of aerosols in urban areas. Recently, a new Multi-Angle Implementation of Atmospheric Correction (MAIAC) algorithm was developed for MODIS which provides AOD at 1 km resolution. Using MAIAC data, the relationship between MAIAC AOD and PM(sub 2.5) as measured by the 27 EPA ground monitoring stations was investigated. These results were also compared to conventional MODIS 10 km AOD retrievals (MOD04) for the same days and locations. The coefficients of determination for MOD04 and for MAIAC are R(exp 2) =0.45 and 0.50 respectively, suggested that AOD is a reasonably good proxy for PM(sub 2.5) ground concentrations. Finally, we studied the relationship between PM(sub 2.5) and AOD at the intra-urban scale (10 km) in Boston. The fine resolution results indicated spatial variability in particle concentration at a sub-10 kilometer scale. A local analysis for the Boston area showed that the AOD-PM(sub 2.5) relationship does not depend on relative humidity and air temperatures below approximately 7 C. The correlation improves for temperatures above 7 - 16 C. We found no dependence on the boundary layer height except when the former was in the range 250-500 m. Finally, we apply a mixed effects model approach to MAIAC aerosol optical depth (AOD) retrievals from MODIS to predict PM(sub 2.5) concentrations within the greater Boston area. With this approach we can control for the inherent day-to-day variability in the AOD-PM(sub 2.5) relationship, which depends on time-varying parameters such as particle optical properties, vertical and diurnal concentration profiles and ground surface reflectance. Our results show that the model-predicted PM(sub 2.5) mass concentrations are highly correlated with the actual observations (out-of-sample R(exp 2) of 0.86). Therefore, adjustment for the daily variability in the AOD-PM(sub 2.5) relationship provides a means for obtaining spatially-resolved PM(sub 2.5) concentrations.

remote sensing↗

NASA NTRS2013

10 Yr Spatial and Temporal Trends of PM2.5 Concentrations in the Southeastern US Estimated Using High-resolution Satellite Data

Long-term PM2.5 exposure has been reported to be associated with various adverse health outcomes. However, most ground monitors are located in urban areas, leading to a potentially biased representation of the true regional PM2.5 levels. To facilitate epidemiological studies, accurate estimates of spatiotemporally continuous distribution of PM2.5 concentrations are essential. Satellite-retrieved aerosol optical depth (AOD) has been widely used for PM2.5 concentration estimation due to its comprehensive spatial coverage. Nevertheless, an inherent disadvantage of current AOD products is their coarse spatial resolutions. For instance, the spatial resolutions of the Moderate Resolution Imaging Spectroradiometer (MODIS) and the Multiangle Imaging SpectroRadiometer (MISR) are 10 km and 17.6 km, respectively. In this paper, a new AOD product with 1 km spatial resolution retrieved by the multi-angle implementation of atmospheric correction (MAIAC) algorithm was used. A two-stage model was developed to account for both spatial and temporal variability in the PM2.5-AOD relationship by incorporating the MAIAC AOD, meteorological fields, and land use variables as predictors. Our study area is in the southeastern US, centered at the Atlanta Metro area, and data from 2001 to 2010 were collected from various sources. The model was fitted for each year individually, and we obtained model fitting R2 ranging from 0.71 to 0.85, MPE from 1.73 to 2.50 g m3, and RMSPE from 2.75 to 4.10 g m3. In addition, we found cross validation R2 ranging from 0.62 to 0.78, MPE from 2.00 to 3.01 g m3, and RMSPE from 3.12 to 5.00 g m3, indicating a good agreement between the estimated and observed values. Spatial trends show that high PM2.5 levels occurred in urban areas and along major highways, while low concentrations appeared in rural or mountainous areas. A time series analysis was conducted to examine temporal trends of PM2.5 concentrations in the study area from 2001 to 2010. The results showed that the PM2.5 levels in the study area followed a generally declining trend from 2001 to 2010 and decreased about 20 during the period. However, there was an exception of an increase in year 2005, which is attributed to elevated sulfate concentrations in the study area in warm months of 2005. An investigation of the impact of wild and prescribed fires on PM2.5 levels in 2007 suggests a positive relationship between them.

particulate matter (PM2.5)↗

NASA NTRS2012

Multi-Angle Implementation of Atmospheric Correction (MAIAC) Algorithm

Multi-Angle Implementation of Atmospheric Correction (MAIAC) is a new algorithm developed for MODIS. MAIAC uses a time series analysis and processing of groups of pixels to perform simultaneous retrievals of aerosol properties and surface bidirectional reflectance without typical assumptions about the surface. It is a generic algorithm which works over both dark and bright land surfaces, including deserts. MAIAC has an internal Cloud Mask, a dynamic land-water-snow classification and a surface change mask which allows it to flexibly choose processing path over different surfaces. A distinct feature of MAIAC is a high 1 km resolution of aerosol retrievals which is required in different applications including the air quality analysis. The novel features of MAIAC include the high quality cloud mask, discrimination of aerosol type, including biomass burning smoke and dust, and detection of surface change - all required for high quality aerosol retrievals. An overview of the algorithm, results of AERONET validation, and examples of comparison with MODIS Collection 5 aerosol product and Deep Blue algorithm for different parts of the world, will be presented.

Lyapustin, A.↗

NASA NTRS2012

The First in situ Observation of Kelvin-Helmholtz Waves at High-Latitude Magnetopause during Strongly Dawnward Interplanetary Magnetic Field Conditions

We report the first in situ observation of high-latitude magnetopause (near the northern duskward cusp) Kelvin-Helmholtz waves (KHW) by Cluster on January 12, 2003, under strongly dawnward interplanetary magnetic field (IMF) conditions. The fluctuations unstable to Kelvin-Helmholtz instability (KHI) are found to propagate mostly tailward, i.e., along the direction almost 90 deg. to both the magnetosheath and geomagnetic fields, which lowers the threshold of the KHI. The magnetic configuration across the boundary layer near the northern duskward cusp region during dawnward IMF is similar to that in the low-latitude boundary layer under northward IMF, in that (1) both magnetosheath and magnetospheric fields across the local boundary layer constitute the lowest magnetic shear and (2) the tailward propagation of the KHW is perpendicular to both fields. Approximately 3-hour-long periods of the KHW during dawnward IMF are followed by the rapid expansion of the dayside magnetosphere associated with the passage of an IMF discontinuity that characterizes an abrupt change in IMF cone angle, Phi = acos (B(sub x) / absolute value of Beta), from approx. 90 to approx. 10. Cluster, which was on its outbound trajectory, continued observing the boundary waves at the northern evening-side magnetopause during sunward IMF conditions following the passage of the IMF discontinuity. By comparing the signatures of boundary fluctuations before and after the IMF discontinuity, we report that the frequencies of the most unstable KH modes increased after the discontinuity passed. This result demonstrates that differences in IMF orientations (especially in f) are associated with the properties of KHW at the high-latitude magnetopause due to variations in thickness of the boundary layer, and/or width of the KH-unstable band on the surface of the dayside magnetopause.

solar wind↗

NASA NTRS2012

Discrimination of Biomass Burning Smoke and Clouds in MAIAC Algorithm

The multi-angle implementation of atmospheric correction (MAIAC) algorithm makes aerosol retrievals from MODIS data at 1 km resolution providing information about the fine scale aerosol variability. This information is required in different applications such as urban air quality analysis, aerosol source identification etc. The quality of high resolution aerosol data is directly linked to the quality of cloud mask, in particular detection of small (sub-pixel) and low clouds. This work continues research in this direction, describing a technique to detect small clouds and introducing the smoke test to discriminate the biomass burning smoke from the clouds. The smoke test relies on a relative increase of aerosol absorption at MODIS wavelength 0.412 micrometers as compared to 0.47-0.67 micrometers due to multiple scattering and enhanced absorption by organic carbon released during combustion. This general principle has been successfully used in the OMI detection of absorbing aerosols based on UV measurements. This paper provides the algorithm detail and illustrates its performance on two examples of wildfires in US Pacific North-West and in Georgia/Florida of 2007.

Lyapustin, A.↗

NASA NTRS2012

Improved Cloud and Snow Screening in MAIAC Aerosol Retrievals Using Spectral and Spatial Analysis

An improved cloud/snow screening technique in the Multi-Angle Implementation of Atmospheric Correction (MAIAC) algorithm is described. It is implemented as part of MAIAC aerosol retrievals based on analysis of spectral residuals and spatial variability. Comparisons with AERONET aerosol observations and a large-scale MODIS data analysis show strong suppression of aerosol optical thickness outliers due to unresolved clouds and snow. At the same time, the developed filter does not reduce the aerosol retrieval capability at high 1 km resolution in strongly inhomogeneous environments, such as near centers of the active fires. Despite significant improvement, the optical depth outliers in high spatial resolution data are and will remain the problem to be addressed by the application-dependent specialized filtering techniques.

Lyapustin, A.↗

NASA NTRS2012

Analysis of Satellite-Derived Arctic Tropospheric BrO Columns in Conjunction with Aircraft Measurements During ARCTAS and ARCPAC

We derive tropospheric column BrO during the ARCTAS and ARCPAC field campaigns in spring 2008 using retrievals of total column BrO from the satellite UV nadir sensors OMI and GOME-2 using a radiative transfer model and stratospheric column BrO from a photochemical simulation. We conduct a comprehensive comparison of satellite-derived tropospheric BrO column to aircraft in-situ observations ofBrO and related species. The aircraft profiles reveal that tropospheric BrO, when present during April 2008, was distributed over a broad range of altitudes rather than being confined to the planetary boundary layer (PBL). Perturbations to the total column resulting from tropospheric BrO are the same magnitude as perturbations due to longitudinal variations in the stratospheric component, so proper accounting of the stratospheric signal is essential for accurate determination of satellite-derived tropospheric BrO. We find reasonably good agreement between satellite-derived tropospheric BrO and columns found using aircraft in-situ BrO profiles, particularly when satellite radiances were obtained over bright surfaces (albedo> 0.7), for solar zenith angle < 80 and clear sky conditions. The rapid activation of BrO due to surface processes (the bromine explosion) is apparent in both the OMI and GOME-2 based tropospheric columns. The wide orbital swath of OMI allows examination of the evolution of tropospheric BrO on about hourly time intervals near the pole. Low surface pressure, strong wind, and high PBL height are associated with an observed BrO activation event, supporting the notion of bromine activation by high winds over snow.

Choi, S.↗

NASA NTRS2012

Multi-Angle Implementation of Atmospheric Correction for MODIS (MAIAC). Part 3: Atmospheric Correction

This paper describes the atmospheric correction (AC) component of the Multi-Angle Implementation of Atmospheric Correction algorithm (MAIAC) which introduces a new way to compute parameters of the Ross-Thick Li-Sparse (RTLS) Bi-directional reflectance distribution function (BRDF), spectral surface albedo and bidirectional reflectance factors (BRF) from satellite measurements obtained by the Moderate Resolution Imaging Spectroradiometer (MODIS). MAIAC uses a time series and spatial analysis for cloud detection, aerosol retrievals and atmospheric correction. It implements a moving window of up to 16 days of MODIS data gridded to 1 km resolution in a selected projection. The RTLS parameters are computed directly by fitting the cloud-free MODIS top of atmosphere (TOA) reflectance data stored in the processing queue. The RTLS retrieval is applied when the land surface is stable or changes slowly. In case of rapid or large magnitude change (as for instance caused by disturbance), MAIAC follows the MODIS operational BRDF/albedo algorithm and uses a scaling approach where the BRDF shape is assumed stable but its magnitude is adjusted based on the latest single measurement. To assess the stability of the surface, MAIAC features a change detection algorithm which analyzes relative change of reflectance in the Red and NIR bands during the accumulation period. To adjust for the reflectance variability with the sun-observer geometry and allow comparison among different days (view geometries), the BRFs are normalized to the fixed view geometry using the RTLS model. An empirical analysis of MODIS data suggests that the RTLS inversion remains robust when the relative change of geometry-normalized reflectance stays below 15%. This first of two papers introduces the algorithm, a second, companion paper illustrates its potential by analyzing MODIS data over a tropical rainforest and assessing errors and uncertainties of MAIAC compared to conventional MODIS products.

Lyapustin, A.↗

NASA NTRS2011

Comparative Analysis of Aerosol Retrievals from MODIS, OMI and MISR Over Sahara Region

MODIS is a wide field-of-view sensor providing daily global observations of the Earth. Currently, global MODIS aerosol retrievals over land are performed with the main Dark Target algorithm complimented with the Deep Blue (DB) Algorithm over bright deserts. The Dark Target algorithm relies on surface parameterization which relates reflectance in MODIS visible bands with the 2.1 micrometer region, whereas the Deep Blue algorithm uses an ancillary angular distribution model of surface reflectance developed from the time series of clear-sky MODIS observations. Recently, a new Multi-Angle Implementation of Atmospheric Correction (MAIAC) algorithm has been developed for MODIS. MAIAC uses a time series and an image based processing to perform simultaneous retrievals of aerosol properties and surface bidirectional reflectance. It is a generic algorithm which works over both dark vegetative surfaces and bright deserts and performs retrievals at 1 km resolution. In this work, we will provide a comparative analysis of DB, MAIAC, MISR and OMI aerosol products over bright deserts of northern Africa.

Lyapustin, A.↗

NASA NTRS2011

Reduction of Aerosol Absorption in Beijing Since 2007 from MODIS and AERONET

An analysis of the time series of MODIS-based and AERONET aerosol records over Beijing reveals two distinct periods, before and after 2007. The MODIS data from both the Terra and Aqua satellites were processed with the new Multi-Angle Implementation of Atmospheric Correction (MAIAC) algorithm. A comparison of MAIAC and AERONET AOT shows that whereas MAIAC consistently underestimated peak AOT values by 10-20% in the prior period, the bias mostly disappears after mid- 2007. Independent analysis of the AERONET dataset reveals little or no change in the effective radii of the fine and coarse fractions and of the Angstrom exponent. At the same time, it shows an increasing trend in the single scattering albedo, by ~0.02 in 9 years. As MAIAC was using the same aerosol model for the entire 2000-2010 period, the decrease in AOT bias after 2007 can be explained only by a corresponding decrease of aerosol absorption caused by a reduction in local black carbon emissions. The observed changes correlate in time with the Chinese government's broad measures to improve air quality in Beijing during preparations for the Summer Olympics of 2008.

Lyapustin, A.↗

NASA NTRS2011

Improving the Automatic Inversion of Digital ISIS-2 Ionogram Reflection Traces into Topside Vertical Electron-Density Profiles

The topside-sounders on the four satellites of the International Satellites for Ionospheric Studies (ISIS) program were designed as analog systems. The resulting ionograms were displayed on 35-mm film for analysis by visual inspection. Each of these satellites, launched between 1962 and 1971, produced data for 10 to 20 years. A number of the original telemetry tapes from this large data set have been converted directly into digital records. Software, known as the TOPside Ionogram Scalar with True-height (TOPIST) algorithm has been produced that enables the automatic inversion of ISIS-2 ionogram reflection traces into topside vertical electron-density profiles Ne(h). More than million digital Alouette/ISIS topside ionograms have been produced and over 300,000 are from ISIS 2. Many of these ISIS-2 ionograms correspond to a passive mode of operation for the detection of natural radio emissions and thus do not contain ionospheric reflection traces. TOPIST, however, is not able to produce Ne(h) profiles from all of the ISIS-2 ionograms with reflection traces because some of them did not contain frequency information. This information was missing due to difficulties encountered during the analog-to-digital conversion process in the detection of the ionogram frame-sync pulse and/or the frequency markers. Of the many digital topside ionograms that TOPIST was able to process, over 200 were found where direct comparisons could be made with Ne(h) profiles that were produced by manual scaling in the early days of the ISIS program. While many of these comparisons indicated excellent agreement (<10% average difference over the entire profile) there were also many cases with large differences (more than a factor of two). Here we will report on two approaches to improve the automatic inversion process: (1) improve the quality of the digital ionogram database by remedying the missing frequency-information problem when possible, and (2) using the above-mentioned comparisons as teaching examples of how to improve the original TOPIST software.

Benson, R. F.↗

NASA NTRS2011

Multiangle Implementation of Atmospheric Correction (MAIAC): 2. Aerosol Algorithm

An aerosol component of a new multiangle implementation of atmospheric correction (MAIAC) algorithm is presented. MAIAC is a generic algorithm developed for the Moderate Resolution Imaging Spectroradiometer (MODIS), which performs aerosol retrievals and atmospheric correction over both dark vegetated surfaces and bright deserts based on a time series analysis and image-based processing. The MAIAC look-up tables explicitly include surface bidirectional reflectance. The aerosol algorithm derives the spectral regression coefficient (SRC) relating surface bidirectional reflectance in the blue (0.47 micron) and shortwave infrared (2.1 micron) bands; this quantity is prescribed in the MODIS operational Dark Target algorithm based on a parameterized formula. The MAIAC aerosol products include aerosol optical thickness and a fine-mode fraction at resolution of 1 km. This high resolution, required in many applications such as air quality, brings new information about aerosol sources and, potentially, their strength. AERONET validation shows that the MAIAC and MOD04 algorithms have similar accuracy over dark and vegetated surfaces and that MAIAC generally improves accuracy over brighter surfaces due to the SRC retrieval and explicit bidirectional reflectance factor characterization, as demonstrated for several U.S. West Coast AERONET sites. Due to its generic nature and developed angular correction, MAIAC performs aerosol retrievals over bright deserts, as demonstrated for the Solar Village Aerosol Robotic Network (AERONET) site in Saudi Arabia.

Lyapustin, A.↗

NASA NTRS2011

Reduction of Aerosol Absorption in Beijing Since 2007 from MODIS and AERONET

An analysis of the time series of MODIS-based and AERONET aerosol records over Beijing reveals two distinct periods, before and after 2007. The MODIS data from both the Terra and Aqua satellites were processed with the new Multi-Angle Implementation of Atmospheric Correction (MAIAC) algorithm. A comparison of MAIAC and AERONET AOT shows that whereas MAIAC consistently underestimated peak AOT values by 10-20% in the prior period, the bias mostly disappears after mid-2007. Independent analysis of the AERONET dataset reveals little or no change in the effective radii of the fine and coarse fractions and of the Angstrom exponent. At the same time, it shows an increasing trend in the single scattering albedo, by approx.0.02 in 9 years. As MAIAC was using the same aerosol model for the entire 2000-2010 period, the decrease in AOT bias after 2007 can be explained only by a corresponding decrease of aerosol absorption caused by a reduction in local black carbon emissions. The observed changes correlate in time with the Chinese government's broad measures to improve air quality in Beijing during preparations for the Summer Olympics of 2008.

Lyapustin, A.↗

NASA NTRS2011

Current Status of Multi-Angle Implementation of Atmospheric Correction (MAIAC) Algorithm

A new Multi-Angle Implementation of Atmospheric Correction (MAIAC) algorithm has been developed for MODIS. MAIAC uses a time series and an image based rather than pixel-based processing to perform simultaneous retrievals of aerosol properties and surface bidirectional reflectance. It is a generic algorithm which works over all land surface types with the exception of snow. MAIAC has an internal Cloud Mask, a dynamic land-water-snow classification and a surface change mask which allows it to flexibly choose processing path over different surfaces. A distinct feature of MAIAC is a high 1 km resolution of aerosol retrievals including optical thickness and fine mode fraction, which is required in different applications including the air quality analysis. An overview of the algorithm, results of AERONET validation, and examples of comparison with MODIS Collection 5 aerosol product, including Deep Blue algorithm, will be presented for different parts of the world including continental USA, Persian Gulf region and India.

Lyapustin, A.↗

NASA NTRS2011

Analysis of MAIAC Dust Aerosol Retrievals from MODIS Over North Africa

An initial comparison of aerosol optical thickness over North Africa for year 2007 was performed between the Deep Blue and Multi-Angle Implementation of Atmospheric Correction (MAIAC) algorithms complimented with MISR and OMI data. The new MAIAC algorithm has a better sensitivity to the small dust storms than the DB algorithm, but it also has biases in the brightest desert regions indicating the need for improvement. The quarterly averaged AOT values in the Bodele depression and western downwind transport region show a good agreement among MAIAC, MISR and OMI data, while the DB algorithm shows a somewhat different seasonality.

Lyapustin, A.↗

NASA NTRS2010

Analysis of Snow Bidirectional Reflectance from ARCTAS Spring-2008 Campaign

The spring 2008 Arctic Research of the Composition of the Troposphere from Aircraft and Satellites (ARCTAS) experiment was one of major intensive field campaigns of the International Polar Year aimed at detailed characterization of atmospheric physical and chemical processes in the Arctic region. A part of this campaign was a unique snow bidirectional reflectance experiment on the NASA P-3B aircraft conducted on 7 and 15 April by the Cloud Absorption Radiometer (CAR) jointly with airborne Ames Airborne Tracking Sunphotometer (AATS) and ground-based Aerosol Robotic Network (AERONET) sunphotometers. The CAR data were atmospherically corrected to derive snow bidirectional reflectance at high 1 degree angular resolution in view zenith and azimuthal angles along with surface albedo. The derived albedo was generally in good agreement with ground albedo measurements collected on 15 April. The CAR snow bidirectional reflectance factor (BRF) was used to study the accuracy of analytical Ross-Thick Li-Sparse (RTLS), Modified Rahman-Pinty-Verstraete (MRPV) and Asymptotic Analytical Radiative Transfer (AART) BRF models. Except for the glint region (azimuthal angles phi less than 40 degrees), the best fit MRPV and RTLS models fit snow BRF to within 0.05. The plane-parallel radiative transfer (PPRT) solution was also analyzed with the models of spheres, spheroids, randomly oriented fractal crystals, and with a synthetic phase function. The latter merged the model of spheroids for the forward scattering angles with the fractal model in the backscattering direction. The PPRT solution with synthetic phase function provided the best fit to measured BRF in the full range of angles. Regardless of the snow grain shape, the PPRT model significantly over-/underestimated snow BRF in the glint/backscattering regions, respectively, which agrees with other studies. To improve agreement with experiment, we introduced a model of macroscopic snow surface roughness by averaging the PPRT solution over the slope distribution function and by adding a simple model of shadows. With macroscopic roughness described by two parameters, the AART model achieved an accuracy of about plus or minus 0.05 with a possible bias of plus or minus 0.03 in the spectral range 0.4-2.2 micrometers. This high accuracy holds at view zenith angles below 55-60 degrees covering the practically important range for remote sensing applications, and includes both glint and backscattering directions.

Lyapustin, A.↗

NASA NTRS2009

Space Technology 5 Multipoint Observations of Temporal and Spatial Variability of Field-Aligned Currents

Space Technology 5 (ST5) is a constellation mission consisting of three microsatellites. It provides the first multipoint magnetic field measurements in low Earth orbit, which enables us to separate spatial and temporal variations. In this paper, we present a study of the temporal variability of field-aligned currents using the ST5 data. We examine the field-aligned current observations during and after a geomagnetic storm and compare the magnetic field profiles at the three spacecraft. The multipoint data demonstrate that mesoscale current structures, commonly embedded within large-scale current sheets, are very dynamic with highly variable current density and/or polarity in approx.10 min time scales. On the other hand, the data also show that the time scales for the currents to be relatively stable are approx.1 min for mesoscale currents and approx.10 min for large-scale currents. These temporal features are very likely associated with dynamic variations of their charge carriers (mainly electrons) as they respond to the variations of the parallel electric field in auroral acceleration region. The characteristic time scales for the temporal variability of mesoscale field-aligned currents are found to be consistent with those of auroral parallel electric field.

Le, G.↗