Stress field control of eruption dynamics at a rift volcano: Nyamuragira, D.R.Congo

Using the record of 30 flank eruptions over the last 110 years at Nyamuragira, we have tested the relationship between the eruption dynamics and the local stress field. There are two groups of eruptions based on their duration (< 80days >) that are also clustered in space and time. We find that the eruptions fed by dykes parallel to the East African Rift Valley have longer durations (and larger volumes) than those eruptions fed by dykes with other orientations. This is compatible with a model for compressible magma transported through an elastic-walled dyke in a differential stress field from an over-pressured reservoir (Woods et al., 2006). The observed pattern of eruptive fissures is consistent with a local stress field modified by a northwest-trending, right lateral slip fault that is part of the northern transfer zone of the Kivu Basin rift segment. We have also re-tested with new data the stochastic eruption models for Nyamuragira of Burt et al. (1994). The time-predictable, pressure-threshold model remains the best fit and is consistent with the typically observed declining rate of sulphur dioxide emission during the first few days of eruption with lava emission from a depressurising, closed, crustal reservoir. The 2.4-fold increase in long-term eruption rate that occurred after 1977 is confirmed in the new analysis. Since that change, the record has been dominated by short-duration eruptions fed by dykes perpendicular to the Rift. We suggest that the intrusion of a major dyke during the 1977 volcano-tectonic event at neighbouring Nyiragongo volcano inhibited subsequent dyke formation on the southern flanks of Nyamuragira and this may also have resulted in more dykes reaching the surface elsewhere. Thus that sudden change in output was a result of a changed stress field that forced more of the deep magma supply to the surface. Another volcano-tectonic event in 2002 may also have changed the magma output rate at Nyamuragira.

Are superhumps good measures of the mass ratio for AM CVn systems?

We extend recent work that included the effect of pressure forces to derive the precession rate of eccentric accretion discs in cataclysmic variables to the case of double degenerate systems. We find that the logical scaling of the pressure force in such systems results in predictions of unrealistically high primary masses. Using the prototype AM CVn as a calibrator for the magnitude of the effect, we find that there is no scaling that applies consistently to all the systems in the class. We discuss the reasons for the lack of a superhump period to mass ratio relationship analogous to that known for SU UMa systems and suggest that this is because these secondaries do not have a single valued mass-radius relationship. We highlight the unreliability of mass-ratios derived by applying the SU UMa expression to the AM CVn binaries.

Superhumps: confronting theory with observation

We review the theory and observations related to the "superhump" precession of eccentric accretion discs in close binary systems. We agree with earlier work, although for different reasons, that the discrepancy between observation and dynamical theory implies that the effect of pressure in the disc cannot be neglected. We extend earlier work that investigates this effect to include the correct expression for the radius at which resonant orbits occur. Using analytic expressions for the accretion disc structure, we derive a relationship between the period excess and mass ratio with the pressure effects included. This is compared to the observed data, recently derived results for detailed integration of the disc equations and the equivalent empirically derived relations and used to predict values for the mass ratio based on measured values of the period excess for 88 systems.

Superhumps, magnetic fields and the mass ratio in AM Canum Venaticorum

We show that the observed K velocities and periodicities of AM CVn can be reconciled given a mass ratio q~0.22 and a secondary star with a modest magnetic field of surface strength B~1 T. We see that the new mass ratio implies that the secondary is most likely semidegenerate. The effect of the field on the accretion disc structure is examined. The theory of precessing discs and resonant orbits is generalized to encompass higher order resonances than 3:2 and shown to retain consistency with the new mass ratio.

The effect of the magnetic field of the secondary star in dwarf novae

We investigate the effect of a secondary star magnetic field on the accretion disc dynamics of dwarf novae. Simulations have been carried out with a particle code and a dipolar magnetic field structure. The magnetic field acts to remove angular momentum from the disc material, increasing the inward mass flow. This makes the accretion disc more centrally condensed, causing a reduction in the recurrence time for dwarf nova outbursts. We have produced Doppler tomograms and light curves which may be compared with observations. These tomograms are significantly different from those produced in the absence of a magnetic field on the secondary. We derive an upper limit to the magnetic moment of the secondary star in UGem of mu_2<2x10^32 A m^2. The magnetic truncation of the accretion disc produces resonance phenomena similar to those seen in the superoutbursts of SUUMa systems. While these have not been observed for systems like UGem, observations of the SUUMa systems provide us with a useful diagnostic of the disc-field interaction. We are able to place an upper limit on the magnetic moment of the secondary in ZCha of mu_2<1x10^30 A m^2.

Spin equilibrium and outflows in young stellar objects

