Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Data depth”

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 469 records · Page 26

Millimeter-wave observations of Euclid Deep Field South using the South Pole Telescope: A data release of temperature maps and catalogs

Context. The South Pole Telescope third-generation camera (SPT-3G) has observed over 10,000 square degrees of sky at 95, 150, and 220 GHz (3.3, 2.0, 1.4 mm, respectively) and will significantly overlap the ongoing 14,000 square-degree Euclid Wide Survey. The Euclid collaboration recently released Euclid Deep Field South (EDF-S) observations of 23 square degrees at wide field depths in the first quick data release (Q1). Aims. With the goal of releasing complementary millimeter-wave data and encouraging legacy science, we performed dedicated observations of a 57-square-degree field overlapping the EDF-S. Methods. The observing time totaled 20 days, and we reached noise depths of 4.3, 3.8, and 13.2 $μ$K-arcmin at 95, 150, and 220 GHz, respectively. Results. In this work we present the temperature maps and two catalogs constructed from these data. The emissive source catalog contains 601 objects (334 inside EDF-S) with 54% synchrotron-dominated sources and 46% thermal dust emission-dominated sources. The 5$σ$ detection thresholds are 1.7, 2.0, and 6.5 mJy in the three bands. The cluster catalog contains 217 cluster candidates (121 inside EDF-S) with median mass $M_{500c}=2.12 \times 10^{14} M_{\odot}/h_{70}$ and median redshift $z$ = 0.70, corresponding to an order-of-magnitude improvement in cluster density over previous tSZ-selected catalogs in this region (3.81 clusters per square degree). Conclusions. The overlap between SPT and Euclid data will enable a range of multiwavelength studies of the aforementioned source populations. This work serves as the first step toward joint projects between SPT and Euclid and provides a rich dataset containing information on galaxies, clusters, and their environments.

