An Approach for Averaging Surface Temperature and Surface Fluxes Over Heterogeneous Terrain
Successful prediction of possible climate change depends on realistic parameterizations of land-surface processes in climate models.
SEARCH · Search NASA
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.
Successful prediction of possible climate change depends on realistic parameterizations of land-surface processes in climate models.
Dust aerosols consist on a variety of minerals with different physic and chemical properties. As such, they interact with short and long wave radiation, potentially form clouds, act as nutrients modulating biogeochemical cycles, or influence atmospheric chemistry, differently. Most current state-of-the-art Earth System Models (ESMs) neglect the complexity in dust composition, mainly due to computational constraints, but also to the existing uncertainties in the size resolved composition of parent soils, the resulting distribution of minerals in airborne dust, and the scarcity of observations to constrain them. Within this work, we assess the variability of global dust composition due to uncertainties in the characterization of the parent soil mineralogy. To that end, we consider two available global soil mineralogy atlases, developed by Claquin et al. (1999) –C1999- and Journet et al. (2014) –J2014-, which represent respectively 8 and 12 relevant minerals for climate (namely: illite, smectite/montmorillonite, kaolinite, calcite, gypsum, hematite, quartz, and feldspars in C1999, and those plus chlorite, vermiculite, goethite, and mica in J2014). Thanks to a recently developed feature of the MONARCH atmospheric-chemistry model, we are able to explicitly resolve the minerals’ atmospheric cycle. Therefore, we define two global experiments to assess changes on airborne dust composition attributed to the soil mineralogy assumptions and provide a measure of their variability. We also perform a preliminary evaluation of the global mineralogy results against available observations of mineral fractions in surface dust concentration. Our results will inform the climate modelling community about the potential variability in dust composition, an aspect that will gain relevance as ESMs continue growing in complexity and new processes to better characterize aerosols’ forcing or biogeochemical cycles are added. Further observational constraints, such as those that will derive from the EMIT NASA mission on soil composition or the FRAGMENT experimental campaigns on airborne dust characterization, will be key in the near future to improve our understanding of the impact of dust mineralogy on fundamental climate features.
This atlas documents the climate characteristics of version 1 of the NASA Seasonal-to-Interannual Prediction Project (NSIPP) Atmospheric General Circulation Model (AGCM). The AGCM includes an interactive land model (the Mosaic scheme), and is part of the NSIPP coupled atmosphere-land-ocean model. The results presented here are based on a 20-year (December 1979-November 1999) "ANIIP-style" integration of the AGCM in which the monthly-mean sea-surface temperature and sea ice are specified from observations. The climate characteristics of the AGCM are compared with the National Centers for Environmental Prediction (NCEP) and the European Center for Medium-Range Weather Forecasting (ECMWF) reanalyses. Other verification data include Special Sensor Microwave/Imager (SSNM) total precipitable water, the Xie-Arkin estimates of precipitation, and Earth Radiation Budget Experiment (ERBE) measurements of short and long wave radiation. The atlas is organized by season. The basic quantities include seasonal mean global maps and zonal and vertical averages of circulation, variance/covariance statistics, and selected physics quantities.
This study focused on the seasonal relationships and interactions of climate parameters such as the surface temperatures, net radiation, long wave flux, short wave flux, and clouds on a global basis. Five years of observations (December 1984 to November 1989) from the Earth Radiation Budget Experiment (ERBE) and the International Satellite Cloud Climatology Program (ISCCP) were used to study both seasonal variations and interannual variations by use of a basic radiation budget equation. In addition, the study was extended to include an analysis of the cloud forcing due El-Nino's impact on the ERBE parameters.
We examine differences among reanalysis highcloud products in the tropics, assess the impacts of these differences on radiation budgets at the top of the atmosphere and within the tropical upper troposphere and lower stratosphere (UTLS), and discuss their possible origins in the context of the reanalysis models. We focus on the ERA5 (fifthgeneration European Centre for Medium-range Weather Forecasts – ECMWF – reanalysis), ERA-Interim (ECMWF Interim Reanalysis), JRA-55 (Japanese 55-year Reanalysis), MERRA-2 (Modern-Era Retrospective Analysis for Research and Applications, Version 2), and CFSR/CFSv2 (Climate Forecast System Reanalysis/Climate Forecast System Version 2) reanalyses. As a general rule, JRA-55 produces the smallest tropical high-cloud fractions and cloud water contents among the reanalyses, while MERRA-2 produces the largest. Accordingly, long-wave cloud radiative effects are relatively weak in JRA-55 and relatively strong in MERRA-2. Only MERRA-2 and ERA5 among the reanalyses produce tropical-mean values of outgoing long-wave radiation (OLR) close to those observed, but ERA5 tends to underestimate cloud effects, while MERRA-2 tends to overestimate variability. ERA5 also produces distributions of longwave, short-wave, and total cloud radiative effects at the top of the atmosphere that are very consistent with those observed. The other reanalyses all exhibit substantial biases in at least one of these metrics, although compensation between the long-wave and short-wave effects helps to constrain biases in the total cloud radiative effect for most reanalyses. The vertical distribution of cloud water content emerges as a key difference between ERA-Interim and other reanalyses. Whereas ERA-Interim shows a monotonic decrease of cloud water content with increasing height, the other reanalyses all produce distinct anvil layers. The latter is in better agreement with observations and yields very different profiles of radiative heating in the UTLS. For example, whereas the altitude of the level of zero net radiative heating tends to be lower in convective regions than in the rest of the tropics in ERAInterim, the opposite is true for the other four reanalyses. Differences in cloud water content also help to explain systematic differences in radiative heating in the tropical lower stratosphere among the reanalyses. We discuss several ways in which aspects of the cloud and convection schemes impact the tropical environment. Discrepancies in the vertical profiles of temperature and specific humidity in convective regions are particularly noteworthy, as these variables are directly constrained by data assimilation, are widely used, and feed back to convective behaviour through their relationships with thermodynamic stability.
Complex numerical weather prediction models incorporate a variety of physical processes, each described by multiple alternative physical schemes with specific parameters. The selection of the physical schemes and the choice of the corresponding physical parameters during model configuration can significantly impact the accuracy of model forecasts. There is no combination of physical schemes that works best for all times, at all locations, and under all conditions. It is therefore of considerable interest to understand the interplay between the choice of physics and the accuracy of the resulting forecasts under different conditions. This paper demonstrates the use of machine learning techniques to study the uncertainty in numerical weather prediction models due to the interaction of multiple physical processes. The first problem addressed herein is the estimation of systematic model errors in output quantities of interest at future times, and the use of this information to improve the model forecasts. The second problem considered is the identification of those specific physical processes that contribute most to the forecast uncertainty in the quantity of interest under specified meteorological conditions. In order to address these questions we employ two machine learning approaches, random forests and artificial neural networks. The discrepancies between model results and observations at past times are used to learn the relationships between the choice of physical processes and the resulting forecast errors. Numerical experiments are carried out with the Weather Research and Forecasting (WRF) model. The output quantity of interest is the model precipitation, a variable that is both extremely important and very challenging to forecast. The physical processes under consideration include various micro-physics schemes, cumulus parameterizations, short wave, and long wave radiation schemes. The experiments demonstrate the strong potential of machine learning approaches to aid the study of model errors.
The objectives were to apply energy balance principles to plant canopies, and to determine which parameters are essential for predicting plant canopy transpiration (E) in controlled environments. Transpiration was accurately measured in a gas-exchange system. Absorbed radiation (R(sub abs)) by the canopy was measured with a net radiometer and calculated from short and long-wave radiation components. Average canopy foliar temperature T(sub L) can be measured with an infrared radiometer, but since T(sub L) is seldom uniform, a weighed average measurement of T(sub L) must be made. The effective canopy temperature T(sub C) is that temperature that balances the energy flux between absorbed radiation and latent heat L(sub E) and sensible heat (H) fluxes. TC should exactly equal air temperature T(sub A) when L(sub E) equals R(sub abs). When unnecessary thermal radiation from the lighting system is removed by a water filter, the magnitude of L(sub E) from transpiration approaches Rabs and T(sub C) is close to T(sub A). Unlike field models, we included the energy used in photosynthesis and found that up to 10% of Rabs was used in photosynthesis. We calculated aerodynamic conductance for H from measurements of wind speed and canopy height using the wind profile equation. Canopy aerodynamic conductance ranged from.03 to.04 m/s for wind speeds from.6 to 1 m/s; thus a 0.1 C canopy to air temperature difference results in a sensible heat flux of about 4 W/sq m, which is only 1% of R(sub abs). We examined the ability of wide angle infrared transducers to accurately integrate T(sub L) from the top to the bottom of the canopy. We measured evaporation from the hydroponic media to be approximately 1 micro mol/sq m s or 10% of R(sub abs). This result indicates that separating evaporation from transpiration is more important than exact measurement of canopy temperature.
Schoeberl et al. (1983) suggested that convective wavebreaking of monochromatic gravity waves might be suppressed by radiative transfer processes if the vertical wavelength waves were sufficiently short. As the vertical wavelength or the gravity wave decreases, radiative transfer between adjacent vertical layers becomes increasingly important. This exchange can increase the radiative relaxation time scale so that the wave will no longer grow with altitude. Thus, very short vertical wavelength waves may dissipate radiatively rather than become convectively unstable. Apruzese and Strobel (1984) have revised the exchange coefficients used in Schoeberl et al. (1983). Also, Chao and Schoeberl (1984) pointed out that the computation made by LINDZEN (1981) of the convective diffusion rate may be a factor of two too low as the convective adjustment processes tends to minimize the thermal transport by the wave. The purpose here is to revise the values given in Schoeberl et al. (1983). These results also suggest that the very thin turbulent layers observed by mesosphere-stratosphere-troposphere (MST) radars (e.g., WOODMAN, 1980) cannot be produced by the convective instability of monochromatic gravity waves with large horizontal scales.
Surface solar radiation reaching the Earth's surface is the primary forcing function of the land surface energy and water cycle. Therefore, there is a need for information on this parameter, preferably, at global scale. Satellite based estimates are now available at accuracies that meet the demands of many scientific objectives. Selection of an approach to estimate such fluxes requires consideration of trade-offs between the use of multi-spectral observations of cloud optical properties that are more difficult to implement at large scales, and methods that are simplified but easier to implement. In this study, an evaluation of such trade-offs will be performed. The University of Maryland Surface Radiation Model (UMD/SRB) has been used to reprocess five years of GOES-8 satellite observations over the United States to ensure updated calibration and improved cloud detection over snow. The UMD/SRB model was subsequently modified to allow input of information on aerosol and cloud optical depth with information from independent satellite sources. Specifically, the cloud properties from the Atmospheric Radiation Measurement (ARM) Satellite Data Analysis Program (Minnis et al., 1995) are used to drive the modified version of the model to estimate surface short-wave fluxes over the Southern Great Plain ARM sites for a twelve month period. The auxiliary data needed as model inputs such as aerosol optical depth, spectral surface albedo, water vapor and total column ozone amount were kept the same for both versions of the model. The estimated shortwave fluxes are evaluated against ground observations at the ARM Central Facility and four satellite ARM sites. During summer, the estimated fluxes based on cloud properties derived from the multi-spectral approach were in better agreement with ground measurements than those derived from the UMD/SRB model. However, in winter, the fluxes derived with the UMD/SRB model were in better agreement with ground observations than those estimated from cloud properties provided by the ARM Satellite Data Analysis Program. During the transition periods, the results were comparable.
We examine the characteristics of the classical radiation emitted by a relativistic electron beam that propagates perpendicularly through a large amplitude relativistic plasma wave. Such a study is useful for evaluating the feasibility of using relativistic plasma waves as extremely short wavelength undulators for generating short wavelength radiation. The electron trajectories in a plasma wave undulator are obtained using perturbation techniques and are then compared to numerical simulation results. The frequency spectrum and angular distribution of the spontaneous radiation emitted by a single electron and the stimulated radiation gain are obtained analytically, and are then compared to 3-D numerical simulations. The characteristics of the plasma wave undulator are compared to the ac FEL undulator and the conventional FEL.
Soil dust aerosols are a key component of the climate system, as they interact with short- and long-wave radiation, alter cloud formation processes, affect atmospheric chemistry and play a role in biogeochemical cycles by providing nutrient inputs such as iron and phosphorus. The influence of dust on these processes depends on its physicochemical properties, which, far from being homogeneous, are shaped by its regionally varying mineral composition. The relative amount of minerals in dust depends on the source region and shows a large geographical variability. However, many state-of-the-art Earth system models (ESMs), upon which climate analyses and projections rely, still consider dust mineralogy to be invariant. The explicit representation of minerals in ESMs is more hindered by our limited knowledge of the global soil composition along with the resulting size-resolved airborne mineralogy than by computational constraints. In this work we introduce an explicit mineralogy representation within the state-of-the-art Multiscale Online Nonhydrostatic AtmospheRe CHemistry (MONARCH) model. We review and compare two existing soil mineralogy datasets, which remain a source of uncertainty for dust mineralogy modeling and provide an evaluation of multiannual simulations against available mineralogy observations. Soil mineralogy datasets are based on measurements performed after wet sieving, which breaks the aggregates found in the parent soil. Our model predicts the emitted particle size distribution (PSD) in terms of its constituent minerals based on brittle fragmentation theory (BFT), which reconstructs the emitted mineral aggregates destroyed by wet sieving. Our simulations broadly reproduce the most abundant mineral fractions independently of the soil composition data used. Feldspars and calcite are highly sensitive to the soil mineralogy map, mainly due to the different assumptions made in each soil dataset to extrapolate a handful of soil measurements to arid and semi-arid regions worldwide. For the least abundant or more difficult-to-determine minerals, such as iron oxides, uncertainties in soil mineralogy yield differences in annual mean aerosol mass fractions of up to ~ 100 %. Although BFT restores coarse aggregates including phyllosilicates that usually break during soil analysis, we still identify an overestimation of coarse quartz mass fractions (above 2 µm in diameter). In a dedicated experiment, we estimate the fraction of dust with undetermined composition as given by a soil map, which makes up ~ 10 % of the emitted dust mass at the global scale and can be regionally larger. Changes in the underlying soil mineralogy impact our estimates of climate-relevant variables, particularly affecting the regional variability of the single-scattering albedo at solar wavelengths or the total iron deposited over oceans. All in all, this assessment represents a baseline for future model experiments including new mineralogical maps constrained by high-quality spaceborne hyperspectral measurements, such as those arising from the NASA Earth Surface Mineral Dust Source Investigation (EMIT) mission.
Clouds interfere with the distribution of short-wave and long-wave radiations over sea ice, and thereby strongly affect the surface energy balance in polar regions. To evaluate the overall effects of clouds on climatic feedback processes in the atmosphere-ice-ocean system, the challenge is to observe sea ice surface thermal states under both clear sky and cloudy conditions. From laboratory experiments, we show that C-band radar (transparent to clouds) backscatter is very sensitive to the surface temperature of first-year sea ice. The effect of sea ice surface temperature on the magnitude of backscatter change depends on the thermal regimes of sea ice thermodynamic states. For the temperature range above the mirabilite (Na2SO4.10H20) crystallization point (-8.2 C), C-band data show sea ice backscatter changes by 8-10 dB for incident angles from 20 to 35 deg at both horizontal and vertical polarizations. For temperatures below the mirabilite point but above the crystallization point of MgCl2.8H2O (-18.0 C), relatively strong backwater changes between 4-6 dB are observed. These backscatter changes correspond to approximately 8 C change in temperature for both cases. The backscattering mechanism is related to the temperature which determines the thermodynamic distribution of brine volume in the sea ice surface layer. The backscatter is positively correlated to temperature and the process is reversible with thermodynamic variations such as diurnal insolation effects. From two different dates in May 1993 with clear and overcast conditions determined by the Advanced Very High Resolution Radiometer (AVHRR), concurrent Earth Resources Satellite 1 (ERS-1) C-band ice observed with increases in backscatter over first-year sea ice, and verified by increases in in-situ sea ice surface temperatures measured at the Collaborative-Interdisciplinary Cryosphere Experiment (C-ICE) site.
Clouds interfere with the distribution of short-wave and long-wave radiations over sea ice, and thereby strongly affect the surface energy balance in polar regions.
Soil dust aerosols are a key component of the climate system, as they interact with short- and long-wave radiation, alter cloud formation processes, affect atmospheric chemistry and play a role in biogeochemical cycles by providing nutrient inputs such as iron and phosphorus. The influence of dust on these processes depends on its physicochemical properties, which, far from being homogeneous, are shaped by its regionally varying mineral composition. The relative amount of minerals in dust depends on the source region and shows a large geographical variability. However, many state-of-the-art Earth system models (ESMs), upon which climate analyses and projections rely, still consider dust mineralogy to be invariant. The explicit representation of minerals in ESMs is more hindered by our limited knowledge of the global soil composition along with the resulting size-resolved airborne mineralogy than by computational constraints. In this work we introduce an explicit mineralogy representation within the state-of-the-art Multiscale Online Nonhydrostatic AtmospheRe CHemistry (MONARCH) model. We review and compare two existing soil mineralogy datasets, which remain a source of uncertainty for dust mineralogy modeling and provide an evaluation of multiannual simulations against available mineralogy observations. Soil mineralogy datasets are based on measurements performed after wet sieving, which breaks the aggregates found in the parent soil. Our model predicts the emitted particle size distribution (PSD) in terms of its constituent minerals based on brittle fragmentation theory (BFT), which reconstructs the emitted mineral aggregates destroyed by wet sieving. Our simulations broadly reproduce the most abundant mineral fractions independently of the soil composition data used. Feldspars and calcite are highly sensitive to the soil mineralogy map, mainly due to the different assumptions made in each soil dataset to extrapolate a handful of soil measurements to arid and semi-arid regions worldwide. For the least abundant or more difficult-to-determine minerals, such as iron oxides, uncertainties in soil mineralogy yield differences in annual mean aerosol mass fractions of up to ∼ 100 %. Although BFT restores coarse aggregates including phyllosilicates that usually break during soil analysis, we still identify an overestimation of coarse quartz mass fractions (above 2 µm in diameter). In a dedicated experiment, we estimate the fraction of dust with undetermined composition as given by a soil map, which makes up ∼ 10 % of the emitted dust mass at the global scale and can be regionally larger. Changes in the underlying soil mineralogy impact our estimates of climate-relevant variables, particularly affecting the regional variability of the single-scattering albedo at solar wavelengths or the total iron deposited over oceans. All in all, this assessment represents a baseline for future model experiments including new mineralogical maps constrained by high-quality spaceborne hyperspectral measurements, such as those arising from the NASA Earth Surface Mineral Dust Source Investigation (EMIT) mission.
Numerical experiments with the NASA finite-volume general circulation model show that heating of the atmosphere by dust and black carbon can lead to widespread enhanced warming over the Tibetan Plateau (TP) and accelerated snow melt in the western TP and Himalayas. During the boreal spring, a thick aerosol layer, composed mainly of dust transported from adjacent deserts and black carbon from local emissions, builds up over the Indo-Gangetic Plain, against the foothills of the Himalaya and the TP. The aerosol layer, which extends from the surface to high elevation (approx.5 km), heats the mid-troposphere by absorbing solar radiation. The heating produces an atmospheric dynamical feedback the so-called elevated-heat-pump (EHP) effect, which increases moisture, cloudiness, and deep convection over northern India, as well as enhancing the rate of snow melt in the Himalayas and TP. The accelerated melting of snow is mostly confined to the western TP, first slowly in early April and then rapidly from early to mid-May. The snow cover remains reduced from mid-May through early June. The accelerated snow melt is accompanied by similar phases of enhanced warming of the atmosphere-land system of the TP, with the atmospheric warming leading the surface warming by several days. Surface energy balance analysis shows that the short-wave and long-wave surface radiative fluxes strongly offset each other, and are largely regulated by the changes in cloudiness and moisture over the TP. The slow melting phase in April is initiated by an effective transfer of sensible heat from a warmer atmosphere to land. The rapid melting phase in May is due to an evaporation-snow-land feedback coupled to an increase in atmospheric moisture over the TP induced by the EHP effect.
Seasonal and regional variations of clouds and their effects on the climatological parameters were studied. The climatological parameters surface temperature, solar insulation, short-wave absorbed, long wave emitted, and net radiation were considered. The data of climatological parameters consisted of about 20 parameters of Earth radiation budget and clouds of 2070 target areas which covered the globe. It consisted of daily and monthly averages of each parameter for each target area for the period, Jun. 1979 - May 1980. Cloud forcing and black body temperature at the top of the atmosphere were calculated. Interactions of clouds, cloud forcing, black body temperature, and the climatological parameters were investigated and analyzed.
Short pulses of millimeter wave (MMW) radiation at 43 GHz create microplasma within a photonic crystal for pressures from 40 to 600 Torr (1.3 × 10 3 –8.0 × 10 4 Pa). Gas breakdown occurs within a photonic crystal, which acts as an electromagnetic resonator to create a strong initial electric field. The time response of the argon metastable density is experimentally determined during the pulse and in the afterglow using laser absorption. The metastable density overshoots the steady-state condition at the beginning of the pulse and during the afterglow. Modeling is presented to understand these observations. The overproduction of argon 1s5 at the beginning of each pulse is due to a concurrent overshoot in the MMW electric field within the photonic crystal. This field overheats the plasma electrons and enhances the production of excited states. The burst of argon metastables observed in the afterglow is due to the pooled energy of the plasma stored in electrons, ions, and excited states of argon. As a result, understanding metastable production is an important intermediate step to ionization and is also critical in the study of diode-pumped rare gas lasers.
The mission of the Data Assimilation Office (DAO) is to advance the state of the art of data assimilation and produce research-quality assimilated data sets which make optimal use of space-based observations. Development efforts over the last few years have focused on delivering a production data assimilation system in support of NASA's Terra platform. That system, called the Goddard Earth Observing System - version 2 or GEOS-2, represents a major upgrade to the baseline GEOS-1 system employed in NASA's first reanalysis effort. GEOS-2 includes a physical-space three dimensional variational analysis algorithm (the Physical-space Statistical Analysis System or PSAS) and numerous improvements to the general circulation model. The latter include a Soil-Vegetation-Atmosphere Transfer (SVAT) land surface scheme, a level 2.5 moist turbulence scheme and new Short Wave (SW) and Long Wave (LW) radiation code. The system also includes an off-line ozone assimilation system, and the capability to assimilate scatterometer surface winds, and TIROS Operational Vertical Sounder (TOVS) and Special Sensor Microwave Imager (SSM/I) moisture data. GEOS-2 is currently run at 1 degree horizontal resolution and 48 levels extending to O.Olmb. Experimental versions of GEOS-2 are run with a global stretched grid allowing enhanced (e.g. 1/4 deg) regional resolution. Other capabilities being developed include, the assimilation of Tropical Rainfall Measuring Mission (TRMM) precipitation and Global Positioning System (GPS) data, an off-line land surface assimilation system, and a retrospective analysis scheme. The DAO is also engaged in a number of collaborative efforts to help accelerate the development of the next generation data assimilation system. These include, a joint modeling effort between the DAO and NCAR/CGDD to develop a new Global Circulation Model (GCM), and a Department of Energy Lawrence Livermore National Laboratory (DOE/LLNL) collaboration on model parallelization. Plans for the next reanalysis will be discussed in the context of current and near term system quality and computing capabilities, and the need for multiple reanalysis products.