Skip to main content
NIHPA Author Manuscripts logoLink to NIHPA Author Manuscripts
. Author manuscript; available in PMC: 2026 Oct 17.
Published before final editing as: Curr Biol. 2026 Oct 15:S0960-9822(26)01267-4. doi: 10.1016/j.cub.2026.09.058

Longitudinal Study of Toba/Qom Communities Reveals Delay and Shortening of Sleep in Real-Time

Leandro P Casiraghi 1,2,*, Ignacio Spiousas 1,2, Laura L Trebucq 1, M Florencia Coldeira 1,2, Camila R Godoy Peirone 1,†, Victor Y Zhang 3, Malen D Moyano 1, Diego A Golombek 1,2, Horacio O de la Iglesia 3,*
PMCID: PMC13674011  NIHMSID: NIHMS2213432  PMID: 42880258

Summary:

While artificial light and digital technologies are widely assumed to delay and curtail sleep, direct, longitudinal evidence documenting these changes is remarkably scarce. From 2012 to 2024, we conducted a longitudinal study of native Toba/Qom communities in northern Argentina, including rural groups newly introduced to electricity and semi-urban groups with earlier access to electricity. Using linear mixed-effects models on a dataset of over 12,000 actigraphically recorded sleep events from 156 participants at different times, we observed striking shifts in sleep dynamics. Remarkably, both groups showed delays of close to one hour in sleep timing, and participants from rural communities lost 35 minutes of sleep within the lapse of a decade. These changes occurred across a period that included, first, the introduction of electricity to rural communities after 2016, and, in more recent years, the spread of the internet and smartphone use among the Toba/Qom. While previous work in these and other populations has shown that both factors are associated with changes in sleep habits, our analyses suggest they cannot fully explain the reported shifts in sleep timing across time by themselves, which may therefore be consistent with concurrent broader processes of modernization and acculturation. These transformations highlight the rapid changes that can take place in communities going through the early stages of the adoption of modern technologies, offering a contemporary parallel to the shifts in human sleep hypothesized to have accompanied global industrialization.

eTOC Blurb:

Casiraghi et al. track sleep over 12 years in Toba/Qom communities in Argentina and find that sleep onset has been delayed by about an hour and sleep duration, in a rural community, shortened by 35 minutes. These shifts parallel electrification and later smartphone/internet access, linking sleep changes to broader modernization.

Introduction

The harnessing of electric power and then the invention of the electric bulb revolutionized human life in ways that are impossible to quantify. Electric light further allowed for the extension of wakefulness and chores into the evening hours and, eventually, the birth of the 24/7 society of work and entertainment in the current hyper-industrialized world. This process has significantly delayed sleep times, but whether it shortened sleep duration is still a matter of debate. A systematic review of studies reporting the duration of sleep between the 1960s and the 2000s failed to find a consistent shortening of sleep,1 and a recent analysis of published sleep data suggests that in industrialized societies sleep may be longer than in non-industrialized ones.2 In contrast, a poll by Gallup3 and a study on self-reported sleep suggest that in the last two decades sleep has shortened in the US.4 Similarly, a systematic review of available children’s data from 1905 to 2008 suggests a decrease in daily sleep of one hour in this period.5

In order to study the ways in which industrialization has affected human behavior, we need to be able to study populations in which this process is still ongoing. Pre-industrial societies offer unique opportunities to investigate the impact of the adoption of modern technologies and their related habits on sleep patterns, as they usually display differential access to electricity, high levels of exposure to the natural environment, weaker scheduling impositions, and rely more on environmental and ecological time references.6–9

Our group is focused on the study of native Western Toba/Qom populations in the northern region of Argentina.10 These people, formerly hunter-gatherers, remain a very isolated and underserved group spread in the Formosa province, relying mostly on government assistance for subsistence. For the last 14 years, we have been studying sleep and circadian rhythms of both rural and semi-urban Toba/Qom communities (referred to as “urban”, for the sake of simplicity; see “Toba/Qom participants” in Methods), allowing us to contrast the behavior of communities in which the most relevant distinguishing characteristic was, until not so long ago, their access to electricity. Our previous work has shown that, indeed, urban Toba/Qom displayed later sleep onsets and shorter sleep duration than rural communities with no access to electricity.11 This correlation between electric light and delayed sleep has also been reported by several other groups studying pre-industrial communities around the world.12–15

Since we started working with Toba/Qom communities in 2012, electricity has been progressively introduced to most rural areas in this region of the country, including the most isolated settlements (Figure 1A). Since such efforts started around 2016, all the communities we have worked with are now connected to the electric grid, and at the time of writing this paper, every single family we have worked with has electricity in their homes. In recent years, television sets, smartphones, and internet connections, while in limited supply, are already available to many families.

Figure 1. Evolution of night light pollution and sleep from 2012 to 2024 in both rural and urban areas.

Figure 1.

A) Radiance DNB levels measured by NASA’s VIIRS instrument in the geographic regions of Ingeniero Juárez (bottom-left arrow) and Vaca Perdida (top-right arrow) over the years, obtained from Jurij Stare’s Light Trends application. B) Summary of sleep start times and duration recorded through ambulatory wrist-actigraphy. Each marker indicates the mean time of sleep start for a participant at a single data collection campaign, with gray lines connecting repeated measurements. Fill colors represent the mean sleep duration, and shapes indicate the participant’s group. See also Table S2.

We have longitudinally monitored sleep patterns through wrist-actigraphy in these Toba/Qom communities at different seasons and time points since 2012, generating a dataset of over 12,000 sleep records. This dataset represents a unique opportunity to assess the impact of the introduction of electricity on the sleep of these communities in real-time. The study we present here provides a close insight into the history of human sleep patterns in developing pre-industrial communities, highlighting the profound impact of modern technology on human behavior.

Results

Our sleep dataset from Toba/Qom communities includes over 12,000 sleep events recorded through wrist-actigraphy between 2012 and 2024 from a total of 156 participants (93 in the rural group, 62 in the urban group, and one who moved from one community to the other after 2018). Table S1 details the number of participants, recorded events, and age and sex composition of all data according to community, and the dates when it was collected. Starting around the year 2016, electricity was progressively introduced in the rural communities in the Vaca Perdida and neighboring areas, including the communities we studied (Figure 1A). Throughout this period of significant transformations for the rural Toba/Qom, we observed a progressive delay and shortening of sleep both at the individual and community level (Figure 1B and Table S2).

Sleep start, midpoint, and end times, sleep duration, and sleep regularity (measured by the Sleep Regularity Index) were analyzed through linear mixed-effects models including groups (urban or rural), date, photoperiod length, type of day (workday or free-day), and sex and age as fixed effects, while including random intercepts and random slopes across time for each participant. Interaction terms between group and date, and between sex and age were also included in the models. Due to differences in the actigraphic devices employed across time, we implemented a non-parametric multiple imputation strategy utilizing empirical error distributions based on a side-by-side comparison of Philips Actiwatch Spectrum Plus and Axivity AX3 devices. The results reported are therefore based on bootstrapping of one hundred independent simulations, effectively incorporating this variability as added imprecision in our study (see Methods for the complete details on the statistical analysis and for additional considerations on the older Actiwatch-L model).

Table 1 displays the full results of models for sleep onset, offset, and midpoint; Table 2 shows the results for sleep duration and sleep regularity. Coefficients are expressed along with their 95% CI. The models revealed photoperiod length as a strong predictor of sleep times, with 0.20 [0.17–0.23] h delays in sleep start times (t = 12.244, p < 0.001), 0.16 [0.13–0.19] h advances in sleep end times (t = 10.557, p < 0.001), and a reduction of 0.36 [0.33–0.40] h of sleep (t = 19.677, p < 0.001) per hour of increase in daylight. A longer photoperiod was also associated with slightly higher sleep regularity (1.0 [0.4–1.7] units/h of photoperiod extension, t = 3.097, p = 0.002). Weekends displayed slightly later sleep start times (0.18 [0.12–0.25] h, t = 5.623, p < 0.001) and later sleep offset times (0.21 [0.15–0.26] h, t = 7.055, p < 0.001), and thus did not affect sleep duration. Accordingly, the midpoint of sleep was minimally delayed in nights before free days (0.19 [0.15–0.24] h, t = 7.907, p < 0.001), indicating that social jet lag is minimal in these communities. Male participants started their sleep significantly later (0.59 [0.33–0.84] h, t = 4.438, p < 0.001), woke up later (0.46 [0.20–0.73] h, t = 3.446, p < 0.001), and had a lower sleep regularity (-2.85 [-4.88 to −0.82] SRI units, t = 2.762, p = 0.006), but slept about the same amount as female participants. Consistent with previous findings in these communities, the rural group showed overall earlier sleep start times than the urban group (0.32 [0.05, 0.6] h, t = 2.303, p = 0.021).

Table 1.

Linear models including random effects of sleep start, end, and midpoint times, expressed as hours from 00:00 on the date of the start of sleep.

