Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Statistical accuracy”

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 865 records · Page 48

Toward the characterization of upper tropospheric clouds using Atmospheric Infrared Sounder and Microwave Limb Sounder observations

We estimate the accuracy of cloud top altitude (Z) retrievals from the Atmospheric Infrared Sounder (AIRS) and Advanced Microwave Sounding Unit (AMSU) observing suite (ZA) on board the Earth Observing System Aqua platform. We compare ZA with coincident measurements of Z derived from the micropulse lidar and millimeter wave cloud radar at the Atmospheric Radiation Measurement (ARM) program sites of Nauru and Manus islands (ZARM) and the inferred Z from vertically resolved Microwave Limb Sounder (MLS) ice water content (IWC) retrievals. The mean difference in ZA minus ZARM plus or minus one standard deviation ranges from −2.2 to 1.6 km ± 1.0 to 4.2 km for all cases of AIRS effective cloud fraction (f A) > 0.15 at Manus Island using the cloud radar only. The range of mean values results from using different approaches to determine ZARM, day/night differences, and the magnitude of f A; the variation about the mean decreases for increasing values of f A. Analysis of ZARM from the micropulse lidar at Nauru Island for cases restricted to 0.05 ≤ f A ≤ 0.15 indicates a statistically significant improvement in ZA − ZARM over the cloud radar-derived values at Manus Island. In these cases the ZA − ZARM difference is −1.1 to 2.1 km ± 3.0 to 4.5 km. These results imply that the operational ZA is quantitatively useful for constraining cirrus altitude despite the nominal 45 km horizontal resolution. Mean differences of cloud top pressure (PCLD) inferred from coincident AIRS and MLS ice water content (IWC) retrievals depend upon the method of defining AIRS PCLD (as with the ARM comparisons) over the MLS spatial scale, the peak altitude and maximum value of MLS IWC, and f A. AIRS and MLS yield similar vertical frequency distributions when comparisons are limited to f A > 0.1 and IWC > 1.0 mg m−3. Therefore the agreement depends upon the opacity of the cloud, with decreased agreement for optically tenuous clouds. Further, the mean difference and standard deviation of AIRS and MLS PCLD are highly dependent on the MLS tangent altitude. For MLS tangent altitudes greater than 146 hPa, the strength of the limb technique, the disagreement becomes statistically significant. This implies that AIRS and MLS “agree” in a statistical sense at lower tangent altitudes and “disagree” at higher tangent altitudes. These results provide important insights on upper tropospheric cloudiness as observed by nadir-viewing AIRS and limb-viewing MLS.

Wu, Dong L.↗

Testing higher-order Lagrangian perturbation theory against numerical simulations. 2: Hierarchical models

We present results showing an improvement of the accuracy of perturbation theory as applied to cosmological structure formation for a useful range of scales. The Lagrangian theory of gravitational instability of Friedmann-Lemaitre cosmogonies is compared with numerical simulations. We study the dynamics of hierarchical models as a second step. In the first step we analyzed the performance of the Lagrangian schemes for pancake models, the difference being that in the latter models the initial power spectrum is truncated. This work probed the quasi-linear and weakly non-linear regimes. We here explore whether the results found for pancake models carry over to hierarchical models which are evolved deeply into the non-linear regime. We smooth the initial data by using a variety of filter types and filter scales in order to determine the optimal performance of the analytical models, as has been done for the 'Zel'dovich-approximation' - hereafter TZA - in previous work. We find that for spectra with negative power-index the second-order scheme performs considerably better than TZA in terms of statistics which probe the dynamics, and slightly better in terms of low-order statistics like the power-spectrum. However, in contrast to the results found for pancake models, where the higher-order schemes get worse than TZA at late non-linear stages and on small scales, we here find that the second-order model is as robust as TZA, retaining the improvement at later stages and on smaller scales. In view of these results we expect that the second-order truncated Lagrangian model is especially useful for the modelling of standard dark matter models such as Hot-, Cold-, and Mixed-Dark-Matter.

Melott, A. L.↗

Boiler Tube Corrosion Characterization with a Scanning Thermal Line

