Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Spatial statistics”

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,495 records · Page 83

Global Ocean Mass Change Derived From Atmospheric Reanalyses

Barystatic sea level change is a significant climate variable and a gauge for water exchanges between land, atmosphere, and ocean, with consequential relevance to stakeholders. The monitoring and attribution of globally averaged sea level variations represents substantial challenges to the climate community. Atmospheric reanalyses provide a level of monitoring for elements of the climate system that are in communication with the ocean, including the atmosphere, terrestrial land, and the surfaces of glaciated land. Global atmospheric reanalyses comprise model output that is statistically adjusted to available observations. Potentially, reanalyses offer a climatic context to observed changes in water mass distribution. Improvement in the representation of the water cycle and in water mass conservation have been stated goals of reanalysis projects. Here, we assess the capability of atmospheric reanalyses to represent variations in the global water budget from the perspective of ocean mass balance. The purpose of this work is to provide a detailed evaluation of current atmospheric reanalyses and their ability to reproduce spatial and temporal variations in atmospheric and terrestrial water storage. A second goal is to identify and characterize recent apparent trends in the water balance. We examine the ERA5, MERRA-2, and JRA-3Q reanalyses with a focus on the radio-occultation observing period from 2007 to the present. For terrestrial water storage we also compare with ERA5-Land, a surface reanalysis. The primary validation is from the NASA/DLR Gravity Recovery and Climate Experiment (GRACE) and its follow-on, GRACE-FO. The Atmospheric Infrared Sounder (AIRS) onboard NASA's Aqua satellite is also used for evaluation of the atmospheric moisture budget. Comparisons reveal strengths and weaknesses in the reanalyses of each component of the climate system. Reanalyses agree closely on atmospheric moisture storage trends, although mean values differ. Reanalyses show significant differences in land water storage. They also apply large corrections to the land water budget. Finally, reanalyses have differing abilities in representing the surface mass balance over glaciated land. Using the observed gravimetry record, changes in the underlying glaciated land ice are examined as a residual.

Richard Cullather↗

Analysis of alternative climate datasets and evapotranspiration methods for the Upper Mississippi River Basin using SWAT within HAWQS

This study reports the application of Soil andWater Assessment Tool (SWAT) within the Hydrologic andWater Quality System(HAWQS) on-line platform, for the UpperMississippi River Basin (UMRB). TheUMRB is an important ecosystem located in the north central U.S. that is experiencing a range of ecological stresses. Specifically, testing of SWAT was performed for: (1) Hargreaves (HG) and Penman-Monteith (PM) PET methods, and (2) Livneh, National Climatic Data Center (NCDC) and Parameter-elevation Regressions on Independent Slopes Model (PRISM) climate datasets. The Livneh-PM combination resulted in the highest average annual water yield of 380.6 mm versus the lowest estimated water yield of 193.9 mm for the Livneh-HG combination, in response to 23-year uncalibrated simulations. Higher annual ET and PET values were predicted with HG method versus the PM method for all three weather datasets in response to the uncalibrated simulations, due primarily to higher HG-based estimates during the growing season. Based on these results, it was found that the HG method is the preferred PET option for the UMRB. Initial calibration of SWAT was performed using the Livneh data and HGmethod for threeMississippi Rivermain stem gauge sites, which was followed by spatial validation at 10 other gauge sites located within the UMRB streamnetwork. Overall satisfactory results were found for the calibration and validation gauge sites, with the majority of R2 values ranging between 0.61 and 0.82, Nash- Sutcliffe modeling efficiency (NSE) values ranging between 0.50 and 0.79, and Kling-Gupta efficiency (KGE) values ranging between 0.61 and 0.84. The results of an additional experimental suite of six scenarios,which represented different combinations of climate data sets and calibrated parameters, revealed that suggested statistical criteria were again satisfied by the different scenario combinations. Overall, the PRISM data exhibited the strongest reliability for the UMRB.

54 ENVIRONMENTAL SCIENCES↗

Dynamics of micrometeoroids

