Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Spatial process”

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 1,765 records · Page 98

Improved Earth Oblateness Rate Reveals Increased Ice Sheet Losses and Mass-Driven Sea Level Rise

Satellite laser ranging (SLR) observations are routinely applied toward the estimation of dynamic oblateness, C(sub 20), which is the largest globally integrated component of Earth's time-variable gravity field. Since 2002, GRACE and GRACE Follow-On have revolutionized the recovery of higher spatial resolution features of global time-variable gravity, with SLR continuing to provide the most reliable estimates of C (sub 20).We quantify the effect of various SLR processing strategies on estimating C(sub 20) and demonstrate better signal recovery with the inclusion of GRACE-derived low-degree gravity information in the forward model. This improved SLR product modifies the Antarctic and Greenland Ice Sheet mass trends by -15.4 and -3.5 Gt/year, respectively, as compared to CSR TN11, and improves global mean sea level budget closure by modifying sea level rise by +0.08 mm/year. We recommend that this new C(sub 20) product be applied to RL06 GRACE data products for enhanced accuracy and scientific interpretation.

Loomis, B. D.↗

Cluster Analysis of Spectroscopic Line Profiles and EUV Emission in RMHD Simulations and Observations of the Solar Atmosphere

Spatially-resolved observations from the IRIS, SDO/AIA, and other space mission and ground-based telescopes, coupled with realistic 3D RMHD simulations, are a powerful tool for analysis of processes in the solar atmosphere. To better understand the dynamical and thermodynamic properties in the simulation data and their connection to observations, it is essential to determine similarities in the behaviors of the synthesized and observed emission. However, the complexity of observational data and physical processes makes comparison of observations and modeling results difficult. In this work, we show the initial results of application of K-Means clustering (unsupervised machine learning) algorithm to two different problems: 1) recognition of the typical spectroscopic line profiles observed by IRIS during solar flares and their typical dynamic behavior; 2) recognition of shocks and heating events in synthetic AIA emission data obtained from StellarBox quiet-Sun simulations. The average silhouette width technique for the KMeans algorithm is utilized in different ways to obtain optimal numbers of clusters. We discuss application of the emission clustering to visualizations of the computational volume, understanding its evolutionary trends and behavior patterns, and inversion (reconstruction) of physical properties of the solar atmosphere from synthesizes emission data.

Sadykov, Viacheslav↗

Development of a “Nature Run” for Observing System Simulation Experiments (OSSEs) for Snow Mission Development

Snow is a fundamental component of global and regional water budgets, particularly in mountainous areas and regions downstream that rely on snowmelt for water resources. Land surface models (LSMs) are commonly used to develop spatially distributed estimates of snow water equivalent (SWE) and runoff. However, LSMs are limited by uncertainties in model physics and parameters, among other factors. In this study, we describe the use of model calibration tools to improve snow simulations within the Noah-MP LSM as the first step in an observing system simulation experiment (OSSE). Noah-MP is calibrated against the University of Arizona (UA) SWE product over a western Colorado domain. With spatially varying calibrated parameters, we run calibrated and default Noah-MP simulations for water years 2010–20. By evaluating both simulations against the UA dataset, we show that calibration decreases domain averaged temporal RMSE and bias for snow depth from 0.15 to 0.13 m and from −0.036 to −0.0023 m, respectively, and improves the timing of snow ablation. Increased snow simulation performance also improves estimates of model-simulated runoff in four of six study basins, though only one has statistically significant improvement. Spatially distributed Noah-MP snow parameters perform better than default uniform values. We demonstrate that calibrating variables related to snow albedo calculations and rain–snow partitioning, among other processes, is a necessary step for creating a nature run that reasonably approximates true snow conditions for the OSSEs. Additionally, the inclusion of a snowfall scaling term can address biases in precipitation from meteorological forcing datasets, further improving the utility of LSMs for generating reliable spatiotemporal estimates of snow.

Melissa L. Wrzesien↗

Hyperspectral Sounder Products: SiFSAP and ClimFiSP

