Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “polar format processing”

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.

664 records · Page 37

Highly Supercooled Riming and Unusual Triple-Frequency Radar Signatures Over McMurdo Station, Antarctica

Riming of ice crystals by supercooled water droplets is an efficient ice growth process, but its basic properties are still poorly known. While it has been shown to contribute significantly to surface precipitation at mid-latitudes, little is known about its occurrence at high latitudes. In Antarctica, two competing effects can influence the occurrence of riming: (i) the scarcity of supercooled liquid water clouds due to the extremely low tropospheric temperatures and (ii) the low aerosol concentration, which may lead to the formation of fewer and larger supercooled drops potentially resulting in an enhanced riming efficiency. In this work, by exploiting the deployment of an unprecedented number of multiwavelength remote sensing systems (including triple-frequency radar measurements) in West Antarctica, during the Atmospheric Radiation Measurements West Antarctic Radiation Experiment (AWARE) field campaign, we evaluate the riming incidence at McMurdo Station and find that riming occurs at much lower temperatures when compared to previous results in the mid-latitudes. This suggests the possible occurrence of a common atmospheric state over Antarctica that includes a rather stable atmosphere inhibiting turbulent mixing, and a high riming efficiency driven by large cloud droplets. We then focus on a peculiar case study featuring a persistent layer with a particularly pronounced riming signature in triple-frequency radar data but only a relatively modest amount of supercooled liquid water. In-depth analysis of the radar observations suggests that such signatures can only be explained by the combined effects of moderately rimed aggregates or similarly shaped florid polycrystals and a narrow particle size distribution (PSD). Simulations of this case study performed with a one-dimensional bin model indicate that similar triple-frequency radar observations can be reproduced when narrow PSDs are simulated. Such narrow PSDs can in turn be explained by two key factors: (i) the presence of a shallow homogeneous droplet or humidified aerosol freezing layer aloft seeding an underlying supercooled liquid layer, and (ii) the absence of turbulent mixing throughout a stable polar atmosphere that sustains narrow PSDs, as hydrometeors grow from the nucleation region aloft to ice particles of several millimeters in size, by vapor deposition and then riming.

Riming↗

The Structure and Dynamics of the Subparsec Jet in M87 Based on 50 VLBA Observations over 17 Years at 43 GHz

The central radio source in M87 provides the best opportunity to study jet formation because it has a large angular size for the gravitational radius of the black hole and has a bright jet that is well resolved by very long baseline interferometry observations. We present intensive monitoring observations from 2007 and 2008, plus roughly annual observations that span 17 years, all made with the the Very Long Baseline Array at 43 GHz with a resolution of about 30 by 60 R S . Our high dynamic range images clearly show the wide opening angle structure and the counterjet. The jet and counterjet are nearly symmetric in the inner 1.5 mas (0.12 pc in projection), with both being edge brightened. Both show deviations from parabolic shape in the form of an initial rapid expansion and subsequent contraction followed by further rapid expansion and, beyond the visible counterjet, subsequent collimation. Proper motions and counterjet/jet intensity ratios both indicate acceleration from apparent speeds of ≲ 0.5c to ≳ 2c in the inner ~2 mas and suggest a helical flow. The jet displays a sideways shift with an approximately 8–10 yr quasi-periodicity. Here, the shift propagates outward nonballistically and significantly more slowly than the flow speed revealed by the fastest-moving components. Polarization data show a systematic structure with magnetic field vectors that suggest a toroidal field close to the core.

79 ASTRONOMY AND ASTROPHYSICS↗

Integrated Geophysical Measurements for Bioremediation Monitoring: Combining Spectral Induced Polarization, Nuclear Magnetic Resonance and Magnetic Methods

This documents contains the final report for the project "Integrated Geophysical Measurements for Bioremediation Monitoring: Combining Spectral Induced Polarization, Nuclear Magnetic Resonance and Magnetic Methods" (DE-SC0007049) Executive Summary: Our research aimed to develop borehole measurement techniques capable of monitoring subsurface processes, such as changes in pore geometry and iron/sulfur geochemistry, associated with remediation of heavy metals and radionuclides. Previous work has demonstrated that geophysical method spectral induced polarization (SIP) can be used to assess subsurface contaminant remediation; however, SIP signals can be generated from multiple sources limiting their interpretation value. Integrating multiple geophysical methods, such as nuclear magnetic resonance (NMR) and magnetic susceptibility (MS), with SIP, could reduce the ambiguity of interpretation that might result from a single method. Our research efforts entails combining measurements from these methods, each sensitive to different mineral forms and/or mineral-fluid interfaces, providing better constraints on changes in subsurface biogeochemical processes and pore geometries significantly improving our understanding of processes impacting contaminant remediation. The Rifle Integrated Field Research Challenge (IFRC) site was used as a test location for our measurements. The Rifle IFRC site is located at a former uranium ore-processing facility in Rifle, Colorado. Leachate from spent mill tailings has resulted in residual uranium contamination of both groundwater and sediments within the local aquifer. Studies at the site include an ongoing acetate amendment strategy, native microbial populations are stimulated by introduction of carbon intended to alter redox conditions and immobilize uranium. To test the geophysical methods in the field, NMR and MS logging measurements were collected before, during, and after acetate amendment. Next, laboratory NMR, MS, and SIP measurements were collected on columns of Rifle sediments during acetate amendment. The laboratory experiments were designed to simulate the field experiments; changes in geophysical signals were expected to correlate with changes in redox conditions and iron speciation. Field MS logging measurements revealed vertically stratified magnetic mineralization, likely the result of detrital magnetic fraction within the bulk alluvium. Little to no change was observed in the MS data suggesting negligible production of magnetic phases (e.g. magnetite, pyrrhotite) as a result of sulfidogenesis. Borehole NMR measurements contained high levels of noise contamination requiring significant signal processing, and analysis suggests that any changes may be difficult to differentiate from simultaneous changes in water content. Laboratory MS and NMR measurements remained relatively stable throughout the course of the acetate amendment experiment, consistent with field measurements. However, SIP measurements changed during the acetate amendment associated with the formation of iron-sulfide mineral phases; a finding that is consistent with chemical analysis of the solid phase materials in the columns.

54 ENVIRONMENTAL SCIENCES↗

GRB and blazar jets shining through their stripes