Wall thinning due to corrosion in utility boiler water wall tubing is a significant operational concern for boiler operators. Historically, conventional ultrasonics has been used for inspection of these tubes. Unfortunately, ultrasonic inspection is very manpower intense and slow. Therefore, thickness measurements are typically taken over a relatively small percentage of the total boiler wall and statistical analysis is used to determine the overall condition of the boiler tubing. Other inspection techniques, such as electromagnetic acoustic transducer (EMAT), have recently been evaluated, however they provide only a qualitative evaluation - identifying areas or spots where corrosion has significantly reduced the wall thickness. NASA Langley Research Center, in cooperation with ThermTech Services, has developed a thermal NDE technique designed to quantitatively measure the wall thickness and thus determine the amount of material thinning present in steel boiler tubing. The technique involves the movement of a thermal line source across the outer surface of the tubing followed by an infrared imager at a fixed distance behind the line source. Quantitative images of the material loss due to corrosion are reconstructed from measurements of the induced surface temperature variations. This paper will present a discussion of the development of the thermal imaging system as well as the techniques used to reconstruct images of flaws. The application of the thermal line source coupled with the analysis technique represents a significant improvement in the inspection speed and accuracy for large structures such as boiler water walls. A theoretical basis for the technique will be presented to establish the quantitative nature of the technique. Further, a dynamic calibration system will be presented for the technique that allows the extraction of thickness information from the temperature data. Additionally, the results of the application of this technology to actual water wall tubing samples and in-situ inspections will be presented.

Cramer, K. Elliott↗

An error covariance model for sea surface topography and velocity derived from TOPEX/POSEIDON altimetry

In order to facilitate the use of satellite-derived sea surface topography and velocity oceanographic models, methodology is presented for deriving the total error covariance and its geographic distribution from TOPEX/POSEIDON measurements. The model is formulated using a parametric model fit to the altimeter range observations. The topography and velocity modeled with spherical harmonic expansions whose coefficients are found through optimal adjustment to the altimeter range residuals using Bayesian statistics. All other parameters, including the orbit, geoid, surface models, and range corrections are provided as unadjusted parameters. The maximum likelihood estimates and errors are derived from the probability density function of the altimeter range residuals conditioned with a priori information. Estimates of model errors for the unadjusted parameters are obtained from the TOPEX/POSEIDON postlaunch verification results and the error covariances for the orbit and the geoid, except for the ocean tides. The error in the ocean tides is modeled, first, as the difference between two global tide models and, second, as the correction to the present tide model, the correction derived from the TOPEX/POSEIDON data. A formal error covariance propagation scheme is used to derive the total error. Our global total error estimate for the TOPEX/POSEIDON topography relative to the geoid for one 10-day period is found tio be 11 cm RMS. When the error in the geoid is removed, thereby providing an estimate of the time dependent error, the uncertainty in the topography is 3.5 cm root mean square (RMS). This level of accuracy is consistent with direct comparisons of TOPEX/POSEIDON altimeter heights with tide gauge measurements at 28 stations. In addition, the error correlation length scales are derived globally in both east-west and north-south directions, which should prove useful for data assimilation. The largest error correlation length scales are found in the tropics. Errors in the velocity field are smallest in midlatitude regions. For both variables the largest errors caused by uncertainty in the geoid. More accurate representations of the geoid await a dedicated geopotential satellite mission. Substantial improvements in the accuracy of ocean tide models are expected in the very near future from research with TOPEX/POSEIDON data.

Tsaoussi, Lucia S.↗

Hierarchical Statistical 3D ' Atomistic' Simulation of Decanano MOSFETs: Drift-Diffusion, Hydrodynamic and Quantum Mechanical Approaches

