Vertical ion flow in the polar ionosphere

Outflowing ions from the polar ionosphere fall into two categories: the classical polar wind and the suprathermal ion flows. The flows in both these categories vary a great deal with altitude. The classical polar wind is supersonic at high altitude: at ∼3 RE geocentric, the observed polar wind is H+ dominated and has a Mach number of 2.5–5.1. At 400–600 km, thermal and suprathermal upward O+ ion fluxes frequently occur at the poleward edge of the nightside auroral oval during magnetically active times. Above 500 km, ions are accelerated transverse to the local geomagnetic field. At 1400 km, transversely accelerated ions are frequently observed in winter nights but rarely appear in the summer. In the dayside cleft above ∼2000 km, ions of all species are transversely heated and upwell with significant number and heat fluxes, forming a cleft ion fountain as they convect across the polar cap. Upwelling ions are observed most (least) frequently in the summer (winter). At yet higher altitudes, energetic (>10 eV to several kiloelectron volts) upflowing H+ and O+ ions are frequently observed, their active time occurrence frequency being as high as 0.7 at auroral latitudes and 0.3 in the polar cap. Their composition, intensity, and angular characteristics vary quantitatively with solar activity, being O+ dominant and more intense near solar maximum. Their resulting ion outflow is dominated by ions below 1 keV and reaches 3.5×10^26 O+ and 7×10^25 H+ ions s^{−1} at magnetically active times (Kp≥5) near solar maximum. In comparison, the estimated polar wind ion outflow at times of moderate solar activity is 7×10^25H+ and 4×10^24 He+ ions s^{−1}. The estimated <10-eV cleft ion fountain flow is 3.8×10^25 O+ and 8.6×10^23 H+ ions s^{−1} near solar maximum.

High-precision measurements of the co-polar correlation coefficient: non-Gaussian errors and retrieval of the dispersion parameter µ in rainfall

The co-polar correlation coefficient (ρhv) has many applications, including hydrometeor classification, ground clutter and melting layer identification, interpretation of ice microphysics and the retrieval of rain drop size distributions (DSDs). However, we currently lack the quantitative error estimates that are necessary if these applications are to be fully exploited. Previous error estimates of ρhv rely on knowledge of the unknown "true" ρhv and implicitly assume a Gaussian probability distribution function of ρhv samples. We show that frequency distributions of ρhv estimates are in fact highly negatively skewed. A new variable: L = -log10(1 - ρhv) is defined, which does have Gaussian error statistics, and a standard deviation depending only on the number of independent radar pulses. This is verified using observations of spherical drizzle drops, allowing, for the first time, the construction of rigorous confidence intervals in estimates of ρhv. In addition, we demonstrate how the imperfect co-location of the horizontal and vertical polarisation sample volumes may be accounted for. The possibility of using L to estimate the dispersion parameter (µ) in the gamma drop size distribution is investigated. We find that including drop oscillations is essential for this application, otherwise there could be biases in retrieved µ of up to ~8. Preliminary results in rainfall are presented. In a convective rain case study, our estimates show µ to be substantially larger than 0 (an exponential DSD). In this particular rain event, rain rate would be overestimated by up to 50% if a simple exponential DSD is assumed.

The counter-propagating Rossby-wave perspective on baroclinic instability. Part III: Primitive-equation disturbances on the sphere

Baroclinic instability of perturbations described by the linearized primitive quations, growing on steady zonal jets on the sphere, can be understood in terms of the interaction of pairs of counter-propagating Rossby waves (CRWs). The CRWs can be viewed as the basic components of the dynamical system where the Hamiltonian is the pseudoenergy and each CRW has a zonal coordinate and pseudomomentum. The theory holds for adiabatic frictionless flow to the extent that truncated forms of pseudomomentum and pseudoenergy are globally conserved. These forms focus attention on Rossby wave activity. Normal mode (NM) dispersion relations for realistic jets are explained in terms of the two CRWs associated with each unstable NM pair. Although derived from the NMs, CRWs have the conceptual advantage that their structure is zonally untilted, and can be anticipated given only the basic state. Moreover, their zonal propagation, phase-locking and mutual interaction can all be understood by ‘PV-thinking’ applied at only two ‘home-bases’—potential vorticity (PV) anomalies at one home-base induce circulation anomalies, both locally and at the other home-base, which in turn can advect the PV gradient and modify PV anomalies there. At short wavelengths the upper CRW is focused in the mid-troposphere just above the steering level of the NM, but at longer wavelengths the upper CRW has a second wave-activity maximum at the tropopause. In the absence of meridional shear, CRW behaviour is very similar to that of Charney modes, while shear results in a meridional slant with height of the air-parcel displacement-structures of CRWs in sympathy with basic-state zonal angular-velocity surfaces. A consequence of this slant is that baroclinically growing eddies (on jets broader than the Rossby radius) must tilt downshear in the horizontal, giving rise to up-gradient momentum fluxes that tend to accelerate the barotropic component of the jet.

