Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Heterogeneity”

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 541 records · Page 30

Photochemical modeling of the Antarctic stratosphere: Observational constraints from the airborne Antarctic ozone experiment and implications for ozone behavior

The rapid decrease in O3 column densities observed during Antarctic spring has been attributed to several chemical mechanisms involving nitrogen, bromine, or chlorine species, to dynamical mechanisms, or to a combination of the above. Chlorine-related theories, in particular, predict greatly elevated concentrations of ClO and OClO and suppressed abundances of NO2 below 22 km. The heterogeneous reactions and phase transitions proposed by these theories could also impact the concentrations of HCl, ClNO3 and HNO3 in this region. Observations of the above species have been carried out from the ground by the National Ozone Expedition (NOZE-I, 1986, and NOZE-II, 1987), and from aircrafts by the Airborne Antarctic Ozone Experiment (AAOE) during the austral spring of 1987. Observations of aerosol concentrations, size distribution and backscattering ratio from AAOE, and of aerosol extinction coefficients from the SAM-II satellite can also be used to deduce the altitude and temporal behavior of surfaces which catalyze heterogeneous mechanisms. All these observations provide important constraints on the photochemical processes suggested for the spring Antarctic stratosphere. Results are presented for the concentrations and time development of key trace gases in the Antarctic stratosphere, utilizing the AER photochemical model. This model includes complete gas-phase photochemistry, as well as heterogeneous reactions. Heterogeneous chemistry is parameterized in terms of surface concentrations of aerosols, collision frequencies between gas molecules and aerosol surfaces, concentrations of HCl/H2O in the frozen particles, and probability of reaction per collision (gamma). Values of gamma are taken from the latest laboratory measurements. The heterogeneous chemistry and phase transitions are assumed to occur between 12 and 22 km. The behavior of trace species at higher altitudes is calculated by the AER 2-D model without heterogeneous chemistry. Calculations are performed for solar illumination conditions typical of 60, 70, and 80 S, from July 15 to October 31.

Rodriguez, Jose M.↗

The chemical and radiative effects of the Mount Pinatubo eruption

To clarify the mechanisms leading to effects on stratospheric ozone, time-dependent stratospheric aerosol and gas experiment II (SAGE II) and cryogenic limb array elaton spectrometer (CLAES) aerosol optical extinction data and SAGE II surface area density are used as parameters in a two-dimensional (2-D) zonally averaged chemical radiative transport model. The model was integrated with time from before the eruption through December 1993. The modeled impact on global ozone results from increased rates of heterogeneous reactions on sulfate aerosols and from the increased radiative heating and scattering caused by these aerosols. When the aerosol heating is allowed to modify the temperature distribution, the maximum change calculated in equatorial column ozone is -1.6%. The calculated equatorial temperature change and peak local ozone change in October 1991 are +6K and -4%, respectively. When aerosol heating perturbs the circulation in the model, the maximum change in equatorial column ozone is -6%. Increased heterogeneous processing on sulfate aerosols is calculated to have changed equatorial column ozone in late 1991 by -1.5%. Global column ozone in the model in 1992 and 1993 changed by -2.8% and -2.4%, respectively. The relationship of ozone-controlling processes in the lower stratosphere is altered as well; HO(x) becomes the most important catalytic cycle, followed by ClO(x) and NO(x). This is driven by significant changes in trace gas concentrations. In October 1991, lower stratospheric, equatorial NO(x) decreased by 40%, ClO(x) increased by 60%, and HO(x) increased by 25%. When the effect of heterogeneous chemical processing on sulfate aerosols is combined with aerosol heating, modifying either circulation or temperature, dramatically different ozone fingerprints with time and latitude are predicted. Model-derived changes in the equatorial region in column ozone best represented the observed data when perturbed circulation was combined with heterogeneous chemical effects. However, at high latitudes, the increased ozone production from the strengthening of the mean circulation tends to cancel the heterogeneous reduction of ozone. This is not in good agreement with observed data, especially in 1992 and 1993. When the circulation is held fixed and the temperature allowed to change, and heterogeneous chemical effects are included, the equatorial ozone decrease predicted was too small for 1991. However, the mid- to high-latitude decrease in 1992 and 1993 is in better agreement with observed data.

Kinneson, Douglas E.↗

Development of a Novel Method for Determination of Residual Stresses in a Friction Stir Weld