Infrared hyperspectral sounders have demonstrated to provide very high value information for a range of weather and climate applications. Traditional retrieval algorithms of hyperspectral sounders use the ‘cloud clearing’ methodology in the retrieval process due to the lack of means to account for cloud scattering effect. The corresponding Level-2 data products are therefore limited by two major factors: the degradation of spatial resolution as compared with the native resolution of the instruments and the lack of ‘radiance closure’ that is critically needed for climate trend studies. The single field-of-view (FOV) sounder atmospheric products (SiFSAP) algorithm has been developed to supplement existing operational products. The SiFSAP algorithm uses the principal component based radiative transfer model (PCRTM) that is ultra-fast and has been extensively validated for the cloud scattering simulation capability. The cloud properties are physically retrieved along with profiles of temperature, moisture, and trace gases of interest, surface properties simultaneously by fitting single FOV spectral measurements. Therefore, SiFSAP ensure the radiometric consistency between the retrieved geophysical properties and the top-of-atmosphere (TOA) spectral measurements. This presentation provides a general introduction of SiFSAP and details about the physical retrieval algorithm. Also introduced here is a derivative data product of SiFSAP, the climate fingerprinting Sounder Products (ClimFiSP). ClimFiSP include the space-time averaged properties of key climate variables that are derived from the space-time averaged radiances via the spectral fingerprinting scheme. The spectral fingerprinting method allows a low latency data update on ClimFiSP following any potential calibration update on satellite measurements. ClimFiSP also include the rigorously defined uncertainty characterization for trends or anomalies of climate variables, being critical for the construction of long-term climate data record. Both SiFSAP and ClimFiSP will be available to users through NASA's Goddard Earth Sciences Data and Information Services Center (GES DISC).

Wan Wu↗

Spatial displays as a means to increase pilot situational awareness

Experiences raise a number of concerns for future spatial-display developers. While the promise of spatial displays is great, the cost of their development will be correspondingly large. The knowledge and skills which must be coordinated to ensure successful results is unprecedent. From the viewpoint of the designer, basic knowledge of how human beings perceive and process complex displays appears fragmented and largely unquantified. Methodologies for display development require prototyping and testing with subject pilots for even small changes. Useful characterizations of the range of differences between individual users is nonexistent or at best poorly understood. The nature, significance, and frequency of interpretation errors associated with complex integrated displays is unexplored and undocumented territory. Graphic displays have intuitive appeal and can achieve face validity much more readily than earlier symbolic displays. The risk of misleading the pilot is correspondingly greater. Thus while some in the research community are developing the tools and techniques necessary for effective spatial-display development, potential users must be educated about the issues so that informed choices can be made. The scope of the task facing all is great. The task is challenging and the potential for meaningful contributions at all levels is high indeed.

Fadden, Delmar M.↗

Introducing Multisensor Satellite Radiance-Based Evaluation for Regional Earth System Modeling

Earth System modeling has become more complex, and its evaluation using satellite data has also become more difficult due to model and data diversity. Therefore, the fundamental methodology of using satellite direct measurements with instrumental simulators should be addressed especially for modeling community members lacking a solid background of radiative transfer and scattering theory. This manuscript introduces principles of multisatellite, multisensor radiance-based evaluation methods for a fully coupled regional Earth System model: NASA-Unified Weather Research and Forecasting (NU-WRF) model. We use a NU-WRF case study simulation over West Africa as an example of evaluating aerosol-cloud-precipitation-land processes with various satellite observations. NU-WRF-simulated geophysical parameters are converted to the satellite-observable raw radiance and backscatter under nearly consistent physics assumptions via the multisensor satellite simulator, the Goddard Satellite Data Simulator Unit. We present varied examples of simple yet robust methods that characterize forecast errors and model physics biases through the spatial and statistical interpretation of various satellite raw signals: infrared brightness temperature (Tb) for surface skin temperature and cloud top temperature, microwave Tb for precipitation ice and surface flooding, and radar and lidar backscatter for aerosol-cloud profiling simultaneously. Because raw satellite signals integrate many sources of geophysical information, we demonstrate user-defined thresholds and a simple statistical process to facilitate evaluations, including the infrared-microwave-based cloud types and lidar/radar-based profile classifications.

Planetary Boundary Layer↗

A Goddard Multi-Scale Modeling System with Unified Physics

