Synchronization among arrays of beating cilia is one of the emergent phenomena in biological processes at meso-scopic scales. Strong inter-ciliary couplings modify the natural beating frequencies, $omega$, of individual cilia to produce a collective motion that moves around a group frequency $omega_m$. Here we study the thermodynamic cost of synchronizing cilia arrays by mapping their dynamics onto a generic phase oscillator model. The model suggests that upon synchronization the mean heat dissipation rate is decomposed into two contributions, dissipation from each ciliums own natural driving force and dissipation arising from the interaction with other cilia, the latter of which can be interpreted as the one produced by a potential with a time-dependent protocol in the framework of our model. The spontaneous phase-synchronization of beating dynamics of cilia induced by strong inter-ciliary coupling is always accompanied with a significant reduction of dissipation for the cilia population, suggesting that organisms as a whole expend less energy by attaining a temporal order. At the level of individual cilia, however, a population of cilia with $|omega|< omega_m$ expend more amount of energy upon synchronization.