Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “parameterization”

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 901 records · Page 50

Synergy of Satellite Radiation, Precipitation, and Other Meteorological Variable Observations for Global Mean Sea Surface Turbulent Heat Flux Estimation

Sea surface turbulent heat flux is one of the key components in global freshwater and energy balances and plays an important role in atmospheric dynamics, thermodynamics and general circulation. This turbulent heat flux is dominantly decided by sea surface latent heat release with some contribution from sensible heat exchange. The flux and its anomaly could significantly affect ocean heat storage and ocean circulation. They are part of and have great impacts on climate variability. Though the turbulent heat flux is extremely important for the climate and weather systems, only very limited ship and buoy observations of the flux are available over open oceans. There are no global, operational, direct measurements of this crucial meteorological variable. Global estimates are, basically, indirectly calculated from bulk turbulent flux parameterization with a combination of satellite column water vapor, sea surface water temperature, and wind speed observations. Critical parameters such as sea surface air temperature and humidity are estimated from empirical relations of theses variables with the water vapor and sea surface water temperature, respectively. Lacking accurate knowledge on surface air temperature and humidity, uncertainties in the parameterization and potential changes in the non-linear turbulent processes with long-term climate variations could cause large errors in estimated long-term turbulent fluxes from the indirect method as shown in current global sea surface turbulent flux datasets. This study uses synergized data of satellite global radiation, precipitation, and other meteorological variable observations to estimate sea surface turbulent heat fluxes. The radiation observations are made by the satellite Clouds and the Earth’s Radiant Energy System (CERES) sensors, while the global precipitation data is from NASA’s satellite Global Precipitation Climatology Project (GPCP). Other data includes satellite sea surface water temperature and wind speed observations. These datasets are obtained from a wide range of space sensors from passive to active instruments and from visible and near infrared to thermal infrared and microwave spectral sounders. One significant feature of these datasets is that they have multi-decades long climate records. Top-of-atmosphere (TOA) radiation and its anomaly represents the net heat energy input to the climate system, and the oceanic precipitation and ocean-land moisture transport can be used to quantify sea surface latent heat release. Based on the synergized datasets and the principle of global water and energy balances, global mean turbulent heat fluxes are estimated. The turbulent heat anomalies for the first two decades of the 21st century are, then, obtained mainly from global CERES radiation and GPCP precipitation anomalies, along with CERES derived ocean-land heat transports. Analysis indicates that the uncertainties in the estimated turbulent flux anomalies may be reduced considerably. The results suggest a strong needs in synchronized and synergized observations of atmospheric radiation, precipitation, and oceanic meteorological variables for long-term climate studies.

Bing Lin↗

Accurate Assessment of Land-Atmosphere Coupling in Climate Models Requires High Frequency Data Output

Land-atmosphere (L-A) interactions are important for understanding convective processes, climate feedbacks, the development and perpetuation of droughts, heatwaves, pluvials, and other land-centred climate anomalies. Local L-A coupling (LoCo) metrics capture relevant L-A processes, highlighting the impact of soil and vegetation states on surface flux partitioning, and the impact of surface fluxes on boundary layer (BL) growth, development, and entrainment of air above the BL. A primary goal of the Climate Process Team on Coupling Land and Atmospheric Subgrid Parameterizations (CLASP) is parameterizing and characterizing the impact of subgrid heterogeneity in global and regional earth system models (ESMs) to improve the connection between land and atmospheric states and processes. A critical step in achieving that aim is the incorporation of L-A metrics, especially LoCo metrics, into climate model diagnostic process streams. However, because land-atmosphere interactions span time scales of minutes (e.g., turbulent fluxes), hours (e.g., BL growth and decay), days (e.g., soil moisture memory), and seasons (e.g., variability of behavioural regimes between soil moisture and latent heat flux), with multiple processes of interest happening in different geographic regions at different times of year, there is not a single metric that captures all the modes, means, and methods of interaction between the land and the atmosphere. And while monthly means of most of the LoCo-relevant variables are routinely saved from ESM simulations, data storage constraints typically preclude routine archival of the hourly data that would enable the calculation of all LoCo metrics. Here we outline a reasonable data request that would allow for adequate characterization of sub-daily coupling processes between the land and the atmosphere, preserving enough sub-daily output to describe, analyse, and better understand L-A coupling in modern climate models. A secondary request involves embedding calculations within the models to determine mean properties in and above the BL to further improve characterization of model behaviour. Higher-frequency model output will (i) allow for more direct comparison with observational field campaigns on process-relevant time scales, (ii) enable demonstration of inter-model spread in L-A coupling processes, and (iii) aid in targeted identification of sources of deficiencies and opportunities for improvement of the models.