Parameter Coefficient SE 95% CI t p
Sleep onset (Intercept) 23.05 (1) 0.11 [22.84, 23.26] 213.178 <0.001
group = Urban 0.32 0.14 [0.05, 0.6] 2.303 0.021
date (years) 0.09 0.02 [0.05, 0.13] 4.473 <0.001
photop. length (h) 0.20 0.02 [0.17, 0.23] 12.244 <0.001
day = weekday −0.18 0.03 [−0.25, −0.12] −5.623 <0.001
sex = male 0.59 0.13 [0.33, 0.84] 4.438 <0.001
age^1 13.69 11.81 [−9.45, 36.83] 1.159 0.246
age^2 0.71 15.65 [−29.97, 31.39] 0.045 0.964
age^3 28.75 15.23 [−1.1, 58.6] 1.888 0.059
(group = Urban):date 0.00 0.03 [−0.06, 0.06] 0.043 0.966
(sex = male):age^1 17.46 14.59 [−11.14, 46.06] 1.197 0.231
(sex = male):age^2 −36.48 16.97 [−69.75, −3.21] −2.149 0.032
(sex = male):age^3 −8.44 16.35 [−40.5, 23.61] −0.516 0.606
Sleep offset (Intercept) 31.21 (2) 0.11 [31, 31.42] 293.569 <0.001
group = Urban 0.03 0.13 [−0.24, 0.29] 0.209 0.835
date (years) 0.03 0.02 [−0.01, 0.07] 1.550 0.121
photop. length (h) −0.16 0.01 [−0.19, −0.13] −10.557 <0.001
day = weekday −0.21 0.03 [−0.26, −0.15] −7.055 <0.001
sex = male 0.46 0.13 [0.2, 0.73] 3.446 <0.001
age^1 −16.85 11.86 [−40.1, 6.4] −1.420 0.156
age^2 14.71 15.00 [−14.7, 44.12] 0.980 0.327
age^3 25.90 14.93 [−3.37, 55.16] 1.735 0.083
(group = Urban):date 0.06 0.03 [−0.01, 0.12] 1.656 0.098
(sex = male):age^1 7.04 14.66 [−21.7, 35.79] 0.480 0.631
(sex = male):age^2 −47.60 16.33 [−79.61, −15.59] −2.915 0.004
(sex = male):age^3 −6.94 16.07 [−38.44, 24.56] −0.432 0.666
Sleep midpoint (Intercept) 27.13 (3) 0.10 [26.93, 27.33] 266.239 <0.001
group = Urban 0.17 0.13 [−0.08, 0.43] 1.313 0.189
date (years) 0.06 0.02 [0.02, 0.1] 3.158 0.002
photop. length (h) 0.02 0.01 [0, 0.05] 1.634 0.102
day = weekday −0.19 0.02 [−0.24, −0.15] −7.907 <0.001
sex = male 0.52 0.13 [0.27, 0.77] 4.085 <0.001
age^1 −2.86 11.23 [−24.86, 19.15] −0.254 0.799
age^2 10.12 14.15 [−17.62, 37.85] 0.715 0.475
age^3 31.71 13.87 [4.53, 58.89] 2.286 0.022
(group = Urban):date 0.03 0.03 [−0.03, 0.09] 1.039 0.299
(sex = male):age^1 14.85 13.92 [−12.42, 42.13] 1.067 0.286
(sex = male):age^2 −45.56 15.41 [−75.77, −15.34] −2.956 0.003
(sex = male):age^3 −10.97 14.92 [−40.21, 18.28] −0.735 0.462

Intercept levels: age = 22 y.o., photop. length = 12 h, group = Rural, date = 2018/7/1, sex = female, type of day = free day. Clock time conversion: (1) 23:03; (2) 7:13; (3) 3:08.

Table 2.

Linear models including random effects of sleep duration (h) and sleep regularity (SRI units)

Parameter Coefficient SE 95% CI t p
Sleep duration (Intercept) 8.14 0.08 [7.99, 8.3] 102.936 <0.001
group = Urban −0.26 0.10 [−0.46, −0.05] −2.459 0.014
date (years) −0.06 0.01 [−0.09, −0.03] −4.132 <0.001
photop. length (h) −0.36 0.02 [−0.4, −0.33] −19.677 <0.001
day = weekday −0.02 0.04 [−0.1, 0.05] −0.600 0.548
sex = male −0.02 0.09 [−0.19, 0.16] −0.216 0.829
age −0.03 0.01 [−0.05, −0.02] −4.599 <0.001
(group = Urban):date 0.06 0.02 [0.01, 0.1] 2.429 0.015
(sex = male):age −0.01 0.01 [−0.03, 0.01] −1.228 0.219
SleepRegularityIndex (Intercept) 77.29 0.76 75.8 to 78.77 102.065 <0.001
group = Urban −2.85 1.03 −4.88 to −0.82 −2.762 0.006
date (years) −0.52 0.18 −0.88 to −0.17 −2.891 0.004
photop. length (h) 1.03 0.33 0.38 to 1.69 3.097 0.002
sex = male −3.55 0.93 −5.38 to −1.73 −3.818 <0.001
age 0.01 0.08 −0.13 to 0.16 0.198 0.843
block length 0.26 0.03 0.2 to 0.32 8.647 <0.001
(group = Urban):date −0.49 0.30 −1.09 to 0.1 −1.635 0.103
(sex = male):age 0.00 0.10 −0.19 to 0.19 0.045 0.964

Intercept levels: age = 22 y.o., photop. length = 12 h, group = Rural, date = 2018/7/1, sex = female, type of day = free day (does not apply to SRI).

Figure 2 and Table S3 present the model-predicted changes in sleep variables for each group under a 12-h photoperiod length for a 22-year-old (the median of ages across all campaigns), while averaging for sex and type of day, for the 2014 to 2024 ten-year span. Models indicate similar delays in the times of start of sleep of 0.88 [0.49–1.27] h (t = 4.444, p < 0.001) for the rural group and of 0.90 [0.40–1.40] h (t = 3.528, p < 0.001) for the urban group. Evolution of sleep end times was asymmetrical between groups: while the rural group showed only a marginal, non-significant change in wake-up times over ten years (0.32 [-0.09 to 0.73] h, t = 1.540, p = 0.124), the urban participants delayed theirs substantially (0.87 [0.36–1.39] h, t = 3.324, p < 0.001). Consistently with these asymmetric shifts, changes in sleep duration were different between communities. While the rural group displayed the loss of 35 minutes of sleep (0.58 [0.30–0.86] h, t = 4.105, p < 0.001), this reduction was absent in the urban group. This asymmetry in the evolution of sleep times also implied that the urban group showed the larger shift in their midpoint of sleep, with a delay of 0.93 [0.45–1.41] h (t = 3.808, p < 0.001), as compared to the 0.61 [0.23–0.99] h delay of the rural group (t = 3.138, p = 0.002). Sleep regularity, measured by the Sleep Regularity Index, showed a significant decrease in both groups across time. Between 2014 and 2024, the model predicts a loss of 5.2 [1.6–8.8] units in the rural population (t = 2.876, p = 0.005) and of 10.2 [5.3–15.0] units in the urban group (t = 4.148, p < 0.001).

Figure 2. Predictions of sleep patterns for 2014 and 2024 based on linear models including random effects.

Figure 2.

Vertical dashed lines were drawn for the start, end, and midpoint of sleep, and horizontal black lines represent the standard error of these predictions. Predicted sleep duration in hours is also indicated within each box. Predicted differences in minutes are indicated for all sleep variables, and asterisks indicate statistically significant changes across the modelled 2014–2024 span (see Table S3 for the full results of the analyses). The yellow vertical line represents the approximate time of dawn. * p<0.05, ** p<0.01, *** p<0.001 (Tukey-corrected for multiple comparisons). Predictions were calculated for a 12 h photoperiod, averaged for sex and type of day, and with age set at 22 y.o., the median across all studies.

The 2012 to 2024 period involved important transformations in these Toba/Qom communities. First, starting around 2016, electricity was progressively introduced in rural communities. Secondly, while to our knowledge not a single participant from either community under study, nor anyone in their families, owned a smartphone by 2018, about half of the participants in our studies had their own smartphone, and about the same had access to a TV and/or a computer at their homes by October 2024. Moreover, even in very remote rural areas, where cell phone coverage is extremely limited, some participants can access the internet through free public Wi-Fi hotspots.