Black hole-driven relativistic astrophysical jets, such as blazars and gamma-ray bursts (GRBs), are powerful sources of electromagnetic radiation. Their emission is powered by some energy dissipation and particle acceleration mechanism operating in situ at large distances from the black hole. We propose that the formation of the dissipative structures in the jet is controlled by the time variability of the accretion disc. We argue that the open magnetic field lines through the black hole, which drive a strongly magnetized jet, may have their polarity reversing over time-scales related to the growth of the magnetorotational dynamo in the disc. Consequently, the jet is characterized by an alternating toroidal field polarity along its propagation axis, i.e. it is a ‘striped jet’. Magnetic reconnection in the current sheets that form between the stripes dissipates the alternating field energy and thus powers further jet acceleration. Here, we consider a jet with a broad distribution of stripe widths l > l min , above a dominant scale l min . We find that the bulk acceleration of the jet, driven by the annihilation of the stripes, is very gradual. The dissipation rate peaks at a distance $z_{\rm peak}\sim 10^6 R_\mathrm{ g}\, (\Gamma _\infty /30)^2\, (l_{\rm min}/1000R_\mathrm{ g})$| , where R g is the black hole’s gravitational radius and Γ ∞ the jet’s asymptotic Lorentz factor, and exhibits a very broad plateau extending by ~4–5 orders of magnitude in distance. The prolonged energy dissipation accounts for the flat-to-inverted long-wavelength spectra commonly observed in jets. The model can also account for the broad range of flaring time-scales of blazars and the fact that their bulk acceleration appears to continue out to ~100 pc scales. In GRB jets, the model predicts comparable power for the photopsheric and Thomson-thin emission components.

79 ASTRONOMY AND ASTROPHYSICS↗

Clementine High Resolution Camera Mosaicking Project

This report constitutes the final report for NASA Contract NASW-5054. This project processed Clementine I high resolution images of the Moon, mosaicked these images together, and created a 22-disk set of compact disk read-only memory (CD-ROM) volumes. The mosaics were produced through semi-automated registration and calibration of the high resolution (HiRes) camera's data against the geometrically and photometrically controlled Ultraviolet/Visible (UV/Vis) Basemap Mosaic produced by the US Geological Survey (USGS). The HiRes mosaics were compiled from non-uniformity corrected, 750 nanometer ("D") filter high resolution nadir-looking observations. The images were spatially warped using the sinusoidal equal-area projection at a scale of 20 m/pixel for sub-polar mosaics (below 80 deg. latitude) and using the stereographic projection at a scale of 30 m/pixel for polar mosaics. Only images with emission angles less than approximately 50 were used. Images from non-mapping cross-track slews, which tended to have large SPICE errors, were generally omitted. The locations of the resulting image population were found to be offset from the UV/Vis basemap by up to 13 km (0.4 deg.). Geometric control was taken from the 100 m/pixel global and 150 m/pixel polar USGS Clementine Basemap Mosaics compiled from the 750 nm Ultraviolet/Visible Clementine imaging system. Radiometric calibration was achieved by removing the image nonuniformity dominated by the HiRes system's light intensifier. Also provided are offset and scale factors, achieved by a fit of the HiRes data to the corresponding photometrically calibrated UV/Vis basemap, that approximately transform the 8-bit HiRes data to photometric units. The sub-polar mosaics are divided into tiles that cover approximately 1.75 deg. of latitude and span the longitude range of the mosaicked frames. Images from a given orbit are map projected using the orbit's nominal central latitude. Polar mosaics are tiled into squares 2250 pixels on a side, which spans approximately 2.2 deg. Two mosaics are provided for each pole: one corresponding to data acquired while periapsis was in the south, the other while periapsis was in the north. The CD-ROMs also contain ancillary data files that support the HiRes mosaic. These files include browse images with UV/Vis context stored in a Joint Photographic Experts Group (JPEG) format, index files ('imgindx.tab' and 'srcindx.tab') that tabulate the contents of the CD, and documentation files.

Source record↗

Triplet-pair spin signatures from macroscopically aligned heteroacenes in an oriented single crystal

The photo-driven process of singlet fission generates coupled triplet pairs (TT) with fundamentally intriguing and potentially useful properties. The quintet 5 TT 0 sublevel is particularly interesting for quantum information because it is highly entangled, is addressable with microwave pulses, and could be detected using optical techniques. Previous theoretical work on a model Hamiltonian and nonadiabatic transition theory, called the JDE model, has determined that this sublevel can be selectively populated if certain conditions are met. Among the most challenging, the molecules within the dimer undergoing singlet fission must have their principal magnetic axes parallel to one another and to an applied Zeeman field. Here, we present time-resolved electron paramagnetic resonance (TR-EPR) spectroscopy of a single crystal sample of a tetracenethiophene compound featuring arrays of dimers aligned in this manner, which were mounted so that the orientation of the field relative to the molecular axes could be controlled. The observed spin sublevel populations in the paired TT and unpaired (T+T) triplets are consistent with predictions from the JDE model, including preferential 5 TT 0 formation at z ‖ B 0 , with one caveat—two 5 TT spin sublevels have little to no population. This may be due to crossings between the 5 TT and 3 TT manifolds in the field range investigated by TR-EPR, consistent with the intertriplet exchange energy determined by monitoring photoluminescence at varying magnetic fields.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Measurement report: Role of organic coating and chemical composition on ice nucleation potential of atmospheric particles in European Arctic

Understanding the ice nucleation (IN) potential of Arctic aerosols is critical for predicting their influence on cloud formation and water cycles in this vulnerable region. This study investigates the role of particle composition, organic coatings, and aerosol sources in modulating ice nucleating particle (INPs) abundance across five aerosol samples collected at the Gruvebadet Observatory Station in Ny-Ålesund, Svalbard. The IN potential of Arctic aerosol particles was studied by investigating chemical, morphological, and INP abundance measurements. Single-particle analyses revealed distinct differences in mixing state, organic volume fraction (OVF), and organic coating morphology across samples. OVF distributions were linked to particle origin, with marine-influenced Na-rich particles often exhibiting thin organic coatings, while long-range transported particles showed thicker organic coatings. Biogenic contributions, though variable, were linked to heat-sensitive INPs, suggesting a role for labile biological macromolecules under certain meteorological conditions. Spearman rank correlation analysis between particle composition and immersion-mode INP concentrations at two freezing temperatures indicated that organic-rich and Na-rich particles were positively associated with enhanced INP abundance. However, discrepancies in INP abundance were observed for particles with thicker organic coatings, where the morphological configuration of the organic material may play a role. The results highlight that Arctic INP variability is governed not only by chemical composition but also by the morphological configuration of organic material, which can either enhance or inhibit ice nucleation depending on its abundance, distribution, thickness, and mixing state. These findings underscore the combined influence of source regions, atmospheric processing, and organic–inorganic interactions in shaping Arctic aerosol freezing behavior.