Kirsten L. Findell↗

GMAO OSSE Framework in Support of PBL Mission Science

The planetary boundary layer is the bottom layer of the troposphere where most of human activities take place. Atmospheric pollutants are capped by the temperature inversion layer in PBL. The thickness of PBL, ranging from tens of meter to several kilometer, affects pollutant dispersion and air quality. The thickness of PBL typically shows a strong diurnal cycle but it is difficult to model and predict its change because complicated dynamic and thermodynamic processes involved with air-surface exchange of temperature and moisture, convective mixing, surface friction, topography, advection and radiation heating and cooling influence the PBL height. Various parameterization schemes for PBL were developed but the underlying processes of PBL are neither clearly understood nor represented in NWP models and adding large uncertainty into weather and climate predictions. There are no systematic global observations to provide information on thermodynamic structure of PBL and efforts to explore new spaceborne instruments and measurement techniques are growing. This study aims to develop a PBL OSSE framework leveraging existing GMAO’s OSSE system that utilizes GEOS data assimilation and global forecast model (1) to evaluate existing and potential new observation types for the PBL structure analysis and prediction and (2) to test sensitivity of PBL parameterization schemes to PBL forecasts.

Min-Jeong Kim↗

Boundary-Layer Cloud Modeling Challenges on the North Slope of Alaska

The accurate modeling and prediction of cloud base heights is critical for energy balance calculations and aviation operations, alike. Low-level (i.e., boundary-layer) Arctic clouds can be difficult to model, making prediction of formation and dissipation challenging. Primarily mixed-phase, these clouds typically contain low quantities of supercooled liquid water and often slowly precipitate relatively small amounts of moderately and heavily rimed snow particles. While this appears to be the predominant cloudy state on the North Slope of Alaska (NSA), the delicate balance of microphysical, dynamical, radiative, surface coupling, and advective processes can rapidly shift to heavy snow (with various degrees of riming) or to a complete dissipation of the cloud layer without any precipitation, depending on the dominant processes. Here we strive to disentangle these various processes. First, we compare the predictive performances of four different numerical weather models in forecasting the presence and base-heights of low-level clouds: the High-Resolution Rapid Refresh - Alaska (HRRR-AK) model, the Polar Weather Research and Forecasting (Polar WRF) model, the Unified Model (UM), and the European Centre for Medium-range Weather Forecasting (ECMWF) model. Initial results comparing model output at two U.S. Department of Energy Atmospheric Radiation Measurement (AMT) NSA sites, during the fall season in 2019 and 2022, show that the UM slightly outperforms the HRRR-AK in terms of accurately forecasting the presence of a low-level cloud layer (89% of the time). All models have a significant bias of 300 to 800 meters in forecasting cloud base height (lower than is observed); however, the UM and ECMWF models have the lowest biases. Finally, a case study for a particularly challenging April 2017 thin-cloud event is presented, wherein we compare the performance of four different bulk microphysical parameterization schemes using a higher-resolution large eddy simulation (LES) model, the WRF-LES. Initial results show that the Thompson scheme was the only one able to reproduce and sustain a substantial supercooled liquid layer, but it was unable to reproduce the transition from a deep, liquid-rich cloud to a thin layer with moderately and heavily rimed precipitation. This is the first step in linking simulated LES-scale riming processes with those parameterized at a coarser mesoscale model scale. This has important implications for forecasting low-level clouds in an operational environment, given the efficiency of the riming process.

cloud base heights↗

A new spectrally-invariant approach to the remote sensing of inhomogeneous clouds

Current operational satellite retrievals of cloud optical and microphysical properties go back to the Nakajima-King technique developed at the end of the 1980s. This technique is based on library calculations for plane-parallel homogeneous clouds. It often works well for overcast skies but leads to substantial errors for inhomogeneous and broken cloud fields where the plane-parallel-geometry assumption is no longer valid. The basic concept of a new technique was first introduced in nuclear physics to quantify the critical conditions of reactors. Based on the eigenvalues of the radiative transfer equation, their approach provides a powerful means to parameterize the structure of 3D media. This parameterization was later successfully applied to relate surface reflectance spectra to 3D canopy structure and is known as the spectrally-invariant approximation. The proposed approach adapts this technique to the remote sensing of cloud properties such as droplet single scattering albedo and the average number of scatterings (which are the fundamental parameters in radiative transfer theory), with an emphasis on quantifying the associated errors and uncertainties. This retrieval is free from the plane-parallel homogeneous cloud assumption.