When MOSFETs are scaled to deep submicron dimensions the discreteness and randomness of the dopant charges in the channel region introduces significant fluctuations in the device characteristics. This effect, predicted 20 year ago, has been confirmed experimentally and in simulation studies. The impact of the fluctuations on the functionality, yield, and reliability of the corresponding systems shifts the paradigm of the numerical device simulation. It becomes insufficient to simulate only one device representing one macroscopical design in a continuous charge approximation. An ensemble of macroscopically identical but microscopically different devices has to be characterized by simulation of statistically significant samples. The aims of the numerical simulations shift from predicting the characteristics of a single device with continuous doping towards estimating the mean values and the standard deviations of basic design parameters such as threshold voltage, subthreshold slope, transconductance, drive current, etc. for the whole ensemble of 'atomistically' different devices in the system. It has to be pointed out that even the mean values obtained from 'atomistic' simulations are not identical to the values obtained from continuous doping simulations. In this paper we present a hierarchical approach to the 'atomistic' simulation of aggressively scaled decanano MOSFETs. A full scale 3D drift-diffusion'atomostic' simulation approach is first described and used for verification of the more economical, but also more restricted, options. To reduce the processor time and memory requirements at high drain voltage we have developed a self-consistent option based on a thin slab solution of the current continuity equation only in the channel region. This is coupled to the Poisson's equation solution in the whole simulation domain in the Gummel iteration cycles. The accuracy of this approach is investigated in comparison with the full self-consistent solution. At low drain voltage only single solution of the nonlinear Poisson equation is sufficient to extract the current with satisfactory accuracy. A pilot version of a hydrodynamic 'atomistic' simulator has been developed in order to study the effect of the nonequilibrium, non local transport in decanano MOSFETs on the random dopant induced current fluctuations. For the first time we have also applied the density gradient approach in 3D to investigate the effect of the quantum confinement on the threshold voltage fluctuations. The developed 'atomistic' simulation techniques have been applied to study various fluctuation resistant MOSFET architectures including epitaxial and delta doped devices.

Asenov, Asen↗

Land-surface temperature measurement from space - Physical principles and inverse modeling

To apply the multiple-wavelength (split-window) method used for satellite measurement of sea-surface temperature from thermal-infrared data to land-surface temperatures, the authors statistically analyze simulations using an atmospheric radiative transfer model. The range of atmospheric conditions and surface temperatures simulated is wide enough to cover variations in clear atmospheric properties and surface temperatures, both of which are larger over land than over sea. Surface elevation is also included in the simulation as the most important topographic effect. Land covers characterized by measured or modeled spectral emissivities include snow, clay, sands, and tree leaf samples. The empirical inverse model can estimate the surface temperature with a standard deviation less than 0.3 K and a maximum error less than 1 K, for viewing angles up to 40 degrees from nadir under cloud-free conditions, given satellite measurements in three infrared channels. A band in the region from 10.2 to 11.0 microns will usually give the most reliable single-band estimate of surface temperature. In addition, a band in either the 3.5-4.0-micron region or in the 11.5-12.6-micron region must be included for accurate atmospheric correction, and a band below the ozone absorption feature at 9.6 microns (e.g., 8.2-8.8 microns) will increase the accuracy of the estimate of surface temperature.

Wan, Zhengming↗

Calibration of the 11x11-Foot Transonic Wind Tunnel at NASA Ames Research Center

A static pipe calibration characterizing the streamwise static pressure distribution was conducted at the 11-By 11-Foot Transonic Wind Tunnel at NASA Ames Research Center. This data is used to determine the local Mach number within the test section and evaluate buoyancy corrections to axial force measurements. The 60 foot long, 6 inch diameter pipe contained 444 static pressure taps spanning the test section and nozzle regions of the tunnel. The forward end of the pipe extends into the settling chamber and is held by four cables mounting to the tunnel shell, and the aft end is fixed on the institutional model support system. A hydraulic cylinder at the aft end of the pipe provides tension on the system to reduce vibration and to keep the pipe as level as possible throughout the test section. The previous calibration was improved upon by using pressure scanners with greater accuracy, ensuring a uniform pressure tube length for each tap to control pneumatic lag, optically tracking any streamwise movement of the pipe, and more tightly controlling the tunnel condition set points. Typically this calibration is conducted with the pipe on tunnel centerline and 33 inches below centerline for sting-mounted models and semi-span (i.e. floor-mounted) models respectively, however schedule demands permitted only the centerline calibration. The semi-span calibration is planned to be completed in the summer of 2020. Immediately following the static pipe calibration, a shorter, 9 foot static pipe used as the calibration check standard was installed to obtain its first post-calibration pressure dataset. This short static pipe consists of 148 static pressure taps distributed along the pipe section and one total pressure tap at the end of an ogive nose. Performing the calibration test and its check standard back-to-back allows this dataset to establish a reliable baseline for future calibration check standard testing. Over time the use of a calibration check standard offers the ability to assess the stability of the calibration through statistical process control in an efficient and costeffective manner, thereby potentially increasing the time required between full tunnel calibrations.

