Play Behavior Varies with Age, Sex, and Socioecological Context in Wild, Immature Orangutans (Pongo spp.)

Play is thought to serve different purposes at different times during ontogeny. The nature and frequency of play are expected to change accordingly over the developmental trajectory and with socio-ecological context. Orangutans offer the opportunity to disentangle the ontogenetic trajectories of solitary and social play with their extended immature phase, and socio-ecological variation among populations and species. We evaluated the frequency of play in 39 immature individuals across two populations (Pongo pygmaeus wurmbii, at Tuanan, Borneo, and P. abelii at Suaq, Sumatra), age (0–11 years), sex, and social context, using more than 11 500 h of full-day focal observation data. We found independent age trajectories of different play types, with solitary object and solitary locomotor peaking before social play. Social play partners changed during ontogeny, and male immatures were more likely to engage in non-mother social play than females. Overall, social play was more frequent at Suaq than Tuanan, linked to the more frequent availability of partners. Furthermore, per time in association with conspecifics, Tuanan immatures were as likely to engage in social play as their peers at Suaq, suggesting similar intrinsic motivation. Increasing fruit availability correlated with both longer associations and increased social play frequency in the less sociable population of Tuanan, but not at Suaq. Our findings on orangutans support evidence from other species that different play types follow different developmental trajectories, vary with sex, social opportunities, and ecological context. Although drawing functional inferences is challenging, the distinct developmental trajectories reflecting adult sociability and behavioral repertoires may indicate that play serves several, non-mutually exclusive functions during ontogeny. Supplementary Information The online version contains supplementary material available at 10.1007/s10764-023-00414-2.


