Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “typical meteorological year”

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 37 records · Page 2

Timing the Flames: Geostationary Satellite Detection of Diurnally Shifting Stubble Burning in Northwestern India

Post-monsoon open-field stubble burning in northwestern (NW) India—a key agricultural region known as the “breadbasket”—is a longstanding practice used to clear fields. Satellite observations spanning over two decades have revealed significant upward trends in crop production, vegetative greenness, and the frequency of post-harvest fires, with this last contributing to hazardous air quality during the peak burning season (mid-October to mid-November). Since 2022, thermal anomaly data from Aqua-MODIS and SNPP-VIIRS sensors have shown a sharp decline in reported fire events—an observation that contrasts starkly with the concurrent rise in regional aerosol loading detected from space. This apparent discrepancy became particularly pronounced in 2024–2025, prompting a closer examination using high-temporal-resolution imagery from the Advanced Meteorological Imager (AMI) on the geostationary satellite GEO-KOMPSAT-2A. These observations revealed a clear spike in fire-related signals occurring around and after 4:00 p.m. local time, i.e., outside the typical noon to 2:00 p.m. detection window of the MODIS and VIIRS. A fire detection algorithm exploiting the fire-sensitive shortwave-infrared 3.8 μm signal and its contrast to 11.2 μm infrared observations is designed to adopt AMI observations and applied to its multi-year observations (2019–2025). The resulting fire dataset unambiguously shows a gradual shift in stubble burning activity toward the late afternoon hours beginning in 2022 which is underreported by polar-orbiting satellites. The orbital drift of NASA’s MODIS sensor on the Aqua platform allows detection of some of the gradually shifting fires during afternoon hours, but the MODIS still misses a large number of fires occurring around and after 4 p.m. The AMI’s relatively coarse spatial resolution (~4 km), a consequence of its slant viewing geometry over NW India, imposes inherent limitations on quantifying the full extent of fire occurrences. The operational air quality forecasting models currently assimilate satellite fire detections predominantly captured during early afternoon overpasses of the MODIS and VIIRS. The temporal shift in fire activity complicates such forecast, leading to a substantial underestimation of emissions. Intense stubble burning and the resulting air pollution highlight the need for effective crop residue management practices for mitigating the frequency of open biomass burning and thereby reducing episodic degradation of air quality and its associated public health and economic impacts.

post-monsoon stubble burning; northwestern India;↗

Interannual Variability of the North Polar Vortex in the Lower Stratosphere During the UARS Mission

Northern winters since the 1991 launch of UARS are compared to earlier years (1978 -1991) with respect to the potential for formation of Polar Stratospheric Clouds and for isolation of the north polar vortex. Daily NMC temperature minima at 465 K late in the winter of 1993-94 and again in December 1994 were the lowest values experienced during these times of the year (since 1978). Northern PV gradients were unusually strong in 1991-92 prior to late January and throughout the winter in both 1992-93 and 1994-95. Of all northern winters since 1978, 1994-95 with its early extended cold spell and persistently strong PV gradients most resembled the Antarctic winter lower stratosphere. Even so, temperatures was never as low, nor was the polar vortex as large, as during a typical southern winter. Judged by daily temperature minima and PV gradients at 465 K, meteorological conditions in the Arctic winter lower stratosphere during the UARS period were more conducive to vortex ozone loss by heterogeneous chemistry than in most previous winters since 1978-79.

Zurek, R. W.↗

Using Total Lightning Observations to Enhance Lightning Safety