Material constitutive properties, which describe the mechanical behavior of a material under loading, are vital to the design and implementation of engineering materials. For homogeneous materials, the standard process for determining these properties is the tensile test, which is used to measure the material stress-strain response. However, a majority of the applications for engineering materials involve the use of heterogeneous materials and structures (i.e. alloys, welded components) that exhibit heterogeneity on a global or local level. Regardless of the scale of heterogeneity, the overall response of the material or structure is dependent on the response of each of the constituents. Therefore, in order to produce materials and structures that perform in the best possible manner, the properties of the constituents that make up the heterogeneous material must be thoroughly examined. When materials exhibit heterogeneity on a local level, such as in alloys or particle/matrix composites, they are often treated as statistically homogenous and the resulting 'effective' properties may be determined through homogenization techniques. In the case of globally heterogeneous materials, such as weldments, the standard tensile test provides the global response but no information on what is Occurring locally within the different constituents. This information is necessary to improve the material processing as well as the end product.

Reynolds, Anthony P.↗

A Proposed Extension to the Soil Moisture and Ocean Salinity Level 2 Algorithm for Mixed Forest and Moderate Vegetation Pixels

The Soil Moisture and Ocean Salinity (SMOS)mission, launched in November 2009, provides global maps of soil moisture and ocean salinity by measuring the L-band (1.4 GHz) emission of the Earth's surface with a spatial resolution of 40-50 km.Uncertainty in the retrieval of soilmoisture over large heterogeneous areas such as SMOS pixels is expected, due to the non-linearity of the relationship between soil moisture and the microwave emission. The current baseline soilmoisture retrieval algorithm adopted by SMOS and implemented in the SMOS Level 2 (SMOS L2) processor partially accounts for the sub-pixel heterogeneity of the land surface, by modelling the individual contributions of different pixel fractions to the overall pixel emission. This retrieval approach is tested in this study using airborne L-band data over an area the size of a SMOS pixel characterised by a mix Eucalypt forest and moderate vegetation types (grassland and crops),with the objective of assessing its ability to correct for the soil moisture retrieval error induced by the land surface heterogeneity. A preliminary analysis using a traditional uniform pixel retrieval approach shows that the sub-pixel heterogeneity of land cover type causes significant errors in soil moisture retrieval (7.7%v/v RMSE, 2%v/v bias) in pixels characterised by a significant amount of forest (40-60%). Although the retrieval approach adopted by SMOS partially reduces this error, it is affected by errors beyond the SMOS target accuracy, presenting in particular a strong dry bias when a fraction of the pixel is occupied by forest (4.1%v/v RMSE,-3.1%v/v bias). An extension to the SMOS approach is proposed that accounts for the heterogeneity of vegetation optical depth within the SMOS pixel. The proposed approach is shown to significantly reduce the error in retrieved soil moisture (2.8%v/v RMSE, -0.3%v/v bias) in pixels characterised by a critical amount of forest (40-60%), at the limited cost of only a crude estimate of the optical depth of the forested area (better than 35% uncertainty). This study makes use of an unprecedented data set of airborne L-band observations and ground supporting data from the National Airborne Field Experiment 2005 (NAFE'05), which allowed accurate characterisation of the land surface heterogeneity over an area equivalent in size to a SMOS pixel.

Panciera, Rocco↗

Crater Morphometry and Scaling in Coarse, Rubble-Like Targets: Insights from Impact Experiments

Spacecraft images reveal that the asteroids Itokawa, Ryugu, and Bennu are covered with coarse, boulder-rich material [1–3]. Impactors that collide with these bodies encounter a target with extreme physical heterogeneity. Other bodies can also possess significant physical heterogeneity (e.g., megaregolith, layering, etc.). Such heterogeneities establish free surfaces and impedance contrasts that can affect shock propagation and attenuation. Therefore, such heterogeneities may also affect crater formation and excavation [4], melt generation [5–7] and crater scaling [4]. As described by [8,9], the extent to which target heterogeneity affects crater formation likely depends on how the length scale, d, of the heterogeneity (e.g., boulder size on a rubble-pile asteroid) compares to the width of the shock, w, generated by impact. Here we further test this hypothesis using impact experiments across a broad range of impact velocities and target grain sizes to systematically vary the ratio between the width of the shock and the diameter of target grains.

Daly, R. T.↗

Continuation of tailored composite structures of ordered staple thermoplastic material