White, Bethany↗

A Dark Target research aerosol algorithm for MODIS observations over eastern China: increasing coverage while maintaining accuracy at high aerosol loading

Satellite aerosol products such as the Dark Target (DT) produced from the MODerate resolution Imaging Spectroradiometer (MODIS) are useful for monitoring the progress of air pollution. Unfortunately, the DT often fails to retrieve during the heaviest aerosol events as well as the more moderate events in winter. Some of the literature at-tributes this lack of retrieval to the cloud mask. However, we found this lack of retrieval is mainly traced to thresholds used for masking of inland water and snow. Modifications to these two masks greatly increase 50 % of the retrievals of aerosol optical depth at 0.55 μm (AOD) greater than 1.0. The “extra”-high-AOD retrievals tend to be biased when com-pared with a ground-based sun photometer (AErosol RObotic NETwork, AERONET). Reducing bias in new retrievals re-quires two additional steps. One is an update to the assumed aerosol optical properties (aerosol model); the haze in this region is both less absorbing and lower in altitude than what is assumed in the global algorithm. The second is account-ing for the scale height of the aerosol, specifically that the heavy-aerosol events in the region are much closer to the surface than what is assumed by the global DT algorithm. The resulting combination of modified masking thresholds, new aerosol model, and lower aerosol layer scale height was applied to 3 months of MODIS observations (January–March2013) over eastern China. After these two additional steps are implemented, the significant increase in new retrievals introduces no overall bias at a high-AOD regime but does degrade other overall validation statistics. We also find that the research algorithm is able to identify additional pollution events that AERONET instruments may not due to different spatial sampling. Mean AOD retrieved from the re-search algorithm increases from 0.11 to 0.18 compared to values calculated from the operational DT algorithm during January to March of 2013 over the study area. But near Beijing, where the severe pollution occurs, the new algorithm increases AOD by as much as 3.0 for each 0.5°grid box over the previous operational-algorithm values.

Dark Target↗

Global Energetics of Solar Flares. VII. Aerodynamic Drag in Coronal Mass Ejections

