The existence of a stably stratified layer underneath the core-mantle boundary (CMB) has been recently revived by corroborating evidences coming from seismic studies, mineral physics and thermal evolution models. Such a layer could find its physical origination either in compositional stratification due to the accumulation of light elements at the top or the core or in thermal stratification due to the heat flux becoming locally sub-adiabatic. The exact properties of this stably-stratified layer, namely its size $\mathcal{H}_S$ and the degree of its stratification characterised by the Brunt-V\"ais\"al\"a frequency $N$, are however uncertain and highly debated. A stable layer underneath the CMB can have crucial dynamical impacts on the geodynamo. Because of the inhibition of the convective motions, a stable layer is expected to primarily act as a low-pass filter on the magnetic field, smoothing out the rapidly-varying and small-scale features by skin effect. To investigate this effect more systematically, we compute 70 global geodynamo models varying the size of the stably-stratified layer from 0 to 300~km and its amplitude from $N/\Omega = 0$ to $N/\Omega \simeq 50$, $\Omega$ being the rotation rate. We show that the penetration of the convective flow in the stably-stratified layer is controlled by the typical size of the convective eddies and by the local variations of the ratio $N/\Omega$. Using quantitative measures of the degree of morphological semblance between the magnetic field obtained in numerical models and the geomagnetic field at the CMB, we establish an upper bound for the stable layer thickness $\mathcal{H}_s < (N/\Omega)^{-1} d_c$, $d_c$ being the horizontal size of the convective flow at the base of the stable layer. This defines a strong geomagnetic constraint on the properties of a stably-stratified layer beneath the CMB.

We present in this work the development of a solar data assimilation method based on an axisymmetric mean field dynamo model and magnetic surface data, our mid-term goal is to predict the solar quasi cyclic activity. Here we focus on the ability of our algorithm to constrain the deep meridional circulation of the Sun based on solar magnetic observations. To that end, we develop a variational data assimilation technique. Within a given assimilation window, the assimilation procedure minimizes the differences between data and the forecast from the model, by finding an optimal meridional circulation in the convection zone, and an optimal initial magnetic field, via a quasi-Newton algorithm. We demonstrate the capability of the technique to estimate the meridional flow by a closed-loop experiment involving 40 years of synthetic, solar-like data. By assimilating the synthetic magnetic proxies annually, we are able to reconstruct a (stochastic) time-varying meridional circulation which is also slightly equatorially asymmetric. We show that the method is robust in estimating a flow whose level of fluctuation can reach 30% about the average, and that the horizon of predictive capability of the method is of the order of 1 cycle length.

We show how magnetic observations of the Sun can be used in conjunction with an axisymmetric flux-transport solar dynamo model in order to estimate the large-scale meridional circulation throughout the convection zone. Our innovative approach rests on variational data assimilation, whereby the distance between predictions and observations (measured by an objective function) is iteratively minimized by means of an optimization algorithm seeking the meridional flow which best accounts for the data. The minimization is performed using a quasi-Newton technique, which requires the knowledge of the sensitivity of the objective function to the meridional flow. That sensitivity is efficiently computed via the integration of the adjoint flux-transport dynamo model. Closed-loop (also known as twin) experiments using synthetic data demonstrate the validity and accuracy of this technique, for a variety of meridional flow configurations, ranging from unicellular and equatorially symmetric to multicellular and equatorially asymmetric. In this well-controlled synthetic context, we perform a systematic study of the behavior of our variational approach under different observational configurations, by varying their spatial density, temporal density, noise level, as well as the width of the assimilation window. We find that the method is remarkably robust, leading in most cases to a recovery of the true meridional flow to within better than 1%. These encouraging results are a first step towards using this technique to i) better constrain the physical processes occurring inside the Sun and ii) better predict solar activity on decadal time scales.

Space weather is a matter of practical importance in our modern society. Predictions of forecoming solar cycles mean amplitude and duration are currently being made based on flux-transport numerical models of the solar dynamo. Interested in the forecast horizon of such studies, we quantify the predictability window of a representative, advection-dominated, flux-transport dynamo model by investigating its sensitivity to initial conditions and control parameters through a perturbation analysis. We measure the rate associated with the exponential growth of an initial perturbation of the model trajectory, which yields a characteristic time scale known as the e-folding time $\tau_e$. The e-folding time is shown to decrease with the strength of the $\alpha$-effect, and to increase with the magnitude of the imposed meridional circulation. Comparing the e-folding time with the solar cycle periodicity, we obtain an average estimate for $\tau_e$ equal to 2.76 solar cycle durations. From a practical point of view, the perturbations analysed in this work can be interpreted as uncertainties affecting either the observations or the physical model itself. After reviewing these, we discuss their implications for solar cycle prediction.

