Shearing Box Simulations Research Articles

TESTING CONVERGENCE FOR GLOBAL ACCRETION DISKS

Global disk simulations provide a powerful tool for investigating accretion and the underlying magnetohydrodynamic turbulence driven by the magneto-rotational instability (MRI). Using them to predict accurately quantities such as stress, accretion rate, and surface brightness profile requires that purely numerical effects, arising from both resolution and algorithm, be understood and controlled. We use the flux-conservative Athena code to conduct a series of experiments on disks having a variety of magnetic topologies to determine what constitutes adequate resolution. We develop and apply several resolution metrics: Qz and Qphi, the ratio of the grid zone size to the characteristic MRI wavelength, alpha_mag, the ratio of the Maxwell stress to the magnetic pressure, and the ratio of radial to toroidal magnetic field energy. For the initial conditions considered here, adequate resolution is characterized by Qz > 15, Qphi > 20, alpha_mag = 0.45, and a field energy ratio of 0.2. These values are associated with > 35 zones per scaleheight, a result consistent with shearing box simulations. Numerical algorithm is also important. Use of the HLLE flux solver or second-order interpolation can significantly degrade the effective resolution compared to the HLLD flux solver and third-order interpolation. Resolution at this standard can be achieved only with large numbers of grid zones, arranged in a fashion that matches the symmetries of the problem and the scientific goals of the simulation.

WIND-DRIVEN ACCRETION IN PROTOPLANETARY DISKS. II. RADIAL DEPENDENCE AND GLOBAL PICTURE

Non-ideal magnetohydrodynamical effects play a crucial role in determining the mechanism and efficiency of angular momentum transport as well as the level of turbulence in protoplanetary disks (PPDs), which are key to understanding PPD evolution and planet formation. It was shown in our previous work that at 1 AU, the magnetorotational instability (MRI) is completely suppressed when both Ohmic resistivity and ambipolar diffusion (AD) are taken into account, resulting in a laminar flow with accretion driven by magnetocentrifugal wind. In this work, we study the radial dependence of the laminar wind solution using local shearing-box simulations. Scaling relation on the angular momentum transport for the laminar wind is obtained, and we find that the wind-driven accretion rate can be approximated as M_dot~0.91x10^(-8)R_AU^(1.21)(B_z/10mG)^(0.93)M_Sun/yr, where B_z is the strength of the large-scale vertical magnetic field threading the disk. The result is independent of disk surface density. Four criteria are outlined for the existence of the laminar wind solution: 1). Ohmic resistivity dominated midplane region; 2). AD dominated disk upper layer; 3). Presence of (not too weak) net vertical magnetic flux. 4). Sufficiently well ionized gas beyond disk surface. All these criteria are likely to be met in the inner region of the disk from ~0.3 AU to about 5-10 AU for typical PPD accretion rates. Beyond this radius, angular momentum transport is likely to proceed due to a combination of the MRI and disk wind, and eventually dominated by the MRI (in the presence of strong AD) in the outer disk. Our simulation results provide key ingredients for a new paradigm on the accretion processes in PPDs.

EFFECTS OF LOCAL DISSIPATION PROFILES ON MAGNETIZED ACCRETION DISK SPECTRA

We present spectral calculations of non-LTE accretion disk models appropriate for high-luminosity stellar mass black hole X-ray binary systems. We first use a dissipation profile based on scaling the results of shearing box simulations of Hirose et al. to a range of annuli parameters. We simultaneously scale the effective temperature, orbital frequency, and surface density with luminosity and radius according to the standard α-model. This naturally brings increased dissipation to the disk surface layers (around the photospheres) at small radii and high luminosities. We find that the local spectrum transitions directly from a modified blackbody to a saturated Compton scattering spectrum as we increase the effective temperature and orbital frequency while decreasing midplane surface density. Next, we construct annuli models based on the parameters of a L/LEdd = 0.8 disk orbiting a 6.62 solar mass black hole using two modified dissipation profiles that explicitly put more dissipation per unit mass near the disk surface. The new dissipation profiles are qualitatively similar to the one found by Hirose et al., but produce strong near power-law spectral tails. Our models also include physically motivated magnetic acceleration support based once again on scaling the Hirose et al. results. We present three full-disk spectra, each based on one of the dissipation prescriptions. Our most aggressive dissipation profile results in a disk spectrum that is in approximate quantitative agreement with certain observations of the steep power-law spectral states from some black hole X-ray binaries.

