Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Global sensitivity analysis”

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 415 records · Page 23

Uncertainty of Passive Imager Cloud Optical Property Retrievals to Instrument Radiometry and Model Assumptions: Examples from MODIS

The optical and microphysical structure of clouds is of fundamental importance for understanding a variety of cloud radiation and precipitation processes. With the advent of MODIS on the NASA EOS Terra and Aqua platforms, simultaneous global-daily 1 km retrievals of cloud optical thickness (COT) and effective particle radius (CER) are provided, as well as the derived water path (WP). The cloud product (MOD06/MYD06 for MODIS Terra and Aqua, respectively) provides separate retrieval datasets for various two-channel retrievals, typically a VISNIR channel paired with a 1.6, 2.1, and 3.7 m spectral channel. The MOD06 forward model is derived from on a homogeneous plane-parallel cloud. In Collection 5 processing (completed in 2007 with a modified Collection 5.1 completed in 2010), pixel-level retrieval uncertainties were calculated for the following non-3-D error sources: radiometry, surface spectral albedo, and atmospheric corrections associated with model analysis uncertainties (water vapor only). The latter error source includes error correlation across the retrieval spectral channels. Estimates of uncertainty in 1 aggregated (Level-3) means were also provided assuming unity correlation between error sources for all pixels in a grid for a single day, and zero correlation of error sources from one day to the next. I n Collection 6 (expected to begin in late summer 2013) we expanded the uncertainty analysis to include: (a) scene-dependent calibration uncertainty that depends on new band and detector-specific Level 1B uncertainties, (b) new model error sources derived from the look-up tables which includes sensitivities associated with wind direction over the ocean and uncertainties in liquid water and ice effective variance, (c) thermal emission uncertainties in the 3.7 m band associated with cloud and surface temperatures that are needed to extract reflected solar radiation from the total radiance signal, (d) uncertainty in the solar spectral irradiance at 3.7 m, and (e) addition of stratospheric ozone uncertainty in visible atmospheric corrections. A summary of the approach and example Collection 6 results will be shown.

Cloud Properties Observations↗

RH1020 Single Event Clock Upset Summary Report

This report summarizes the testing and analysis of "single event clock upset' in the RH1020. Also included are SEU-rate predictions and design recommendations for risk analysis and reduction. The subject of "upsets" in the RH1020 is best understood by using a model consisting of a global clock buffer and a D-type flip-flop as the basic memory unit. The RH1020 is built on the ACT 1 family architecture. As such, it has one low-skew global clock buffer with a TTL-level input threshold that is accessed via a single dedicated pin. The clock signal is driven to full CMOS levels, buffered, and sent to individual row buffers with one buffer per channel. For low-skew performance, the outputs of all of the RH1020 row buffers are shorted together via metal lines, as is done in the A1020B. All storage in the RH1020 consists of routed flip-flops, constructed with multiplexors and feedback through the routing segments. A simple latch can be constructed from a single (combinatorial or C) module; an edge-triggered flip-flop is constructed using two concatenated latches. There is no storage in the I/O modules. The front end of the clock buffering circuitry, at a common point relative to the row buffer, is a sub-circuit that was determined to be the most susceptible to heavy ions. This is due, in part, to its smaller transistors compared to the rest of the circuitry. This conclusion is also supported by SPICE simulations and an analysis of the heavy ion data, described in this report. The edge triggered D flip-flop has two single-event-upset modes. Mode one, called C-module upset, is caused by a heavy ion striking the C-module's sensitive area on the silicon and produces a soft single bit error at the output of the flip-flop. Mode two, called clock upset, is caused by a heavy ion strike on the clock buffer, generating a runt pulse interpreted as a false clock signal and consequently producing errors at the flip-flop outputs. C-module upset sensitivity in the RH1020 is essentially the same as that of its ACT 1 siblings (A1020, A1020A and A1020B), which were well tested, analyzed, and documented in the literature.

Katz, Richard B.↗

Validation Assessment for the Soil Moisture Active Passive (SMAP) Level 4 Carbon (L4_C) Data Product Version 5