clouds↗

The Sensitivity of the Equatorial Pacific ODZ to Particulate Organic Matter Remineralization in A Climate Model Under Pre-Industrial Conditions

Marine oxygen plays a fundamental role in regulating the transfer of organic carbon and nutrients to their dissolved inorganic forms, serving as the terminal electron acceptor for heterotrophic respiration. Oxygen can become limiting to these processes in coastal and open-ocean oxygen deficient zones (ODZs). The maintenance of ODZs depends on the balance between physical processes such as ventilation and biogeochemical processes such as remineralization. These two processes act in opposing directions on ODZs, with ventilation responsible for transporting oxygen from the surface, where it is near saturation with the partial pressure of O2 in the atmosphere, to deeper oxygen-depleted layers, and remineralization acting to consume oxygen through biogeochemical processes such as microbial degradation throughout the water column. Remineralization is represented in all CMIP6 models, but the actual parameterizations, as well as their magnitude, are widely varying. In this study, we examine the sensitivity to remineralization of the equatorial Pacific ODZ O2 ol/kg) in a model; the NASA GISS Ocean Biogeochemical Model (NOBM-G) embedded in the NASA-GISS coupled atmosphere-ocean model. We find that increasing the remineralization rate shoals the ODZ onset depth (i.e. the regionally averaged upper bound of the ODZ), and decreases ODZ volume and regional-averaged thickness. Changes in biological consumption and vertical convergence of oxygen are identified as the primary contributors to changes in ODZ onset depth. ODZ thickness is mainly influenced by the shoaling of the bottom boundary of the ODZ, which decreases with maximum remineralization rate due to decreasing oxygen consumption and increasing horizontal oxygen convergence in deep waters. On the other hand, vertically-averaged ODZ area has a more complex, non-monotonic relationship with maximum remineralization rate. While these findings suggest an important role for remineralization in determining the shape of the ODZ in our model, the relative importance of remineralization vs. other physical parameterizations remains to be established. While our results reflect the structure of our ocean biogeochemical model, the relationship of remineralization to ODZ characteristics in other models should be examined to better inform model inter-comparisons.

Oxygen-deficient zone↗

Version 8 of the SMAP Level-4 Soil Moisture Data Assimilation Product

The NASA Soil Moisture Active Passive (SMAP) mission Level-4 Soil Moisture (L4_SM) product provides global, 9-km resolution, 3-hourly surface (0-5 cm) and root-zone (0-100 cm) soil moisture from April 2015 to present with a mean latency of 2.5 days from the time of observation. The product is based on the assimilation of SMAP L-band (1.4 GHz) brightness temperature (Tb) observations into the NASA Catchment land surface model as the model is driven with observations-based precipitation forcing. This presentation discusses the improvements in the forthcoming Version 8 of L4_SM, including updates in the precipitation forcing, the Catchment model parameters, and the L-band microwave radiative transfer model (mwRTM). The precipitation observations used in L4_SM Version 8 outside of North America and the high latitudes are from the latest (Version 7) NASA Integrated Multi-satellitE Retrievals for the Global Precipitation Measurement mission (IMERG) products. Moreover, occasionally excessive precipitation rates in earlier versions of L4_SM along certain longitudes in North America were eliminated by a bug fix in the precipitation corrections algorithm. The Catchment model in L4_SM Version 8 uses climatological snow albedo values based on observations from the Moderate Resolution Imaging Spectroradiometer, replacing the look-up table parameterization of earlier versions. Additionally, corrected soil parameters were implemented for a small region in Argentina that had erroneously been classified as peat because of an error in the ancillary soil data. Finally, the mwRTM in L4_SM Version 8 uses the Mironov soil mixing approach and updated values of the L-band scattering albedo, soil roughness, and vegetation opacity climatology obtained from the latest (Release 19) SMAP Level-2 dual-channel soil moisture retrieval product. During the development of L4_SM Version 8, the change in the mwRTM parameterization resulted in a reduced unbiased RMSE of surface soil moisture when verified against in situ measurements. It also reduced the standard deviation of the Tb observation-minus-forecast residuals by ~0.15 K, highlighting the importance of the mwRTM for successful data assimilation. The bug fixes in the precipitation corrections algorithm and the Catchment model soil parameters in the Argentina region result in locally large improvements of the simulated land surface states. In summary, the ongoing refinements of the L4_SM product continue to improve its science quality and performance for global soil moisture monitoring

