Global observed changes in daily climate extremes of temperature and precipitation

A suite of climate change indices derived from daily temperature and precipitation data, with a primary focus on extreme events, were computed and analyzed. By setting an exact formula for each index and using specially designed software, analyses done in different countries have been combined seamlessly. This has enabled the presentation of the most up-to-date and comprehensive global picture of trends in extreme temperature and precipitation indices using results from a number of workshops held in data-sparse regions and high-quality station data supplied by numerous scientists world wide. Seasonal and annual indices for the period 1951-2003 were gridded. Trends in the gridded fields were computed and tested for statistical significance. Results showed widespread significant changes in temperature extremes associated with warming, especially for those indices derived from daily minimum temperature. Over 70% of the global land area sampled showed a significant decrease in the annual occurrence of cold nights and a significant increase in the annual occurrence of warm nights. Some regions experienced a more than doubling of these indices. This implies a positive shift in the distribution of daily minimum temperature throughout the globe. Daily maximum temperature indices showed similar changes but with smaller magnitudes. Precipitation changes showed a widespread and significant increase, but the changes are much less spatially coherent compared with temperature change. Probability distributions of indices derived from approximately 200 temperature and 600 precipitation stations, with near-complete data for 1901-2003 and covering a very large region of the Northern Hemisphere midlatitudes (and parts of Australia for precipitation) were analyzed for the periods 1901-1950, 1951-1978 and 1979-2003. Results indicate a significant warming throughout the 20th century. Differences in temperature indices distributions are particularly pronounced between the most recent two periods and for those indices related to minimum temperature. An analysis of those indices for which seasonal time series are available shows that these changes occur for all seasons although they are generally least pronounced for September to November. Precipitation indices show a tendency toward wetter conditions throughout the 20th century.

Perturbed Rankine vortices in surface quasi-geostrophic dynamics

An analytical dispersion relation is derived for linear perturbations to a Rankine vortex governed by surface quasi-geostrophic dynamics. Such a Rankine vortex is a circular region of uniform anomalous surface temperature evolving under quasi-geostrophic dynamics with uniform interior potential vorticity. The dispersion relation is analysed in detail and compared to the more familiar dispersion relation for a perturbed Rankine vortex governed by the Euler equations. The results are successfully verified against numerical simulations of the full equations. The dispersion relation is relevant to problems including wave propagation on surface temperature fronts and the stability of vortices in quasi-geostrophic turbulence.

Associations between stratospheric variability and tropospheric blocking

There is widely believed to be a link between stratospheric flow variability and stationary, persistent “blocking” weather systems, but the precise nature of this link has proved elusive. Using data from the ERA-40 Reanalysis and an atmospheric general circulation model (GCM) with a well-resolved stratosphere (HadGAM), it is shown that there are in fact several different highly significant associations, with blocking in different regions being related to different patterns of stratospheric variability. This is true in both hemispheres and in both data sets. The associations in HadGAM are shown to be very similar to those in ERA-40, although the model has a tendency to underestimate both European blocking and the wave number 2 stratospheric variability to which this is related. Although the focus is on stratospheric variability in general, several of the blocking links are seen to occur in association with the major stratospheric sudden warmings. In general, the direction of influence appears to be upward, as blocking anomalies are shown to modify the planetary stationary waves, leading to an upward propagation of wave activity into the stratosphere. However, significant correlations are also apparent with the zonal mean flow in the stratosphere leading the occurrence of blocking at high latitudes. Finally, the underestimation of blocking is an enduring problem in GCMs, and an example has recently been given in which improving the resolution of the stratosphere improved the representation of blocking. Here, however, another example is given, in which increasing the stratospheric resolution unfortunately does not lead to an improvement in blocking.

Galactic cosmic ray hazard in the unusual extended solar minimum between solar cycle 23 and 24

Galactic cosmic rays (GCRs) are extremely difficult to shield against and pose one of the most severe long-term hazards for human exploration of space. The recent solar minimum between solar cycles 23 and 24 shows a prolonged period of reduced solar activity and low interplanetary magnetic field strengths. As a result, the modulation of GCRs is very weak, and the fluxes of GCRs are near their highest levels in the last 25 years in the fall of 2009. Here we explore the dose rates of GCRs in the current prolonged solar minimum and make predictions for the Lunar Reconnaissance Orbiter (LRO) Cosmic Ray Telescope for the Effects of Radiation (CRaTER), which is now measuring GCRs in the lunar environment. Our results confirm the weak modulation of GCRs leading to the largest dose rates seen in the last 25 years over a prolonged period of little solar activity.

