Study of earth's atmosphere final report
Earth and planetary atmospheres, isopycnic levels, temperature determination from diffusion data, and transition model atmospheres
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.
Earth and planetary atmospheres, isopycnic levels, temperature determination from diffusion data, and transition model atmospheres
The ability of the Weather Research and Forecasting model with chemistry (WRF-Chem) version 3.7 and the Community Atmosphere Model version 5.3 (CAM5) in simulating profiles of aerosol properties is quantified using extensive in situ and remote sensing measurements from the Two Column Aerosol Project (TCAP) conducted during July of 2012. TCAP was supported by the U.S. Department of Energy’s Atmospheric Radiation Measurement program and was designed to obtain observations within two atmospheric columns; one fixed over Cape Cod, Massachusetts and the other several hundred kilometers over the ocean. The performance is quantified using most of the available aircraft and surface measurements during July, and two days are examined in more detail to identify the processes responsible for the observed aerosol layers. The higher resolution WRF-Chem model produced more aerosol mass in the free troposphere than the coarser resolution CAM5 model so that the fraction of aerosol optical thickness above the residual layer from WRF-Chem was more consistent with lidar measurements. We found that the free troposphere layers are likely due to mean vertical motions associated with synoptic-scale convergence that lifts aerosols from the boundary layer. The vertical displacement and the time period associated with upward transport in the troposphere depend on the strength of the synoptic system and whether relatively high boundary layer aerosol concentrations are present where convergence occurs. While a parameterization of subgrid scale convective clouds applied in WRF-Chem modulated the concentrations of aerosols aloft, it did not significantly change the overall altitude and depth of the layers.
Detecting changes in the Earth's environment using satellite images of ocean and land surfaces must take into account atmospheric effects. As a result, major programs are underway to develop algorithms for image retrieval of atmospheric aerosol properties and atmospheric correction. However, because of the temporal and spatial variability of atmospheric transmittance it is very difficult to model atmospheric effects and implement models in an operational mode. For this reason, simultaneous in situ ground measurements of atmospheric optical properties are vital to the development of accurate atmospheric correction techniques. Presented in this paper is a spectroradiometer system that provides an optimized set of surface measurements for the calibration and validation of atmospheric correction algorithms. The Portable Ground-based Atmospheric Monitoring System (PGAMS) obtains a comprehensive series of in situ irradiance, radiance, and reflectance measurements for the calibration of atmospheric correction algorithms applied to multispectral. and hyperspectral images. The observations include: total downwelling irradiance, diffuse sky irradiance, direct solar irradiance, path radiance in the direction of the north celestial pole, path radiance in the direction of the overflying satellite, almucantar scans of path radiance, full sky radiance maps, and surface reflectance. Each of these parameters are recorded over a wavelength range from 350 to 1050 nm in 512 channels. The system is fast, with the potential to acquire the complete set of observations in only 8 to 10 minutes depending on the selected spatial resolution of the sky path radiance measurements
We examine the changes in clouds and cloud feedback between Phase 5 (AMIP5) and Phase 6 (AMIP6) of the Atmospheric Model Intercomparison Project. Each model is perturbed by uniformly increasing the sea surface temperature by 4 K. The simulated cloud fraction, the perturbed states and cloud radiative kernels are used to derive cloud feedback in the shortwave (SW), longwave (LW) and their sum (Net). Compared to AMIP5, the cloud fraction in AMIP6 increases by 9.1%, while the perturbation leads to a 0.25% decrease. The Net cloud feedback at the top of the atmosphere (TOA) is almost double (174%). Statistical tests support that this change is mainly due to an increase in the surface SW cloud feedback caused by optically thick, middle and low clouds. The contribution of the atmospheric Net component (12%) stems from the increase in the atmospheric LW cloud feedback, likely to play a role in weakening (strengthening) the northward (southward) meridional atmospheric energy transport, while the opposite is true for the surface LW and Net cloud feedback in the meridional oceanic energy transport. The substantial increase in cloud feedback at the TOA primarily contributes to the higher climate sensitivity. The cloud feedback spread in AMIP6 is comparable to that in AMIP5.
Temperature correction procedure for atmosphere model calculation
We analyze new/archival VLT/NaCo and Gemini/NICI high-contrast imaging of the young, self-luminous planet Beta Pictoris b in seven near-to-mid IR photometric filters, using advanced image processing methods to achieve high signal-to-noise, high precision measurements. While Beta Pic b's near-IR colors mimic those of a standard, cloudy early-to-mid L dwarf, it is overluminous in the mid-infrared compared to the field L/T dwarf sequence. Few substellar/planet-mass objects-i.e., κ And b and 1RXJ 1609B-match Beta Pic b's JHKsL photometry and its 3.1 micron and 5 micron photometry are particularly difficult to reproduce. Atmosphere models adopting cloud prescriptions and large (approx. 60 micron)dust grains fail to reproduce the Beta Pic b spectrum. However, models incorporating thick clouds similar to those found forHR8799 bcde, but also with small (a fewmicrons) modal particle sizes, yield fits consistent with the data within the uncertainties. Assuming solar abundance models, thick clouds, and small dust particles (a = 4 micron), we derive atmosphere parameters of log(g) = 3.8 +/- 0.2 and Teff = 1575-1650 K, an inferred mass of 7+4 −3 MJ, and a luminosity of log(L/L) approx. −3.80 +/- 0.02. The best-estimated planet radius, is approx. equal to 1.65 +/- 0.06 RJ, is near the upper end of allowable planet radii for hot-start models given the host star's age and likely reflects challenges constructing accurate atmospheric models. Alternatively, these radii are comfortably consistent with hot-start model predictions if Beta Pic b is younger than is approx. equal to 7 Myr, consistent with a late formation well after its host star's birth approx. 12+8 −4 Myr ago.
We analyze new/archival VLT/NaCo and Gemini/NICI high-contrast imaging of the young, self-luminous planet β Pictoris b in seven near-to-mid IR photometric filters, using advanced image processing methods to achieve high signal-to-noise, high precision measurements. While β Pic b's near-IR colors mimic those of a standard, cloudy early-to-mid L dwarf, it is overluminous in the mid-infrared compared to the field L/T dwarf sequence. Few substellar/planet-mass objects—i.e., κ And b and 1RXJ 1609B—match β Pic b's JHK{sub s}L' photometry and its 3.1 μm and 5 μm photometry are particularly difficult to reproduce. Atmosphere models adopting cloud prescriptions and large (∼60 μm) dust grains fail to reproduce the β Pic b spectrum. However, models incorporating thick clouds similar to those found for HR 8799 bcde, but also with small (a few microns) modal particle sizes, yield fits consistent with the data within the uncertainties. Assuming solar abundance models, thick clouds, and small dust particles ((a) = 4 μm), we derive atmosphere parameters of log (g) = 3.8 ± 0.2 and T{sub eff} = 1575-1650 K, an inferred mass of 7{sup +4}{sub -3} M{sub J} , and a luminosity of log(L/L{sub ☉}) ∼–3.80 ± 0.02. The best-estimated planet radius, ≈1.65 ± 0.06 R{sub J} , is near the upper end of allowable planet radii for hot-start models given the host star's age and likely reflects challenges constructing accurate atmospheric models. Alternatively, these radii are comfortably consistent with hot-start model predictions if β Pic b is younger than ≈7 Myr, consistent with a late formation well after its host star's birth ∼12{sup +8}{sub -4} Myr ago.
The first young giant planet to be discovered by the Gemini Planet Imager was the (is) approximately 2MJ planet 51 Eri b. This approximately 20 Myr old young Jupiter is the first directly imaged planet to show unmistakable methane in H band. To constrain the planet's mass, atmospheric temperature, and composition, the GPI J and H band spectra as well as some limited photometric points were compared to the predictions of substellar atmosphere models. The best fitting models reported in the discovery paper (Macintosh et al. 2015) relied upon a combination of clear and cloudy atmospheric columns to reproduce the data. However for an object as cool as 700 K, the origin of the cloud coverage is somewhat puzzling, as the global silicate and iron clouds would be expected to have sunk well below the photosphere by this effective temperature. While strong vertical mixing in these low gravity atmospheres remains a plausible explanation, we have explored whether atmospheric photochemistry, driven by the UV flux from the primary star, may yield hazes that also influence the observed spectrum of the planet. To explore this possibility we have modeled the atmospheric photochemistry of 51 Eri b using two state-of-the-art photochemical models, both capable of predicting yields of complex hydrocarbons under various atmospheric conditions. In our presentation we will summarize the modeling approach employed to characterize 51 Eri b, explaining constraints on the planet's effective temperature, gravity, and atmospheric composition and also present results of our studies of atmospheric photochemistry. We will discuss whether photochemical hazes could indeed be responsible for the particulate opacity that apparently sculpts the spectrum of the planet.
Mars atmosphere models at three density levels
We present a deep near-infrared image of the newly discovered brown dwarf WISE J085510.83-071442.5 (W0855) using the FourStar imager at Las Campanas Observatory. Our detection of J3 = 24.8{sub −0.35}{sup +0.53} (J {sub MKO} = 25.0{sub −0.35}{sup +0.53}) at 2.6σ—or equivalently an upper limit of J3 > 23.8 (J {sub MKO} > 24.0) at 5σ makes W0855 the reddest brown dwarf ever categorized (J {sub MKO} – W2 = 10.984{sub −0.35}{sup +0.53} at 2.6σ—or equivalently an upper limit of J {sub MKO} – W2 > 9.984 at 5σ) and refines its position on color-magnitude diagrams. Comparing the new photometry with chemical equilibrium model atmosphere predictions, we demonstrate that W0855 is 2.7σ from models using a cloudless atmosphere and well reproduced by partly cloudy models (50%) containing sulfide and water ice clouds. Non-equilibrium chemistry or non-solar metallicity may change predictions, however using currently available model approaches, this is the first candidate outside our own solar system to have direct evidence for water clouds.
T Coronae Borealis is a recurrent, symbiotic nova system currently in quiescence between its periodic ≈80 yr cycle of eruptions. Observations during inter-outburst epochs provide an opportunity to study properties of the accretion disk and the M red giant. Here we present new irradiated (blackbody veiling) models, incorporating modern molecular opacities and line lists, of spectra derived from high-resolution (22,000 ≲ R ≲ 120,000) optical echelle observations obtained at two epochs, one prior to and one post the 2015 rebrightening event at similar spectroscopic system phase. We find a lithium abundance in the secondary at both epochs to be comparable. The non-irradiated (classical) model atmospheres yield a lithium abundance, A(Li) = 1.3 ± 0.1. The irradiated model (veiled) atmospheres, which are likely a better representation of the system in which the white dwarf and accretion disk illuminate the red giant, give A(Li) = 2.4 ± 0.1.
The Modern-Era Retrospective Analysis for Research and Applications, version 2 (MERRA-2), provides global, 1-hourly estimates of land surface conditions for 1980-present at 50-km resolution. Outside of the high latitudes, MERRA-2 uses observations-based precipitation data products to correct the precipitation falling on the land surface. This paper describes the precipitation correction method and evaluates the MERRA-2 land surface precipitation and hydrology. Compared to monthly GPCPv2.2 observations, the corrected MERRA-2 precipitation (M2CORR) is better than the precipitation generated by the atmospheric models within the cyclingMERRA-2 system and the earlier MERRA reanalysis. Compared to 3-hourlyTRMM observations, the M2CORR diurnal cycle has better amplitude but less realistic phasing than MERRA-2 model-generated precipitation. Because correcting the precipitation within the coupled atmosphere-land modeling system allows the MERRA-2 near-surface air temperature and humidity to respond to the improved precipitation forcing, MERRA-2 provides more self-consistent surface meteorological data than were available from the earlier, offline MERRA-Land reanalysis. Overall, MERRA-2 land hydrology estimates are better than those of MERRA-Land and MERRA. A comparison against GRACE satellite observations of terrestrial water storage demonstrates clear improvements in MERRA-2 over MERRA in South America and Africa but also reflects known errors in the observations used to correct the MERRA-2 precipitation. The MERRA-2 and MERRA-Land surface and root zone soil moisture skill vs. in situ measurements is slightly higher than that of ERA-Interim Land and higher than that of MERRA (significantly for surface soil moisture). Snow amounts from MERRA-2 have lower bias and correlate better against reference data than do those of MERRA-Land and MERRA, with MERRA-2 skill roughly matching that of ERA-Interim Land. Seasonal anomaly R values against naturalized stream flow measurements in the United States are, on balance, highest for MERRA-2 and ERA-Interim Land, somewhat lower for MERRA-Land, and lower still for MERRA.
CO 2.6-mm wavelength observations of Mars during March 29-April 1, 1980, were analyzed using a model atmosphere based on Viking measurements. A least-squares fit of the model to the observed line profile yielded an average CO mixing ratio of 0.0032 + or - 0.0011, a value four times larger than that obtained by Kaplan (1969) from the analysis of an infrared spectrum obtained in 1967. Martian atmospheric chemistry models indicate that the temporal variation implied could exist, being due primarily to variations in H2O abundance.
Abstract Ural blocking (UB) is a prominent mode of variability of the Northern Hemisphere atmospheric circulation, particularly in fall. It can persist for several days and exert a lagged influence on the wintertime NH circulation, providing predictability at the subseasonal time scale. Using two atmospheric models, we explore how the early winter atmospheric circulation responds to a 2‐week persistent UB anomaly imposed in early November. Experiments are carried out with two different configurations of Barents‐Kara (BK) sea‐ice concentration to examine whether it plays a role in how UB impacts atmospheric variability. In both models, the UB anomaly is followed by a weakening of the stratospheric polar vortex and a negative phase of the North Atlantic Oscillation (NAO), which lasts up to 2 months after the forcing is released. Interestingly, the response is more persistent under low BK sea‐ice conditions, that is, BK sea ice modulates the atmospheric response to UB. Additional experiments with prescribed sea ice concentration anomalies alone suggest that BK sea ice exerts a limited influence on early winter NH atmospheric variability. The response to UB involves a weakening of the polar vortex that persists longer under low BK sea ice, which explains the more persistent response in that configuration. Our study highlights that UB variability in November is a robust precursor for early winter NAO/polar stratosphere anomalies, and this may be more relevant in the context of declining Arctic sea‐ice extent. Provided that climate models accurately capture this teleconnection, it has the potential to improve subseasonal predictions of the NH wintertime climate.
An atmospheric radiation model is used here to illustrate several features associated with modeling the diurnal cycle of the planetary albedo. It is found that even for clear regions there appear to be deficiencies in our knowledge of how to model this quantity. The diurnal amplitude factor, defined as the ratio of the diurnally averaged planetary albedo to that at noon, between two GCMs and measurements made from a geostationary satellite. While reasonable consistency is found, the comparisons underscore difficulties associated with converting local-time albedo measurements, as made from sun-synchronous satellites, to diurnally averaged albedos.
During its emergence from the Saturnian shadow on June 28, 1980, Titan was observed simultaneously in the visible and IR ranges (6000-9000 A, and 11.8 and 20 microns). The eclipse observation indicates that thin-atmosphere models with a modest greenhouse effect are likely to be incorrect. A thin atmosphere of pure methane in equilibrium with an ice surface of methane is in marginal agreement with the present data, but improved VLA temperature measurements would be a crucial test of this model. Both the eclipse observation and VLA measurements agree with a nitrogen-dominant atmospheric model, as indicated by preliminary Voyager results.
Given the importance of offshore wind energy development to the U.S. clean energy targets, it is vital to be able to characterize the wind resource in that environment accurately. Toward that end, two Bureau of Ocean Energy Management buoys equipped with Doppler lidar are being maintained by Pacific Northwest National Laboratory on behalf of the Department of Energy and deployed to regions of potential offshore wind development. In addition to standard meteorological and oceanographic measurements, the buoys document the wind profile between about 40 m and 250 m above the sea surface through Doppler lidar retrievals. After a multiyear deployment of two buoys along the U.S. East Coast, the buoys were redeployed to the U.S. West coast from 2020 – 2022 to locations near the Humboldt and Morro Bay lease areas. The buoys provide nearly continuous, multiyear datasets that can be used to evaluate predictions of hub-height (~100 m) wind speed for standard atmospheric models in the region. In the absence of measurements at the study site, offshore wind developers rely on model-based data to assess site conditions. Potential sources of model error in this environment include under-resolution or misrepresentation of coastal topographically forced flows, marine boundary layer dynamics and the evolution of their associated cloud and turbulence fields, the role of upwelling and other currents on surface heat fluxes into the boundary layer, and the impact of wave fields on surface momentum fluxes and thus the wind speed profile. In particular, most predictive models of wind speed do not predict wave fields at all, relying on parameterizations to represent their effects. In thus study, we focus on evaluating the role of wind / wave interactions on modeled hub-height wind speed and error by using the Coupled Ocean–Atmosphere–Wave–Sediment–Transport Modeling System to capture two-way interactions between an atmospheric model (Weather Research and Forecasting (WRF)) and a wave model (WAVEWATCHIII (WW3)) and compare to both stand-alone WRF and one-way coupled WRF / WW3 configurations. Our approach is similar to that used in Gaudet et al. (2022) to evaluate wind / wave coupling over the U.S. East Coast, but applied to the very different environment of the U.S. West Coast. We show examples for two cases, a cold-season frontal case and a warm-season low-level jet case. We find that wind / wave coupling makes little impact on model error for these cases at the location of the lidar buoys, for which other misrepresentations of model physics seems to be responsible for model-observation discrepancies. However, domain-wide evaluations, which also make use of the National Buoy Data Center network, show that a two-way coupling approach is less prone to systematic errors in the hub-height wind field than the one-way coupled approach. WRF resolution of kilometer-scale or less is needed to properly capture the sharp wind speed gradients that can be found along the coastline, and WW3 simulations driven by the downscaled WRF produce better bulk and spectral wave fields when compared to observations. Implications of the results for wind resource characterization are then discussed.
This work is an intercomparison of four tropospheric-stratospheric climate models, the Unified Model (UM) of the U.K. Meteorological Office (UKMO), the model of the Free University in Berlin (FUB). the ARPEGE-climat model of the National Center for Meteorological Research (CNRM), and the Extended UGAMP GCM (EUGCM) of the Center for Global Atmospheric Modelling (CGAM), against the UKMO analyses. This comparison has been made in the framework of the "GSM-Reality Intercomparison Project for SPARC" (GRIPS). SPARC (Stratospheric Processes and their Role in Climate) aims are to investigate the effects of the middle atmosphere on climate and the GRIPS purpose is to organized a comprehensive assessment of current Middle Atmosphere-Climate Models (MACMs). The models integrations were made without identical contraints e.g. boundary conditions, incoming solar radiation). All models are able to represent the dominant features of the extratropical circulation. In this paper, the structure of the tropical winds and the strengths of the Kelvin waves are examined. Explanations for the differences exhibited. between the models. as well as between models and analyses, are also proposed. In the analyses a rich spectrum of waves (eastward and westward) is present and contributes to drive the SAO (SemiAnnual Oscillation) and the QBO (Quasi-Biennal Oscillation). The amplitude of the Kelvin waves is close to the one observed in UARS (Upper Atmosphere Research Satellite) data. In agreement with observations, the Kelvin waves generated in the models propagate into the middle atmosphere as wave packets which underlines convective forcing origin. In most models, slow Kelvin waves propagate too high and are hence overestimated in the upper stratosphere and in the mesosphere, except for the UM which is more diffusive. These waves are not sufficient to force realistic westerlies of the QBO or SAO westerly phases. If the SAO is represented by all models only two of them are able to generate westerlies between 10 hPa and 50 hPa. The importance of the role played by subgrided gravity waves is more and more recognized. Actually, the EUGCM which includes a parametrization of gravity waves with a non-zero phase speed is able to simulate. with however some unrealistic features, clear easterly to westerly transitions as well as westerlies downward propagations. Thermal damping is also important in the westerlies forcing in the stratosphere. The model ARPEGE-climat shows more westerlies in the stratosphere than tile other three models probably due to the use of a simplified scheme to predict the ozone distribution in the middle atmosphere.