Introduction
Play can be observed in most vertebrates, and in some invertebrates (Burghardt, 2005(Burghardt, , 2015;;Fagen, 1981;Zylinski, 2015).Because of its omnipresence and the costs that come with it, play is expected to have an adaptive value, rather than being a byproduct or even a maladaptation (Bekoff & Byers, 1981).However, defining play, let alone evaluating its fitness consequences, remains notoriously difficult (Miller, 2017).We use the definition that "play is repeated, seemingly non-functional behavior differing from more adaptive versions structurally, contextually, or developmentally, and initiated when the animal is in a relaxed, unstimulating, or low stress setting" (p.91 Burghardt, 2014).Play is commonly divided into three types: social, solitary object, and locomotor-rotational play (Bekoff & Byers, 1981), and these three play types can be combined within one single bout.Being a heterogenous phenomenon, it is debated whether play types have evolved independently and serve different functions (Burghardt, 1998;Graham & Burghardt, 2010).Empirical studies supporting the numerous theories about the evolution of play behavior remain scarce (Sharpe, 2019).
Supporting the argument that play serves multiple functions, play behavior varies with age, sex, and environmental and social context between and within species 1 3 Play Behavior Varies with Age, Sex, and Socioecological Context… and may reflect context-and development-dependent functions and constraints.First, the three play types show different age trajectories with consecutive peaks of solitary locomotor, solitary object, and social play (chimpanzees [Pan troglodytes]: Cordoni & Palagi, 2011; vervet monkeys [Chlorocebus pygerythrus]: Fairbanks, 2000; gorillas [Gorilla gorilla]: Maestripieri & Ross, 2004;meerkats [Suricata suricatta]: Sharpe, 2005).Second, sex differences reflect adult sociality and skill repertoires, in that male immatures exhibit higher social play frequencies, particularly rough-and-tumble play, compared to female immatures (humans [Homo sapiens]: Lew-Levy et al., 2022;primates: Lonsdorf, 2017;mammals: Marley et al., 2022).Further, immature chimpanzee males have more social partners (Lonsdorf et al., 2014a, b) and show different types of object play compared to females (Kahlenberg & Wrangham, 2010).Third, elevated rates and more sophisticated forms of object play and manipulation are associated with the occurrence of habitual tool use (e.g., chimpanzees versus bonobos: Koops et al., 2015a, b;Myowa-Yamakoshi & Yamakoshi, 2011) and with specific (complex) foraging niches (Bock & Johnson, 2004;O'Hara & Auersperg, 2017).Fourth, play is a costly behavior, and among the first activities to cease when an individual is exposed to stressors (Held & Špinka, 2011;Martin & Caro, 1985).In wild populations, social play frequency varies accordingly with food abundance in some species (Burghardt, 2005; e.g., mantled howler monkeys [Alouatta palliata]: del Toro et al., 2019;squirrel monkeys [Saimiri sciureus]: Stone, 2008).Finally, the frequency and particularly the quality of social play interactions among immature individuals have been associated with more socially complex societies, and adult social and/or fighting competence (e.g., domestic pigs [Sus scrofa domesticus]: Turner et al., 2020).Where this is important, mothers may modify their own behavior to facilitate the socialization of their offspring by associating with suitable partners.This is most apparent in species with high fission-fusion dynamics.Chimpanzee mothers with juvenile offspring, at the age when most play occurs, have been reported to associate more compared to mothers with infants (Williams et al., 2002), and mothers with sons provided them with more play opportunities compared to mothers with daughters (Murray et al., 2014), reflecting male immatures' need to set up social bonds in the light of male-male coalitions for territory defense and hunting in a male philopatric society.In sum, while inferring functions to play behavior in empirical studies may be challenging, play behavior and its different types vary with individual and socio-ecological context, and reflect differences in adult behavioral repertoires in various taxa.
Play behavior of orangutans has been studied predominantly in captivity (Jantschke, 1971;Zucker & Thibaut, 1995), and likely differs from behavior in the wild due to (i) the housing in permanent groups, which offers more opportunities for social play, and (ii) the safety of the captive environment, resulting in more terrestrial activity (e.g., Fröhlich et al., 2021).In wild populations, unweaned offspring spent about 15-45% of their time in solitary play but only 0.5-1.5% in social play on average (van Noordwijk et al., 2009).Reflecting their natural lifestyle, orangutans in zoos also engage more in solitary play than the more gregarious great ape species, such as chimpanzees (Miller & Nadler, 1981).Despite lower interaction rates, orangutans show high social competence in the context-dependent use of play faces (Davila Ross et al., 2008;Waller et al., 2015) and play-solicitation gestures (Cartmill & Byrne, 2010;Fröhlich et al., 2021).Furthermore, like other great apes (Cordoni & Palagi, 2011), Sumatran orangutans (Pongo abelii) show different ontogenetic trajectories for social and solitary play (van Adrichem et al., 2006) in one very gregarious population (Mitra Setia et al., 2009;van Schaik, 1999).
Due to their semi-solitary lifestyle, dependent immature orangutans only have the opportunity for social play with individuals other than the mother when the mother associates with potential play partners.These occasions are rare, however, particularly in Bornean populations where females with unweaned offspring associate with others less than 10% of their active time (Kunz et al., 2021;Mitra Setia et al., 2009;van Noordwijk et al., 2012).The orangutans' semi-solitariness is likely caused by limited food availability, as higher orangutan density, sociability and social tolerance in North-West Sumatran compared to North-East Sumatran and Bornean populations (Husson et al., 2009;Roth et al., 2020;Schuppli et al., 2017;van Schaik, 1999) align with higher forest productivity in these regions (Marshall et al., 2009;Wich et al., 2011).Yet, providing opportunities for offspring socialization was proposed to be the main benefit of female-female associations in Bornean orangutans (Pongo pygmaeus) (van Noordwijk et al., 2012).Within populations, mothers spend less time in association with each other during lean food periods which likely affects their offsprings' play opportunities (Ashbury et al., 2022).Moreover, orangutan females are philopatric (Arora et al., 2012) and preferably associate with their maternal relatives in a Bornean population (van Noordwijk et al., 2012), which may imply an even more limited number of potential play partners.
We aimed to evaluate how play behavior relates to individual and socio-ecological factors in two populations of wild orangutans, one at Suaq, South Aceh, North-West Sumatra (Pongo abelii) and one at Tuanan, Central Kalimantan, Borneo (P.pygmaeus wurmbii).We assessed the three play types (social, solitary object and solitary locomotor play) separately, and evaluated variation with age, sex, study site, local fruit availability and social context.If the three different play types serve different developmental functions or reflect distinct levels of skill acquisition, we predicted distinct ontogenetic trajectories (age) (P1) and no compensation for a lack of social play by solitary play (P2).Moreover, if energetic costs limit play, we predicted its frequency would be positively correlated with local fruit availability (P3).If play behavior reflects adult skill repertoires, we predicted object play would be more prevalent in the Suaq population, where feeding techniques are more complex (Schuppli et al., 2017) and include frequent, sophisticated stick tool use (van Schaik 1 3 Play Behavior Varies with Age, Sex, and Socioecological Context… et al., 1996), than at Tuanan, where tool use is absent (P4).Finally, if play behavior reflects adult sociability and serves socialization, we predicted that orangutans in the more sociable Sumatran population (Suaq) would show more social play than immature individuals of the less gregarious Bornean population (Tuanan) (P5).
In a semi-solitary species, variation in social play behavior may either be explained by association rates (i.e., opportunities to play), or interaction rates during these associations (i.e., the motivation to interact).Accordingly, if social play was mediated by association rates, we predicted that social play frequency would be higher for individuals from larger matrilines or with older siblings (P6) or in contexts (study sites) where association frequency is overall high (P7).If intrinsic motivation to play reflects higher adult sociability, we predicted the motivation to play to be higher at the more sociable population of Suaq compared to Tuanan (P8).Further, if play had a socializing function during ontogeny, we predicted a shift in play partners throughout immaturity (P9): from the mother and older sibling to additional association partners.Also following the socialization hypothesis, we predicted that mothers and/or siblings would not compensate for limited play (opportunities) with other associates (P10), but rather that association frequency with potential play partners would change over immatures' age (possibly regulated by the mother).Finally, if play interactions reflect adult sociability and social tolerance, we predicted more play with related individuals at Tuanan compared to Suaq, as a larger range of (unrelated) association partners may qualify as play partners in a setting of higher social tolerance (P11).

Study Site and Subjects
We conducted this study at two field sites in Indonesia: Suaq Balimbing (03°02'N; 97°25'E) in the Kluet region of Gunung Leuser National Park, South Aceh, Sumatra, and Tuanan (02°15'S; 114°44'E) in the Mawas Conservation Area, Central Kalimantan, Borneo.The study area in Suaq Balimbing consists of 500 ha primary, coastal peat-swamp forest.Suaq has the highest known orangutan density of 7 individuals per square kilometer (Singleton & van Schaik, 2001).The study area in Tuanan, Borneo, is situated in 750 ha of regenerating, formerly selectively logged peat-swamp forest.Tuanan has the highest reported orangutan density on Borneo with 4.5 individuals per square kilometer (van Schaik et al., 2005).Year-round, long-term studies on orangutan behavior are conducted at both field sites.For this study, we collected detailed behavioral data between June 2007 and March 2018 at Suaq and from July 2003 to June 2018 at Tuanan.
At both field sites, the studied orangutans are individually recognized.Based on demographic (van Schaik, 1999) and genetic data (Arora et al., 2012;Lenzi, 2014;Noordwijk et al., 2023), two major clusters of maternally related females (from here on referred to as matrilines) are known in Suaq.In Tuanan one central matriline of six females lives in the core of the study area and several mother-offspring units with fewer known close maternal relatives range more in the periphery of the study area (Table 1).At both sites, adult females have ranges that overlap with both closely related and non-related females (Arora et al., 2012;Ashbury et al., 2020;van Noordwijk et al., 2012).We accounted for the different matriline sizes in all analyses by adding a matriline size term, assigning immatures of females without known adult female relatives as small matrilines and all others to large matrilines (Table 1).While we analyzed play trajectories using absolute age (Table 1), we defined three different age classes of immature individuals to categorize social play partners.Dependent immatures comprise all not yet weaned individuals.Semi-dependent immatures do not nurse anymore, but still spend most of their time in association with their mothers.Independent immatures range independently of their mothers.We refer to these three categories collectively as immatures, and to both semidependent and independent immatures as weaned immatures.We knew the sex of all immatures, but we estimated some ages to the nearest half year (ages of young infants that were born in recent years were estimated to the closest month) (Table 1).At Suaq, we followed 12 immatures (5 females and 7 males, born between 2001 and 2016) regularly.At Tuanan, we followed 27 immatures (11 females and 16 males, born between 1999 and 2016) regularly (Table 1).

Behavioral Data
We used instantaneous behavioral sampling, following an established protocol (https:// www.aim.uzh.ch/ en/ resea rch/ orang utann etwork/ sfm.html) during full-day focal follows (mean 11.11 ± SD 1.10 focal hours per day), from the time when individuals left their morning nest to when they entered their evening nests.We recorded the focal individual's activity, visibility, and distance to other individuals at 2-min intervals.We defined three play categories as follows: Solitary Object Play: Manipulation of objects without an apparent immediate goal, often repetitive, e.g., swinging twigs, ripping off and discarding leaves, shaking, and manipulating sticks.We included explorative object manipulation (cf.Schuppli et al., 2016a, b) in this category.Solitary Locomotor Play: Repetitive movement of body parts or the whole body, like twirling up-side-down or swinging arms and/or legs.Solitary locomotor play can be in one spot, but also includes moving around with no obvious directed displacement.Social Play: Non-aggressive interaction between two or more individuals that does not serve any apparent, immediate purpose, like repetitive tagging, wrestling, and chasing.Social play includes body contact in most instances and is often accompanied by play faces or play vocalizations.
For each 2-min bout, we noted only one behavioral category.If the activity was a mix of categories at the 2-min scan, observers documented the behavior in a hierarchical order (social > solitary object > solitary locomotor play).Because the two solitary play types were often combined within a single play bout (mean 11.3 ± SD 6% of solitary play bouts is a combination, Fig. S1), we also analyzed overall solitary play frequency (Table S1).We labelled all individuals (including the own mother, younger or older siblings and any other individual) within a 50 m radius of the focal as association partners.Well-trained observers collected data, and passed an interobserver reliability test with an experienced observer with a Cohen's kappa of k ≥ 0.8 after an intense training phase.Several observers collected data at both study sites, further ensuring consistency in classification of play categories.

3
Play Behavior Varies with Age, Sex, and Socioecological Context…

Statistical Analysis
We used two main datasets for this study.First, we analyzed ontogenetic play trajectories and association frequency using longitudinal, 2-min 'instantaneous' scan data from full-day focal follows of immatures for whom at least 6 observation days were available over the study period (N = 1037 full-day focal follows; Table 1; Fig. S2).We included only data for dependent and semi-dependent immatures aged 11 years and younger, as males start to disperse around the age of 11 years (unpubl.data).Second, we used data on dyadic associations during full-day focal follows to analyze the occurrence of social play with specific partner classes (mother, sibling [older and younger] and associates [dependent, semi-dependent and independent immatures, and unflanged males]).We included data for dependent and semi-dependent immatures, who still ranged with their mother most of their active time, in these analyses (N = 2576 dyadic associations; Table 1).
We conducted data analysis using R version 4.2.1 (R Core Team, 2022).We performed Generalized Linear Mixed-effect Models (GLMMs) using the package 'glmmTMB' (Brooks et al., 2017).We checked models for heteroscedasticity and influential cases, also on the level of all random intercepts, using the package 'DHARMa' (Hartig & Lohse, 2021) and the model evaluation functions provided by the 'glmmTMB' package developers (https:// glmmt mb.github.io/ glmmT MB/ artic les/ model_ evalu ation.html).If we detected influential cases, we reran analyses excluding them to determine whether we found the same patterns.We examined models assuming a Poisson error distribution for overdispersion (Bolker, 2021), and used negative binomial GLMMs instead if overdispersion appeared to be an issue.We generated all figures using 'ggplot2' (Wickham, 2016) and 'ggeffects' (Lüdecke, 2018(Lüdecke, , 2020) ) to illustrate model predictions.We conducted likelihood ratio tests to compare the full models (including all fixed effects) to their corresponding null model (including only offsets and random effects).We only report second order interaction terms when adding them further improved model fit based on likelihood ratio tests and lower AIC values.We report model estimates and standard errors, exponentiated estimates and 95% confidence intervals from the confint() function in the result tables.We used the alpha value of P < 0.05 as cut-off for significance.

Ontogenetic Play Trajectories (P1-P6)
First, we tested if the three play types follow different age trajectories using a negative binomial GLMM.We pooled daily 2-min scan counts of each play type in one analysis and evaluated if the second order interaction terms between play type and the linear and quadratic term of age explained a significant amount of variation in the occurrence of play.We included the total count of daily 2-min scans with known activity as an offset term using a log-link function (we report details of accounting for visibility constraints in the supplementary materials).To avoid pseudo-replication resulting from having the same individual during consecutive days in the analyses, we included the follow period (FP), i.e., the calendar year and calendar month of focal observation, nested in individual identity (ID) as random intercepts in all the models (~ 1|ID/FP).An individual follow period comprised a mean of 3 ± 2.3 fullday focal follows (range: 1-11; Fig. S2).Additionally, we included observer identity and observation day as a random intercept.Second, we evaluated the total count of daily 2-min scans spent in solitary object, solitary locomotor, or social play during full-day focal follows in three separate analyses (N = 1037).We assessed variation in daily play scan counts in relation to study site, age (linear and quadratic term), sex, and zFAI in negative binomial GLMMs.We additionally assessed if time devoted to social play varied with matriline size (large vs. small) or with the presence of an older sibling (no vs. yes).To assess if patterns differ between populations, we tested for second order interactions of study site with any other fixed factor in the analyses of all three play types.We included follow period nested in individual identity and observer identity as random intercepts, and daily known activity scans as offset term using a log-link function, as for the first analysis.

Play Opportunities -Association Frequency (P7)
We analyzed (i) the number of and (ii) cumulative association time (in 2-min scans) with potential play partners based on full-day focal follows with a subset of the data for which information on association duration and partners were available (N = 731 follow days).We labelled dependent, semi-dependent and independent immatures, and unflanged males as potential play partners, because we did not observe the other sex-age classes (flanged males and adult females) playing with immatures.We did not include the individual's mother, with whom a dependent immature is always in association, or older siblings in these analyses.We included study site, age (linear and quadratic term), sex, matriline size, and zFAI as fixed factors and tested for any second order interaction terms of study site with the other fixed factors.As in the analyses of the play trajectories, we included follow period nested in focal identity as a random intercept.We included the total number of daily activity scans as an offset term using a log-link function when analyzing the duration of association.

Motivation to Play with Different Partners (P8-P11)
To investigate the variation in social play in detail, we first evaluated play occurrence during dyadic associations with different partner classes.We used a binomial GLMM to analyze the occurrence of social play (0/1) in a given dyadic association (N = 2576 dyadic associations [own mother: 1203; older sibling: 237; younger sibling: 241; other associate: 895]) with the total time in association (hours per day) as an offset term using a logit-link function.We included (i) focal identity nested in the follow period and (ii) the partner identity as crossed random intercept to avoid pseudo-replication.We added study site, age (linear and quadratic terms), sex and the partner class (mother, younger sibling, older sibling and associate) and zFAI as fixed effects.Further, we tested for second order interactions including either study site or age (either or both the linear and quadratic term) and all the other fixed factors to evaluate (i) potential site differences and (ii) age trajectory variation depending on other factors.

3
Play Behavior Varies with Age, Sex, and Socioecological Context… Second, we assessed the total count of observed social play scans in a given dyadic association per day.We evaluated play interactions with mothers, older siblings, and other associates in three different analyses.We set-up all GLMMs with a negative binomial error structure, included the total time in association as an offset term using a log-link function, and set the focal and partner identity as crossed random intercepts (~ 1|FocalID + 1|Partn-erID).In all three analyses, we assessed social play bout variation in relation to study site, age (linear and quadratic term), sex, matriline size and zFAI as fixed factors.When evaluating play frequency with mothers, we also added the presence/absence of an older sibling as a fixed effect.In analyses of play with associates (non-sibling and non-mothers), we added the category of the partner (unweaned [dependent] immatures, weaned immatures [including semi-dependent and independent] and unflanged males) as fixed effect.In all analyses, we tested for second order interactions including either site, age or sex and any other factor.

Ethical note
Behavioral data collection was non-invasive and exclusively observational.The data collection protocol for this study adheres to legal requirements of Indonesia and was approved by the Indonesian State Ministry for Research and Technology (RISTEK), the Directorate General of Natural Resources and Ecosystem Conservation-Ministry of Environment and Forestry of Indonesia (KSDAE-KLHK), the Ministry of Internal Affairs, Indonesia, the Nature Conservation Agency of Central Kalimantan (BKSDA) and Balai Besar Taman Nasional Gunung Leuser (BBTNGL).object play varied significantly with age (Fig. S3), but not with zFAI (Fig. S4), sex or the occurrence of social play on the same day (Table 2b).Solitary object play followed a quadratic age trajectory, increasing until around the age of 2 years and then decreasing towards the age of 6-8 years (Fig. 1).The age trajectories differed slightly but significantly between the two study sites, with a faster decline in solitary object play at Tuanan than at Suaq (Fig. 1), as shown by the significant interaction term between study site and the linear and quadratic age term (Table 2b).No other second order interaction term between study site and any of the fixed factors further enhanced the model fit.

Solitary Locomotor Play
The full model for solitary locomotor play explained significantly more variation than the corresponding null model (Table 2c).Solitary locomotor play showed a quadratic age trajectory which differed between the two study sites (Fig. S3), with a much higher peak in Tuanan than in Suaq (Table 2c; Fig. 1).Solitary locomotor play frequency did not vary significantly with local fruit availability (Fig. S4), sex or the occurrence of social play (Table 2c).No other interaction terms between study site and any of the fixed factors further improved the model fit.

Social Play
The full model explained significantly more variation in social play scans than the corresponding null model (Table 2d).Age trajectories in social play differed between the two study sites (Fig. S3).Social play followed a quadratic age trajectory at Suaq, the more sociable population, and peaked around the age of 4-6 years, whereas at Tuanan this peak was less pronounced, and social play frequency remained low(er) at Tuanan (Table 2d; Fig. 1).There was a significant interaction between zFAI and study site (Table 2d): Daily social play increased with increasing zFAI at Tuanan, whereas there was no such correlation at Suaq (Fig. S4), where absolute fruit availability is generally higher.Daily social play scans did not differ significantly between sexes, matriline size or if the immature had an older sibling (Table 2d).

Play Opportunities-Association Frequency (P7)
As predicted, both the mean number and the mean cumulative time spent in associations with potential play partners, including dependent, semi-dependent und independent immatures, and unflanged males, were higher at Suaq (1.25 ± SD 1.44 potential play partners per full-day focal follow and 4.45 ± SD 7.98 cumulative association hours) than at Tuanan (0.78 ± SD 0.84 partners and 1.22 ± SD 2.92 cumulative hours).While both number and duration of associations with potential play partners increased with local fruit availability at Tuanan, they remained constant or even decreased at Suaq, as the interaction between site and zFAI shows (Table 3; Fig. S5).The overall number of association partners did 1 3 Play Behavior Varies with Age, Sex, and Socioecological Context… not vary with age or sex of the immature (Table 3).However, in separate analyses of association partner types, the number of dependent immatures and unflanged males increased and the number of weaned immatures decreased with the age of immatures (Table S2; Fig. S6).

Probability of Play During Dyadic Association
The full model explained significantly more variation in social play occurrence during associations than the corresponding null model (Table 4).During an association, the probability of social play between the dyadic association partners varied with (i) partner type dependent on immatures' age (quadratic effect) (Fig. 2), (ii) sex, and (iii) study site depending on local fruit availability (Table 4).First, the probability of social play in an association with the older sibling peaked at an earlier age and decreased sooner than with other partners.Moreover, the probability of social play with the individual's mother decreased at a younger age than with other (non-sibling) associates (Table 4).Second, male immatures were more likely to engage in social play than female immatures in a given dyadic association (Table 4).The second order interaction between sex and either (i) study site, (ii) partner type, or (iii) age did not further improve the model fit.Finally, the probability of social play in a dyadic association increased with increasing zFAI at Tuanan, but did not vary with zFAI at Suaq (Table 4).

Play Frequency with Own Mother
We observed immature focal individuals playing with their mother on 17.3 ± SD 37.9% (Suaq: 21.6 ± SD 41.2%; Tuanan: 16.3 ± SD 36.9%) of full-day focal follows and in a mean of 0.75 ± SD 2.51 (Suaq: 1.01 ± SD 2.62; Tuanan 0.68 ± SD 2.47) daily 2-min scans.The age trajectory of total daily play scans with the mother followed a quadratic function (Table 5a).While social play with their own mother was recorded for only a few scans per day at both sites (on average), the frequency increased somewhat with increasing zFAI at Tuanan, but remained constant with zFAI at Suaq, as indicated by the significant interaction term (Table 5a; Fig. S8).Social play with the own mother did not vary with sex or study site.

Play Frequency with Older Sibling
When they were in association, immature focal individuals played with their older siblings during 56.9 ± 49.9% (mean ± SD) of full-day focal follows at Suaq, and 33.7 ± 47.4% at Tuanan.Immature focal individuals spent more 2-min scans in social play per time in association with their older siblings at Suaq (4.0 ± SD 6.8) than at Tuanan (3.1 ± SD 7.1) (Table 5b).Daily play frequency (2-min scans) with older siblings followed a quadractic age trajectory but did not vary significantly with the sex of either the focal individual or the older sibling or with zFAI (Table 5b).

Play Frequency with Associates
The amount of social play with other associates, including dependent, semi-dependent, and independent immatures, and unflanged males, did not vary significantly with zFAI, study site, matriline size or if the focal had an older sibling, but did vary significantly as (i) a quadratic function of age, (ii) an interaction of focal sex and partner class, and (iii) an interaction of study site and if play partners were maternal relatives (Table 5c).First, male and female immatures did not differ significantly in their time playing with dependent immatures (male focal: 10.0 ± SD 24.4 scans; female focal: 14.8 ± SD 25.3 scans), but male immatures spent relatively more time in social play with weaned immatures (male focal: 5.4 ± SD 10.5 scans; female focal: 2.5 ± SD 8.8 scans) and unflanged males (male focal: 1.6 ± SD 5.5 scans; female focal: 0.1 ± SD 0.5), compared to female immatures (Fig. 3; Fig. S7 for distribution of raw data).Second, immatures' play frequency was higher with maternally related association partners than unrelated association partners at Tuanan, but not at Suaq, where play frequency did not depend on maternal relatedness (Fig. S9).

Discussion
Play behavior in orangutans varies with age, sex, local fruit availability and socioecological context.As predicted, our results showed that solitary object, solitary locomotor and social play follow different age trajectories (P1), whereas we could not find any evidence for compensation effects between solitary and social play (P2).While we did not account for the combination of different play types, particularly of solitary play, our results may reflect developmental constraints on the expression of    1.97 play or indicate that different play types serve different developmental functions.As predicted, fewer opportunities (P7) appear to cause the lower social play frequency at the less sociable population of Tuanan compared to Suaq (P5).Supporting our predictions and the socialization hypothesis, play partners changed during ontogeny (P9), and mothers did not compensate for the limited number of play opportunities with other play partners (P10).While solitary play was not affected by local fruit availability, social play decreased in all dyads at Tuanan at lower fruit availabilities (P3).Altogether, our findings support hypotheses that play behavior reflects adult behavioral repertoire and sociability, while being constrained by energy.

Solitary Locomotor and Solitary Object Play Trajectories
Solitary play might be linked to skill acquisition during the first few years of life: Both solitary object and solitary locomotor play were most frequent during the first 4 years of life, and ceased almost completely around the age of weaning at 6 to 7 years.First, the drop in solitary locomotor play coincides with the start of consistent independent travel when the offspring is rarely carried by the mother during  travel (Chappell et al., 2015;van Noordwijk et al., 2009).The trajectory of solitary locomotor play differed between study sites, with a more pronounced peak at Tuanan than at Suaq.We did not predict a site difference in solitary locomotor play, but it may be related to the different forest structure at the two sites, which may require  (Chappell et al., 2015;Manduell et al., 2012).The trees at Suaq tend to be large and thus the canopy is more 'rugged' than at Tuanan, where many small and still flexible trees fill most gaps.Future studies comparing locomotor play across study populations with varying forest structures are needed to further elucidate if locomotor play is more frequent in habitats with more challenging forest structures and may reflect increased need for motor skill training.
Second, solitary object play remained prevalent at older ages at Suaq (P4), where complex feeding techniques, including habitual, sophisticated (stick) tool use, are regularly observed (Schuppli et al., 2017;van Schaik et al., 1996), unlike at Tuanan where feeding techniques appear less complex and tool use is absent.However, the difference in solitary object play was modest.Moreover, immatures become proficient tool users well after the drop in solitary object play, around the age of 8 years at the earliest (Meulman et al., 2013).Solitary object play may thus be a first step in manipulation skills and understanding the physical concepts needed for tool use, combined with social learning (Schuppli et al., 2016a, b;sensu van Schaik & Burkart, 2011), and more fine-tuned practice in the form of exploratory activities (Schuppli, van Cauwenberghe et al., 2021b).More sophisticated adult behavior may further be foreshadowed by increasing complexity, particular manipulation, or specific types of objects during play rather than by absolute rates of play (Schuppli, van Cauwenberghe et al., 2021b).In the genus Pan, chimpanzees in populations with frequent stick tool use showed a higher use of detached objects with increasing age, and more varied object manipulation than bonobos, where such stick tool use is absent (Koops et al., 2015a, b;Koops et al., 2015a, b;Myowa-Yamakoshi & Yamakoshi, 2011).Object manipulation quality also foreshadowed particular skill sets in a chimpanzee population without habitual tool use (Lamon et al., 2018).Indonesia, 2003Indonesia, -2018) ) and Sumatran orangutans (P.abelii at Suaq Balimbing, South Aceh, Indonesia, 2007Indonesia, -2018)).Symbols indicate model predictions and error bars show 95% confidence intervals obtained from the model (Table 5c)  1 3 Play Behavior Varies with Age, Sex, and Socioecological Context… Future studies will show whether the nature of object play, such as the objects used, or distinct manipulative actions performed, rather than frequencies per se further reflect the differences between populations with or without complex feeding techniques in orangutans.Object play may generally serve motor skill acquisition (e.g., geladas [Theropithecus gelada]: Cangiano & Palagi, 2020).The slightly extended period of solitary object play at Suaq compared to Tuanan might be linked to slower developmental trajectories overall, as overall slower development and a slower acquisition of feeding skills has been suggested for Sumatran compared to Bornean orangutans (Schuppli, et al., 2016a, b;van Noordwijk et al., 2009).For now, we conclude that there is only a minor predisposition or motivational difference for increased overall play frequencies with objects in early life phases between populations with different levels of feeding technique complexity in orangutans.

Social Play
As predicted, social play frequency reflected population sociability (P5).At the less sociable population of Tuanan, social play remained constant throughout immaturity and decreased towards adolescence, with large variation between days.The time spent on social play was generally higher and peaked around the age of 4 to 6 years at Suaq, the more sociable population.The higher individual social play frequency at Suaq compared to Tuanan, however, is the direct result of the increased association frequency, and hence more play opportunities (P7), rather than the immatures' intrinsic motivation to play (P8).At Tuanan, immatures (i) were as likely to play and (ii) spent as much time in social play within a given association as at Suaq (cf.Fröhlich et al., 2020).The only exceptions were older siblings with whom play frequency was generally higher at Suaq than at Tuanan, which can likely be explained by the energetic costs of play.Indeed, the lower association frequencies and hence fewer play opportunities at Tuanan compared to Suaq may increase the immatures' motivation to play once they get the chance, and their mother associates with potential play partners.Thus, we may see a compensation effect at Tuanan, which is not necessary at Suaq due to the higher availability of play partners.Associations also lasted longer at Suaq and there may have been a saturation effect on motivation or energy to play.In sum, these findings support the basic rule 'play when you get the opportunity' unless limited by energy.However, we cannot infer the directionality of the relationship between amount of social play and adult social competence posited by the social skills hypothesis (Baldwin & Baldwin, 1974) or the socialization hypothesis (Fagen, 1981).

No Compensation for Limited Play Opportunities by the Mother
The availability of play partners is a limiting factor for social play in orangutans given their semi-solitary lifestyle.Mothers do not appear to make up for the limited play opportunities for their offspring (P10).First, when associations with other partners were rare, during extremely low fruit availability, play frequency with mothers was also low.Second, neither mothers nor siblings in small matrilines, with fewer observed social play opportunities (van Noordwijk et al., 2012), played more with their offspring or younger sibling, respectively.Third, play partner identity transitioned from the mother to other individuals with increasing age supporting our prediction (P9).From the age of 4 years onwards, immatures hardly played with their own mothers.Our finding of changing play partners during ontogeny supports the socialization hypothesis and is in accordance with reports on released orangutans (Mendonça et al., 2017;Rijksen, 1978), chimpanzees (Cordoni & Palagi, 2011), bottlenose dolphins (Tursiops truncatus: Mackey et al., 2014) and humans (Pellegrini, 2009).It remains to be investigated whether play quality in wild orangutans also changes with age, or play partner and context, as for example in chimpanzees (Cordoni & Palagi, 2011;Flack et al., 2004).

Do Mother Drive Socialization in their Offspring
Social play facilitates the acquisition of social competence and rank (e.g., chimpanzees: Shimada & Sueur, 2014), but social competence is also required to maintain social play bouts (e.g., van Leeuwen et al., 2013).In species with high fission-fusion dynamics, such as orangutans or chimpanzees, associations tend to be costly for nursing mothers (Kunz et al., 2021), which limits the access their dependent offspring have to potential playmates.Active play partner choice by immatures may therefore rarely be feasible, unless mothers drive socialization via selective associations (Murray et al., 2014;Williams et al., 2002).We found only limited evidence for this: while overall association frequency did not vary with the age or sex of the immature, association frequency with different partner types varied with the age of immatures, which might indicate some selection of association (and play) partners by the mothers.Further, mothers' preference for associating with maternal relatives has been suggested to lead to more social play opportunities for immatures in larger matrilines (van Noordwijk et al., 2012).We did not observe overall differences in social play frequency in immatures from larger versus smaller matrilines but we did find a difference in social play frequency with maternally related associates compared to unrelated associates at Tuanan, although not at Suaq.This finding indicates that maternal relatives are indeed preferred play partners, at least in one less sociable population, where social tolerance may be a limiting factor for play to occur.

Sex Differences in Social Play Probability
While overall social play frequency did not differ between male and female immatures, immature male orangutans were more likely to engage in social play during dyadic associations and showed higher play frequencies with other weaned immatures and unflanged males than immature females did.Similar sex differences in orangutan social play behavior have been found in rehabilitation centers (Descovich et al., 2011), and they match sex differences in adult orangutan skill sets and their acquisition (Ehmann et al., 2021;Schuppli et al., 2021a, b).Immature male orangutans eventually leave their natal home range (Arora et al., 2012), and show relatively high levels of sociability post-dispersal and even regular playful 1 3 Play Behavior Varies with Age, Sex, and Socioecological Context… interactions (Fox, 2002;Galdikas, 1985;Mörchen et al., 2023, n.d;Utami Atmoko et al., 2009).Whether social play serves social skill acquisition, to practice fighting skills or to reduce tension (e.g., Verreaux's sifaka [Propithecus verreauxi]: Antonacci et al., 2010;Palagi, 2018) in (male) orangutans remains to be tested with more detailed data on social play quality.Unfortunately, longitudinal studies evaluating play behavior in relation to social and skill competence in adolescent and unflanged male orangutans are hardly feasible in the wild due to their natal dispersal (Arora et al., 2012).Immature female orangutans, in contrast, appeared to play slightly more with dependent immatures than with weaned immatures and unflanged males.Although our data set is one of the largest available on wild orangutan populations, we are still limited in accounting for potential confounding variables such as demography, dyad composition, and age and sex differences between play partners due to orangutans' slow life history (van Noordwijk et al., 2018) and semi-solitary lifestyle (van Schaik, 1999), which make associations and social interactions rare.Thus, the sex differences we found should be treated with caution, although they are in line with earlier studies of orangutans and other more gregarious great ape species (chimpanzees: Kahlenberg & Wrangham, 2010;Koops et al., 2015a, b;Lonsdorf et al., 2014a, b;Pusey, 1990;gorillas: Maestripieri & Ross, 2004).

Energetic Constraints to Play
Play is presumably only feasible when animals have surplus energy (Berghänel et al., 2015;Held & Špinka, 2011;Palagi, 2018;Spencer, 1872).Accordingly, social play frequency increased with increasing local fruit availability at Tuanan, although this increase was a few minutes per day at most, whereas it did not vary with local fruit availability at Suaq.Although these differences may be in part attributed to association patterns, and hence, decreased opportunities to play at lower local fruit availability at Tuanan, the motivation, or rather energy, to play appears to decrease with decreasing local fruit availability at Tuanan.Moreover, social play frequency with the mother also decreased at lower fruit abundance at Tuanan.While this effect was small and social play with the mother is usually very gentle, it suggests that mothers might be particularly reluctant to play when only low-value nutrition is available.This pattern fits the notion that mothers oversee social play opportunities for their dependents for which they have to provide the additional energy via milk and transport (carrying) (van Noordwijk et al., 2013).Importantly, social play did not cease completely even during extended periods of very low fruit availabilities (Ashbury et al., 2022).This finding suggests that the benefits of low but persistent rates of social play outweigh the costs of social play even during lean periods.
Contra our prediction (P3), solitary play frequency did not vary with local fruit availability and thus, immatures appear to be cost-insensitive.Mothers have less control over solitary play despite acting as an energy buffer for their offspring.Therefore, it is not surprising that solitary play ceases by the time of weaning.The energy expenditure for solitary play may thus represent a hidden mother-offspring conflict (sensu Trivers, 1974).

Conclusion
We show that solitary object, solitary locomotor, and social play follow different ontogenetic trajectories in orangutans.Play seems crucial during the ontogeny of immature orangutans, as the intrinsic motivation to play appears largely independent of the socio-ecological context.However, social play frequency is shaped by variation in sociability and ecology across two populations.Assessing the potential functions of play and specifically, whether individual play experience translates into adult competence, survival, and reproductive success, will only be feasible through extended longitudinal studies, and will be easier in the philopatric females.Our study supports growing evidence that both solitary and social play behavior changes during ontogeny and its variation reflects adult behavioral repertoires.While play behavior may be limited by energy and skill development to some extent, the intrinsic motivation of immature individuals to play suggests that it has developmental functions.

Fig. 1
Fig. 1 Play trajectories: Predicted number of daily 2-min scans spent in play behavior by immature orangutans by age (years) (x-axis), play type and study site (left: Sumatran orangutans (P.abelii) at Suaq Balimbing, South Aceh, Indonesia, 2007-2018; right: Bornean orangutans (Pongo pygmaeus wurmbii) at Tuanan, Central Kalimantan, Indonesia, 2003-2018).Lines indicate model predictions based on the full model for each study site and shaded areas illustrate 95% confidence intervals Effects of individual, social and ecological variables on the occurrence of social play in a given dyadic association, derived using a GLMM with binomial error structure (N = 2576 dyadic associations during full-day focal follows of immature focal individuals [Bornean (Pongo pygmaeus wurmbii at Tuanan, Central Kalimantan, Indonesia, 2003-2018) and Sumatran orangutans (P.abelii at Suaq Balimbing, South Aceh, Indonesia, 2007-2018)]).Comparisons to the null model based on a likelihood ratio test are indicated below the response variables.Significant fixed effects (P < 0.05) are indicated in bold

Fig. 2
Fig. 2 Probability of social play (0/1) during associations with different partners (associates, the immatures' own mother and their older sibling) by the age of immature (years) in Bornean (Pongo pygmaeus wurmbii at Tuanan, Central Kalimantan, Indonesia, 2003-2018) and Sumatran orangutans (P.abelii at Suaq Balimbing, South Aceh, Indonesia, 2007-2018).Lines indicate the predicted probability of social play during a given dyadic association.'Rugs' (i.e., the short vertical dashes at the bottom and top of the panel) indicate the distribution of associations (bottom: association, but no play observed; top: association and play observed), and shaded areas the 95% confidence intervals

Fig. 3
Fig. 3 Predicted number of 2-min social play scans during association with different partners by the sex of immature focal individuals in Bornean (Pongo pygmaeus wurmbii at Tuanan, Central Kalimantan, Indonesia, 2003-2018) and Sumatran orangutans (P.abelii at Suaq Balimbing, South Aceh, Indonesia, 2007-2018).Symbols indicate model predictions and error bars show 95% confidence intervals obtained from the model (Table5c)

Table 1
(continued) Arora et al., 2012;Lenzi, 2014;van Noordwijk et al., 2023 et al., 2023† only selected observers (researchers focusing on immature development) # focal data between 2010 and 2018 from dependent and semi-dependent immatures included, all association partners counted here (own mother, other dependent, semi-dependent and independent immatures, adult females, and unflanged and flanged males) k known birth date Birth months and years in italics indicate estimates for immatures present at the start of the study period

Table 2 (
Play Behavior Varies with Age,Sex, and Socioecological Context…

Table 3
Effects of individual, social and ecological variables on the number of play partners and cumulative duration of association with them (includes dependent, semidependent and independent immatures, and unflanged males) during full-day focal follows on Bornean (Pongo pygmaeus wurmbii

Table 5
Effects of individual, social and ecological variables on the count of 2-min social play scans per full-day focal follows with a) the mother, b) the older sibling, and c) association partners (unweaned and weaned immatures, and unflanged males) in Bornean (Pongo pygmaeus wurmbii at Tuanan, CentralKalimantan, Indonesia, 2003-2018)and Sumatran (P.abelii at Suaq Balimbing, South Aceh, Indonesia, 2007-2018) orangutans.Comparisons to the null models (fixed effect = 1, offset and random intercepts) based on likelihood ratio tests and the sample size (dyadic associations with the partner category on full-day focal follows) are indicated below the response variable.Significant fixed effects (P < 0.05) are indicated in bold Play Behavior Varies with Age,Sex, and Socioecological Context…