Where do winds drive the Antarctic Circumpolar Current?

The strength of the Antarctic Circumpolar Current (ACC) is believed to depend on the westerly wind stress blowing over the Southern Ocean, although the exact relationship between winds and circumpolar transport is yet to be determined. Here we show, based on theoretical arguments and a hierarchy of numerical modeling experiments, that the global pycnocline depth and the baroclinic ACC transport are set by an integral measure of the wind stress over the path of the ACC, taking into account its northward deflection. Our results assume that the mesoscale eddy diffusivity is independent of the mean flow; while the relationship between wind stress and ACC transport will be more complicated in an eddy-saturated regime, our conclusion that the ACC is driven by winds over the circumpolar streamlines is likely to be robust.

Excess open solar magnetic flux from satellite data: 1. Analysis of the third perihelion Ulysses pass

We use the third perihelion pass by the Ulysses spacecraft to illustrate and investigate the “flux excess” effect, whereby open solar flux estimates from spacecraft increase with increasing heliocentric distance. We analyze the potential effects of small-scale structure in the heliospheric field (giving fluctuations in the radial component on timescales smaller than 1 h) and kinematic time-of-flight effects of longitudinal structure in the solar wind flow. We show that the flux excess is explained by neither very small-scale structure (timescales < 1 h) nor by the kinematic “bunching effect” on spacecraft sampling. The observed flux excesses is, however, well explained by the kinematic effect of larger-scale (>1 day) solar wind speed variations on the frozen-in heliospheric field. We show that averaging over an interval T (that is long enough to eliminate structure originating in the heliosphere yet small enough to avoid cancelling opposite polarity radial field that originates from genuine sector structure in the coronal source field) is only an approximately valid way of allowing for these effects and does not adequately explain or account for differences between the streamer belt and the polar coronal holes.

Excess open solar magnetic flux from satellite data: 2. A survey of kinematic effects

We investigate the “flux excess” effect, whereby open solar flux estimates from spacecraft increase with increasing heliocentric distance. We analyze the kinematic effect on these open solar flux estimates of large-scale longitudinal structure in the solar wind flow, with particular emphasis on correcting estimates made using data from near-Earth satellites. We show that scatter, but no net bias, is introduced by the kinematic “bunching effect” on sampling and that this is true for both compression and rarefaction regions. The observed flux excesses, as a function of heliocentric distance, are shown to be consistent with open solar flux estimates from solar magnetograms made using the potential field source surface method and are well explained by the kinematic effect of solar wind speed variations on the frozen-in heliospheric field. Applying this kinematic correction to the Omni-2 interplanetary data set shows that the open solar flux at solar minimum fell from an annual mean of 3.82 × 1016 Wb in 1987 to close to half that value (1.98 × 1016 Wb) in 2007, making the fall in the minimum value over the last two solar cycles considerably faster than the rise inferred from geomagnetic activity observations over four solar cycles in the first half of the 20th century.

Space physics for graduate students: an activities-based approach