The free energy that is dissipated in a magnetic reconnection process of a solar flare, generally accompanied by a coronal mass ejection (CME), has been considered as the ultimate energy source of the global energy budget of solar flares in previous statistical studies. Here we explore the effects of the aerodynamic drag force on CMEs, which supplies additional energy from the slow solar wind to a CME event, besides the magnetic energy supply. For this purpose, we fit the analytical aerodynamic drag model of Cargill and Vršnak et al. to the height–time profiles r(t) of LASCO/SOHO data in 14,316 CME events observed during the first 8 yr (2010–2017) of the Solar Dynamics Observatory era (ensuring EUV coverage with AIA). Our main findings are (1) a mean solar wind speed of w = 472 ± 414 km s(exp −1), (2) a maximum drag-accelerated CME energy of E(drag) <~2 × 10(exp32) erg, (3) a maximum flare-accelerated CME energy of E(flare <~1.5 × 10(exp33) erg, (4) the ratio of the summed kinetic energies of all flare accelerated CMEs to the drag-accelerated CMEs amounts to a factor of 4, (5) the inclusion of the drag force slightly lowers the overall energy budget of CME kinetic energies in flares from ≈7% to ≈4%, and (6) the arrival times of CMEs at Earth can be predicted with an accuracy of ≈23%.

Aschwanden, Markus J.↗

Ionospheric corrections to precise time transfer using GPS

The free electrons in the earth's ionosphere can retard the time of reception of GPS signals received at a ground station, compared to their time in free space, by many tens of nanoseconds, thus limiting the accuracy of time transfer by GPS. The amount of the ionospheric time delay is proportional to the total number of electrons encountered by the wave on its path from each GPS satellite to a receiver. This integrated number of electrons is called Total Electron Content, or TEC. Dual frequency GPS receivers designed by Allen Osborne Associates, Inc. (AOA) directly measure both the ionospheric differential group delay and the differential carrier phase advance for the two GPS frequencies and derive from this the TEC between the receiver and each GPS satellite in track. The group delay information is mainly used to provide an absolute calibration to the relative differential carrier phase, which is an extremely precise measure of relative TEC. The AOA Mini-Rogue ICS-4Z and the AOA TurboRogue ICS-4000Z receivers normally operate using the GPS P code, when available, and switch to cross-correlation signal processing when the GPS satellites are in the Anti-Spoofing (A-S) mode and the P code is encrypted. An AOA ICS-Z receiver has been operated continuously for over a year at Hanscom AFB, MA to determine the statistics of the variability of the TEC parameter using signals from up to four different directions simultaneously. The 4-channel ICS-4Z and the 8-channel ICS-4000Z, have proven capabilities to make precise, well calibrated, measurements of the ionosphere in several directions simultaneously. In addition to providing ionospheric corrections for precise time transfer via satellite, this dual frequency design allows full code and automatic codeless operation of both the differential group delay and differential carrier phase for numerous ionospheric experiments being conducted. Statistical results of the data collected from the ICS-4Z during the initial year of ionospheric time delay in the northeastern U.S., and initial results with the ICS-4000Z, will be presented.

Snow, Robert W.↗

Star Confusion Effect on SIM PlanetQuest Astrometric Performance

SIM PlanetQuest will measure star positions to an accuracy of a few microarcseconds using precise white light fringe measurements. One challenge for SIM observation scenario is "star confusion," where multiple stars are present in the instrument field of view. This is especially relevant for observing dim science targets because the density of number of stars increases rapidly with star magnitude. We study the effect of star confusion on the SIM astrometric performance due to systematic fringe errors caused by the extra photons from the confusion star(s}. Since star confusion from multiple stars may be analyzed as a linear superposition of the effect from single star confusion, we quantify the astrometric errors due to single star confusion surveying over many spectral types, including AOV, FOV, K5III, and MOV, and for various visual magnitude differences. To the leading order, the star confusion effect is characterized by the magnitude difference, spectral difference, and the angular separation between the target and confusion stars.Strategies for dealing with star confusion are presented. For example, since the presence of additional sources in the field of view leads to inconsistent delay estimates from different channels, with sufficient signal to noise ratio, the star confusion can be detected using chi-square statistics of fringe measurements from multiple spectral channels. An interesting result is that the star confusion can be detected even though the interferometer cannot resolve the separation between the target and confusion stars when their spectra are sufficiently different. Other strategies for mitigating the star confusion effect are also discussed.

crowded field↗

Quantitative Daily Maps of PM 2.5 Episodes for California and Other Regions: Satellite Column Water and Optical Depth as Allied Tracers of Dilution

The Western US and many regions globally present daunting difficulties in understanding PM 2.5 episodes. We evaluate extensions of a method independent of modeled source-description and transport/transformation and using several satellite remote sensing products from imaging spectrometers. The San Joaquin Valley (SJV) especially suffers few-day episodes due to shallow mixing; PM 2.5 retrieval suffers low satellite AOT (Aerosol Optical Thickness) and bright surfaces.Nevertheless, we find residual errors in our maps of of typically 5-8 micrograms per cubic meter. Episodes in the Valley reaching 60-100 micrograms per cubic meter. These maps detail pollution from Interstate 5 at the scale of a few kilometers. The maps are based on NASA's MODerate resolution Imaging Spectrometer (MODIS) data at circa 1 kilometer as processed with the Multi-Angle Implementation of Atmospheric Correction. The Bay Area Air Quality Management District has requested that we test our methods in their challenging environment characterized by multiple sub-basins defined by complex topography. Our tests suggest that nearly similar precision may be expected for wintertime conditions with high PM 2.5 . We note difficulties when measured PM 2.5 is less than 8-10 micrograms per cubic meter, but good relative precision when PM 2.5 rises above 20; i.e. in episodes of concern for morbidity and mortality. Our method stresses physically meaningful functions of MODIS-MAIAC (Multi-Angle Implementation of Atmospheric Correction)-derived AOD (Aerosol Optical Depth) and total water vapor column. A mixed-effects statistical model exploiting existing station data works powerfully to allow us daily AOT-to-PM 2.5 relationships that allow a calibration of the map. In those cases where water vapor and particles have generally similar surface sources, using the ratio of AOT / Column_water can improve the daily calibrations so as to reach our quoted precision. We briefly present some cartoon idealizations that explain this success and also the likely reasons that our mixed effects model (or "daily calibration") works; also when it should not work. The combined satellite/mixed-effects model works best for wintertime San Joaquin Valley episodes, where the meteorology of particle and H2O(v) dilution is quite appropriate. We extended and tested the methodology (a) for the Bay Area wintertime situations and (b) for smoke plume events (e.g. the October 2017 fire events of the Sonoma area). Our SJV work was evaluated using NASA's DISCOVER-AQ (Deriving Information on Surface conditions from Column and Vertically Resolved Observations Relevant to Air Quality) airborne measurements, and by season- long measurements in Fresno. If the composition and size distribution of the aerosols can be assessed for the regions we describe, retrievals should have improved accuracy.

Chatfield, Robert B.↗

Accuracy Assessment of Global Internal-Tide Models Using Satellite Altimetry Models Using Satellite Altimetry

Altimeter measurements are corrected for several geophysical parameters in order to access ocean signals of interest, like mesoscale or sub-mesoscale variability. The ocean tide is one of the most critical corrections due to the amplitude of the tidal elevations and to the aliasing phenomena of high-frequency signals into the lower-frequency band, but the internal-tide signatures at the ocean surface are not yet corrected globally. Internal tides can have a signature of several centimeters at the surface with wavelengths of about 50–250 km for the first mode and even smaller scales for higher-order modes. The goals of the upcoming Surface Water Ocean Topography (SWOT) mission and other high-resolution ocean measurements make the correction of these small-scale signals a challenge, as the correction of all tidal variability becomes mandatory to access accurate measurements of other oceanic signals. In this context, several scientific teams are working on the development of new internal-tide models, taking advantage of the very long altimeter time series now available, which represent an unprecedented and valuable global ocean database. The internal-tide models presented here focus on the coherent internal-tide signal and they are of three types: empirical models based upon analysis of existing altimeter missions, an assimilative model and a three-dimensional hydrodynamic model. A detailed comparison and validation of these internal-tide models is proposed using existing satellite altimeter databases. The analysis focuses on the four main tidal constituents: M2, K1, O1 and S2. The validation process is based on a statistical analysis of multi-mission altimetry including Jason-2 and Cryosphere Satellite-2 data. The results show a significant altimeter variance reduction when using internal-tide corrections in all ocean regions where internal tides are generating or propagating. A complementary spectral analysis also gives some estimation of the performance of each model as a function of wavelength and some insight into the residual non-stationary part of internal tides in the different regions of interest. This work led to the implementation of a new internal-tide correction (ZARON'one) in the next geophysical data records version-F (GDR-F) standards.

Loren Carrere↗

Statistical and Hydrological Evaluation of TRMM-Based Multi-Satellite Precipitation Analysis over the Wangchu Basin of Bhutan: Are the Latest Satellite Precipitation Products 3B42V7 Ready for Use in Ungauged Basins?

The objective of this study is to quantitatively evaluate the successive Tropical Rainfall Measuring Mission (TRMM) Multi-satellite Precipitation Analysis (TMPA) products and further to explore the improvements and error propagation of the latest 3B42V7 algorithm relative to its predecessor 3B42V6 using the Coupled Routing and Excess Storage (CREST) hydrologic model in the mountainous Wangchu Basin of Bhutan. First, the comparison to a decade-long (2001-2010) daily rain gauge dataset reveals that: 1) 3B42V7 generally improves upon 3B42V6s underestimation both for the whole basin (bias from -41.15 to -8.38) and for a 0.250.25 grid cell with high-density gauges (bias from -40.25 to 0.04), though with modest enhancement of correlation coefficients (CC) (from 0.36 to 0.40 for basin-wide and from 0.37 to 0.41 for grid); and 2) 3B42V7 also improves its occurrence frequency across the rain intensity spectrum. Using the CREST model that has been calibrated with rain gauge inputs, the 3B42V6-based simulation shows limited hydrologic prediction NSCE skill (0.23 in daily scale and 0.25 in monthly scale) while 3B42V7 performs fairly well (0.66 in daily scale and 0.77 in monthly scale), a comparable skill score with the gauge rainfall simulations. After recalibrating the model with the respective TMPA data, significant improvements are observed for 3B42V6 across all categories, but not as much enhancement for the already well-performing 3B42V7 except for a reduction in bias (from -26.98 to -4.81). In summary, the latest 3B42V7 algorithm reveals a significant upgrade from 3B42V6 both in precipitation accuracy (i.e., correcting the underestimation) thus improving its potential hydrological utility. Forcing the model with 3B42V7 rainfall yields comparable skill scores with in-situ gauges even without recalibration of the hydrological model by the satellite precipitation, a compensating approach often used but not favored by the hydrology community, particularly in ungauged basins.

