Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “weather normalization”

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 217 records · Page 12

Multi-Sensor Cloud and Aerosol Retrieval Simulator and Remote Sensing from Model Parameters : Aerosols - Part 2

The Multi-sensor Cloud Retrieval Simulator (MCRS) produces a simulated radiance product from any high-resolution general circulation model with interactive aerosol as if a specific sensor such as the Moderate Resolution Imaging Spectroradiometer (MODIS) were viewing a combination of the atmospheric column and land ocean surface at a specific location. Previously the MCRS code only included contributions from atmosphere and clouds in its radiance calculations and did not incorporate properties of aerosols. In this paper we added a new aerosol properties module to the MCRS code that allows users to insert a mixture of up to 15 different aerosol species in any of 36 vertical layers. This new MCRS code is now known as MCARS (Multi-sensor Cloud and Aerosol Retrieval Simulator). Inclusion of an aerosol module into MCARS not only allows for extensive, tightly controlled testing of various aspects of satellite operational cloud and aerosol properties retrieval algorithms, but also provides a platform for comparing cloud and aerosol models against satellite measurements. This kind of two-way platform can improve the efficacy of model parameterizations of measured satellite radiances, allowing the assessment of model skill consistently with the retrieval algorithm. The MCARS code provides dynamic controls for appearance of cloud and aerosol layers. Thereby detailed quantitative studies of the impacts of various atmospheric components can be controlled. In this paper we illustrate the operation of MCARS by deriving simulated radiances from various data field output by the Goddard Earth Observing System version 5 (GEOS-5) model. The model aerosol fields are prepared for translation to simulated radiance using the same model sub grid variability parameterizations as are used for cloud and atmospheric properties profiles, namely the ICA technique. After MCARS computes modeled sensor radiances equivalent to their observed counterparts, these radiances are presented as input to operational remote-sensing algorithms. Specifically, the MCARS-computed radiances are input into the processing chain used to produce the MODIS Data Collection 6 aerosol product (MOYD04). TheMOYD04 product is of course normally produced from MOYD021KM MODIS Level-1B radiance product directly acquired by the MODIS instrument. MCARS matches the format and metadata of a MOYD021KM product. The resulting MCARS output can be directly provided to MODAPS (MODIS Adaptive Processing System) as input to various operational atmospheric retrieval algorithms. Thus the operational algorithms can be tested directly without needing to make any software changes to accommodate an alternative input source. We show direct application of this synthetic product in analysis of the performance of the MOD04 operational algorithm. We use biomass-burning case studies over Amazonia employed in a recent Working Group on Numerical Experimentation (WGNE)-sponsored study of aerosol impacts on numerical weather prediction (Freitas et al., 2015). We demonstrate that a known low bias in retrieved MODIS aerosol optical depth appears to be due to a disconnect between actual column relative humidity and the value assumed by the MODIS aerosol product.

aerosol retrieval↗

Variables Influencing RNAV STAR Adherence

In this study we investigated how variables in the aviation domain impact adherence levels of aircraft flying area navigation arrivals with optimized profile descents (RNAV OPDs) (RNAV STARs: aRea NAVigation Standard Terminal Arrival Routes). Variable categories were: weather, aircraft, procedure, and traffic. Non-adherence events analyzed were: miss above, miss below, skip before merge, and skip after merge. Miss below and miss above describe when a flight does not comply vertically with a procedure. Skips refer to a flight leaving a procedure, then returning. Findings of this work reveal that vertical events are most impacted by altitude restriction size, steepness of flight paths, and merging routes. Lateral events were impacted by merging flight conflicts, number of speed restrictions, and the flow rate of the arrival traffic. This study helps increase understanding of how the system is functioning and identifies where procedures are not flexible enough to handle the variability in normal operations.major airports, procedure design, and recommendations for future work.

RNAV OPD↗

Midwest Food Security & Agriculture II: Leveraging NASA Earth Observations to Analyze and Display Crop Phenology Data and Weather Conditions to Support Expansion of Small Grain Crops in the Midwest