We set out to test whether either of these two events (i.e. introduction of electricity, and access to digital connectivity) could support the observed changes in sleep timing across this timespan, based on their potential for extending activity in the evening hours and thus delaying the start of sleep. To do this, we fitted a new linear model in which every sleep event was labeled as recorded either under conditions of access to electricity or not, and whether it was recorded during the “smartphones and internet era,” an umbrella term for all data from 2023 onwards (see “Test of candidate mechanisms” in Methods for further details on this analysis). This model revealed that neither the introduction of electricity into rural areas nor the “smartphone era” can explain by themselves the longitudinal changes in sleep onset observed in these communities, which were still consistent at an estimated delay of 0.12 [0.05–0.19] h per year (t = 3.321, p < 0.001; Table S4). Because both electricity access and the “smartphone era” indicators are, by construction, correlated with calendar date (see “Test of candidate mechanisms” in Methods), we cannot claim that this analysis fully rules out a contribution of either factor to the observed changes in sleep, particularly since both are independently plausible drivers of sleep timing, and a non-significant coefficient under any degree of shared variance with date does not, by itself, demonstrate the absence of an effect. Our diagnostics indicate, however, that this correlation is mild rather than severe, and that the residual date effect does not attenuate once electricity and smartphone use are included in the model (if anything, it remains of a similar, even slightly larger magnitude, and stays significant). We therefore interpret these results as evidence that a longitudinal trend in sleep persists beyond what can be attributed to these two milestones in the “modernization” of Toba/Qom lifestyle.

Discussion

Our analysis reveals the rapid evolution of sleep habits in two Toba/Qom communities undergoing early stages of modernization, marked by the sequential introduction of electricity, mobile devices, and internet access. Drawing on over 12,000 recorded sleep events spanning twelve years, our models document significant shifts in sleep timing and duration: in the span of a single decade, sleep onset was delayed by close to one hour, and total daily sleep was reduced by 35 minutes in rural communities. These findings echo patterns reported in historical records16 and self-reported sleep studies4,5 tracking the transition from preindustrial to industrial (and more recently digital or postindustrial) environments. Yet, longitudinal investigations using objective sleep measures in communities where modernization, understood as a set of co-occurring environmental and social changes, unfolds incrementally have remained absent from the literature. Our results are consistent with the hypothesis that such transformations delay and shorten human sleep, and provide the longitudinal, objectively measured evidence of change over time that this hypothesis has so far lacked. Importantly, our findings imply there are other factors driving sleep changes beyond the availability of electricity and the use of electric light, and even smartphones, computers, or TV sets. Together, these results highlight the need to monitor how the introduction of modern amenities reshapes the lifestyle and biology of communities in transition.

In the time we studied these communities, we detected similar delays in the time of sleep start across groups, but only the urban community displayed a comparable delay in the offset of sleep. Rural communities, who live under even more limited conditions than their counterparts in the city of Ingeniero Juárez, may not be able to extend their sleep much later in the morning either due to being more exposed to natural elements (sunlight, temperature), to the need to start with sustenance activities (e.g. harvesting firewood, managing livestock, gathering), or to other necessities related to their isolated rural environments (e.g. distance to health centers or convenience shops, lack of transportation options). This contrast may explain the reduction in sleep duration in the rural group in the face of similar shifts in the timing of sleep onset. That is, the large shift in sleep start times appears to be the most relevant culprit driving changes in overall sleep behavior in these communities.

Internet access and devices such as smartphones, however limited in supply, allow Toba/Qom communities to connect to a world beyond their local environment, for example, through social media and video platforms. The adoption of these technologies in these communities has been very rapid. To the best of our knowledge, smartphones were virtually unavailable by 2018, yet by 2024, half of the Toba/Qom participants owned a smartphone, and about the same had access to a TV or a computer in their homes, reflecting a swift shift toward digital connectivity. A public internet antenna was installed as recently as 2022 in Isla García, providing free internet access to the community; one family in Vaca Perdida even had a satellite internet connection. These technological advances likely favor the delay of sleep and, whenever the time for waking up cannot be further delayed, the reduction of sleeping time by providing almost unlimited entertainment for the night hours. However, our models indicate that the introduction of smartphones and the internet cannot explain the observed delays in the start of sleep by itself. Moreover, this delay in sleep times could not be driven merely by late-night smartphone users. A recent study on these communities revealed that users of smartphones display only a 13-minute delay of sleep compared to non-users,17 making it impossible for just these users to explain the population-level changes reported here. Electronic devices are hypothesized to impact sleep via several simultaneous mechanisms: direct reduction of sleep time, effects of emitted light, promotion of arousal, and causing sleep disruptions.18 For the Toba/Qom population, an additional factor may contribute beyond these established pathways. Electronic devices such as smartphones, computers, and TVs, along with access to the internet, can open a door into the “modern” lifestyle, of which nightlife is a central element. This is particularly relevant in Argentina, where, for instance, the TV “prime time” extends from 21:00 to midnight and football matches can start as late as 22:00.19 This phenomenon merits further research incorporating Toba/Qom self-perceived experiences regarding access to electricity and electric light sources, as well as the roles of the internet, the TV, and other electronic devices.

Sleep loss is considered a modern world epidemic with significant consequences for public health, performance and safety, learning and memory, and overall economic impact.20–23 Lack of sleep can affect health at several levels, including, but not limited to, cardiovascular, metabolic, immune, and mental health,24,25 and is a strong predictor of all-cause mortality.26 Although our models predict an average sleep duration still within the consensus recommended sleep duration range for these ages,27 we note that these recommendations are aimed at people in modern developed societies. More research is needed to assess whether this reduction of sleep time in this population is associated with changes in physical health, psychological wellbeing, and overall quality of life.

The observed reduction in sleep regularity in both communities may be more concerning. Sleep regularity has gained attention recently as a marker of sleep health,28 and has been proposed as a more precise predictor of the negative impacts of poor sleep than sleep duration itself.29 Delayed sleep has already been associated with low sleep regularity in hyper-industrialized societies,30 and in a study combining data from adolescents in the Toba/Qom and in rural communities of Southern Brazil, we already reported that sleep regularity was reduced in students exposed to electric light in the evening hours.14 Further investigation should determine whether this change is driven by the mere access to electric light and other electric appliances and devices at home, and/or if it is correlated with changes in evening habits that involve the whole community.

Numerous studies of preindustrial communities support the notion that electricity and other postindustrial environmental factors impact sleep timing.7 However, the hypothesis that industrialization has necessarily shortened sleep remains under debate. A recent meta-analysis, for instance, suggests that people in industrialized societies sleep longer than those in preindustrial settings, and attributes modern sleep-related health issues not to shorter sleep but to heightened circadian disruption.2 However insightful, these observations need to be carefully weighed, considering the inherent diversity of pre-industrial and industrialized groups, without overlooking the unique environmental and sociocultural factors that can shape sleep. Our longitudinal study of the Toba/Qom communities provides a rare opportunity to assess the real-time impact of certain features of modernization, including the acquisition of electricity and related technologies. Our data reveals that in these Toba/Qom communities, this process correlates with a significant reduction in sleep duration and an increased circadian disruption —marked by a reduced sleep regularity and the further decoupling of sleep onset from sunset. We should also note that the transition from harsh preindustrial conditions to an electrically lit industrial environment can also have beneficial effects on sleep, for example, by improving comfort and sense of safety at bedtime.31 This underscores the intricate ways in which the environment and technology interact to shape human sleep patterns.32

Some limitations in our study need to be addressed. First, it should be noted that the recordings in our dataset have been obtained using different devices and algorithms. Because of the long timespan of our study, this was an obstacle hard to overcome. Of note, the Actiwatch devices used from 2012 to 2018 do not provide raw accelerometric data, as the Axivity AX3 devices do, precluding us from using a common algorithm to extract sleep timing. These different methods are all based on longitudinal wrist-accelerometry and the sleep estimation in each of them was performed following well-tested protocols, based on similar criteria of activity/inactivity and aided by sleep logs to detect unreliable estimations.

Our side-by-side tests indicate that average sleep time estimations from AX3 and Actiwatch Spectrum devices are comparable within a margin of less than 20 minutes of difference between them (see Figure S2). However, because a discrepancy of this magnitude is biologically relevant, and because differences between actigraphic methods may even vary between individuals with different activity profiles,33 we followed a rigorous procedure to integrate the available data by using multiple imputation of errors to account for the empirical disparities between devices. By presenting results pooled from one hundred simulated models for each sleep variable, this approach provides strong confidence that the detected longitudinal differences that survive the addition of this imprecision into our statistical analysis reflect a genuine temporal phenomenon rather than an artifact of changing instrumentation. Our correction method resamples from the empirical distribution of 206 paired differences obtained from non-Toba/Qom participants who wore both devices simultaneously, thus propagating the location, dispersion, and skew of this discrepancy into our uncertainty estimates. The assumption, nevertheless, is that such error can be comparable across groups. This is a limitation that cannot be overcome with the data available to us, and that demands further precautions in the interpretation of our results. Finally, because we could not perform a specific side-by-side validation of the older Actiwatch-L devices, we performed an additional sensitivity analysis over a reduced dataset excluding all data obtained with these devices, which revealed essentially the same results (Table S5).