Numerical cloud resolving models (CRMs), which are based the non-hydrostatic equations of motion, have been extensively applied to cloud-scale and mesoscale processes during the past four decades. Recent GEWEX Cloud System Study (GCSS) model comparison projects have indicated that CRMs agree with observations in simulating various types of clouds and cloud systems from different geographic locations. Cloud resolving models now provide statistical information useful for developing more realistic physically based parameterizations for climate models and numerical weather prediction models. It is also expected that Numerical Weather Prediction (NWP) and regional scale model can be run in grid size similar to cloud resolving model through nesting technique. Current and future NASA satellite programs can provide cloud, precipitation, aerosol and other data at very fine spatial and temporal scales. It requires a coupled global circulation model (GCM) and cloud-scale model (termed a szrper-parameterization or multi-scale modeling -framework, MMF) to use these satellite data to improve the understanding of the physical processes that are responsible for the variation in global and regional climate and hydrological systems. The use of a GCM will enable global coverage, and the use of a CRM will allow for better and more sophisticated physical parameterization. NASA satellite and field campaign can provide initial conditions as well as validation through utilizing the Earth Satellite simulators. At Goddard, we have developed a multi-scale modeling system with unified physics. The modeling system consists a coupled GCM-CRM (or MMF); a state-of-the-art weather research forecast model (WRF) and a cloud-resolving model (Goddard Cumulus Ensemble model). In these models, the same microphysical schemes (2ICE, several 3ICE), radiation (including explicitly calculated cloud optical properties), and surface models are applied. In addition, a comprehensive unified Earth Satellite simulator has been developed at GSFC, which is designed to fully utilize the multi-scale modeling system. A brief review of the multi-scale modeling system with unified physics/simulator and examples is presented in this article.

Tao, W.K.↗

Pyroscopegridding: Geo-Leo Aerosol Data Fusion (an Open-Source Package)

The retrieval of aerosol optical depths (AODs) from sun-synchronous polar orbiting satellites, such as MODISs, VIIRSs, OMI, TROPOMI, etc., has been widely adopted as a method for obtaining information regarding particulate matter (PM) and related atmospheric processes. However, the advent of recently launched geostationary satellites, such as GOES-16/17/18, Himawari-8/9, and Meteosat Third Generation (MTG), has led to an increased temporal resolution of AOD observations (order of 10 minutes), resulting in typically one or more images per hour during daylight hours, compared to the once-per-day observations obtained from LEO satellites. By integrating these observations, the diurnal cycle of global AOD can be characterized at local, regional, and global scales. The scientific community is still examining the novel data from geostationary satellite observations and evaluating methods for effectively merging these observations with differing spatial and temporal resolutions. This presents a significant ""Big Data"" challenge, encompassing not only data storage, but also data discoverability, accessibility, and migration within cloud computing environments. This study presents our attempts at fusing Level 2 aerosol data from six satellites, three of which are geostationary (GOES-16/17 and Himawari-8) and three of which are polar orbiting (TERRA/MODIS, AQUA/MODIS, and SNPP-VIIRS), using the Dark Target aerosol retrieval algorithm. The ability to fuse remote sensing products on demand into desired temporal and spatial domains empowers researchers and practitioners to more efficiently work with satellite and sensor data. It is our hope that through making our open-source package and accompanying functionality available, the scientific community will have improved access to aerosol data processing resources.

Jennifer Wei↗

The cold land processes experiment (CLPX) local scale observatin site (LSOS)

The Local Scale Observation Site (LSOS) is the smallest study site of the Cold LandProcesses Experiment (CLPX) and is located within the Fraser Meso-cell Study Area (MSA), near the Fraser Experimental Forest Headquarters Facility, in Fraser, CO USA.The 100-m x 100-m site consists of a small open field, a managed dense canopy and an open, mixed age canopy. Unlike the other components of the experiment, which focus on spatial distributions at relatively brief snapshots in time, measurements at the local scale site focused on the temporal domain.

cold land processes ground-based remote sensing va↗

Interstellar Isomers: The Importance of Bonding Energy Differences