Agriculture in the Midwest is dominated by monoculture systems that strip the soil of nutrients, decrease yields, and worsen water quality. Crop diversification and cultivating small grains is economically and ecologically advantageous, but limited in practice due to a lack of data about small-grain crop performance as well as monoculture-favoring insurance coverage, economic incentives, and sociocultural traditions. Practical Farmers of Iowa (PFI) strives to create more sustainable agriculture and supports Midwestern farmers in adding small grains into their crop rotations. The DEVELOP team partnered with the USDA Agricultural Research Service and the USA National Phenology Network to assist PFI in applying NASA Earth observation data to track small-grain crop performance. This project developed a Google Earth Engine user interface that enables PFI farmers to compare satellite-assessed vegetation production and climatic conditions of their farm with other PFI farmers. The team assessed the crop performance through Normalized Difference Vegetation Index (NDVI) values derived from Landsat 7 Enhanced Thematic Mapper Plus (ETM+) and Landsat 8 Operational Land Imager (OLI). Additionally, the team created a regression model to analyze the predictive properties of vegetation indices and crop yield, ultimately establishing a relationship with mean NDVI from the early growing season. PFI farmers can utilize these products to make informed decisions about their crop production as well as advocate for more expansive insurance policies to protect small-grain crops.

Sophie Barrowman↗

Measurements of the earth radiation budget from satellites during the first GARP global experiment

Radiation budget data (which will aid in climate model development) and solar constant measurements (both to be used for the study of long term climate change and interannual seasonal weather variability) are presented, obtained during Nimbus-6 and Nimbus-7 satellite flights, using wide-field-of-view, scanner, and black cavity detectors. Data on the solar constant, described as a function of the date of measurement, are given. The unweighed mean amounts to 1377 + or - 20 per sq Wm, with a standard deviation of 8 per sq Wm. The new solar data are combined with earlier measurements, and it is suggested that the total absolute energy output of the sun is a minimum at 'solar maximum' and vice versa. Attention is given to the measurements of the net radiation budget, the planetary albedo, and the infrared radiant exitance. The annual and semiannual cycles of normal variability explain most of the variance of energy exchange between the earth and space. Examination of separate ocean and atmospheric energy budgets implies a net continent-ocean region energy exchange.

Vonder Haar, T. H.↗

Estimates of Long Term Surface Soil Moisture in the Midwestern U.S. Derived from Satellite Microwave Observations

Soil moisture is a key component of the water and energy balances of the Earth's surface, and has been identified as a parameter of significant potential for improving the accuracy of large-scale land surface-atmosphere interaction models. However, soil moisture is often somewhat difficult to measure accurately in both space and time, especially at large spatial scales. Soil moisture is highly variable, and while point measurements are typically quite accurate, subsequent areal averaging of these measurements often leads to large errors. Since remotely sensed land surface observations are already a spatially averaged or areally integrated value, they are a logical input parameter to regional or larger scale land process models. A database of long-term soil moisture was compared to satellite microwave observations over test sites in the Midwestern United States. Ground measurements of average volumetric surface soil moisture in the top ten cm were made bimonthly at 19 locations throughout the state of Illinois. Nighttime microwave brightness temperatures were observed at a frequency of 6.6 GHz, by the Scanning Multichannel Microwave Radiometer (SMMR), onboard the Nimbus 7 satellite. The life of the SMMR instrument spanned from Nov. 1978 to Aug. 1987. At 6.6 GHz, the instrument provided a spatial resolution of approximately 150 km, and a temporal frequency over the test area of about 3 nighttime orbits per week. Vegetation radiative transfer characteristics, such as the canopy transmissivity, were estimated from vegetation indices such as the Normalized Difference Vegetation Index (NDVI) and the 37 GHz Microwave Polarization Difference Index (MPDI). Because the time of satellite coverage does not always coincide with the ground measurements of soil moisture, the existing ground data were used to calibrate a water balance for the top IO cm surface layer in order to interpolate daily surface moisture values. Such a climate-based approach is often more appropriate for estimating large-area average soil moisture because meteorological data are generally more spatially representative than isolated point measurements of soil moisture, Passive microwave remote sensing presents the greatest potential for providing regular spatially representative estimates of surface soil moisture at global scales. Real time estimates should improve weather and climate modelling efforts, while the development of historical data sets will provide necessary information for simulation and validation of long-term climate and global change studies.

Owe, M.↗

The necessity for a new parameterization of an empirical model for wind/ocean scatterometry