WIND-DRIVEN ACCRETION IN PROTOPLANETARY DISKS. I. SUPPRESSION OF THE MAGNETOROTATIONAL INSTABILITY AND LAUNCHING OF THE MAGNETOCENTRIFUGAL WIND

We perform local, vertically stratified shearing-box MHD simulations of protoplanetary disks (PPDs) at a fiducial radius of 1 AU that take into account the effects of both Ohmic resistivity and ambipolar diffusion (AD). The magnetic diffusion coefficients are evaluated self-consistently from a look-up table based on equilibrium chemistry. We first show that the inclusion of AD dramatically changes the conventional picture of layered accretion. Without net vertical magnetic field, the system evolves into a toroidal field dominated configuration with extremely weak turbulence in the far-UV ionization layer that is far too inefficient to drive rapid accretion. In the presence of a weak net vertical field (plasma beta~10^5 at midplane), we find that the MRI is completely suppressed, resulting in a fully laminar flow throughout the vertical extent of the disk. A strong magnetocentrifugal wind is launched that efficiently carries away disk angular momentum and easily accounts for the observed accretion rate in PPDs. Moreover, under a physical disk wind geometry, all the accretion flow proceeds through a strong current layer with thickness of ~0.3H that is offset from disk midplane with radial velocity of up to 0.4 times the sound speed. Both Ohmic resistivity and AD are essential for the suppression of the MRI and wind launching. The efficiency of wind transport increases with increasing net vertical magnetic flux and the penetration depth of the FUV ionization. Our laminar wind solution has important implications on planet formation and global evolution of PPDs.

Local outflows from turbulent accretion disks

The aim of this paper is to investigate the properties of accretion disks threaded by a weak vertical magnetic field, with a particular focus on the interplay between MHD turbulence driven by the magnetorotational instability (MRI) and outflows that might be launched from the disk. For that purpose, we use a set of numerical simulations performed with the MHD code RAMSES in the framework of the shearing box model. We concentrate on the case of a rather weak vertical magnetic field such that the initial ratio beta0 of the thermal and magnetic pressures in the disk midplane equals 10^4. As reported recently, we find that MHD turbulence drives an efficient outflow out of the computational box. We demonstrate a strong sensitivity of that result to the box size: enlargements in the radial and vertical directions lead to a reduction of up to an order of magnitude in the mass-loss rate. Such a dependence prevents any realistic estimates of disk mass-loss rates being derived using shearing-box simulations. We find however that the flow morphology is robust and independent of the numerical details of the simulations. Its properties display some features and approximate invariants that are reminiscent of the Blandford & Payne launching mechanism, but differences exist. For the magnetic field strength considered in this paper, we also find that angular momentum transport is most likely dominated by MHD turbulence, the saturation of which scales with the magnetic Prandtl number, the ratio of viscosity and resistivity, in a way that is in good agreement with expectations based on unstratified simulations. This paper thus demonstrates for the first time that accretion disks can simultaneously exhibit MRI-driven MHD turbulence along with magneto-centrifugally accelerated outflows.

FROM DUST TO PLANETESIMALS: CRITERIA FOR GRAVITATIONAL INSTABILITY OF SMALL PARTICLES IN GAS

Dust particles sediment toward the midplanes of protoplanetary disks, forming dust-rich sublayers encased in gas. What densities must the particle sublayer attain before it can fragment by self-gravity? We describe various candidate threshold densities. One of these is the Roche density, which is that required for a strengthless satellite to resist tidal disruption by its primary. Another is the Toomre density, which is that required for de-stabilizing self-gravity to defeat the stabilizing influences of pressure and rotation. We show that for sublayers containing aerodynamically well-coupled dust, the Toomre density exceeds the Roche density by many (up to about 4) orders of magnitude. We present 3D shearing box simulations of self-gravitating, stratified, dust-gas mixtures to test which of the candidate thresholds is relevant for collapse. All our simulations indicate that the larger Toomre density is required for collapse. This result is sensible because sublayers are readily stabilized by pressure. Sound-crossing times for thin layers are easily shorter than free-fall times, and the effective sound speed in dust-gas suspensions decreases only weakly with the dust-to-gas ratio (as the inverse square root). Our findings assume that particles are small enough that their stopping times in gas are shorter than all other timescales. Relaxing this assumption may lower the threshold for gravitational collapse back down to the Roche criterion. In particular, if the particle stopping time becomes longer than the sound-crossing time, sublayers may lose pressure support and become gravitationally unstable.