Lightning is often the underrated threat faced by the public when it comes to dangerous weather phenomena. Typically, larger scale events such as floods, hurricanes, and tornadoes receive the vast majority of attention by both the general population and the media. This comes from the fact that these phenomena are large, longer lasting, can impact a large swath of society at one time, and are dangerous events. The threat of lightning is far more isolated on a case by case basis, although millions of cloud-to-ground lightning strikes hit this United States each year. While attention is given to larger meteorological events, lightning is the second leading cause of weather related deaths in the United States. This information raises the question of what steps can be taken to improve lightning safety. Already, the meteorological community s understanding of lightning has increased over the last 20 years. Lightning safety is now better addressed with the National Weather Service s access to the National Lightning Detection Network data and enhanced wording in their severe weather warnings. Also, local groups and organizations are working to improve public awareness of lightning safety with easy phrases to remember, such as "When Thunder Roars, Go Indoors." The impacts can be seen in the greater array of contingency plans, from airports to sports stadiums, addressing the threat of lightning. Improvements can still be made and newer technologies may offer new tools as we look towards the future. One of these tools is a network of sensors called a lightning mapping array (LMA). Several of these networks exist across the United States. NASA s Short-term Prediction Research and Transition Center (SPoRT), part of the Marshall Spaceflight Center, has access to three of these networks from Huntsville, Alabama, the Kennedy Space Center, and Washington D.C. The SPoRT program s mission is to help transition unique products and observations into the operational forecast environment. SPoRT has been collaborating with the Huntsville National Weather Service (NWS) Office since 2003 and has since included several other offices to better implement LMA observations into real-time applications. Much of that work has focused on the LMA s ability to detect intra-cloud lightning in addition to cloud-to-ground lightning strikes. Combined, these observations are called total lightning. With total lightning observations, NWS offices can enhance their situational awareness and improve severe weather warnings. Just as importantly, the observed intra-cloud flashes often precede the first cloud-to-ground strike by a few minutes. SPoRT and its partner NWS offices are working to develop visualizations and applications to better utilize these data. However, there is a drawback. The LMAs have a short range of no more than 200 km. This is being addressed with the next generation geostationary satellite, GOES-R, which will boast the Geostationary Lightning Mapper (GLM). SPoRT, in conjunction with NOAA s GOES-R Proving Ground, is working to prepare the end user community for the GLM era using the LMA observations as a demonstration tool. Working collaboratively with our NWS partners, SPoRT is working to determine how best to integrate these future observations to improve both severe storm warnings and lightning safety.

Stano, Geoffrey T.↗

The effect of latent heat release on synoptic-to-planetary wave interactions and its implication for satellite observations: Theoretical modeling

The project objectives are to develop process models to investigate the interaction of planetary and synoptic-scale waves including the effects of latent heat release (precipitation), nonlinear dynamics, physical and boundary-layer processes, and large-scale topography; to determine the importance of latent heat release for temporal variability and time-mean behavior of planetary and synoptic-scale waves; to compare the model results with available observations of planetary and synoptic wave variability; and to assess the implications of the results for monitoring precipitation in oceanic-storm tracks by satellite observing systems. Researchers have utilized two different models for this project: a two-level quasi-geostrophic model to study intraseasonal variability, anomalous circulations and the seasonal cycle, and a 10-level, multi-wave primitive equation model to validate the two-level Q-G model and examine effects of convection, surface processes, and spherical geometry. It explicitly resolves several planetary and synoptic waves and includes specific humidity (as a predicted variable), moist convection, and large-scale precipitation. In the past year researchers have concentrated on experiments with the multi-level primitive equation model. The dynamical part of that model is similar to the spectral model used by the National Meteorological Center for medium-range forecasts. The model includes parameterizations of large-scale condensation and moist convection. To test the validity of results regarding the influence of convective precipitation, researchers can use either one of two different convective schemes in the model, a Kuo convective scheme or a modified Arakawa-Schubert scheme which includes downdrafts. By choosing one or the other scheme, they can evaluate the impact of the convective parameterization on the circulation. In the past year researchers performed a variety of initial-value experiments with the primitive-equation model. Using initial conditions typical of climatological winter conditions, they examined the behavior of synoptic and planetary waves growing in moist and dry environments. Surface conditions were representative of a zonally averaged ocean. They found that moist convection associated with baroclinic wave development was confined to the subtropics.

Branscome, Lee E.↗

Bias Reduction in Short Records of Satellite Soil Moisture

Although surface soil moisture data from different sources (satellite retrievals, ground measurements, and land model integrations of observed meteorological forcing data) have been shown to contain consistent and useful information in their seasonal cycle and anomaly signals, they typically exhibit very different mean values and variability. These biases pose a severe obstacle to exploiting the useful information contained in satellite retrievals through data assimilation. A simple method of bias removal is to match the cumulative distribution functions (cdf) of the satellite and model data. However, accurate cdf estimation typically requires a long record of satellite data. We demonstrate here that by wing spatial sampling with a 2 degree moving window we can obtain local statistics based on a one-year satellite record that are a good approximation to those that would be derived from a much longer time series. This result should increase the usefulness of relatively short satellite data records.