Difficulties related to the paucity of weather observation data regarding oceans were potentially alleviated for three summer months in 1978 when NASA's Seasat telemetered data from three wind-measuring instruments. The present study is concerned with one of these instruments, the Seasat A Scatterometer System (SASS). Attention is given to an internal consistency check of the SASS 1 model, comparisons of SASS 1 and other model-predicted winds with in situ winds, and a brief summary of the principal findings. It is found that a new wind retrieval system is required if SASS wind data are to be globally applicable and, at the same time, are to meet the required performance specifications. The sum-of-squares (SOS) technique for inverting SASS NRCS (normalized radar cross section) measurements results in the discarding of valuable data in low-speed areas and for higher incidence angles.

Woiceshyn, P. M.↗

The value of NDBC observations during March 1993's 'Storm of the Century'

National Data Buoy Center (NDBC) observations from the western Gulf of Mexico Provided several early indications of more rapid cyclogenesis on 12 March 1993 than was forecast by numerical weather prediction (NWP) models. Observations demostrated a closed circulation with gale-force winds shortly after the storm entered the Gulf of Mexico. Pressure measurements at two buoys off the Texas coast were 4-6 hPa lower than the 12-h NWP forecasts, a significant forecast error. Observations from NDBC's moored buoys and Coastal-Marine Automated Network stations revealed that the developing storm was significantly deeper than was indicated on the National Meteorological Center's automated surface analyses. Ocean wave observations reveal some of the steepest waves NDBC has ever measured, indicating phenomenonal wave growth and a high potential for damage to vessels and structures. A warm eddy caused sea surface temperatures (SST) to be several degrees above normal under the track of the storm, creating a strong SST gradient to the north. This provided ample energy and strengthened the baroclinity. NDBC observations showed the eddy to be somewhat larger and warmer than indicated by the most recent National Hurricane Center analysis. This event demonstrates the tremendous value of NDBC marine observations for accurately detecting the occurrence and strength of coastal cyclogenesis events.

Gilhousen, David B.↗

Radar-anomalous, high-altitude features on Venus

Over nearly all of the surface of Venus the reflectivity and emissivity at centimeter wavelengths are about 0.15 and 0.85 respectively. These values are consistent with moderately dense soils and rock populations, but the mean reflectivity is about a factor of 2 greater than that for the Moon and other terrestrial planets. Pettingill and Ford, using Pioneer Venus reflectivities and emissivities, found a number of anomalous features on Venus that showed much higher reflectivities and much lower emissivities with both values approaching 0.5. These include Maxwell Montes, a number of high regions in Aphrodite Terra and Beta Regio, and several isolated mountain peaks. Most of the features are at altitudes above the mean radius by 2 to 3 km or more. However, such features have been found in the Magellan data at low altitudes and the anomalies do not exist on all high structures, Maat Mons being the most outstanding example. A number of papers have been written that attempt to explain the phenomena in terms of the geochemistry balance of weathering effects on likely surface minerals. The geochemists have shown that the fundamentally basaltic surface would be stable at the temperatures and pressures of the mean radius in the form of magnetite, but would evolve to pyrite and/or pyrrhotite in the presence of sulfur-bearing compounds such as SO2. Pyrite will be stable at altitudes above 4 or 5 km on Venus. Although the geochemical arguments are rather compelling, it is vitally important to rationally look at other explanations for radar and radio emission measurements such as that presented by Tryka and Muhleman. The radar reflectivity values are retrieved from the raw Magellan backscatter measurements by fitting the Hagfors' radar scattering model in which a surface roughness parameters and a normal incidence electrical reflectivity are estimated. The assumptions of the theory behind the model must be considered carefully before the results can be believed. These include that the surface roughness exists only at horizontal scales large compared to the wavelength, the vertical deviations are gaussianly distributed, there is no shadowing, and that the reflection occurs at the interface of two homogeneous dielectric half-spaces. Probably all these conditions are violated at the anomalous features under discussion. The most important of these is the homogeneity of the near surface of Venus, particularly in highlands. Under the assumptions of the theory, all of the radio energy is reflected by the impedance jump at the very boundary. However, in heterogeneous soil some fraction of the illuminating energy is propagated into the soil and then scattered back out by impedance discontinuities such as rock, voids, and cracks. In light soils, the latter effect can overwhelm the scattering effects of the true surface and greatly enhance the backscatter power, suggesting a much higher value of an effective dielectric constant that would be estimated from Hagfors' model.

Muhleman, Duane O.↗

Analysis of Aircraft Clusters to Measure Sector-Independent Airspace Congestion