Second, thanks to the improvement of fieldwork conditions, in terms of the available equipment and our relationship with the communities, our dataset is much denser in recent years, compared to the first four years of work. While this may pose a limitation in the model’s precision for those early years, the use of linear models, including random effects accounting for sex, age, and repeated measures, is robust enough to assess the dynamics of sleep across the studied period in these communities. Importantly, the lower number of participants in the rural community before the appearance of electricity may have contributed to the lack of a statistically significant effect of electricity on sleep onset in our analysis. Obviously, this result does not mean that the appearance of electricity was not critical for the sleep delay and shortening we observed; electricity not only allows for artificially lit evenings but also for internet hotspots, cell phones, TVs, and other sleep disruptors.

Finally, it is worth addressing the five-year interruption in data collection between 2018 and 2023. This covers the time of the COVID-19 pandemic, which of course impacted on the possibility of continuing fieldwork. The pandemic was associated with substantial changes in sleep in populations worldwide, including in Argentina, generally involving delayed and lengthened sleep together with reduced social jetlag.34–37 However, the mechanisms proposed to drive pandemic-related sleep change (i.e., confinement indoors and consequent reduction of natural light exposure, the relaxation or pausing of fixed work and school schedules, the elimination of commuting, and the adoption of “working from home” routines) are largely tied to urban, densely populated, and wage-employed ways of life. These conditions are substantially attenuated, if not absent, in the Toba/Qom communities we study, where subsistence activities are largely outdoors, commuting is very infrequent, and formal employment is rare. Moreover, our results show the opposite trend in sleep duration to that reported during the pandemic, so the only dimension on which a pandemic legacy could plausibly act in the direction we observe is the delay in sleep onset. An additional sub analysis of sleep start times in the 2012 to 2018 data determined that the trend of sleep delay was already present at a comparable, albeit numerically larger, rate with overlapping 95% confidence intervals (Table S6), making it unlikely that the longitudinal trend we report is driven by an acute change during those years.

Our work suggests changes in sleep habits as one of the most plausible consequences and multifaceted impacts of technological innovation on human behavior and health. These impacts could differentially affect underserved communities like the Toba/Qom, who have finally attained the benefits of the access to electricity, smartphones, and the internet, but remain significantly behind in the access to many of the advantages of modern, developed societies in other dimensions such as infrastructure and health services. This disparity may exacerbate any potential negative consequences of the observed delayed and shortened sleep in these communities, which are commonly reported features of modern lifestyle. In this sense, the Toba/Qom case we present exemplifies what has been described as “cultural arrhythmia”, i.e., an imbalance between the pace of changes in human lives and their surrounding habitats.38 Consequently, it is imperative that such transformative processes consider the unique features and conditions of each community, ensuring they enhance health and overall wellbeing while proactively addressing any potential challenges. Our results highlight the need to monitor sleep health in groups amid such rapid environmental transitions, the full implications of which remain to be established.

Resource availability

Lead contact:

Requests for further information and resources should be directed to and will be fulfilled by the lead contact, Horacio de la Iglesia (horaciod@uw.edu).

Material availability:

This study did not generate new unique reagents.

Data and code availability:

  • Due to the vulnerable nature of these communities, the longitudinal sleep data from Toba/Qom participants will be available to researchers upon reasonable request to the lead contact. Requests for de-identified, individual-level data should detail the proposed uses of the data and research questions, procedures that will be employed to ensure data security, any potential benefits to the study communities, and precautions for avoiding stigmatizing interpretations of the results, and its use may require a data use agreement.

  • Data used in the cross-correlation of AX3 and Actiwatch devices along with a set of scripts in the R language that allows users to replicate the linear model analysis including the multiple imputation of errors, with the application of Rubin’s Rules, are available at https://github.com/LeandroCasiraghi/ax3-aw-multiple-error-imputation.

Any additional information required to reanalyze the data reported in this paper is available from the lead contact upon proper request.

STAR Methods

EXPERIMENTAL MODEL AND STUDY PARTICIPANT DETAILS

Toba/Qom participants