Reichle, Rolf H.↗

Use of Polarimetric Radar Measurements to Constrain Simulated Convective Cell Evolution: A Pilot Study with Lagrangian Tracking

To probe the potential value of a radar-driven field campaign to constrain simulation of isolated convection subject to a strong aerosol perturbation, convective cells observed by the operational KHGX weather radar in the vicinity of Houston, Texas, are examined individually and statistically. Cells observed in a single case study of onshore flow conditions during July 2013 are first examined and compared with cells in a regional model simulation. Observed and simulated cells are objectively identified and tracked from observed or calculated positive specific differential phase (K(sub DP)) above the melting level, which is related to the presence of supercooled liquid water. Several observed and simulated cells are subjectively selected for further examination. Below the melting level, we compare sequential cross sections of retrieved and simulated raindrop size distribution parameters. Above the melting level, we examine time series of KDP and radar differential reflectivity (Z(sub DR)) statistics from observations and calculated from simulated supercooled rain properties, alongside simulated vertical wind and supercooled rain mixing ratio statistics. Results indicate that the operational weather radar measurements offer multiple constraints on the properties of simulated convective cells, with substantial value added from derived K(sub DP) and retrieved rain properties. The value of collocated three-dimensional lightning mapping array measurements, which are relatively rare in the continental US, supports the choice of Houston as a suitable location for future field studies to improve the simulation and understanding of convective updraft physics. However, rapid evolution of cells between routine volume scans motivates consideration of adaptive scan strategies or radar imaging technologies to amend operational weather radar capabilities. A 3-year climatology of isolated cell tracks, prepared using a more efficient algorithm, yields additional relevant information. Isolated cells are found within the KHGX domain on roughly 40 % of days year-round, with greatest concentration in the northwest quadrant, but roughly 5-fold more cells occur during June through September. During this enhanced occurrence period, the cells initiate following a strong diurnal cycle that peaks in the early afternoon, typically follow a south-to-north flow, and dissipate within 1 h, consistent with the case study examples. Statistics indicate that ∼ 150 isolated cells initiate and dissipate within 70 km of the KHGX radar during the enhanced occurrence period annually, and roughly 10 times as many within 200 km, suitable for multi-instrument Lagrangian observation strategies. In addition to ancillary meteorological and aerosol measurements, robust vertical wind speed retrievals would add substantial value to a radar-driven field campaign.

Fridlind, Ann M.↗

The latitude dependence of the variance of zonally averaged quantities

Geometric characteristics of the spherical earth are shown to be responsible for the increase of variance with latitude of zonally averaged meteorological statistics. An analytic model is constructed to display the effect of a spherical geometry on zonal averages, employing a sphere labeled with radial unit vectors in a real, stochastic field expanded in complex spherical harmonics. The variance of a zonally averaged field is found to be expressible in terms of the spectrum of the vector field of the spherical harmonics. A maximum variance is then located at the poles, and the ratio of the variance to the zonally averaged grid-point variance, weighted by the cosine of the latitude, yields the zonal correlation typical of the latitude. An example is provided for the 500 mb level in the Northern Hemisphere compared to 15 years of data. Variance is determined to increase north of 60 deg latitude.

North, G. R.↗

The MST Radar Technique

The past ten year have witnessed the development of a new radar technique to examine the structure and dynamics of the atmosphere between roughly 1 to 100 km on a continuous basis. The technique is known as the MST (for Mesosphere-Stratosphere-Troposphere) technique and is usable in all weather conditions, being unaffected by precipitation or cloud cover. MST radars make use of scattering from small scale structure in the atmospheric refractive index, with scales of the order of one-half the radar wavelength. Pertinent scale sizes for middle atmospheric studies typically range between a fraction of a meter and a few meters. The structure itself arises primarily from atmospheric turbulence. The technique is briefly described along with the meteorological parameters it measures.

Balsley, B. B.↗