The Distributed Air/Ground Traffic Management (DAG-TM) concept of operations* permits appropriately equipped aircraft to conduct Free Maneuvering operations. These independent aircraft have the freedom to optimize their trajectories in real time according to user preferences; however, they also take on the responsibility to separate themselves from other aircraft while conforming to any local Traffic Flow Management (TFM) constraints imposed by the air traffic service provider (ATSP). Examples of local-TFM constraints include temporal constraints such as a required time of arrival (RTA), as well as spatial constraints such as regions of convective weather, special use airspace, and congested airspace. Under current operations, congested airspace typically refers to a sector(s) that cannot accept additional aircraft due to controller workload limitations; hence Dynamic Density (a metric that is indicative of controller workload) can be used to quantify airspace congestion. However, for Free Maneuvering operations under DAG-TM, an additional metric is needed to quantify the airspace congestion problem from the perspective of independent aircraft. Such a metric would enable the ATSP to prevent independent aircraft from entering any local areas of congestion in which the flight deck based systems and procedures may not be able to ensure separation. This new metric, called Gaggle Density, offers the ATSP a mode of control to regulate normal operations and to ensure safety and stability during rare-normal or off-normal situations (e.g., system failures). It may be difficult to certify Free Maneuvering systems for unrestricted operations, but it may be easier to certify systems and procedures for specified levels of Gaggle Density that could be monitored by the ATSP, and maintained through relatively minor flow-rate (RTA type) restrictions. Since flight deck based separation assurance is airspace independent, the challenge is to measure congestion independent of sector boundaries. Figure 1 , reproduced from Ref. 1, depicts an example traffic situation. When the situation is analyzed by sector boundaries (left side of figure), a Dynamic Density metric would identify excessive congestion in the central sector. When the same traffic situation is analyzed independent of sector boundaries (right side of figure), a Gaggle Density metric would identify congestion in two dynamically defined areas covering portions of several sectors. The first step towards measuring airspace-independent congestion is to identify aircraft clusters, i.e., groups of closely spaced aircraft. The objective of this work is to develop techniques to detect and classify clusters of aircraft.

Bilimoria, Karl D.↗

The Transition of Atmospheric Infrared Sounder Total Ozone Products to Operations

The National Aeronautics and Space Administration Short-term Prediction Research and Transition Center (NASA SPoRT) has transitioned a total column ozone product from the Atmospheric Infrared Sounder (AIRS) retrievals to the Weather Prediction Center and Ocean Prediction Center. The total column ozone product is used to diagnose regions of warm, dry, ozone-rich, stratospheric air capable of descending to the surface to create high-impact non-convective winds. Over the past year, forecasters have analyzed the Red, Green, Blue (RGB) Air Mass imagery in conjunction with the AIRS total column ozone to aid high wind forecasts. One of the limitations of the total ozone product is that it is difficult for forecasters to determine whether elevated ozone concentrations are related to stratospheric air or climatologically high values of ozone in certain regions. During the summer of 2013, SPoRT created an AIRS ozone anomaly product which calculates the percent of normal ozone based on a global stratospheric ozone mean climatology. With the knowledge that ozone values 125 percent of normal and greater typically represent stratospheric air; the anomaly product can be used with the total column ozone product to confirm regions of stratospheric air. This paper describes the generation of these products along with forecaster feedback concerning the use of the AIRS ozone products in conjunction with the RGB Air Mass product to access the utility and transition of the products.

Berndt, E. B.↗

The Transition of Atmospheric Infrared Sounder Total Ozone Products to Operations

The National Aeronautics and Space Administration Short-term Prediction Research and Transition Center (NASA SPoRT) has transitioned a total column ozone product from the Atmospheric Infrared Sounder (AIRS) retrievals to the Weather Prediction Center and Ocean Prediction Center. The total column ozone product is used to diagnose regions of warm, dry, ozone-rich, stratospheric air capable of descending to the surface to create high-impact non-convective winds. Over the past year, forecasters have analyzed the Red, Green, Blue (RGB) Air Mass imagery in conjunction with the AIRS total column ozone to aid high wind forecasts. One of the limitations of the total ozone product is that it is difficult for forecasters to determine whether elevated ozone concentrations are related to stratospheric air or climatologically high values of ozone in certain regions. During the summer of 2013, SPoRT created an AIRS ozone anomaly product which calculates the percent of normal ozone based on a global stratospheric ozone mean climatology. With the knowledge that ozone values 125 percent of normal and greater typically represent stratospheric air; the anomaly product can be used with the total column ozone product to confirm regions of stratospheric air. This paper describes the generation of these products along with forecaster feedback concerning the use of the AIRS ozone products in conjunction with the RGB Air Mass product to access the utility and transition of the products.

