Thin-Shell Magnetohydrodynamic Equations for the Solar Tachocline

We present a new system of equations designed to study global-scale dynamics in the stably-stratified portion of the solar tachocline. This system is derived from the 3D equations of magnetohydrodynamics in a rotating spherical shell under the assumption that the shell is thin and stably-stratified (subadiabatic). The resulting thin-shell model can be regarded as a magnetic generalization of the hydrostatic primitive equations often used in meteorology. It is simpler in form than the more general anelastic or Boussinesq equations, making it more amenable to analysis and interpretation and more computationally efficient. However, the thin-shell system is still three-dimensional and as such represents an important extension to previous 2D and shallow-water approaches. In this paper we derive the governing equations for our thin-shell model and discuss its underlying assumptions, its context relative to other models, and its application to the solar tachocline. We also demonstrate that the dissipationless thin-shell system conserves energy, angular momentum and magnetic helicity.     

An iterative scheme for use in MT 2D forward modelling — to speed up production of the Jacobian and to improve accuracy of the inverse matrix

An important stage in two-dimensional magnetotelluric modelling is the calculation of the Earth's response functions for an assumed conductivity model and the calculation of the associated Jacobian relating those response functions to the model parameters. The efficiency of the calculation of the Jacobian will affect the efficiency of the inversion modelling. Rodi (1976) produced all the Jacobian elements by inverting a single matrix and using an approximate first-order algorithm. Since only one inverse matrix required calculation the procedure speeded up the inversion. An iterative scheme to improve the approximation to the Jacobian information is presented in this paper. While this scheme takes a little longer than Rodi's algorithm, it enables a more accurate determination of the Jacobian information. It is found that the Jacobian elements can be produced in 10% of the time required to calculate an inverse matrix or to calculate a 2D starting model. A modification of the algorithm can further be used to improve the accuracy of the original inverse matrix calculated in a 2D finite difference program and hence the solution this program produces. The convergence of the iteration scheme is found to be related both to the originally calculated inverse matrix and to the change in the newly formed matrix arising from perturbation of the model parameter. A ridge regression inverse algorithm is used in conjunction with the iterative scheme for forward modelling described in this paper to produce a 2D conductivity section from field data.     

Photometric Inference of the Corotation Radii of Ten Barred Galaxies

We have used photometric images of ten barred galaxies in the B and I bands to infer the geometrical and dynamical parameters of their bars: Their length, (R _bar), their strengths, (S _b), and their corotation radii, (R _CR).These parameters have been obtained studying azimuthal profiles from the B and I band images of the galaxies. We find that R _CR is in all cases slightly larger than R _bar, with mean values close to 1.2 R _bar. We have also found a dependence of the ratio R _CR/R _bar on S _bar. This revised version was published online in July 2006 with corrections to the Cover Date.     

Persistent 22-year cycle in sunspot activity: Evidence for a relic solar magnetic field

We use the recently presented group sunspot number series to show that a persistent 22-year cyclicity exists in sunspot activity throughout the entire period of about 400 years of direct sunspot observations. The amplitude of this cyclicity is about 10% of the present sunspot activity level. A 22-year cyclicity in sunspot activity is naturally produced by the 22-year magnetic polarity cycle in the presence of a relic dipole magnetic field. Accordingly, a persistent 22-year cyclicity in sunspot activity gives an evidence for the existence of such a relic magnetic field in the Sun. The stable phase and the roughly constant amplitude of this cyclicity during times of very different sunspot activity level strongly support this interpretation.     

Investigation of Barred Galaxies. VII. Comparative Statistics of SB and SA Galaxies. Near and Mid-IR Regions

The near and mid-IR properties of barred and unbarred spiral galaxies are discussed on the basis of complete samples that were compiled earlier. The two types of spirals are shown not to differ from one another in emission power in the near and mid-IR ranges. Multiple regression and principal component analysis have been applied to investigate near and mid-IR properties of SB and SA galaxies, particularly their relation to X-ray and radio continuum emissions. There are definite differences between SB and SA spirals in the near and mid-IR. In the case of SB galaxies, the compactness of 10 m emission is closely related to the J - H color index, and the redder J - H color corresponds to relatively more extended emission at 10 m. It is assumed that these are caused by the bar itself, which stimulates enhanced star formation in a barred galaxy with respect to unbarred spiral.     

Structure through colour: a pixel approach towards understanding galaxies

We present a study of pixel colour–magnitude diagrams (pCMDs) for a sample of 69 nearby galaxies chosen to span a wide range of Hubble types. Our goal is to determine how useful a pixel approach is for studying galaxies according to their stellar light distributions and content. The galaxy images were analysed on a pixel-by-pixel basis to reveal the structure of the individual pCMDs. We find that the average surface brightness (or projected mass density) in each pixel varies according to galaxy type. Early-type galaxies exhibit a clear 'prime sequence' and some pCMDs of face-on spirals reveal 'inverse-L' structures. We find that the colour dispersion at a given magnitude is found to be approximately constant in early-type galaxies but this quantity varies in the mid and late types. We investigate individual galaxies and find that the pCMDs can be used to pick out morphological features. We discuss the discovery of 'Red Hooks' in the pCMDs of six early-type galaxies and two spirals and postulate their origins. We develop quantitative methods to characterize the pCMDs, including measures of the blue-to-red light ratio and colour distributions of each galaxy and we organize these by morphological type. We compare the colours of the pixels in each galaxy with the stellar population models of Bruzual & Charlot to calculate star formation histories for each galaxy type and compare these to the stellar mass within each pixel. Maps of pixel stellar mass and mass-to-light ratio are compared to galaxy images. We apply the pCMD technique to three galaxies in the Hubble Ultra Deep Field to test the usefulness of the analysis at high redshift. We propose that these results can be used as part of a new system of automated classification of galaxies that can be applied at high redshift.     

Identification of some of the unidentified EIC-1 sources

We have identified 20 of the 75 unidentified Equatorial Infrared Catalogue Number-1 sources with stars in the BD Catalogue. Two of these are also found in the HD Catalogue.     

Excitation of whistler waves by reflected auroral electrons

Excitation of electron cyclotron waves and whistlers by reflected auroral electrons which possess a loss-cone distribution is investigated. Based on a given magnetic field and density model, the instability problem is studied over a broad region along the auroral field lines. This region covers altitudes ranging from one quarter of an Earth radius to five Earth radii. It is found that the growth rate is significant only in the region of low altitude, say below the source region of the auroral kilometric radiation. In the high altitude region the instability is insignificant either because of low refractive indices or because of small loss cone angles.