The post-launch Cal/Val phase of the SMAP mission is guided by two primary objectives for each science product team: 1) to calibrate, verify, and improve the performance of the science algorithms, and 2) validate accuracies of the science data products as specified in the SMAP Level-1 mission science requirements. Algorithm science and product maintenance activities during the SMAP extended mission phase have also involved periodic algorithm calibration and product refinements to maintain or enhance product consistency and performance as well as science utility. This report provides an assessment of the latest (Version 5) SMAP Level 4 Carbon (L4_C) product. The L4_C Version 5 (v5) global record now spans more than six years (March 2015 – present) of SMAP operations and has benefited from five major reprocessing updates to the operational product. These reprocessing events and L4_C product release updates have incorporated various algorithm refinements and calibration adjustments to account for similar refinements to the upstream GEOS land model assimilation system, SMAP brightness temperatures, and MODIS vegetation inputs used for L4_C processing. The SMAP L4_C algorithms utilize a terrestrial carbon flux model informed by daily surface and root zone soil moisture information contributed from the SMAP Level 4 Soil Moisture (L4_SM) product along with optical remote sensing-based (e.g. MODIS-based) land cover and canopy fractional photosynthetic active radiation (fPAR), and other ancillary biophysical data. The carbon flux model estimates global daily net ecosystem CO2 exchange (NEE) and the component carbon fluxes, namely, vegetation gross primary production (GPP) and soil heterotrophic respiration (Rh). Other L4_C product elements include surface (~0-5 cm depth) soil organic carbon (SOC) stocks and associated environmental constraints to these processes, including soil moisture-related controls on GPP and ecosystem respiration (Kimball et al. 2014, Jones et al. 2017). The L4_C product addresses SMAP carbon cycle science objectives by: 1) providing a direct link between terrestrial carbon fluxes and underlying freeze/thaw and soil moisture-related constraints to these processes, 2) documenting primary connections between terrestrial water, energy and carbon cycles, and 3) improving understanding of terrestrial carbon sink activity. The SMAP L4_C algorithms and operational product are mature and at a CEOS Validation Stage 4 level (Jackson et al. 2012) based on extensive validation of the multi-year record against a diverse array of independent benchmarks, well characterized global performance, and systematic refinements gained from five major reprocessing events. There are no Level-1 mission science requirements for the L4_C product; however, self-imposed requirements have been established focusing on NEE as the primary product field for validation, and on demonstrating L4_C accuracy and success in meeting product science requirements (Jackson et al. 2012). The other L4_C product fields also have strong utility for carbon science applications (e.g., Liu et al. 2019, Endsley et al. 2020); however, analysis of these other fields is considered secondary relative to primary validation activities focusing on NEE. The L4_C targeted accuracy requirements are to meet or exceed a mean unbiased root-mean-square error (ubRMSE, or standard deviation of the error) for NEE of 1.6 g C m-2 d-1 and 30 g C m-2 yr-1, emphasizing northern (≥45°N) boreal and arctic ecosystems; this accuracy is similar to that of tower eddy covariance measurement-based observations (Baldocchi 2008). Methods used for the latest v5 L4_C product performance and validation assessment have been established from the SMAP Cal/Val plan and previous studies (Jackson et al. 2012, Jones et al. 2017) and include: 1) consistency evaluations of the product fields against earlier product releases (version 4 or earlier); 2) comparisons of daily carbon flux estimates with independent tower eddy 4covariance measurement-based daily carbon (CO2) flux observations from core tower validation sites (CVS); and 3) consistency checks against other global carbon products, including soil carbon inventory records, global GPP records derived from tower observation upscaling methods, and satellite-based observations of canopy solar induced chlorophyll fluorescence (SIF) as a surrogate for GPP. Metrics used to evaluate relative agreement between L4_C product fields and observational benchmarks include correlation (r-value), RMSE differences, bias and model sensitivity diagnostics. Following these validation criteria, the present report provides a validation assessment of the latest L4_C product release (v5). Detailed descriptions of the L4_C algorithm and additional global product accuracy and performance results are given elsewhere (Jones et al. 2017, Endsley et al. 2020). The v5 L4_C product replaces earlier product versions and continues to show: (i) accuracy and performance levels meeting or exceeding SMAP L4_C science requirements; (ii) improvement over the previous product version (version 4); and (iii) suitability for a diversity of science applications. Example L4_C applications from the recent literature include clarifying environmental trends and controls on the northern terrestrial carbon sink (Liu et al. 2019), diagnosing drought-related impacts on ecosystem productivity (Li et al. 2020), and regional monitoring of cropland conditions for projecting annual yields (Wurster et al. 2020). 2EXPECTED L4_C ALGORITHM AND PRODUCTPERFORMANCE The L4_C algorithm performance, including variance and uncertainty estimates of model outputs, was determined during the mission pre-launch phase through spatially explicit model sensitivity studies using available model inputs similar to those currently being used for operational production and evaluating the resulting model simulations over the observed range of northern (≥45 °N) and global conditions (Kimball et al. 2012, Entekhabi et al. 2014). The L4_C algorithm options were also evaluated during the mission prelaunch phase, including deriving canopy fPAR from lower order NDVI (Normalized Difference Vegetation Index) inputs in lieu of using MODIS (MOD15) fPAR; and including an explicit model representation of boreal fire disturbance recovery impacts. These results indicated that the L4_C accuracy requirements (i.e., NEE ubRMSE ≤ 30 g C m-2 yr-1or ≤ 1.6 g C m-2 d-1) could be met from the baseline algorithms over more than 82% and 89% of global and northern vegetated land areas, respectively (Yi et al. 2013, Kimball et al. 2014). The global L4_C algorithm error budget for NEE derived during the mission prelaunch phase indicated that the estimated NEE ubRMSE uncertainty is proportional to GPP and is therefore larger in higher biomass productivity areas, including forests and croplands (Kimball et al. 2014). Likewise, NEE ubRMSE uncertainty is expected to be lower in less-productive areas, including grasslands and shrublands. Expected model NEE ubRMSE levels were also generally within targeted accuracy levels for characteristically less-productive boreal and Arctic biomes, even though relative model error as a proportion of total productivity (NEE RMSE / GPP) may be large in these areas. The estimated NEE uncertainty was lower than expected in some warmer tropical high biomass productivity areas (e.g. Amazon rainforest) because of reduced low temperature and moisture constraints to the L4_C respiration calculations so that the bulk of model uncertainty is contributed by GPP in these areas. Model NEE uncertainty in the African Congo was estimated to