Shearing box simulations of accretion disk winds

The launching process of a magnetically driven outflow from an accretion disk is investigated in a local, shearing box model which allows a study of the feedback between accretion and angular momentum loss. The mass-flux instability found in previous linear analyses of this problem is recovered in a series of 2D (axisymmetric) simulations in the MRI-stable (high magnetic field strength) regime. At low field strengths that are still sufficient to suppress MRI, the instability develops on a short radial length scale and saturates at a modest amplitude. At high field strengths, a long-wavelength "clump" instability of large amplitude is observed, with growth times of a few orbits. As speculated before, the unstable connection between disk and outflow may be relevant for the time dependence observed in jet-producing disks. The success of the simulations is due in a large part to the implementation of an effective wave-transmitting upper boundary condition.

INTERDEPENDENCE OF ELECTRIC DISCHARGE AND MAGNETOROTATIONAL INSTABILITY IN PROTOPLANETARY DISKS

We study how the magnetorotational instability (MRI) in protoplanetary disks is affected by the electric discharge caused by the electric field in the resistive MHD. We have performed three-dimensional shearing box simulations with various values of plasma beta and electrical breakdown models. We find the self-sustainment of the MRI in spite of the high resistivity. The instability gives rise to the large electric field that causes the electrical breakdown, and the breakdown maintains the high ionization degree required for the instability. The condition for this self-sustained MRI is set by the balance between the energy supply from the shearing motion and the energy consumed by the Ohmic dissipation. We apply the condition to various disk models and study where the active, self-sustained, and dead zones of MRI are. In the fiducial minimum-mass solar nebula (MMSN) model, the newly-found sustained zone occupies only the limited volume of the disk. In the late-phase gas-depleted disk models, however, the sustained zone occupies larger volume of the disk.

Outward migration of a super-Earth in a disc with outward propagating density waves excited by a giant planet

In this paper we consider a new mechanism for stopping the inward migration of a low-mass planet embedded in a gaseous protoplanetary disc. It operates when a low-mass planet (for example a super-Earth), encounters outgoing density waves excited by another source in the disc. This source could be a gas giant in an orbit interior to that of the low-mass planet. As the super-Earth passes through the wave field, angular momentum is transferred to the disc material and then communicated to the planet through coorbital dynamics, with the consequence that its inward migration can be halted or even reversed. We illustrate how the mechanism we consider works in a variety of different physical conditions employing global two-dimensional hydrodynamical calculations. We confirm our results by performing local shearing box simulations in which the super-Earth interacts with density waves excited by an independent harmonically varying potential. Finally, we discuss the constraints arising from the process considered here, on formation scenarios for systems containing a giant planet and lower mass planet in an outer orbit with a 2:1 commensurability such as GJ876.

DENSITY WAVES EXCITED BY LOW-MASS PLANETS IN PROTOPLANETARY DISKS. II. HIGH-RESOLUTION SIMULATIONS OF THE NONLINEAR REGIME

We investigate numerically the propagation of density waves excited by a low-mass planet in a protoplanetary disk in the nonlinear regime, using 2D local shearing box simulations with the grid-based code Athena at high spatial resolution (256 grid points per scale height h). The nonlinear evolution results in the wave steepening into a shock, causing damping and angular momentum transfer to the disk. On long timescales this leads to spatial redistribution of the disk density, causing migration feedback and potentially resulting in gap opening. Previous numerical studies concentrated on exploring these secondary phenomena as probes of the nonlinear wave evolution. Here we focus on exploring the evolution of the basic wave properties, such as its density profile evolution, shock formation, post-shock wave behavior, and provide comparison with analytical theory. The generation of potential vorticity at the shock is computed analytically and is subsequently verified by simulations and used to pinpoint the shock location. We confirm the theoretical relation between the shocking length and the planet mass (including the effect of the equation of state), and the post-shock decay of the angular momentum flux carried by the wave. The post-shock evolution of the wave profile is explored, and we quantitatively confirm its convergence to the theoretically expected N-wave shape. The accuracy of various numerical algorithms used to compute the nonlinear wave evolution is also investigated: we find that higher order spatial reconstruction and high resolution are crucial for capturing the shock formation correctly.

ASSESSING QUANTITATIVE RESULTS IN ACCRETION SIMULATIONS: FROM LOCAL TO GLOBAL