Berndt, Emily↗

Monitoring Airspace Complexity and Determining Contributing Factors

The national airspace has evolved over many years to accommodate increased traffic demand [1] while simultaneously maintaining one of the safest forms of transportation [2], [3]. One of the reasons for this success is the ability of the system and the operators to adapt and accommodate to situations that routinely disrupt optimal operations. These situations may include: adverse weather, delays, early arrivals, equipment outages, and other factors that are outside the operators ability to control. These factors can lead to states where automation is unable to properly handle these issues and therefore air traffic controllers and pilots have to intervene, ultimately increasing communication between operators resulting in higher workload. As controller workload increases to handle sub-optimal operating conditions this can be viewed as an increase in complexity. The reasoning for this is because humans are now required to make tactical decisions in response to external factors, resulting in a departure from the strategic plan where operations would be more efficiently managed. Human operators control airspace complexity under rigid regulations that are constantly changing. The airspace is divided into sectors and the number of aircraft assigned to each controller is limited for safe handling. There has been past work that devised airspace complexity metrics in commercial aviation and related these metrics to controller workload (e.g., [4],[5]). The upper bounds on the system load are pre-determined. Such bounds on complexity make for a safe system, but the system cannot scale and adapt to autonomous, dense, and heterogeneous traffic, including the many types of Unmanned Aerial Vehicles (UAVs) envisioned to be added to the operations. We hypothesize that, as traffic density and heterogeneity grow, and other key metrics change, there will be phase transitions at which the way traffic should be managed changes significantly [6]. We offer a method for in-time detection of contributing factors that lead to phase transitions, characterized by increased complexity. To the best of our knowledge, there is no tool similar to our proposed effort that identifies such contributing factors or precursor patterns. To define the scope we are proposing to measure complexity from the viewpoint of the Terminal Radar Approach Control Facilities (TRACON) controller’s perspective. In particular we are analyzing arrivals into KSFO. With safety as the top concern for airspace operators, it is important to recognize that as density and heterogeneity grow, the focus of the system will change. Times of the day when the airspace has low density and heterogeneity, the flights will follow more efficient paths where the aircraft move on established routes that are more or less directly to the destination. However, when density and heterogeneity increases, the system will begin changing focus to avoiding conflicts and collisions and route the flights in a more flexible way. Higher flexibility requires more communication and coordination between controllers and pilots which the current automation is unable to handle. This paper proposes a novel approach that monitors airspace complexity at multiple scales, uses a Machine Learning-based tool that predicts when operations will transition to a regime of greater complexity, and identifies actions that can reduce the complexity while still maintaining efficient and safe operations. We demonstrate our proposed approach using data from multiple complementary sources. This includes, but is not limited to: historical aircraft surveillance data from NASA’s Sherlock Data Warehouse [7], METAR weather data, and airport configuration data from Aviation System Performance Metrics (ASPM). The surveillance data flight paths are sampled at a variable sample rate — increasing as the aircraft approaches the airport. This is due to how Sherlock manages flight track stitching between different radar facilities which have different sampling rates. The weather and performance data are logged at defined intervals throughout the day at a courser refresh rate. In addition to the logged data and metrics, we leverage pre-defined Standard Terminal Arrival Routes (STARs) procedures to characterize the path of each flight. Each flight files for one of these routes in the flight plan well before entering the terminal airspace, and approximately follows the route until it leaves the STAR, typically on the final fix of a runway transition. However, most flights do not always fly the full STAR procedure to completion [8], but the majority do adhere to the fixes within the common route of the procedure. Our approach leverages fixes in the common route of each of the STARs to build a reference path to the airport. This allows us to characterize the flight paths in what we are defining as the “maneuvering area” (the airspace between the STAR and before the flight is lined up on the runway’s final approach) to determine how off nominal the flights are to calculate its complexity score. Determining airspace complexity is a concept that does not have a concrete answer. In designing this metric, we consider what increases the workload for the air traffic controllers. Consequently more specialized vectoring maneuvers results in higher workload. Accordingly, we start with a theory: each flight has a direct path it takes from the STAR’s common route to the final approach’s outer marker fix for the flight’s landing runway. It is important to note that the direct path is only used as a reference. If the majority of the flights have a large consistent offset as compared to other routes it does not necessarily mean that those flights have higher complexity. We are merely building a distribution based on this direct path for that particular STAR and runway pair to determine the normal mode of operations for that route. Flights that are in the upper tail of these distributions will result in higher complexity scores and flights that fly in the median will represent the normal mode of operations and therefore will have lower complexity scores. Since flights following each STAR route take different paths to the airport, we have a different distribution for each STAR route and therefore can model these distributions to compute a complexity score from their respective normalized distributions. To evaluate the effectiveness of our proposed airspace complexity metric we will compare against an established approach based on trajectory clustering [9]. This unsupervised learning technique consists of the following steps: (1) identify the general maneuvering areas (waypoints) by performing $\kappa$-means or DBSCAN clustering on locations where aircraft frequently turn based on the surveillance radar track data, (2) map flight trajectories onto sequences of waypoints, and (3) cluster the sequences based on their common subsequences. From a high-level perspective, this baseline model learns nominal operations in the airspace through the sequence of waypoints that are representative of where aircraft change direction and defines deviations from the nominal operations as “complex.” Therefore, more deviations from the nominal operations correspond to higher complexity values. For our validation, we re-implemented this technique and tune model hyper-parameters to correctly detect waypoints for the arrival traffic into the San Francisco bay area. We will compute the complexity measure over a one-year period using our proposed technique as well as the baseline. Our validation will be based on each technique’s ability to detect a set of undesirable outcomes (e.g., go-arounds, holding patterns, average time in the airspace, etc.). Since our current complexity metric is derived from the offset from the direct reference path, it’s important to understand what causes these offsets. In many of the flights with high offset distance, flights performing holding patterns and S turns can be observed. These maneuvering tactics are utilized to add distance between the aircraft and the destination runway to prevent multiple flights from having conflicting arrival times. In order to predict a rise in complexity (or the precursor to complexity), it’s necessary to be able to identify these potential conflicts (which in turn, result in higher offsets). To do this, we define a “representative flight” for each STAR route and runway pair. This flight is approximately the path the flight would take if there was a clear path with no other flights in the airspace — including the time remaining to the airport. We first identify the flights for a given STAR runway pair using the offset to the reference path distributions that fall between the 44-55 percentiles. This yields the flights that conform to the most normal mode of operation. Each of these flights is partitioned based on the percent complete from the entry point into the maneuvering areas from 0\% – 100\% complete. Then for each percent “bin”, we take the median value of the flight’s latitude/longitude coordinates, airspeed, and (non causal) time remaining to the airport to construct a lookup table for each percent complete bin on a given route. As a flight enters the maneuvering area, we can find the estimated arrival time of a flight to the airport by finding the closest point to the representative path’s percent complete bin (relative to the flight’s current position at any snapshot in the airspace) and therefore retrieve the corresponding remaining time left on the “representative path”. We assume that the flight will follow the representative path to completion when deriving these estimates. We can then compare these estimated arrival times against other flights for the same snapshot in time to identify potential conflicts. If more flights are estimated to arrive within a tolerance window than there are runways available, then we have a potential conflict. We can use this derived measure along with other factors expected to add disruption to the operation such as weather and runway configuration changes as an input to machine learning tools to detect precursors that increases in our complexity measure. This novel method will assist in uncovering insights into the contributing factors that lead to increased complexity that may allow for in-time responses to avoid reaching a high complexity state in the airspace.