Many photometric time-domain surveys are driven by specific goals, such as searches for supernovae or transiting exoplanets, which set the cadence with which fields are re-imaged. In the case of the Palomar Transient Factory (PTF), several sub-surveys are conducted in parallel, leading to non-uniform sampling over its $\sim$$20,000 \mathrm{deg}^2$ footprint. While the median $7.26 \mathrm{deg}^2$ PTF field has been imaged $\sim$40 times in \textit{R}-band, $\sim$$2300 \mathrm{deg}^2$ have been observed $>$100 times. We use PTF data to study the trade-off between searching for microlensing events in a survey whose footprint is much larger than that of typical microlensing searches, but with far-from-optimal time sampling. To examine the probability that microlensing events can be recovered in these data, we test statistics used on uniformly sampled data to identify variables and transients. We find that the von Neumann ratio performs best for identifying simulated microlensing events in our data. We develop a selection method using this statistic and apply it to data from fields with $>$10 $R$-band observations, $1.1\times10^9$ light curves, uncovering three candidate microlensing events. We lack simultaneous, multi-color photometry to confirm these as microlensing events. However, their number is consistent with predictions for the event rate in the PTF footprint over the survey's three years of operations, as estimated from near-field microlensing models. This work can help constrain all-sky event rate predictions and tests microlensing signal recovery in large data sets, which will be useful to future time-domain surveys, such as that planned with the Large Synoptic Survey Telescope.

Dynamo models relying on the Babcock-Leighton mechanism are successful in reproducing most of the solar magnetic field dynamical characteristics. However, considering that such models operate only above a lower magnetic field threshold, they do not provide an appropriate magnetic field regeneration process characterizing a self-sustainable dynamo. In this work we consider the existence of an additional \alpha-effect to the Babcock-Leighton scenario in a mean-field axisymmetric kinematic numerical model. Both poloidal field regeneration mechanisms are treated with two different strength-limiting factors. Apart from the solar anti-symmetric parity behavior, the main solar features are reproduced: cyclic polarity reversals, mid-latitudinal equatorward migration of strong toroidal field, poleward migration of polar surface radial fields, and the quadrature phase shift between both. Long-term variability of the solutions exhibits lengthy periods of minimum activity followed by posterior recovery, akin to the observed Maunder Minimum. Based on the analysis of the residual activity during periods of minimum activity, we suggest that these are caused by a predominance of the \alpha-effect over the Babcock-Leighton mechanism in regenerating the poloidal field.

We present tomography of the circum-galactic metal distribution at redshift 1.7 to 4.5 derived from echellete spectroscopy of binary quasars. We find CIV systems at similar redshifts in paired sightlines more often than expected for sightline-independent redshifts. As the separation of the sightlines increases from 36 kpc to 907 kpc, the amplitude of this clustering decreases. At the largest separations, the CIV systems cluster similar to Lyman-break galaxies (Adelberger et al. 2005a). The CIV systems are significantly less correlated than these galaxies, however, at separations less than R_1 ~ 0.42 +/- 0.15 h-1 comoving Mpc. Measured in real space, i.e., transverse to the sightlines, this length scale is significantly smaller than the break scale estimated from the line-of-sight correlation function in redshift space (Scannapieco et al. 2006a). Using a simple model, we interpret the new real-space measurement as an indication of the typical physical size of enriched regions. We adopt this size for enriched regions and fit the redshift-space distortion in the line-of-sight correlation function. The fitted velocity kick is consistent with the peculiar velocity of galaxies as determined by the underlying mass distribution and places an upper limit on the outflow (or inflow) speed of metals. The implied time scale for dispersing metals is larger than the typical stellar ages of Lyman-break galaxies (Shapley et al. 2001), and we argue that enrichment by galaxies at z > 4.3 played a greater role in dispersing metals. To further constrain the growth of enriched regions, we discuss empirical constraints on the evolution of the CIV correlation function with cosmic time. This study demonstrates the potential of tomography for measuring the metal enrichment history of the circum-galactic medium.

Azimuthal magnetorotational instability is a mechanism that generates nonaxisymmetric field pattern. Nonlinear simulations in an infinite Taylor-Couette system with current-free external field show, that not only the linearly unstable mode m=1 appears, but also an inverse cascade transporting energy into the axisymmetric field is possible. By varying the Reynolds number of the flow and the Hartmann number for the magnetic field, we find that the ratio between axisymmetric (m=0) and dominating nonaxisymmetric mode (m=1) can be nearly free chosen. On the surface of the outer cylinder this mode distribution appears similarly, but with weaker axisymmetric fields. We do not find significant differences in the case that a constant current within the flow is added.

We present numerical simulations of the flow in a rapidly rotating cylindrical annulus. We show that at the rotation rates relevant to the magneto-rotational instability, the flow is strongly constrained by the Taylor-Proudman theorem. As a result, it is controlled almost entirely by the end-plates. We then consider two possible options for minimizing these end-effects, namely (i) simply taking a very long cylinder, and (ii) splitting the end-plates into a series of differentially rotating rings. Regarding option (i), we show that the cylinder would have to be hundreds of times as long as it is wide before end-effects become unimportant in the interior. Since this is clearly not feasible, we turn to option (ii), and show that in order to obtain a smooth angular velocity profile, the end-plates would have to be split into around ten rings. If the end-plates are split into fewer rings, perhaps 3-5, the angular velocity profile will not be smooth, but will instead consist of a series of Stewartson layers at the boundaries from one ring to the next. We suggest therefore that the instabilities one obtains in this system will be the familiar Kelvin-Helmholtz instabilities of these Stewartson layers, rather than the magneto-rotational instability. At best, one might hope to obtain the MRI superimposed on these Kelvin-Helmholtz modes. Any subsequent interpretation of results is thus likely to be quite complicated.