NRLMSIS® 2.0 is an empirical atmospheric model that extends from the ground to the exobase and describes the average observed behavior of temperature, eight species densities, and mass density via a parametric analytic formulation. The model inputs are location, day of year, time of day, solar activity, and geomagnetic activity. NRLMSIS 2.0 is a major, reformulated upgrade of the previous version, NRLMSISE‐00. The model now couples thermospheric species densities to the entire column, via an effective mass profile that transitions each species from the fully mixed region below ~70 km altitude to the diffusively separated region above ~200 km. Other changes include the extension of atomic oxygen down to 50 km and the use of geopotential height as the internal vertical coordinate. We assimilated extensive new lower and middle atmosphere temperature, O, and H data, along with global average thermospheric mass density derived from satellite orbits, and we validated the model against independent samples of these data. In the mesosphere and below, residual biases and standard deviations are considerably lower than NRLMSISE‐00. The new model is warmer in the upper troposphere and cooler in the stratosphere and mesosphere. In the thermosphere, N₂ and O densities are lower in NRLMSIS 2.0; otherwise, the NRLMSISE‐00 thermosphere is largely retained. Future advances in thermospheric specification will likely require new in situ mass spectrometer measurements, new techniques for species density measurement between 100 and 200 km, and the reconciliation of systematic biases among thermospheric temperature and composition data sets, including biases attributable to long‐term changes.

[1] The variability of stratospheric aerosol loading between 1985 and 2010 is explored with measurements from SAGE II, CALIPSO, GOMOS/ENVISAT, and OSIRIS/Odin space‐based instruments. We find that, following the 1991 eruption of Mount Pinatubo, stratospheric aerosol levels increased by as much as two orders of magnitude and only reached “background levels” between 1998 and 2002. From 2002 onwards, a systematic increase has been reported by a number of investigators. Recently, the trend, based on ground‐based lidar measurements, has been tentatively attributed to an increase of SO₂ entering the stratosphere associated with coal burning in Southeast Asia. However, we demonstrate with these satellite measurements that the observed trend is mainly driven by a series of moderate but increasingly intense volcanic eruptions primarily at tropical latitudes. These events injected sulfur directly to altitudes between 18 and 20 km. The resulting aerosol particles are slowly lofted into the middle stratosphere by the Brewer‐Dobson circulation and are eventually transported to higher latitudes.

We report the first measurement of nighttime atmospheric temperature and horizontal wind profiles in the lower thermosphere up to 140 km with the Na lidar at Andes Lidar Observatory in Cerro Pachón, Chile (30.25°S, 70.74°W), when enhanced thermospheric Na was observed. Temperature and horizontal wind were derived up to 140 km using various resolutions, with the lowest resolution of about 2.7 hr and 15 km above 130 km. Thus, the measurements span 60 km in vertical, more than double the traditional 25 km. On the night of 17 April 2015, the horizontal wind magnitude in the thermosphere exceeds 150 ms⁻¹, consistent with past rocket measurements. The meridional wind shows a clear transition from the diurnal‐tide‐dominant mesopause to the semidiurnal‐tide‐dominant lower thermosphere. A lidar with a 100 times the power aperture product will be able to measure wind and temperature above 160 km and cover longer time span, providing key measurements for the study of atmosphere‐space interactions in this region.

