We employed a surface-based laminar analysis pipeline to extract vertex-based resting state fMRI time series for 68 cortical regions separately from six layers of the neocortex (Fig. 2). Then, we performed vertex-by-vertex blind deconvolution [58] to get each vertex’s latent neural response and HRF. Please refer to the methods section for further details on these analyses.