Estimation of Surface Turbulent Fluxes From Land Surface Moisture and Temperature Via a Variational Data Assimilation Framework
Explore the source record for details and available documents.
SEARCH · Engineering Papers
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.
Explore the source record for details and available documents.
Results are presented from a series of forecast experiments which were conducted to assess the importance of large-scale dynamical processes, diabatic heating, and initial data to the prediction of the President's Day cyclone. The synoptic situation and NMC model forecasts for this case are summarized, and the analysis/forecast system and experiments are described. The GLAS Model forecast from the GLAS analysis at 0000 GMT 18 February is found to have correctly predicted intense coastal cyclogenesis and heavy precipitation. A forecast with surface heat and moisture fluxes eliminated failed to predict any cyclogenesis while a similar forecast with only the surface moisture flux excluded showed weak development. Diabatic heating resulting from oceanic fluxes significantly contributed to the generation of low-level cyclonic vorticity and the intensification and slow rate of movement of an upper level ridge over the western Atlantic.
A model is developed for the estimation of the surface fluxes of momentum, heat, and moisture of the cloud topped marine atmospheric boundary layer by use of satellite remotely sensed parameters. The parameters chosen for the problem are the integrated liquid water content, q sub li, the integrated water vapor content, q sub vi, the cloud top temperature, and either a measure of the 10 meter neutral wind speed or the friction velocity at the surface. Under the assumption of a horizontally homogeneous, well-mixed boundary layer, the model calculates the equivalent potential temperature and total water profiles of the boundary layer along with the boundary layer height from inputs of q sub li, q sub vi, and cloud top temperature. These values, along with the 10m neutral wind speed or friction velocity and the sea surface temperature are then used to estimate the surface fluxes. The development of a scheme to parameterize the integrated water vapor outside of the boundary layer for the cases of cold air outbreak and California coastal stratus is presented.
A series of forecast experiments was conducted to assess the accuracy of the GLAS model, and to determine the importance of large scale dynamical processes and diabatic heating to the cyclogenesis. The GLAS model correctly predicted intense coastal cyclogenesis and heavy precipitation. Repeated without surface heat and moisture fluxes, the model failed to predict any cyclone development. An extended range forecast, a forecast from the NMC analysis interpolated to the GLAS grid, and a forecast from the GLAS analysis with the surface moisture flux excluded predicted weak coastal low development. Diabatic heating resulting from oceanic fluxes significantly contributed to the generation of low level cyclonic vorticity and the intensification and slow rate of movement of an upper level ridge over the western Atlantic. As an upper level short wave trough approached this ridge, diabatic heating associated with the release of latent heat intensified, and the gradient of vorticity, vorticity advection and upper level divergence in advance of the trough were greatly increased, providing strong large scale forcing for the surface cyclogenesis.
The Land Surface, Snow and Soil Moisture Model Intercomparison Project (LS3MIP) is designed to provide a comprehensive assessment of land surface, snow, and soil moisture feedbacks on climate variability and climate change, and to diagnose systematic biases in the land modules of current Earth System Models (ESMs). The solid and liquid water stored at the land surface has a large influence on the regional climate, its variability and predictability, including effects on the energy, water and carbon cycles. Notably, snow and soil moisture affect surface radiation and flux partitioning properties, moisture storage and land surface memory. They both strongly affect atmospheric conditions, in particular surface air temperature and precipitation, but also large-scale circulation patterns. However, models show divergent responses and representations of these feedbacks as well as systematic biases in the underlying processes. LS3MIP will provide the means to quantify the associated uncertainties and better constrain climate change projections, which is of particular interest for highly vulnerable regions (densely populated areas, agricultural regions, the Arctic, semi-arid and other sensitive terrestrial ecosystems).The experiments are subdivided in two components, the first addressing systematic land biases in offline mode (LMIP, building upon the 3rd phase of Global Soil Wetness Project; GSWP3) and the second addressing land feedbacks attributed to soil moisture and snow in an integrated framework (LFMIP, building upon the GLACE-CMIP blueprint).
Convective instability is affected by both surface heat fluxes into near-surface air parcels and the atmospheric profile encountered as these air parcels are lifted aloft. In this study, we investigate the evolution of buoyancy experienced by surface-based parcels, taking interstorm soil moisture drydowns as fundamental time periods. We use profile observations of atmospheric temperature and humidity from ground-based passive remote sensing and in situ observations available at the U.S. Department of Energy Atmospheric Radiation Measurement Southern Great Plains site in Oklahoma, United States, and examine drydown events during the warm seasons from 2017 to 2022. Results reveal new insights on the emergence and growth of convective available potential energy (CAPE) during drydowns. We observe that the peak convective potential appears in the afternoon and grows in the midtroposphere (68% of 66 drydown events), with the largest CAPE values toward the end of the drydown event. We also find that the near-surface parcel contributes more to the evolution of CAPE than the atmospheric environment (86% of drydown events). These findings highlight the importance of surface heat fluxes and moisture availability from surface-based parcels in modulating the evolution of atmospheric conditions favorable for convection, underscoring the interplay between surface processes and atmospheric dynamics in driving convective instability.
In order to understand how the climate responds to variations in forcing, one necessary component is to understand the full distribution of variability of exchanges of heat and moisture between the atmosphere and ocean. Surface heat and moisture fluxes are critical to the generation and decay of many coupled air-sea phenomena. These mechanisms operate across a number of scales and contain contributions from interactions between the anomalous (i.e. non-mean), often extreme-valued, flux components. Satellite-derived estimates of the surface turbulent and radiative heat fluxes provide an opportunity to assess results from modeling systems. Evaluation of only time mean and variability statistics, however only provides limited traceability to processes controlling what are often regime-dependent errors. This work will present an approach to evaluate the representation of the turbulent fluxes at the air-sea interface in the current and evolving Goddard Earth Observing System (GEOS) model. A temperature and moisture vertical profile-based clustering technique is used to identify robust weather regimes, and subsequently intercompare the turbulent fluxes and near-surface parameters within these regimes in both satellite estimates and GEOS-driven data sets. Both model reanalysis (MERRA) and seasonal-to-interannual coupled GEOS model simulations will be evaluated. Particular emphasis is placed on understanding the distribution of the fluxes including extremes, and the representation of near-surface forcing variables directly related to their estimation. Results from these analyses will help identify the existence and source of regime-dependent biases in the GEOS model ocean surface turbulent fluxes. The use of the temperature and moisture profiles for weather-state clustering will be highlighted for its potential broad application to 3-D output typical of model simulations.
NASA funded Mesoscale Environmental Simulations and Operations (MESO), Inc. to develop a version of the Mesoscale Atmospheric Simulation System (MASS). The model has been modified specifically for short-range forecasting in the vicinity of KSC/CCAS. To accomplish this, the model domain has been limited to increase the number of horizontal grid points (and therefore grid resolution) and the model' s treatment of precipitation, radiation, and surface hydrology physics has been enhanced to predict convection forced by local variations in surface heat, moisture fluxes, and cloud shading. The objective of this paper is to (1) provide an overview of MASS including the real-time initialization and configuration for running the data pre-processor and model, and (2) to summarize the preliminary evaluation of the model's forecasts of temperature, moisture, and wind at selected rawinsonde station locations during February 1994 and July 1994. MASS is a hydrostatic, three-dimensional modeling system which includes schemes to represent planetary boundary layer processes, surface energy and moisture budgets, free atmospheric long and short wave radiation, cloud microphysics, and sub-grid scale moist convection.
The planetary boundary layer (PBL) is the lowest layer of the atmosphere in contact with the Earth's surface. It plays an important role in atmospheric circulation and dynamics by influencing surface fluxes of moisture, heat, and momentum. The PBL is characterized by its turbulence structure. Organized convection within the PBL may generate gravity waves within the free troposphere. These waves can propagate vertically, modifying global circulation. Remote sensing of the PBL using a downward-looking lidar is providing new insight into dynamic processes of this important layer and its connection with wave generation. An airborne lidar is an ideal tool to study the PBL. The lidar is installed on the NASA Electra aircraft. It consists of a frequency doubled Nd:YAG laser aligned with a 40 cm diameter telescope. As the laser beam propagates downward from the aircraft toward the surface, it is scattered by the aerosols and molecules in the atmosphere. Generally, there is a sharp gradient of aerosol scattering associated with the PBL top, with high aerosol scattering from within the PBL, and generally very low scattering from the free atmosphere above. It is this gradient in aerosol scattering that is used to visualize the structure of the PBL. A small portion of the scattered radiation is collected by the telescope, filtered, detected by a photomultiplier, and digitized at 100 nsec rate. The horizontal resolution is typically around 13 m, using a 10 Hz laser and a nominal aircraft moves, the resolution of features in the atmosphere below the aircraft is about 15 x 13 m. Results from several studies are presented and discussed.
Hydrological models have been calibrated and validated using catchment streamflows. However, using a point measurement does not guarantee correct spatial distribution of model computed heat fluxes, soil moisture and surface temperatures. With the advent of satellites in the late 70s, surface temperature is being measured two to four times a day from various satellite sensors and different platforms. The purpose of this paper is to demonstrate use of satellite surface temperature in (a) validation of model computed surface temperatures and (b) assimilation of satellite surface temperatures into a hydrological model in order to improve the prediction accuracy of soil moistures and heat fluxes. The assimilation is carried out by comparing the satellite and the model produced surface temperatures and setting the "true"temperature midway between the two values. Based on this "true" surface temperature, the physical relationships of water and energy balance are used to reset the other variables. This is a case of nudging the water and energy balance variables so that they are consistent with each other and the true" surface temperature. The potential of this assimilation scheme is demonstrated in the form of various experiments that highlight the various aspects. This study is carried over the Red-Arkansas basin in the southern United States (a 5 deg X 10 deg area) over a time period of a year (August 1987 - July 1988). The land surface hydrological model is run on an hourly time step. The results show that satellite surface temperature assimilation improves the accuracy of the computed surface soil moisture remarkably.
Soil moisture impacts the biosphere–atmosphere exchange of CO2 and CH4 and plays an important role in the terrestrial carbon cycle. A better representation of soil moisture would improve coupled carbon–water dynamics in terrestrial ecosystem models and could potentially improve model estimates of large-scale carbon fluxes and climate feedbacks. Here, we investigate using soil moisture observations from the Soil Moisture Active Passive (SMAP) satellite mission to inform simulated carbon fluxes in the global terrestrial ecosystem model LPJ-wsl. Results suggest that the direct insertion of SMAP reduces the bias in simulated soil moisture at in situ measurement sites by 40%, with a greater improvement at temperate sites. A wavelet analysis between the model and measurements from 26 FLUXNET sites suggests that the assimilated run modestly reduces the bias of simulated carbon fluxes for boreal and subtropical sites at 1–2-month time scales. At regional scales, SMAP soil moisture can improve the estimated responses of CO2 and CH4 fluxes to extreme events such as the 2018 European drought and the 2019 rainfall event in the Sudd (Southern Sudan) wetlands. The simulated improvements to land–surface carbon fluxes using the direct insertion of SMAP are shown across a variety of timescales, which suggests the potential of SMAP soil moisture in improving the model representation of carbon–water coupling.
In the Energy Exascale Earth System Model (E3SM) Multi-scale Modeling Framework (MMF), where parameterizations of convection and turbulence are replaced by a 2-D cloud-resolving model (CRM), there are multiple options to represent land–atmosphere interactions. Here, we propose three different coupling strategies, namely the (1) coupling of a single land surface model to the global grid (MMF), (2) coupling a single land copy directly to the embedded CRM (SFLX2CRM), and (3) coupling a single copy of land model to each column of the CRM grid (MAML). In the MAML (Multi-Atmosphere Multi-Land) framework, a land model is coupled to CRM at the CRM-grid scale by coupling an individual copy of a land model to each CRM grid. Therefore, we can represent intra-CRM heterogeneity in the land–atmosphere interaction processes. There are 5-year global simulations run using these three coupling strategies, and we find some regional differences but overall small changes with respect to whether a land model is coupled to CRM or a global atmosphere. In MAML, the spatial heterogeneity within CRM induces stronger turbulence, which leads to the changes in soil moisture, surface heat fluxes, and precipitation. However, the differences in the MAML from the other two cases are rather weak, suggesting that the impact of using MAML does not justify the increase in cost.
The atmosphere can be defined as the relatively thin gaseous envelope surrounding the entire planet Earth. It possesses a number of properties related to its physical state and chemical composition, and it undergoes a variety of internal processes and external interactions that can either maintain or alter these properties. Whereas descriptions of the atmosphere's chemical properties form much of the remaining chapters of this book, the present chapter will highlight the atmosphere's gases, and these define its temperature structure. In contrast, the larger-scale motions comprise the winds, the global organization of which is often referred to as the general circulation. The framework of the dynamical and thermodynamical laws, including the three principles of conversation of mass, momentum, and energy, are fundamental in describing both the internal processes of the atmosphere and its external interactions. The atmosphere is not a closed system, because it exchanges all three of these internally conservative quantities across the atmosphere's boundary below and receives input from regions outside it. Thus surface fluxes of moisture, momentum, and heat occur to and from the underlying ocean and land. The atmosphere exchanges very little mass and momentum with space, though it absorbs directly a portion of the solar radiational energy received from above.
The National Climate Assessment-Land Data Assimilation System (NCA-LDAS) is an Integrated Terrestrial Water Analysis, and is one of NASAs contributions to the NCA of the United States. The NCA-LDAS has undergone extensive development, including multi-variate assimilation of remotely-sensed water states and anomalies as well as evaluation and verification studies, led by the Goddard Space Flight Centers Hydrological Sciences Laboratory (HSL). The resulting NCA-LDAS data have recently been released to the general public and include those from the Noah land-surface model (LSM) version 3.3 (Noah-3.3) and the Catchment LSM version Fortuna-2.5 (CLSM-F2.5). Standard LSM output variables including soil moistures temperatures, surface fluxes, snow cover depth, groundwater, and runoff are provided, as well as streamflow using a river routing system. The NCA-LDAS data are archived at and distributed by the NASA Goddard Earth Sciences Data and Information Services Center (GES DISC). The data can be accessed via HTTP, OPeNDAP, Mirador search and download, and NASA Earth data Search. To further facilitate access and use, the NCA-LDAS data are integrated into the NASA Giovanni, for quick visualization and analysis, and into the Data Rods system, for retrieval of time series of long time periods. The temporal and spatial resolutions of the NCA-LDAS data are, respectively, daily-averages and 0.125x0.125 degree, covering North America (25N 53N; 125W 67W) and the period January 1979 to December 2015. The data files are in self-describing, machine-independent, CF-compliant netCDF-4 format.
Regional surface fluxes are essential lower boundary conditions for large scale numerical weather and climate models and are the elements of global budgets of important trace gases. Surface properties affecting the exchange of heat, moisture, momentum and trace gases vary with length scales from one meter to hundreds of km. A classical difficulty is that fluxes have been measured directly only at points or along lines. The process of scaling up observations limited in space and/or time to represent larger areas was done by assigning properties to surface classes and combining estimated or calculated fluxes using an area weighted average. It is not clear that a simple area weighted average is sufficient to produce the large scale from the small scale, chiefly due to the effect of internal boundary layers, nor is it known how important the uncertainty is to large scale model outcomes. Simultaneous aircraft and tower data obtained in the relatively simple terrain of the western Alaska tundra were used to determine the extent to which surface type variation can be related to fluxes of heat, moisture, and other properties. Surface type was classified as lake or land with aircraft borne infrared thermometer, and flight level heat and moisture fluxes were related to surface type. The magnitude and variety of sampling errors inherent in eddy correlation flux estimation place limits on how well any flux can be known even in simple geometries.
Regional surface fluxes are essential lower boundary conditions for large scale numerical weather and climate models and are the elements of global budgets of important trace gases. Surface properties affecting the exchange of heat, moisture, momentum and trace gases vary with length scales from one meter to hundreds of km. A classical difficulty is that fluxes have been measured directly only at points or along lines. The process of scaling up observations limited in space and/or time to represent larger areas was done by assigning properties to surface classes and combining estimated or calculated fluxes using an area weighted average. It is not clear that a simple area weighted average is sufficient to produce the large scale from the small scale, chiefly due to the effect of internal boundary layers, nor is it known how important the uncertainty is to large scale model outcomes. Simultaneous aircraft and tower data obtained in the relatively simple terrain of the western Alaska tundra were used to determine the extent to which surface type variation can be related to fluxes of heat, moisture, and other properties. Surface type was classified as lake or land with aircraft borne infrared thermometer, and flight level heat and moisture fluxes were related to surface type. The magnitude and variety of sampling errors inherent in eddy correlation flux estimation place limits on how well any flux can be known even in simple geometries.
Abstract Flows in the atmospheric boundary layer are turbulent, characterized by a large Reynolds number, the existence of a roughness sublayer and the absence of a well-defined viscous layer. Exchanges with the surface are therefore dominated by turbulent fluxes. In numerical models for atmospheric flows, turbulent fluxes must be specified at the surface; however, surface fluxes are not known a priori and therefore must be parametrized. Atmospheric flow models, including global circulation, limited area models, and large-eddy simulation, employ Monin–Obukhov similarity theory (MOST) to parametrize surface fluxes. The MOST approach is a semi-empirical formulation that accounts for atmospheric stability effects through universal stability functions. The stability functions are determined based on limited observations using simple regression as a function of the non-dimensional stability parameter representing a ratio of distance from the surface and the Obukhov length scale (Obukhov in Trudy Inst Theor Geofiz AN SSSR 1:95–115, 1946), $$z/L$$ z / L . However, simple regression cannot capture the relationship between governing parameters and surface-layer structure under the wide range of conditions to which MOST is commonly applied. We therefore develop, train, and test two machine-learning models, an artificial neural network (ANN) and random forest (RF), to estimate surface fluxes of momentum, sensible heat, and moisture based on surface and near-surface observations. To train and test these machine-learning algorithms, we use several years of observations from the Cabauw mast in the Netherlands and from the National Oceanic and Atmospheric Administration’s Field Research Division tower in Idaho. The RF and ANN models outperform MOST. Even when we train the RF and ANN on one set of data and apply them to the second set, they provide more accurate estimates of all of the fluxes compared to MOST. Estimates of sensible heat and moisture fluxes are significantly improved, and model interpretability techniques highlight the logical physical relationships we expect in surface-layer processes.
Flooding and drought are two key forecasting challenges for the Kenya Meteorological Department (KMD). Atmospheric processes leading to excessive precipitation and/or prolonged drought can be quite sensitive to the state of the land surface, which interacts with the boundary layer of the atmosphere providing a source of heat and moisture. The development and evolution of precipitation systems are affected by heat and moisture fluxes from the land surface within weakly-sheared environments, such as in the tropics and sub-tropics. These heat and moisture fluxes during the day can be strongly influenced by land cover, vegetation, and soil moisture content. Therefore, it is important to represent the land surface state as accurately as possible in numerical weather prediction models. Enhanced regional modeling capabilities have the potential to improve forecast guidance in support of daily operations and high-end events over east Africa. KMD currently runs a configuration of the Weather Research and Forecasting (WRF) model in real time to support its daily forecasting operations, invoking the Nonhydrostatic Mesoscale Model (NMM) dynamical core. They make use of the National Oceanic and Atmospheric Administration / National Weather Service Science and Training Resource Center's Environmental Modeling System (EMS) to manage and produce the WRF-NMM model runs on a 7-km regional grid over eastern Africa. Two organizations at the National Aeronautics and Space Administration Marshall Space Flight Center in Huntsville, AL, SERVIR and the Short-term Prediction Research and Transition (SPoRT) Center, have established a working partnership with KMD for enhancing its regional modeling capabilities. To accomplish this goal, SPoRT and SERVIR will provide experimental land surface initialization datasets and model verification capabilities to KMD. To produce a land-surface initialization more consistent with the resolution of the KMD-WRF runs, the NASA Land Information System (LIS) will be run at a comparable resolution to provide real-time, daily soil initialization data in place of interpolated Global Forecast System soil moisture and temperature data. Additionally, real-time green vegetation fraction data from the Visible Infrared Imaging Radiometer Suite will be incorporated into the KMD-WRF runs, once it becomes publicly available from the National Environmental Satellite Data and Information Service. Finally, model verification capabilities will be transitioned to KMD using the Model Evaluation Tools (MET) package, in order to quantify possible improvements in simulated temperature, moisture and precipitation resulting from the experimental land surface initialization. The transition of these MET tools will enable KMD to monitor model forecast accuracy in near real time. This presentation will highlight preliminary verification results of WRF runs over east Africa using the LIS land surface initialization.