Discretized numerical simulations are a powerful tool for investigation of nonlinear MHD turbulence in accretion disks. However, confidence in their quantitative predictions requires a demonstration that further refinement of the spatial gridscale would not result in any significant change. This has yet to be accomplished, particularly for global disk simulations. In this paper, we combine data from previously published stratified shearing box simulations and new global disk simulations to calibrate several quantitative diagnostics by which one can estimate progress toward numerical convergence. Using these diagnostics, we find that the established criterion for an adequate numerical description of linear growth of the magneto-rotational instability (the number of cells across a wavelength of the fastest-growing vertical wavenumber mode) can be extended to a criterion for adequate description of nonlinear MHD disk turbulence, but the standard required is more stringent. We also find that azimuthal resolution, which has not often been extensively examined in previous studies, can significantly affect the evolution of the poloidal magnetic field. We further analyze the comparative resolution requirements of a small sample of initial magnetic field geometries; not surprisingly, more complicated initial field geometries require higher spatial resolution. Otherwise, they tend to evolve to qualitatively similar states if evolved for sufficient time. Applying our quantitative resolution criteria to a sample of previously published global simulations, we find that, with perhaps a single exception, they are significantly under-resolved, and therefore underestimate the magnetic turbulence and resulting stress levels throughout the accretion flow.

EFFECT OF AMBIPOLAR DIFFUSION ON THE NONLINEAR EVOLUTION OF MAGNETOROTATIONAL INSTABILITY IN WEAKLY IONIZED DISKS

We study the role of ambipolar diffusion (AD) on the non-linear evolution of the MRI in protoplanetary disks using the strong coupling limit, which applies when the electron recombination time is much shorter than the orbital time. The effect of AD in this limit is characterized by the dimensionless number Am, the frequency of which neutral particles collide with ions normalized to the orbital frequency. We perform three-dimensional unstratified shearing-box simulations of the MRI over a wide range of Am as well as different magnetic field strengths and geometries. The saturation level of the MRI turbulence depends on the magnetic geometry and increases with the net magnetic flux. There is an upper limit to the net flux for sustained turbulence, corresponding to the requirement that the most unstable vertical wavelength be less than the disk scale height. Correspondingly, at a given Am, there exists a maximum value of the turbulent stress alpha_max. For Am<1, the largest stress is associated with a field geometry that has both net vertical and toroidal flux. In this case, we confirm the results of linear analyses that show the fastest growing mode has a non-zero radial wave number with growth rate exceeding the pure vertical field case. We find there is a very tight correlation between the turbulent stress (alpha) and the plasma beta=P_gas/P_mag~1/(2alpha) at the saturated state of the MRI turbulence regardless of field geometry, and alpha_max rapidly decreases with decreasing Am. In particular, we quote alpha_max~0.007 for Am=1 and alpha_max~0.0006 for Am=0.1.

LOW-FREQUENCY OSCILLATIONS IN GLOBAL SIMULATIONS OF BLACK HOLE ACCRETION

We have identified the presence of large-scale, low-frequency dynamo cycles in a long-duration, global, magnetohydrodynamic (MHD) simulation of black hole accretion. Such cycles had been seen previously in local shearing box simulations, but we discuss their evolution over 1,500 inner disk orbits of a global pi/4 disk wedge spanning two orders of magnitude in radius and seven scale heights in elevation above/below the disk midplane. The observed cycles manifest themselves as oscillations in azimuthal magnetic field occupying a region that extends into a low-density corona several scale heights above the disk. The cycle frequencies are ten to twenty times lower than the local orbital frequency, making them potentially interesting sources of low-frequency variability when scaled to real astrophysical systems. Furthermore, power spectra derived from the full time series reveal that the cycles manifest themselves at discrete, narrow-band frequencies that often share power across broad radial ranges. We explore possible connections between these simulated cycles and observed low-frequency quasi-periodic oscillations (LFQPOs) in galactic black hole binary systems, finding that dynamo cycles have the appropriate frequencies and are located in a spatial region associated with X-ray emission in real systems. Derived observational proxies, however, fail to feature peaks with RMS amplitudes comparable to LFQPO observations, suggesting that further theoretical work and more sophisticated simulations will be required to form a complete theory of dynamo-driven LFQPOs. Nonetheless, this work clearly illustrates that global MHD dynamos exhibit quasi-periodic behavior on timescales much longer than those derived from test particle considerations.

DISSIPATION AND VERTICAL ENERGY TRANSPORT IN RADIATION-DOMINATED ACCRETION DISKS