SMAP↗

The Effects of Normal Metal Stripes on TES Performance

Exploring the effects of size and geometry of normal metal features on the transition shapes and performance of transition-edge sensor microcalorimeters. The spectral resolution of transition-edge sensor (TES) microcalorimeters is very sensitive to the specific dependencies of the resistance R in the superconducting transition on the current I, magnetic field B, and temperature T. In particular, it has been shown that transitions that are very steep in (R,T) space lead to a significant noise term, in excess of conventional expectations. This so-called unexplained noise is known to be reduced by the addition of normal metal stripes across the TES perpendicular to the direction of current flow. These normal metal stripes have been shown to drastically alter the oscillatory patterns seen in measurements of the critical current as a function of magnetic field. However, there are many remaining questions about the exact impact of the stripes on current distributions within the TES, the Fraunhofer pattern and, therefore, the shape of the R(I, B, T) surface. Through measurements of the resistance under DC bias of TES devices of various sizes, with different stripe patterns and dimensions, we will discuss how these stripes can affect the R(I, B, T) surface. In addition, using measurements and analysis of the noise spectra of various devices we will present how these changes to the stripe pattern may affect the performance of the TES. In particular, we will discuss strategies to reduce the presence of localized discontinuities in the derivative of R, associated with increased noise, while maintaining the globally low levels of unexplained noise currently achieved with conventional metal stripe patterns. Implementing these strategies is a path towards producing large arrays with highly uniform transitions and high spectral resolution. These large uniform arrays will be required for future x-ray astronomy applications, such as the X-IFU on ATHENA.

TES↗

Negative extreme events in gross primary productivity and their drivers in China during the past three decades

Climate extremes have remarkable impacts on ecosystems and are expected to increase with future global warming. However, only few studies have focused on the ecological extreme events and their drivers in China. In this study, we carried out an analysis of negative extreme events in gross primary productivity (GPP) in China and the sub-regions during 1982–2015, using monthly GPP simulated by 12 process-based models (TRENDYv6) and an observation-based model (Yao-GPP). Extremes were defined as the negative 5th percentile of GPP anomalies, which were further merged into individual extreme events using a three-dimensional contiguous algorithm. Spatio-temporal patterns of negative GPP anomalies were analyzed by taking the 1000 largest extreme events into consideration. Results showed that the effects of extreme events decreased annual GPP by 2.8% (i.e. 208 TgC/year) in TRENDY models and 2.3% (i.e. 151 TgC/year) in Yao-GPP. Hotspots of extreme GPP deficits were mainly observed in North China (−53 gC/sq.m year) in TRENDY models and Northeast China (−42 gC/sq.m year) in Yao-GPP. For China as a whole, attribution analyses suggested that extreme low precipitation was associated with 40%–50% of extreme negative GPP events. Most events in northern and western China could be explained by meteorological droughts (i.e. low precipitation) while GPP extreme events in southern China were more associated with temperature extremes, in particular with cold spells. GPP was revealed to be much more sensitive to heat/drought than to cold/wet extreme events. Combined with projected changes in climate extremes in China, GPP negative anomalies caused by drought events in northern China and by temperature extremes in southern China might be more prominent in the future.

Climate change↗

Improvement of Automated POST Case Success Rate Using Support Vector Machines