We present strong detections of methyl cyanide (CH3CN), vinyl cyanide (CH2CHCN), ethyl cyanide (CH3CH2CN) and cyanodiacetylene (HC4CN) molecules with the Green Bank Telescope (GBT) toward the Sgr B2(N) molecular cloud. Attempts to detect the corresponding isocyanide isomers were only successful in the case of methyl isocyanide (CH3NC) for its J(sub K) = 1(sub 0) - 0(sub 0) transition, which is the first interstellar report of this line. To determine the spatial distribution of CH3NC, we used archival Berkeley-Illinois-Maryland Association (BIMA) array data for the J(sub K) = 1(sub 0) - 0(sub 0) transitions but no emission was detected. From ab initio calculations, the bonding energy difference between the cyanide and isocyanide molecules is greater than 8500 per centimeter (greater than 12,000 K). Thus, cyanides are the more stable isomers and would likely be formed more preferentially over their isocyanide counterparts. That we detect CH3NC emission with a single antenna (Gaussian beamsize(omega(sub B))=1723 arcsec(sup 2)) but not with an interferometer (omega(sub b)=192 arcsec(sup 2)), strongly suggests that CH3NC has a widespread spatial distribution toward the Sgr B2(N) region. Other investigators have shown that CH3CN is present both in the LMH hot core of Sgr B2(N) and in the surrounding medium, while we have shown that CH3NC appears to be deficient in the LMH hot core. Thus, largescale, non-thermal processes in the surrounding medium may account for the conversion of CH3CN to CH3NC while the LMH hot core, which is dominated by thermal processes, does not produce a significant amount of CH3NC. Ice analog experiments by other investigators have shown that radiation bombardment of CH3CN can produce CH3NC, thus supporting our observations. We conclude that isomers separated by such large bonding energy differences are distributed in different interstellar environments, making the evaluation of column density ratios between such isomers irrelevant unless it can be independently shown that these species are co-spatial.

Remijan, Anthony J.↗

Low-Level Jets and the Convergence of Mars Data Assimilation Algorithms

Data assimilation is an increasingly popular technique in Mars atmospheric science, but its effect on the mean states of the underlying atmosphere models has not been thoroughly examined. The robustness of results to the choice of model and assimilation algorithm also warrants further study. We investigate these issues using two Mars general circulation models (MGCMs), with particular emphasis on zonal wind and temperature fields. When temperature retrievals from the Mars Global Surveyor Thermal Emission Spectrometer (TES) are assimilated into the U.K.-Laboratoire de M ́et ́eorologie Dynamique (UK-LMD) MGCM to create the Mars Analysis Correction Data Assimilation (MACDA) reanalysis, low-level zonal jets in the winter northern hemisphere shift equatorward and weaken relative to a free-running control simulation from the same MGCM. The Ensemble Mars Atmosphere Reanalysis System (EMARS) reanalysis, which is also based on TES temperature retrievals, also shows jet weakening (but less if any shifting) relative to a control simulation performed with the underlying Geophysical Fluid Dynamics Laboratory (GFDL) MGCM. Examining higher levels of the atmosphere, monthly mean three-dimensional temperature and zonal wind fields are in generally better agreement between the two reanalyses than between the two control simulations. In conjunction with information about the MGCMs’ physical parametrizations, intercomparisons between the various reanalyses and control simulations suggest that overall the EMARS control run is plausibly less biased (relative to the true state of the Martian atmosphere) than the MACDA control run. Implications for future observational studies are discussed. Plain Language Summary An increasingly popular way to study Martian weather and climate is to combine atmospheric temperature observations with a computer model (specifically, a global climate model). The process of combining model and observations is called “data assimilation”, and the resulting merged data set is called a “reanalysis”. One advantage of reanalyses is that they include variables (such as wind) that are not directly observed. For scientific and practical applications we want these variables to be reasonably accurate— however, it is not clear how well data assimilation algorithms compute them. Our study investigates this issue using two Mars reanalyses and two model simulations that do not assimilate temperature data. We focus on slowly-varying atmospheric phenomena (timescales from 10 Mars days to a season). Assimilating temperature data into two different global climate models changes the strength and/or spatial pattern of east-west winds at low altitudes. Furthermore, monthly mean three-dimensional temperature and east-west wind fields agree better between reanalyses than between non-assimilating model simulations. This suggests that the data assimilation process is basically successful. One non-assimilating model simulation has less realistic representations of atmospheric physical processes than the other—we argue that this plausibly gives it larger biases relative to the true state of the atmosphere.

Todd A Mooring↗

Motion sickness and development of synergy within the spatial orientation system. A hypothetical unifying concept

