Adarsh Kuruvanthodi, Daniel Schaerer, Rui Marques-Chaves, Andrea Weibel, Damien Korber, Corinne Charbonnel

The Balmer break (BB) is a key spectral feature for constraining stellar population ages, star formation histories, and redshifts of high-redshift sources. The redshift evolution and distribution of BB strength, together with the properties of BB galaxies, constrain stellar population characteristics and the nature of star formation across cosmic epochs. However, a systematic and unbiased characterization of BB strengths across the full galaxy population remains limited with the James Webb Space Telescope (JWST). We aim to characterize the redshift evolution of BB strength over $z=3.5$-10 and its distribution across different epochs using photometry. We also examine correlations between BB strength and key physical parameters within these redshift intervals. We further assess the implications of BB galaxies for the nature of star formation and stellar populations at $z>3.5$. We used the JWST NIRCam photometric observations taken as part of various programs, including CEERS, JADES, FRESCO, and PRIMER. We estimated the BB strength of the objects with two adjacent broadband filters in various redshift windows between redshifts 3 and 10, which exclude strong line contamination. We employ the SED-fitting code CIGALE for both SED fitting and the generation of mock galaxy simulations. We find that the median Balmer break strength (expressed as a flux ratio) increases from cosmic dawn to cosmic noon, from 1.1 to 1.5, primarily driven by the age of the stellar population. These estimations are in agreement with the latest spectroscopic estimations in the literature. We identify objects with extremely large BB strengths (BB$>3.0$) at $z=3.5$-4 and $z=7$-10, indicating strong extinction combined with an old stellar population and the presence of Little Red Dots (LRDs) in the former, and predominantly LRDs in the latter.[abridged]