Standard models of radiation supported accretion disks generally assume that diffusive radiation flux is solely responsible for vertical heat transport. This requires that heat must be generated at a critical rate per unit volume if the disk is to be in hydrostatic and thermal equilibrium. This raises the question of how heat is generated and how energy is transported in MHD turbulence. By analysis of a number of radiation/MHD stratified shearing-box simulations, we show that the divergence of the diffusive radiation flux is indeed capped at the critical rate, but deep inside the disk, substantial vertical energy flux is also carried by advection of radiation. Work done by radiation pressure is a significant part of the energy budget, and much of this work is dissipated later through damping by radiative diffusion. We show how this damping can be measured in the simulations, and identify its physical origins. Radiative damping accounts for as much as tens of percent of the total dissipation, and is the only realistic physical mechanism for dissipation of turbulence that can actually be resolved in numerical simulations of accretion disks. Buoyancy associated with dynamo-driven, highly magnetized, nearly-isobaric nonlinear slow magnetosonic fluctuations is responsible for the radiation advection flux, and also explains the persistent periodic magnetic upwelling seen at all values of the radiation to gas pressure ratio. The intimate connection between radiation advection and magnetic buoyancy is the first example we know of in astrophysics in which a dynamo has direct impact on the global energetics of a system.

Black Hole Accretion in Low States: Electron Heating

Plasmas in an accretion flow are heated by magneto-hydrodynamic turbulence generated through the magneto- rotational instability (MRI). The viscous stress driving the accretion is intimately connected to the microscopic processes of turbulence dissipation. We show that, in a few well-observed black hole accretion systems, there is compelling observational evidence of efficient electron heating by turbulence or collective plasma effects in low accretion states, when Coulomb collisions are not efficient enough to establish a thermal equilibrium between electrons and ions at small radii. We consider a Keplerian two-temperature accretion flow with a constant mass accretion rate in the pseudo- Newtonian gravitational potential and take into account the bremsstrahlung, synchrotron, and inverse Comptonization cooling processes. The balance of gravitational energy dissipation and turbulence energy cascade requires that the viscous stress be proportional to the product of the turbulence kinetic energy density and the total turbulence energy density, which may contradict the result of some shearing box simulations that the viscous stress is proportional to the magnetic field energy density. The critical mass accretion rate, below which the two-temperature solution may exist, is determined by the cooling processes and the collisional energy exchanges between electrons and ions and has very weak dependence on the collision-less heating of electrons by turbulence, which becomes more important at lower accretion rates. If the collision-less heating is dominated by the transit-time damping processes, small scale waves propagating obliquely with respect to the large scale magnetic field are prohibited, which may affect the saturate state of the MRI driven turbulence significantly. The plasma also needs to be strongly magnetized with the magnetic field and proton energy densities comparable so that electrons can share more of the dissipated gravitational energy. The heating of relativistic electrons is efficient since the heating rate is proportional to the mean momentum of the particles, and the electron heating may also be enhanced by their resonant scattering with small scale nearly parallel propagating waves.

Non-linear energy transfers in accretion discs MRI turbulence

The magnetorotational instability (MRI) is believed to be responsible for most of the angular momentum transport in accretion discs. However, molecular dissipation processes may drastically change the efficiency of MRI turbulence in realistic astrophysical situations. The physical origin of this dependency is still poorly understood as linear and quasi linear theories fail to explain it. In this paper, we look for the link between molecular dissipation processes and MRI transport of angular momentum in non stratified shearing box simulations including a mean vertical field. We show that magnetic helicity is unimportant in the model we consider. We perform a spectral analysis on the simulations tracking energy exchanges in spectral space when turbulence is fully developed. We find that the energy exchanges are essentially direct (from large to small scale) whereas some non linear interactions appear to be non local in spectral space. We speculate that these non local interactions are responsible for the correlation between turbulent transport and molecular dissipation. We argue that this correlation should then disappear when a significant scale separation is achieved and we discuss several methods by which one can test this hypothesis.

CHARACTERIZING THE TIME VARIABILITY IN MAGNETIZED NEUTRINO-COOLED ACCRETION DISKS: SIGNATURES OF THE GAMMA-RAY BURST CENTRAL ENGINE

The central engine of Gamma Ray Bursts is hidden from direct probing with photons mainly due to the high densities involved. Inferences on their properties are thus made from their cosmological setting, energetics, low-energy counterparts and variability. If GRBs are powered by hypercritical accretion onto compact objects, on small spatial scales the flow will exhibit fluctuations, which could in principle be reflected in the power output of the central engine and ultimately in the high energy prompt emission. Here we address this issue by characterizing the variability in neutrino cooled accretion flows through local shearing box simulations with magnetic fields, and then convolving them on a global scale with large scale dynamical simulations of accretion disks. The resulting signature is characteristic, and sensitive to the details of the cooling mechanism, providing in principle a discriminant for GRB central engine properties.