Adaptation to research paradigms such as rotating rooms and optical alteration of visual feedback during movement results in development of perceptual-motor programs that provide the reflexive assistance that is necessary to skilled control of movement and balance. The discomfort and stomach awareness that occur during the adaptation process has been attributed to conflicting sensory information about the state of motion. Vestibular signals depend on the kinematics of head movements irrespective of the presence or absence of signals from other senses. We propose that sensory conflict when vestibular signals are at least one component of the conflict are innately disturbing and unpleasant. This innate reaction is part of a continuum that operates early in life to prevent development of inefficient perceptual-motor programs. This reaction operates irrespective of and in addition to reward and punishment from parental guidance or goal attainment to yield efficient control of whole body movement in the operating environment of the individual. The same mechanism is involved in adapting the spatial orientation system to strange environments. This conceptual model "explains" why motion sickness is associated with adaptation to novel environments and is in general consistent with motion sickness literature.

Review, Tutorial↗

How Cities Make Their Own Weather

Urbanization is one of the extreme cases of land use change. Most of world's population has moved to urban areas. Although currently only 1.2% of the land is considered urban, the spatial coverage and density of cities are expected to rapidly increase in d e near future. It is estimated that by the year 2025, 60% of the world's population will live in cities. Human activity in urban environments also alters weather and climate processes. However, our understanding of urbanization on the total Earth-weather-climate system is incomplete. Recent literature continues to provide evidence that anomalies in precipitation exist over and downwind of major cities. Current and future research efforts are actively seeking to verify these literature findings and understand potential cause-effect relationships. The novelty of this study is that it utilizes rainfall data from multiple satellite data sources (e.g. TRMM precipitation radar, TRMM-geosynchronous-rain gauge merged product, and SSM/I) and ground-based measurements to identify spatial anomalies and temporal trends in precipitation for cities around the world. We will also present results from experiments using a regional atmospheric-land surface modeling system. Early results will be presented and placed within the context of weather prediction, climate assessment, and societal applications.

Shepherd, J. Marshall↗

Mesoscale Eddies, Satellite Altimetry, and New Production in the Sargasso Sea

Satellite altimetry and hydrographic observations are used to characterize the mesoscale eddy field in the Sargasso Sea near Bermuda and to address the role of physical processes on the supply of new nutrients to the euphotic zone. The observed sea level anomaly (SLA) field is dominated by the occurrence of westward propagating features with SLA signatures as large as 25 cm, Eulerian temporal scales of roughly a month, lifetimes of several months, spatial scales of approximately 200 km, and a propagation of approximately 5 cm/s . Hydrographic estimates of dynamic height anomaly (referenced to 4000 dbar) are well correlated with satellite SLA (r(sup 2) = 0.65), and at least 85% of the observed dynamic height variability is associated with the first baroclinic mode of motion. This allows us to apply the satellite observations to remotely estimate isopycnal displacements and the flux of nutrients into the euphotic zone due to eddy pumping. Eddy pumping is the process by which mesoscale eddies induce isopycnal displacements that lift nutrient- replete waters into the euphotic zone, driving new primary production. A kinematic approach to the estimation of the eddy pumping results in a flux of 0.24+/-0.1 mol N/sq m/yr (including a scale estimate for the small contribution due to 18 deg water eddies). This flux is more than an order of magnitude larger than the diapycnal diffusive flux as well as scale estimates for the vertical transport due to isopycnal mixing along sloping isopycnal surfaces. Eddy pumping and wintertime convection are the two dominant mechanisms transporting new nutrients into the euphotic zone, and the sum of all physical new nutrient supply fluxes effectively balances previous geochemical estimates of annual new production for this site. However, if biological transports (e.g., nitrogen fixation, etc.) are significant, the new nitrogen supply budget will be in excess of geochemical new production estimates. This suggests that the various physical and biological transport fluxes, as well as geochemical inferences of new production, still need to be reconciled and many outstanding questions remain.

Siegel, David A.↗

Mesoscale Eddies, Satellite Altimetry, and New Production in the Sargasso Sea