The counter-propagating Rossby-wave perspective on baroclinic instability. Part IV: Nonlinear life cycles

Pairs of counter-propagating Rossby waves (CRWs) can be used to describe baroclinic instability in linearized primitive-equation dynamics, employing simple propagation and interaction mechanisms at only two locations in the meridional plane—the CRW ‘home-bases’. Here, it is shown how some CRW properties are remarkably robust as a growing baroclinic wave develops nonlinearly. For example, the phase difference between upper-level and lower-level waves in potential-vorticity contours, defined initially at the home-bases of the CRWs, remains almost constant throughout baroclinic wave life cycles, despite the occurrence of frontogenesis and Rossby-wave breaking. As the lower wave saturates nonlinearly the whole baroclinic wave changes phase speed from that of the normal mode to that of the self-induced phase speed of the upper CRW. On zonal jets without surface meridional shear, this must always act to slow the baroclinic wave. The direction of wave breaking when a basic state has surface meridional shear can be anticipated because the displacement structures of CRWs tend to be coherent along surfaces of constant basic-state angular velocity, U. This results in up-gradient horizontal momentum fluxes for baroclinically growing disturbances. The momentum flux acts to shift the jet meridionally in the direction of the increasing surface U, so that the upper CRW breaks in the same direction as occurred at low levels

Extension of spherical nonparametric estimators to nonisotropic kernels: An oceanographic application

In a recently published paper. spherical nonparametric estimators were applied to feature-track ensembles to determine a range of statistics for the atmospheric features considered. This approach obviates the types of bias normally introduced with traditional estimators. New spherical isotropic kernels with local support were introduced. Ln this paper the extension to spherical nonisotropic kernels with local support is introduced, together with a means of obtaining the shape and smoothing parameters in an objective way. The usefulness of spherical nonparametric estimators based on nonisotropic kernels is demonstrated with an application to an oceanographic feature-track ensemble.

Spherical nonparametric estimators applied to the UGAMP model integration for AMIP

The aim of this paper is essentially twofold: first, to describe the use of spherical nonparametric estimators for determining statistical diagnostic fields from ensembles of feature tracks on a global domain, and second, to report the application of these techniques to data derived from a modern general circulation model. New spherical kernel functions are introduced that are more efficiently computed than the traditional exponential kernels. The data-driven techniques of cross-validation to determine the amount elf smoothing objectively, and adaptive smoothing to vary the smoothing locally, are also considered. Also introduced are techniques for combining seasonal statistical distributions to produce longer-term statistical distributions. Although all calculations are performed globally, only the results for the Northern Hemisphere winter (December, January, February) and Southern Hemisphere winter (June, July, August) cyclonic activity are presented, discussed, and compared with previous studies. Overall, results for the two hemispheric winters are in good agreement with previous studies, both for model-based studies and observational studies.

An improved algorithm for generating global window brightness temperatures from multiple satellite infrared imagery

An improved algorithm for the generation of gridded window brightness temperatures is presented. The primary data source is the International Satellite Cloud Climatology Project, level B3 data, covering the period from July 1983 to the present. The algorithm rakes window brightness, temperatures from multiple satellites, both geostationary and polar orbiting, which have already been navigated and normalized radiometrically to the National Oceanic and Atmospheric Administration's Advanced Very High Resolution Radiometer, and generates 3-hourly global images on a 0.5 degrees by 0.5 degrees latitude-longitude grid. The gridding uses a hierarchical scheme based on spherical kernel estimators. As part of the gridding procedure, the geostationary data are corrected for limb effects using a simple empirical correction to the radiances, from which the corrected temperatures are computed. This is in addition to the application of satellite zenith angle weighting to downweight limb pixels in preference to nearer-nadir pixels. The polar orbiter data are windowed on the target time with temporal weighting to account for the noncontemporaneous nature of the data. Large regions of missing data are interpolated from adjacent processed images using a form of motion compensated interpolation based on the estimation of motion vectors using an hierarchical block matching scheme. Examples are shown of the various stages in the process. Also shown are examples of the usefulness of this type of data in GCM validation.

A Tropospheric Assessment of the ERA-40, NCEP, and JRA-25 Global Reanalyses in the Polar Regions