The O(¹S) (green line) night airglow emission in response to atmospheric gravity wave (AGW) perturbations was simulated with a linear, one‐dimensional model. The results were combined with previously modeled O₂(b, 0–1) atmospheric band and OH Meinel band emission response (Liu and Swenson, 2003) to derive amplitude and phase relations among multiple airglow layers in response to gravity waves with various intrinsic parameters and damping rates (β). The simulations show that the vertical profile of the standard deviation of the perturbed green line volume emission rate (VER) has a centroid altitude that is 3 km lower and a full‐width‐half‐maximum 2.1 km smaller than the unperturbed VER profile, similar to findings for the OH and O₂(b) band layers. Relative phase differences and amplitudes of vertically propagating waves can be deduced from zenith observations of the layers. Airglow weighted responses to waves are related through a cancellation factor (CF) for both layer intensity and temperature. The vertical wavelength can be deduced from relative phase information of three airglow layers separated in altitude. The vertical flux of horizontal momentum associated with gravity waves is deduced from intrinsic wave parameters. Wave damping versus altitude is used to deduce the flux divergence and local accelerations resulting from dissipative waves. The simulations are useful in calculating wave information and wave effects on the atmosphere from multiwavelength, zenith airglow observations.

The Aerospace Corporation's Nightglow Imager (ANI) observes nighttime OH emission (near 1.6 μm) every 2 s over an approximate 73° field of view. ANI had previously been used to study instability features seen over Maui. Here we describe observations of instabilities seen from 5 to 8 UT on 24 March 2012 over Cerro Pachon, Chile, and compare them with previous results from Maui, with theory, and with Direct Numerical Simulations (DNS). The atmosphere had reduced stability because of the large negative temperature gradients measured by a Na lidar. Thus, regions of dynamical and convective instabilities are expected to form, depending on the value of the Richardson number. Bright primary instabilities are formed with a horizontal wavelength near 9 km and showed the subsequent formation of secondary instabilities, rarely seen over Maui, consistent with the primaries being dynamical instabilities. The ratio of the primary to secondary horizontal wavelength was greater over Chile than over Maui. After dissipation of the instabilities, smaller‐scale features appeared with sizes in the buoyancy subrange between 1.5 and 6 km. Their size spectra were consistent with the model of Weinstock (1978) if the turbulence is considered to be increasing. The DNS results produce secondary instabilities with sizes comparable to what is seen in the images although their spectra are somewhat steeper than is observed. However, the DNS results also show that after the complete decay of the primary features, scale sizes considerably smaller than 1 km are produced and these cannot be seen by the ANI instrument.

We provide here an overview of, and a summary of results arising from, an extensive experimental campaign (the Spread F Experiment, or SpreadFEx) performed from September to November 2005, with primary measurements in Brazil. The motivation was to define the potential role of neutral atmosphere dynamics, specifically gravity wave motions propagating upward from the lower atmosphere, in seeding Rayleigh-Taylor instability (RTI) and plasma bubbles extending to higher altitudes. Campaign measurements focused on the Brazilian sector and included ground-based optical, radar, digisonde, and GPS measurements at a number of fixed and temporary sites. Related data on convection and plasma bubble structures were also collected by GOES 12, and the GUVI instrument aboard the TIMED satellite. Initial results of our SpreadFEx analyses are described separately by Fritts et al. (2009). Further analyses of these data provide additional evidence of 1) gravity wave (GW) activity near the mesopause apparently linked to deep convection predominantly to the west of our measurement sites, 2) small-scale GWs largely confined to lower altitudes, 3) larger-scale GWs apparently penetrating to much higher altitudes, 4) substantial GW amplitudes implied by digisonde electron densities, and 5) apparent influences of these perturbations in the lower F-region on the formation of equatorial spread F, RTI, and plasma bubbles extending to much higher altitudes. Other efforts with SpreadFEx data have also yielded 6) the occurrence, locations, and scales of deep convection, 7) the spatial and temporal evolutions of plasma bubbles, 8) 2-D (height-resolved) structures in electron density fluctuations and equatorial spread F at lower altitudes and plasma bubbles above, and 9) the occurrence of substantial tidal perturbations to the large-scale wind and temperature fields extending to bottomside F-layer and higher altitudes. Collectively, our various SpreadFEx analyses suggest direct links between deep tropical convection and large GW perturbations at large spatial scales at the bottomside F-layer and their likely contributions to the excitation of RTI and plasma bubbles extending to much higher altitudes.