complexity↗

Initial Results from the Variable Intensity Sonic Boom Propagation Database

An extensive sonic boom propagation database with low- to normal-intensity booms (overpressures of 0.08 lbf/sq ft to 2.20 lbf/sq ft) was collected for propagation code validation, and initial results and flight research techniques are presented. Several arrays of microphones were used, including a 10 m tall tower to measure shock wave directionality and the effect of height above ground on acoustic level. A sailplane was employed to measure sonic booms above and within the atmospheric turbulent boundary layer, and the sailplane was positioned to intercept the shock waves between the supersonic airplane and the ground sensors. Sailplane and ground-level sonic boom recordings were used to generate atmospheric turbulence filter functions showing excellent agreement with ground measurements. The sonic boom prediction software PCBoom4 was employed as a preflight planning tool using preflight weather data. The measured data of shock wave directionality, arrival time, and overpressure gave excellent agreement with the PCBoom4-calculated results using the measured aircraft and atmospheric data as inputs. C-weighted acoustic levels generally decreased with increasing height above the ground. A-weighted and perceived levels usually were at a minimum for a height where the elevated microphone pressure rise time history was the straightest, which is a result of incident and ground-reflected shock waves interacting.

Haering, Edward A., Jr.↗

Imager for Mars Pathfinder (IMF)