During early conceptual design of complex systems, concept down selection can have a large impact upon program life-cycle cost. Therefore, any concepts selected during early design will inherently commit program costs and affect the overall probability of program success. For this reason it is important to consider as large a design space as possible in order to better inform the down selection process. For conceptual design of launch vehicles, trajectory analysis and optimization often presents the largest obstacle to evaluating large trade spaces. This is due to the sensitivity of the trajectory discipline to changes in all other aspects of the vehicle design. Small deltas in the performance of other subsystems can result in relatively large fluctuations in the ascent trajectory because the solution space is non-linear and multi-modal [1]. In order to help capture large design spaces for new launch vehicles, the authors have performed previous work seeking to automate the execution of the industry standard tool, Program to Optimize Simulated Trajectories (POST). This work initially focused on implementation of analyst heuristics to enable closure of cases in an automated fashion, with the goal of applying the concepts of design of experiments (DOE) and surrogate modeling to enable near instantaneous throughput of vehicle cases [2]. Additional work was then completed to improve the DOE process by utilizing a graph theory based approach to connect similar design points [3]. The conclusion of the previous work illustrated the utility of the graph theory approach for completing a DOE through POST. However, this approach was still dependent upon the use of random repetitions to generate seed points for the graph. As noted in [3], only 8% of these random repetitions resulted in converged trajectories. This ultimately affects the ability of the random reps method to confidently approach the global optima for a given vehicle case in a reasonable amount of time. With only an 8% pass rate, tens or hundreds of thousands of reps may be needed to be confident that the best repetition is at least close to the global optima. However, typical design study time constraints require that fewer repetitions be attempted, sometimes resulting in seed points that have only a handful of successful completions. If a small number of successful repetitions are used to generate a seed point, the graph method may inherit some inaccuracies as it chains DOE cases from the non-global-optimal seed points. This creates inherent noise in the graph data, which can limit the accuracy of the resulting surrogate models. For this reason, the goal of this work is to improve the seed point generation method and ultimately the accuracy of the resulting POST surrogate model. The work focuses on increasing the case pass rate for seed point generation.

Zwack, Matthew R.↗

NASA Team 2 Sea Ice Concentration Algorithm Retrieval Uncertainty

Satellite microwave radiometers are widely used to estimate sea ice cover properties (concentration, extent, and area) through the use of sea ice concentration (IC) algorithms. Rare are the algorithms providing associated IC uncertainty estimates. Algorithm uncertainty estimates are needed to assess accurately global and regional trends in IC (and thus extent and area), and to improve sea ice predictions on seasonal to interannual timescales using data assimilation approaches. This paper presents a method to provide relative IC uncertainty estimates using the enhanced NASA Team (NT2) IC algorithm. The proposed approach takes advantage of the NT2 calculations and solely relies on the brightness temperatures (TBs) used as input. NT2 IC and its associated relative uncertainty are obtained for both the Northern and Southern Hemispheres using the Advanced Microwave Scanning Radiometer for the Earth Observing System (AMSR-E) TB. NT2 IC relative uncertainties estimated on a footprint-by-footprint swath-by-swath basis were averaged daily over each 12.5-km grid cell of the polar stereographic grid. For both hemispheres and throughout the year, the NT2 relative uncertainty is less than 5%. In the Southern Hemisphere, it is low in the interior ice pack, and it increases in the marginal ice zone up to 5%. In the Northern Hemisphere, areas with high uncertainties are also found in the high IC area of the Central Arctic. Retrieval uncertainties are greater in areas corresponding to NT2 ice types associated with deep snow and new ice. Seasonal variations in uncertainty show larger values in summer as a result of melt conditions and greater atmospheric contributions. Our analysis also includes an evaluation of the NT2 algorithm sensitivity to AMSR-E sensor noise. There is a 60% probability that the IC does not change (to within the computed retrieval precision of 1%) due to sensor noise, and the cumulated probability shows that there is a 90% chance that the IC varies by less than +/-3%. We also examined the daily IC variability, which is dominated by sea ice drift and ice formation/melt. Daily IC variability is the highest, year round, in the MIZ (often up to 20%, locally 30%). The temporal and spatial distributions of the retrieval uncertainties and the daily IC variability is expected to be useful for algorithm intercomparisons, climate trend assessments, and possibly IC assimilation in models.

Sea Ice↗

Building a standardized Observing System Simulation Experiment (OSSE) framework for Mars

We advocate that the Decadal Survey recommends the NASA Science Mission Directorate to develop a rigorous Observing System Simulation Experiment (OSSE) framework for Mars, to optimize future atmospheric observations. Atmospheric conditions on Mars are a potential hazard source for landing missions. Errors in the estimates of atmospheric density profiles, inadequate knowledge of wind vertical structure and dust concentration as a function of height are likely causes of uncertainty at the landing site on the order of kilometers. An operational real-time weather forecasting capability for Mars would reduce such uncertainties, carrying enormous benefits to future robotic missions, and would be an invaluable prerequisite for human missions.A real-time forecasting capability relies upon three fundamental components: a critical mass of observing systems, a data assimilation system (DAS), and a global forecast model. The DAS allows the model to ingest the data effectively, optimizing the observational information content,and transforming them into a gridded representation of the atmosphere at a given time, called an ‘analysis’. The analysis is the best estimate of the atmospheric state for that time, and also represents a set of ‘initial conditions’ from which a global model can be initialized, to predict a future state of the atmosphere. The connection between analysis and forecast represents the foundation of modern weather forecasting. However, from the point of view of a forecast system,not all observations are equally impactful, partially because of the problem of “observational error correlation”, one important research topic in data assimilation development. For the Earth, partly due to the spontaneous and deregulated development of observations and forecast capabilities worldwide for more than half a century,the use of observations in contemporary operational forecast systems is suboptimal, with many potentially useful data being underutilized. On the contrary, Mars atmospheric scientists are in the unique situation of designing the next-generation observing systems by learning from the experience gathered on the Earth, so as to assure that the future instruments are specifically optimized to give the maximum benefit to a future weather forecast capability.An immensely powerful tool that has been firmly established by atmospheric scientists on the Earth is represented by a properly designed OSSE framework. A realistic OSSE framework cannot only quantify the benefit of future data types, be them surface based or space borne, but can also help design and optimize an entire observational network. Furthermore, OSSEs can provide deep insights into an atmosphere’s behavior, by addressing conceptual problems of its intrinsic predictability and delineating the regions or features of the atmosphere which are more sensitive to additional data and would benefit from a denser sampling. The difficulties posed by OSSEs are fundamentally different for Earth and Mars. For Earth, the enormous data volume imposes a tremendous constraint on any innovation in the observing systems: it is very hard for a single sensor to impact the skill. For Mars, the problem is the opposite: almost any additional instrument will exert some impact. However, OSSEs can help to evaluate the cost/benefit for every sensor and suggest optimal data configuration and density.The purpose of this white paper is to provide an introduction to a rigorously designed OSSE framework, explain the underlying problems and challenges, and engage the Mars community to collaborate with Earth Atmospheric scientists in order to develop a joint-OSSE framework for Mars with the largest consensual basis possible. An OSSE infrastructure would increase the understanding of the Martian atmosphere, would help NASA to optimize instrument specifications and orbit choice, providing the maximium benefit for a given expenditure of resources, and could even help establishing a roadmap for a future real-time weather forecasting capability.