All the activities including Toba/Qom participants described in this work, past and recent, were reviewed and approved by the Internal Review Board from the University of Washington’s Human Subject Division (Study IDs: #43876, #52422, #11222, and #15158) and followed the Declaration of Helsinki. The research did not involve any risk to participants or researchers.

The participants (n = 156, median age at the time of inclusion = 20 y.o., median age throughout the study = 22 y.o., 52% female; see Table S1) belonged to three Toba/Qom communities in Formosa, a northern province in Argentina, all of which share a recent historical past and are closely related culturally and ethnically:10

  1. Barrio Toba: a community located adjacent to Ingeniero Juárez (23°47′ S, 61°48′ W), a town with close to 20,000 inhabitants, including the Toba/Qom people. This community has had 24/7 access to electricity at their homes throughout the time we surveyed them. People in this community settled adjacent to Ingeniero Juárez after migrating from the northern Toba/Qom communities in the 1990s.

  2. Vaca Perdida: a rural community of a few hundred members located ~60 km north of Ingeniero Juárez (23°29′ S, 61°38′ W). This community slowly gained access to electricity when the local government started introducing it in 2016 and has since had regular access to electricity and, more recently, paid internet access.

  3. Isla García: this small rural community (50–70 people) is named after the family that founded it. Since we started working with them, the García family and their relatives have moved through three different locations in the vicinity of Vaca Perdida. The first year they were able to consistently access electricity was 2017, and since 2018, they have remained at their current location with full access to electricity. In 2023, the local government installed a public internet antenna on the community grounds.

Consent was obtained orally from every participant after detailing the study procedures in Spanish, with the help of a Toba/Qom translator whenever necessary. For subjects under 18 years of age, consent as well as assent by one of the parents were obtained. Participants were not monetarily compensated for their participation in any of these studies. In appreciation of their participation, they received different gifts, usually clothing and daily-use items. At sign-up, participants completed a demographic questionnaire with date of birth, sex, and fixed scheduled activities such as school, adult education, or work. Gender identity was not elicited as a separate variable: the distinction between sex and gender as separately reportable categories presupposes a framework that is not locally salient, and eliciting it, even if through an interpreter, would have imposed external categories rather than measured a locally meaningful one. The sex term in our models should therefore be understood to index both biological sex and the socially patterned activities that covary with it in these communities; the two cannot be separated in these data, and no gender-based analysis is possible.

It is important to note that these communities engage in very similar social behavior and cultural activities. Children and adolescents attend school during daylight hours, although the time for each age group may differ slightly between communities. Schools are typically located within walking distance of the houses. Most adults are formally unemployed and rely economically on occasional gigs and government assistance. Of note, substance abuse disorders are rare among these communities, in part due to their Anglican religious practices; indeed, selling alcohol is forbidden in Vaca Perdida. These groups are among the most neglected in Argentina. Recently, the National Institute of Statistics and Census (INDEC) classified communities in this area as “very low income/poor.” Socioeconomic differences between families and/or across time are negligible and thus were not included as a covariate in analyses.

Classifying rural and urban communities includes considering factors such as access to electricity, built housing, and market and labor integration. Beyond their history of access to electricity, rural communities are far from health, continued education, and administrative centers, grocery shops, and banks, and usually lack transportation mediums of their own to travel the challenging ~60 km road necessary to access these in Ingeniero Juárez. The “rural” group in this paper is represented by the aggregated data of both Vaca Perdida and Isla García participants. The Toba/Qom community in Barrio Toba will be assigned the “urban” label in this paper, due to being settled within an urbanized area, even if their access to these urban amenities may be limited due to their economic situation.

Device validation participants

The comparison between actigraphic devices was performed in a separate sample of 22 adult volunteers who were not members of the Toba/Qom communities (13 female, mean age = 35, range 21 to 68 y.o.). These data were collected as part of a separate study on the comparison of actigraphic methods. Each volunteer wore an Actiwatch Spectrum Plus and an Axivity AX3 concurrently on the non-dominant wrist for 7 to 14 days during their normal daily routine, yielding 206 paired sleep events.

That study was approved by the Ethics Committee of Universidad Nacional de Quilmes (CEUNQ No. 6/2023), where L.P.C., I.S. and D.A.G. held teaching appointments at the time; approval was sought from this committee because the study predates the establishment of Universidad de San Andrés Ethics Committee. All volunteers gave written informed consent.

METHOD DETAILS

Study design and dataset composition

Data was collected at different points between 2012 and 2014, 2016 and 2018, and 2023 and 2024, and comprises a total of 12,307 recorded sleep events from 156 unique participants (93 from the rural group, 62 from the urban group, and one participant who moved from the urban group to the rural group after 2018). Data was recorded mostly in the late fall or early winter, and in the late spring or early summer, to better capture seasonal differences. Of note, not every participant was included at every time point. Table S1 details the number of participants, their age and sex composition, and the total and the average number of recorded events per participant according to community, year, and season of collection. Figure S1 shows the study trajectories for each single individual in the study, from their first to their last participation.

Electricity and devices access survey

Since efforts to introduce electricity into different rural settlements in the region started circa 2016, participants were asked upon sign-up whether they had received access to the electric grid at their homes. Every participant included in the study from 2017 to 2024 answered affirmatively.

Starting in the 2023 Spring campaign, when signing up, each participant was asked whether they owned a smartphone. We followed up by asking participants if they had access to a TV, a computer or tablet device, and/or an internet connection at home.

Sunlight and light pollution data

Solar data for the Ingeniero Juárez and Vaca Perdida geographical regions was obtained using the suncalc package.39 For every date present in the dataset, we collected the sunrise and sunset times. These were used to calculate the natural photoperiod for each date, as well as for the filtering of nightly sleep episodes that are described below.

Satellite images through the years of the regions of interest superimposed with data from the Day/Night Band of NASA’s Black Marble Group’s Visible Infrared Imaging Radiometer Suite instrument were obtained from Jurij Stare’s Light Trends application (https://lighttrends.lightpollutionmap.info).

Actigraphy and sleep estimation

Sleep was assessed by wrist actigraphy, and participants also completed daily sleep logs throughout each recording session. Participants wore the actigraphs on their non-dominant wrists, and locomotor activity data was always collected in one-minute bins. Individual recordings were visually examined before proceeding to the sleep estimation to confirm data integrity and discard low-quality recordings (e.g. fragmented data, frequent missing periods of data, frequent saturated counts of activity).

Three different wrist-actigraphic devices were used across the years, and consequently, the methods used for sleep detection also differed:

  • 1

    Participants in the studies before 2014 wore the Actiwatch-L device from Mini Mitter (Bend, OR); data was analyzed for sleep detection, aided by sleep logs, using the Sleep Analysis package of the Mini Mitter Software (version 3.4).11

  • 2

    Participants recorded between 2014 and 2018 were equipped with the Actiwatch Spectrum Plus loggers from Respironics (Bend, OR); sleep detection was performed using the Philips Respironics Actiware software V.6.0.9. These data were included in the analysis of a previous study.40

For both the Actiwatch-L and Actiwatch Spectrum Plus loggers, the software scores an epoch as sleep (as opposed to wake) if the comparison of the activity in that epoch vs. the epochs immediately surrounding the epoch falls below a threshold. The “Medium” setting was used, which sets a threshold for wake at 40 activity counts/min. The length of the interval of immobility for the detection of sleep onset was 10 min.

  • 3

    The newly reported data in this paper, collected since 2023, were obtained with the AX3 devices from Axivity (Newcastle upon Tyne, UK), which are equipped with a tri-axial accelerometer.41 Devices were set at 50hz with a range of 8g. Analysis of the actigraphic data was performed using the GGIR software (version 3.2.6) for R.42 Raw accelerometer data was summarized using the default Euclidean Norm Minus One (ENMO) metric. Sleep was detected using the implemented Van Hees algorithm,41 which identifies sustained immobility bouts based on arm angle changes. A minimum of 16 hours of wear-time was required for each day of recording.

Algorithm-based sleep detection from data cohorts 2 and 3 above was contrasted with individual sleep logs non-exhaustively to assess the overall rationality of estimations. Participants from these cohorts displayed a loss of data from a median of 4% (IQR: 0 – 25) of days due to failure to wear the actimeters for at least 16 hours per day. On those days considered in the analyses, the median time of non-wear was 0.32 hours per day (IQR: 0.07 – 0.79).

Cross-device comparison

Our study followed the recommendations of the Society of Behavioral Sleep Medicine (SBSM) on the use of actigraphy for sleep studies43 (e.g., use of validated devices, placement on the non-dominant wrist, multi-day recordings, and artifact removal). However, a specific recommendation that could not be met due to the study’s nature was the use of the same device model throughout the entire monitoring period. To address this, we performed a comparison of the sleep/wake estimations between the Actiwatch Spectrum Plus and the Axivity AX3 devices in the separate sample of 22 volunteers described above. We performed a Bland-Altman analysis44 based on 206 sleep events, selected with the same criteria used for the Toba/Qom data (described in the next subsection). Due to the non-normality of the differences, we performed a non-parametric analysis of the results. Differences between methods, expressed as AX3/GGIR minus Actiwatch/Actiware, were: for sleep start, mean = 16.5 min, median = 1.4 min; for sleep end, mean = −0.9 min, median = 2.6 min; for sleep duration, mean = −17.5 min, median = 1.2 min. See Figure S2 for more details. These differences in the mean estimations from Actiwatch and AX3 devices are of a biologically relevant magnitude. In consequence, they needed to be addressed at the statistical analysis stage. The method used to account for these differences across data collection methods is described further below.

Data treatment and removal of outlier data

The whole 2012–2024 dataset of sleep events was filtered first to keep only nightly sleep events. Every sleep event starting outside the hours from dusk till dawn on any date was removed from the dataset. This step removed only 181 sleep events (1.5% of the total sleep events). While actigraphy performs poorly for the detection of naps, daytime sleep is rare in these Toba/Qom communities. According to the analysis of sleep logs, the proportion of Toba/Qom that reported naps decreased from 75% in 2012 to 39% in 2024. For those who reported napping, the frequency of reported siestas per person was very low and changed from a median of 2.3 naps per week in 2012 to only one nap per week in 2024.

We then proceeded to detect events that could be considered outliers based on their duration (i.e. too long or too short events). This approach allowed us to remove data points that could represent either artifacts of the sleep detection method, or sleep events from “non-representative” nights (e.g., during travel or while sick). The outlier analysis was performed using the median absolute deviation (MAD) method45 through the Routliers package46, setting a threshold of 2.5 MADs (lower limit = 3.49 h, upper limit = 12.44 h). This step removed an additional 506 (4.1%) sleep events from the dataset. The final dataset on which the analysis was performed had a size of 11,620 sleep records.

Sleep regularity index

Sleep regularity was evaluated using the Sleep Regularity Index (SRI). The SRI, proposed by Phillips et al. 2017,30 represents a re-scaled measure of the probability that a participant is at an equal state of wake or sleep at two time points separated by 24 hours on any series of at least seven consecutive night sleep events. A person who always wakes up and falls asleep at the same time every day will score 100, while a person who randomly falls asleep and wakes up at any time of the day will score 0.

First, we detected blocks of consecutive nightly sleep events for each participant. Subjects were considered to be asleep within the recorded sleep bout and awake outside of it. Starting at the time of the first recorded sleep start and in one-minute bins, we calculated the SRI for each individual block. Only scores for sets of at least seven consecutive recordings were included in the analysis. It is important to note that for the SRI analysis, the filter step of outliers based on sleep duration described above was purposefully not performed, as we preferred to prevent the loss of the continuity of recordings, even if that meant including potentially inaccurate records of sleep.

QUANTIFICATION AND STATISTICAL ANALYSIS

Statistical reporting

Unless otherwise stated, n denotes the number of recorded sleep events, and the corresponding number of participants is given alongside. The final analyzed dataset comprised 11,620 sleep events from 156 participants. Model results are reported as fixed-effect coefficients with their standard errors and 95% confidence intervals; post-hoc contrasts are reported as estimated differences with 95% confidence intervals, t values and p values. Significance was defined as p < 0.05. Full fixed-effects outputs for every model reported in this paper are given in Tables 1, 2, S4, and S6.

Linear mixed-effects models

To evaluate the evolution of night sleep timing in these communities, we fit linear models including mixed-effects using the lme4 package.47 Fixed factors included a binary rural/urban factor, a continuous date variable, photoperiod length, type of day (before a workday or before the weekend or a holiday; excluded in the analysis of the SRI), sex, and age. Interaction terms for group and date, as well as sex and age, were also included. Exclusively for the analysis of the sleep start, midpoint, and end variables, the age factor was included as a 3-grade polynomial term to account for the known non-linear progression of sleep timing across adolescence and adulthood.48 Individual identities were included as a random factor to account for repeated measures. Below is an example of a lme4 model equation in this study, illustrating the interaction terms and general syntax:

Sleep_onset~Group*Date+Photoperiod_length+Sex*Age+Week_day+1+Date|participant_ID

Intercept levels for all linear models were set as follows: age = 22 y.o. (the median age across all participants), photoperiod length = 12 h, group = Rural, date = July 1st, 2018 (i.e. 2018.5), sex = female, type of day = free day.

Sex was included as a fixed effect in every model and was allowed to interact with age, so that sex differences in the age trajectory of sleep timing are explicitly modelled. Ancestry and ethnicity are constant across the sample by design, and socioeconomic differences within these communities are negligible and therefore were not modelled.

Post-hoc contrasts

To evaluate post-hoc contrasts between years or groups, we used the emmeans package.49 The tools in this package return the values of the dependent variable predicted by the models based on given independent variables, and test contrasts between levels (p-values are Tukey-corrected for multiple comparisons). A ten-year span between 2014 and 2024 was used when comparing between years to allow for an easier appreciation of the magnitude of the results. We decided to set the modelled age at 22 years old, which is the median age of participants in the dataset across all recording campaigns. It is important to note that age was not included as an interacting factor with groups, years, or photoperiod length in our models, and thus while the absolute predicted values can change according to their dynamics through age, the differences in the contrasts between the aforementioned factors do not change regardless of the set age value. In other words, had we applied our model to a 40-year-old participant, the delays or advances in sleep timing between the two communities would be the same, even when absolute sleep times may be different.

Correction for the change in devices

To account for the systematic differences and variability in sleep estimations between the Actiwatch and AX3 devices, we implemented a non-parametric multiple imputation strategy utilizing empirical error distributions. In essence, we pooled the results of 100 simulated analyses, each of them using a different dataset created to reflect the values that would have likely been obtained had the older Actiwatch data been collected with the newer AX3 devices. This process consisted of the following steps:

Step 1: Derivation of empirical errors. For each sleep parameter, we extracted the exact paired differences observed between the two devices during the simultaneous calibration phase. For each sleep variable, these consisted of 206 observations. Because these differences exhibited a skewed, non-normal distribution, utilizing the raw empirical distribution rather than parametric assumptions provides a more accurate representation of the measurement error.

Step 2: Data imputation via bootstrapping. In the longitudinal dataset, every sleep observation recorded with an Actiwatch device was shifted by adding an error value randomly sampled with replacement from the empirical distribution of differences derived in Step 1.

Step 3: Model fitting and contrasts. The linear mixed-effects models (including all relevant covariates and random effects) were fitted to this newly adjusted dataset. Additionally, estimated marginal means were computed to extract the 2014–2024 post-hoc temporal contrasts.

Step 4: Iteration and pooling. Steps 2 and 3 were then repeated independently 99 times, generating a final set of 100 distinct model estimations. Finally, Rubin’s rules50 were applied to pool the variable coefficients, contrast estimates, standard errors, 95% confidence intervals, and t and p-values, thereby incorporating the uncertainty introduced by the change in instrumentation into the final statistical inference.

A standalone, documented implementation of this imputation and pooling procedure, together with an example of use, is available at https://github.com/LeandroCasiraghi/ax3-aw-multiple-error-imputation.

Model assumptions and software

All the data manipulation and analysis were carried out in the R language51 (version 4.5.1) on the RStudio platform (version 2026.08.2)52 unless specifically indicated. The assumptions of the models were evaluated using the check_model function from the performance package.53 We tested the following features: 1) posterior predictive check (are real values far from modelled values?), 2) linearity, 3) homogeneity of variance, 4) normality of residuals, 5) normality of random effects, and 6) influential observations (i.e., are there too many outliers?).