The spatial and orbital distributions of meteoroids of mass less than 10 to the -8 g are discussed. Pioneer 8 and 9, HEOS-2 and Helios 1 observations of the smallest meteoroids leaving the solar system in hyperbolic orbits due to radiation pressure (beta particles), intermediate mass particles arriving from the heliocentric spacecraft apex direction (alpha meteoroids) and the largest micrometeoroids observed with no preferred directionality (sporadic meteoroids) are reviewed, and evidence for the three groups of dynamically different micrometeoroids in lunar microcrater statistics is indicated. The collisional model of the origin of dust particles from the dispersion of sporadic meteoroids of cometary origin due to the Poynting-Robertson effect and collisions is then applied to account for the dynamics of the beta and alpha micrometeoroids.

Gruen, G.↗

Full-Coverage High-Resolution Daily PM(sub 2.5) Estimation using MAIAC AOD in the Yangtze River Delta of China

Satellite aerosol optical depth (AOD) has been used to assess population exposure to fine particulate matter (PM (sub 2.5)). The emerging high-resolution satellite aerosol product, Multi-Angle Implementation of Atmospheric Correction(MAIAC), provides a valuable opportunity to characterize local-scale PM(sub 2.5) at 1-km resolution. However, non-random missing AOD due to cloud snow cover or high surface reflectance makes this task challenging. Previous studies filled the data gap by spatially interpolating neighboring PM(sub 2.5) measurements or predictions. This strategy ignored the effect of cloud cover on aerosol loadings and has been shown to exhibit poor performance when monitoring stations are sparse or when there is seasonal large-scale missngness. Using the Yangtze River Delta of China as an example, we present a Multiple Imputation (MI) method that combines the MAIAC high-resolution satellite retrievals with chemical transport model (CTM) simulations to fill missing AOD. A two-stage statistical model driven by gap-filled AOD, meteorology and land use information was then fitted to estimate daily ground PM(sub 2.5) concentrations in 2013 and 2014 at 1 km resolution with complete coverage in space and time. The daily MI models have an average R(exp 2) of 0.77, with an inter-quartile range of 0.71 to 0.82 across days. The overall Ml model 10-fold cross-validation R(exp 2) (root mean square error) were 0.81 (25 gm(exp 3)) and 0.73 (18 gm(exp 3)) for year 2013 and 2014, respectively. Predictions with only observational AOD or only imputed AOD showed similar accuracy.Comparing with previous gap-filling methods, our MI method presented in this study performed bette rwith higher coverage, higher accuracy, and the ability to fill missing PM(sub 2.5) predictions without ground PM(sub 2.5) measurements. This method can provide reliable PM(sub 2.5)predictions with complete coverage that can reduce biasin exposure assessment in air pollution and health studies.

PM2.5↗

Extension of DAMAS Phased Array Processing for Spatial Coherence Determination (DAMAS-C)

The present study reports a new development of the DAMAS microphone phased array processing methodology that allows the determination and separation of coherent and incoherent noise source distributions. In 2004, a Deconvolution Approach for the Mapping of Acoustic Sources (DAMAS) was developed which decoupled the array design and processing influence from the noise being measured, using a simple and robust algorithm. In 2005, three-dimensional applications of DAMAS were examined. DAMAS has been shown to render an unambiguous quantitative determination of acoustic source position and strength. However, an underlying premise of DAMAS, as well as that of classical array beamforming methodology, is that the noise regions under study are distributions of statistically independent sources. The present development, called DAMAS-C, extends the basic approach to include coherence definition between noise sources. The solutions incorporate cross-beamforming array measurements over the survey region. While the resulting inverse problem can be large and the iteration solution computationally demanding, it solves problems no other technique can approach. DAMAS-C is validated using noise source simulations and is applied to airframe flap noise test results.

Brooks, Thomas F.↗

Tracer-field cross-correlations with k -nearest neighbour distributions

