Moonstruck sleep: Synchronization of human sleep with the moon ...

Share
RESOURCES
As a library, NLM provides access to scientific literature. Inclusion in an NLM database does not imply endorsement of, or agreement with,
the contents by NLM or the National Institutes of Health.
Learn more:
PMC Disclaimer
|
PMC Copyright Notice
This is an open-access article distributed under the terms of the Creative Commons Attribution-NonCommercial license, which permits use, distribution, and reproduction in any medium, so long as the resultant use is not for commercial advantage and provided the original work is properly cited.
Associated Data
Sleep is delayed and shortened preceding the full moon phase, both in absence of electric light and in modern urban societies.
Abstract
Before the availability of artificial light, moonlight was the only source of light sufficient to stimulate nighttime activity; still, evidence for the modulation of sleep timing by lunar phases is controversial. Here, we use wrist actimetry to show a clear synchronization of nocturnal sleep timing with the lunar cycle in participants living in environments that range from a rural setting with and without access to electricity in indigenous Toba/Qom communities in Argentina to a highly urbanized postindustrial setting in the United States. Our results show that sleep starts later and is shorter on the nights before the full moon when moonlight is available during the hours following dusk. Our data suggest that moonlight likely stimulated nocturnal activity and inhibited sleep in preindustrial communities and that access to artificial light may emulate the ancestral effect of early-night moonlight.
INTRODUCTION
Timing and duration of sleep have changed vastly throughout human evolution and history, following changes in social organization and subsistence. Human beings, with reduced vision capabilities in low-lit environments, are mostly diurnal, and it is believed that nomadic groups timed their sleep onset to the time after dusk when it became too dark to be safe hunting and gathering (1). The establishment of industrial societies, with widespread availability of artificial light sources, allowed humans to accommodate their sleep and wake patterns to modern society demands by creating well-lit—or darkened—environments that isolated them considerably from natural cycles. These artificially lit environments, which can acutely inhibit sleep, also entrain the central body clock in the brain that controls the timing of sleep leading to a delayed onset of sleep and a shorter nocturnal sleep bout (1–4).
While the sun is the most important source of light and synchronizer of circadian rhythms for almost all species, moonlight also modulates nocturnal activity in organisms ranging from invertebrate larvae to primates (5). Moonlight is so bright to the human eye that it is entirely reasonable to imagine that, in the absence of other sources of light, this source of nocturnal light could have had a role in modulating human nocturnal activity and sleep. However, whether the moon cycle can modulate human nocturnal activity and sleep remains a matter of controversy. Some authors have argued against strong effects of moon phase on human behavior and biological rhythms (6–8), but recent studies have reported that human sleep and cortical activity under strictly controlled laboratory conditions are synchronized with lunar phases (9, 10). The controversy generated by these studies has underscored the need for longitudinal studies that can assess the potential effects of moon cycle on sleep (11).
To examine the hypothesis that the lunar phase affects sleep, we conducted a study with three Western Toba/Qom communities of the Argentinian province of Formosa. Once exclusively hunter-gatherers, these geographically spread indigenous communities share a recent historical past and live under very different levels of urbanization (12). We tested the prediction that in communities without access to electricity, moonlit nights would be associated with increased nocturnal activity and decreased sleep. We worked with three Toba/Qom communities (see detailed information in Materials and Methods and table S1): one in an urban setting with full access to electricity (Ur) and two rural communities, one with access to limited electric light (Ru-LL) and another with no access to electric light at all (Ru-NL).
RESULTS
Consistent with previous studies (3), shorter sleep duration and a delayed onset of sleep were associated with increased access to electric light (fig. S1 and table S2). Moreover, both the duration and the time of sleep onset showed a clear modulation throughout the moon cycle that was evident in the whole population, as well as in the individual communities. The peak in sleep onset time and trough of sleep duration took place 3 to 5 days before the night of full moon (Fig. 1). Times of sleep offset suggested a negligible variation across the moon cycle (fig. S2).
(A and B) Double plots of the average duration and onset of sleep in the Toba/Qom population across the moon cycle expressed as average z scores (±SEM; N = 69 participants). Solid lines represent the best fit of the complete dataset to sinusoidal curves with a 30-day period from a nonlinear least squares fit (see Materials and Methods), and the vertical dashed lines indicate the trough of sleep duration (i.e., the shorter sleep events) and the acrophase of sleep onset (i.e., the latest sleeping times). Best-fit equations are indicated for each variable. Fitted sine wave amplitudes mean and 95% confidence intervals (CIs): duration: 0.31 [0.25 to 0.37]; onset: 0.46 [0.40 to 0.51]. (C and D) Double plots of the average values (±SEM) of duration (in minutes) and onset of sleep (in minutes after the astronomical dusk) by community. Solid lines represent best fits for each community data subset. Fitted sine wave amplitudes mean and 95% CIs: duration: Ru-NL, 8.8 [4.9 to 12.8]; Ru-LL, 7.5 [4.0 to 11.0]; Ur, 9.4 [4.6 to 14.2]; onset: Ru-NL, 10.0 [6.2 to 13.7]; Ru-LL, 12.1 [9.1 to 15.1]; Ur, 6.4 [2.5 to 10.3]. Amplitude and phase parameters for all fits are summarized in table S9. Individual data series for participants with records for at least 80% of the moon cycle were filtered through moving average with a window of 7 days before summarizing the data. Number of participants: Ru-LL, 20; Ru-NL, 23; Ur, 26.
The modulation of sleep duration and onset across the moon cycle was evident at the individual level for most participants in the three communities (Fig. 2 and figs. S3 and S4). We fitted individual data to sine waves with a 30-day period through a nonlinear least squares approach and analyzed the parameters of the best-fitting participants (i.e., the three best quartiles according to the standard error of the regression, S; n = 51). The individual phases of sleep duration and onset showed a consistent clustering of the troughs of sleep duration and the peaks of sleep onset on the days before the full moon (Rayleigh z tests, mean phase in days before the full moon [fiducial limits]: duration, 2.8 [5.0 to 2.5], P = 6 × 10−4; onset: 3.3 [4.4 to 2.2], P = 3 × 10−7; fig. S5 and Supplementary Text). Changes in each participant’s sleep duration across the lunar cycle ranged from 20 to more than 90 min and did not differ considerably between groups {mean duration change in minutes [95% confidence interval (CI)]: Ru-NL, 46 [36 to 56]; Ru-LL, 52 [41 to 63]; Ur, 58 [50 to 67]}. Changes in the onset of sleep varied from 30 to 80 min (Ru-NL, 29 [17 to 41]; Ru-LL, 32 [20 to 43]; Ur, 32 [24 to 40]). Thus, moon phases were associated with predictable and biologically relevant changes in daily sleep timing.
Double plots of (A) sleep duration and (B) onset of sleep, expressed as z scores, of the 35 best-fitting participants (as evaluated by the standard error of regressions, S) in the study for each variable. Dots indicate the data on a given night in the cycle, and colored lines represent the best fit of a sine wave with a 30-day period through a nonlinear least squares approach. Individual data series for participants with records for at least 80% of the moon cycle were filtered through moving average with a window of 7 days before summarizing the data. The numbers on the bottom of each plot identify the participant (left) and the S value for the fit (right). The data for the complete set of participants are presented in figs. S3 and S4.
The moon is responsible for several environmental cycles, but its lighting power during the night is arguably the most relevant cycle to humans in natural conditions. Humans typically start their daily sleep bout some hours after dusk but rarely wake up before dawn, a pattern we also documented earlier among the Toba/Qom (3). Moonlight intensity is sufficient to allow outdoor activities, and it is likely to prevent sleep initiation; in contrast, and importantly, it is unlikely to wake somebody who has already fallen asleep. In this context, it is primarily moonlight available during the first hours of the night that is more likely to drive changes in the onset of sleep. In contrast, moonlight late in the night, when most individuals are typically asleep, should have little influence on sleep onset or duration. The hours of available moonlight change predictably through the moon cycle according to the time of moonrise, approximately 50 min later every day (fig. S6); under a mostly symmetrical photoperiod like that of the spring season, moonlight becomes less available during the early night on the nights that follow the full moon night. We hypothesized that the adaptive value of the synchronization between sleep and the moon cycles is to stimulate wakefulness on nights when moonlight is available during the early night, which are the nights that precede the night of full moon but not the ones that follow it, when the moon rises much later than dusk.
To test predictions derived from this hypothesis, we determined the availability of moonlight during the first 6 hours of every night recorded in the study and classified them into three categories of “moonlight phases” (see Materials and Methods and fig. S6): full moonlight (F-ML), no moonlight (No-ML), and waning/waxing moonlight (W-ML). Figure 3 shows the change in sleep timing between F-ML and No-ML for the different communities. Within-subject averaged data comparing the two phases are shown in Fig. 3 and the summarized data are presented in table S3. Linear mixed-effects models (LMEMs) considering community, moonlight phase and the interaction between these, age and sex as fixed effects, and subject identity as a random effect were fit to analyze the associations between changes in sleep variables and moonlight phase [the full description of the estimates and effects sizes are presented in table S4, and type III analysis of variance (ANOVA) evaluation of fixed factors is shown in table S5; see Supplementary Materials for a description of the model and R syntax]. According to the models, the duration of sleep for the three communities was 25 [95% CI: 13 to 37], 19 [7 to 32], and 11 [0 to 21] min longer during No-ML than F-ML nights for the Ru-NL, Ru-LL, and Ur groups, respectively. Accordingly, participants fell asleep 22 [13 to 30], 22 [12 to 31], and 9 [2 to 17] min earlier on No-ML nights for the Ru-NL, Ru-LL, and Ur groups, respectively (table S6). Time of sleep offset did not show consistent differences between phases. Contrasts suggested a differential effect of the moonlight availability on the sleep variables, with larger effects for groups with less access to electric light (table S6). It is worth noting that even a change of about half hour in the duration and time of sleep, which was the lesser change across the cycle found for individuals, is biologically relevant. Both field studies using daily sleep logs during 6 weeks (13) and a retrospective analysis of sleep recordings from a sleep laboratory (9) found that people slept about 20 min less on full moon than on new moon nights, an effect size similar to the one found in our study. Earlier sleep offsets did not coincide with nights in which moonlight was available during the late night, before dawn (fig. S2 and tables S7 and S8), suggesting that late-night moonlight did not drive changes in sleep offset.
Analysis of changes in sleep variables according to the No-ML (no moonlight) and F-ML (full moonlight) phases.
(A) Duration of sleep (hours), (B) sleep onset (hours after dusk), and (C) sleep offset (hours after dawn). Smaller dots represent each individual’s mean, while the bigger dots represent the means for each community and are connected by color lines. Vertical black error bars represent the 95% CI of the mean. The reported P values correspond to the main effect of moonlight for the type III analysis of variance (ANOVA); Cohen’s d values are for Tukey contrasts between moonlight phases in each community in the linear mixed-effects models. Only participants with at least four nights recorded in each phase were considered. Number of participants: Ru-NL, 25; Ru-LL, 25; and Ur, 32.
Beyond moonlight phases within the ~29.5-day synodic month, the gravitational pull of the moon on the Earth surface is maximal twice as frequently (every ~14.75 days), during full or new moons. Recent evidence indicates that bipolar patients can show mood (14, 15) and sleep [one reported case (16)] changes synchronized with these ~15-day semilunar phases. If a similar synchronization pattern exists with the sleep parameters of our Toba/Qom participants, we would expect that some individuals have bimodal peaks or the opposite phase relationship in their sleep parameters. The distribution of sleep duration and onset in the Toba/Qom communities (Fig. 1) and the inspection of the individual recordings (Fig. 2 and figs. S3 and S4) suggest that the two parameters show a second trough and peak, respectively, within the lunar month. To further examine this possibility, we fitted sleep duration and sleep onset data with a periodic function with two sinusoidal components (30- and 15-day periods; see Materials and Methods). This approach led to more informative models than those of single 30-day sine waves for 68 of 69 and 62 of 66 participants for duration and onset, respectively [figs. S7 to S9; differences between Akaike information criterion (AIC) values of ΔAICmedian = −61.6 for duration and ΔAICmedian = −49.2 for onset]. This result indicates that, for the vast majority of participants, sleep timing is associated with both the ~30-day (lunar) and ~ 15-day (semilunar) lunar phases. While the 15-day period component was evident for single participants in every community, we were not able to detect it at the population level in the urban community for neither sleep duration nor sleep onset patterns (table S10). This result indicates that the 15-day component in the urban community is weaker than in the other communities.
All three communities of Toba/Qom, including those in the urbanized setting, showed a strong association of sleep timing with the moon cycle. To explore whether a similar modulation of sleep across the moon cycle occurs in people living in large modern urban environments, we analyzed sleep recordings obtained from 464 University of Washington undergraduate students. Unexpectedly, the changes in sleep duration and onset throughout the moon cycle resembled those of the Toba/Qom people, with sleep events starting later and becoming shorter in the week before the full moon (Fig. 4). Even with the limitations inherent to nonlongitudinal studies, the data suggest that sleep changes across the moon cycle may still be present in completely urbanized environments, where individuals may have little awareness of the synodic moon phases. Light pollution measurements in the highly urbanized areas of Seattle where students typically live reveal values that are above our full moonlight measurements in the Toba/Qom rural environments. Our results are also in line with two retrospective analyses of electroencephalographic recordings of sleep, obtained in the controlled conditions of sleep laboratories, which found that both polysomnographic sleep and cortical activity were synchronized with lunar phase (8, 9); these results point to the importance of longitudinal studies to determine the extent to which the modulation of sleep by the moon cycle prevails under modern living conditions. Studies of this nature may also reveal a semilunar component as the one we report for the Toba/Qom.
Association of sleep duration and onset with the moon cycle in a highly urban setting.
Double plots of (A) sleep duration and (B) sleep onset expressed as z scores (±SEM) on weeknights recorded on 463 college students in different quarters from 2015 to 2018. The differences between individual data points and the mean values in each season were calculated for over ~4300 sleep events. The solid lines represent the best fits to sine waves with a 30-day period to the data from nonlinear least squares fits (see Materials and Methods); wave equations are printed at the bottom left. The dashed vertical lines indicate the phase of shorter (A) and later (B) sleep events. Fitted sine wave amplitudes and 95% CI: duration, 0.34 [0.13 to 0.55]; onset, 0.32 [0.10 to 0.53]. The average data summaries were filtered through a moving average with a window of 7 days. Participants/sleep events per quarter: spring, 173/1729; summer, 66/619; fall, 136/1240; and winter, 88/796.
DISCUSSION
Our results show that sleep timing is synchronized with the moon cycle under a range of living environments. Toba/Qom participants slept less and stayed up later on the days previous to full moon nights, when moonlight is available during the early night. This pattern could represent a response to the availability of moonlight during the first half of the night for communities with limited or no access to electric light. The amplitude of the lunar phase effect on sleep parameters appears to be stronger the more limited the access to electric light is. However, we were able to corroborate this modulation both in a Toba/Qom community living with full access to electricity and in a sample of college students living in a modern city. Together, these results strongly suggest that human sleep is synchronized with lunar phases regardless of ethnic and sociocultural background, and of the level of urbanization.
Increased level of access to artificial light in the Toba/Qom communities correlated with later sleep onsets and shorter duration of sleep. These findings are consistent with previous work from several laboratories including ours (3, 17–19), as well as with the hypothesis that electric light allowed humans to extend their evening activity and push sleep times later into the night, therefore reducing the total amount of night sleep (20). The availability of electric light during the evening mimics the sleep-inhibiting effects of moonlight. This is particularly evident during nights with high moonlight availability during the early hours of the night, in which the timing of sleep is most similar across the Toba/Qom communities with and without access to electric light (Fig. 3). This finding may indicate that the effect that electric light has—delaying sleep onset and shortening sleep—could be emulating an ancestral effect of moonlit evenings, although the light intensities we are typically exposed to in our artificially lit environments are much higher.
While some studies have found minimal, or no association, between the moon cycle and sleep parameters (6, 8, 21, 22), they compared sleep during nights around the full moon to sleep during the nights of both the waxing/waning phases and the new moon, which may not correspond to the peaks and troughs of sleep duration and onset. Inspection of the college students’ data according to full and new moon phases, as these previous studies did, shows no clear association with sleep duration or onset (Fig. 4). In contrast to these studies and in line with our findings, two other studies found an association between sleep parameters and moon phases. Röösli et al. found an effect of moon phase on subjective sleep duration as measured by sleep diaries with shorter sleep durations around full moon nights. Similarly, in a retrospective analysis of polysomnographic sleep in a sleep laboratory, Cajochen et al. (8) found a pattern of sleep duration through moon phases similar to what we found for our population, as well as on the percentage of time spent on specific sleep stages. Last, a recent study reported that nocturnal activity appears to be higher during new moon nights in a group of hunter-gatherer Hadza, but this result is based on a small sample, a sampling window of less than a full moon cycle, and is confounded by the fact that the group has rituals during nights without moonlight (23).
What could be the potential adaptive value of increased activity during moonlit nights? Our interviews with Toba/Qom individuals indicate that moonlit nights are particularly rich in social activities. Toba/Qom elders report that, at times when food was obtained from the forest, moonlit nights had particularly high hunting and fishing activity. Furthermore, mythological stories associate the moon with the female reproductive cycle and sexual relations. The moon in the Toba culture is represented as a man who has sexual relations with women, it induces the first menstruation and regulates the timing of the following menstruations (24). Interestingly, stories told by elder Toba/Qom point to moonlit nights as nights of higher sexual activity. These latter stories suggest the possibility that ancestrally moonlight-associated encounters could have synchronized reproductive activity with women’s fertility (25). Although the true adaptive value of human activity during moonlit nights remains to be determined, our data seem to show that humans—in a variety of environments—are more active and sleep less when moonlight is available during the early hours of the night. This finding, in turn, suggests that the effect of electric light on modern humans may have tapped into an ancestral regulatory role of moonlight on sleep.
Although our results point to moonlight availability during the early night as a likely determinant of later sleep onset and shorter sleep duration, the presence of similar lunar rhythms in sleep parameters in Seattle college students who may not be aware of the availability of moonlight, together with the presence of semilunar (~15-day long) components on the sleep parameters of the Toba/Qom communities, suggests that other physical phenomena associated with the moon cycle could influence sleep. It is thus conceivable that although the ultimate cause—which confers adaptive value—for nocturnal activity in synchrony with the moon cycle is to display activity during moonlit nights, the proximal cause—which induces changes in sleep parameters—for sleep modulation by the moon cycle is the gravitational pull by the moon, which is a more reliable indicator of moon phase than its associated nocturnal illuminance. This hypothesis predicts that both semilunar oscillations— respectively associated with new and full moons—should have the same amplitude. However, although the new and full moon gravitational pulls may be indistinguishable, they consistently occur at different times of the solar day; namely, full moons but not new moons exert their gravitational pull during the night. This raises the possibility that moon gravity could have a time-of-day specific effect on sleep.
A limitation of our observational study is that we cannot establish causality. It would be difficult to manipulate human exposure to the light the moon reflects and virtually impossible to manipulate the exposure to the gravitational pull it exerts on Earth. Nevertheless, it is hard to conceive that the conserved synchronization between sleep and the moon cycle that we report occurred by chance.
MATERIALS AND METHODS
Participants and study groups
All described study procedures were approved by the Internal Review Board of University of Washington’s Human Subjects Division and were in agreement with the Declaration of Helsinki. Oral consent was obtained from every participant from the Toba/Qom communities after a verbal explanation of all procedures in Spanish. All participants were bilingual (Toba/Qom/Spanish). Parental oral consent was obtained for participants under 18 years old, who also gave their assent to participate. University of Washington participants (undergraduate students) were all over 18 years of age and provided written consent.
Toba/Qom participants were aware that we were interested in studying the relationship between the moon cycle and sleep but were unaware of any of our specific predictions, e.g., that moonlit nights would be associated with less sleep. University of Washington participants were unaware of any relationship between their sleep study and moon phases.
Toba/Qom participants
Toba/Qom participants [N = 98, females 56%, mean age (range) = 24.1 (12 to 75)] lived in one of three Toba/Qom communities in the Formosa province, north of Argentina. Each community had different levels of access to electric light:
1)
A community located in the outskirts of Ingeniero Juárez (23°47′ S, 61°48′ W), a town with 19,000 inhabitants. All the participants in this community had 24-hour access to electric light at home and outside through streetlights, as well as to other urban features (consolidated roads, public leisure spaces, and commercial premises). This community is referred to as the Urban group (Ur, n = 40).
2)
A community located in Vaca Perdida (23°29′ S, 61°38′ W), a small rural settlement of approximately 300 people 50 km north of Ingeniero Juárez. These participants had 24-hour access to electricity at their homes, and light fixtures were limited at most to one light bulb per room. In contrast to the Urban community, this community had no electric light poles in outside areas during dark hours. This community is referred to as the Rural, limited light group (Ru-LL, n = 33).
3)
A group of participants living in sparsely distributed houses in a region known as Isla García, approximately 3 km away from Vaca Perdida. These participants lived in small extended-family groups, without any organized settlement features and no access to electric light. While children and some adults may have been exposed to artificial light at school or other settings away from their home location, this could have only occurred during natural daylight times. This community is referred to as the Rural, no light group (Ru-NL, n = 25).
The three communities share the same ethnic and historical past (12). The community at Ingeniero Juárez originated from a group of Toba people who migrated from the northern region in the 1990s. Housing, daily chores, and social behavior are very similar between the communities, and the vast majority of adults are typically unemployed and rely on government subsidies.
Age and sex characteristics of the three study communities are presented in table S1. The three study communities were similar with regard to sex [χ2(2) = 0.686, P = 0.710] and age [ANOVA F(2,96) = 0.774, P = 0.464]. As expected, women displayed onset times that were 15 [95% CI: 0 to 30] min earlier than men (Cohen’s d = 0.432) (26). Age was associated with shorter sleep duration (d = −1.411) and later sleep onset times (d = 0.429), which would be expected for the age range of our participants (table S2) (27, 28).
Data were recorded during field campaigns in three consecutive years: September to October 2016 and 2017 and October to November 2018. The three campaigns were carried out during the spring (September to November in the Southern Hemisphere) to keep environmental factors (including weather, sunrise and sunset times, and sunlight intensity) as stable as possible.
University of Washington participants
We also analyzed sleep data from 464 college students [mean age (range) = 21.5 (18 to 38), 62% females; 41% Caucasian, 26% Asian, 8% Hispanic, and 15% others/nonidentifying] at the University of Washington (Seattle, WA) that were recorded across different quarters between 2015 and 2018 as part of a separate study. The data from one subject who did not wear the watch consistently were discarded.
Locomotor activity and sleep recording
Participants were equipped with Actiwatch Spectrum Plus wrist locomotor-activity loggers (Respironics, OR) for 1 to 2 months in the case of the Toba/Qom participants and from 1 to 3 weeks for the college students. The data acquisition interval was set to 1 min. Recorded data were downloaded and exported using the Philips Respironics Actiware software V.6.0.9. Participants also completed a sleep log throughout their participation, indicating times and locations of sleep events (including naps), and whether they left their home community on any given day. These logs were used for the validation of the data and to discard data from time ranges when the subjects were under different conditions to those in their main study group. The Actiware software was set to determine sleep onsets whenever 10 consecutive minute bins were classified as of immobility (<4 activity counts) and sleep offsets after the last 10 consecutive minutes of immobility in a sleep bout. Minute bins were labeled as “wake” when they surpassed the 40 activity counts.
Data treatment and analysis
Preparation and treatment of raw recordings, as well as the statistical analysis and plotting of data, were performed in R (29) unless otherwise indicated.
Astronomical data
Sun and moon data for the Ingeniero Juárez region through the dates of recording were obtained from NASA’s Jet Propulsion Laboratory HORIZONS Web-Interface (https://ssd.jpl.nasa.gov/horizons.cgi) with a 1-min precision. The synodic moon cycle that determines the well-known “full,” “new,” and “waning” or “waxing” moon phases has an approximate average duration of 29.5 days. All nights in our study were numbered according to their position in the synodic cycle: Full moon nights were considered as “night zero” with all previous and subsequent nights ordered from 1 to 29 or 30, depending on the day of the next full moon, when nights numbering was restarted. In short, a night in position “1” in the cycle represents the night immediately after a full moon.
We use a second classification of phases through the moon cycle according to the availability of moonlight during the early night hours, which we refer to as “moonlight phases.” We calculated the time in the first 6 hours after astronomical dusk—the time that marks the disappearance of natural daylight—during which the moon was above the horizon for every night of recording, rounded this time to the nearest hour, and classified nights into three moonlight phases: (i) nights in which there was no moonlight through the first 6 hours of the evening (no moonlight, No-ML); (ii) nights during which there was moonlight available throughout the whole 6 hours (full moonlight, F-ML); and (iii) all other nights in which the moon either rose after astronomical dusk or set before 6 hours after astronomical dusk (waning/waxing moonlight, W-ML). The distribution of these moonlight phases shows correspondence with the standard moon phases, with the F-ML nights preceding the full moon by 1 week (fig. S6).
Sleep measurements
Actograms from individual recordings were visually inspected to check for data integrity before any analysis. Low-quality criteria included fragmented recording or very frequent missing-data periods, sustained saturated Actiwatch counts, or total lack of activity counts during extended periods of time. Seven Toba/Qom participants who displayed very low quality recordings, due to noncompliance wearing the watch or to Actiwatch malfunction, were removed from the sample at this point.
To detect and discard very long and very short sleep events (artifacts from Actiwatch recording or the software algorithm), as well as irregular sleep events (long daytime siestas, naps, all-nighters), we analyzed our sleep databases for outliers in two steps: first, looking at values of duration, and then at times of sleep onset according to the time of dusk. We performed the analyses through the median absolute deviation (MAD) method setting a threshold of three MADs, using the Routliers package for R (30).
To explore the patterns of sleep duration and timing through the moon cycle in the Toba/Qom communities, we selected the participants with data for at least 80% of the nights in the cycle. Sixty-nine participants met this condition: 20 in Ru-LL, 23 in Ru-NL, and 26 in Ur. Each participant’s data points were first averaged by night in the cycle (for those with data for over a whole cycle), and data for any missing nights were completed by linear interpolation. Data were then run through a moving-average filter with a seven-night window to eliminate low-frequency noise. The data were fitted through a nonlinear least squares approach to the best possible sine curve with a 30-day period using the nls tool from the stats package (29). Because recent literature has pointed to the possibility of semilunar month (~15 days) rhythms in mood and sleep (14–16), alternative fits were performed considering the combination of a sine wave with a 30-day period along a 15-day period wave. The equations for these curves were
y=A30sin[(2π/30)*(night–P30)]+M
(1)
y=A30sin[(2π/30)*(night–P30)]+A15sin[(2π/15)*(night–P15)]+M
(2)
where night represents the night order in the moon cycle, A30 and A15 represent the sine amplitudes, P30 and P15 represent the phase angles of the 30- and 15-day components, respectively, and M represents the center value (which was not included in fits of z-scored data). For the evaluation of population data fits, 95% CIs of the estimated parameters were calculated from complete datasets to determine nonzero amplitude values (i.e., a nonnegligible sinusoidal pattern).
Goodness of fit of individual data fits was measured by the standard error of the regression (S), a proper measure for nonlinear equation fits. When analyzing group differences in parameters obtained from individual fits, we only included subjects within the best three quantiles of S values for their z-score fits for each variable (that is, 51 of the abovementioned 69 participants with enough data). For comparing modeling alternatives, we calculated the AIC for the individual fits and considered a more informative model when it was reduced by more than 2 units. For sleep duration patterns, in 68 of 69 participants, the two-component models were more informative, while for the sleep onset patterns, that was the case for 62 of 66 individuals. The single-component model did not converge in a fit for the sleep onsets of one participant, while the same happened with the two-component model for the sleep onsets of two participants in the sample.
Rayleigh z tests were used to analyze phase clustering and were performed and plotted with El Temps software v1.311 (University of Barcelona, Spain). Fiducial limits were used as a measure of CIs of the mean phases, due to the nonnormal nature of the circular distribution. Circular distributions were compared by Watson-Wheeler tests for homogeneity of phases using the circular package for R (31).
Because we did not count with longitudinal recordings for college students throughout the moon phases, these data were analyzed as follows. We first calculated the average and SD of sleep duration and sleep onset on school nights (Sunday to Monday, excluding nights before a holiday) for each season. Then, we normalized each sleep variable by subtracting the season average and dividing it by the season SD. We then averaged these population data according to the day on the moon cycle and smoothed it by running a moving average with a seven-night window. The data were then fitted to a sine wave with a 30-day period through the nonlinear least squares approach described above.
Statistical analysis
To estimate the associations between the moonlight phases and demographic variables on sleep features, we applied LMEMs using the lme4 package (32). Only participants who presented data for at least four nights within each moonlight phase were considered for these analyses. QQ plots for every LMEM fit were visually inspected to check normality of residues; these are presented in the “Statistical model definitions and diagnostics” section in the Supplementary Text. All the random effects (intercepts by subject) of the best-fitted models were normally distributed. To estimate the predicted difference of sleep duration, onset, and offset with moonlight phase, we calculated Tukey contrasts using the estimated marginal means of each group [via the emmeans R library (33)]. The syntax used for LMEMs is also presented in the Supplementary Text.
Cohen’s d effect size for each fixed effect of the best-fitted models and the above-defined contrasts was calculated using the lme.dscore function of the EMATools R package (34). This package calculates the Cohen’s d using the following equation
d=2β/Std.Error(β)df
The degrees of freedom of the model (df) were calculated using the Satterthwaite approximation [via the lmerTest R library (35)]. CIs for the fixed effect estimations were calculated using the confint function of the stats R library with alpha level set at 5% and using the Wald method. For the CIs of the effect size, we followed the lme.dscore methodology but using the beta estimates limits of the CI instead of the beta estimate itself.
Acknowledgments
We are thankful to A. García, M. Pérez, K. Ortiz, M. Rotundo, E. R. Guiñazú, and N. West for field assistance and the Toba/Qom people for the tremendous support and collaboration. We thank B. Schwartz and C. Helfrich Foster and her team for comments on the manuscript. Funding: This study was supported by NSF RAPID award #1743364 to E.F.-D., C.V., and H.O.d.l.I. and by Leakey Foundation grant 1266 to H.O.d.l.I. Author contributions: L.C. designed the study, collected Toba/Qom data, analyzed data, and wrote the manuscript. I.S. analyzed data and wrote the manuscript. G.P.D. collected Toba/Qom and UW students’ data and revised the manuscript. K.M. organized and processed Toba/Qom recordings and revised the manuscript. E.F.-D. and C.V. provided resources for the study and wrote the manuscript. H.O.d.l.I designed the study, provided resources for the study, collected Toba/Qom and UW students’ data, analyzed data, and wrote the manuscript. Competing interests: The authors declare that they have no competing interests. Data and materials availability: All data needed to evaluate the conclusions in the paper are present in the paper and/or the Supplementary Materials. Additional data related to this paper may be requested from the authors or accessed at https://osf.io/nxtvk/?view_only=6baffd62ae8a485c9e64cd777f9a128c.