Supplementary analyses and tests

To discard any spurious or artifactual effects that may inflate or artificially generate the observed results, we conducted a series of additional sensitivity analyses and checks. These consist of modified models for sleep onset and/or sleep duration (what we consider to be the most important variables under study), adding potentially problematic variables or combinations of variables. The models were run exactly as described above for the main results, including the detailed multiple imputation of errors through 100 iterations whenever the data included those from both AX3 and Actiwatch devices.

Test of candidate mechanisms

To evaluate whether access to electricity or to smartphones/internet could account for the longitudinal trend in sleep onset, we re-fit the model including two additional binary covariates: electricity access, and a “smartphone/internet era” indicator. Because both covariates are, by construction, correlated with the continuous date term (i.e. each represents a binary partition of calendar time) we assessed collinearity directly in the fitted model using the adjusted generalized variance-inflation factor (GVIF1/(2·Df))54 with the “vif” function from the car package,55 and the condition number of the scaled fixed-effect design matrix. The “electricity access” variable reported a value of 1.43, the “smartphone/internet era” factor a value of 1.76, both suggesting low collinearity. The date variable returned a value of 2.27, which warrants careful consideration but that is well below the threshold for severe collinearity (i.e. 3.16, equivalent to a traditional VIF value of 10).56 The condition number (12.0) remained below the conventional threshold for even moderate collinearity (<30)57, and the date coefficient remained significant despite an expectedly larger standard error (Table S4).

Older Actiwatch-L devices

Since we no longer have access to the old Actiwatch-L devices, we decided to treat them in the same way as the Actiwatch Plus devices during the multiple imputation of errors phase of our analysis. While these devices shared some similar features, they differed in both manufacturing as well as in the analysis pipeline and software. In order to test whether the inclusion of these older data in the analysis produced inflated results, we ran the analyses on a dataset excluding all data collected with the Actiwatch-L devices until 2014 (a total of 328 sleep events). Even after this reduction in the dataset, all longitudinal effects are preserved in direction, magnitude, and significance (Table S5).

Sleep delays before the COVID-19 pandemic

To assess whether the detected longitudinal delay in sleep start times was driven by an acute change in sleep habits during the COVID-19 pandemic, we fit a model analyzing sleep onsets only in the 2012 to 2018 span, before the pandemic. The model was fitted over 3,975 sleep events from 95 participants, and due to the absence of data from AX3 devices, no multiple imputation was required. The model revealed a delay of sleep onsets that does not differ significantly from that detected for the longitudinal dataset (Table S6).

Short- or long-term participation

This analysis aimed to determine if the longitudinal changes described truly represent a within-person change or merely a between-person turnover. To test this, we classified each participant according to whether they participated in the study for a short time (across less than two years) or for a long time (across two years or more). Then, we re-estimated the changes in sleep with data only from participants in the “long-window” classification. We recovered the full-sample effects: both communities showed significant delays in sleep onset, while only the rural community showed a significant reduction in sleep duration.

Year of entry into the study

This analysis aimed to assess whether the detected longitudinal effects represent an artifact of cohort composition (that is, of different participants entering the study in different years), rather than a genuine within-person change. To do this, we replicated the models including the “entry year” variable for each participant, as well as its interaction term with the longitudinal chronological variable. The longitudinal trend remained robust to these additions (0.12 ± 0.03 h/yr for sleep onset, −0.05 ± 0.02 h/yr for sleep duration).

Supplementary Material

1

Document S1. Figures S1–S2 and Tables S1-S6.

KEY RESOURCES TABLE

REAGENT or RESOURCE SOURCE IDENTIFIER
Deposited data
Longitudinal actigraphic sleep data from Toba/Qom communities 2012–2024 This paper 10.5281/zenodo.22216371
Concurrent Actiwatch Spectrum Plus / Axivity AX3 device comparison dataset + Original code for data analysis including multiple imputation of errors This paper 10.5281/zenodo.22697596
Satellite images including VIIRS Day/Night Band night-time radiance, 2012–2024 NASA Black Marble Group; Radiance Light Trends https://lighttrends.lightpollutionmap.info/
Software and algorithms
R version 4.5.1 R Core Team52 https://cran.r-project.org/bin/windows/base/old/4.5.1/
RStudio version 2026.08.2 Posit team53 https://dailies.rstudio.com/version/2026.08.2+200/
GGIR version 3.2.6 van Hees and Migueles43 https://doi.org/10.32614/CRAN.package.GGIR; https://cran.r-project.org/src/contrib/Archive/GGIR/GGIR_3.2-6.tar.gz
van Hees sleep detection algorithm van Hees et al.42 https://doi.org/10.1038/s41598-018-31266-z
Philips Respironics Actiware version 6.0.9 Philips Respironics https://www.usa.philips.com/healthcare/sites/actigraphy/solutions/actiware; discontinued
Mini Mitter Sleep Analysis software version 3.4 Mini Mitter Co. https://web.archive.org/web/20040220223248/http://minimitter.com/Products/ActiWatch/software.html; discontinued
lme4 Bates et al.48 https://doi.org/10.32614/CRAN.package.lme4
emmeans Lenth50 https://doi.org/10.32614/CRAN.package.emmeans
performance Lüdecke et al.54 https://doi.org/10.32614/CRAN.package.performance
Routliers Delacre and Klein47 https://doi.org/10.32614/CRAN.package.Routliers
suncalc Thieurmel and Elmarhraoui39 https://doi.org/10.32614/CRAN.package.suncalc
car Fox et al.56 https://doi.org/10.32614/CRAN.package.car
Sleep Regularity Index Phillips et al.30 https://doi.org/10.1038/s41598-017-03171-4
Other
Actiwatch-L actigraph Mini Mitter Co. (Bend, OR, USA) https://web.archive.org/web/20070815181903/http://www.minimitter.com/products/actiwatch/awl.html; discontinued
Actiwatch Spectrum Plus actigraph Philips Respironics (Bend, OR, USA) https://web.archive.org/web/20190922023022/http://www.actigraphy.respironics.com/solutions/actiwatch/actiwatch-plus.html; discontinued
Axivity AX3 tri-axial accelerometer Axivity Ltd (Newcastle upon Tyne, UK) https://axivity.com/product/ax3

Highlights:

  1. A 12-year actigraphy study tracked sleep in Toba/Qom rural and urban communities.

  2. Sleep onset was delayed by about an hour in both groups over a decade within this span.

  3. Rural participants lost 35 min of sleep over that decade; urban duration held stable.

  4. These shifts paralleled electrification and later smartphone/internet access.