Analysis of Rapidly Developing Low Cloud Ceilings in a Stable Environment

This report describes the work done by the Applied Meteorology Unit (AMU) in developing a database of days that experienced rapid (< 90 minutes) low cloud formation in a stable atmosphere, resulting in ceilings at the Shuttle Landing Facility (TTS) that violated Space Shuttle Flight Rules (FR). The meteorological conditions favoring the rapid formation of low ceilings include the presence of any inversion below 8000 ft, high relative humidity beneath the inversion, and a clockwise turning of the winds from the surface to the middle troposphere (approx. 15000 ft). The AMU compared and contrasted the atmospheric and thermodynamic conditions between days with rapid low ceiling formation and days with low ceiling resulting from other mechanism. The AMU found that the vertical wind profile is the probable discerning factor between the rapidly-forming ceiling days and other low ceiling days at TTS. Most rapidly-developing low ceiling days had a clockwise turning of the winds with height, whereas other low ceiling days typically had a counter-clockwise turning of the winds with height or negligible vertical wind shear. Forecasters at the Space Meteorology Group (SMG) issue 30 to 90 minute forecasts for low cloud ceilings at TTS to support Space Shuttle landings. Mission verification statistics have shown ceilings to be the number one forecast challenge. More specifically, forecasters at SMG are concerned with any rapidly developing clouds ceilings below 8000 ft in a stable, capped thermodynamic environment, Therefore, the AMU was tasked to examine archived events of rapid stable cloud formation resulting in ceilings below 8000 ft, and document the atmospheric regimes favoring this type of cloud development. The AMU examined the cool season months of November to March during the years of 1993-2003 for days that had low-level inversions and rapid, stable low cloud formation that resulted in ceilings violating the Space Shuttle FR. The AMU wrote and modified existing code to identify inversions from the morning Cape Canaveral, FL rawinsonde (XMR) during the cool season and output pertinent sounding information. They parsed all days with cloud ceilings below 8000 ft at TTS, forming a database of possible rapidly-developing low ceiling events. Days with precipitation or noticeable fog bum-off situations were excluded from the database. Only the daytime hours were examined for possible ceiling development events since low clouds are easier to diagnose with visible satellite imagery. Follow-on work would expand the database to include nighttime cases, using a special enhancement of the infrared imagery for identifying areas of low clouds. The report presents two sample cases of rapidly-developing low cloud ceilings. These cases depict the representative meteorological and thermodynamic characteristics of such events. The cases also illustrate how quickly the cloud decks can develop, sometimes forming in 30 minutes or less. The report also summarizes the composite meteorological conditions for 20 event days with rapid low cloud ceiling formation and 48 non-events days consisting of advection or widespread low cloud ceilings. The meteorological conditions were quite similar for both the event and non-event days, since both types of days experienced low cloud ceilings. Both types of days had a relatively moist environment beneath the inversion based below 8000 ft. In the 20 events identified, de onset of low ceilings occurred between 1200-1800 UTC in every instance. The distinguishing factor between the event and non-event days appears to be the vertical wind profile in the XMR sounding. Eighty-five percent of the event days had a clockwise turning of the winds with height in the lower to middle troposphere whereas 83% of the non-events had a counter-clockwise turning of the winds with height or negligible vertical wind shear. A clockwise turning of the winds with height indicates a warm advection regime, which supports large-scale rising motn and possible cloud formation. Meanwhile, a counter-clockwise turning of the winds with height indicates cold advection or sinking motion in a post-cold frontal environment.

Wheeler, Mark M.↗

Global Modeling of Tropospheric Chemistry with Assimilated Meteorology: Model Description and Evaluation