ABSTRACT In astronomy and cosmology significant effort is devoted to characterizing and understanding spatial cross-correlations between points – e.g galaxy positions, high energy neutrino arrival directions, X-ray and AGN sources, and continuous fields – e.g. weak lensing meiand Cosmic Microwave Background maps. Recently, we introduced the k-nearest neighbour (kNN) formalism to better characterize the clustering of discrete (point) data sets. Here, we extend it to the point – field cross-correlations analysis. It combines kNN measurements of the point data set with measurements of the field smoothed at many scales. The resulting statistics are sensitive to all orders in the joint clustering of the points and the field. We demonstrate that this approach, unlike the 2-pt cross-correlation, can measure the statistical dependence of two data sets even when there are no linear (Gaussian) correlations between them. We further demonstrate that this framework is far more effective than the two point function in detecting cross-correlations when the continuous field is contaminated by high levels of noise. For a particularly high level of noise, the cross-correlation between haloes and the underlying matter field in a cosmological simulation, between 10 h−1 Mpc and 30 h−1 Mpc, is detected at >5σ significance using the technique presented here, when the two-point cross-correlation significance is ∼1σ. Finally, we show that kNN cross-correlations of haloes and the matter field can be well modelled on quasi-linear scales using the Hybrid Effective Field Theory (HEFT) framework, with the same set of bias parameters as are used for 2-pt cross-correlations. The substantial improvement in the statistical power of detecting cross-correlations using this method makes it a promising tool for various cosmological applications.

79 ASTRONOMY AND ASTROPHYSICS↗

Data from: 'Abiotic influences on continuous conifer forest structure across a subalpine watershed'

This package archives the core data used for analysis and inference in 'Abiotic influences on continuous conifer forest structure across a subalpine watershed' (Worsham et al., 2025). All data were collected in the East River, Washington Gulch, Slate River, and Coal Creek watersheds of Colorado. In the paper, we quantified the relative influence of climate, topographic, edaphic, and geologic factors on conifer stand structure and composition, and their functional relationships, at the watershed scale. We used waveform LiDAR data to derive spatially continuous stand structure metrics. We fused these with a species-level classification map to estimate tree species abundance. We applied generalized additive and generalized boosted models to evaluate the covariability of structural and compositional metrics with abiotic variables. The package contains the essential products required for reproducing our analysis and the tables and figures reported in the publication. The products comprise four classes: (1) geospatial data, (2) tabular data used for inferential analysis, (3) tabular data describing analytical results and performance statistics, and (4) a data user guide. (1) includes discretized waveform LiDAR data, locations and attributes of individual tree crowns, sampling locations and domain boundaries, a canopy height model, and raster files of estimated forest structural and compositional metrics at 100 m grid scale. (2) includes all response and explanatory variable values applied in inferential models. Response variables include conifer forest stand density, basal area, 95th percentile height, quadratic mean diameter, and others. Explanatory variables include climatic water deficit, actual evapotranspiration, elevation, heat load, soil available water content, and others. (3) includes results of training and testing several individual tree detection (ITD) algorithms, as well as inferential modeling results. (4) is a PDF user guide for this data package, including detailed descriptions and data dictionaries for all files. The data package root contains 17 assets: 8 compressed tape archive (.tar.gz) files, 5 comma-separated values (.csv) files, 3 Geographic Tagged Image File Format (GeoTIFF) (.tif) files, and 1 Portable Document Format (.pdf) file. The compressed .tar.gz archives contain ESRI shapefiles (.shp) .tif, compressed LASer (.laz), and .csv files. The archives must first be decompressed using the widely distributed command-line software utility TAR. All other files, including constituent files within the .tar.gz archives, can be opened in the open-source R statistical computing environment. Alternatively, .csv files may also be read in any simple text editor software or Microsoft Excel. Geospatial files including .shp and .tif files can also be opened in GIS software, such as QGIS (open-source) or ESRI ArcGIS (proprietary). The .pdf Data User Guide can be read with Adobe Acrobat Reader or other compatible readers.

2018 NEON and 2025 CHESS Campaigns↗

Reflectance and Thermal Infrared Spectroscopy of Mars: Relationship Between ISM and TES for Compositional Determinations

