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 721 records · Page 40

All sources NASA NTRS DOE OSTI

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.↗

DOE OSTI2025

Standardising the “Gregory method” for calculating equilibrium climate sensitivity

The equilibrium climate sensitivity (ECS) – the equilibrium global mean temperature response to a doubling of atmospheric CO 2 – is a high-profile metric for quantifying the Earth system's response to human-induced climate change. A widely applied approach to estimating the ECS is the “Gregory method” (Gregory et al., 2004), which uses an ordinary least squares (OLS) regression between the net radiative flux, N, and surface air temperature anomalies, ΔT, from a 150 year experiment in which atmospheric CO 2 concentrations are quadrupled. The ECS is determined by extrapolating the linear fit to N=0, i.e. the ΔT-intercept, indicating the point at which the system is back in equilibrium. This method has been used to compare ECS estimates across the CMIP5 and CMIP6 ensembles and will likely be a key diagnostic for CMIP7. Despite its widespread application, there is little consistency or transparency between studies in how the climate model data is processed prior to the regression, leading to potential discrepancies in ECS estimates. We identify 32 alternative data processing pathways, varying by differences in global mean weighting, net radiative flux variable, anomaly calculation method, and linear regression fit. Using 44 CMIP6 models, we systematically assess the impact of these choices on ECS estimates and calculate uncertainty ranges using two bootstrap approaches. While the inter-model ECS range is insensitive to the data processing pathway, individual outlier models exhibit notable differences. Approximating a model's native grid cell area (if irregular) with cosine of the latitude can decrease the ECS by 11 %, the choice of N-variable can change the ECS by 6 %, and some anomaly calculation methods can introduce spurious temporal correlations in the processed data. Beyond data processing choices, we also evaluate an alternative linear regression method – total least squares (TLS) – which has a more statistically robust basis than OLS. However, for consistency with previous literature, and given TLS may reduce the ECS compared to OLS (by up to 24 %), thereby making a known bias in the Gregory method worse, we do not feel there is sufficient clarity to recommend a transition to TLS in all cases. To improve reproducibility and comparability in future studies, we recommend a standardised Gregory method: weighting the global mean by cell area, using the top of the atmosphere (as opposed to the top of model) N-variable, and calculating anomalies by first applying a rolling average to the preindustrial control timeseries then subtracting from the raw CO 2 quadrupling experiment. This approach accounts for model drift while reducing noise in the data to best meet the pre-conditions of the linear regression. While CMIP6 results of the multi-model mean ECS appear insensitive to these processing choices, similar assumptions may not hold for CMIP7, underscoring the need for standardised data preparation in future climate sensitivity assessments.

Geosciences↗

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↗

DOE OSTI2023

Volcanic Imprints in Last-Millennium Land Summer Temperatures in the Circum-North Atlantic Area

Summer cooling is one of the most direct consequences of explosive volcanic eruptions that can affect ecosystems and human societies. Recent studies revealed a multiyear cooling impact on hemispheric and global summer temperatures after tropical eruptions, yet the volcanic responses appear to vary on regional scales. Here, in this work, we revisit volcano-induced summer cooling in eastern Canada and northern and central Europe by applying superposed epoch analysis on CMIP6-PMIP4 simulations and millennial temperature reconstructions based on tree-ring density. We then examine potential causes modulating region-specific volcanic impact. While confirming that, on average, tropical eruptions over the last millennium have induced a longer cooling (>4 yr) than eruptions from extratropical Northern Hemisphere in all three North Atlantic regions, we show that the peak magnitude of cooling is stronger in eastern Canada. We also find that the detected volcanic temperature anomalies can be strongly affected by the selection and number of volcanic events and nonvolcanic signals embedded in the climate time series. This study highlights the risks of using highly noisy proxy records to investigate volcanic impacts, especially in regions with strong unforced climate variability. The CMIP6-PMIP4 simulations generally agree with the three reconstructions on the average response to tropical eruptions, but their performance is poorer regarding the production of significant cooling after extratropical eruptions. Our results further suggest that the particular sensitivity to tropical eruptions in eastern Canada is likely related to increased sea ice surrounding Quebec–Labrador associated with the positive Arctic Oscillation and North Atlantic Oscillation formed during the first posteruption winter.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Evaluation of aerosol–cloud interactions in E3SM using a Lagrangian framework

A Lagrangian framework is used to evaluate aerosol–cloud interactions in the U.S. Department of Energy's Energy Exascale Earth System Model (E3SM) version 1 (E3SMv1) for measurements taken at Graciosa Island in the Azores where a U.S. Department of Energy Atmosphere Radiation Measurement (ARM) site is located. This framework uses direct measurements of cloud condensation nuclei (CCN) concentration (instead of relying on satellite retrievals of aerosol optical depth) and incorporates a suite of ground-based ARM measurements, satellite retrievals, and meteorological reanalysis products that when applied to over a 1500 trajectories provides key insights into the evolution of low-level clouds and aerosol radiative forcing that is not feasible from a traditional Eulerian analysis framework. Significantly lower concentrations (40 %) of surface CCN concentration are measured when precipitation rates in 48 h back trajectories average above 1.2 mm d -1 in the Integrated Multi-satellitE Retrievals for Global Precipitation Measurement (IMERG) product. The depletion of CCN concentration when precipitation rates are elevated is nearly twice as large in the ARM observations compared to E3SMv1 simulations. The model CCN concentration bias remains significant despite modifying the autoconversion and accretion rates in warm clouds. As the clouds in trajectories associated with larger surface-based CCN concentration advect away from Graciosa Island, they maintain higher values of droplet number concentrations (N d ) over multiple days in observations and E3SM simulations compared to trajectories that start with lower CCN concentrations. The response remains robust even after controlling for meteorological factors such as lower troposphere stability, the degree of cloud coupling with the surface, and island wake effects. E3SMv1 simulates a multi-day aerosol effect on clouds and a Twomey radiative effect that is within 30 % of the ARM and satellite observations. However, the mean cloud droplet concentration is more than 2–3 times larger than in the observations. While Twomey radiative effects are similar amongst autoconversion and accretion sensitivity experiments, the liquid water path and cloud fraction adjustments are positive when using a regression model as opposed to negative when using the present-day minus pre-industrial aerosol emissions approach. This result suggests that tuning the autoconversion and accretion alone is unlikely to produce the desired aerosol susceptibilities in E3SMv1.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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↗