The geospace environment is controlled largely by events on the Sun, such as solar flares and coronal mass ejections, which generate significant geomagnetic and upper atmospheric disturbances. The study of this Sun-Earth system, which has become known as space weather, has both intrinsic scientific interest and practical applications. Adverse conditions in space can damage satellites and disrupt communications, navigation, and electric power grids, as well as endanger astronauts. The Center for Integrated Space Weather Modeling (CISM), a Science and Technology Center (STC) funded by the U.S. National Science Foundation (see http://www.bu.edu/cism/), is developing a suite of integrated physics-based computer models that describe the space environment from the Sun to the Earth for use in both research and operations [Hughes and Hudson, 2004, p. 1241]. To further this mission, advanced education and training programs sponsored by CISM encourage students to view space weather as a system that encompasses the Sun, the solar wind, the magnetosphere, and the ionosphere/thermosphere. This holds especially true for participants in the CISM space weather summer school [Simpson, 2004].

Conservation of open solar magnetic flux and the floor in the heliospheric magnetic field

The near-Earth heliospheric magnetic field intensity, |B|, exhibits a strong solar cycle variation, but returns to the same ``floor'' value each solar minimum. The current minimum, however, has seen |B| drop below previous minima, bringing in to question the existence of a floor, or at the very least requiring a re-assessment of its value. In this study we assume heliospheric flux consists of a constant open flux component and a time-varying contribution from CMEs. In this scenario, the true floor is |B| with zero CME contribution. Using observed CME rates over the solar cycle, we estimate the ``no-CME'' |B| floor at ~4.0 +/- 0.3 nT, lower than previous floor estimates and below |B| observed this solar minimum. We speculate that the drop in |B| observed this minimum may be due to a persistently lower CME rate than the previous minimum, though there are large uncertainties in the supporting observational data.

Suprathermal electron evolution in a Parker spiral magnetic field

Suprathermal electrons (>70 eV) form a small fraction of the total solar wind electron density but serve as valuable tracers of heliospheric magnetic field topology. Their usefulness as tracers of magnetic loops with both feet rooted on the Sun, however, most likely fades as the loops expand beyond some distance owing to scattering. As a first step toward quantifying that distance, we construct an observationally constrained model for the evolution of the suprathermal electron pitch-angle distributions on open field lines. We begin with a near-Sun isotropic distribution moving antisunward along a Parker spiral magnetic field while conserving magnetic moment, resulting in a field-aligned strahl within a few solar radii. Past this point, the distribution undergoes little evolution with heliocentric distance. We then add constant (with heliocentric distance, energy, and pitch angle) ad-hoc pitch-angle scattering. Close to the Sun, pitch-angle focusing still dominates, again resulting in a narrow strahl. Farther from the Sun, however, pitch-angle scattering dominates because focusing is effectively weakened by the increasing angle between the magnetic field direction and intensity gradient, a result of the spiral field. We determine the amount of scattering required to match Ulysses observations of strahl width in the fast solar wind, providing an important tool for inferring the large-scale properties and topologies of field lines in the interplanetary medium. Although the pitch-angle scattering term is independent of energy, time-of-flight effects in the spiral geometry result in an energy dependence of the strahl width that is in the observed sense although weaker in magnitude.

Ambient solar wind's effect on ICME transit times

Most empirical and numerical models of Interplanetary Coronal Mass Ejection (ICME) propagation use the initial CME velocity as their primary, if not only, observational input. These models generally predict a wide spread of 1 AU transit times for ICMEs with the same initial velocity. We use a 3D coupled MHD model of the corona and heliosphere to determine the ambient solar wind's effect on the propagation of ICMEs from 30 solar radii to 1 AU. We quantitatively characterize this deceleration by the velocity of the upstream ambient solar wind. The effects of varying solar wind parameters on the ICME transit time are quantified and can explain the observed spread in transit times for ICMEs of the same initial velocity. We develop an adjustment formula that can be used in conjunction with other models to reduce the spread in predicted transit times of Earth-directed ICMEs.

Metrics for solar wind prediction models: Comparison of empirical, hybrid, and physics-based schemes with 8 years of L1 observations

Space weather effects on technological systems originate with energy carried from the Sun to the terrestrial environment by the solar wind. In this study, we present results of modeling of solar corona-heliosphere processes to predict solar wind conditions at the L1 Lagrangian point upstream of Earth. In particular we calculate performance metrics for (1) empirical, (2) hybrid empirical/physics-based, and (3) full physics-based coupled corona-heliosphere models over an 8-year period (1995–2002). L1 measurements of the radial solar wind speed are the primary basis for validation of the coronal and heliosphere models studied, though other solar wind parameters are also considered. The models are from the Center for Integrated Space-Weather Modeling (CISM) which has developed a coupled model of the whole Sun-to-Earth system, from the solar photosphere to the terrestrial thermosphere. Simple point-by-point analysis techniques, such as mean-square-error and correlation coefficients, indicate that the empirical coronal-heliosphere model currently gives the best forecast of solar wind speed at 1 AU. A more detailed analysis shows that errors in the physics-based models are predominately the result of small timing offsets to solar wind structures and that the large-scale features of the solar wind are actually well modeled. We suggest that additional “tuning” of the coupling between the coronal and heliosphere models could lead to a significant improvement of their accuracy. Furthermore, we note that the physics-based models accurately capture dynamic effects at solar wind stream interaction regions, such as magnetic field compression, flow deflection, and density buildup, which the empirical scheme cannot.