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.     

Damped Oscillations of Coronal Loops

A mechanism of damped oscillations of a coronal loop is investigated. The loop is treated as a thin toroidal flux rope with two stationary photospheric footpoints, carrying both toroidal and poloidal currents. The forces and the flux-rope dynamics are described within the framework of ideal magnetohydrodynamics (MHD). The main features of the theory are the following: i) Oscillatory motions are determined by the Lorentz force that acts on curved current-carrying plasma structures and ii) damping is caused by drag that provides the momentum coupling between the flux rope and the ambient coronal plasma. The oscillation is restricted to the vertical plane of the flux rope. The initial equilibrium flux rope is set into oscillation by a pulse of upflow of the ambient plasma. The theory is applied to two events of oscillating loops observed by the Transition Region and Coronal Explorer (TRACE). It is shown that the Lorentz force and drag with a reasonable value of the coupling coefficient (c _d ) and without anomalous dissipation are able to accurately account for the observed damped oscillations. The analysis shows that the variations in the observed intensity can be explained by the minor radial expansion and contraction. For the two events, the values of the drag coefficient consistent with the observed damping times are in the range c _d ≈2 – 5, with specific values being dependent on parameters such as the loop density, ambient magnetic field, and the loop geometry. This range is consistent with a previous MHD simulation study and with values used to reproduce the observed trajectories of coronal mass ejections (CMEs).     

Large-scale cosmic homogeneity from a multifractal analysis of the PSCz catalogue

We investigate the behaviour of galaxy clustering on large scales using the PSC z catalogue. In particular, we ask whether there is any evidence of large-scale fractal behaviour in this catalogue. We find the correlation dimension in this survey varies with scale, consistent with other analyses. For example, our results on small and intermediate scales are consistent those obtained from the QDOT sample, but the larger PSC z sample allows us to extend the analysis out to much larger scales. We find firm evidence that the sample becomes homogeneous at large scales; the correlation dimension of the sample is D ₂=2.992±0.003 for r >30  h ⁻¹ Mpc. This provides strong evidence in favour of a universe that obeys the cosmological principle.     

The rotation of Janus and Epimetheus

Epimetheus, a small moon of Saturn, has a rotational libration (an oscillation about synchronous rotation) of 5.9°±1.2°, placing Epimetheus in the company of Earth’s Moon and Mars’ Phobos as the only natural satellites for which forced rotational libration has been detected. The forced libration is caused by the satellite’s slightly eccentric orbit and non-spherical shape.Detection of a moon’s forced libration allows us to probe its interior by comparing the measured amplitude to that predicted by a shape model assuming constant density. A discrepancy between the two would indicate internal density asymmetries. For Epimetheus, the uncertainties in the shape model are large enough to account for the measured libration amplitude. For Janus, on the other hand, although we cannot rule out synchronous rotation, a permanent offset of several degrees between Janus’ minimum moment of inertia (long axis) and the equilibrium sub-Saturn point may indicate that Janus does have modest internal density asymmetries.The rotation states of Janus and Epimetheus experience a perturbation every 4 years, as the two moons “swap” orbits. The sudden change in the orbital periods produces a free libration about synchronous rotation that is subsequently damped by internal friction. We calculate that this free libration is small in amplitude (<0.1°) and decays quickly (a few weeks, at most), and is thus below the current limits for detection using Cassini images.     

Assimilation of thermal emission spectrometer atmospheric data during the Mars Global Surveyor aerobraking period

The Thermal Emission Spectrometer aboard the Mars Global Surveyor spacecraft has produced an extensive atmospheric data set, beginning during aerobraking and continuing throughout the extended scientific mapping phase. Temperature profiles for the atmosphere below about 40 km, surface temperatures and total dust and water ice opacities, can be retrieved from infrared spectra in nadir viewing mode. This paper describes assimilation of nadir retrievals from the spacecraft aerobraking period, L_S=190°–260°, northern hemisphere autumn to winter, into a Mars general circulation model. The assimilation scheme is able to combine information from temperature and dust optical depth retrievals, making use of a model forecast containing information from the assimilation of earlier observations, to obtain a global, time-dependent analysis. Given sufficient temperature retrievals, the assimilation procedure indicates errors in the a priori dust distribution assumptions even when lacking dust observations; in this case there are relatively cold regions above the poles compared to a model which assumes a horizontally-uniform dust distribution. One major reason for using assimilation techniques is in order to investigate the transient wave behavior on Mars. Whilst the data from the 2-h spacecraft mapping orbit phase is much more suitable for assimilation, even the longer (45–24 h) period aerobraking orbit data contain useful information about the three-dimensional synoptic-scale martian circulation which the assimilation procedure can reconstruct in a consistent way. Assimilations from the period of the Noachis regional dust storm demonstrate that the combined assimilation of temperature and dust retrievals has a beneficial impact on the atmospheric analysis.     

The topography of Iapetus' leading side