DOE OSTI2022

M2CT-22IN1202096 Design for Carbon Conversion Product Pathways with Nuclear Power Plant Integration

Coal is a globally abundant resource that historically has been used for power generation via combustion. As the power industry replaces coal with cleaner methods of generation, energy-rich coal could be used in other chemical and fuel applications. This study presents a coal utilization option in which coal combustion is replaced with a carbon-free nuclear power plant and the coal is upgraded to valuable products for a variety of markets. Coal is prepared for conversion first by the pyrolysis process, which will optimize solid, liquid, and gaseous products based on the market size and potential product value, maximizing the monetary value of coal. This process is designed using bituminous coal from the Appalachian region as a basis to provide a pathway to preserve or transition coal-related jobs and create new jobs associated with the clean energy transition. Process modeling in the AspenOne Suite will be used to determine each component’s sensitivities, costs, inputs, and outputs. Dispatch modeling in the FORCE toolset will optimize the entire system and calculate the NPV for the refinery lifetime. Advanced and light-water reactors are considerations to supply the heat, steam, and electricity to the process. This paper focuses on the technical and market analysis used to determine the optimal processes and product pathways for the carbon refinery. Product pathways are on activated carbon, formic acid synthesis, and methanol synthesis for further upgrading to marketable chemical and polymer products.

01 COAL, LIGNITE, AND PEAT↗

DOE OSTI2023

Spatio-Temporal Denoising Graph Autoencoders with Data Augmentation for Missing Photovoltaic Data Imputation

The integration of the global Photovoltaic (PV) market with real time data-loggers has enabled large scale PV data analytical pipelines for power forecasting and long-term reliability assessment of PV fleets. Nevertheless, the performance of PV data analysis heavily depends on the quality of PV timeseries data. This paper proposes a novel Spatio-Temporal Denoising Graph Autoencoder (STD-GAE) framework to impute missing PV Power Data. STDGAE exploits temporal correlation, spatial coherence, and value dependencies from domain knowledge to recover missing data. It is empowered by two modules. (1) To cope with sparse yet various scenarios of missing data, STD-GAE incorporates a domain-knowledge aware data augmentation module that creates plausible variations of missing data patterns. This generalizes STD-GAE to robust imputation over different seasons and environment. (2) STD-GAE nontrivially integrates spatiotemporal graph convolution layers (to recover local missing data by observed “neighboring” PV plants) and denoising autoencoder (to recover corrupted data from augmented counterpart) to improve the accuracy of imputation accuracy at PV fleet level. We have evaluated our proposed model on two realworld PV datasets. Experimental results show that STD-GAE can achieve a gain of 43.14% in imputation accuracy and remains less sensitive to missing rate, different seasons, and missing scenarios, compared with state-of-the-art data imputation methods such as MIDA and LRTC-TNN.

Fan, Yangxin↗

DOE OSTI2022

M2CT-22IN1202096 Design for Carbon Conversion Product Pathways with Nuclear Power Plant Integration

01 COAL, LIGNITE, AND PEAT↗

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.↗

DOE OSTI2022

Application of a Novel Heat Pump Model for Estimating Economic Viability and Barriers of Heat Pumps in Dairy Applications in the United States

Heat pumps represent an important opportunity for energy savings and decarbonization. This work investigates the techno-economic performance of high-temperature heat pumps (HTHPs) for use in the U.S. dairy industry. The studied heat pump performs a 50 °C temperature lift on a waste heat stream of cleaning water and applies the upgraded heat stream to a fluid milk pasteurization process. This work involved the creation of a HTHP model that estimated the coefficient of performance (COP), internal rate of return (IRR), net present value (NPV), and payback period (PBP), and emissions saved for a heat pump replacing a natural gas boiler. Capital costs, operations, and maintenance (O&M) cost, heat pump lifetime, electricity prices, natural gas prices, and a cost of carbon were varied to perform a parametric study on the factors affecting the break-even price of HTHPs. The results show that HTHP economics are highly sensitive to COP and energy price environment, and less sensitive to capital and O&M cost variance, leading to a large scatter of positive and negative NPVs based on U.S. location. PBPs demonstrate a defined threshold, based on energy price environment, below which favorable two-to-three-year PBPs predominate. This work is focused on the U.S. dairy industry, but international application in relation to fossil vs. electricity price regimes. Heat pumps have seen wider adoption in regions with a high ratio of fossil energy to electricity prices ($/MMBTU vs. $/kWh). The U.S. has plentiful natural gas resulting in lower fossil energy prices which has reduced heat pump adoption. This paper identifies potential first mover industries for HTHP adoption and their associated price regimes even in regions with lower ratios of fossil energy to electricity prices that exist many places globally.

ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATION,↗

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↗