Lata, Nurun Nahar [Pacific Northwest National Labo↗

Formation and fate of freshwater on an ice floe in the Central Arctic

The melt of snow and sea ice during the Arctic summer is a significant source of relatively fresh meltwater. The fate of this freshwater, whether in surface melt ponds or thin layers underneath the ice and in leads, impacts atmosphere–ice–ocean interactions and their subsequent coupled evolution. Here, we combine analyses of datasets from the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition (June–July 2020) for a process study on the formation and fate of sea ice freshwater on ice floes in the Central Arctic. Our freshwater budget analyses suggest that a relatively high fraction (58 %) is derived from surface melt. Additionally, the contribution from stored precipitation (snowmelt) outweighs by 5 times the input from in situ summer precipitation (rain). The magnitude and rate of local meltwater production are remarkably similar to those observed on the prior Surface Heat Budget of the Arctic Ocean (SHEBA) campaign, where the cumulative summer freshwater production totaled around 1 m during both. A relatively small fraction (10 %) of freshwater from melt remains in ponds, which is higher on more deformed second-year ice (SYI) compared to first-year ice (FYI) later in the summer. Most meltwater drains laterally and vertically, with vertical drainage enabling storage of freshwater internally in the ice by freshening brine channels. In the upper ocean, freshwater can accumulate in transient meltwater layers on the order of 0.1 to 1 m thick in leads and under the ice. The presence of such layers substantially impacts the coupled system by reducing bottom melt and allowing false bottom growth; reducing heat, nutrient, and gas exchange; and influencing ecosystem productivity. Regardless, the majority fraction of freshwater from melt is inferred to be ultimately incorporated into the upper ocean (75 %) or stored internally in the ice (14 %). Terms such as the annual sea ice freshwater production and meltwater storage in ponds could be used in future work as diagnostics for global climate and process models. For example, the range of values from the CESM2 climate model roughly encapsulate the observed total freshwater production, while storage in melt ponds is underestimated by about 50 %, suggesting pond drainage terms as a key process for investigation.

54 ENVIRONMENTAL SCIENCES↗

Landsat at 45: How it Changed the Way We See the Earth

On October 24, 1946, more than 10 years before the launch of the first artificial satellite Sputnik, scientists at the White Sands Missile Range in New Mexico placed a camera on top of a captured German V-2 ballistic missile. As the rocket flew to an altitude of about 65 miles - just above the generally recognized border of outer space - the 35-millimeter motion picture camera snapped a frame every one and a half seconds. Minutes later, the missile came crashing back down and slammed into the ground at more than 340 mph, but the film survived and gave us our first glimpse of Earth from space. Earth Resources Technology Satellite aka Landsat It was images like those first grainy black and white pictures and later those taken by America's first astronauts in the 1960's that inspired the development of the Earth Resources Technology Satellite (ERTS). From the unique vantage point of space, we could now observe Earth using a variety of different instruments to monitor changes over time. The ERTS-1 satellite, wisely renamed Landsat-1, was launched aboard a Delta rocket on July 23, 1972, into a Sun-synchronous polar orbit at an altitude of about 560 miles. In this unique orbit, Landsat could observe the same point on the Earth every 18 days, always with the same solar illumination, allowing for precise monitoring of changes on the ground over time. Landsat-1, derived from the highly successful Nimbus weather satellites, carried two instruments that allowed it to take images not only in visible light but also in infrared, well-suited to track changes in vegetation over time. Designed to last only one year, Landsat-1 actually operated for nearly three years, by which time it had been joined in space by Landsat-2, a near identical copy of the original. Since then, ever more sophisticated instruments were flown aboard Landsat-3 through -8, with Landsat-9 planned for launch in 2020, acquiring millions of images of Earth over more than four decades. At first, images from Landsat were processed by NASA and hardcopies sold to users for a fee, a somewhat tedious process. Since 2008, images have been made available to all interested users by the US Geological Survey (USGS) at no cost via the Internet in near real time. So, how can Landsat help? In short, Landsat looks. And looks. And proves Yogi right. Space-based images from Landsat and other similar satellites offer a unique and critical capability to observe land use over time by providing repetitive observations of the Earth otherwise unavailable. The data provided by the images can be used by scientists and politicians to inform wise decisions in areas such as agriculture, climate, ecosystems and biodiversity, energy, forest management, human health, fire, natural disasters, urban growth and water management. This overview article doesn't allow for examples from each of these disciplines, but details can be found at the following website: https://landsat.gsfc.nasa.gov/how_landsat_helps/. One third of the US economy is influenced by climate, weather and natural hazards, providing strong economic incentives to sustain a healthy space-based Earth observation program. One example, however, may be illustrative of how Landsat and other space-based observations can be helpful in monitoring and documenting some effects of global climate change. It had been noted since the 1970's that permafrost in subarctic areas like Siberia was melting at an accelerating rate. It also became apparent that this led to the formation of hundreds of melt lakes and the liberation of large quantities of methane, a gas that in the short-term has more potent greenhouse effects than carbon dioxide. The methane turned the water in these lakes blue, making them easy to track over time by satellite. Additionally, the newly released methane has been carbon dated to tens of thousands of years ago, meaning that it had remained frozen since the last Ice Age and therefore hadn't been accounted for in models of the Earth's overall carbon balance. Release of sufficient methane by an increasingly warming climate could actually cause a positive feedback loop in global warming, melting more permafrost and releasing yet more methane. Further monitoring by Landsat and other platforms will reveal whether this process is reversible or whether we've passed the tipping point. In addition to the highly successful Landsat series of satellites, NASA and other agencies such as the National Oceanic and Atmospheric Administration (NOAA) operate a fleet of other Earth observing platforms, many with more specific research goals such as monitoring sea ice levels or atmospheric carbon content. Several instruments aboard the International Space Station also contribute to this overall effort to better understand short- and long-term changes to the Earth. Many of these missions are guided by the 2009 Decadal Survey published by the National Research Council of the National Academy of Sciences. Europe, Russia and China all see the value of space-based Earth observation by deploying their own fleet of satellites.

Uri, John↗

Aldol Condensation and Esterification over Ti-Substituted *BEA Zeolite: Mechanisms and Effects of Pore Hydrophobicity

Aldol condensation and esterification reactions provide paths to upgrade ethanol and acetaldehyde to higher- value molecules useful as fuels or intermediates for the synthesis of polymers. Transition-metal-substituted BEA zeolites (M-BEA) catalyze these reactions; however, the mechanisms for these processes in M-BEA and the effects of incidental or purposefully included silanol groups are not reported. Here, we combine kinetic and spectroscopic measurements obtained during catalytic reactions of acetaldehyde (CH 3 CHO), ethanol (C 2 H 5 OH), and hydrogen (H 2 ) mixtures over a series of Ti-BEA catalysts that possess a known range of silanol group densities to examine the kinetic relevance of intervening steps and the impact of silanol groups on catalytic rates. Across all Ti-BEA, rates for aldol condensation and esterification increase with the pressure of CH 3 CHO; however, C 2 H 5 OH and H 2 O weakly inhibit the rates of these reactions. The substitution of CD 3 CDO for CH 3 CHO decreases aldol condensation rates slightly (~10%) but leads to greater esterification rates (2- to 5-fold). The kinetic isotope effects together with the measured dependence of rates on reactant pressures suggest that aldol condensation and esterification occur on unoccupied Ti sites and involve multiple kinetically relevant steps. CH 3 CHO deprotonates irreversibly, and the kinetically relevant nucleophilic attack of the enolate to CH 3 CHO* (i.e., adsorbed CH 3 CHO on Ti sites) leads to aldol products, while the nucleophilic attack of the enolate to C 2 H 5 OH* gives esters. Selectivities toward aldol condensation increase with the ratio of CH 3 CHO to C 2 H 5 OH pressure and with increases in the silanol density of the as-synthesized Ti-BEA. During catalysis, in situ infrared spectroscopy demonstrates that these silanol groups react with C 2 H 5 OH to form ethoxysilane groups (i.e., SiOC 2 H 5 ) that modify the polarity of the environment near Ti active sites. As initial silanol densities increase, steady-state turnover rates for aldol condensation and esterification increase by factors of 5 and 2, respectively. The changes in rates and selectivities among Ti-BEA catalysts likely reflect changes in excess free energies of transition states for enolization and nucleophilic attack of the enolate to adsorbed coreactants. The differences in excess stability report on the interactions among reactive intermediates at framework Ti atoms and the ethoxysilane and remaining silanol groups present. Furthermore, the in situ modification of these pore environments confers changes in the stability of reactive species in a manner that contradicts intuition when considering the initial state of the catalyst but can be reconciled after accounting for the formation of persistent alkoxy surface moieties in the pores.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

GC/MS Method Development for Separating Lunar Volatile Ice Simulant Headspace Gases

Various investigators propose the lunar surface contains widely distributed volatiles, especially water- like species, i.e. OH and H2O. Surface volatiles are theorized to exist as a hydrated regolith layer, concentrated in extremely cold polar permanently shadowed regions (PSR), and/or solar wind implantation reservoirs in lunar glasses. The proposed sources of lunar surface volatiles range from cometary impacts, solar wind, or a supply present during moon formation. Future Artemis missions aim to collect and return the samples containing volatiles collected near lunar polar craters or PSRs. We, as advanced curation scientists, are responsible for developing techniques and methodologies for preserving returned sample integrity as much as possible. Pristine volatile-bearing samples are invaluable to the scientific community seeking to unravel the history of the solar system. Realistically, a sample will experience alteration during collection, transportation back to earth, and storage. The Planetary Exploration and Astromaterials Research Lab (PEARL) seeks to understand temperature and pressure effects on high-fidelity volatile-containing regolith simulants, the foundation for the future of cold curation. This abstract outlines the separation, identification, and quantification of headspace gases over volatile ice feed stock material using gas chromatography/mass spectrometry (GC/MS). Preliminary objectives concentrated on sample handling, reproducibility, and understanding the elution characteristics for each analyte. Initial GC/MS method development experiments utilized diluted static headspace sample preparation. Diluted samples were used because sampling headspace gases directly from a vial containing liquid analyte resulted in overloading of the column and detector. Overloading is evident based on chromatogram peak shapes and instrument contamination, or carry over, between experiments. A mixture of three alcohols were used for a majority of the sample handling and reproducibility studies. Reproducibility was tested via multiple users, calibration curves, and check standards. Stock solutions of condensed lunar volatile analytes included methanol, ammonia in methanol, hydrogen sulfide in water, and an equal volume mixture of methanol, ethanol, and isopropanol. Current samples use room air as the headspace sample matrix, however future experiments will incorporate an inert purge gas, such as argon or nitrogen. Three mL of each analyte solution were capped in separate 20 mL crimp top GC vials. Dilutions were carried out by removing an aliquot of headspace gases with a calibrated 1 mL gastight syringe and immediately transferring to a 20 mL capped crimp top vial. The GC/MS is a Thermo Fisher Trace 1310/ISQ 7000 with a TriPlus RSH autosampler and split/splitless injector module. The experiments outlined in this abstract use the following hardware: a 2.5 mL gastight headspace syringe tool, 1 mm ID x 78.5 mm length ultra-inert straight injection liner, and a TG-BondQ 30 m × 0.32 mm × 10 μm column. Various parameters, such as hardware selection and the temperature, pressure, and split ratio set points, continue to evolve as the overall experiment is refined. Diluted headspace chromatograms were collected for the individual stock solutions. Retention times, peak shapes, and mass spectra were evaluated and added to the data processing method for each molecule of interest. Figure 1 shows the total ion chromatograms for the three major lunar volatile simulant stock solutions: methanol, 7 N ammonia in methanol, and 0.4% hydrogen sulfide in water. Tailing peak shapes for ammonia (2.98 min rt) and water (4.06 min rt) indicate the molecules are not properly eluting from the selected column with the current separation method. Additionally, hydrogen sulfide and ammonia have overlapping peak windows, which could impact quantification. Ongoing experiments aim to address the peak shape and overlapping via the separation method and hardware selection. Sample preparation reproducibility experiments used stock solution containing equal volumes of a non- interactive mixture of methanol, ethanol, and isopropanol. Mass spectrum ion traces were used to identify and quantify all three alcohols. Peaks were automatically detected, identified, and integrated through the mass spectra detection and processing parameters. Calibration response curves and check standards were used to evaluate the validity of the sample preparation procedure. Figure 2 shows the methanol chromatogram peak area versus total headspace dilution volume transferred from the alcohol mixture vial. The calibration response curves and check standards validate sample preparation procedure. Continuing data analysis efforts are working towards correlating the peak area and instrument response factor to the headspace analyte concentration and condensed phase composition. Static headspace gas chromatography theory relies on Dalton’s law, Raoult’s law, Henry’s Law, and the Kolb and Ettre equation to associate peak area to the analyte composition in a non-ideal solution. Equation 1 is a simplified expression derived from the aforementioned theories. Future experiments involve liquid injections of the individual stock solutions, liquid and headspace analysis of various stock solution combinations, and the addition of regolith simulants to the mixtures. Temperature is another variable expected to affect reaction rates and will be explored.

Cecilia L. Amick↗

Final DOE-ASR Report for the Project “Advancing our Understanding and the Remote Sensing of Ice Clouds”

This project has evolved during its execution, and what follows are the key project findings. This project has arguably provided the first global view of how cirrus cloud (defined as having cloud base temperature T < 235 K) nucleation physics (evaluated through satellite retrievals of ice particle number concentration Ni, effective diameter De and ice water content IWC) evolves with the seasons for a given temperature, latitude zone and surface type (e.g. ocean vs. land), based on a new satellite remote sensing method developed for this project. The retrieval method is unique in that it is very sensitive to the small ice crystals that govern the number concentration Ni, allowing Ni to be retrieved. The method currently samples single-layer cirrus clouds having visible optical depth ranging from about 0.3 to 3.0, using co-located observations from the Infrared Imaging Radiometer (IIR) and from the CALIOP (Cloud and Aerosol Lidar with Orthogonal Polarization) lidar aboard the CALIPSO (Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation) polar orbiting satellite, employing IIR channels at 10.6 μm and 12.05 μm. Retrievals of Ni are primarily used to estimate the cirrus cloud formation mechanism; that is, either homo- or heterogeneous ice nucleation (henceforth hom and het). This is possible since, in general, hom produces more than an order of magnitude more ice crystals than does het. Thus the retrievals provide insight on how these mechanisms change with the seasons for a given latitude zone or region, based on the years 2008 and 2013. Using a conservative criterion for hom cirrus, on average, the sampled cirrus clouds formed through hom occur about 43% of the time in the Arctic and 50% of the time in the Antarctic, and during winter at mid-latitudes in the Northern Hemisphere, hom cirrus occur 37% of the time. Elsewhere (and during other seasons in the Northern Hemisphere mid-latitudes), this hom cirrus fraction is lower, and it is lowest in the tropics. Thus, the microphysical properties of cirrus clouds in the Polar Regions are much different than they are in the tropics; something unknown prior to this study. Moreover, the frequency of cirrus cloud occurrence in the Polar Regions varies strongly with season, peaking during winter in the Arctic and during spring in the Antarctic. Considering these seasonal changes in microphysics and inferred cloud coverage, this leads us to speculate that the buildup of Arctic cirrus during winter may significantly contribute to tropospheric heating in that region, possibly affecting winter jet-stream dynamics and mid-latitude weather patterns through the thermal-wind balance relationship. This cirrus cloud research provides essential guidance for realistically representing cirrus clouds in climate models; guidance previously unavailable. For example, mid-latitude hom cirrus were widespread during winter over or nearby mountainous terrain, evidently due to mountain-induced waves that produce strong updrafts at cirrus cloud levels. The treatment of turbulent mountain stress and gravity waves will likely need to be improved in climate models in order to adequately represent cirrus clouds outside the tropics. Another goal of this project was to develop a ground-based 94-GHz radar retrieval for winter snowstorms, based on (1) an improved analytical framework describing the interaction of radiation from radar with snowfall and (2) the development of a steady-state snow growth model that predicts the height-evolution of the ice particle size distribution through ice particle growth by vapor diffusion, aggregation and riming (i.e. the growth of snow through collisions with supercooled cloud droplets). Although activities (1) and (2) were completed, there was insufficient time to test and finalize the radar retrieval scheme. However, activity (2) provided a new method for relating ice particle mass “m” and projected area “A” to the ice particle maximum dimension “D”. The ice cloud microphysical processes (which determine ice cloud radiative properties) in climate models are parameterized in terms of these m-D and A-D relationships. By improving these relationships, the ice cloud radiative properties in Community Atmosphere Model version 5, or CAM5 (an atmosphere global climate model, or GCM) were improved. Student funding from the University of Nevada, Reno, was combined with funds from this project to conduct some basic research on the mechanism of the North American monsoon, or NAM. Federal research on the NAM has dwindled since 2006, but atmospheric soundings taken during research vessel cruises in the Gulf of California (GC) during the North American Monsoon Experiment (NAME) were used to reveal a likely mechanism that explains the relationship between an intrusion of tropical warm water into the GC during late spring-early summer and the onset of relatively heavy NAM rainfall in northwest Mexico and the southwestern United States. These soundings, combined with reanalysis data, satellite sea surface temperatures and satellite measurements of outgoing longwave radiation were used to develop and provide evidence for a planetary-scale NAM mechanism. As far as we know, no other physical explanation has been offered for the spring-summer evolution of the NAM system.

54 ENVIRONMENTAL SCIENCES↗

Investigating Low-Altitude Constellations of Ad-Hoc Lunar PNT System for Distributed Spacecraft Autonomy

In this study, we examine a low-altitude Lunar Position, Navigation, and Timing (LPNT) constellations and the localization performance of Centralized Extended Kalman Filter (CEKF) and Decentralized Extended Kalman Filter (DEKF) algorithms. The primary investigation involves a 100-node swarm operating at a 100 km altitude, in contrast to previous studies that examined a 21-node asset in a frozen-orbit at 5,500 km. The autonomous operation of large-scale swarm is based on two-way Inter-Satellite Link (ISL) measurements, which involve pseudoranges and relative velocities among swarm nodes. We perform a numerical assessment of the two filtering approaches, utilizing ‘fully sampled’ measurements from all available assets as well as ‘two ISL’ measurements where each spacecraft is restricted to only two antennas. This research includes an analysis of CEKF under 2-ISL constraints and evaluates the performance of DEKF in a 100-node swarm, which has not been explored in previous studies. In addition, we examine the impact of increasing the sampling frequency for DEKF, showing that the update cycle can be shortened from a 10-minute interval. A novel approach for ‘2-ISL limited’ DEKF will also be introduced, using a matching formulation that exhaustively enumerates all potential matches. This study provides valuable insights into large-scale distributed swarm operations, considering various filter configurations, sampling frequencies, matching strategies, and scalability of CEKF and DEKF for low-altitude LPNT applications. The Lunar PNT technology plays a key role in providing reliable and robust navigation services on the Moon's surface and the South pole, where the primary Lunar missions are planned. To support upcoming Lunar missions, including small satellites from NASA's Commercial Lunar Payload Services program, the Lunar PNT system must be adaptable to smaller platforms like CubeSats. Driven by the growing involvement of public and private exploration partnerships, the traditional low Earth orbit missions are shifting to beyond geosynchronous orbit [1]. These upcoming missions aim to foster a sustainable and innovative exploration program, in collaboration with commercial and international partners, to facilitate human expansion throughout the solar system and return new knowledge and opportunities to Earth [2]. As part of this trend, there are increasing efforts to utilize science missions in Lunar orbit to develop a non-dedicated and ad-hoc PNT network system. Two traditional approaches, the Deep Space Network (DSN) and the weak signal Global Positioning System (GPS), are established deep-space navigation technologies for missions beyond the geosynchronous orbit. Beginning in 1958, the DSN was developed to communicate with the Explorer 1 spacecraft based on the use of radiometric tracking in spacecraft navigation [3]. The DSN is capable of providing nearly unfettered coverage to spacecraft beyond low-Earth orbit (LEO), however, increased space mission volume has created concerns about future expectations of DSN usage for spacecraft navigation [4]. For cislunar mission applications, the position accuracy using DSN achieves 100 m (3σ) with at least three geometrically diverse ground stations when using radiometric tracking alone [5]. The DSN's dependence on Earth-based ground stations restricts its operational capabilities to periods of Earth visibility. This limitation, coupled with its poor localization performance, renders the DSN unsuitable for future lunar missions that demand continuous tracking and precise positioning. To satisfy the increasing requirements of DSN in Lunar applications, spacecrafts are also required to improve their onboard antenna power and efficiency of the transmission. However, there is an important aggregate cost trade between adding capabilities to every spacecraft and adding to a capacity on the ground that serves multiple spacecraft [6]. A weak GPS system can provide PNT service while the user spacecraft is bound to the Moon, leveraging a single, steerable high gain antenna with the relatively narrow beam which includes all the sources in its field of view [7]. However, the higher the altitude the receiver is above the GPS constellations, the poorer and the weaker are the relative geometry and the received signal powers, respectively, leading to a significant navigation accuracy reduction [8]. The transmitted power becomes weaker with increasing distance from the Earth as well as signals tracked from one of the side lobes of the GPS antenna pattern. As a results, the number of visible satellites and relative geometric condition of the GPS satellites at very high altitude drops dramatically and reduces the navigation solution accuracy. Therefore, the weak GPS system is also not an ideal way to provide PNT service to upcoming Lunar missions when considering its limited geometric condition and the recued navigation accuracy. Another navigation approach on the Moon is being developed, similar to the Global Navigation Satellite System (GNSS) on Earth, aiming to offer navigation service with continuous 24/7 coverage across the entire Lunar surface. For example, lunar communications relay and navigation systems (LCRNS) by NASA and Lunar navigation satellite systems (LNSS) by JAXA are designed to serve as dedicated Position, Navigation, and Timing (PNT) systems for the Moon. However, designing a dedicated LNSS and PNT service involves additional challenges, which are unique to the lunar environment, including limited payload capacity for the CubeSat platform, i.e., the size, weight, and power (SWaP) of the onboard clock, limited lunar ground monitoring stations, and limited financial investment as compared to the legacy Earth-GPS [9]. NASA’s focus on utilizing CubeSat platforms on the Moon leads to an alternative Lunar navigation platform that leverages the existing Lunar science and exploration assets. The small satellites used in Lunar missions can be used to create a low-cost, autonomous, ad-hoc, and on-demand mission-centric Lunar PNT swarm capable of providing PNT services to these low-cost lunar missions [10]. As upcoming Lunar missions will often operate at low-altitude about 30 km to 100 km for scientific observations and mapping purposes, the low-altitude orbital constellations could be employed to create an ad-hoc Lunar PNT system. However, several issues must be addressed, such as the instability of these orbits, which often require maintenance or are only suitable for short-duration missions, operating for fewer than 90 days. Additionally, at an altitude of 100 km, the satellites have a limited period during which they are above the horizon and capable of providing PNT service to users. The implementation of a non-dedicated, ad-hoc Lunar navigation constellation facilitates on-demand PNT services. A preliminary study of ad-hoc Lunar PNT system was conducted using 21 spacecraft in 5,5000 km altitude frozen orbits to test its feasibility and a basic performance of orbital asset localization among ad-hoc Lunar constellations in small satellites format [10]. These swarm assets are designed for autonomous localization with minimal Earth interaction, reducing dependency on bandwidth and ground resources. The design in [10] demonstrated the feasibility of a decentralized PNT approach, specifically employing a DEKF approach for state estimation, which helps minimize onboard operating costs. The DEKF method distributes computation across individual satellites, which lightens the computational load while maintaining accuracy in orbit ephemeris and clock offsets, similar to centralized systems [11]. In a follow-on study [12], each spacecraft was limited to 2 communications antennae, forcing the selection of measurements and scheduling spacecraft activities to perform the measurements. A matching algorithm is implemented to select the best measurements and schedule position estimation updates. The decentralized localization performance is also investigated with increasing levels of network degradation for swarm assets considering the impact of intermittent and permanent communication failure, to demonstrate the robustness and fidelity of the decentralized Lunar PNT service [13]. This study confirmed that the ad-hoc PNT constellations in frozen orbit are highly robust and resilient to communication failures. However, unlike frozen orbit swarm assets, the low-altitude satellites have a limited ground view at an altitude of 100 km, where the ad-hoc Lunar constellation consists of 98 low-altitude satellites, evenly distributed across seven circular polar orbital planes, alongside two satellites in a frozen orbit at an altitude of 5,500 km (Figure 1). Therefore, the number of satellites visible to ground users is significantly limited in low-altitude orbit constellations. As each visibility of a spacecraft remains intact for only a few ticks before it moves out of the field of view, the ground user encounters challenges in maintaining continuous navigation service, resulting in sparse availability and provision of Lunar PNT system. Consequently, service availability is primarily restricted to the Lunar South Pole region (Figure 2). Given these limitations and concerns, the localization performance of low-altitude swarm assets will be assessed in this study. We focus on the investigation of the localization performance of low-altitude swarm assets and ground users near the Lunar South Pole. The overall flow of the Lunar PNT simulation incorporates the DEKF approach of asset localization and the weighted least-squares approach in user localization (Figure 3). The autonomous Lunar PNT simulation is primarily implemented in MATLAB, where the DEKF based on the matching scheduler is implemented with Google’s OR-tools as a model builder and Gurobi optimization tool as a backend solver. The General Mission Analysis Tool (GMAT) is utilized to generate ephemeris data for swarm assets, and accounts for satellite orbital details, mass, and perturbations like solar radiation pressure and drag coefficients. Each ephemeris dataset is produced in the Moon International Celestial Reference Frame (ICRF) inertial coordinate system. For state estimation, the distributed swarm assets rely on two-way Inter-Satellite Link (ISL) measurements, which involve tracking pseudoranges and relative velocities between visible satellites and anchor nodes during each observation. Numerical evaluations of the decentralized localization process are conducted to demonstrate the feasibility of the low-altitude PNT system in providing reliable navigation services. The main approach involves using DEKF and CEKF to localize 100 satellites in low-altitude constellations, where the CEKF is implemented to serve as a baseline for comparing the performance of distributed algorithms. In both cases, we evaluate ‘fully sampled’ measurements from all available assets, and ‘two ISL’ measurements when spacecraft are constrained to have only two antennas. We test four estimation techniques: CEKF fully sampled, CEKF two ISL, DEKF fully sampled, and DEKF two ISL filters. As the DEKF update cycle is comprised of network setup, communication, and computations, a global broadcast network and 2-way ISL network setup will take from 4 to 6 minutes as maximum [12]. In this simulation, the DEKF update cycle is set to 10 minutes, including a 4-minute latency for obtaining and computing the actual measurement updates. We experiment an increased update cycle to demonstrate the feasibility and evaluate the impact on localization performance using various tuning values for measurement noise covariances (Figures 4 and 5). By comparing centralized and decentralized approaches using a matching algorithm, we analyze the influence of cross-correlation factors in the covariance matrix, assuming 100% reliability of all assets and measurements. The increased frequency and the adjustments of tuning parameters reveal distinct error patterns between the two scenarios. The localization accuracy of the swarm assets and ground users is assessed by taking the median error across 100 assets and one ground user (84.9°S, 137.5°E) over 7-day simulation period (Table 1). Since the user localization accuracy is significantly affected by the performance of the swarm assets, it is crucial to maintain high localization accuracy within the swarm. This study will continue to explore decentralized filtering for autonomous LPNT operations, with further investigation of an 'iterative' matching approach which enumerates every valid matching pair, planned for the following month.

Yeji Kim↗

Tuning Organic Semiconductor Packing and Morphology through Non-equilibrium Solution Processing (Final Report)

Organic semiconductors (OSCs) are a promising candidate to produce low-cost, large area, and flexible electronics. There has previously been successful development of organic field effect transistors (OFETs), photovoltaics (OPV), and bioelectronics using OSCs. Solution processing of these materials allows for the fabrication of large-area devices in the kinetic crystallization regime. The charge transport capabilities have been shown to depend on the morphology and molecular packing of the OSC within thin films. Our hypothesis is that solution processing conditions may significantly impact OSC morphology and as a result the charge transport. Therefore, this work has focused on obtaining a fundamental understanding of the morphology and molecular packing of OSCs for optimal device performance. Through our work, we have gained a better understanding of how these structural conditions were influenced by the solution processing conditions. Our studies have led to a more systematic understanding of the various parameters that impact OSC morphology. Solution processing leads to non-equilibrium films, thus allowing for the formation of diverse morphologies that are inaccessible by other fabrication methods. Previously, our group has focused on tuning the morphology of small-molecular OSCs [53, 55-57]. However, little work had been devoted to polymer OSCs. There was a lack of detailed studies characterizing the solution-state of polymer OSCs in terms of their conformation, degree of entanglement, polymer aggregation, and chain relaxation dynamics. The solution state properties are also likely to be influenced by the rigidity and molecular weight of the polymer OSCs investigated. Thus, we have focused our attention on gaining a better insight of polymer OSC films and solution processing methods. Through this proposal, our approach is to investigate the correlation between solution-state properties and the morphology of the resulting polymer OSC films. We worked three specific aims: investigate the effects of 1) polymer OSC solution-state properties on final film morphology; 2) molecular additives on the solution-state and final film properties; 3) controlled pre-aggregation in the solution-state in the final film morphology. More rigid and planar polymer backbones should promote interchain charge transport and more efficient interchain hopping between polymer chains. Through tailoring the polymer backbone, polymer sidechains, and molecular weight, we expected the altered solution-state properties to affect the final film morphology. In addition, reducing the entanglements and promoting chain alignment will likely prevent charge carrier trapping through conformation disorders. We thus studied different mechanisms, such as molecular additives, to reduce entanglements in solution. Devices fabricated from these solutions were expected to have improved charge transport abilities. In addition to tailoring the solution-state characteristics of the polymer OSCs, we investigated the effect of solution processing methods on the molecular packing and morphology of polymer OSC films. Two main solution processing techniques, e.g., spin-coating and solution shearing, were employed to fabricate OFETs. Spin coating was employed to prepare OSC films, which produce isotropic films. This technique creates several parameters to tune such as spin coating speed, acceleration, and time. Additionally, our research group developed the solution shearing method, which consists of the solution initially sandwiched between two plates. By sliding the top plate, the solution front is exposed and drying begins. This technique allows for the creation of aligned large crystalline domains. Solution shearing consists of different processing parameters which can affect the final morphology, such as shearing speed, substrate temperature and temperature gradient, distance between both plates, and the tilt angle of the top plate. Due to the complexity of the polymer systems investigated, numerous techniques were employed to characterize the polymer OSC solution state and films. All the materials were characterized using the DOE supported synchrotron X-ray scattering facilities at the Stanford Synchrotron Radiation Lightsource (SSRL). Using grazing incidence X-ray scattering (GIXS) and Near Edge X-ray Absorption Fine Structure (NEXAFS) techniques, the crystalline structure and molecular orientations of the thin films were measured. Optical absorption (UV-Vis) spectroscopy was employed to determine the aggregation state in solution and films of the polymer systems. Polarized UV-Vis also allowed for the determination of the relative degree of polymer chain alignment for solution sheared films. Additionally, various other techniques were used to investigate other properties within the film, such as atomic force microscopy (AFM) and solution rheology. Finally, the device performance is quantified through the fabrication and characterization of OFETs, which will highlight the effects of morphology on charge transport.