Rolf Reichle↗

A Neural Network Parametrization of Volumetric Cloud Fraction Profiles Using Satellite Observations and MERRA-2 Reanalysis Meteorological Data

Clouds play a crucial role in regulating the hydrologic cycle and Earth's radiative energy budget, yet they are often poorly represented in global climate models (GCMs). This study applies deep machine learning techniques to develop a physical parameterization of volumetric cloud fraction (VCF), the fraction of a 3-D grid volume occupied by clouds using satellite lidar-radar measurements. The neural network (NN) captures the complicated relationships between observed VCF profiles and collocated meteorological variables from MERRA-2 reanalysis data. Our results show that the NN model, particularly a sequence-to-sequence long short-term memory (LSTM) network with a sixfactor loss function, effectively learns the underlying cloud physical processes. The NN model outperforms MERRA-2 reanalysis in representing low-level clouds in tropical and subtropical regions and low- and middle-level clouds over midlatitude storm-track regions, and also improves VCF histograms. These improvements are reflected in the vertical distributions of zonally, meridionally, and globally averaged VCFs, geographic distributions of low-, middle-, and high-level clouds, and seasonal variations in monthly-mean VCF. Furthermore, the NN predictions effectively capture the El Niño-Southern Oscillation (ENSO) effects and other interannual variations. The NN parameterization is further evaluated through a sensitivity analysis, in which a single predictor is perturbed at a time. This reveals that relative humidity (RH) is the dominant factor influencing variations in globally averaged VCF at low and middle altitudes, followed by temperature. At higher altitudes, temperature becomes the primary driver of VCF through its effect on RH. Changes in wind components had minimal impact on globally averaged VCF.

Shan Zeng↗

A Single-column Model Ensemble Approach Applied to the TWP-ICE Experiment

Single-column models (SCM) are useful test beds for investigating the parameterization schemes of numerical weather prediction and climate models. The usefulness of SCM simulations are limited, however, by the accuracy of the best estimate large-scale observations prescribed. Errors estimating the observations will result in uncertainty in modeled simulations. One method to address the modeled uncertainty is to simulate an ensemble where the ensemble members span observational uncertainty. This study first derives an ensemble of large-scale data for the Tropical Warm Pool International Cloud Experiment (TWP-ICE) based on an estimate of a possible source of error in the best estimate product. These data are then used to carry out simulations with 11 SCM and two cloud-resolving models (CRM). Best estimate simulations are also performed. All models show that moisture-related variables are close to observations and there are limited differences between the best estimate and ensemble mean values. The models, however, show different sensitivities to changes in the forcing particularly when weakly forced. The ensemble simulations highlight important differences in the surface evaporation term of the moisture budget between the SCM and CRM. Differences are also apparent between the models in the ensemble mean vertical structure of cloud variables, while for each model, cloud properties are relatively insensitive to forcing. The ensemble is further used to investigate cloud variables and precipitation and identifies differences between CRM and SCM particularly for relationships involving ice. This study highlights the additional analysis that can be performed using ensemble simulations and hence enables a more complete model investigation compared to using the more traditional single best estimate simulation only.

numerical weather forecasting↗

Modeling and Evaluation of the Global Sea-Salt Aerosol Distribution: Sensitivity to Emission Schemes and Resolution Effects at Coastal/Orographic Sites

One of the major sources of uncertainty in model estimates of the global sea-salt aerosol distribution is the emission parameterization. We evaluate a new sea-salt aerosol life cycle module coupled to the online multi-scale chemical transport model NMMB/BSC-CTM. We compare 5 year global simulations using five state-of-the-art sea-salt open-ocean emission schemes with monthly averaged coarse aerosol optical depth (AOD) from selected AERONET sun photometers, surface concentration measurements from the University of Miami's Ocean Aerosol Network, and measurements from two NOAA/PMEL cruises (AEROINDOEX and ACE1). Model results are highly sensitive to the introduction of sea-surface-temperature (SST)-dependent emissions and to the accounting of spume particles production. Emission ranges from 3888 teragrams per year to 8114 teragrams per year, lifetime varies between 7.3 hours and 11.3 hours, and the average column mass load is between 5.0 teragrams and 7.2 teragrams. Coarse AOD is reproduced with an overall correlation of around 0.5 and with normalized biases ranging from +8.8 percent to +38.8 percent. Surface concentration is simulated with normalized biases ranging from minus 9.5 percent to plus 28 percent and the overall correlation is around 0.5. Our results indicate that SST-dependent emission schemes improve the overall model performance in reproducing surface concentrations. On the other hand, they lead to an overestimation of the coarse AOD at tropical latitudes, although it may be affected by uncertainties in the comparison due to the use of all-sky model AOD, the treatment of water uptake, deposition and optical properties in the model and/or an inaccurate size distribution at emission.