We present a first description and evaluation of GEOS-CHEM, a global three-dimensional (3-D) model of tropospheric chemistry driven by assimilated meteorological observations from the Goddard Earth Observing System (GEOS) of the NASA Data Assimilation Office (DAO). The model is applied to a 1-year simulation of tropospheric ozone-NOx-hydrocarbon chemistry for 1994, and is evaluated with observations both for 1994 and for other years. It reproduces usually to within 10 ppb the concentrations of ozone observed from the worldwide ozonesonde data network. It simulates correctly the seasonal phases and amplitudes of ozone concentrations for different regions and altitudes, but tends to underestimate the seasonal amplitude at northern midlatitudes. Observed concentrations of NO and peroxyacetylnitrate (PAN) observed in aircraft campaigns are generally reproduced to within a factor of 2 and often much better. Concentrations of HNO3 in the remote troposphere are overestimated typically by a factor of 2-3, a common problem in global models that may reflect a combination of insufficient precipitation scavenging and gas-aerosol partitioning not resolved by the model. The model yields an atmospheric lifetime of methylchloroform (proxy for global OH) of 5.1 years, as compared to a best estimate from observations of 5.5 plus or minus 0.8 years, and simulates H2O2 concentrations observed from aircraft with significant regional disagreements but no global bias. The OH concentrations are approximately 20% higher than in our previous global 3-D model which included an UV-absorbing aerosol. Concentrations of CO tend to be underestimated by the model, often by 10-30 ppb, which could reflect a combination of excessive OH (a 20% decrease in model OH could be accommodated by the methylchloroform constraint) and an underestimate of CO sources (particularly biogenic). The model underestimates observed acetone concentrations over the South Pacific in fall by a factor of 3; a missing source from the ocean may be implicated.

Bey, Isabelle↗

Estimating Future Changes of Energy Demand for Heating and Cooling Buildings at NASA Centers GC23J-1198