36 MATERIALS SCIENCE↗

Principles of water oxidation and O2-based hydrocarbon transformation by multinuclear catalytic sites

Abstract The central thrust of this integrated experimental and computational research program was to obtain an atomistic-level understanding of the structural and dynamic factors underlying the design of catalysts for water oxidation and selective reductant-free O2-based transformations. The focus was on oxidatively robust polyoxometalate (POM) complexes in which a catalytic active site interacts with proximal metal centers in a synergistic manner. Thirty five publications in high-impact journals arose from this grant. I. Developing an oxidatively and hydrolytically stable and fast water oxidation catalyst (WOC), a central need in the production of green fuels using water as a reductant, has proven particularly challenging. During this grant period we have designed and investigated several carbon-free, molecular (homogenous), oxidatively and hydrolytically stable WOCs, including the Rb8K2[{Ru4O4(OH)2(H2O)4}(γ-SiW10O36)2]·25H2O (1) and [Co4(H2O)2(α-PW9O34)2]10- (2). Although complex 1 is fast, oxidatively and hydrolytically stable WOC, Ru is neither abundant nor inexpensive. Therefore, development of a stable and fast carbon-free homogenous WOC, based on earth-abundant elements became our highest priority. In 2010, we reported the first such catalyst, complex 2. This complex is substantially faster than 1 and stable under homogeneous conditions. Recently, we have extended our efforts and reported a V2-analog of the complex 2, i.e. [Co4(H2O)2(α-VW9O34)2]10- (3), which shows an even greater stability and reactivity. We succeeded in: (a) immobilizing catalysts 1 and 2 on the surface of various electrodes, and (b) elucidating the mechanism of O2 formation and release from complex 1, as well as the Mn4O4L6 “cubane” cluster. We have shown that the direct O-O bond formation is the most likely pathway for O2 formation during water oxidation catalyzed by 1. II. Oxo transfer catalysts that contain two proximal and synergistically interacting redox active metal centers in the active site form another part of considerable interest of our grant because species with such sites [including methane monooxygenase (MMO) and more] are some of the most effective oxygenase catalysts known. Our team conducted the following research on γ-M2-Keggin complexes: (a) investigated stability of the trimer [{Fe3(OH)3(H2O)2}3(γ-SiW10O36)3]15-, 4, in water, and developed the chemistry and catalysis of the di-iron centered POM, [γ(1,2)-SiW10{Fe(OH)}2O38]6-, 5, in organic solvents (Figure 2). We also study the thermodynamic and structural stability of γ-M2-Keggin in aqueous media for different M’s (d-electron metals). We have defined two structural classes of POMs with proximally bound d-electron metal centers. We refer to these structural isomers of the {γ-M2SiW10} family of POMs as “in-pocket” and “out-of pocket”. We have elucidated the factors controlling the structure and stability of the V, Fe, Ru, Tc, Mo and Rh derivatives of [(SiO4)M2(OH)2W10O32]4- using a range of computational tools. We have: (a) demonstrated that heteroatom X in these polyanions may function as an “internal switch” for defining the ground electronic states and, consequently, the reactivity of the γ-M2-Keggin POM complexes; (b) elucidated reactivity of divacant lacunary species and polyperoxotungstates (PPTs), {Xn+O4[WO(O2)2]4}n-, which could be degradation products of γ-M2-Keggin complexes in aqueous media; (c) elucidated the role of the POM ligand in stabilization of {Ru2} and {(Ru-oxo)2} fragments in the reactant and product of the reaction of {γ-[(Xn+O4)Ru2(OH)2W10O32]}(8-n)- (where X = Si4+, P5+ and S6+) with O2, and (d) the mechanisms of olefin epoxidation catalyzed by these di-d-transition metal substituted and divacant lacunary γ-M2-Keggin complexes. III. Complementing the efforts presented above was the development of less time-consuming but reasonably accurate computational methods allowing one to explore more deeply large catalytic systems. We developed Reactive Force Field (ReaxFF) to study interaction of the targeted POMs with water, proton and hydroxide ions in the liquid phase. We tested our ReaxFF parameters on the Lindqvist POMs, M6O19n-, where M = Nb and Ta. These parameters are made available as part of the ReaxFF code. In addition, we have developed parameters for Sc, Ti, Fe, Co and Ni in combination with H, C, N, O, as well as the same metal (M-M) for the spin-polarized self-consistent-charge density-functional tight-binding (DFTB) method. Test calculations showed that the DFTB method with the present parameters in most cases reproduces structural properties very well. These parameters are made available as part of the DFTB code. Thus, this DOE BES funded research project has clarified several key areas impacting (a) water oxidation and O2-based hydrocarbon transformation, (b) stabilization of key structures and catalytic intermediates in such processes, (c) immobilization of molecular catalysts on metal oxide surfaces, and (d) application of optimal computational methods to study reaction dynamics in large systems.