TRMM↗

Evaluation of MODIS and VIIRS Cloud-Gap-Filled Snow-Cover Products for Production of an Earth Science Data Record

MODerate resolution Imaging Spectroradiometer (MODIS) cryosphere products have been available since 2000 – following the 1999 launch of the Terra MODIS and the 2002 launch of the Aqua MODIS – and include global snow-cover extent (SCE) (swath, daily, and 8 d composites) at 500 m and ∼5 km spatial resolutions. These products are used extensively in hydrological modeling and climate studies. Reprocessing of the complete snow-cover data record, from Collection 5 (C5) to Collection 6 (C6) and Collection 6.1 (C6.1), has provided improvements in the MODIS product suite. Suomi National Polar-orbiting Partnership (S-NPP) Visible Infrared Imaging Radiometer Suite (VIIRS) Collection 1 (C1) snow-cover products at a 375 m spatial resolution have been available since 2011 and are currently being reprocessed for Collection 2 (C2). Both the MODIS C6.1 and the VIIRS C2 products will be available for download from the National Snow and Ice Data Center beginning in early 2020 with the complete time series available in 2020. To address the need for a cloud-reduced or cloud-free daily SCE product for both MODIS and VIIRS, a daily cloud-gap-filled (CGF) snow-cover algorithm was developed for MODIS C6.1 and VIIRS C2 processing. MOD10A1F (Terra) and MYD10A1F (Aqua) are daily, 500 m resolution CGF SCE map products from MODIS. VNP10A1F is the daily, 375 m resolution CGF SCE map product from VIIRS. These CGF products include quality-assurance data such as cloud-persistence statistics showing the age of the observation in each pixel. The objective of this paper is to introduce the new MODIS and VIIRS standard CGF daily SCE products and to provide a preliminary evaluation of uncertainties in the gap-filling methodology so that the products can be used as the basis for a moderate-resolution Earth science data record (ESDR) of SCE. Time series of the MODIS and VIIRS CGF products have been developed and evaluated at selected study sites in the US and southern Canada. Observed differences, although small, are largely attributed to cloud masking and differences in the time of day of image acquisition. A nearly 3-month time-series comparison of Terra MODIS and S-NPP VIIRS CGF snow-cover maps for a large study area covering all or parts of 11 states in the western US and part of southwestern Canada reveals excellent correspondence between the Terra MODIS and S-NPP VIIRS products, with a mean difference of 11 070 sqkm, which is ∼0.45 % of the study area. According to our preliminary validation of the Terra and Aqua MODIS CGF SCE products in the western US study area, we found higher accuracy of the Terra product compared with the Aqua product. The MODIS CGF SCE data record beginning in 2000 has been extended into the VIIRS era, which should last at least through the early 2030s.