Acknowledgments

The authors extend their gratitude to the participants from the Toba/Qom communities of Ingeniero Juárez, Vaca Perdida, and Isla García. Their welcoming generosity has been invaluable to our research team. We would like to especially thank community leaders whose assistance on-site has been essential for our work. Additionally, we are grateful to Hotel Indiano in Ingeniero Juárez for their hospitality and support. We are thankful to Professors Claudia Valeggia and Eduardo Fernández-Duque (Yale University), who introduced our group and were key to establishing our ties to these communities. Marcelo Rotundo and Fundación ECO (Formosa) provided invaluable help during field campaigns throughout the years. Finally, the authors thank all the researchers who participated in fieldwork across all these years. This work has been supported by funding from these agencies and institutions: National Institutes of Health, Grant # R01HL162311; Program of Support to Research from Universidad de San Andrés; NSF RAPID, Award #1743364; Leakey Foundation grant #1266.

Footnotes

Publisher's Disclaimer: This is a PDF file of an unedited manuscript that has been accepted for publication. As a service to our customers we are providing this early version of the manuscript. The manuscript will undergo copyediting, typesetting, and review of the resulting proof before it is published in its final form. Please note that during the production process errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain.

Declaration of interests:

The authors declare no competing interests.

References

  • 1.Bin YS, Marshall NS, and Glozier N. (2012). Secular trends in adult sleep duration: a systematic review. Sleep Med. Rev. 16, 223–230. 10.1016/j.smrv.2011.07.003. [DOI] [PubMed] [Google Scholar]
  • 2.Samson DR, and McKinnon L. (2025). Are humans facing a sleep epidemic or enlightenment? Large-scale, industrial societies exhibit long, efficient sleep yet weak circadian function. Proc. R. Soc. B Biol. Sci. 292, 20242319. 10.1098/rspb.2024.2319. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 3.Inc, G. (2024). Americans Sleeping Less, More Stressed. Gallup.com. https://news.gallup.com/poll/642704/americans-sleeping-less-stressed.aspx. [Google Scholar]
  • 4.Tubbs AS, Ghani SB, Valencia D, Jean-Louis G, Killgore WDS, Fernandez F-X, and Grandner MA. (2022). Racial/ethnic minorities have greater declines in sleep duration with higher risk of cardiometabolic disease: An analysis of the U.S. National Health Interview Survey. Sleep Epidemiol. 2, 100022. 10.1016/j.sleepe.2022.100022. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 5.Matricciani L, Olds T, and Petkov J. (2012). In search of lost sleep: Secular trends in the sleep time of school-aged children and adolescents. Sleep Med. Rev. 16, 203–211. 10.1016/j.smrv.2011.03.005. [DOI] [PubMed] [Google Scholar]
  • 6.Bulbulshoev U, Kassam K-A, and Ruelle M. (2011). Ecology of Time: Calendar of the Human Body in the Pamir Mountains. J. Persianate Stud. 4, 146–170. 10.1163/187471611X600369. [DOI] [Google Scholar]
  • 7.Casiraghi LP, and de la Iglesia HO. (2022). Sleep Under Preindustrial Conditions: What We Can Learn from It. In Circadian Regulation Methods in Molecular Biology., Solanas G and Welz P-S, eds. (Springer; US: ), pp. 1–14. 10.1007/978-1-0716-2249-0_1. [DOI] [PubMed] [Google Scholar]
  • 8.Roenneberg T, Kumar CJ, and Merrow M. (2007). The human circadian clock entrains to sun time. Curr. Biol. 17, R44–R45. 10.1016/j.cub.2006.12.011. [DOI] [PubMed] [Google Scholar]
  • 9.Yetish G, and McGregor R. (2019). Hunter-Gatherer Sleep and Novel Human Sleep Adaptations. In Handbook of Behavioral Neuroscience (Elsevier; ), pp. 317–331. 10.1016/B978-0-12-813743-7.00021-9. [DOI] [Google Scholar]
  • 10.Miller ES. (2001). Peoples of the Gran Chaco (Greenwood Publishing Group; ). [Google Scholar]
  • 11.de la Iglesia HO, Fernández-Duque E, Golombek DA, Lanza N, Duffy JF, Czeisler CA, and Valeggia CR. (2015). Access to Electric Light Is Associated with Shorter Sleep Duration in a Traditionally Hunter-Gatherer Community. J. Biol. Rhythms 30, 342–350. 10.1177/0748730415590702. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 12.Beale AD, Pedrazzoli M, Gonçalves B. da S.B., Beijamini F, Duarte NE, Egan KJ, Knutson KL, Schantz M. von, and Roden LC. (2017). Comparison between an African town and a neighbouring village shows delayed, but not decreased, sleep during the early stages of urbanisation. Sci. Rep. 7, 5697. 10.1038/s41598-017-05712-3. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 13.Moreno CRC, Vasconcelos S, Marqueze EC, Lowden A, Middleton B, Fischer FM, Louzada FM, and Skene DJ. (2015). Sleep patterns in Amazon rubber tappers with and without electric light at home. Sci. Rep. 5, 14074. 10.1038/srep14074. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 14.Peixoto CAT, da Silva AGT, Carskadon MA, and Louzada FM. (2009). Adolescents Living in Homes Without Electric Lighting Have Earlier Sleep Times. Behav. Sleep. Med. 7, 73–80. 10.1080/15402000902762311. [DOI] [PubMed] [Google Scholar]
  • 15.Pilz LK, Levandovski R, Oliveira MAB, Hidalgo MP, and Roenneberg T. (2018). Sleep and light exposure across different levels of urbanisation in Brazilian communities. Sci. Rep. 8, 11389. 10.1038/s41598-018-29494-4. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 16.Ekirch AR. (2001). Sleep We Have Lost: Pre-industrial Slumber in the British Isles. Am. Hist. Rev. 106, 343–386. 10.1086/ahr/106.2.343. [DOI] [PubMed] [Google Scholar]
  • 17.Trebucq LL, Moyano MD, Coldeira MF, Golombek DA, De La Iglesia HO, Spiousas I, and Casiraghi LP. (2026). Smartphone use is associated with later and shorter actigraphic sleep, and circadian clock delays in an indigenous Toba/Qom community. SLEEP, zsag039. 10.1093/sleep/zsag039. [DOI] [PubMed] [Google Scholar]
  • 18.Bauducco S, Pillion M, Bartel K, Reynolds C, Kahn M, and Gradisar M. (2024). A bidirectional model of sleep and technology use: A theoretical review of How much, for whom, and which mechanisms. Sleep Med. Rev. 76, 101933. 10.1016/j.smrv.2024.101933. [DOI] [PubMed] [Google Scholar]
  • 19.Tassino B, and Leone MJ. (2025). Night owls of Rio de la Plata region: Real-life scenarios to understand the biological clock. Neuroscience, S0306452225001368. 10.1016/j.neuroscience.2025.02.022. [DOI] [PubMed] [Google Scholar]
  • 20.Anauati MV, Gómez Seeber M, Campanario S, Sosa Escudero W, and Golombek DA. (2024). The economic costs and consequences of (insufficient) sleep: a case study from Latin America. Eur. J. Health Econ. 10.1007/s10198-024-01733-8. [DOI] [PubMed] [Google Scholar]
  • 21.Golombek DA, Eyre H, Spiousas I, Casiraghi LP, Hartikainen KM, Partonen T, Pyykkö M, Reynolds CF, Hynes WM, Bassetti CLA, et al. (2025). Sleep Capital: Linking Brain Health to Wellbeing and Economic Productivity Across the Lifespan. Am. J. Geriatr. Psychiatry Off. J. Am. Assoc. Geriatr. Psychiatry 33, 92–106. 10.1016/j.jagp.2024.07.011. [DOI] [PubMed] [Google Scholar]
  • 22.Hafner M, Stepanek M, Taylor J, Troxel WM, and van Stolk C. (2017). Why Sleep Matters-The Economic Costs of Insufficient Sleep: A Cross-Country Comparative Analysis. Rand Health Q. 6, 11. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 23.Streatfeild J, Smith J, Mansfield D, Pezzullo L, and Hillman D. (2021). The social and economic cost of sleep disorders. Sleep 44, zsab132. 10.1093/sleep/zsab132. [DOI] [PubMed] [Google Scholar]
  • 24.Chaput J-P, Dutil C, Featherstone R, Ross R, Giangregorio L, Saunders TJ, Janssen I, Poitras VJ, Kho ME, Ross-White A, et al. (2020). Sleep duration and health in adults: an overview of systematic reviews. Appl. Physiol. Nutr. Metab. 45, S218–S231. 10.1139/apnm-2020-0034. [DOI] [PubMed] [Google Scholar]
  • 25.Grandner MA. (2017). Sleep, Health, and Society. Sleep Med. Clin. 12, 1–22. 10.1016/j.jsmc.2016.10.012. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 26.Cappuccio FP, D’Elia L, Strazzullo P, and Miller MA. (2010). Sleep duration and all-cause mortality: a systematic review and meta-analysis of prospective studies. Sleep 33, 585–592. 10.1093/sleep/33.5.585. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 27.Hirshkowitz M, Whiton K, Albert SM, Alessi C, Bruni O, DonCarlos L, Hazen N, Herman J, Adams Hillard PJ, Katz ES, et al. (2015). National Sleep Foundation’s updated sleep duration recommendations: final report. Sleep Health 1, 233–243. 10.1016/j.sleh.2015.10.004. [DOI] [PubMed] [Google Scholar]
  • 28.Fischer D, Klerman EB, and Phillips AJK. (2021). Measuring sleep regularity: theoretical properties and practical usage of existing metrics. Sleep 44, zsab103. 10.1093/sleep/zsab103. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 29.Windred DP, Burns AC, Lane JM, Saxena R, Rutter MK, Cain SW, and Phillips AJK. (2024). Sleep regularity is a stronger predictor of mortality risk than sleep duration: A prospective cohort study. SLEEP 47, zsad253. 10.1093/sleep/zsad253. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 30.Phillips AJK, Clerx WM, O’Brien CS, Sano A, Barger LK, Picard RW, Lockley SW, Klerman EB, and Czeisler CA. (2017). Irregular sleep/wake patterns are associated with poorer academic performance and delayed circadian and sleep/wake timing. Sci. Rep. 7, 3216. 10.1038/s41598-017-03171-4. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 31.McKinnon L, Samson DR, Nunn CL, Rowlands A, Salvante KG, and Nepomnaschy PA. (2022). Technological infrastructure, sleep, and rest-activity patterns in a Kaqchikel Maya community. PLOS ONE 17, e0277416. 10.1371/journal.pone.0277416. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 32.Samson DR. (2021). The Human Sleep Paradox: The Unexpected Sleeping Habits of Homo sapiens. Annu. Rev. Anthropol. 50, 259–274. 10.1146/annurev-anthro-010220-075523. [DOI] [Google Scholar]
  • 33.Maczák B, Hordós AZ, and Vadai G. (2025). The Impact of Quantifying Human Locomotor Activity on Examining Sleep–Wake Cycles. Sensors 25, 7659. 10.3390/s25247659. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 34.Leone MJ, Sigman M, and Golombek DA. (2020). Effects of lockdown on human sleep and chronotype during the COVID-19 pandemic. Curr. Biol. 30, R930–R931. 10.1016/j.cub.2020.07.015. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 35.Blume C, Schmidt MH, and Cajochen C. (2020). Effects of the COVID-19 lockdown on human sleep and rest-activity rhythms. Curr. Biol. 30, R795–R797. 10.1016/j.cub.2020.06.021. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 36.Wright KP, Linton SK, Withrow D, Casiraghi L, Lanza SM, Iglesia H. de la, Vetter C, and Depner CM. (2020). Sleep in university students prior to and during COVID-19 Stay-at-Home orders. Curr. Biol. 30, R797–R798. 10.1016/j.cub.2020.06.022. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 37.Rice A, Sather O, Wright KP, Vetter C, Martin MA, and de la Iglesia HO. (2023). COVID-19 stay-at-home restrictions increase the alignment in sleep and light exposure between school days and weekends in university students. Sleep 46, zsad059. 10.1093/sleep/zsad059. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 38.Iparraguirre G. (2022). Cultural Rhythmics: Applied Anthropology and Global Development from Latin America (Emerald Publishing Limited; ) 10.1108/9781803828237. [DOI] [Google Scholar]
  • 39.Thieurmel B, and Elmarhraoui A. (2017). suncalc: Compute Sun Position, Sunlight Phases, Moon Position and Lunar Phase. 10.32614/CRAN.package.suncalc. [DOI] [Google Scholar]
  • 40.Casiraghi LP, Spiousas I, Dunster GP, McGlothlen K, Fernández-Duque E, Valeggia C, and De La Iglesia HO. (2021). Moonstruck sleep: Synchronization of human sleep with the moon cycle under field conditions. Sci. Adv. 7, eabe0465. 10.1126/sciadv.abe0465. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 41.Van Hees VT, Sabia S, Anderson KN, Denton SJ, Oliver J, Catt M, Abell JG, Kivimäki M, Trenell MI, and Singh-Manoux A. (2015). A Novel, Open Access Method to Assess Sleep Duration Using a Wrist-Worn Accelerometer. PLOS ONE 10, e0142533. 10.1371/journal.pone.0142533. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 42.Van Hees VT, and Migueles JH. (2013). GGIR: Raw Accelerometer Data Analysis. 10.32614/CRAN.package. [DOI] [Google Scholar]
  • 43.Ancoli-Israel S, Cole R, Alessi C, Chambers M, Moorcroft W, and Pollak CP. (2003). The role of actigraphy in the study of sleep and circadian rhythms. Sleep 26, 342–392. 10.1093/sleep/26.3.342. [DOI] [PubMed] [Google Scholar]
  • 44.Bland JM, and Altman DG. (1999). Measuring agreement in method comparison studies. Stat. Methods Med. Res. 8, 135–160. 10.1177/096228029900800204. [DOI] [PubMed] [Google Scholar]
  • 45.Leys C, Ley C, Klein O, Bernard P, and Licata L. (2013). Detecting outliers: Do not use standard deviation around the mean, use absolute deviation around the median. J. Exp. Soc. Psychol. 49, 764–766. 10.1016/j.jesp.2013.03.013. [DOI] [Google Scholar]
  • 46.Delacre M, and Klein O. (2019). Routliers: Robust Outliers Detection. 10.32614/CRAN.package.Routliers. [DOI] [Google Scholar]
  • 47.Bates D, Maechler M, Bolker B, Walker S, Christensen RHB, Singmann H, Dai B, Scheipl F, Grothendieck G, Green P, et al. (2025). lme4: Linear Mixed-Effects Models using “Eigen” and S4. Version 1.1–36. [Google Scholar]
  • 48.Roenneberg T, Kuehnle T, Pramstaller PP, Ricken J, Havel M, Guth A, and Merrow M. (2004). A marker for the end of adolescence. Curr. Biol. 14, R1038–R1039. 10.1016/j.cub.2004.11.039. [DOI] [PubMed] [Google Scholar]
  • 49.Lenth RV. (2017). emmeans: Estimated Marginal Means, aka Least-Squares Means. 10.32614/CRAN.package.emmeans. [DOI] [Google Scholar]
  • 50.Rubin DB. (1987). Multiple Imputation for Nonresponse in Surveys 1st ed. (Wiley; ) 10.1002/9780470316696. [DOI] [Google Scholar]
  • 51.R Core Team (2024). R: A Language and Environment for Statistical Computing (R Foundation for Statistical Computing; ). [Google Scholar]
  • 52.Posit team (2024). RStudio: Integrated Development Environment for R (Posit Software, PBC; ). [Google Scholar]
  • 53.Lüdecke D, Makowski D, Ben-Shachar MS, Patil I, Waggoner P, Wiernik BM, and Thériault R. (2019). performance: Assessment of Regression Models Performance. 10.32614/CRAN.package.performance. [DOI] [Google Scholar]
  • 54.Fox J, and Monette G. (1992). Generalized Collinearity Diagnostics. J. Am. Stat. Assoc. 87, 178–183. 10.2307/2290467. [DOI] [Google Scholar]
  • 55.Fox J, Weisberg S, and Price B. (2001). car: Companion to Applied Regression. 10.32614/CRAN.package.car. [DOI] [Google Scholar]
  • 56.Hair J, Anderson R, Babin B, and Black W. (2013). Multivariate Data Analysis (Pearson International; ). [Google Scholar]
  • 57.Belsley DA, Kuh E, and Welsch RE. (1980). Regression Diagnostics: Identifying Influential Data and Sources of Collinearity 1st ed. (Wiley; ) 10.1002/0471725153. [DOI] [Google Scholar]

Associated Data

This section collects any data citations, data availability statements, or supplementary materials included in this article.

Supplementary Materials

1

Data Availability Statement

  • Due to the vulnerable nature of these communities, the longitudinal sleep data from Toba/Qom participants will be available to researchers upon reasonable request to the lead contact. Requests for de-identified, individual-level data should detail the proposed uses of the data and research questions, procedures that will be employed to ensure data security, any potential benefits to the study communities, and precautions for avoiding stigmatizing interpretations of the results, and its use may require a data use agreement.

  • Data used in the cross-correlation of AX3 and Actiwatch devices along with a set of scripts in the R language that allows users to replicate the linear model analysis including the multiple imputation of errors, with the application of Rubin’s Rules, are available at https://github.com/LeandroCasiraghi/ax3-aw-multiple-error-imputation.

Any additional information required to reanalyze the data reported in this paper is available from the lead contact upon proper request.

RESOURCES