58 GEOSCIENCES↗

Multi-resolution Arctic Shrub Cover Dataset Derived from UAS and Airborne SfM and LiDAR (2013-2025)

We synthesized 177 unoccupied aerial system flights and 77 airborne flights across the Arctic and created a multi-resolution benchmark data of low-to-tall shrub fractional cover leveraging Structure-from-Motion and Light Detection and Ranging. The resulting dataset covered a total of 1899 km2 across Alaska, Western Canada, Sweden, and Siberian Arctic, including key sites from the Oro Arctic to the High Arctic. The dataset is organized into 6 primary data collection directories (“Abisko,” “AWI,” “ERE,” “Fairbanks,” “NGEE,” “Toolik”), each containing site and flight subdirectories. Flight directories include shrub cover rasters (*.tifs) at 1 m, 5 m, and 30 m resolution, the canopy height model at 1 m resolution (*.tifs), and a bounding box *.kml file. For the AWI, Abisko, NGEE, and Fairbanks collections, we also include the GCC raster at 1 m resolution (*.tif). Files are organized by Collection > Site > Flight Name > Data Files. Flight rasters are in the local UTM zone and the .kml files are in the geographic coordinate system EPSG 4326. We also include a .csv file that details the source datasets for every flight. The Next-Generation Ecosystem Experiments in the Arctic (NGEE Arctic) project is a research effort to reduce uncertainty in the Department of Energy’s Energy Exascale Earth System Model (E3SM) by developing a predictive understanding of Arctic tundra ecosystems underlain by permafrost and to quantify feedbacks from the Arctic tundra to the Earth system. NGEE Arctic is supported by the Department of Energy's Office of Biological and Environmental Research. Over Phases 1–3, observations made by the NGEE Arctic team across a gradient of permafrost landscapes in Arctic Alaska improved the representation of tundra processes in the land surface component of E3SM (the E3SM Land Model, ELM). Model improvements emphasized unique aspects of permafrost environments and explored reductions in model complexity while retaining predictive power. The Arctic-informed ELM developed by NGEE Arctic has been used to make novel predictions on processes ranging from permafrost thaw to soil biogeochemical cycling to Earth system feedbacks associated with the unique characteristics of tundra plants. In Phase 4, the NGEE Arctic team is evaluating our new predictive understanding under novel conditions across the Arctic domain. In collaboration with partners at long-term pan-Arctic research sites we are examining whether an Arctic-informed ELM can faithfully simulate interactions among surface and subsurface processes at site, regional, and pan-Arctic scales. In turn, we are using variety of tools to dynamically extend and evaluate ELM inference, with an emphasis on data synthesis and pan-Arctic model evaluation, reintegration of code with an evolving E3SM, scaling across heterogeneous Arctic landscapes, and the appropriate representation of the impacts of increasingly frequent Arctic disturbances.

canopy height model↗