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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 397 records · Page 22

All sources NASA NTRS DOE OSTI

NASA NTRS2025

Mind the Gap: Addressing Data Gaps and Assessing Noise Mismodeling in LISA

Due to the sheer complexity of the Laser Interferometer Space Antenna (LISA) space mission, data gaps arising from instrumental irregularities and/or scheduled maintenance are unavoidable. Focusing on merger-dominated massive black hole binary signals, we test the appropriateness of the Whittle-likelihood on gapped data in a variety of cases. From first principles, we derive the likelihood valid for gapped data in both the time and frequency domains. Cheap-to-evaluate proxies to p-p plots are derived based on a Fisher-based formalism, and verified through Bayesian techniques. Our tools allow to predict the altered variance in the parameter estimates that arises from noise mismodeling, as well as the information loss represented by the broadening of the posteriors. The result of noise mismodeling with gaps is sensitive to the characteristics of the noise model, with strong low-frequency (red) noise and strong high-frequency (blue) noise giving statistically significant fluctuations in recovered parameters. We demonstrate that the introduction of a tapering window reduces statistical inconsistency errors, at the cost of less precise parameter estimates. We also show that the assumption of independence between inter-gap segments appears to be a fair approximation even if the data set is inherently coherent. However, if one instead assumes fictitious correlations in the data stream, when the data segments are actually independent, then the resultant parameter recoveries could be inconsistent with the true parameters. The theoretical and numerical practices that are presented in this work could readily be incorporated into global-fit pipelines operating on gapped data.

LISA↗

NASA NTRS2017

Importance of a Priori Vertical Ozone Profiles for TEMPO Air Quality Retrievals