Oreste Reale↗

Understanding the Uncertainty and Variability of Earth’s Radiation Budget Components Over the Past 3 Decades from the NASA’s GEWEX Surface Radiation Budget Release 4 Integrated Product

The NASA/GEWEX Surface Radiation Budget (SRB) project produced a 30+ year time series of 3-hourly Top-of-Atmosphere (TOA) and surface radiative flux components at a 1x1 degree resolution. This SRB Release 4 Integrated Productare derived using imager data products from the cross-calibrated from ISCCP (International Satellite Cloud Climatology Project) “HSeries” radiance, cloud and surface property data products. The ISCCP HXS product replaces ISCCP DX with a ninefold increase in pixel count (10km instead of 30km). ISCCP also provides an atmospheric temperature and moisture dataset known as nnHIRS which we use and discuss radiative flux sensitivities to in this presentation. Total solar irradiance is now variable consistent with SORCE measurements. Sea surface (SeaFlux v2.0) and land surface products (LandFlux) are utilized to achieve an radiative flux estimatesintegrated with global data sets from these GEWEX projects. The variable aerosol composition are specified using a detailed aerosol history from the Max Planck Institute Aerosol Climatology (MAC) also used for background aerosol in the ISCCP HX. Numerous other improvements were made to the SW and LW algorithms including improved cloud optical property treatments and specification of surface spectral albedo and emissivity for various surface conditions. Here we analyze the new SWand LWSRB datasets, comparing them to the previous Release 3, BSRN, GEBA and PMEL surface measurements, and ERBE and CERES satellite datasets. For surface flux validation besides ensemble comparisons, we show the variability of SRB vs surface measurements from BSRN beginning in 1992 and GEBA SW from 1983. We assess the variability of the surface flux anomalies, particularly the net fluxes and the cloud radiative effects, over the entire time period including the period containing CERES overlap. For the TOA fluxes, the anomalies of the fluxes in the early period are compared to the MEaSURES ERBE and the latter period to CERES EBAF/SYN1Deg datasets for a benchmark. Long-term changes in the surface radiation budget components and cloud radiative effects are shown and discussed relative to CERES and surface measurements. This includes a change point analysis that helps associate various changes in the ISCCP GEO/LEO configurations and other inputs to variations in the radiative budget components.

Earth Radiation Budget↗

Air Pollution Over the Ganges Basin and Northwest Bay of Bengal in the Early Postmonsoon Season Based on NASA MERRAero Data

The MERRA Aerosol Reanalysis (MERRAero) has been recently developed at NASA's Global Modeling Assimilation Office. This reanalysis is based on a version of the Goddard Earth Observing System-5 (GEOS-5) model radiatively coupled with Goddard Chemistry, Aerosol, Radiation, and Transport aerosols, and it includes assimilation of bias-corrected aerosol optical thickness (AOT) from the Moderate Resolution Imaging Spectroradiometer (MODIS) sensor on both Terra and Aqua satellites. In October over the period 2002-2009, MERRAero showed that AOT was lower over the east of the Ganges basin than over the northwest of the Ganges basin: this was despite the fact that the east of the Ganges basin should have produced higher anthropogenic aerosol emissions because of higher population density, increased industrial output, and transportation. This is evidence that higher aerosol emissions do not always correspond to higher AOT over the areas where the effects of meteorological factors on AOT dominate those of aerosol emissions. MODIS AOT assimilation was essential for correcting modeled AOT mainly over the northwest of the Ganges basin, where AOT increments were maximal. Over the east of the Ganges basin and northwest Bay of Bengal (BoB), AOT increments were low and MODIS AOT assimilation did not contribute significantly to modeled AOT. Our analysis showed that increasing AOT trends over northwest BoB (exceeding those over the east of the Ganges basin) were reproduced by GEOS-5, not because of MODIS AOT assimilation butmainly because of the model capability of reproducing meteorological factors contributing to AOT trends. Moreover, vertically integrated aerosol mass flux was sensitive to wind convergence causing aerosol accumulation over northwest BoB.