With its unique and trusted earth observations, NASA is a critical source in informing decisions that will help achieve the U.S. goal of Net-Zero Greenhouse Gas (GHG) Emissions by 2050. NASA’s Prediction of Worldwide Energy Resource (POWER) project facilitates the use of NASA Earth Science data holdings within the energy, agricultural, and building heating/cooling design industries. POWER packages solar and meteorological data from several NASA projects in a user friendly GIS-enabled web services system (https://power.larc.nasa.gov). As part of the development of new data products to support the energy and building heating/cooling design communities, we estimate the changes in energy required to heat and cool buildings in the future climate at 14 different NASA site locations spread throughout the continental United States, as projected by CMIP6 climate models under different emissions scenarios. Bias-corrected downscaled time series of meteorological variables are taken from NASA Earth Exchange (NEX) Global Daily Downscaled Projections (GDDP-CMIP6) downscaled climate model data. The spread between the different model projections is accounted for by analyzing both the ensemble average of 22 CMIP6 models and 6 representative models with different climate sensitivities and different interannual variability. Changes in energy use are estimated in two ways. First, changes in the total annual heating and cooling degree days (HDD and CDD, respectively) are calculated relative to the current climate. This is done at all 14 sites. Second, the downscaled time series are used as inputs into RETScreen(R), a clean energy management decision tool, to give an estimate of heating/cooling energy use for a typical office building. We use this estimation method with model data at Langley Research Center. In the next 50 years, the annual total of HDD (CDD) is projected to decrease by 8-38% (increase by 5-28%) at all sites, with the increase in CDD typically a larger magnitude the decrease in HDD. For a typical small office building at Langley Research Center, the amount of energy needed to cool increases by 33-54% and the amount of energy to heat decreases by 29-40%. POWER is working to develop long term climate data services based on these results to include in future data products to provide to users.

Bradley M. Hegyi↗

Seasonal Soil Freeze/Thaw Variability Across North America via Ensemble Land Surface Modeling

Land surface modeling provides the opportunity to investigate winter soil processes such as soil freezing and thawing over large domains. However, the variability in simulated winter soil characteristics among land surface models and forcing datasets is not well understood. In this study, a nine-member ensemble was employed to characterize the spatial and inter-annual variability of three winter soil characteristics, annual number of frozen days, annual minimum temperature and annual number of freeze-thaw cycles over North America during the water years 2010–2016. The ensemble, including three land surface model (JULES, Noah2.7.1 and Noah-MP) and three forcing datasets (ECMWF, GDAS and MERRA2), was developed through the Snow Ensemble Uncertainty Project (SEUP). In many regions, there was remarkably good agreement across the ensemble for the winter soil temperatures. However, the variability among the ensemble's annual number of frozen days, as quantified with standard deviation, exceeded 150 days at the northern Pacific coastline. While the differences among the ensemble members' annual minimum temperature were typically <3 °C, the differences exceeded 6 °C north of 50°N. The ensemble members generally agreed within one to three freeze-thaw cycles in mid latitude regions, Alaska, and the south and west coasts of the USA. High variability among the ensemble, more than six freeze-thaw cycles, occurred in the Great Plains, northern Pacific coastline, and along the Appalachian Mountains. Differences in winter soil temperature characteristics were more apparent among the LSMs rather than the meteorological forcing datasets. Except for maritime regions, the Noah2.7.1 members had the fewest annual number of frozen days and the least number of freeze-thaw cycles. Noah2.7.1 also had the coldest annual minimum temperatures except for ephemeral regions. Comparisons between in-situ observations and the SEUP estimates of winter soil characteristics revealed that the modeled frozen period was much longer than observed, that the modeled annual minimum temperatures were much colder than observed, and the modeled freeze-thaw cycles occurred more frequently than observed. Excluding the high latitude sites, the observed frozen period is less than two months with the minimum temperature above − 5 °C at most of the studied sites, while the ensemble members simulated, on average, a four month frozen period with − 10 °C minimum temperature. Errors in capturing snow during the accumulation period appear to impact differences in modeled versus observed soil temperature throughout the entire winter.

Soil freeze/thaw↗

Evaluation of MERRA-2-Based Ozone Profile Simulations with the Global Ozonesonde Network

Chemical transport model (CTM) hindcasts of ozone (O3) are useful for filling in observational gaps and providing context for observed O3 variability and trends. We use global networks of ozonesonde stations to evaluate the O3 profiles in two simulations running versions of the NASA Global Modeling Initiative (GMI) chemical mechanism. Both simulations are tied to the NASA Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2) meteorological reanalysis: 1) The GMI CTM, and 2) The MERRA-2 GMI Replay (M2 GMI). Both simulations start in 1980, and are compared against >50,000 ozonesonde profiles from 37 global stations from the tropics to the poles. The comparisons allow us to evaluate how the Replay technique affects modeled O3 distribution, how an updated chemical mechanism in the GMI CTM affects simulated tropospheric O3 amounts, and how observed O3 distributions compare to the full set of model output. In general, M2 GMI O3 is ~10% higher than in the GMI CTM, and shows global near-surface and tropical upper troposphere/lower stratosphere (UT/LS) high biases. The updated chemical mechanism in the GMI CTM reduces these high biases. Both simulations show similar negative biases in tropical free-tropospheric O3, especially during typical biomass burning seasons. The simulations are highly-correlated with ozonesonde measurements, particularly in the UT/LS (r > 0.8), showing the ability of MERRA-2 to capture tropopause height variations. Both simulations show improved correlations with ozonesonde data and smaller O3 biases in recent years. We expect to use the sonde/model comparisons to diagnose causes of disagreement and to gauge the feasibility of calculating multidecadal O3 trends from the model output.

Stauffer, Ryan↗

Evaluation of MERRA-2-Based Ozone Profile Simulations with the Global Ozonesonde Network

Chemical transport model (CTM) hindcasts of ozone (O3) are useful for filling in observational gaps and providing context for observed O3 variability and trends. We use global networks of ozonesonde stations to evaluate the O3 profiles in two simulations running versions of the NASA Global Modeling Initiative (GMI) chemical mechanism. Both simulations are tied to the NASA Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2) meteorological reanalysis: 1) The GMI CTM, and 2) The MERRA-2 GMI Replay (M2 GMI). Both simulations start in 1980, and are compared against >50,000 ozonesonde profiles from 37 global stations from the tropics to the poles. The comparisons allow us to evaluate how the Replay technique affects modeled O3 distribution, how an updated chemical mechanism in the GMI CTM affects simulated tropospheric O3 amounts, and how observed O3 distributions compare to the full set of model output. In general, M2 GMI O3 is ~10% higher than in the GMI CTM, and shows global near-surface and tropical upper troposphere/lower stratosphere (UT/LS) high biases. The updated chemical mechanism in the GMI CTM reduces these high biases. Both simulations show similar negative biases in tropical free-tropospheric O3, especially during typical biomass burning seasons. The simulations are highly-correlated with ozonesonde measurements, particularly in the UT/LS (r > 0.8), showing the ability of MERRA-2 to capture tropopause height variations. Both simulations show improved correlations with ozonesonde data and smaller O3 biases in recent years. We expect to use the sonde/model comparisons to diagnose causes of disagreement and to gauge the feasibility of calculating multidecadal O3 trends from the model output.