Hall, Dorothy K.↗

Fused Bead Analysis of Diogenite Meteorites

Bulk rock chemistry is an essential dataset in meteoritics and planetary science [1]. A common method used to obtain the bulk chemistry of meteorites is ICP-MS. While the accuracy, precision and low detection limits of this process are advantageous [2], the sample size used for analysis (approx.70 mg) can be a problem in a field where small and finite samples are the norm. Fused bead analysis is another bulk rock analytical technique that has been used in meteoritics [3]. This technique involves forming a glass bead from ~10 mg of sample and measuring its chemistry using a defocused beam on a microprobe. Though the ICP-MS has lower detection limits than the microprobe, the fused bead method destroys a much smaller sample of the meteorite. Fused bead analysis was initially designed for samples with near-eutectic compositions and low viscosities. Melts generated of this type homogenize at relatively low temperatures and produce primary melts near the sample s bulk composition [3]. The application of fused bead analysis to samples with noneutectic melt compositions has not been validated. The purpose of this study is to test if fused bead analysis can accurately determine the bulk rock chemistry of non-eutectic melt composition meteorites. To determine this, we conduct two examinations of the fused bead. First, we compare ICP-MS and fused bead results of the same samples using statistical analysis. Secondly, we inspect the beads for the presence of crystals and chemical heterogeneity. The presence of either of these would indicate incomplete melting and quenching of the bead.

