(1) Aims: Atrial numerical modelling has generally represented the organ as either a surface or tissue with thickness. While surface models have significant computational advantages over tissue models, they cannot fully capture propagation patterns seen in vivo, such as dissociation of activity between endo-and epicardium. We introduce an intermediate representation, a bilayer model of the human atria, which is capable of recreating recorded activation patterns. (2) Methods: We simultaneously solved two surface monodomain problems by formalizing an optimization method to set a coupling coefficient between them. Two different asymptotically equivalent numerical implementations of the model are presented. We then built a geometrically and electrophysiologically detailed model of the human atria based on CT data, including two layers of fibre directions, major muscle bundles, and discrete atrial coupling. We adjusted parameters to recreate clinically measured activation times. Activation was compared to a monolayer model. (3) Results: Activation was fit to the physiological range measured over the entire atria. The crista terminalis and pectinate muscles were important for local right atrial activation, but did not significantly affect total activation time. Propagation in the bilayer model was similar to that of a monolayer, but with noticeable difference, due to three-dimensional propagationwhere fibre direction changed abruptly across the wall, resulting in a slight dissociation of activity. (4) Conclusion: Atrial structure plays the dominant role in determining activation. A bilayer model is able to take into account transmural heterogeneities, while maintaining thelow computational load associated with surface models.Read less <