Bias↗

Ionospheric Mapping with the GPS/MET

The GPS/MET experiment, which placed a GPS receiver in a low-Earth orbit tracking the GPS in an occultation geometry, has collected thousands of occultations since its launch in April of 1995. Each occultation can be inverted to give an electron density profile in the ionosphere, temperature and pressure profiles in the lower mesosphere, stratosphere and upper troposphere, and water vapor density profiles in the lower troposphere. This paper gives a summary of the ionospheric effects on the GPS/MET signal and examines some of the retrieved electron density profiles. We examine the bending induced by the ionosphere on the occulting signal and the resulting separation of the two GPS links corresponding to the L1 and L2 phase signals. We also examine the amplitude scintillation caused by sharp layers at the bottom of the ionosphere. We briefly describe the Abel inversion method for obtaining the index of refraction and show several examples of electron density profiles and comparisons derived from the Parameterized Ionospheric Model (PIM) and to incoherent scatter radar measurements.

GPS↗

A Fast Vector Radiative Transfer Model for the Atmosphere-Ocean Coupled System

To infer atmospheric and oceanic constituent properties from polarimetric observations, an efficient and accurate retrieval algorithm is desirable. In-line radiative transfer calculations are indispensable if a large state vector, including both atmospheric profiles and surface properties, is used to improve retrieval accuracy. However, in-line radiative transfer calculations are usually not computationally efficient for remote sensing applications. Therefore, there is a pressing need to develop an accurate and fast vector radiative transfer model (RTM) to fully utilize satellite polarimetric observations. This paper reports on a fast vector RTM, referred to as TAMU-VRTM, in support of polarimetric remote sensing, which is capable of simulating the Stokes vector values observed at the top of the atmosphere and at the surface by fully considering absorption, scattering, and emission in the atmosphere and ocean. Gaseous absorption is parameterized with respect to gas concentration, temperature, and pressure, by using a regression method applicable to an inhomogeneous atmospheric path. An efficient two-component approach combining the small angle approximation and the adding-doubling method is utilized to solve the vector radiative transfer equation (RTE). The thermal emission component of the RTE solution is obtained by an efficient doubling process. The air-sea interface is treated as a wind-ruffled rough surface in the model to mimic a realistic ocean surface. Several oceanic optical property models are introduced to model ocean inherent optical properties. To demonstrate the applicability of the TAMUVRTM, simulations are compared with satellite observations, and results from other vector radiative transfer methods including benchmarks.

radiative transfer in coupled atmosphere-ocean sys↗

A model-independent data assimilation (MIDA) module and its applications in ecology

Models are an important tool to predict Earth system dynamics. An accurate prediction of future states of ecosystems depends on not only model structures but also parameterizations. Model parameters can be constrained by data assimilation. However, applications of data assimilation to ecology are restricted by highly technical requirements such as model-dependent coding. To alleviate this technical burden, we developed a model-independent data assimilation (MIDA) module. MIDA works in three steps including data preparation, execution of data assimilation, and visualization. The first step prepares prior ranges of parameter values, a defined number of iterations, and directory paths to access files of observations and models. The execution step calibrates parameter values to best fit the observations and estimates the parameter posterior distributions. The final step automatically visualizes the calibration performance and posterior distributions. MIDA is model independent, and modelers can use MIDA for an accurate and efficient data assimilation in a simple and interactive way without modification of their original models. We applied MIDA to four types of ecological models: the data assimilation linked ecosystem carbon (DALEC) model, a surrogate-based energy exascale earth system model: the land component (ELM), nine phenological models and a stand-alone biome ecological strategy simulator (BiomeE). The applications indicate that MIDA can effectively solve data assimilation problems for different ecological models. Additionally, the easy implementation and model-independent feature of MIDA breaks the technical barrier of applications of data–model fusion in ecology. MIDA facilitates the assimilation of various observations into models for uncertainty reduction in ecological modeling and forecasting.

Earth system dynamics↗