Mittlefehldt, D.W.↗

Ocean Wave Studies with Applications to Ocean Modeling and Improvement of Satellite Altimeter Measurements

Combining analysis of satellite data (altimeter, scatterometer, high-resolution visible and infrared images, etc.) with mathematical modeling of non-linear wave processes, we investigate various ocean wave fields (on scales from capillary to planetary), their role in ocean dynamics and turbulent transport (of heat and biogeochemical quantities), and their effects on satellite altimeter measuring accuracy. In 1998 my attention was focused on long internal gravity waves (10 to 1000 km), known also as baroclinic inertia-gravity (BIG) waves. We found these waves to be a major factor of altimeter measurements "noise," resulting in a greater uncertainty [up to 10 cm in terms of sea surface height (SSH) amplitude] in the measured SSH signal than that caused by the sea state bias variations (up to 5 cm or so). This effect still remains largely overlooked by the satellite altimeter community. Our studies of BIG waves address not only their influence on altimeter measurements but also their role in global ocean dynamics and in transport and turbulent diffusion of biogeochemical quantities. In particular, in collaboration with Prof Peter Weichman, Caltech, we developed a theory of turbulent diffusion caused by wave motions of most general nature. Applied to the problem of horizontal turbulent diffusion in the ocean, the theory yielded the effective diffusion coefficient as a function of BIG wave parameters obtainable from satellite altimeter data. This effort, begun in 1997, has been successfully completed in 1998. We also developed a theory that relates spatial fluctuations of scalar fields (such as sea surface temperature, chlorophyll concentration, drifting ice concentration, etc.) to statistical characteristics of BIG waves obtainable from altimeter measurements. A manuscript is in the final stages of preparation. In order to verify the theoretical predictions and apply them to observations, we are now analyzing Sea-viewing Wide Field of view Sensor (SeaWiFS) and Field of view Sensor (SeaWiFS) and Advanced Very High-Resolution Radiometer (AVHRR) data on sea surface temperature (SST) and chlorophyll concentration jointly with TOPEX/POSEIDON data on SSH variations.

Glazman, Roman E.↗

Trend Change Detection in NDVI Time Series: Effects of Inter-Annual Variability and Methodology

Changing trends in ecosystem productivity can be quantified using satellite observations of Normalized Difference Vegetation Index (NDVI). However, the estimation of trends from NDVI time series differs substantially depending on analyzed satellite dataset, the corresponding spatiotemporal resolution, and the applied statistical method. Here we compare the performance of a wide range of trend estimation methods and demonstrate that performance decreases with increasing inter-annual variability in the NDVI time series. Trend slope estimates based on annual aggregated time series or based on a seasonal-trend model show better performances than methods that remove the seasonal cycle of the time series. A breakpoint detection analysis reveals that an overestimation of breakpoints in NDVI trends can result in wrong or even opposite trend estimates. Based on our results, we give practical recommendations for the application of trend methods on long-term NDVI time series. Particularly, we apply and compare different methods on NDVI time series in Alaska, where both greening and browning trends have been previously observed. Here, the multi-method uncertainty of NDVI trends is quantified through the application of the different trend estimation methods. Our results indicate that greening NDVI trends in Alaska are more spatially and temporally prevalent than browning trends. We also show that detected breakpoints in NDVI trends tend to coincide with large fires. Overall, our analyses demonstrate that seasonal trend methods need to be improved against inter-annual variability to quantify changing trends in ecosystem productivity with higher accuracy.

Inter-Annual↗