The reliability of the global reanalyses in the polar regions is investigated. The overview stems from an April 2006 Scientific Committee on Antarctic Research (SCAR) workshop on the performance of global reanalyses in high latitudes held at the British Antarctic Survey. Overall, the skill is much higher in the Arctic than the Antarctic, where the reanalyses are only reliable in the summer months prior to the modern satellite era. In the Antarctic, large circulation differences between the reanalyses are found primarily before 1979, when vast quantities of satellite sounding data started to be assimilated. Specifically for ERA-40, this data discontinuity creates a marked jump in Antarctic snow accumulation, especially at high elevations. In the Arctic, the largest differences are related to the reanalyses depiction of clouds and their associated radiation impacts; ERA-40 captures the cloud variability much better than NCEP1 and JRA-25, but the ERA-40 and JRA-25 clouds are too optically thin for shortwave radiation. To further contrast the reanalyses skill, cyclone tracking results are presented. In the Southern Hemisphere, cyclonic activity is markedly different between the reanalyses, where there are few matched cyclones prior to 1979. In comparison, only some of the weaker cyclones are not matched in the Northern Hemisphere from 1958-2001, again indicating the superior skill in this hemisphere. Although this manuscript focuses on deficiencies in the reanalyses, it is important to note that they are a powerful tool for climate studies in both polar regions when used with a recognition of their limitations.

An objective climatology of the dynamical forcing of polar lows in the Nordic seas

A new objective climatology of polar lows in the Nordic (Norwegian and Barents) seas has been derived from a database of diagnostics of objectively identified cyclones spanning the period January 2000 to April 2004. There are two distinct parts to this study: the development of the objective climatology and a characterization of the dynamical forcing of the polar lows identified. Polar lows are an intense subset of polar mesocyclones. Polar mesocyclones are distinguished from other cyclones in the database as those that occur in cold air outbreaks over the open ocean. The difference between the wet-bulb potential temperature at 700 hPa and the sea surface temperature (SST) is found to be an effective discriminator between the atmospheric conditions associated with polar lows and other cyclones in the Nordic seas. A verification study shows that the objective identification method is reliable in the Nordic seas region. After demonstrating success at identifying polar lows using the above method, the dynamical forcing of the polar lows in the Nordic seas is characterized. Diagnostics of the ratio of mid-level vertical motion attributable to quasi-geostrophic forcing from upper and lower levels (U/L ratio) are used to determine the prevalence of a recently proposed category of extratropical cyclogenesis, type C, for which latent heat release is crucial to development. Thirty-one percent of the objectively identified polar low events (36 from 115) exceeded the U/L ratio of 4.0, previously identified as a threshold for type C cyclones. There is a contrast between polar lows to the north and south of the Nordic seas. In the southern Norwegian Sea, the population of polar low events is dominated by type C cyclones. These possess strong convection and weak low-level baroclinicity. Over the Barents and northern Norwegian seas, the well-known cyclogenesis types A and B dominate. These possess stronger low-level baroclinicity and weaker convection.

A new look at stratospheric sudden warmings. Part III. Polar vortex evolution and vertical structure

The evolution of the Arctic polar vortex during observed major mid-winter stratospheric sudden warmings (SSWs) is investigated for the period 1957-2002, using European Centre for Medium-Range Weather Forecasts (ECMWF) ERA-40 Ertel’s potential vorticity (PV) and temperature fields. Time-lag composites of vertically weighted PV, calculated relative to the SSW onset time, are derived for both vortex displacement SSWs and vortex splitting SSWs, by averaging over the 15 recorded displacement and 13 splitting events. The evolving vertical structure of the polar vortex during a typical SSW of each type is clearly illustrated by plotting an isosurface of the composite PV field, and is shown to be very close to that observed during representative individual events. Results are verified by comparison with an elliptical diagnostic vortex moment technique. For both types of SSW, little variation is found between individual events in the orientation of the developing vortex relative to the underlying topography, i.e. the location of the vortex during SSWs of each type is largely fixed in relation to the Earth’s surface. During each type of SSW, the vortex is found to have a distinctive vertical structure. Vortex splitting events are typically barotropic, with the vortex split occurring near-simultaneously over a large altitude range (20-40 km). In the majority of cases, of the two daughter vortices formed, it is the ‘Siberian’ vortex that dominates over its ‘Canadian’ counterpart. In contrast, displacement events are characterized by a very clear baroclinic structure; the vortex tilts significantly westward with height, so that the top and bottom of the vortex are separated by nearly 180◦ longitude before the upper vortex is sheared away and destroyed.