The IMP camera is a near-surface sensing experiment with many capabilities beyond those normally associated with an imager. It is fully pointable in both elevation and azimuth with a protected, stowed position looking straight down. Stereo separation is provided with two optical paths; each has a 12-position filter wheel. The primary function of the camera, strongly tied to mission success, is to take a color panorama of the surrounding terrain. IMP requires approximately 120 images to give a complete downward hemisphere from the deployed position. IMP provides the geologist, and everyone else, a view of the local morphology with millimeter-tometer-scale resolution over a broad area. In addition to the general morphology of the scale, IMP has a large compliment of specially chosen filters to aid in both the identification of the mineral types and their degree of weathering.

Smith, Peter H.↗

A Data-driven, Physics-based Transport Model of Solar Energetic Particles Accelerated by Coronal Mass Ejection Shocks Propagating through the Solar Coronal and Heliospheric Magnetic Fields

In an effort to develop computational tools for predicting radiation hazards from solar energetic particles (SEPs), we have created a data-driven physics-based particle transport model to calculate the injection, acceleration, and propagation of SEPs from coronal mass ejection (CME) shocks traversing through the solar corona and interplanetary magnetic fields. The model runs on an input of corona and heliospheric plasma and magnetic field configuration from a magnetohydrodynamic model driven by solar photospheric magnetic field measurements superposed with observed CME shocks determined from coronagraph images. SEP source particles are injected at the shock using the result of diffusive shock acceleration formulation from a characteristic obliquity-dependent injection from a heated solar wind thermal tail population. With several advanced computation techniques involving stochastic simulation and integration, the model obtains the particle intensity at any location in interplanetary space through the rigorous solution to the time-dependent 5D focus transport equation in the phase space that includes perpendicular diffusion. We apply the model to the 2011 November 3 CME event. The calculation results reproduce multispacecraft SEP observations at Earth and STEREO-B reasonably well without normalization of particle flux. The observations at STEREO-A can be reproduced by rescaling particle energy or modified energy dependence of particle diffusion coefficients. This circumsolar SEP event seen by spacecraft at Earth, STEREO-A, and STEREO-B at widely separated longitudes can be explained by diffusive shock acceleration by a single CME shock with a moderate speed.

Solar energetic particles↗

Prediction of Coronal Mass Ejections From Vector Magnetograms: Quantitative Measures as Predictors

We derived two quantitative measures of an active region's global nonpotentiality from the region's vector magnetogram, 1) the net current (I(sub N)), and 2) the length of strong-shear, strong-field main neutral line (Lss), and used these two measures in a pilot study of the CME productivity of 4 active regions. We compared the global nonpotentiality measures to the active regions' CME productivity determined from GOES and Yohkoh/SXT observations. We found that two of the active regions were highly globally nonpotential and were CME productive, while the other two active regions had little global nonpotentiality and produced no CMEs. At the Fall 2000 AGU, we reported on an expanded study (12 active regions and 17 magnetograms) in which we evaluated four quantitative global measures of an active region's magnetic field and compared these measures with the CME productivity. The four global measures (all derived from MSFC vector magnetograms) included our two previous measures (I(sub N) and L(sub ss)) as well as two new ones, the total magnetic flux (PHI) (a measure of an active region's size), and the normalized twist (alpha (bar)= muIN/PHI). We found that the three quantitative measures of global nonpotentiality (I(sub N), L(sub ss), alpha (bar)) were all well correlated (greater than 99% confidence level) with an active region's CME productivity within plus or minus 2 days of the day of the magnetogram. We will now report on our findings of how good our quantitative measures are as predictors of active-region CME productivity, using only CMEs that occurred after the magnetogram. We report the preliminary skill test of these quantitative measures as predictors. We compare the CME prediction success of our quantitative measures to the CME prediction success based on an active region's past CME productivity. We examine the cases of the handful of false positive and false negatives to look for improvements to our predictors. This work is funded by NSF through the Space Weather Program and by NASA through the Solar Physics Supporting Research and Technology Program.

Falconer, D. A.↗

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↗