Accretion across scales: streamers, surface-layer transport, and rapid replenishment in young protoplanetary discs

Protoplanetary discs evolve around newly-formed stars through an interplay of infall from surrounding turbulent cloud material, accretion towards the young star, and outflow driven mass-loss. It has been challenging to determine if discs are fed predominantly through infall along the disc midplane, or along the poles, and if accretion occurs in a steady or burst-like fashion. Here, we present a suite of 3D ideal magnetohydrodynamical simulations of protoplanetary disc formation and evolution in a dynamic, large-scale molecular cloud environment using the adaptive mesh refinement framework DISPATCH. We focus on nine stellar systems, where we resolve discs down to a scale of 0.8 au. Across the sample, stellar accretion proceeds at rates of about 10-5 M yr-1 over 10 5 yr, with significant variability. Discs grow to 100 au scales and remain gravitationally stable in time, with disc-to-star mass ratios below 10 %. Transient high-density streamers, with 10 kyr infall times, can drive anisotropic mass delivery at rates comparable to the background accretion flow. Their interaction with discs typically results in a temporary reduction of the disc size by half, and disc mass by 40 %. During later quiescent disc evolution stages (t 50 kyr), accretion predominantly occurs through the midplane and disc surface layers. This is associated with the development of a toroidal magnetic field morphology, which includes field reversals across both disc surfaces. In this way, the full disc mass reservoir is replenished on 10 kyr-timescales. These findings support that the outer parts of very young discs, when well-ionised and close to the ideal MHD regime, are not yet conducive to planet formation, due to high replenishment rates, strong turbulence, and disruptive streamer infall events.