Reflectance spectroscopy has demonstrated that high albedo surfaces on Mars contain heavily altered materials with some component of hematite, poorly crystalline ferric oxides, and an undefined silicate matrix. The spectral properties of many low albedo regions indicate crystalline basalts containing both low and high calcium pyroxene, a mineralogy consistent with the basaltic SNC meteorites. The Thermal Emission Spectrometer (TES) experiment on the Mars Geochemical Surveyor has acquired critical new data relevant to surface composition and mineralogy, but in a wavelength region that is complementary to reflectance spectroscopy. The essence of the completed research was to analyze TES data in the context of reflectance data obtained by the French ISM imaging spectrometer experiment in 1989. This approach increased our understanding of the complementary nature of these wavelength regions for mineralogic determinations using actual observations of the martian surface. The research effort focused on three regions of scientific importance: Syrtis Major-Isidis Basin, Oxia Palus-Arabia, and Valles Marineris. In each region distinct spatial variations related to reflectance, and in derived mineralogic information and interpreted compositional units were analyzed. In addition, specific science questions related to the composition of volcanics and crustal evolution, soil compositions and pedogenic processes, and the relationship between pristine lithologies and weathering provided an overall science-driven framework for the work. The detailed work plan involved colocation of TES and ISM data, extraction of reflectance and emissivity spectra from areas of known reflectance variability, and quantitative analysis using factor analysis and statistical techniques to determine the degree of correspondence between these different wavelength regions. Identified coherent variations in TES spectroscopy were assessed against known atmospheric effects to validate that the variations are due to surface properties. With this new understanding of reflectance and emission spectroscopy, mineralogic interpretations were derived and applied to the science objectives of the three regions.

Boyce, Joseph↗

Modification of Eccentric Gaze-Holding

