We develop a theoretical description based on an existent mean-field model for the transient dynamics prior to steady flow of yielding materials. The mean-field model not only reproduces the experimentally observed non-linear time dependence of the shear-rate response to an external stress, but also allows for the determination of the different physical processes involved in the onset of the re-acceleration phase after the initial slowing down and a distinct fluidization phase. The fluidization time displays a power-law dependence on the distance of the applied stress to an age dependent yield stress, which is not universal but strongly dependent on initial conditions.Yield stress fluids (YSF), such as dense emulsions or pastes, display a rich rheological behavior that has attracted considerable attention recently [1,2]. Stationary flow is typically described by a nonlinear flow curve σ(γ), where σ is the stress andγ the deformation rate. However, the flow curve is far from accounting for the full complexity of these systems, that involves an interplay between external driving and internal aging, potentially leading to complex thixotropic behavior. In recent years, many experiments and molecular simulations have tried to reveal this complexity using creep experiments, in which the flow rate is measured in response to a fixed stress σ applied at a given waiting time t w after sample preparation [3][4][5][6][7][8][9][10]. These experiments, that lead, for σ larger than a yield stress σ Y , to flow or failure, reveal an intriguing behavior, with two salient features: (i) the strain-rateγ(t) in response to a stress larger than the yield stress is strongly non-linear and nonmonotonous, with a so called "s-shaped" dependence ofγ(t) [4,11,12], including a nontrivial "primary creep regime" often described by a power law t −µ .(ii) The fluidization time scale τ f diverges when approaching yield stress, however in a non-universal manner.In this work, we develop an approach that explains these features in athermal systems, in which thermal fluctuations have essentially no influence on the flow. These systems are a large subset of YSF [2], including e.g. foams, emulsions, physical gels, or granular media. In spite of the irrelevance of thermal fluctuations, the creep dynamics will depend on the initial condition determined by the preparation process and the subsequent waiting period, during which slow processes such as coarsening or compaction can alter the level of relaxation. Our approach contrasts previous attempts addressing systems in which thermal fluctuations are important, based on the soft glassy rheology model [13][14][15][16] or mode-coupling theory [17,18], Our description is based on a mean-field version of the elasto-plastic scenario, that describes the flow as resulting from interactions among local plastic events triggered by the external driving, and accounts for the flow properties of athermal YSF [19][20][21][22][23][24][25]. It extends a previous formulation for imposed shear-rate [26] to a system subjec...