AOT↗

New Mastcam Multispectral Rock Classes in Sulfate-Bearing Strata, Gale Crater, Mars

Orbital data of Gale Crater, Mars has identified a transition in the stratigraphy of Mt. Sharp indicating an environmental change from a wetter one that accommodated clay mineral formation to a drier environment that led to the precipitation of sulfates. Mt. Sharp is not the only location on Mars for which this transition has been observed which suggests a global environmental change around the Noachian-Hesperian transition. Also visible from orbit is a dark-toned marker band that lies within the sulfate-bearing strata hypothesized to be a volcanic ash deposit. The Mars Science Laboratory (MSL) Curiosity rover is currently exploring this region of Mt. Sharp. In-situ analyses by the rover show changes in morphology and chemistry between the sulfate-bearing layers below and above the marker band and the marker band itself. The MSL instrument payload includes the Mast Camera (Mastcam), which collects multispectral images to provide context for other instruments and broad mineralogic interpretations. In addition to mafic Fe-bearing minerals, Mastcam is sensitive to oxides and some hydrated minerals, important indicators of alteration. In our previous work, we compiled a database of Mastcam spectra through Curiosity’s exploration of Vera Rubin ridge (sols 0-2302), from which 9 spectral classes of rocks were identified (Fig. 1D). Each class has diagnostic spectral features that reflect a common mineralogic interpretation. More recently, as Curiosity traversed through the clay-rich region, Glen Torridon (GT), there was some variation in the rock spectral classes potentially due to nontronite and other phyllosilicates, but not enough to warrant new spectral classes (Fig. 1B). We hypothesize that new classes will appear in accordance with sulfate-bearing strata and the marker band, specifically appearing as a ‘downturn’ in the last two wavelength filters that is inherent, but not unique, to polyhydrated Mg-sulfates. Here, we expand upon our analysis of Mastcam multispectral observations through Curiosity’s second encounter with the marker band (up to sol 3672), comparing the spectral diversity within the sulfate-bearing unit and the marker band to the spectral classes encountered previously in the traverse. These analyses will provide a basis for comparison and aid in selecting targets for in-situ investigations as Curiosity continues its ascent of Mt. Sharp.

A M Eng↗

ICESat Observations of Arctic Sea Ice: A First Look

Analysis of near-coincident ICESat and RADARSAT imagery shows that the retrieved elevations from the laser altimeter are sensitive to new openings (containing thin ice or open water) in the sea ice cover as well as to surface relief of old and first-year ice. The precision of the elevation estimates, measured over relatively flat sea ice, is approx. 2 cm. Using the thickness of thin-ice in recent openings to estimate sea level references, we obtain the sea-ice freeboard along the altimeter tracks. This step is necessitated by the large uncertainties in the sea surface topography compared to that required for accurate determination of freeboard. Unknown snow depth introduces the largest uncertainty in the conversion of freeboard to ice thickness. Surface roughness is also derived, for the first time, from the variability of successive elevation estimates along the altimeter track. Overall, these ICESat measurements provide an unprecedented view of the Arctic Ocean ice cover at length scales at and above the spatial dimension of the altimeter footprint of approx. 70 m.

Arctic Ocean↗

Radiative Susceptibility of Cloudy Atmospheres to Droplet Number Perturbations: 1. Theoretical Analysis and Examples from MODIS

Theoretical and satellite-based assessments of the sensitivity of broadband shortwave radiative fluxes in cloudy atmospheres to small perturbations in the cloud droplet number concentration (N) of liquid water clouds under constant water conditions are performed. Two approaches to study this sensitivity are adopted: absolute increases in N, for which the radiative response is referred to as absolute cloud susceptibility, and relative increases in N or relative cloud susceptibility. Estimating the former is more challenging as it requires an assumed value for either cloud liquid water content or geometrical thickness; both susceptibilities require an assumed relationship between the droplet volume and effective radius. Expanding upon previous susceptibility studies, present radiative calculations include the effect of AN perturbations on droplet asymmetry parameter and single-scattering albedo, in addition to extinction. Absolute cloud susceptibility has a strong nonlinear dependence on the droplet effective radius as expected, while relative cloud susceptibility is primarily dependent on optical thickness. Molecular absorption and reflecting surfaces both reduce the relative contribution of the cloud to the top-of-atmosphere (TOA) flux and therefore also reduce the TOA albedo susceptibility. Transmittance susceptibilities are negative with absolute values similar to albedo susceptibility, while atmospheric absorptance susceptibilities are about an order of magnitude smaller than albedo susceptibilities and can be either positive or negative. Observation-based susceptibility calculations are derived from MODIS pixel-level retrievals of liquid water cloud optical thickness, effective radius, and cloud top temperature; two data granule examples are shown. Susceptibility quantifies the aerosol indirect effect sensitivity in a way that can be easily computed from model fields. As such, susceptibilities derived from MODIS observations provide a higher-order test of model cloud properties used for indirect effect studies. MODIS-derived global distributions of cloud susceptibility and radiative forcing calculations are presented in a companion paper.