Clear vision and accurate localization of objects in the environment are prerequisites for reliable performance of motor tasks. Space flight confronts the crewmember with a stimulus rearrangement that requires adaptation to function effectively with the new requirements of altered spatial orientation and motor coordination. Adaptation and motor learning driven by the effects of cerebellar disorders may share some of the same demands that face our astronauts. One measure of spatial localization shared by the astronauts and those suffering from cerebellar disorders that is easily quantified, and for which a neurobiological substrate has been identified, is the control of the angle of gaze (the "line of sight"). The disturbances of gaze control that have been documented to occur in astronauts and cosmonauts, both in-flight and postflight, can be directly related to changes in the extrinsic gravitational environment and intrinsic proprioceptive mechanisms thus, lending themselves to description by simple non-linear statistical models. Because of the necessity of developing robust normal response populations and normative populations against which abnormal responses can be evaluated, the basic models can be formulated using normal, non-astronaut test subjects and subsequently extended using centrifugation techniques to alter the gravitational and proprioceptive environment of these subjects. Further tests and extensions of the models can be made by studying abnormalities of gaze control in patients with cerebellar disease. A series of investigations were conducted in which a total of 62 subjects were tested to: (1) Define eccentric gaze-holding parameters in a normative population, and (2) explore the effects of linear acceleration on gaze-holding parameters. For these studies gaze-holding was evaluated with the subjects seated upright (the normative values), rolled 45 degrees to both the left and right, or pitched back 30 and 90 degrees. In a separate study the further effects of acceleration on gaze stability was examined during centrifugation (+2 G (sub x) and +2 G (sub z) using a total of 23 subjects. In all of our investigations eccentric gaze-holding was established by having the subjects acquire an eccentric target (+/-30 degrees horizontal, +/- 15 degrees vertical) that was flashed for 750 msec in an otherwise dark room. Subjects were instructed to hold gaze on the remembered position of the flashed target for 20 sec. Immediately following the 20 sec period, subjects were cued to return to the remembered center position and to hold gaze there for an additional 20 sec. Following this 20 sec period the center target was briefly flashed and the subject made any corrective eye movement back to the true center position. Conventionally, the ability to hold eccentric gaze is estimated by fitting the natural log of centripetal eye drifts by linear regression and calculating the time constant (G) of these slow phases of "gaze-evoked nystagmus". However, because our normative subjects sometimes showed essentially no drift (tau (sub c) = m), statistical estimation and inference on the effect of target direction was performed on values of the decay constant theta = 1/(tau (sub c)) which we found was well modeled by a gamma distribution. Subjects showed substantial variance of their eye drifts, which were centrifugal in approximately 20 % of cases, and > 40% for down gaze. Using the ensuing estimated gamma distributions, we were able to conclude that rightward and leftward gaze holding were not significantly different, but that upward gaze holding was significantly worse than downward (p<0.05). We also concluded that vertical gaze holding was significantly worse than horizontal (p<0.05). In the case of left and right roll, we found that both had a similar improvement to horizontal gaze holding (p<0.05), but didn't have a significant effect on vertical gaze holding. For pitch tilts, both tilt angles significantly decreased gaze-holding ility in all directions (p<0.05). Finally, we found that hyper-g centrifugation significantly decreased gaze holding ability in the vertical plane. The main findings of this study are as follows: (1) vertical gaze-holding is less stable than horizontal, (2) gaze-holding to upward targets is less stable than to downward targets, (3) tilt affects gaze holding, and (4) hyper-g affects gaze holding. This difference between horizontal and vertical gaze-holding may be ascribed to separate components of the velocity-to-position neural integrator for eye movements, and to differences in orbital mechanics. The differences between upward and downward gaze-holding may be ascribed to an inherent vertical imbalance in the vestibular system. Because whole body tilt and hyper-g affects gaze-holding, it is implied that the otolith organs have direct connections to the neural integrator and further studies of astronaut gaze-holding are warranted. Our statistical method for representing the range of normal eccentric gaze stability can be readily applied to normals who maybe exposed to environments which may modify the central integrator and require monitoring, and to evaluate patients with gaze-evoked nystagmus by comparing to the above established normative criteria.

Reschke, M. F.↗

A Tracking Sun Photometer Without Moving Parts

This innovation is small, lightweight, and consumes very little electricity as it measures the solar energy attenuated by gases and aerosol particles in the atmosphere. A Sun photometer is commonly used on the Earth's surface, as well as on aircraft, to determine the solar energy attenuated by aerosol particles in the atmosphere and their distribution of sizes. This information is used to determine the spatial and temporal distribution of gases and aerosols in the atmosphere, as well as their distribution sizes. The design for this Sun photometer uses a combination of unique optics and a charge coupled device (CCD) array to eliminate moving parts and make the instrument more reliable. It could be selfcalibrating throughout the year. Data products would be down-welling flux, the direct-diffuse flux ratio, column abundance of gas phase constituents, aerosol optical depth at multiple-wavelengths, phase functions, cloud statistics, and an estimate of the representative size of atmospheric particles. These measurements can be used to obtain an estimate of aerosol size distribution, refractive index, and particle shape. Incident light is received at a light-reflecting (inner) surface, which is a truncated paraboloid. Light arriving from a hemispheric field of view (solid angle 2 steradians) enters the reflecting optic at an entrance aperture at, or adjacent to, the focus of the paraboloid, and is captured by the optic. Most of this light is reflected from an inner surface. The light proceeds substantially parallel to the paraboloid axis, and is detected by an array detector located near an exit aperture. Each of the entrance and exit apertures is formed by the intersection of the paraboloid with a plane substantially perpendicular to the paraboloid axis. Incident (non-reflected) light from a source of limited extent (the Sun) illuminates a limited area on the detector array. Both direct and diffuse illumination may be reflected, or not reflected, before being received on the detector array. As the Sun traverses a path in the sky over some time interval, the track of the Sun can be traced on the detector array. A suitably modified Sun photometer might be used to study the dynamics of an environment on another planet or satellite with an atmosphere.

Strawa, Anthony W.↗

Metric Learning to Enhance Hyperspectral Image Segmentation

Unsupervised hyperspectral image segmentation can reveal spatial trends that show the physical structure of the scene to an analyst. They highlight borders and reveal areas of homogeneity and change. Segmentations are independently helpful for object recognition, and assist with automated production of symbolic maps. Additionally, a good segmentation can dramatically reduce the number of effective spectra in an image, enabling analyses that would otherwise be computationally prohibitive. Specifically, using an over-segmentation of the image instead of individual pixels can reduce noise and potentially improve the results of statistical post-analysis. In this innovation, a metric learning approach is presented to improve the performance of unsupervised hyperspectral image segmentation. The prototype demonstrations attempt a superpixel segmentation in which the image is conservatively over-segmented; that is, the single surface features may be split into multiple segments, but each individual segment, or superpixel, is ensured to have homogenous mineralogy.

Thompson, David R.↗

Characterization of Extreme Wave Conditions for Wave Energy Converter Design and Project Risk Assessment

Best practices and international standards for determining n-year return period extreme wave (sea states) conditions allow wave energy converter designers and project developers the option to apply simple univariate or more complex bivariate extreme value analysis methods. The present study compares extreme sea state estimates derived from univariate and bivariate methods and investigates the performance of spectral wave models for predicting extreme sea states at buoy locations within several regional wave climates along the US East and West Coasts. Two common third-generation spectral wave models are evaluated, a WAVEWATCH III® model with a grid resolution of 4 arc-minutes (6–7 km), and a Simulating WAves Nearshore model, with a coastal resolution of 200–300 m. Both models are used to generate multi-year hindcasts, from which extreme sea state statistics used for wave conditions characterization can be derived and compared to those based on in-situ observations at National Data Buoy Center stations. Comparison of results using different univariate and bivariate methods from the same data source indicates reasonable agreement on average. Discrepancies are predominantly random. Large discrepancies are common and increase with return period. There is a systematic underbias for extreme significant wave heights derived from model hindcasts compared to those derived from buoy measurements. This underbias is dependent on model spatial resolution. However, simple linear corrections can effectively compensate for this bias. A similar approach is not possible for correcting model-derived environmental contours, but other methods, e.g., machine learning, should be explored.

16 TIDAL AND WAVE POWER↗

Robust Measurement of Stellar Streams around the Milky Way: Correcting Spatially Variable Observational Selection Effects in Optical Imaging Surveys

Observations of density variations in stellar streams are a promising probe of low-mass dark matter substructure in the Milky Way. However, survey systematics such as variations in seeing and sky brightness can also induce artificial fluctuations in the observed densities of known stellar streams. These variations arise because survey conditions affect both object detection and star–galaxy misclassification rates. To mitigate these effects, we use Balrog synthetic source injections in the Dark Energy Survey (DES) Y3 data to calculate detection rate variations and classification rates as functions of survey properties. We show that these rates are nearly separable with respect to survey properties and can be estimated with sufficient statistics from the synthetic catalogs. Applying these corrections reduces the standard deviation of relative detection rates across the DES footprint by a factor of 5, and our corrections significantly change the inferred linear density of the Phoenix stream when including faint objects. Additionally, for artificial streams with DES-like survey properties we are able to recover density power spectra with reduced bias. We also find that uncorrected power-spectrum results for Legacy Survey of Space and Time (LSST)-like data can be around 5 times more biased, highlighting the need for such corrections in future ground-based surveys.

79 ASTRONOMY AND ASTROPHYSICS↗

Measurement of high-degree solar oscillation frequencies

We present m-averaged solar p- and f-mode oscillation frequencies over the frequency range nu greater than 1.8 and less than 5.0 mHz and the spherical harmonic degree range l greater than or equal to 100 and less than or equal to 1200 from full-disk, 1000 x 1024 pixel, Ca II intensity images collected 1993 June 22-25 with a temporal cadence of 60 s. We itemize the sources and magnitudes of statistical and systematic uncertainties and of small frequency corrections, and we show that our frequencies represent an improvement in accuracy and coverage over previous measurements. Our frequencies agree at the 2 micro Hz level with Mount Wilson frequencies determined for l less than or equal to 600 from full-disk images, and we find systematic offsets of 10-20 micro Hz with respect to frequencies measured from Big Bear and La Palma observations. We give evidence that these latter offsets are indicative of spatial scaling uncertainties associated with the analysis of partial-disk images. In comparison with theory, our p-mode frequencies agree within 10 micro Hz of frequencies predicted by the Los Alamos model but are as much as 100 micro Hz smaller than frequencies predicted by the Denmark and Yale models at degrees near 1000. We also find systematic differences between our n = 0 frequencies and the frequencies closely agreed upon by all three models.

Bachmann, K. T.↗

The EPIC-MOS Particle-Induced Background Spectra

In order to analyse diffuse emission that fills the field of view, one must accurately characterize the instrumental backgrounds. For the XMM-Newton EPIC instrument these backgrounds include a temporally variable "quiescent" component. as well as the strongly variable soft proton contamination. We have characterized the spectral and spatial response of the EPIC detectors to these background components and have developed tools to remove these backgrounds from observations. The "quiescent" component was characterized using a combination of the filter-wheel-closed data and a database of unexposed-region data. The soft proton contamination was characterized by differencing images and spectra taken during flared and flare-free intervals. After application of our modeled backgrounds, the differences between independent observations of the same region of "blank sky" are consistent with the statistical uncertainties except when there is clear spectral evidence of solar wind charge exchange emission. Using a large sample of blank sky data, we show that strong magnetospheric SWCX emission requires elevated solar wind fluxes; observations through the densest part of the magnetosheath are not necessarily strongly contaminated with SWCX emission.

Kuntz, K. D.↗

Comparison of Surface and Column Variations of CO2 Over Urban Areas for Future Active Remote CO2 Sensors

High resolution in-situ CO2 measurements were recorded onboard the NASA P-3B during the DISCOVER-AQ (Deriving Information on Surface Conditions from Column and Vertically Resolved Observations Relevant to Air Quality) Field Campaign, to investigate the ability of space-based observations to accurately assess near surface conditions related to air quality. This campaign includes, Washington DC/Baltimore, MD (July 2011), San Joaquin Valley, CA (January - February 2013), Houston, TX (September 2013), and Denver, CO (July-August 2014). Each of these campaigns consisted of missed approaches and approximately two hundred vertical soundings of CO2 within the lower troposphere (surface to about 5 km). In this study, surface (0 - 1 km) and column-averaged (0 - 3.5 km) CO2 mixing ratio values from the vertical soundings in the four geographically different urban areas are used to investigate the temporal and spatial variability of CO2 within the different urban atmospheric emission environments. Tracers such as CO, CH2O, NOx, and NMHCs are used to identify the source of CO2 variations in the urban sites. Additionally, we apply nominal CO2 column weighting functions for potential future active remote CO2 sensors operating in the 1.57-microns and 2.05-microns measurement regions to convert the in situ CO2 vertical mixing ratio profiles to variations in CO2 column optical depths, which is what the active remote sensors actually measure. Using statistics calculated from the optical depths at each urban site measured during the DISCOVER-AQ field campaign and for each nominal weighting function, we investigate the natural variability of CO2 columns in the lower troposphere; relate the CO2 column variability to the urban surface emissions; and show the measurement requirements for the future ASCENDS (Active Sensing of CO2 Emissions over Nights, Days, and Seasons) in the continental U.S. urban areas.

Choi, Yonghoon↗

Global Gridded Crop Model Evaluation: Benchmarking, Skills, Deficiencies and Implications.

Crop models are increasingly used to simulate crop yields at the global scale, but so far there is no general framework on how to assess model performance. Here we evaluate the simulation results of 14 global gridded crop modeling groups that have contributed historic crop yield simulations for maize, wheat, rice and soybean to the Global Gridded Crop Model Intercomparison (GGCMI) of the Agricultural Model Intercomparison and Improvement Project (AgMIP). Simulation results are compared to reference data at global, national and grid cell scales and we evaluate model performance with respect to time series correlation, spatial correlation and mean bias. We find that global gridded crop models (GGCMs) show mixed skill in reproducing time series correlations or spatial patterns at the different spatial scales. Generally, maize, wheat and soybean simulations of many GGCMs are capable of reproducing larger parts of observed temporal variability (time series correlation coefficients (r) of up to 0.888 for maize, 0.673 for wheat and 0.643 for soybean at the global scale) but rice yield variability cannot be well reproduced by most models. Yield variability can be well reproduced for most major producing countries by many GGCMs and for all countries by at least some. A comparison with gridded yield data and a statistical analysis of the effects of weather variability on yield variability shows that the ensemble of GGCMs can explain more of the yield variability than an ensemble of regression models for maize and soybean, but not for wheat and rice. We identify future research needs in global gridded crop modeling and for all individual crop modeling groups. In the absence of a purely observation-based benchmark for model evaluation, we propose that the best performing crop model per crop and region establishes the benchmark for all others, and modelers are encouraged to investigate how crop model performance can be increased. We make our evaluation system accessible to all crop modelers so that other modeling groups can also test their model performance against the reference data and the GGCMI benchmark.

wheat↗