The search for the cost effective composite structure has motivated the investigation of several approaches to develop composite structure from innovative material forms. Among the promising approaches is the conversion of a planar sheet to components of complex curvature through sheet forming or stretch forming. In both cases, the potential for material stretch in the fiber direction appears to offer a clear advantage in formability over continuous fiber systems. A framework was established which allows the simulation of the anisotropic mechanisms of deformation of long discontinuous fiber laminates wherein the matrix phase is a viscous fluid. Predictions for the effective viscosities of a hyper-anisotropic medium consisting of collimated, discontinuous fibers suspended in viscous matrix were extended to capture the characteristics of typical polymers including non-Newtonian behavior and temperature dependence. In addition, the influence of fiber misorientation was also modeled by compliance averaging to determine ensemble properties for a given orientation distribution. A design tool is presented for predicting the effect of material heterogeneity on the performance of curved composite beams such as those used in aircraft fuselage structures. Material heterogeneity can be induced during manufacturing processes such as sheet forming and stretch forming of thermoplastic composites. This heterogeneity can be introduced in the form of fiber realignment and spreading during the manufacturing process causing radial and tangential gradients in material properties. Two analysis procedures are used to solve the beam problems. The first method uses separate two-dimensional elasticity solutions for the stresses in the flange and web sections of the beam. The separate solutions are coupled by requiring that forces and displacements match section boundaries. The second method uses an approximate Rayleigh-Ritz technique to find the solutions for more complex beams. Analyses are performed for curved beams of various cross-sections loaded in pure bending and with a uniform distributed load. Preliminary results show that the geometry of the beam dictates the effect of heterogeneity on performance. The role of heterogeneity is larger in beams with a small average radius-to-depth ration, R/t, where R is the average radius of the beam and t is the difference between the inside and outside radii. Results of the anlysis are in the form of stresses and displacements and are compared to both mechanics of materials and numerical solutions obtained using finite element analysis.

Santare, Michael H.↗

Impact of stratospheric aircraft emissions on ozone: A two dimensional model study

Atmospheric perturbations caused by the emission of nitrogen oxides from a projected fleet of stratospheric aircraft are studied with a two dimensional chemistry, transport model. Photochemistry of the lower stratosphere, the region where these aircraft may fly, is now known to be influenced by heterogeneous reactions involving sulfuric acid aerosols. This study examines the sensitivity of the atmospheric effects of aircraft to heterogeneous reactions. Information of background aerosols based on the SAGE 2 measurements have been used in the parameterization of the heterogeneous conversion rates. It is found that heterogeneous reactions make the lower stratospheric ozone less sensitive to perturbations in the odd nitrogen level. The calculated reduction in global ozone due to NO(x) injection from a fleet of Mach 2.4 aircraft is 1.28 percent if gas phase reactions only are considered in the model, and 0.06 percent if heterogeneous reactions are included.

Natarajan, M.↗

Efficient Use of Distributed Systems for Scientific Applications

Distributed computing has been regarded as the future of high performance computing. Nationwide high speed networks such as vBNS are becoming widely available to interconnect high-speed computers, virtual environments, scientific instruments and large data sets. One of the major issues to be addressed with distributed systems is the development of computational tools that facilitate the efficient execution of parallel applications on such systems. These tools must exploit the heterogeneous resources (networks and compute nodes) in distributed systems. This paper presents a tool, called PART, which addresses this issue for mesh partitioning. PART takes advantage of the following heterogeneous system features: (1) processor speed; (2) number of processors; (3) local network performance; and (4) wide area network performance. Further, different finite element applications under consideration may have different computational complexities, different communication patterns, and different element types, which also must be taken into consideration when partitioning. PART uses parallel simulated annealing to partition the domain, taking into consideration network and processor heterogeneity. The results of using PART for an explicit finite element application executing on two IBM SPs (located at Argonne National Laboratory and the San Diego Supercomputer Center) indicate an increase in efficiency by up to 36% as compared to METIS, a widely used mesh partitioning tool. The input to METIS was modified to take into consideration heterogeneous processor performance; METIS does not take into consideration heterogeneous networks. The execution times for these applications were reduced by up to 30% as compared to METIS. These results are given in Figure 1 for four irregular meshes with number of elements ranging from 30,269 elements for the Barth5 mesh to 11,451 elements for the Barth4 mesh. Future work with PART entails using the tool with an integrated application requiring distributed systems. In particular this application, illustrated in the document entails an integration of finite element and fluid dynamic simulations to address the cooling of turbine blades of a gas turbine engine design. It is not uncommon to encounter high-temperature, film-cooled turbine airfoils with 1,000,000s of degrees of freedom. This results because of the complexity of the various components of the airfoils, requiring fine-grain meshing for accuracy. Additional information is contained in the original.

Taylor, Valerie↗

GCSS Cirrus Parcel Model Comparison Project