Platnick, Steven↗

Analysis of Present Day and Future OH and Methane Lifetime in the ACCMIP Simulations

Results from simulations performed for the Atmospheric Chemistry and Climate Modeling Intercomparison Project (ACCMIP) are analysed to examine how OH and methane lifetime may change from present day to the future, under different climate and emissions scenarios. Present day (2000) mean tropospheric chemical lifetime derived from the ACCMIP multi-model mean is 9.8+/-1.6 yr (9.3+/-0.9 yr when only including selected models), lower than a recent observationally-based estimate, but with a similar range to previous multi-model estimates. Future model projections are based on the four Representative Concentration Pathways (RCPs), and the results also exhibit a large range. Decreases in global methane lifetime of 4.5 +/- 9.1% are simulated for the scenario with lowest radiative forcing by 2100 (RCP 2.6), while increases of 8.5+/-10.4% are simulated for the scenario with highest radiative forcing (RCP 8.5). In this scenario, the key driver of the evolution of OH and methane lifetime is methane itself, since its concentration more than doubles by 2100 and it consumes much of the OH that exists in the troposphere. Stratospheric ozone recovery, which drives tropospheric OH decreases through photolysis modifications, also plays a partial role. In the other scenarios, where methane changes are less drastic, the interplay between various competing drivers leads to smaller and more diverse OH and methane lifetime responses, which are difficult to attribute. For all scenarios, regional OH changes are even more variable, with the most robust feature being the large decreases over the remote oceans in RCP8.5. Through a regression analysis, we suggest that differences in emissions of non-methane volatile organic compounds and in the simulation of photolysis rates may be the main factors causing the differences in simulated present day OH and methane lifetime. Diversity in predicted changes between present day and future OH was found to be associated more strongly with differences in modelled temperature and stratospheric ozone changes. Finally, through perturbation experiments we calculated an OH feedback factor (F) of 1.24 from present day conditions (1.50 from 2100 RCP8.5 conditions) and a climate feedback on methane lifetime of 0.33+-0.13 yr/K, on average. Models that did not include interactive stratospheric ozone effects on photolysis showed a stronger sensitivity to climate, as they did not account for negative effects of climate-driven stratospheric ozone recovery on tropospheric OH, which would have partly offset the overall OH/methane lifetime response to climate change.

atmospheric composition↗

Aquarius Instrument Science Calibration During the Risk Reduction Phase

This final report presents the results of work performed under NASA Grant NAG512726 during the period 15 January 2003 through 30 June 2004. An analysis was performed of a possible vicarious calibration method for use by Aquarius to monitor and stabilize the absolute and relative calibration of its microwave radiometer. Stationary statistical properties of the brightness temperature (T(sub B)) measured by a low Earth orbiting radiometer operating at 1.4135 GHz are considered as a means of validating its absolute calibration. The global minimum, maximum, and average T(sub B) are considered, together with a vicarious cold reference method that detects the presence of a sharp lower bound on naturally occurring values for T(sub B). Of particular interest is the reliability with which these statistics can be extracted from a realistic distribution of T(sub B) measurements that would be observed by a typical sensor. Simulations of measurements are performed that include the effects of instrument noise and variable environmental factors such as the global water vapor and ocean surface temperature, salinity and wind distributions. Global minima can vary widely due to instrument noise and are not a reliable calibration reference. Global maxima are strongly influenced by several environmental factors as well as instrument noise and are even less stationary. Global averages are largely insensitive to instrument noise and, in most cases, to environmental conditions as well. The global average T(sub B) varies at only the 0.1 K RMS level except in cases of anomalously high winds, when it can increase considerably more. The vicarious cold reference is similarly insensitive to instrument effects and most environmental factors. It is not significantly affected by high wind conditions. The stability of the vicarious reference is, however, found to be somewhat sensitive (at the several tenths of Kelvins level) to variations in the background cold space brightness, T(sub c). The global average is much less sensitive to this parameter and so using two approaches together can be mutually beneficial.

Ruf, Christopher S.↗

The Impact of Assimilating Large Volumes of GNSS Radio Occultation Observations from Spire’s Commercial Constellation into NASA’s Global Earth Observing System

