X-ray reverberation echoes are assumed to be produced in the strongly distorted spacetime around accreting supermassive black holes. This signal allows us to spatially map the geometry of the inner accretion flow - a region which cannot yet be spatially resolved by any telescope - and provides a direct measure of the black hole mass and spin. The reverberation timescale is set by the light travel path between the direct emission from a hot X-ray corona and the reprocessed emission from the inner edge of the accretion disc. However, there is an inherent degeneracy in the reverberation signal between black hole mass, inner disc radius and height of the illuminating corona above the disc. Here, we use a long X-ray observation of the highly-variable active galaxy, IRAS 13224-3809, to track the reverberation signal as the system evolves on timescales of a day. With the inclusion of all the relativistic effects, modelling reveals that the height of the X-ray corona increases with increasing luminosity, providing a dynamic view of the inner accretion region. This simultaneous modelling allows us to break the inherent degeneracies and obtain an independent timing-based estimate for the mass and spin of the black hole. The uncertainty on black hole mass is comparable to the leading optical reverberation method, making X-ray reverberation a powerful technique, particularly for sources with low optical variability.

We present a detailed X-ray timing analysis of the highly variable NLS1 galaxy, IRAS 13224-3809. The source was recently monitored for 1.5 Ms with XMM-Newton which, combined with 500 ks archival data, makes this the best studied NLS1 galaxy in X-rays to date. We apply standard time- and Fourier-domain in order to understand the underlying variability process. The source flux is not distributed lognormally, as would be expected for accreting sources. The first non-linear rms-flux relation for any accreting source in any waveband is found, with $\mathrm{rms} \propto \mathrm{flux}^{2/3}$. The light curves exhibit significant strong non-stationarity, in addition to that caused by the rms-flux relation, and are fractionally more variable at lower source flux. The power spectrum is estimated down to $\sim 10^{-7}$ Hz and consists of multiple peaked components: a low-frequency break at $\sim 10^{-5}$ Hz, with slope $\alpha < 1$ down to low frequencies; an additional component breaking at $\sim 10^{-3}$ Hz. Using the high-frequency break we estimate the black hole mass $M_\mathrm{BH} = [0.5-2] \times 10^{6} M_{\odot}$, and mass accretion rate in Eddington units, $\dot m_{\rm Edd} \gtrsim 1$. The non-stationarity is manifest in the PSD with the normalisation of the peaked components increasing with decreasing source flux, as well as the low-frequency peak moving to higher frequencies. We also detect a narrow coherent feature in the soft band PSD at $0.7$ mHz, modelled with a Lorentzian the feature has $Q \sim 8$ and an $\mathrm{rms} \sim 3$ %. We discuss the implication of these results for accretion of matter onto black holes.