The Cirrus Parcel Model Comparison Project, a project of GCSS Working Group on Cirrus Cloud Systems (WG2), involves the systematic comparison of current models of ice crystal nucleation and growth for specified, typical, cirrus cloud environments. The goal of this project is to document and understand the factors resulting in significant inter-model differences. The intent is to foment research leading to model improvement and validation. In Phase 1 of the project reported here, simulated cirrus cloud microphysical properties are compared for situations of "warm" (-40 C) and "cold" (-60 C) cirrus subject to updrafts of 4, 20 and 100 cm/s, respectively. Five models participated. These models employ explicit microphysical schemes wherein the size distribution of each class of particles (aerosols and ice crystals) is resolved into bins. Simulations are made including both homogeneous and heterogeneous ice nucleation mechanisms. A single initial aerosol population of sulfuric acid particles is prescribed for all simulations. To isolate the treatment of the homogeneous freezing (of haze drops) nucleation process, the heterogeneous nucleation mechanism is disabled for a second parallel set of simulations. Qualitative agreement is found for the homogeneous-nucleation-only simulations, e.g., the number density of nucleated ice crystals increases with the strength of the prescribed updraft. However, non-negligible quantitative differences are found. Detailed analysis reveals that the homogeneous nucleation formulation, aerosol size, ice crystal growth rate (particularly the deposition coefficient), and water vapor uptake rate are critical components that lead to differences in predicted microphysics. Systematic bias exists between results based on a modified classical theory approach and models using an effective freezing temperature approach to the treatment of nucleation. Each approach is constrained by critical freezing data from laboratory studies, but each includes assumptions that can only be justified by further laboratory data. Consequently, it is not yet clear if the two approaches can be made consistent. Large haze particles may deviate considerably from equilibrium size in moderate to strong updrafts (20-100 cm/s) at -60 C when the commonly invoked equilibrium assumption is lifted. The resulting difference in particle-size-dependent solution concentration of haze particles may significantly affect the ice nucleation rate during the initial nucleation interval. The uptake rate for water vapor excess by ice crystals is another key component regulating the total number of nucleated ice crystals. This rate, the product of ice number concentration and ice crystal diffusional growth rate, which is sensitive to the deposition coefficient when ice particles are small, partially controls the peak nucleation rate achieved in an air parcel and the duration of the active nucleation time period. The effects of heterogeneous nucleation are most pronounced in weak updraft situations. Vapor competition by the nucleated (heterogeneous) ice crystals limits the achieved ice supersaturation and thus suppresses the contribution of homogeneous nucleation. Correspondingly, ice crystal number density is markedly reduced. Definitive laboratory and atmospheric benchmark data are needed for the heterogeneous nucleation process. Inter-model differences are correspondingly greater than in the case of the homogeneous nucleation process acting alone.

Lin, Ruei-Fong↗

Cirrus Parcel Model Comparison Project: The Critical Components to Simulate Cirrus Initiation Explicitly - Phase 1

The Cirrus Parcel Model Comparison Project, a project of the GCSS (GEWEX Cloud System Studies) Working Group on Cirrus Cloud Systems, involves the systematic comparison of current models of ice crystal nucleation and growth for specified, typical, cirrus cloud environments. In Phase I of the project reported here, simulated cirrus cloud microphysical properties are compared for situations of "warm" (40 C) and "cold" (-60 C) cirrus, both subject to updrafts of 4, 20 and 100 centimeters per second. Five models participated. The various models employ explicit microphysical schemes wherein the size distribution of each class of particles (aerosols and ice crystals) is resolved into bins or treated separately. Simulations are made including both the homogeneous and heterogeneous ice nucleation mechanisms. A single initial aerosol population of sulfuric acid particles is prescribed for all simulations. To isolate the treatment of the homogeneous freezing (of haze droplets) nucleation process, the heterogeneous nucleation mechanism is disabled for a second parallel set of simulations. Qualitative agreement is found for the homogeneous-nucleation- only simulations, e.g., the number density of nucleated ice crystals increases with the strength of the prescribed updraft. However, significant quantitative differences are found. Detailed analysis reveals that the homogeneous nucleation rate, haze particle solution concentration, and water vapor uptake rate by ice crystal growth (particularly as controlled by the deposition coefficient) are critical components that lead to differences in predicted microphysics. Systematic bias exists between results based on a modified classical theory approach and models using an effective freezing temperature approach to the treatment of nucleation. Each approach is constrained by critical freezing data from laboratory studies, but each includes assumptions that can only be justified by further laboratory research. Consequently, it is not yet clear if the two approaches can be made consistent. Large haze particles may deviate considerably from equilibrium size in moderate to strong updrafts (20-100 centimeters per second) at -60 C when the commonly invoked equilibrium assumption is lifted. The resulting difference in particle-size- dependent solution concentration of haze particles may significantly affect the ice particle formation rate during the initial nucleation interval. The uptake rate for water vapor excess by ice crystals is another key component regulating the total number of nucleated ice crystals. This rate, the product of particle number concentration and ice crystal diffusional growth rate, which is particularly sensitive to the deposition coefficient when ice particles are small, modulates the peak particle formation rate achieved in an air parcel and the duration of the active nucleation time period. The effects of heterogeneous nucleation are most pronounced in weak updraft situations. Vapor competition by the heterogeneously nucleated ice crystals may limit the achieved ice supersaturation and thus suppresses the contribution of homogeneous nucleation. Correspondingly, ice crystal number density is markedly reduced. Definitive laboratory and atmospheric benchmark data are needed for the heterogeneous nucleation process. Inter-model differences are correspondingly greater than in the case of the homogeneous nucleation process acting alone.