The upcoming retrospective analysis for the 21st century (R21C) reanalysis product from NASA will include the full dataset of GNSS Radio Occultation (RO) observations collected by Spire with their constellation of smallsats. The Spire RO dataset, which was purchased by NASA for its Commercial Smallsat Data Acquisition (CSDA) archive, has global coverage and includes approximately 6 thousand RO profiles per day during 2019 increasing to nearly 20 thousand in 2022-2023. This more than doubles the volume of RO profiles from all other routinely assimilated RO missions combined during 2023. The increase in the number of RO profiles available is particularly important in the extratropics as the next largest RO constellation, the state-of-the-science Formosa Satellite Mission 7 (FORMOSAT-7)/Constellation Observing System for Meteorology, Ionosphere and Climate (COSMIC-2) mission (hereafter COSMIC-2), is focused on the Tropics and only provides approximately 4 thousand RO profiles per day. The large number of observations in the Spire RO dataset have great potential to improve analyses and forecasts of the Earth’s atmosphere produced by numerical weather prediction systems. Their impact is assessed through numerical experiments using NASA’s Global Earth Observing System (GEOS) Atmospheric Data Assimilation System (ADAS) with and without the observations from the COSMIC-2 mission. The ability of the Spire RO observations to make up for the omission of the RO observations from COSMIC-2 is directly assessed over the tropics and the impact of Spire over the extratropics is compared to that of the available RO missions that sample over this region. Finally, the quality of the Spire RO observations is compared to that of the routinely assimilated RO missions using Forecast Sensitivity-based Observation Impact (FSOI).

Michael J. Murphy, Jr.↗

The Impact of Assimilating Large Volumes of GNSS Radio Occultation Observations from Spire’s Commercial Constellation into NASA’s GEOS

The upcoming retrospective analysis for the 21st century (R21C) reanalysis product from NASA will include the full dataset of GNSS Radio Occultation (RO) observations collected by Spire with their constellation of smallsats. The Spire RO dataset, which was purchased by NASA for its Commercial Smallsat Data Acquisition (CSDA) archive, has global coverage and includes approximately 6 thousand RO profiles per day during 2019 increasing to nearly 20 thousand in 2022-2023. This more than doubles the volume of RO profiles from all other routinely assimilated RO missions combined during 2023. The increase in the number of RO profiles available is particularly important in the extratropics as the next largest RO constellation, the state-of-the-science Formosa Satellite Mission 7 (FORMOSAT-7)/Constellation Observing System for Meteorology, Ionosphere and Climate (COSMIC-2) mission (hereafter COSMIC-2), is focused on the Tropics and only provides approximately 4 thousand RO profiles per day. The large number of observations in the Spire RO dataset have great potential to improve analyses and forecasts of the Earth’s atmosphere produced by numerical weather prediction systems. Their impact is assessed through numerical experiments using NASA’s Global Earth Observing System (GEOS) Atmospheric Data Assimilation System (ADAS) with and without the observations from the COSMIC-2 mission. The ability of the Spire RO observations to make up for the omission of the RO observations from COSMIC-2 is directly assessed over the tropics and the impact of Spire over the extratropics is compared to that of the available RO missions that sample over this region. Finally, the quality of the Spire RO observations is compared to that of the routinely assimilated RO missions using Forecast Sensitivity-based Observation Impact (FSOI).

Michael J. Murphy↗

The Impact of Assimilating Large Volumes of GNSS Radio Occultation Observations from Spire’s Commercial Constellation into NASA’s GEOS

The upcoming retrospective analysis for the 21st century (R21C) reanalysis product from NASA will include the full dataset of GNSS Radio Occultation (RO) observations collected by Spire with their constellation of smallsats. The Spire RO dataset, which was purchased by NASA for its Commercial Smallsat Data Acquisition (CSDA) archive, has global coverage and increases in number from approximately 6 thousand RO profiles per day during 2019 to approximately 20 thousand in 2022-2023. This more than doubles the volume of RO profiles from all other routinely assimilated RO missions combined during 2022. The increase in the number of RO profiles available is particularly important in the extratropics as the next largest RO constellation, the state-of-the-science Formosa Satellite Mission 7 (FORMOSAT-7)/Constellation Observing System for Meteorology, Ionosphere and Climate (COSMIC-2) mission (hereafter COSMIC-2), is focused on the Tropics and only provides approximately 4 thousand RO profiles per day. The large number of observations in the Spire RO dataset have great potential to improve analyses and forecasts of the Earth’s atmosphere produced by numerical weather prediction systems. Their impact is assessed through numerical experiments using NASA’s Global Earth Observing System (GEOS) Atmospheric Data Assimilation System (ADAS) with and without the observations from the COSMIC-2 mission. The ability of the Spire RO observations to make up for the omission of the RO observations from COSMIC-2 is directly assessed over the tropics and the impact of Spire over the extratropics is compared to that of the available RO missions that sample over this region. Finally, the quality of the Spire RO observations is compared to that of the routinely assimilated RO missions using Forecast Sensitivity-based Observation Impact (FSOI).

Michael J. Murphy↗