We present N-body simulations of accretion discs about young stellar objects (YSOs). The simulation includes the presence of a magnetic loop structure on the central star which interacts with the particles by means of a magnetic drag force. We find that an equilibrium spin rate is achieved when the corotation radius coincides with the edge of the loop. This spin rate is consistent with observed values for TTauri stars, being an order of magnitude less than the breakup value. The material ejected from the system by the rotating loop has properties consistent with the observed molecular outflows, given the presence of a suitable containing cavity.

Rapid Keck Spectroscopy of Cataclysmic Variables

We present an analysis of Rapid Keck Spectroscopy of the CVs AM Her (polar) and SS Cyg (dwarf nova). We decompose the spectra into constant and variable components and identify different types of variability in AM Her with different characteristic timescales. The variable flickering component of the accretion disc flux and the observational characteristics of a small flare in SS Cyg are isolated.

Climate at high-obliquity

The question of climate at high obliquity is raised in the context of both exoplanet studies (e.g. habitability) and paleoclimates studies (evidence for low-latitude glaciation during the Neoproterozoic and the ”Snowball Earth” hypothesis). States of high obliquity, φ, are distinctive in that, for φ ≥54◦, the poles receive more solar radiation in the annual mean than the Equator, opposite to the present day situation. In addition, the seasonal cycle of insolation is extreme, with the poles alternatively “facing” the sun and sheltering in the dark for months. The novelty of our approach is to consider the role of a dynamical ocean in controlling the surface climate at high obliquity, which in turn requires understanding of the surface winds patterns when temperature gradients are reversed. To address these questions, a coupled ocean-atmosphere-sea ice GCM configured on an aquaplanet is employed. Except for the absence of topography and modified obliquity, the set-up is Earth-like. Two large obliquities φ, 54◦ and 90◦, are compared to today’s Earth value, φ=23.5◦. Three key results emerge at high obliquity: 1) despite reversed temper- ature gradients, mid-latitudes surface winds are westerly and trade winds exist at the equator (as for φ=23.5◦) although the westerlies are confined to the summer hemisphere, 2) a habitable planet is possible with mid-latitude temperatures in the range 300-280 K and 3) a stable climate state with an ice cap limited to the equatorial region is unlikely. We clarify the dynamics behind these features (notably by an analysis of the potential vorticity structure and conditions for baroclinic instability of the atmosphere). Interestingly, we find that the absence of a stable partially glaciated state is critically linked to the absence of ocean heat transport during winter, a feature ultimately traced back to the high seasonality of baroclinic instability conditions in the atmosphere.

A regime-dependent parametrization of subgrid-scale cloud water content variability

The subgrid-scale spatial variability in cloud water content can be described by a parameter f called the fractional standard deviation. This is equal to the standard deviation of the cloud water content divided by the mean. This parameter is an input to schemes that calculate the impact of subgrid-scale cloud inhomogeneity on gridbox-mean radiative fluxes and microphysical process rates. A new regime-dependent parametrization of the spatial variability of cloud water content is derived from CloudSat observations of ice clouds. In addition to the dependencies on horizontal and vertical resolution and cloud fraction included in previous parametrizations, the new parametrization includes an explicit dependence on cloud type. The new parametrization is then implemented in the Global Atmosphere 6 (GA6) configuration of the Met Office Unified Model and used to model the effects of subgrid variability of both ice and liquid water content on radiative fluxes and autoconversion and accretion rates in three 20-year atmosphere-only climate simulations. These simulations show the impact of the new regime-dependent parametrization on diagnostic radiation calculations, interactive radiation calculations and both interactive radiation calculations and in a new warm microphysics scheme. The control simulation uses a globally constant f value of 0.75 to model the effect of cloud water content variability on radiative fluxes. The use of the new regime-dependent parametrization in the model results in a global mean which is higher than the control's fixed value and a global distribution of f which is closer to CloudSat observations. When the new regime-dependent parametrization is used in radiative transfer calculations only, the magnitudes of short-wave and long-wave top of atmosphere cloud radiative forcing are reduced, increasing the existing global mean biases in the control. When also applied in a new warm microphysics scheme, the short-wave global mean bias is reduced.

Spatial variability of liquid cloud and rain: observations and microphysical effects

Spatial variability of liquid cloud water content and rainwater content is analysed from three different observational platforms: in situ measurements from research aircraft, land-based remote sensing techniques using radar and lidar, and spaceborne remote sensing from CloudSat. The variance is found to increase with spatial scale, but also depends strongly on the cloud or rain fraction regime, with overcast regions containing less variability than broken cloud fields. This variability is shown to lead to large biases, up to a factor of 4, in both the autoconversion and accretion rates estimated at a model grid scale of ≈40 km by a typical microphysical parametrization using in-cloud mean values. A parametrization for the subgrid variability of liquid cloud and rainwater content is developed, based on the observations, which varies with both the grid scale and cloud or rain fraction, and is applicable for all model grid scales. It is then shown that if this parametrization of the variability is analytically incorporated into the autoconversion and accretion rate calculations, the bias is significantly reduced.