Satellite altimetry and hydrographic observations are used to characterize the mesoscale eddy field in the Sargasso Sea near Bermuda and to address the role of physical processes on the supply of new nutrients to the euphotic zone. The observed sea level anomaly (SLA) field is dominated by the occurrence of westward propagating features with SLA signatures as large as 25 cm, Eulerian temporal scales of roughly a month, lifetimes of several months, spatial scales of approximately 200 km, and a propagation of approximately 5 cm/s. Hydrographic estimates of dynamic height anomaly (referenced to 4000 dbar) are well correlated with satellite SLA (r(exp 2) = 0.65), and at least 85% of the observed dynamic height variability is associated with the first baroclinic mode of motion. This allows us to apply the satellite observations to remotely sensed estimate isopycnal displacements and the flux of nutrients into the euphotic zone due to eddy pumping. Eddy pumping is the process by which mesoscale eddies induce isopycnal displacements that lift nutrient-replete waters into the euphotic zone, driving new primary production. A kinematic approach to the estimation of the eddy pumping results in a flux of 0.24 +/- 0.1 mol N/sq m (including a scale estimate for the small contribution due to 18 deg water eddies). This flux is more than an order of magnitude larger than the diapycnal diffusive flux as well as scale estimates for the vertical transport due to isopycnal mixing along sloping isopycnal surfaces. Eddy pumping and wintertime convection are the two dominant mechanisms transporting new nutrients into the euphotic zone, and the sum of all physical new nutrient supply fluxes effectively balances previous geochemical estimates of annual new production for this site. However, if biological transports (e.g., nitrogen fixation, etc.) are significant, the new nitrogen supply budget will be in excess of geochemical new production estimates. This suggests that the various physical and biological transport fluxes, as well as geochemical inferences of new production, still need to be reconciled and many outstanding questions remain.

Siegel, David A.↗

An Analysis of Gravity-Field Estimation Based on Intersatellite Dual-1-Way Biased Ranging

The GRACE (Gravity Recovery And Climate Experiment) mission is designed to make global, highly accurate measurements of the Earth's gravity field with high spatial resolution. Ancillary GPS occultation measurements are also to be carried out for atmospheric monitoring. In the dual-1-way biased ranging of this mission, the range between two satellites separated by 100 to 200 km in nearly polar, coplanar, circular orbits, is measured to very high precision, to within an additive constant, through the exchange of K- and Ka-band sinusoidal signals. Such biased ranging data, along with GPS L-band range and phase data, can be processed and fit over successive multiday intervals to obtain accurate estimates of the Earth's gravity field. This report approximately models and analyzes this process, from the generation of the RF signals at the two satellites through the extraction of the geopotential. The steps include generation of the transmitted signals, processing the received signals to extract high-rate baseband phase, carrying out a dual-1-way combination of baseband phase to extract high-rate biased range for each band, combining K- and Ka-band ranges to correct for the ionosphere effect, and processing the resulting high-rate biased range values to extract three types of reduced-rate observables: biased range, range rate and range acceleration. The version of dual-1-way biased ranging developed by this report improves upon previous versions in a number of ways: highly accurate satellite-timetag corrections derived from concurrent GPS data, better baseband phase extraction using highly digital processing, highly accurate USO-rate calibration derived from concurrent GPS data, an improved method for extracting high-rate biased range from baseband phase, improved filtering for extracting reduced- rate observables from high-rate biased range, and parallel extraction of three observable types.

SINUSOIDAL SIGNALS↗

GMI-IPS: Python Processing Software for Aircraft Campaigns