Archipley, M. [Chicago U., Astron. Astrophys. Ctr.↗

A Methodology for Surface Soil Moisture and Vegetation Optical Depth Retrieval Using the Microwave Polarization Difference Index

A methodology for retrieving surface soil moisture and vegetation optical depth from satellite microwave radiometer data is presented. The procedure is tested with historical 6.6 GHz brightness temperature observations from the Scanning Multichannel Microwave Radiometer over several test sites in Illinois. Results using only nighttime data are presented at this time, due to the greater stability of nighttime surface temperature estimation. The methodology uses a radiative transfer model to solve for surface soil moisture and vegetation optical depth simultaneously using a non-linear iterative optimization procedure. It assumes known constant values for the scattering albedo and roughness. Surface temperature is derived by a procedure using high frequency vertically polarized brightness temperatures. The methodology does not require any field observations of soil moisture or canopy biophysical properties for calibration purposes and is totally independent of wavelength. Results compare well with field observations of soil moisture and satellite-derived vegetation index data from optical sensors.

Owe, Manfred↗

Role of Forcing Uncertainty and Background Model Error Characterization in Snow Data Assimilation

Accurate specification of the model error covariances in data assimilation systems is a challenging issue. Ensemble land data assimilation methods rely on stochastic perturbations of input forcing and model prognostic fields for developing representations of input model error covariances. This article examines the limitations of using a single forcing dataset for specifying forcing uncertainty inputs for assimilating snow depth retrievals. Using an idealized data assimilation experiment, the article demonstrates that the use of hybrid forcing input strategies (either through the use of an ensemble of forcing products or through the added use of the forcing climatology) provide a better characterization of the background model error, which leads to improved data assimilation results, especially during the snow accumulation and melt-time periods. The use of hybrid forcing ensembles is then employed for assimilating snow depth retrievals from the AMSR2 (Advanced Microwave Scanning Radiometer 2) instrument over two domains in the continental USA with different snow evolution characteristics. Over a region near the Great Lakes, where the snow evolution tends to be ephemeral, the use of hybrid forcing ensembles provides significant improvements relative to the use of a single forcing dataset. Over the Colorado headwaters characterized by large snow accumulation, the impact of using the forcing ensemble is less prominent and is largely limited to the snow transition time periods. The results of the article demonstrate that improving the background model error through the use of a forcing ensemble enables the assimilation system to better incorporate the observational information.

assimilation↗

Occurrences of REE and CM in volcanic ash beds associated with coals in SW Wyoming & NW Colorado

Rare earth elements (REEs) and critical minerals (CMs) are vital components in the nation’s economy, national security, and the further strengthening of our renewable energy capabilities and infrastructure. Recent work has revealed more insight into previously unknown modes of occurrence for REEs, particularly in coal seams and adjacent noncoal sediments, but many data gaps persist. Temporal relationships between coal beds and adjacent tonstiens, ash beds, or other partings are widely recognized but not fully understood. Variations in REE concentration and their direct or indirect association with coal beds are not predictable. This study presents findings from examining lateral and potentially stratigraphic variation of REE and CM occurrences in coals of late Cretaceous and Paleocene age. In this study we collected and examined samples of differing ages both core and coal mine exposures for qualitative x-ray fluorescence (XRF), and multi-element geochemistry. Preliminary XRF data indicates variable levels of light REEs in a core from the Yampa coalfield of Moffat County in northwest Colorado. An absence of similar enrichments in slightly older coals in the Adaville Formation of the Fold and Thrust belt in southwest Wyoming is described from XRF and ICP-MS evaluation. Previous work has shown an association with the Yampa Bed and REE enrichment in coals of the Upper Cretaceous Williams Fork Formation. While not common and widespread, this association suggests that leaching or reworking of the ash bed could lead to elevated REEs in the adjacent coal. We investigate differences between the two data occurringas geologic compositions and post-depositional histories. In-depth petrography and precise geochemistry should reveal data that at the very least will help improve some aspects of the geologic models of coal basins and the effect of an influx of non-coal sediments.

REE, CM, ash, volcanic ash, coal↗

NASA Tech Briefs, March 2013

Topics covered include: Remote Data Access with IDL Data Compression Algorithm Architecture for Large Depth-of-Field Particle Image Velocimeters Vectorized Rebinning Algorithm for Fast Data Down-Sampling Display Provides Pilots with Real-Time Sonic-Boom Information Onboard Algorithms for Data Prioritization and Summarization of Aerial Imagery Monitoring and Acquisition Real-time System (MARS) Analog Signal Correlating Using an Analog-Based Signal Conditioning Front End Micro-Textured Black Silicon Wick for Silicon Heat Pipe Array Robust Multivariable Optimization and Performance Simulation for ASIC Design; Castable Amorphous Metal Mirrors and Mirror Assemblies; Sandwich Core Heat-Pipe Radiator for Power and Propulsion Systems; Apparatus for Pumping a Fluid; Cobra Fiber-Optic Positioner Upgrade; Improved Wide Operating Temperature Range of Li-Ion Cells; Non-Toxic, Non-Flammable, -80 C Phase Change Materials; Soft-Bake Purification of SWCNTs Produced by Pulsed Laser Vaporization; Improved Cell Culture Method for Growing Contracting Skeletal Muscle Models; Hand-Based Biometric Analysis; The Next Generation of Cold Immersion Dry Suit Design Evolution for Hypothermia Prevention; Integrated Lunar Information Architecture for Decision Support Version 3.0 (ILIADS 3.0); Relay Forward-Link File Management Services (MaROS Phase 2); Two Mechanisms to Avoid Control Conflicts Resulting from Uncoordinated Intent; XTCE GOVSAT Tool Suite 1.0; Determining Temperature Differential to Prevent Hardware Cross-Contamination in a Vacuum Chamber; SequenceL: Automated Parallel Algorithms Derived from CSP-NT Computational Laws; Remote Data Exploration with the Interactive Data Language (IDL); Mixture-Tuned, Clutter Matched Filter for Remote Detection of Subpixel Spectral Signals; Partitioned-Interval Quantum Optical Communications Receiver; and Practical UAV Optical Sensor Bench with Minimal Adjustability.

Source record↗

Planck 2018 results: III. High frequency instrument data processing and frequency maps

This paper presents the High Frequency Instrument (HFI) data processing procedures for the Planck 2018 release. Major improvements in mapmaking have been achieved since the previous Planck 2015 release, many of which were used and described already in an intermediate paper dedicated to the Planck polarized data at low multipoles. These improvements enabled the first significant measurement of the reionization optical depth parameter using Planck-HFI data. This paper presents an extensive analysis of systematic effects, including the use of end-to-end simulations to facilitate their removal and characterize the residuals. The polarized data, which presented a number of known problems in the 2015 Planck release, are very significantly improved, especially the leakage from intensity to polarization. Calibration, based on the cosmic microwave background (CMB) dipole, is now extremely accurate and in the frequency range 100-353 GHz reduces intensity-to-polarization leakage caused by calibration mismatch. The Solar dipole direction has been determined in the three lowest HFI frequency channels to within one arc minute, and its amplitude has an absolute uncertainty smaller than 0.35 μK, an accuracy of order 10 -4 . This is a major legacy from the Planck HFI for future CMB experiments. The removal of bandpass leakage has been improved for the main high-frequency foregrounds by extracting the bandpass-mismatch coefficients for each detector as part of the mapmaking process; these values in turn improve the intensity maps. This is a major change in the philosophy of "frequency maps", which are now computed from single detector data, all adjusted to the same average bandpass response for the main foregrounds. End-to-end simulations have been shown to reproduce very well the relative gain calibration of detectors, as well as drifts within a frequency induced by the residuals of the main systematic effect (analogue-to-digital convertor non-linearity residuals). Using these simulations, we have been able to measure and correct the small frequency calibration bias induced by this systematic effect at the 10 -4 level. There is no detectable sign of a residual calibration bias between the first and second acoustic peaks in the CMB channels, at the 10 -3 level.

79 ASTRONOMY AND ASTROPHYSICS↗

Oregon Rivers Characterization Tool

DHS’s Regional Resiliency Assessment Program (RRAP) undertook the Oregon State Transportation Systems RRAP project between 2018-2021 to assess the seismic vulnerabilities of the state’s transportation system to a Cascadia Subduction Zone (CSZ) earthquake (CISA 2021). The project primarily focused on the effects of a CSZ earthquake on the mobility of emergency response efforts into the region. In the event that bridges located at crucial river crossings fail during a CSZ earthquake, rivers will become natural barriers affecting the ability to support response and recovery efforts. The characteristics of rivers (e.g., width, depth, and flow) are critical data to inform planning for the recovery of transportation systems across river barriers and determine reopening times following a CSZ earthquake. As part of the transportation system assessment, an Argonne National Laboratory (Argonne) project team developed a River Characterization Tool (RCT) to quantify the physical characteristics of the river system in Oregon, which can be used to inform a system-level analysis of transportation recovery and reopening times. This report mainly discusses approaches, methods, development, and applications of the RCT. The use of the river data generated from this tool for a system-level analysis of transportation recovery is discussed in the Oregon Transportation Systems RRAP Resiliency Assessment report (CISA 2021). The primary goal of the RCT is to generate statewide dataset of widths, depths, and flows for all rivers across Oregon to assess the effects of river barriers on emergency response and recovery efforts following a CSZ earthquake. The RCT leverages remote sensing data that multispectral sensors collect from satellites, employs a deep learning model to characterize river predictions, and computes river channel geometric dimensions with hydraulic relationships and river predictions. The entire process of river characterization is automated and implemented through the RCT.

58 GEOSCIENCES↗

Crustal interpretation of the MAGSAT data in the continental United States

The processing of MAGSAT scalar data to construct a crustal magnetic anomaly map over the continental U.S. involves removal of the reference field model, a path-by-path subtraction of a low order polynomial through a least-squares fit to reduce orbital offset errors, and a two dimensional spectral filtering to mitigate the spectral bias induced by the path-by-path orbital correction scheme. The resultant anomaly map shows reasonably good correlations with an aeromagnetic map derived from the project MAGNET. Prominent satellite magnetic anomalies are identified in terms of geological provinces and age boundaries. An inversion method was applied to MAGSAT data which produces both the Curie depth topography and laterally varying magnetic susceptibility of the crust. A contoured Curie depth map thus derived shows general agreements with a crustal thickness map based on seismic data.

Won, I. J.↗

Voyager 2 photopolarimeter observations of Titan

Observations of Titan's whole-disk polarization at 2460 and 7500 A are presented and analyzed in terms of model scattering atmospheres. If the Titan aerosols are spherical or nearly spherical, no single combination of refractive index and size distribution is able to fit data at both wavelengths. However, a vertically inhomogeneous distribution suggested by Tomasko and Smith (1980), characterized by a size gradient with altitude, fits the data at 2640 A moderately well but must be modified at intermediate and large optical depths to fit the 7500-A data. Results for synthetic phase functions indicate that the single-scattering polarization must be 70 percent or larger in the UV and 78 percent or larger in the near-IR at 90-deg phase angle, depending on the phase function. If the correct phase function is similar to that for 0.5-micron-radius spheres, the UV single-scattered polarization must be 84 percent and the near-IR single-scattered polarization must be over 90 percent. Such large polarizations are impossible for 0.5-micron-radius spheres but may be possible for nonspherical particles with effective radii near 0.5 micron, although the existence of nonspherical particles with the scattering properties required by these and other observations has not been demonstrated.

West, R. A.↗

Methodology to Define Delivery Accuracy Under Current Day ATC Operations

In order to enable arrival management concepts and solutions in a NextGen environment, ground- based sequencing and scheduling functions have been developed to support metering operations in the National Airspace System. These sequencing and scheduling algorithms as well as tools are designed to aid air traffic controllers in developing an overall arrival strategy. The ground systems being developed will support the management of aircraft to their Scheduled Times of Arrival (STAs) at flow-constrained meter points. This paper presents a methodology for determining the undelayed delivery accuracy for current day air traffic control operations. This new method analyzes the undelayed delivery accuracy at meter points in order to understand changes of desired flow rates as well as enabling definition of metrics that will allow near-future ground automation tools to successfully achieve desired separation at the meter points. This enables aircraft to meet their STAs while performing high precision arrivals. The research presents a possible implementation that would allow delivery performance of current tools to be estimated and delivery accuracy requirements for future tools to be defined, which allows analysis of Estimated Time of Arrival (ETA) accuracy for Time-Based Flow Management (TBFM) and the FAA's Traffic Management Advisor (TMA). TMA is a deployed system that generates scheduled time-of-arrival constraints for en- route air traffic controllers in the US. This new method of automated analysis provides a repeatable evaluation of the delay metrics for current day traffic, new releases of TMA, implementation of different tools, and across different airspace environments. This method utilizes a wide set of data from the Operational TMA-TBFM Repository (OTTR) system, which processes raw data collected by the FAA from operational TMA systems at all ARTCCs in the nation. The OTTR system generates daily reports concerning ATC status, intent and actions. Due to its availability, ease of use, and vast collection of data across several airspaces it was determined that the OTTR data set would be the best method to utilize moving forward with this analysis. The particular variables needed for further analysis were determined along with the necessary OTTR reports, by working closely with the repository team additional analysis reports were developed that provided key ETA and STA information at the freeze horizon. One major benefit of the OTTR data is that using the correct reports the data across several airports could be analyzed over large periods of time. The OTTR data processes the TBFM data daily and is stored in various formats across several airspaces. This allowed us to develop our own parsing methods and raw data processing that would not rely on other computationally expensive tools that perform more in depth analysis of similar sets of data. The majority of this work consisted of the development of the ability to filter flights to create a subset of flights that could be considered undelayed, which is defined as a flight at the freeze horizon with an ETA and STA difference that was minimal or close to zero. This was a broad method that allowed the consideration of a large data set which consisted of all the traffic across a two month period in 2013, the hottest and coldest months, arriving into four airports: George Bush Intercontinental, Denver International, Los Angeles International, and Phoenix Sky Harbor.

delivery accuracy↗

Challenges and Lessons Learned from an Analysis of Three Zero Energy Buildings

Zero energy buildings, or zero energy ready buildings, which are designed and operated by public or private commercial property owners can play an important role in reducing carbon emissions. This paper discusses lessons learned and key takeaways from an in-depth analysis of three zero energy buildings that took an integrated design and construction approach to significantly reduce energy use. Two of these projects are new construction and one is a retrofit to zero energy. Findings are based on project literature review, data analysis, and in-depth interviews with the building design teams and staff who were involved during the design, construction, and operation of the buildings. The paper addresses value proposition and cost data in such a way that other building owners can replicate the strategies and technology solutions in response to regulatory mandates or organizational goals. It also provides details on operational improvements taken at each zero energy building to enhance energy performance and increase the potential for load flexibility, and discusses challenges and lessons gleaned from design teams and building staff. The findings serve as a reference for building owners, designers, engineers, contractors, or others interested in, or involved with, the design, construction, or use of new or existing buildings. The paper also includes recommended pathways for widespread adoption of zero energy strategies that can be applied in various locations.

ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATION↗

Revised depolarization corrections for atmospheric extinction

Possible sources of errors in depolarization data are analyzed with special reference to the rotational Raman contribution to the optical depth. It is shown that the depolarization data now available are so accurate that their uncertainty affects the calculated Rayleigh optical depth by less than 1%, and consequently, the anomalously low extinctions occasionally reported must be due to systematic errors of measurement rather than errors in the Rayleigh optical depths.

Young, A. T.↗

GIS-Based Modeling of Contaminated Soil Volumes at Multiple Sites in the Formerly Utilized Sites Remedial Action Program - 20149

The remediation of hazardous, toxic, and radioactive waste (HTRW) sites produces cost-related risks associated with the estimation of contaminated soil or debris volumes. Historical risk-management techniques include cost contingencies to cover volume uncertainties that affect project budgeting and decision-making. The Buffalo District teamed with project partners to lessen volume uncertainty and reduce project risks at multiple HTRW sites managed under the Formerly Utilized Sites Remedial Action Program (FUSRAP). Historical remedial investigations under FUSRAP commonly identified the presence of radiological material in site media, the associated human health risk, and then areas of remediation. To manage remedial execution and reduce risk, pre-design or remediation-phase sampling essentially 'chased' contamination, which was not conducive to efficient predictive budgeting derived from Feasibility Study (FS) cost analyses. The Buffalo District first optimized their approach to better understand volume uncertainty by utilizing the Argonne National Laboratory's Bayesian Approaches for Adaptive Spatial Sampling (BAASS) software [1]. BAASS processed soft data (e.g., gamma walk-over data) and spatial sampling data to estimate the lateral extent of contaminated soil irrespective of depth (i.e., gross contamination extent) and define areas of contaminant uncertainty. The software performed a binary transformation of contaminant concentrations at all sampling points based upon remedial action goals or a sum of ratios approach (i.e., clean, impacted, or range of impacts in soil). The model produced two-dimensional (horizontal) contaminant probability contours and statistical uncertainty in the sampling coverage and resulting contaminant extents. This method was translated vertically by partitioning the sampling data into depth brackets that produced a stacked representation of contaminant extents and uncertainty in the subsurface (i.e., similar to construction lifts). The results commonly led to a better understanding of project uncertainty and the need for sampling strategies that produce high-confidence soil volumes, which control costs. The BAASS-based delineations were eventually replaced by Empirical Bayesian Kriging (EBK) methods available in ArcGIS Spatial or 3D Analysts [2]. The EBK method calculates contaminant probability zones derived from user-controlled semivariograms of the spatial datasets. The resulting probability zones (e.g., 50% or 80% of contaminant probability) represent the two-dimensional surface delineation of the overall horizontal remedial area, similarly to BAASS. However, unlike BAASS, the vertical sampling data within these probability zones became vertical control points to contour a subterranean surface that connects subsurface points to the land-surface delineations of contamination. The resulting representation of horizontal and vertical impacts within an enclosed envelop (volume) of soil included uncertainty distributions that are used to plan uncertainty-reduction sampling. These data-driven and math-based models of three-dimensional sampling results produced well-bounded remedial volumes for project planning and better uncertainty predictions during project budgeting. The EBK method was applied to several FUSRAP sites managed by the Buffalo District and compared to less rigorously modeled sites previously remediated by the District. The comparison of modeled to actual remediated volumes provide a basis for validating the volume-estimation method. This comparison is important to ensure modeled volumes match physical boundaries of site remediation. FUSRAP sites with denser investigative sampling and lesser volume uncertainty proved useful in remedial planning and contracting. The Buffalo District noted that historical sites with sparser sampling arrays had greater disparity between estimated volumes and final remedial volumes. The benefit achieved over the cost of detailed soil sampling appears positive for FUSRAP projects, especially where impacts vary widely and appear unbounded by investigation-phase sampling. The subsequent Empirical Bayesian Kriging of contamination coupled with vertical contouring for soil estimations reduces uncertainty in soil volumes or indicates where sampling is required to reduce uncertainty, which together optimize remedial planning and budgeting. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Thaw depth, soil moisture, and vegetation height, Teller and Kougarok sites, Seward Peninsula, Alaska, 2022

Thaw depth, soil moisture, and vegetation height sampled from locations on the Teller MM27, Kougarok MM80, and Kougarok MM82 NGEE-Arctic sites, Seward Peninsula, Alaska. These data were collected in support of the ongoing NGEE-Arctic and NASA ABoVE data synthesis work. Samples were collected in July 2022, including 6 transects covering the entire Teller MM27 watershed and 4 transects covering 4 thaw ponds at Kougarok MM80 and MM82. This data package includes sample information and thaw depth, soil moisture, vegetation height data (.csv). Metadata files include data descriptions (_dd.csv) for tabular data and file level metadata (.csv). The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research. The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska. Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗

Mitigating the optical depth degeneracy using the kinematic Sunyaev-Zel’dovich effect with CMB-S4 data

The epoch of reionization is one of the major phase transitions in the history of the universe, and is a focus of ongoing and upcoming cosmic microwave background (CMB) experiments with improved sensitivity to small-scale fluctuations. Reionization also represents a significant contaminant to CMB-derived cosmological parameter constraints, due to the degeneracy between the Thomson-scattering optical depth, $\tau$, and the amplitude of scalar perturbations, $A_s$. This degeneracy subsequently hinders the ability of large-scale structure data to constrain the sum of the neutrino masses, a major target for cosmology in the 2020s. In this work, we explore the kinematic Sunyaev-Zel'dovich (kSZ) effect as a probe of reionization, and show that it can be used to mitigate the optical depth degeneracy with high-sensitivity, high-resolution data from the upcoming CMB-S4 experiment. We discuss the dependence of the kSZ power spectrum on physical reionization model parameters, as well as on empirical reionization parameters, namely $\tau$ and the duration of reionization, $\Delta z$. We show that by combining the kSZ two-point function and the reconstructed kSZ four-point function, degeneracies between $\tau$ and $\Delta z$ can be strongly broken, yielding tight constraints on both parameters. We forecast $\sigma(\tau) = 0.003$ and $\sigma(\Delta z) = 0.25$ for a combination of CMB-S4 and Planck data, including detailed treatment of foregrounds and atmospheric noise. The constraint on $\tau$ is nearly identical to the cosmic-variance limit that can be achieved from large-angle CMB polarization data. The kSZ effect thus promises to yield not only detailed information about the reionization epoch, but also to enable high-precision cosmological constraints on the neutrino mass.

79 ASTRONOMY AND ASTROPHYSICS↗

Grail Refinements to Lunar Seismic Structure

To probe a planet's interior, seismology provides the most direct constraints on the variables that govern the dynamic properties of the body. However, the GRAIL (Gravity Recovery and Interior Laboratory) mission's high-resolution measurements of the lunar gravity field provide constraints on crustal thickness, mantle structure, core radius and stratification, and core state (solid vs. molten). These data complement seismic investigations, and joint interpretation permits improved constraints on the Moon's internal structure. Joint interpretation of disparate geophysical datasets helps reduce drawbacks that can result from analyzing them individually. The Apollo seismic network was situated on the lunar nearside surface in a roughly equilateral triangle having sides approximately 1000 km long, with stations 12/14 nearly co-located at one corner. Due to this limited geographical extent, near-surface ray coverage from moonquakes is low, but increase with depth. In comparison, gravity surveys and their resulting gravity anomaly maps have traditionally offered optimal resolution at crustal depths. Gravimetric maps and seismic data sets are therefor well suited to joint inversion, since the complementary information reduces inherent model ambiguity. We will perform a joint inversion of Apollo seismic delay times and gravity data collected by GRAIL lunar gravity mission, in order to recover seismic velocity and density as a function of latitude, longitude and depth within the Moon. We will relate density (rho) to seismic velocity (v) using a linear relationship that is allowed to be depth-dependent. The corresponding coefficient (B) can reflect a variety of material properties that vary with depth, including temperature and composition. The inversion seeks to recover the set of rho, v, and B perturbations that minimize (in a least-squares sense) the difference between the observed and calculated data.

Weber, Renee C.↗

Investigations of Desert Dust and Smoke in the North Atlantic in Support of the TOMS Instrument

During the initial period of the work we concentrated on Saharan dust storms and published a sequence of papers (Colarco et a1 2002,2003a,b, Toon, 2004). The U.S. Air Force liked the dust model so well that they appropriated it for operational dust storm forecasting (Barnum et al., 2004). The Air Force has used it for about 5 yrs in the Middle East where dust storms cause significant operational problems. The student working on this project, Peter Colarco, has graduated and is now a civil servant at Goddard where he continues to interact with the TOMS team. This work helped constrain the optical properties of dust at TOMS wavelengths, which is useful for climate simulations and for TOMS retrievals of dust properties such as optical depth. We also used TOMS data to constrain the sources of dust in Africa and the Middle East, to determine the actual paths taken by Saharan dust storms, to learn more about the mechanics of variations in the optical depths, and to learn more about the mechanisms controlling the altitudes of the dust. During the last two years we have been working on smoke from fires. Black carbon aerosols are one of the leading factors in radiative forcing. The US Climate Change Science Program calls this area out for specific study. It has been suggested by Jim Hansen, and Mark Jacobsen among others, that by controlling emissions of black carbon we might reduce greenhouse radiative forcing in a relatively painless manner. However, we need a greatly improved understanding of the amount of black carbon in the atmosphere, where it is located, where it comes from, how it is mixed with other particles, what its actual optical properties are, and how it evolves. In order to learn about these issues we are using a numerical model of smoke. We have applied this model to the SAFARI field program data, and used the TOMS satellite observations in that period (Sept. 2000). Our goal is to constrain source function estimates for black carbon, and smoke optical properties.

Toon, Owen B.↗

Melt Inclusions in Zircon: A Window to Understanding the Structure and Evolution of the Magmatic System Beneath the Laguna del Maule Volcanic Field

Explosive silicic eruptions pose a significant threat to society, yet the development and destabilization of the underlying silicic magmatic systems are still controversial. Zircons provide simultaneous information on the trace element composition and age of silicic magmatic systems, while melt inclusions in quartz and plagioclase yield important constraints on their volatile content as well as storage depth. Melt inclusions in zircons (MIZs) combine these data from a single mineral grain, recording the age, storage depth, temperature, and composition of magmas, and thus provide unique constraints on the structure and evolution of silicic magmatic systems. We studied MIZs from the Laguna del Maule (LdM) volcanic field in the southern Andes that is among the most active Pleistocene-Holocene rhyolitic volcanic centers worldwide and a potentially hazardous system displaying inflation rates in excess of 25 cm/yr. The host zircon ages suggest that the LdM MIZ record extends to ~30 kyr before eruption, in contrast to the melt inclusions in LdM plagioclase and quartz crystals that formed only decades to centuries before eruption. The major element compositions of MIZs are minimally affected by post-entrapment crystallization, and agree well with the LdM rhyolitic whole rock data. The more evolved major element composition of rle MIZs than rdm MIZs, suggests a long-term deeper connection of the rdm crystal mush to a more primitive magma body than that of the rle. The evidence of slow H diffusion observed in MIZs suggest that their H 2 O contents are not significantly affected by diffusion of H through the host zircon. The storage depths of 1.1 to 2.8 kbars recorded by the H 2 O contents of rdm and rle MIZs are consistent with the optimal emplacement window (2.0 ± 0.5 kbar) of silicic magma reservoir growth, storage, and eruptibility based on thermomechanical modeling (Huber et al. 2019).

Kei Shimizu↗