We propose two ensembles of geomagnetic field models spanning the last three millennia: COV-ARCH is calculated using all available archeological artefacts and volcanic lava flows; lake and marine ...sediment records are added to this data set to build COV-LAKE. Given the sparse distribution of archeomagnetic observations and their associated large uncertainties, the recovery of magnetic field changes from such data is an ill-posed inverse problem that requires assuming some prior information. This is usually performed by imposing arbitrary regularizations in space and time. Instead, we construct the prior knowledge entering the objective function to be minimized from spatial and temporal statistics of the geomagnetic field, as available from satellites, ground-based observatories and paleomagnetic measurements, and validated by numerical geodynamo simulations. Our approach relies on the projection of model coefficients onto temporal cross-covariance functions. We show with synthetic experiments that the dispersion within the ensemble of solutions provides a coherent measure of the model uncertainties. Gauss coefficients inverted from geophysical records compare satisfactorily with those deduced from the independent database built upon historical and observatory records. A posteriori model errors are reduced when incorporating sediment records; they logically increase towards decreasing length-scales, indicating that a partial information is available up to a spherical degree 4-5. Such models and their associated uncertainties are suited to be used as observations in geomagnetic data assimilation studies. Our results advocate for an approximately constant dipole decay since ≈1700 AD, preceded by an era (going back to 1000 AD) where the dipole trend is weak, possibly slightly positive. We observe in both hemispheres, at both low- and high-latitude, persistent patches over the past 3000 yr. We also confirm a westward drift of flux patches at the core-mantle boundary at a speed of about 0.20 to 0.25° yr-1. Despite the sparse data distribution in the southern hemisphere, the South Atlantic Anomaly appears in both ensembles of models around 1800 AD. A similar low-intensity event seems to have appeared below the Indian Ocean over 600-1400 AD. Both global models are in general good agreement with regional master curves, though filtering out some of the centennial oscillations.
Animals require the ability to ignore sensory stimuli that have no consequence yet respond to the same stimuli when they become useful. However, the brain circuits that govern this flexibility in ...sensory processing are not well understood. Here we show in mouse primary auditory cortex (A1) that daily passive sound exposure causes a long-lasting reduction in representations of the experienced sound by layer 2/3 pyramidal cells. This habituation arises locally in A1 and involves an enhancement in inhibition and selective upregulation in the activity of somatostatin-expressing inhibitory neurons (SOM cells). Furthermore, when mice engage in sound-guided behavior, pyramidal cell excitatory responses to habituated sounds are enhanced, whereas SOM cell responses are diminished. Together, our results demonstrate the bidirectional modulation of A1 sensory representations and suggest that SOM cells gate cortical information flow based on the behavioral relevance of the stimulus.
•Passive sound experience causes habituation of sensory representations in A1•Habituation involves an increase in inhibition of layer 2/3 pyramidal cells•Habituation reflects the selective upregulation of SOM interneuron activity•Sound-guided behavior decreases SOM cell activity and rapidly reverses habituation
Kato et al. demonstrate that daily sound experience upregulates SOM interneuron activity in A1 and causes “habituation” of sound representations. Sound-guided behavior reverses these effects, indicating that sensory representations are bidirectionally modified based on the behavioral relevance of sensory stimuli.
ABSTRACT
The process of heating and reionization of the Universe at high redshifts links small-scale structures/galaxy formation and large-scale intergalactic medium (IGM) properties. Even if the ...first one is difficult to observe, an observation window is opening on the second one, with the promising development of current and future radio telescopes. They will permit to observe the 21 cm brightness temperature global signal and fluctuations. The need for large-scale simulations is therefore strong to understand the properties of the IGM that will be observed. However, at the same time, the urge to resolve the structures responsible for those processes is important. We introduce in this study a simulation framework of the cosmic dawn and reionization, based on hydrodynamics and radiative transfer code and a simple subgrid galaxy formation process for 1 Mpc-resolution simulations. Here, this model is calibrated on the state-of-the-art simulation CoDaII. This scheme permits us to follow consistently dark matter, hydrodynamics, and radiative transfer evolution on large scales, while the subgrid model deals with the galaxy formation scale, in particular, taking into account the different feedback on the star formation. We process the simulation to produce a simulated 21 cm signal as close as possible to the observations.
Summary
We invert for motions at the surface of Earth’s core under spatial and temporal constraints that depart from the mathematical smoothings usually employed to ensure spectral convergence of the ...flow solutions. Our spatial constraints are derived from geodynamo simulations. The model is advected in time using stochastic differential equations coherent with the occurrence of geomagnetic jerks. Together with a Kalman filter, these spatial and temporal constraints enable the estimation of core flows as a function of length and time-scales. From synthetic experiments, we find it crucial to account for subgrid errors to obtain an unbiased reconstruction. This is achieved through an augmented state approach. We show that a significant contribution from diffusion to the geomagnetic secular variation should be considered even on short periods, because diffusion is dynamically related to the rapidly changing flow below the core surface. Our method, applied to geophysical observations over the period 1950–2015, gives access to reasonable solutions in terms of misfit to the data. We highlight an important signature of diffusion in the Eastern equatorial area, where the eccentric westward gyre reaches low latitudes, in relation with important up/downwellings. Our results also confirm that the dipole decay, observed over the past decades, is primarily driven by advection processes. Our method allows us to provide probability densities for forecasts of the core flow and the secular variation.
SUMMARY
We make use of recent geodynamo simulations to propose a reduced stochastic model of the dynamics at the surface of Earth’s core. On decadal and longer periods, this model replicates the most ...energetic eigen directions of the geodynamo computation. Towards shorter timescales, it proposes a compensation for weaknesses of these simulations. This model furthermore accounts for the signature, in the geomagnetic secular variation, of errors of representativeness associated with unresolved processes. We incorporate the reduced stochastic model into a geomagnetic data assimilation algorithm—an augmented state ensemble Kalman filter—and apply it to re-analyse magnetic field changes over the period 1880–2015. Errors of representativeness appear to be responsible for an important fraction of the observed changes in the secular variation, as it is the case in the dynamo simulation.
Recovered core surface motions are primarily symmetric with respect to the equator. We observe the persistence of the eccentric westward gyre over the whole studied era and vortices that partly follow isocontours of the radial magnetic field at the core surface. Our flow models provide a good fit to decadal changes in the length-of-day and predict its interannual variations over the period 1940–2005. The largest core flow acceleration patterns are found in an equatorial belt below 10° in latitude and are associated with non-axisymmetric features. No systematic longitudinal drift of acceleration patterns is found, even over the past decades where satellite data are available. The acceleration of the high-latitude westward jet in the Pacific hemisphere is, during the satellite era, a factor 5 smaller than previously reported and its structure shows some evidence for equatorial asymmetry. The era of continuous satellite records provides enhanced contrast on the rapid core flow variations. The proposed assimilation algorithm offers the prospect of evaluating Earth-likeness of geodynamo simulations.
Inferring the core dynamics responsible for the observed geomagnetic secular variation requires knowledge of the magnetic field at the core‐mantle boundary together with its associated model ...covariances. However, most currently available field models have been built using regularization conditions, which force the expansions in the spatial and time domains to converge but also hinder the calculation of reliable second‐order statistics. To tackle this issue, we propose a stochastic approach that integrates, through time covariance functions, some prior information on the time evolution of the geomagnetic field. We consider the time series of spherical harmonic coefficients as realizations of a continuous and differentiable stochastic process. Our specific choice of process, such that it is not twice differentiable, mainly relies on two properties of magnetic observatory records (time spectra, existence of geomagnetic jerks). In addition, the required characteristic times for the low degree coefficients are obtained from available models of the magnetic field and its secular variation based on satellite data. We construct the new family COV‐OBS of field models spanning the observatory and satellite era of 1840–2010. These models include the external dipole and permit sharper time changes of the internal field compared to previous regularized reconstructions. The a posteriori covariance matrix displays correlations in both space and time, which should be accounted for through the secular variation error model in core flow inversions and geomagnetic data assimilation studies.
Key PointsModeling the geomagnetic secular vartiationUse of the stochastic processes frameworkObtain realistic model covariances
SUMMARY
We present an update of the geomagnetic data assimilation tool pygeodyn, use it to analyse ground and satellite-based geomagnetic data sets, and report new findings on the dynamics of the ...Earth’s outer core on interannual to decadal timescales. Our results support the idea that quasi-geostrophic Magneto-Coriolis waves, recently discovered at a period of 7 yr, also operate on both shorter and longer timescales, specifically in period bands centred around 3.5 and 15 yr. We revisit the source of interannual variations in the length of day and argue that both geostrophic torsional Alfvén waves and quasi-geostrophic Magneto-Coriolis waves can possibly contribute to spectral lines that have been isolated around 8.5 and 6 yr. A significant improvement to our ensemble Kalman filter algorithm comes from accounting for cross-correlations between variables of the state vector forecast, using the ‘Graphical lasso’ method to help stabilize the correlation matrices. This allows us to avoid spurious shrinkage of the model uncertainties while (i) conserving important information contained in off-diagonal elements of the forecast covariance matrix, and (ii) considering a limited number of realizations, thus reducing the computational cost. Our updated scheme also permits us to use observations either in the form of Gauss coefficient data or more directly as ground-based and satellite-based virtual observatory series. It is thanks to these advances that we are able to place global constraints on core dynamics even at short periods.
Geomagnetic observations from satellites have highlighted interannual variations in the rate of change of the magnetic field originating from Earth’s core. Downward continued to the core surface, ...these variations primarily show up in the equatorial belt. First, we recall the main characteristics of these patterns, addressing their spatio-temporal resolution, as seen from field models. We then review the several dynamical frameworks proposed so far to understand and model these observations, which populate the frequency spectrum on time scales close to the Alfvén time
τ
A
≈
2
yr, much shorter than the vortex turnover time
τ
U
≈
150
yr in Earth’s core. Magnetic–Archimedes–Coriolis (MAC) waves in a stratified layer below the core surface constitute a first possibility in the case of a sub-adiabatic heat flux at the top of the core. Their period may reach the interannual range for a layer thickness less than
≈
30
km, for a buoyancy frequency of the order of the Earth’s rotation rate. An alternative has been proposed in a context where the Coriolis force dominates the momentum balance, rendering transient motions almost invariant along the rotation axis (quasi-geostrophy, QG). Torsional Alfvén waves, consisting of axisymmetric QG motions, operate at periods similar to the Alfvén time, but are not sufficient to explain the interannual field changes, which require non-axisymmetric motions. QG Alfvén waves (involving the Coriolis and magnetic forces) constitute another possibility, with inertia playing an important role. They have been detected in the latest generation of geodynamo simulations, propagating in a ubiquitous manner at a speed slightly less than the Alfvén velocity. They are localized in longitude and as a result their description requires high azimuthal wavenumber. But the branch of QG waves with large extent in azimuth is also worth considering, as it reaches interannual periods as their radial wavenumber is increased. The excitation of such high-frequency dynamics is discussed with respect to the temporal spectrum of the core field, which presents a slope
∼
f
-
4
for periods approximately between
τ
A
and
τ
U
. We finally summarize the main geophysical implications of the existence of this interannual dynamics on core and lower mantle structure, properties, and dynamics.
SUMMARY
The inference of fluid motion below the core–mantle boundary from geomagnetic observations presents a highly non-unique inverse problem. We propose a new method that provides a unique local ...estimate of the velocity field, assuming quasi-geostrophic flow in the core interior (which implies equatorial mirror symmetry) and negligible magnetic diffusion. These assumptions remove the theoretical underdetermination, enabling us to invert for the flow at each point of a spherical grid representing the core surface. The unreliable reconstruction of small-scale flows, which arises because only large-scale observations are available, is mitigated by smoothing the locally estimated velocity field using a Gaussian process regression. Application of this method to synthetic data provided by a state-of-the-art geodynamo simulation suggests that using this approach, the large-scale flow pattern of the core surface flow can be well reconstructed, while the flow amplitude tends to be underestimated. We compare these results with a core flow inversion using a Bayesian framework that incorporates statistics from numerical geodynamo models as prior information. We find that whether the latter method provides a more accurate recovery of the reference flow than the local estimation depends heavily on how realistic/relevant the chosen prior information is. Application to real geomagnetic data shows that both methods are able to reproduce the main features found in previous core flow studies.