We have used Cassini stereo images to study the topography of Iapetus' leading side. A terrain model derived at resolutions of 4-8 km reveals that Iapetus has substantial topography with heights in the range of −10 km to +13 km, much more than observed on the other middle-sized satellites of Saturn so far. Most of the topography is older than 4 Ga [Neukum, G., Wagner, R., Denk, T., Porco, C.C., 2005. Lunar Planet. Sci. XXXVI. Abstract 2034] which implies that Iapetus must have had a thick lithosphere early in its history to support this topography. Models of lithospheric deflection by topographic loads provide an estimate of the required elastic thickness in the range of 50-100 km. Iapetus' prominent equatorial ridge [Porco, C.C., and 34 colleagues, 2005. Science 307, 1237-1242] reaches widths of 70 km and heights of up to 13 km from their base within the modeled area. The morphology of the ridge suggests an endogenous origin rather than a formation by collisional accretion of a ring remnant [Ip, W.-H., 2006. Geophys. Res. Lett. 33, doi:10.1029/2005GL025386. L16203]. The transition from simple to complex central peak craters on Iapetus occurs at diameters of 11±3 km. The central peaks have pronounced conical shapes with flanking slopes of typically 11° and heights that can rise above the surrounding plains. Crater depths seem to be systematically lower on Iapetus than on similarly sized Rhea, which if true, may be related to more pronounced crater-wall slumping (which widens the craters) on Iapetus than on Rhea. There are seven large impact basins with complex morphologies including central peak massifs and terraced walls, the largest one reaches 800 km in diameter and has rim topography of up to 10 km. Generally, no rings are observed with the basins consistent with a thick lithosphere but still thin enough to allow for viscous relaxation of the basin floors, which is inferred from crater depth-to-diameter measurements. In particular, a 400-km basin shows up-domed floor topography which is suggestive of viscous relaxation. A model of complex crater formation with a viscoplastic (Bingham) rheology [Melosh, H.J., 1989. Impact Cratering. Oxford Univ. Press, New York] of the impact-shocked icy material provides an estimate of the effective cohesion/viscosity at . The local distribution of bright and dark material on the surface of Iapetus is largely controlled by topography and consistent with the dark material being a sublimation lag deposit originating from a bright icy substrate mixed with the dark components, but frost deposits are possible as well.     

Global solar dynamo models: Simulations and predictions

Flux-transport type solar dynamos have achieved considerable success in correctly simulating many solar cycle features, and are now being used for prediction of solar cycle timing and amplitude. We first define flux-transport dynamos and demonstrate how they work. The essential added ingredient in this class of models is meridional circulation, which governs the dynamo period and also plays a crucial role in determining the Sun’s memory about its past magnetic fields. We show that flux-transport dynamo models can explain many key features of solar cycles. Then we show that a predictive tool can be built from this class of dynamo that can be used to predict mean solar cycle features by assimilating magnetic field data from previous cycles.     

Search for extinct aluminum‐26 and titanium‐44 in nanodiamonds from the Allende CV3 and Murchison CM2 meteorites

Abstract– We report Mg‐Al and Ca‐Ti isotopic data for meteoritic nanodiamonds separated from the Allende CV3 and Murchison CM2 meteorites. The goal of this study was to search for excesses in ²⁶Mg and ⁴⁴Ca, which can be attributed to the in situ decay of radioactive and now extinct ²⁶Al and ⁴⁴Ti, respectively. Previous work on presolar SiC and graphite had shown that ²⁶Al/²⁷Al and ⁴⁴Ti/⁴⁸Ti ratios in presolar grains can be used to discriminate between different types of stellar sources. Aluminum and Ti concentrations are low in the meteoritic nanodiamonds of this study. Murchison nanodiamonds have higher Al and Ti concentrations than the Allende nanodiamonds. This can be attributed to contamination and the presence of presolar SiC in the Murchison nanodiamond samples. ²⁶Mg/²⁴Mg and ⁴⁴Ca/⁴⁰Ca ratios are close to normal in Allende nanodiamonds with upper limits on the initial ²⁶Al/²⁷Al and ⁴⁴Ti/⁴⁸Ti ratios of approximately 1 × 10^?3. These ratios are factors of 10–1000 and, respectively, 1–1000 lower than those of presolar SiC and graphite grains from supernovae. The ²⁶Al/²⁷Al and ⁴⁴Ti/⁴⁸Ti data for nanodiamonds are compatible with an asymptotic giant branch star or solar system origin, but not with a supernova origin of a major fraction of meteoritic nanodiamonds. The latter possibility cannot be excluded, though, as the diamond separates may contain significant amounts of contaminating Al and Ti, which would lower the inferred ²⁶Al/²⁷Al and ⁴⁴Ti/⁴⁸Ti ratios considerably.     

The Evolution of Sunspot Magnetic Fields Associated with a Solar Flare

Solar flares occur due to the sudden release of energy stored in active-region magnetic fields. To date, the precursors to flaring are still not fully understood, although there is evidence that flaring is related to changes in the topology or complexity of an active-region’s magnetic field. Here, the evolution of the magnetic field in active region NOAA 10953 was examined using Hinode/SOT-SP data over a period of 12 hours leading up to and after a GOES B1.0 flare. A number of magnetic-field properties and low-order aspects of magnetic-field topology were extracted from two flux regions that exhibited increased Ca ii H emission during the flare. Pre-flare increases in vertical field strength, vertical current density, and inclination angle of ≈ 8° toward the vertical were observed in flux elements surrounding the primary sunspot. The vertical field strength and current density subsequently decreased in the post-flare state, with the inclination becoming more horizontal by ≈ 7°. This behavior of the field vector may provide a physical basis for future flare-forecasting efforts.