Stauffer, Ryan M.↗

Diurnal Waves Forced By Horizontal Convergence of Near-Surface Winds on Mars

Boundary-layer winds near steeply sloped topography excite a distinct type of atmospheric wave. Basic properties of the waves were determined through analysis of 7 Mars years of temperature profiles from the Mars Climate Sounder on Mars Reconnaissance Orbiter. A more complete understanding of the waves was derived from a simulation with the NASA Ames Mars General Circulation Model. This type of wave appears in the tropics in the vicinity of conspicuous surface features, such as large volcanoes and the rim of Isidis Planitia. The amplitude is typically 5–10 K, the vertical wavelength is about 40 km, the zonal wavelength is about 1400 km, and the zonal phase speed relative to the surface is about 25 m s−1. The wave forcing comes from horizontal convergence of the near-surface winds, which generates a daily surge in the vertical winds. These convergence zones were not previously recognized as an important source of atmospheric waves on Mars. Strong forcing occurs only within a limited range of local time, producing wave packets (diurnal pulses) rather than steady oscillations. Waves are present throughout the year and their properties vary with season. Near the solstices, the wave amplitude is largest in the tropics of the summer hemisphere, where the convergence zones are stronger than in the winter hemisphere. The diurnal-mean zonal winds, which vary with latitude and season, control the direction of propagation of the predominant waves through filtering by critical levels.

Mars↗

Assessing the Ability of Instantaneous Aircraft and Sonde Measurements to Characterize Climatological Means and Long-Term Trends in Tropospheric Composition

Over four decades of measurements exist that sample the 3-D composition of reactive trace gases in the troposphere from approximately weekly ozone sondes, instrumentation on civil aircraft, and individual comprehensive aircraft field campaigns. An obstacle to using these data to evaluate coupled chemistry-climate models (CCMs)the models used to project future changes in atmospheric composition and climateis that exact space-time matching between model fields and observations cannot be done, as CCMs generate their own meteorology. Evaluation typically involves averaging over large spatiotemporal regions, which may not reflect a true average due to limited or biased sampling. This averaging approach generally loses information regarding specific processes. Here we aim to identify where discrete sampling may be indicative of long-term mean conditions, using the GEOS-Chem global chemical-transport model (CTM) driven by the MERRA reanalysis to reflect historical meteorology from 2003 to 2012 at 2o by 2.5o resolution. The model has been sampled at the time and location of every ozone sonde profile available from the Would Ozone and Ultraviolet Radiation Data Centre (WOUDC), along the flight tracks of the IAGOSMOZAICCARABIC civil aircraft campaigns, as well as those from over 20 individual field campaigns performed by NASA, NOAA, DOE, NSF, NERC (UK), and DLR (Germany) during the simulation period. Focusing on ozone, carbon monoxide and reactive nitrogen species, we assess where aggregates of the in situ data are representative of the decadal mean vertical, spatial and temporal distributions that would be appropriate for evaluating CCMs. Next, we identically sample a series of parallel sensitivity simulations in which individual emission sources (e.g., lightning, biogenic VOCs, wildfires, US anthropogenic) have been removed one by one, to assess where and when the aggregated observations may offer constraints on these processes within CCMs. Lastly, we show results of an additional 31-year simulation from 1980-2010 of GEOS-Chem driven by the MACCity emissions inventory and MERRA reanalysis at 4o by 5o. We sample the model at every WOUDC sonde and flight track from MOZAIC and NASA field campaigns to evaluate which aggregate observations are statistically reflective of long-term trends over the period.

Atmospheric composition↗

Ice surface temperature retrieval from AVHRR, ATSR, and passive microwave satellite data: Algorithm development and application