Ozone (O3) is a toxic pollutant which plays a major role in air quality. Typically, monitoring of surface air quality and O3 mixing ratios is conducted using in situ measurement networks. This is partially due to high-quality information related to air quality being limited from space-borne platforms due to coarse spatial resolution, limited temporal frequency, and minimal sensitivity to lower tropospheric and surface-level O3. The Tropospheric Emissions: Monitoring of Pollution (TEMPO) satellite is designed to address the limitations of current space-based platforms and to improve our ability to monitor North American air quality. TEMPO will provide hourly data of total column and vertical profiles of O3 with high spatial resolution to be used as a near-real-time air quality product. TEMPO O3 retrievals will apply the Smithsonian Astrophysical Observatory profile algorithm developed based on work from GOME (Global Ozone Monitoring Experiment), GOME-2, and OMI (Ozone Monitoring Instrument). This algorithm is suggested to use a priori O3 profile information from a climatological data-base developed from long-term ozone-sonde measurements (tropopause-based (TB-Clim) O3 climatology). This study evaluates the TB-Clim dataset and model simulated O3 profiles, which could potentially serve as a priori O3 profile information in TEMPO retrievals, from near-real-time data assimilation model products (NASA GMAO's (Global Modeling and Assimilation Office) operational GEOS-5 (Goddard Earth Observing System, Version 5) FP (Forecast Products) model and reanalysis data from MERRA2 (Modern-Era Retrospective analysis for Research and Applications, Version 2)) and a full chemical transport model (CTM), GEOS-Chem. In this study, vertical profile products are evaluated with surface (0-2 kilometers) and tropospheric (0-10 kilometers) TOLNet (Tropospheric Ozone Lidar Network) observations and the theoretical impact of individual a priori profile sources on the accuracy of TEMPO O3 retrievals in the troposphere and at the surface are presented. Results indicate that while the TB-Clim climatological dataset can replicate seasonally-averaged tropospheric O3 profiles, model-simulated profiles from a full CTM resulted in more accurate tropospheric and surface-level O3 retrievals from TEMPO when compared to hourly and daily-averaged TOLNet observations. Furthermore, it is shown that when large surface O3 mixing ratios are observed, TEMPO retrieval values at the surface are most accurate when applying CTM a priori profile information compared to all other data products.

Priori↗

NASA NTRS2014

Diagnosis of Middle Atmosphere Climate Sensitivity by the Climate Feedback Response Analysis Method

We present a new method to diagnose the middle atmosphere climate sensitivity by extending the Climate Feedback-Response Analysis Method (CFRAM) for the coupled atmosphere-surface system to the middle atmosphere. The Middle atmosphere CFRAM (MCFRAM) is built on the atmospheric energy equation per unit mass with radiative heating and cooling rates as its major thermal energy sources. MCFRAM preserves the CFRAM unique feature of an additive property for which the sum of all partial temperature changes due to variations in external forcing and feedback processes equals the observed temperature change. In addition, MCFRAM establishes a physical relationship of radiative damping between the energy perturbations associated with various feedback processes and temperature perturbations associated with thermal responses. MCFRAM is applied to both measurements and model output fields to diagnose the middle atmosphere climate sensitivity. It is found that the largest component of the middle atmosphere temperature response to the 11-year solar cycle (solar maximum vs. solar minimum) is directly from the partial temperature change due to the variation of the input solar flux. Increasing CO2 always cools the middle atmosphere with time whereas partial temperature change due to O3 variation could be either positive or negative. The partial temperature changes due to different feedbacks show distinctly different spatial patterns. The thermally driven globally averaged partial temperature change due to all radiative processes is approximately equal to the observed temperature change, ranging from 0.5 K near 70 km from the near solar maximum to the solar minimum.

MCFRAM↗

NASA NTRS2003

Global Changes of the Water Cycle Intensity

In this study, we evaluate numerical simulations of the twentieth century climate, focusing on the changes in the intensity of the global water cycle. A new diagnostic of atmospheric water vapor cycling rate is developed and employed, that relies on constituent tracers predicted at the model time step. This diagnostic is compared to a simplified traditional calculation of cycling rate, based on monthly averages of precipitation and total water content. The mean sensitivity of both diagnostics to variations in climate forcing is comparable. However, the new diagnostic produces systematically larger values and more variability than the traditional average approach. Climate simulations were performed using SSTs of the early (1902-1921) and late (1979- 1998) twentieth century along with the appropriate C02 forcing. In general, the increase of global precipitation with the increases in SST that occurred between the early and late twentieth century is small. However, an increase of atmospheric temperature leads to a systematic increase in total precipitable water. As a result, the residence time of water in the atmosphere increased, indicating a reduction of the global cycling rate. This result was explored further using a number of 50-year climate simulations from different models forced with observed SST. The anomalies and trends in the cycling rate and hydrologic variables of different GCMs are remarkably similar. The global annual anomalies of precipitation show a significant upward trend related to the upward trend of surface temperature, during the latter half of the twentieth century. While this implies an increase in the hydrologic cycle intensity, a concomitant increase of total precipitable water again leads to a decrease in the calculated global cycling rate. An analysis of the land/sea differences shows that the simulated precipitation over land has a decreasing trend while the oceanic precipitation has an upward trend consistent with previous studies and the available observations. The decreasing continental trend in precipitation is located primarily over tropical land regions, with some other regions, such as North America experiencing an increasing trend. Precipitation trends are diagnosed further using the water tracers to delineate the precipitation that occurs because of continental evaporation, as opposed to oceanic evaporation. These diagnostics show that over global land areas, the recycling of continental moisture is decreasing in time. However, the recycling changes are not spatially uniform so that some regions, most notably over the United States, experience continental recycling of water that increases in time.

Bosilovich, Michael G.↗

NASA NTRSDate not supplied

Assessing the Variability of Radiation, Water and Energy in the Deep Tropics Over the Last 3 Decades

Global surface energy closure and its variability depends heavily upon the surface radiative energy budget in the deep tropics. A key process observed is the change in convection regimes from disorganized convection to organized deep convection. During these periods an imbalance is observed between energy and water fluxes in terms of energy closure in some atmospheric reanalysis systems, implying that processes are not simulated well. A recent paper by Hsaio et al (2023, pre-published version) found that the transition process has wind shear and longwave cloud radiative feedback signatures. In this presentation, we assess the surface radiative budget in the deep tropics by comparing multiple data products (GEWEX SRB Rel4IP, CERES Ed4.1, ISSCP FH, etc.) and describing the variability across the deep tropics for the period from 1988 to near present. We assess this variability of the radiative flux anomalies (including the net TOA, surface and the atmospheric divergence fluxes) against water vapor divergence, cloud properties (include ISCCP “Weather States, Tselioudis et al. 2021) and larger scale wind shear. We include further analysis contrasting 6 key tropical oceanic regions (Indian Ocean east and west, tropical western, central and eastern Pacific Ocean, and Tropical Atlantic Ocean) and the sensitivity of the fluxes as a function of fluctuations in cloud types in response to various larger scale atmospheric processes (e.g. El Niño, Indian Ocean Dipole). These are contrasted against the modeled radiative fluxes from ERA-5 and MERRA-2. Ocean buoy measurements of radiative fluxes are utilized to help assess data radiative flux uncertainty over the nearly 40 year period. In general, the variability for overlap periods of these various data products agrees well, but there are significant differences in the net fluxes that vary according to the rendering of surface, atmospheric, aerosol and cloud properties. We conclude with recommendations for continuing work.

Radiation budget↗

NASA NTRSDate not supplied

Simultaneous lidar profiling of water vapor, aerosol, and wind over the eastern Pacific and implications for observational challenges on clouds, circulation, and climate sensitivity

Moist convective processes and related circulations are key components of the Earth’s climate system. Improving our understanding of moisture and cloud processes over oceans is especially challenging due to observational limitations, but also important due to the large radiative effects of cloud over ocean as well as the ties to local and global-scale circulations. A unique dataset for addressing these topics was collected over 5 flights with the NASA HALO (water vapor DIAL and aerosol HSRL) and DAWN (Doppler wind) lidar systems onboard the NASA DC-8 aircraft during the spring 2019 ADM-Aeolus Cal/Val test flight campaign. High-resolution water vapor, aerosol/cloud, and wind profiles were simultaneously collected and complimented with in situ high-resolution dropsonde measurements of temperature, relative humidity, and wind profiles. Airborne lidar measurements of these key variables captured both large-scale spatial variability across different cloud and atmospheric conditions as well as finer features, providing insight into processes from large-scale circulations down to aerosol/cloud interactions. The 5 flights spanned the northeast Pacific and southwest U.S. from 6 N to 52 N and 112 W to 157 W and captured several features of interest including moist layers from the ascending region of the Hadley circulation, very dry layers above midlatitude marine stratocumulus and cloud-free planetary boundary layers, and cloud top height and spatial variability. High covariance of aerosol and water vapor in the free troposphere was also observed for the first time. Preliminary analysis of these data sets will be presented as they apply to clouds and circulations, supplemented by model reanalysis to inform how regional water vapor structure and large-scale motion may contribute to observed shallow cloud formations and other meteorological processes that underpin our current understanding of how these state variables couple to impact the Earth’s weather and climate systems. The potential of complementary measurements in future field campaigns will also be mentioned.

Brian James Carroll↗

NASA NTRS2007

Performance of the Lidar Design and Data Algorithms for the GLAS Global Cloud and Aerosol Measurements

The Geoscience Laser Altimeter System (GLAS) launched in early 2003 is the first polar orbiting satellite lidar. The instrument design includes high performance observations of the distribution and optical scattering cross sections of atmospheric clouds and aerosol. The backscatter lidar operates at two wavelengths, 532 and 1064 nm. For the atmospheric cloud and aerosol measurements, the 532 nm channel was designed for ultra high efficiency with solid state photon counting detectors and etalon filtering. Data processing algorithms were developed to calibrate and normalize the signals and produce global scale data products of the height distribution of cloud and aerosol layers and their optical depths and particulate scattering cross sections up to the limit of optical attenuation. The paper will concentrate on the effectiveness and limitations of the lidar channel design and data product algorithms. Both atmospheric receiver channels meet and exceed their design goals. Geiger Mode Avalanche Photodiode modules are used for the 532 nm signal. The operational experience is that some signal artifacts and non-linearity require correction in data processing. As with all photon counting detectors, a pulse-pile-up calibration is an important aspect of the measurement. Additional signal corrections were found to be necessary relating to correction of a saturation signal-run-on effect and also for daytime data, a small range dependent variation in the responsivity. It was possible to correct for these signal errors in data processing and achieve the requirement to accurately profile aerosol and cloud cross section down to 10-7 llm-sr. The analysis procedure employs a precise calibration against molecular scattering in the mid-stratosphere. The 1064 nm channel detection employs a high-speed analog APD for surface and atmospheric measurements where the detection sensitivity is limited by detector noise and is over an order of magnitude less than at 532 nm. A unique feature of the GLAS is a full acquisition of the surface return pulse, which has important application to the atmospheric transmission retrieval.

Spinhirne, James D.↗

NASA NTRS2017

The PMIP4 Contribution to CMIP6-Part 4: Scientific Objectives and Experimental Design of the PMIP4-CMIP6 Last Glacial Maximum Experiments and PMIP4 Sensitivity Experiments

The Last Glacial Maximum (LGM, 21,000 years ago) is one of the suite of paleoclimate simulations included in the current phase of the Coupled Model Intercomparison Project (CMIP6). It is an interval when insolation was similar to the present, but global ice volume was at a maximum, eustatic sea level was at or close to a minimum, greenhouse gas concentrations were lower, atmospheric aerosol loadings were higher than today, and vegetation and land-surface characteristics were different from today. The LGM has been a focus for the Paleoclimate Modelling Intercomparison Project (PMIP) since its inception, and thus many of the problems that might be associated with simulating such a radically different climate are well documented. The LGM state provides an ideal case study for evaluating climate model performance because the changes in forcing and temperature between the LGM and pre-industrial are of the same order of magnitude as those projected for the end of the 21st century. Thus, the CMIP6 LGM experiment could provide additional information that can be used to constrain estimates of climate sensitivity. The design of the Tier 1 LGM experiment (lgm) includes an assessment of uncertainties in boundary conditions, in particular through the use of different reconstructions of the ice sheets and of the change in dust forcing. Additional (Tier 2) sensitivity experiments have been designed to quantify feedbacks associated with land-surface changes and aerosol loadings, and to isolate the role of individual forcings. Model analysis and evaluation will capitalize on the relative abundance of paleoenvironmental observations and quantitative climate reconstructions already available for the LGM.

Kageyama, Masa↗

NASA NTRS2022

Soil Moisture Estimation in South Asia via Assimilation of SMAP Retrievals

A soil moisture retrieval assimilation framework is implemented across South Asia in an attempt to improve regional soil moisture estimation as well as to provide a consistent regional soil moisture dataset. This study aims to improve the spatiotemporal variability of soil moisture estimates by assimilating Soil Moisture Active Passive (SMAP) near-surface soil moisture retrievals into a land surface model. The Noah-MP (v4.0.1) land surface model is run within the NASA Land Information System software framework to model regional land surface processes. NASA Modern-Era Retrospective Analysis for Research and Applications (MERRA2) and Global Precipitation Measurement (GPM) Integrated Multi-satellitE Retrievals (IMERG) provide the meteorological boundary conditions to the land surface model. Assimilation is carried out using both cumulative distribution function (CDF)-corrected (DA-CDF) and uncorrected SMAP retrievals (DA-NoCDF). CDF matching is applied to correct the statistical moments of the SMAP soil moisture retrieval relative to the land surface model. Comparison of assimilated and model-only soil moisture estimates with publicly available in situ measurements highlights the relative improvement in soil moisture estimates by assimilating SMAP retrievals. Across the Tibetan Plateau, DA-NoCDF reduced the mean bias and RMSE by 8.4 % and 9.4 %, even though assimilation only occurred during less than 10 % of the study period due to frozen (or partially frozen) soil conditions. The best goodness-of-fit statistics were achieved for the IMERG DA-NoCDF soil moisture experiment. The general lack of publicly available in situ measurements across irrigated areas limited a domain-wide direct model validation. However, comparison with regional irrigation patterns suggested correction of biases associated with an unmodeled hydrologic phenomenon (i.e., anthropogenic influence via irrigation) as a result of SMAP soil moisture retrieval assimilation. The greatest sensitivity to assimilation was observed in cropland areas. Improvements in soil moisture potentially translate into improved spatiotemporal patterns of modeled evapotranspiration, although limited influence from soil moisture assimilation was observed on modeled processes within the carbon cycle such as gross primary production. Improvement in fine-scale modeled estimates by assimilating coarse-scale retrievals highlights the potential of this approach for soil moisture estimation over data-scarce regions.

Jawairia Ahmad↗

NASA NTRS2013

Modeling of 2008 Kasatochi Volcanic Sulfate Direct Radiative Forcing: Assimilation of OMI SO2 Plume Height Data and Comparison with MODIS and CALIOP Observations

Volcanic SO2 column amount and injection height retrieved from the Ozone Monitoring Instrument (OMI) with the Extended Iterative Spectral Fitting (EISF) technique are used to initialize a global chemistry transport model (GEOS-Chem) to simulate the atmospheric transport and lifecycle of volcanic SO2 and sulfate aerosol from the 2008 Kasatochi eruption, and to subsequently estimate the direct shortwave, top-of-the-atmosphere radiative forcing of the volcanic sulfate aerosol. Analysis shows that the integrated use of OMI SO2 plume height in GEOS-Chem yields: (a) good agreement of the temporal evolution of 3-D volcanic sulfate distributions between model simulations and satellite observations from the Moderate Resolution Imaging Spectroradiometer (MODIS) and Cloud-Aerosol Lidar with Orthogonal Polarisation (CALIOP), and (b) an e-folding time for volcanic SO2 that is consistent with OMI measurements, reflecting SO2 oxidation in the upper troposphere and stratosphere is reliably represented in the model. However, a consistent (approx. 25 %) low bias is found in the GEOS-Chem simulated SO2 burden, and is likely due to a high (approx.20 %) bias of cloud liquid water amount (as compared to the MODIS cloud product) and the resultant stronger SO2 oxidation in the GEOS meteorological data during the first week after eruption when part of SO2 underwent aqueous-phase oxidation in clouds. Radiative transfer calculations show that the forcing by Kasatochi volcanic sulfate aerosol becomes negligible 6 months after the eruption, but its global average over the first month is -1.3W/sq m, with the majority of the forcing-influenced region located north of 20degN, and with daily peak values up to -2W/sq m on days 16-17. Sensitivity experiments show that every 2 km decrease of SO2 injection height in the GEOS-Chem simulations will result in a approx.25% decrease in volcanic sulfate forcing; similar sensitivity but opposite sign also holds for a 0.03 m increase of geometric radius of the volcanic aerosol particles. Both sensitivities highlight the need to characterize the SO2 plume height and aerosol particle size from space. While more research efforts are warranted, this study is among the first to assimilate both satellite-based SO2 plume height and amount into a chemical transport model for an improved simulation of volcanic SO2 and sulfate transport.

Wang, J.↗

NASA NTRS2015

An Assessment of Actual and Potential Building Climate Zone Change and Variability From the Last 30 Years Through 2100 Using NASA's MERRA and CMIP5 Simulations

Background: In the US, residential and commercial building infrastructure combined consumes about 40% of total energy usage and emits about 39% of total CO2 emission (DOE/EIA "Annual Energy Outlook 2013"). Building codes, as used by local and state enforcement entities are typically tied to the dominant climate within an enforcement jurisdiction classified according to various climate zones. These climate zones are based upon a 30-year average of local surface observations and are developed by DOE and ASHRAE. Establishing the current variability and potential changes to future building climate zones is very important for increasing the energy efficiency of buildings and reducing energy costs and emissions in the future. Objectives: This paper demonstrates the usefulness of using NASA's Modern Era Retrospective-analysis for Research and Applications (MERRA) atmospheric data assimilation to derive the DOE/ASHRAE building climate zone maps and then using MERRA to define the last 30 years of variability in climate zones for the Continental US. An atmospheric assimilation is a global atmospheric model optimized to satellite, atmospheric and surface in situ measurements. Using MERRA as a baseline, we then evaluate the latest Climate Model Inter-comparison Project (CMIP) climate model Version 5 runs to assess potential variability in future climate zones under various assumptions. Methods: We derive DOE/ASHRAE building climate zones using surface and temperature data products from MERRA. We assess these zones using the uncertainties derived by comparison to surface measurements. Using statistical tests, we evaluate variability of the climate zones in time and assess areas in the continental US for statistically significant trends by region. CMIP 5 produced a data base of over two dozen detailed climate model runs under various greenhouse gas forcing assumptions. We evaluate the variation in building climate zones for 3 different decades using an ensemble and quartile statistics to provide an assessment of potential building climate zone changes relative to the uncertainties demonstrated using MERRA. Findings and Conclusions: These results show that there is a statistically significant increase in the area covered by warmer climate zones and a tendency for a reduction of area in colder climate zones in some limited regions. The CMIP analysis shows that models vary from relatively little building climate zone change for the least sensitive and conservation assumptions to a warming of at most 3 zones for certain areas, particularly the north central US by the end of the 21st century.

Stackhouse, Paul W., Jr.↗

NASA NTRS2017

Mitigating Satellite-Based Fire Sampling Limitations in Deriving Biomass Burning Emission Rates: Application to WRF-Chem Model over the Northern sub-Saharan African Region

Largely used in several independent estimates of fire emissions, fire products based on MODIS sensors aboard the Terra and Aqua polar-orbiting satellites have a number of inherent limitations, including (a) inability to detect fires below clouds, (b) significant decrease of detection sensitivity at the edge of scan where pixel sizes are much larger than at nadir, and (c) gaps between adjacent swaths in tropical regions. To remedy these limitations, an empirical method is developed here and applied to correct fire emission estimates based on MODIS pixel level fire radiative power measurements and emission coefficients from the Fire Energetics and Emissions Research (FEER) biomass burning emission inventory. The analysis was performed for January 2010 over the northern sub-Saharan African region. Simulations from WRF-Chem model using original and adjusted emissions are compared with the aerosol optical depth (AOD) products from MODIS and AERONET as well as aerosol vertical profile from CALIOP data. The comparison confirmed an 30-50% improvement in the model simulation performance (in terms of correlation, bias, and spatial pattern of AOD with respect to observations) by the adjusted emissions that not only increases the original emission amount by a factor of two but also results in the spatially continuous estimates of instantaneous fire emissions at daily time scales. Such improvement cannot be achieved by simply scaling the original emission across the study domain. Even with this improvement, a factor of two underestimations still exists in the modeled AOD, which is within the current global fire emissions uncertainty envelope. Plain Language Summary Polar-orbiting satellites sensors, such as MODIS, have limitations in detecting fires under clouds or when viewing angles are large or in the gaps among satellites' different ground swaths. Here we developed an empirical method to mitigate the effect of these limitations in fire emission estimate. The method is applied to a fire emission inventory (FEER) based on MODIS. We show that, with our method, the adjusted emission inventory improves WRF-Chem simulation of smoke transport and distribution.

emissions properties↗

NASA NTRS2022

A Numerical Study of the Liquid Motion in Titan’s Subsurface Ocean

An ocean filled with liquid water lies beneath the icy surface of several Jovian and Saturnian moons. In such an ocean, the currents are driven by various phenomena such as the tidal forcing, the deformation of the ice shell lying at its top, the temperature gradient resulting from the surface and bottom heat fluxes…The flow induced by the first two forcings can be modelled by means of a 2D depth-averaged model, while the third one generates horizontal and vertical density variations whose effects can only be captured by a 3D baroclinic model. We study the tides of Titan’s subsurface ocean and the impact of the ice shell on the liquid motion by means of the Second-generation Louvain-la-Neuve Ice–ocean Model, SLIM (https://www.climate.be/slim). The impact of the ice shell lying at the top of the ocean is modelled by a surface friction term and surface pressure terms. The latter are a function of the difference between the ocean elevation and the vertical displacement of the shell and the time derivative of this difference. Because of Titan’s appreciable obliquity (0.306°), the tidal motion expected (and found) is similar to the Europa tidal scenario described by Tyler (2008): the surface elevation consists of two bulges rotating around Titan and the associated depth-averaged velocity field consists of two gyres, separated by an area of high speed flow, whose centre follows a sinusoidal path centred on the equator. The ice shell damps the surface motion, thus slowing down the flow, without significantly modifying the spatial patterns of these fields. The depth of the ocean and the mechanical characteristics of the ice shell being poorly constrained, a sensitivity analysis is conducted. The depth-averaged flow slows down when the depth is increased and a lag appears in the tidal phase but the tidal range remains similar. The ice shell mechanical characteristics influences both the elevation and depth-averaged velocity fields in terms of magnitude but does not modify the spatial patterns of these fields. The influence of the surface heat flux is studied by means of the 3D baroclinic version of SLIM. The heat flux derived from Titan’s topography by Kvorka et al. (2018) is used as surface boundary condition for the temperature equation while a uniform bottom heat flux is implemented. Its value is computed assuming that the heat budget of the ocean is at equilibrium. These boundary conditions cause density variations, which impact the hydrodynamics of the ocean. While the flow velocity induced by these variations is two orders of magnitude smaller than the tidal flow, its orientation is time-independent, hence impacting the orientation of the velocity field. Although the variations of ocean surface elevation and speed with respect to the shell mechanical properties can be larger than those induced by the surface heat flux, taking into account the latter results in large variations of the velocity field global patterns, which was not observed when modifying the shell mechanical properties. Future studies should therefore focus on modelling the surface and bottom heat fluxes while uncertainties about the mechanical characteristics of the shell can be tolerated.

Titan↗

NASA NTRS2009

Alternative Representations of Convective Processes in the NASA GEOS-5 AGCM

The gap in explicit resolution of phenomena between global climate models and cloud resolving models is shrinking at a steady pace with global integrations of several km in resolution now practical for at least some time scales. In moving toward finer resolution the long standing problem of convective parameterization is being examined along with the assumption of convective quasi-equilibrium and how this can be reconciled with the stochastic and intermittent nature of convection. In this context we examine the nature of parameterized convection in the NASA Goddard Earth Observing System (GEOS-5) Atmospheric General Circulation Model. Our analysis uses both coarse (2.5 degree) and fine scale 0.25 degree spatial resolution integrations. Two basic formulations are compared: the default option is the Relaxed Arakawa-Schubert (RAS) scheme which invokes a sequence of linearly entraining plumes and quasi-equilibrium closure. An optional modification of this method (the "Stochastic Tokioka" constraint) places a random lower limit on plume entrainment. An alternative representation is the Kain-Fritsch parameterization which was originally developed for mesoscale numerical modeling strategies. Here entrainment is determined by a crude buoyancy sorting approach that allows the plume spectrum to be more responsive to ambient vertical stratification of moisture. Diagnostics of the model behavior are referenced to recent observational evidence of continuous phase transition behavior. In particular we examine the relationship between column water vapor and probablility of convective presence and intensity. Sensitivity of the statistics of convective behavior to parcel mixing/entrainment formulations and parcel initial thermodynamics are considered. Observational statistics from A-Train and TRMM sensors provide validation of the model integrations.

Robertson, Franklin↗

NASA NTRS2009

Climate and Lightning: An updated TRMM-LIS Analysis

The TRMM Lightning Imaging Sensor (LIS) has sampled global tropical and sub-tropical lightning flash densities for approximately 11 years. These data were originally processed and results presented by the authors in the 3rd AMS MALD Conference held in 2007 using both pre and post TRMM-boost lightning data. These data were normalized for the orbit boost by scaling the pre-boost data by a fixed constant based on the different swath areas for the pre and post-boost years (post-boost after 2001). Inevitably, one must question this simple approach to accounting for the orbit boost when sampling such a noisy quantity. Hence we are in the process of reprocessing the entire 11-year TRMM LIS dataset to reduce the orbit swath of the post-boost era to that of the pre-boost in order to eliminate sampling bias in the dataset. Study of the diurnal/seasonal/annual sampling suggests that those biases are already minimal and should not contribute to error in examination of annual trends. We will present new analysis of the 11-year annual trends in total lightning flash density for all latitudinal belts and select regions/regimes of the tropics as related to conventional climate signals and precipitation contents in the same period. The results should enable us to address, in some fashion, the sensitivity of the lightning flash density to subtle changes in climate.

Petersen, Walter A.↗

NASA NTRS2013

Inter-Calibration and Concatenation of Climate Quality Infrared Cloudy Radiances from Multiple Instruments

A change in climate is not likely captured from any single instrument, since no single instrument can span decades of time. Therefore, to detect signals of global climate change, observations from many instruments on different platforms have to be concatenated. This requires careful and detailed consideration of instrumental differences such as footprint size, diurnal cycle of observations, and relative biases in the spectral brightness temperatures. Furthermore, a common basic assumption is that the data quality is independent of the observed scene and therefore can be determined using clear scene data. However, as will be demonstrated, this is not necessarily a valid assumption as the globe is mostly cloudy. In this study we highlight challenges in inter-calibration and concatenation of infrared radiances from multiple instruments by focusing on the analysis of deep convective or anvil clouds. TRMM/VIRS is potentially useful instrument to make correction for observational differences in the local time and foot print sizes, and thus could be applied retroactively to vintage instruments such as AIRS, IASI, IRIS, AVHRR, and HIRS. As the first step, in this study, we investigate and discuss to what extent AIRS and VIRS agree in capturing deep cloudy radiances at the same local time. The analysis also includes comparisons with one year observations from CrIS. It was found that the instruments show calibration differences of about 1K under deep cloudy scenes that can vary as a function of land type and local time of observation. The sensitivity of footprint size, view angle, and spectral band-pass differences cannot fully explain the observed differences. The observed discrepancies can be considered as a measure of the magnitude of issues which will arise in the comparison of legacy data with current data.

cloudy radiance↗

NASA NTRS2011

Passive Microwave Remote Sensing of Falling Snow and Associated GPM Field Campaigns

Retrievals of falling snow from space represent one of the next important challenges for the atmospheric, hydrological, and energy budget scientific communities. Historically, retrievals of falling snow have been difficult due to the relative insensitivity of satellite rain-based channels as used in the past. We emphasize the use of high frequency passive microwave channels (85-200 GHz) since these are more sensitive to the ice in clouds and have been used to estimate falling snow from space. While satellite-based remote sensing provides global coverage of falling snow events and the science is relatively new, retrievals are still undergoing development with challenges remaining. There are several current satellite sensors, though not specifically designed for estimating falling snow, are capable of measuring snow from space. These include NOAA's AMSU-B, the MHS sensors, and CloudSat radar. They use high frequency (greater than 85 GHz) passive and active microwave and millimeter-wave channels that are sensitive to the scattering from ice and snow particles in the atmosphere. Sensors with water vapor channels near 183 GHz center line provide opaqueness to the Earth's surface features that can contaminate the falling snow signatures, especially over snow covered surface. In addition, the Global Precipitation Measurement (GPM) mission scheduled for launch in 2013 is specifically designed to measure both liquid rain and frozen snow precipitation. Since falling snow from space is the next precipitation measurement challenge from space, information must be determined in order to guide retrieval algorithm development for these current and future missions. This information includes thresholds of detection for various sensor channel configurations, snow event system characteristics, and surface types. For example, can a lake effect snow system with low cloud tops having an ice water content (IWC) at the surface of 1.0 gram per cubic meter be detected? If this information is known, we can focus retrieval efforts on detectable storms and concentrate advances on achievable results. In this work, the focus is to determine thresholds of detection for falling snow for various snow conditions over land and lake surfaces. The results rely on simulated Weather Research Forecasting (WRF) simulations of falling snow cases [9] since simulations provide all the information to determine the measurements from space and the ground truth. This analysis relies on data from the Canadian CloudSat/CALIPSO Validation Program (C3VP) field campaign held from October 31, 2006 through March 1, 2007. In January 2012 GPM will return to the C3VP area for the GPM Cold Precipitation Experiment (GCPEx). This presentation will describe the thresholds-of-detection procedure and results, as well as the field campaign details.

Skofronick-Jackson, Gail↗

NASA NTRS2013

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↗