Lin, Ruei-Fong↗

Sensitivity Studies of Dust Ice Nuclei Effect on Cirrus Clouds with the Community Atmosphere Model CAM5

In this study the effect of dust aerosol on upper tropospheric cirrus clouds through heterogeneous ice nucleation is investigated in the Community Atmospheric Model version 5 (CAM5) with two ice nucleation parameterizations. Both parameterizations consider homogeneous and heterogeneous nucleation and the competition between the two mechanisms in cirrus clouds, but differ significantly in the number concentration of heterogeneous ice nuclei (IN) from dust. Heterogeneous nucleation on dust aerosol reduces the occurrence frequency of homogeneous nucleation and thus the ice crystal number concentration in the Northern Hemisphere (NH) cirrus clouds compared to simulations with pure homogeneous nucleation. Global and annual mean shortwave and longwave cloud forcing are reduced by up to 2.0+/-0.1Wm (sup−2) (1 uncertainty) and 2.4+/-0.1Wm (sup−2), respectively due to the presence of dust IN, with the net cloud forcing change of −0.40+/-0.20W m(sup−2). Comparison of model simulations with in situ aircraft data obtained in NH mid-latitudes suggests that homogeneous ice nucleation may play an important role in the ice nucleation at these regions with temperatures of 205-230 K. However, simulations overestimate observed ice crystal number concentrations in the tropical tropopause regions with temperatures of 190- 205 K, and overestimate the frequency of occurrence of high ice crystal number concentration (greater than 200 L(sup-1) and underestimate the frequency of low ice crystal number concentration (less than 30 L(sup-1) at NH mid-latitudes. These results highlight the importance of quantifying the number concentrations and properties of heterogeneous IN (including dust aerosol) in the upper troposphere from the global perspective.

Liu, Xiaohong↗

Evaluation of Satellite Remote Sensing Albedo Retrievals over the Ablation Area of the Southwestern Greenland Ice Sheet

MODerate resolution Imaging Spectroradiometer (MODIS) albedo products have been validated over spatially uniform, snow-covered areas of the Greenland ice sheet (GrIS) using the so-called single 'point-to-pixel' method. This study expands on this methodology by applying a 'multiple-point-to-pixel' method and examination of spatial autocorrelation (here using semivariogram analysis) by using in situ observations, high-resolution World- View-2 (WV-2) surface reflectances, and MODIS Collection V006 daily blue-sky albedo over a spatially heterogeneous surfaces in the lower ablation zone in southwest Greenland. Our results using 232 ground-based samples within two MODIS pixels, one being more spatial heterogeneous than the other, show little difference in accuracy among narrow and broad band albedos (except for Band 2). Within the more homogenous pixel area, in situ and MODIS albedos were very close (error varied from −4% to +7%) and within the range of ASD standard errors. The semivariogram analysis revealed that the minimum observational footprint needed for a spatially representative sample is 30 m. In contrast, over the more spatially heterogeneous surface pixel, a minimum footprint size was not quantifiable due to spatial autocorrelation, and far exceeds the effective resolution of the MODIS retrievals. Over the high spatial heterogeneity surface pixel, MODIS is lower than ground measurements by 4-7%, partly due to a known in situ undersampling of darker surfaces that often are impassable by foot (e.g., meltwater features and shadowing effects over crevasses). Despite the sampling issue, our analysis errors are very close to the stated general accuracy of the MODIS product of 5%. Thus, our study suggests that the MODIS albedo product performs well in a very heterogeneous, low-albedo, area of the ice sheet ablation zone. Furthermore, we demonstrate that single 'point-to-pixel' methods alone are insufficient in characterizing and validating the variation of surface albedo displayed in the lower ablation area. This is true because the distribution of in situ data deviations from MODIS albedo show a substantial range, with the average values for the 10th and 90th percentiles being −0.30 and 0.43 across all bands. Thus, if only single point is taken for ground validation, and is randomly selected from either distribution tails, the error would appear to be considerable. Given the need for multiple in-situ points, concurrent albedo measurements derived from existing AWSs, (low-flying vehicles (airborne or unmanned) and high-resolution imagery (WV-2)) are needed to resolve high sub-pixel variability in the ablation zone, and thus, further improve our characterization of Greenland's surface albedo.

Ablation zone↗

Limitations in Representation of Physical Processes Prevents Successful Simulation of PM2.5 During KORUS-AQ

High levels of fine particulate matter (PM2.5) pollution in East Asia often exceed local air quality standards. Observations from the Korea United States-Air Quality (KORUS-AQ) field campaign in May and June 2016 showed that development of extreme pollution (haze) occurred through a combination of long-range transport and favorable meteorological conditions that enhanced local production of PM2.5. Atmospheric models often have difficulty simulating PM2.5 chemical composition during haze, which is of concern for the development of successful control measures. We use observations from KORUS-AQ to examine the ability of the GEOS-Chem chemical transport model to simulate PM2.5 composition throughout the campaign and identify the mechanisms driving the pollution event. At the surface, the model underestimates sulfate by -64% but overestimates nitrate by +36%. The largest underestimate in sulfate occurs during the pollution event, where models typically struggle to generate elevated sulfate concentrations due to missing heterogeneous chemistry in aerosol liquid water in the polluted boundary layer. Hourly surface observations show that the model nitrate bias is driven by an overestimation of the nighttime peak. In the model, nitrate formation is limited by the supply of nitric acid, which is biased by +100% against aircraft observations. We hypothesize that this is due to a large missing sink, which we implement here as a factor of five increase in dry deposition. We show that the resulting increased deposition velocity is consistent with observations of total nitrate as a function of photochemical age. The model does not account for factors such as the urban heat island effect or the heterogeneity of the built-up urban landscape resulting in insufficient model turbulence and surface area over the study area that likely results in insufficient dry deposition. Other species such as NH3 could be similarly affected but were not measured during the campaign. Nighttime production of nitrate is driven by NO2 hydrolysis in the model, while observations show that unexpectedly elevated nighttime ozone (not present in the model) should result in N2O5 hydrolysis as the primary pathway. The model is unable to represent nighttime ozone due to an overly rapid collapse of the afternoon mixed layer and excessive titration by NO. We attribute this to missing nighttime heating driving deeper nocturnal mixing that would be expected to occur in a city like Seoul. This urban heating is not considered in air quality models run at large enough scales to treat both local chemistry and long-range transport. Key model failures in simulating nitrate, mainly overestimated daytime nitric acid, incorrect representation of nighttime chemistry, and an overly shallow and insufficiently turbulent nighttime mixed layer, exacerbate the model’s inability to simulate the buildup of PM2.5 during haze pollution. To address the underestimate in sulfate most evident during the haze event, heterogeneous aerosol uptake of SO2 is added to the model which previously only considered aqueous production of sulfate from SO2 in cloud water. Implementing a simple parameterization of this chemistry improves the model abundance of sulfate but degrades the SO2 simulation implying that emissions are underestimated. We find that improving model simulations of sulfate has direct relevance to determining local vs. transboundary contributions to PM2.5. During the haze pollution event, the inclusion of heterogeneous aerosol uptake of SO2 decreases the fraction of PM2.5 attributable to long-range transport from 66% to 54%. Locally-produced sulfate increased from 1% to 25% of locally-produced PM2.5, implying that local emissions controls would have a larger effect than previously thought. However, this additional uptake of SO2 is coupled to the model nitrate prediction which affects the aerosol liquid water abundance and chemistry driving sulfate-nitrate-ammonium partitioning. An additional simulation of the haze pollution with heterogeneous uptake of SO2 to aerosol and simple improvements to the model nitrate simulation results in 30% less sulfate due to 40% less nitrate and aerosol water, and results in an underestimate of sulfate during the haze event. Future studies need to better consider the impact of model physical processes such as dry deposition and nighttime boundary layer mixing on the simulation of nitrate and the effect of improved nitrate simulations on the overall simulation of secondary inorganic aerosol (sulfate+nitrate+ammonium) in East Asia. Foreign emissions are rapidly changing, increasing the need to understand the impact of local emissions on PM2.5 in South Korea to ensure continued air quality improvements.

Katherine R. Travis↗

Oxygen Isotopic Composition of Refractory Inclusions from the Miller Range (MIL) 090019 Carbonaceous Chondrite

Introduction: Primitive meteorites contain Calcium-Aluminum-rich Inclusions (CAIs)that preserve the records of the earliest times in the Solar System. CAIs record O isotopic variations with time and/or location from solar (16O-rich) to planetary (16O-poor) compositions [e.g., 1, 2].However, CAIs were also affected by parent body processing to varying degrees. The goal of this work was to identify the processes that resulted in the fine-scale spatial variations in oxygen isotopic compositions within the most ancient Solar System solids. We report oxygen isotopic composition of CAIs from MIL 090019 CO3.1 carbonaceous chondrite. MIL 090019 contains a high abundance of CAIs showing mineralogical as well as textural variations [3]. Methods: We analyzed 10 CAIs with different mineralogical assemblages, classified as corundum-bearing, grossite-bearing, hibonite-bearing, melilite-bearing, anorthite-bearing inclusions and amoeboid olivine aggregates(AOA’s). The mineralogical and petrological characterization of CAIs was performed using the JOEL Hyperprobe 8530 electron microprobe at NASA JSC. Oxygen isotopic imaging of the CAIs was done using the Cameca NanoSIMS 50L ion microprobe at NASA JSC. We followed the analytical protocol described in [4]. O-isotopic maps of the CAI were acquired by rastering a ~3 pA primary Cs+beam at 16 keVover an area of 20 ×20μm2over a period of ~7 hours. Negative secondary ions of 16O−, 17O−, 18O−, 28Si−, 24Mg16O−, 27Al16O−, and 40Ca16O−were simultaneously acquired using electron multiplier detectors at a mass resolving power sufficient to resolve the 16OH−interference from the 17O−peak, where the contribution of 16OH−was <0.1%.An electron flood gun was used to mitigate sample charging during the analyses. We used San Carlos olivine and Madagascar hibonite as isotopic standards to correct for the instrumental mass fractionation. The O-isotopic ratios were corrected for the quasi simultaneous arrival (QSA) effect and the detector dead time. All reported errors are 1 sigma. Results: The MIL 090019 carbonaceous chondrite hosts refractory inclusions varying in their mineralogies and textures[2, 3]. We found the oxygen isotopic compositions of mineral phases in MIL 090019 CAIs show large variations with Δ17O values varying from -27 to 0‰. A corundum-bearing inclusion records isotopic variations with corundum, melilite, perovskite, and anorthite showing Δ17Ovalues of -20 ±2.5, -10.2 ± 6.5, -10.7 ± 1.3,and -1.2 ± 3.2‰ respectively. A grossite-hibonite-bearing inclusion also records O-isotopic heterogeneity, where grossite, hibonite, Mg-rich spinel, and melilite show -3.7 ± 2.7, -22.9 ± 2.5, -21.4 ± 2.6, and -17.3 ± 2.8‰ respectively, whereas perovskite grains record a heterogenous isotopic composition ranging from -22.6 to -13.9‰. A hibonite-bearing inclusion shows average Δ17Ο values of -17.6 ± 5.3 ‰ and -15.1 ± 4.8 ‰ for hibonite and spinel respectively. The two perovskite grains also show relatively 16O-poor composition (Δ17Ο = -11.8 ± 6.4 ‰ and -8.6 ± 5.2 ‰). Anorthite has ~0.1 wt% Na2O and ~0.52 wt% FeO and a Δ17Ο value of -7.2 ± 5.0 ‰, whereas diopside records a Δ17Ο value of -3.5 ± 5.0 ‰. A hibonite-spinel-rich inclusion contains spinel and hibonite that are16O-rich (Δ17Οvalues -21.3 ± 2.8‰ and -19.4 ± 3.0‰, respectively). Melilite shows heterogenous oxygen isotopic composition, with an inner 16O-rich region and an outer 16O-poor region separated by a sharp boundary, with a bulk Δ17Ο= -16.2 ± 2.1‰, whereas two perovskite grains are relatively 16O-poor (Δ17Ο= -13.2 ± 3.7‰ and -7.1 ± 3.4‰).A spinel-rich inclusion contains perovskite grains with Δ17Οranging from -27.0 to -22.8 ‰and spinel with average Δ17Οvalue of -21.4 ± 2.8 ‰. In a melilite and spinel-rich CAI, perovskite, spinel, diopside, melilite, and Al-Ti-rich pyroxene record Δ17Ο values of-26.9 ± 4.6‰, -22.4 ± 4.3‰,-23.4 ± 4.9, -19.3 ± 3.4, and -19.0 ± 4.3 ‰ respectively. Olivine and spinel components of 2 AOA’s record a 16O-rich composition with Δ17O~-23‰. Discussion: The coexistence of 16O-rich and 16O-poor minerals within CAIs has been attributed to differing degrees of O isotopic exchange, reflecting a wide range of O diffusion rates[5]. However, the O isotopic heterogeneity observed in these CAIs cannot be explained by oxygen diffusion in nebular or parent body settings alone. The isotopic imaging of CAIs from the MIL 090019 meteorite shows variations in the O isotopic composition recorded in their mineral components, suggesting that these CAIs record the oxygen isotopic heterogeneity in the nebular gas from which they condensed. The coexistence of 16O-rich and 16O-poor mineral components in the CAI-forming region suggests that the gas in this region was not well-mixed and consisted of distinct 16O-rich and 16O-poor gaseous reservoirs.

P Mane↗

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↗

Characterizing the subglacial environment of lower Thwaites Glacier using radar modeling

Understanding the spatial heterogeneity beneath Thwaites Glacier, West Antarctica, is vital to projecting its impact on future sea levels. Radar-echo sounding (RES) is commonly used to infer subglacial conditions, but these data can be challenging to interpret. We assess basal heterogeneity across Thwaites Glacier by comparing RES returns to a radar backscattering simulator for over 400 km of RES data. The modeled variations in bed returned power exhibited a strong correlation with actual RES data in 40% of our simulated flight segments, which we consider evidence for a relatively homogeneous glacier bed. Other sites (40%) demonstrated improved fit quality when hydrology or substrate transitions were introduced in the bed material model. The remaining simulated segments (20%) were diagnosed as having more complex basal heterogeneity. The spatial distribution of complex heterogeneity appears to coincide with asymmetric patterns in the RES specularity content, which has been interpreted in previous studies as a signature for channelized hydrology. Conversely, the homogeneous substrate locations coincide with areas of fast-moving ice in western Thwaites. Our simulation method can isolate power variations induced by material heterogeneity vs topography, which is an important limitation of existing RES analysis methods.

Antarctic glaciology↗

On the Susceptibility of Cold Tropical Cirrus to Ice Nuclei Abundance

Numerical simulations of cirrus formation in the tropical tropopause layer (TTL) are used to evaluate the impact heterogeneous ice nuclei (IN) abundance on cold cloud microphysical properties. The model includes homogeneous and heterogeneous ice nucleation, deposition growthsublimation, and sedimentation. Reanalysis temperature and wind fields are used to force the simulations, with addition of high-frequency wave temperature variability based on in situ measurements. The model results are constrained by recent in situ observations of TTL cirrus and relative humidity, as well as satellite measurements of cirrus bulk properties and occurrence frequencies. Temperature variability driven by high-frequency waves has a dominant influence on TTL cirrus microphysical properties and occurrence frequencies, and inclusion of these waves is required to produce agreement between the simulated and observed cloud properties as well as the abundance of TTL cirrus. With a composition-independent supersaturation threshold for homogeneous freezing of aqueous aerosols, the model produces excessive ice concentrations compared with in situ observations. Inclusion of relatively numerous heterogeneous nuclei ( 100L) in the simulations improves the agreement with observed ice concentrations. However, when IN contribute significantly to TTL cirrus ice nucleation, the occurrence frequency of large supersaturations with respect to ice is less than indicated by in situ measurements. The simulated TTL cirrus extinction statistics agree with observations (within uncertainties) from both in situ and remote-sensing measurements. We find that the sensitivity of TTL cirrus extinction and ice water content statistics to heterogeneous ice nuclei abundance is relatively weak. Likewise, the simulated occurrence frequencies of TTL cirrus are insensitive to ice nuclei abundance, both in terms of cloud frequency height distribution and regional distribution throughout the tropics.

occurences↗

Effects of Topography on Tropical Forest Structure Depend on Climate Context

Topography affects abiotic conditions which can influence the structure, function, and dynamics of ecological communities. An increasing number of studies have demonstrated biological consequences of fine-scale topographic heterogeneity but we have a limited understanding of how We merged high-resolution (1 sq. meter) data on topography and canopy height derived from airborne lidar with ground-based data from 15 forest plots in Puerto Rico distributed along a precipitation gradient spanning ca. 800 to 3,500 mm yr(exp -1). Ground-based data included species composition, estimated above-ground biomass (AGB), and two key functional traits (wood density and leaf mass per area, LMA) that reflect resource-use strategies and a trade-off between hydraulic safety and hydraulic efficiency. We used hierarchical Bayesian models to evaluate how the interaction between topography climate is related to metrics of forest structure (i.e., canopy height and AGB), as well as taxonomic and functional alpha- and beta-diversity. Fine-scale topography (characterized with the topographic wetness index, TWI) significantly affected forest structure and the strength (and in some cases direction) of these effects varied across the precipitation gradient. In all plots, canopy height increased with topographic wetness but the effect was much stronger in dry compared to wet forest plots. In dry forest plots, topographically wetter microsites also had higher levels of AGB but in wet forest plots, topographically drier microsites had higher AGB. Fine-scale topography influenced functional composition but had only weak or non-significant effects on taxonomic and functional alpha- and beta-diversity. For instance, community-weighted wood density followed a similar pattern to AGB across plots. We also found a marginally significant association between variation of wood density and topographic heterogeneity that depended on climate context. Synthesis: The effects of fine-scale topographic heterogeneity on tropical forest structure and composition depend on the climate context. Our study demonstrates how a stronger integration of topographic heterogeneity across precipitation gradients could improve estimates of forest structure and biomass, and may provide insight to the ways that topography might mediate species responses to drought and climate change.

Tropical dry forest↗