One essential parameter used in the estimation of radiative and turbulent heat fluxes from satellite data is surface temperature. Sea and land surface temperature (SST and LST) retrieval algorithms that utilize the thermal infrared portion of the spectrum have been developed, with the degree of success dependent primarily upon the variability of the surface and atmospheric characteristics. However, little effort has been directed to the retrieval of the sea ice surface temperature (IST) in the Arctic and Antarctic pack ice or the ice sheet surface temperature over Antarctica and Greenland. The reason is not one of methodology, but rather our limited knowledge of atmospheric temperature, humidity, and aerosol vertical, spatial and temporal distributions, the microphysical properties of polar clouds, and the spectral characteristics of snow, ice, and water surfaces. Over the open ocean the surface is warm, dark, and relatively homogeneous. This makes SST retrieval, including cloud clearing, a fairly straightforward task. Over the ice, however, the surface within a single satellite pixel is likely to be highly heterogeneous, a mixture of ice of various thicknesses, open water, and snow cover in the case of sea ice. Additionally, the Arctic is cloudy - very cloudy - with typical cloud cover amounts ranging from 60-90 percent. There are few observations of cloud cover amounts over Antarctica. The goal of this research is to increase our knowledge of surface temperature patterns and magnitudes in both polar regions, by examining existing data and improving our ability to use satellite data as a monitoring tool. Four instruments are of interest in this study: the AVHRR, ATSR, SMMR, and SSM/I. Our objectives are as follows. Refine the existing AVHRR retrieval algorithm defined in Key and Haefliger (1992; hereafter KH92) and applied elsewhere. Develop a method for IST retrieval from ATSR data similar to the one used for SST. Further investigate the possibility of estimating surface temperature from passive microwave data (in conjunction with AVHRR clear sky samples) through the use of 'effective emissivities' and physical relationships between skin temperature and subsurface temperature. Use the general method outlined in MK93 to calculate a 12-year record of clear sky equivalent surface temperatures, or possibly all-sky snow-ice interface physical temperatures, from SMMR and SSM/I, compare these temperatures to climatologies, ECMWF modeled surface temperatures, and surface temperatures predicted by a 2-D ice model. And intercompare several ice surface retrieval methods and validate them against ground measurements from the Swiss Camp on the Greenland ice sheet. Additionally, we intend to develop a surface temperature product based on AVHRR data and possibly blended with drifting buoy and meteorological station temperatures.

Key, Jeff↗

Tropical Pacific moisture variability

The objectives are to describe synoptic scale variability of moisture over the tropical Pacific Ocean and the systems leading to this variability; implement satellite analysis procedures in support of this effort, and to incorporate additional satellite information into operational analysis forecast systems at the National Meteorological Center (NMC). Composite satellite radiance patterns describe features detectable well before the development of synoptic scale tropical plumes. These typical features were extracted from historical files of Tiros Operational Vertical Sounder (TOVS) radiance observations for a pair of tropical plumes which developed during January 1989. Signals were inserted into the NMC operational medium range forecast model and a suite of model integrations were conducted. Many of the 48 h model errors of the historical forecasts were eliminated by the inclusion of more complete satellite observations. Three studies in satellite radiance analysis progressed. An analysis which blended TOVS moisture channels, OLR observations and European Center for Medium Weather Forecasts (ECMWF) model analysis to generate fields of total precipitable water comparable to those estimated from Scanning Multichannel Microwave Radiometer (SMMR) mu-wave observations. This study demonstrated that a 10 y climatology of precipitable water over the oceans is feasible, using available infrared observations (OLR and TOVS) and model analysis (ECMWF, NMC or similar quality). The estimates are sensitive to model quality and the estimating model must be updated with operational model changes. Coe developed a set of tropical plume and ITCZ composites from TOVS observations, and from NMC and ECMWF analyses which had been passed through a radiative transfer model to simulate TOVS radiances. The composites have been completed as well as many statistical diagnostics of individual TOVS channels. Analysis of the computations is commencing. Chung has initiated a study of the differences between TOVS observed vapor structure during El Nino Southern Oscillation (ENSO) (1983) and non-ENSO (1984) years. Preliminary diagnosis demonstrates gross moisture changes between warm and cold sea surface temperature episodes.

Mcguirk, James P.↗