NASA's Atmospheric Tomography Mission (ATom) seeks to understand the impact of anthropogenic air pollution on gases in the Earth's atmosphere. Four flight campaigns are being deployed on a seasonal basis to establish a continuous global-scale data set intended to improve the representation of chemically reactive gases in global atmospheric chemistry models. The Global Modeling Initiative (GMI), is creating chemical transport simulations on a global scale for each of the ATom flight campaigns. To meet the computational demands required to translate the GMI simulation data to grids associated with the flights from the ATom campaigns, the GMI ICARTT Processing Software (GMI-IPS) has been developed and is providing key functionality for data processing and analysis in this ongoing effort. The GMI-IPS is written in Python and provides computational kernels for data interpolation and visualization tasks on GMI simulation data. A key feature of the GMI-IPS, is its ability to read ICARTT files, a text-based file format for airborne instrument data, and extract the required flight information that defines regional and temporal grid parameters associated with an ATom flight. Perhaps most importantly, the GMI-IPS creates ICARTT files containing GMI simulated data, which are used in collaboration with ATom instrument teams and other modeling groups. The initial main task of the GMI-IPS is to interpolate GMI model data to the finer temporal resolution (1-10 seconds) of a given flight. The model data includes basic fields such as temperature and pressure, but the main focus of this effort is to provide species concentrations of chemical gases for ATom flights. The software, which uses parallel computation techniques for data intensive tasks, linearly interpolates each of the model fields to the time resolution of the flight. The temporally interpolated data is then saved to disk, and is used to create additional derived quantities. In order to translate the GMI model data to the spatial grid of the flight path as defined by the pressure, latitude, and longitude points at each flight time record, a weighted average is then calculated from the nearest neighbors in two dimensions (latitude, longitude). Using SciPya's Regular Grid Interpolator, interpolation functions are generated for the GMI model grid and the calculated weighted averages. The flight path points are then extracted from the ATom ICARTT instrument file, and are sent to the multi-dimensional interpolating functions to generate GMI field quantities along the spatial path of the flight. The interpolated field quantities are then written to a ICARTT data file, which is stored for further manipulation. The GMI-IPS is aware of a generic ATom ICARTT header format, containing basic information for all flight campaigns. The GMI-IPS includes logic to edit metadata for the derived field quantities, as well as modify the generic header data such as processing dates and associated instrument files. The ICARTT interpolated data is then appended to the modified header data, and the ICARTT processing is complete for the given flight and ready for collaboration. The output ICARTT data adheres to the ICARTT file format standards V1.1. The visualization component of the GMI-IPS uses Matplotlib extensively and has several functions ranging in complexity. First, it creates a model background curtain for the flight (time versus model eta levels) with the interpolated flight data superimposed on the curtain. Secondly, it creates a time-series plot of the interpolated flight data. Lastly, the visualization component creates averaged 2D model slices (longitude versus latitude) with overlaid flight track circles at key pressure levels. The GMI-IPS consists of a handful of classes and supporting functionality that have been generalized to be compatible with any ICARTT file that adheres to the base class definition. The base class represents a generic ICARTT entry, only defining a single time entry and 3D spatial positioning parameters. Other classes inherit from this base class; several classes for input ICARTT instrument files, which contain the necessary flight positioning information as a basis for data processing, as well as other classes for output ICARTT files, which contain the interpolated model data. Utility classes provide functionality for routine procedures such as: comparing field names among ICARTT files, reading ICARTT entries from a data file and storing them in data structures, and returning a reduced spatial grid based on a collection of ICARTT entries. Although the GMI-IPS is compatible with GMI model data, it can be adapted with reasonable effort for any simulation that creates Hierarchical Data Format (HDF) files. The same can be said of its adaptability to ICARTT files outside of the context of the ATom mission. The GMI-IPS contains just under 30,000 lines of code, eight classes, and a dozen drivers and utility programs. It is maintained with GIT source code management and has been used to deliver processed GMI model data for the ATom campaigns that have taken place to date.

Damon, M. R.↗

Quantifying Libya-4 Surface Reflectance Heterogeneity With WorldView-1, 2 and EO-1 Hyperion

The land surface imaging (LSI) virtual constellation approach promotes the concept of increasing Earth observations from multiple but disparate satellites. We evaluated this through spectral and spatial domains, by comparing surface reflectance from 30-m Hyperion and 2-m resolution WorldView-2 (WV-2) data in the Libya-4 pseudoinvariant calibration site. We convolved and resampled Hyperion to WV-2 bands using both cubic convolution and nearest neighbor (NN) interpolation. Additionally, WV-2 and WV-1 same-date imagery were processed as a cross-track stereo pair to generate a digital terrain model to evaluate the effects from large (>70 m) linear dunes. Agreement was moderate to low on dune peaks between WV-2 and Hyperion (R2 <; 0.4) but higher in areas of lower elevation and slope (R2 > 0.6). Our results provide a satellite sensor intercomparison protocol for an LSI virtual constellation at high spatial resolution, which should start with geolocation of pixels, followed by NN interpolation to avoid tall dunes that enhance surface reflectance differences across this internationally utilized site.

Hyperion↗