On the dynamics of planetesimals embedded in turbulent protoplanetary discs

(abridged) Angular momentum transport and accretion in protoplanetary discs are generally believed to be driven by MHD turbulence via the magneto-rotational instability (MRI). The dynamics of solid bodies embedded in such discs (dust grains, boulders, planetesimals and planets) may be strongly affected by the turbulence, such that the formation pathways for planetary systems are determined in part by the strength and spatial distribution of the turbulent flow. We examine the dynamics of planetesimals, with radii between 1m \^a 10 km, embedded in turbulent protoplanetary discs, using three dimensional MHD simulations. The planetesimals experience gas drag and stochastic gravitational forces due to the turbulent disc. We use, and compare the results from, local shearing box simulations and global models in this study. The main aims of this work are to examine: the growth, and possible saturation, of the velocity dispersion of embedded planetesimals as a function of their size and disc parameters; the rate of radial migration and diffusion of planetesimals; the conditions under which the results from shearing box and global simulations agree. We find good agreement between local and global simulations when shearing boxes of dimension 4H x 16H x 2H are used (H being the local scale height). The magnitude of the density fluctuations obtained is sensitive to the box size, due to the excitation and propagation of spiral density waves. This affects the stochastic forcing experienced by planetesimals. [...] Our models show that fully developed MHD turbulence in protoplanetary discs would have a destructive effect on embedded planetesimals. Relatively low levels of turbulence are required for traditional models of planetesimal accretion to operate, this being consistent with the existence of a dead zone in protoplanetary discs.

Trapping solids at the inner edge of the dead zone: 3-D global MHD simulations

The poorly-ionized interior of the protoplanetary disk is the location where dust coagulation processes may be most efficient. However even here, planetesimal formation may be limited by the loss of solid material through radial drift, and by collisional fragmentation of the particles. Our aim is to investigate the possibility that solid particles are trapped at local pressure maxima in the dynamically evolving disk. We perform the first 3-D global non-ideal MHD calculations of the disk treating the turbulence driven by the magneto-rotational instability. The domain contains an inner MRI-active region near the young star and an outer midplane dead zone, with the transition between the two modeled by a sharp increase in the magnetic diffusivity. The azimuthal magnetic fields generated in the active zone oscillate over time, changing sign about every 150 years. We thus observe the radial structure of the `butterfly pattern' seen previously in local shearing-box simulations. The mean magnetic field diffuses from the active zone into the dead zone, where the Reynolds stress nevertheless dominates. The greater total accretion stress in the active zone leads to a net reduction in the surface density, so that after 800 years an approximate steady state is reached in which a local radial maximum in the midplane pressure lies near the transition radius. We also observe the formation of density ridges within the active zone. The dead zone in our models possesses a mean magnetic field, significant Reynolds stresses and a steady local pressure maximum at the inner edge, where the outward migration of planetary embryos and the efficient trapping of solid material are possible.

A mean-field approach to the propagation of field patterns in stratified magnetorotational turbulence

Local shearing box simulations of stratified magneto rotational turbulence invariably exhibit cyclic field patterns which propagate away from the disc midplane. A common explanation for this is magnetic buoyancy. The recent analysis by Shi et al. however shows that the flow is buoyantly stable below one disc scale height H, necessitating an alternative explanation in this region. We here conduct and analyse direct numerical simulations to explain the observed behaviour by means of a mean-field description. Apart from the mean radial and azimuthal field, we monitor the small-scale current helicity, which we propose as a key indicator for saturation. Reconstructing the horizontally averaged field, we demonstrate that the problem can be reduced to a one-dimensional induction equation. By means of the so-called test field method, we then determine the underlying closure parameters. Our analysis shows that, apart from a possible direct MRI dynamo, two distinct indirect dynamo mechanisms operate in the disc. This resolves the issue of the "wrong" sign of the MRI dynamo effect. Finally, we use the obtained closure parameters to run a dynamically quenched dynamo model. This model approximately recovers the observed field patterns in the mean fields. Moreover, the model reproduces the prevailing parity and the distinct phase pattern in the small-scale current helicity. The latter property might open a potential route to understand the saturation of MRI induced turbulence.