Introduction
Parental care often requires substantial investment of time and energy, and strongly impacts the fitness of the individual carers that provide it (Saether, 1994). Previous studies have shown that care behavior is influenced by individual characteristics such as sex (Liker et al., 2015), age (Ortega et al., 2017), condition (Dearborn, 2001), and personality (Westneat et al., 2011). However, many studies only measure the amount contributed by each carer within one behavioral dimension of care (e.g., food delivery) and during one stage of offspring development (e.g., provisioning nestlings). We currently know relatively little about how carers contribute across multiple dimensions of care behavior, or how the distribution of contributions impacts outcomes for carers and offspring. Similarly, while many studies have explored how carers change the amount they contribute according to the contributions of others (reviewed in Hatchwell, 1999; Harrison et al., 2009), the fine-scale behavioral rules underpinning carer interactions have only recently attracted serious attention (Johnstone et al., 2014).
In this mini-review, we discuss how care behavior can be quantified, classifying measurements into three broad groups: how parental care is distributed (i.e., when care occurs along the timeline), the characteristics of care (i.e., variation among different instances of similar care behavior), and the interactions between carers (i.e., whether one carer's behavior is associated with the behavior of others). Our aim is to illustrate the variety of questions that can be explored using datasets on parental care, and some of the statistical and technical considerations that arise when doing so. We review existing studies and analysis methods that have addressed these different aspects of care behavior, and briefly discuss potential future research directions.
Distribution of Care
Rate and Variance
When individuals deliver discrete, relatively brief care events to offspring (e.g., provisioning, defense), care/visit rates or the mean and variance of carer inter-visit intervals (IVIs) are useful metrics to quantify behavior. Most literally, IVIs refer to the periods between an individual leaving a nest or den and its next arrival (Santema et al., 2017), but IVIs are also commonly characterized as the time between consecutive arrivals (e.g., Johnstone et al., 2014); for clarity we refer to this latter case as the inter-arrival interval (IAI). When discrete care events are somewhat longer (e.g., nest maintenance) it can also be informative to characterize the within (or intra-) visit intervals (WVI) of carers. Passive integrated transponder (PIT) tags are increasingly used to collect large amounts of visit data on provisioning behavior, particularly in cooperative systems with many carers (Browning et al., 2012), but are less valuable when care occurs away from fixed locations like a nest.
The distributions of intervals (of all types) can be compared between individuals and contexts to understand variation in care. Intervals are often approximately gamma- or inverse-gamma distributed, as they are bounded at zero and often have a (soft) minimum duration that depends on the type of care delivered. They can hence be defined using two independent parameters, scale and shape, that reflect their rate and skewness/variance (Lejeune et al., 2019). In biparental and cooperative systems, the distribution of IAIs by the entire care group is more likely to predict breeding success than those of individual contributors, as the overall amount and distribution of care is what determines outcomes for offspring. Outlier IVIs and IAIs may represent carers taking breaks from caring (e.g., due to self-foraging or disturbance), and hence are useful for characterizing and partitioning longer sample periods to avoid applying inappropriate analyses. When analyzing samples from longer periods of care, intervals (and hence analyses based on them) can be biased as the beginning and end of the sample periods are more likely to cut longer intervals; where possible studies should use either naturally bounded periods, or ensure their sample contains many events and acknowledge the bias (Baldan et al., 2019).
For care delivered over substantial periods (e.g., incubation, babysitting) the proportion of active time carers spend on care, or the proportion of opportunities during which care occurs, are more suitable metrics than the intervals between care events. A typical way to model effects on proportional care is a logistic (Bambini et al., 2018) or binomial (Clutton-Brock et al., 2000) regression when the proportion is derived from counts, or a beta/Dirichlet regression when it is based on continuous numbers (Douma and Weedon, 2019).
Trends
The rates at which individuals deliver care may vary across a sample period, driven by environmental variation (e.g., weather) or the states of parents or offspring (e.g., hunger). Such variation will affect parental care over the same period, and can limit the usefulness of randomizations used to infer interactions between carers (Baldan et al., 2019). To quantify trends, one simple metric is how strongly intervals are ordered in time (Schlicht et al., 2016), which will identify a linear increase or decrease in rate. A more detailed picture can be obtained by explicitly fitting a model of interval length, with linear and higher-order time terms as predictors. For more complex trends, especially those with periodicity, one could investigate temporal patterns of care using methods developed for time series analysis that have been previously applied to other aspects of behavior, such as cross-correlations (Hall et al., 2014) or wavelet analysis (Zhang et al., 2017). An alternative approach is to group care events by hour or by day and then fit Poisson-based mixed models to the counts of care behavior, with environmental metrics as covariates (Nomano et al., submitted); the best approach will depend on the study system and length of time analyzed.
Repeatability
In addition to measuring care variation and trends within an observation period, one can also evaluate whether carer behavior is individually repeatable between observations, and hence infer whether that behavior can be regarded as an individual trait. Repeatability is typically defined as the proportion of variance attributable to the differences among groups of observations, before or after controlling for the effect of confounding factors on the response variable (“adjusted repeatability” in the latter case). Detailed guidance on how to measure and interpret repeatability is beyond the scope of this review; see Nakagawa and Schielzeth (2010) for extensive discussion, and the R package “rptR” for useful analysis methods (Stoffel et al., 2017). While there have been several studies on the repeatability of provisioning (Potti et al., 1999; Nakagawa et al., 2007) and other care behaviors including babysitting (English et al., 2010; Sanderson et al., 2015) there is substantial scope for further study.
Multi-Stage Investment
Empirical studies often focus on a single stage of parental care, or assume each stage is broadly independent. However, theory suggests that investment during earlier stages (e.g., egg-laying) can influence investment rules in later stages (e.g., provisioning), particularly if the breeding female can control the number (Smith and Härdling, 2000; Savage et al., 2013) or quality (Savage et al., 2015) of offspring. These predictions are supported by empirical studies (Russell et al., 2007; Canestrari et al., 2011), but our understanding of multi-stage investment dynamics is still incomplete, particularly how these are influenced by environmental conditions (Langmore et al., 2016). As with any form of adaptive plasticity, for between-stage strategies to evolve the environment must be both variable and predictable on the timescale of the care periods, a concept familiar from research on transgenerational effects (Proulx and Teotónio, 2017).
Characteristics of Care
Variation Among Events
Not all instances of a particular care behavior are equivalent from the perspectives of either parents or offspring. For example, provisioned food can vary in mass or nutritional content, and hence carers can vary in contributions without differences in provisioning rate. The size of prey delivered can vary with carer sex (Colombelli-Négrel and Kleindorfer, 2010), these differences may vary with offspring age (Wiebe and Slagsvold, 2009), and males and females may differ in the prey type delivered to offspring (Fraser et al., 2006). Similar considerations apply to other forms of parental care; for example mobbing behavior can vary not only in its frequency but also in the intensity of each event (e.g., contact vs. non-contact) (Strnad et al., 2012), and incubation may be costly to younger but not older carers (Heinsohn and Cockburn, 1994). Characteristics of care events can also affect IVIs, for example larger food items being brought to offspring after parents have been away for longer (Grieco, 2002), altering inferences about carer investment.
Favoritism
When parents deliver care to multiple offspring, the amount each offspring receives is important for the overall outcome of the breeding attempt. Certain offspring may be more dominant or beg more intensively (Drummond, 2006), or carers may preferentially feed some offspring over others due to expected returns (Jeon, 2008) or favor different offspring if their costs or benefits differ (Lessells, 2002). For example, male offspring may receive more food if mothers preferentially provision sons over daughters (Mainwaring et al., 2011), and parents may adjust which offspring they provision based on offspring age and perceived quality (Avilés et al., 2011).
One metric to characterize this variation is a “skew index” (Pamilo and Crozier, 1996; Shen et al., 2010), which varies from 0 (complete equality) to 1 (one offspring receives all the food). When care to individual offspring can be quantified precisely, for example using video cameras deployed inside nests, comparing between offspring any of the above metrics for the distribution of care can also reveal differences in carer delivery behavior.
Non-care and Deception
Carers sometimes visit offspring without providing care. For example, individuals bringing food to dependent offspring may consume it themselves in so-called “false feeding” behavior, perhaps representing a deceptive strategy to lower the costs of care (Boland et al., 1997), or a non-deceptive mediation of carer need against offspring need (Canestrari et al., 2010). Identifying deception is non-trivial, as apparent false-feeding can occur when carers visit offspring that are fully satiated; approaches to identify deception include using remote video monitoring to score offspring begging behavior or attempts to feed offspring (Young et al., 2013), measuring latency between arrival and self-consumption (or departure with food item) compared to normal visits, and testing whether false-feeds are less likely to occur when the provisioner is observed (Boland et al., 1997; Young et al., 2013). In some species such as the bell miner (Manorina melanophrys) carers may also only partially deliver food items; whether these events are treated as false-feeds or not should depend whether they can be more parsimoniously explained by (for example) difficulties in prey transfer (McDonald et al., 2007).
Visits to offspring in which carers arrive without food (and provide no other care) could also be classified as false-feeding, but might instead represent carers updating information about offspring hunger. When carers cannot easily monitor each other's contributions, and visiting offspring is much less costly than finding and delivering food, this additional information is especially valuable to correctly distribute care.
Interactions During Care
Negotiation
As the benefits of care are shared but the costs personal, carers have a conflict of interest over how much each contributes, and this should affect their investment decisions (Trivers, 1972). Theory suggests that individuals in biparental species should respond to changes in the contributions of others by incompletely compensating, both over evolutionary (Houston and Davies, 1985) and behavioral (McNamara et al., 1999) timescales. Further models suggest that incomplete compensation should also occur in cooperative systems (Johnstone, 2011), and that high responsiveness (McNamara et al., 2003), asymmetric information (Johnstone and Hinde, 2006) or threshold effects (Jones et al., 2002) can modify predictions. Empirical work on biparental species largely supports incomplete compensation as the usual response to changes in partner investment, albeit with substantial variation (Harrison et al., 2009) and often sex differences (e.g., Iserbyt et al., 2015). However, cooperative species adopt more diverse investment rules (Hatchwell, 1999), including responding to both the composition and size of the care group (Brouwer et al., 2014) potentially due to greater variation in care during later stages (Savage et al., 2013).
Testing theoretical predictions about negotiation requires careful experiments to manipulate offspring demand (actual or perceived) or carer costs, and monitoring parental responses. The mechanisms through which individuals negotiate are still poorly understood; vocal communication is likely to play a major role (Bell et al., 2010; Boucaud et al., 2016), and as negotiations could also be mediated indirectly through offspring need (Lessells and McNamara, 2012) negotiation behavior is likely to be highly system-specific.
Alternation
Theoretical work suggests that “turn-taking”—carers alternating contributions—can resolve sexual conflict efficiently (Johnstone et al., 2014), and provided the costs and benefits of care are time-dependent this does not require individuals to monitor each other perfectly (Johnstone and Savage, 2019). Empirical work has suggested that several species indeed alternate more than expected by chance, however questions remain over the mechanism of interaction and how strongly this turn-taking is driven by environmental variation vs. individual responsiveness (Ihle et al., 2019).
Turn-taking can be quantified using the proportion of alternated visits (Bebbington and Hatchwell, 2016; Iserbyt et al., 2017) or the log-odds of the deviation between observed and expected number of alternated visits (Baldan et al., in press). Alternatively, the runs test (Wald and Wolfowitz, 1940) can investigate whether two carers alternate more than expected (Johnstone et al., 2014), and a modified version of the test is also applicable to cooperative species (Sheskin, 2011; Khwaja et al., 2017).
Incorporating time information as well as visit sequences, continuous-time Markov models can be used to investigate patterns of alternation in biparental (Johnstone et al., 2014) and cooperative (Savage et al., 2017) species. These analyses can be applied in R using packages such as “msm” (Jackson, 2011), or more simply calculated directly from visit times if within-observation covariates are not required (Savage et al., in review). If enough data exist to characterize the distribution of care intervals precisely, an alternative approach is to fit a semi-Markov model explicitly using the relevant distribution. Such models can be implemented using (e.g.,) the “SemiMarkov” package in R (Król and Saint-Pierre, 2015), although to our knowledge this method has yet to be applied to care behavior.
Synchrony
Depending on the system and behavior in question, pairs or groups synchronizing care activities might have either a positive or negative impact on the success of a breeding attempt. Visiting offspring can increase predation risk (Martin et al., 2000), leading to groups that synchronize visits having increased breeding success (Raihani et al., 2010). Similarly, synchronizing visits might reduce sibling competition by providing resources to more offspring simultaneously (Shen et al., 2010). In contrast, if visits do not increase predation risk, carers deliver multiple (or divisible) food items, and offspring satiate quickly, then carers should deliberately separate their care contributions (anti-synchrony). Beyond the impacts on offspring, individuals might also benefit from synchronizing their visits to advertise their contributions to or monitor other group members (Doutrelant and Covas, 2007), particularly in a “pay-to-stay” cooperative system (Gaston, 1978; Kokko et al., 2002).
One method to quantify synchrony is to characterize a particular visit as synchronous when another individual also provides care within a certain window (Mariette and Griffith, 2012) and then use the square root arcsine–transformed proportion of synchronous care events as a measure of overall synchrony (Mariette and Griffith, 2015). This is appropriate in systems with relatively low care rates and brief care events, but can be sensitive to the window chosen. Alternatively, for more frequent or longer care behaviors one can cross-correlate the time series of care contributions by each individual (Savage et al., 2017), and for both methods randomizations can be used to generate expected levels of synchrony. Potentially useful analysis methods have also been developed in neurobiology, where quantifying the relationships between a number of neural spike trains is a common problem (Oram et al., 2001; Shimazaki et al., 2012), however these have yet to be applied to care behavior.
Alternation and synchrony together provide a good picture of individual interactions (Koenig and Walters, 2016), and investigating both is also important because the interpretation of each metric depends partly on the other. For example, a strict pattern of alternation may suggest that individuals are adopting a turn-taking rule under low synchrony, but under high synchrony an alternative explanation would be that individuals forage together and individual differences (e.g., from state or personality) result in one consistently visiting before the other.
Task Specialization
In many species care occurs simultaneously across multiple behavioral dimensions. For example, in an altricial bird carers may need to feed and brood offspring, maintain the nest, remove fecal sacs, and mob nest predators. Pairs and groups of carers may be comprised of individuals that specialize in particular behaviors and/or generalist individuals, driven by differences in the costs and benefits of each behavior (Arnold et al., 2005). The degree of specialization can change over time (Iserbyt et al., 2017), and within activities carers may sub-specialize (e.g., by food type) or partition roles in time (e.g., helpers rearing first broods while breeders re-nest; Ridley and Raihani, 2008). Comparing parental investment across modalities can be challenging as costs are often accrued in a different “currency” for different care activities and contexts (e.g., mortality via predation risk when mobbing vs. condition via lost self-foraging time during provisioning), but such comparisons are important as these behaviors trade-off against each other (Mutzel et al., 2013).
To investigate factors influencing (e.g.,) the type of prey being delivered to offspring, one approach is to fit the proportional abundance of each prey type as response terms in a (mixed-effect) multinomial logistic regression, to avoid conflating variation in the proportion of each prey type with that of the others (Browning et al., 2012). These models can be most precisely fitted using Markov chain Monte Carlo (MCMC) Bayesian methods, for example with the R packages “RStan” and “MCMCglmm” (Hadfield, 2010; Stan Development Team, 2018). Correlations among the random effects of such models are potentially informative for elucidating individual trade-offs among care behaviors; for discussion of these effects and a detailed treatment of methods around the multinomial analysis of behavior see Koster and McElreath (2017).
Discussion
In this mini-review we have illustrated that how care is distributed, how care events vary, and how carers interact, each have important consequences for carers and offspring. Our review also illuminates a number of questions that remain despite the vast literature on parental care. Firstly, aside from visit rates, we still know relatively little about how patterns of care are influenced by the ecological (predation, food distribution, etc.), environmental (temperature, weather, etc.), physiological (hormone, individual condition), and behavioral (foraging paradigm, pair stability, etc.) contexts of care. Secondly, while provisioning has been well-studied, other dimensions and stages of care—and how these influence each other—require far more attention. Both theory and empirical work suggests these can strongly impact carer behavior to the point that simply measuring one stage and dimension is insufficient.
We advocate both for more in-depth analyses of care behavior, and for raw parental care data from existing studies to be deposited alongside relevant publications, published as data papers, and shared with those interested in applying further analyses where feasible. Collectively, unpublished care data has the potential to greatly advance our understanding of how individuals provide for their offspring.
Our review is restricted to the quantification of carer behavior, but this is inextricably linked to the overall care paradigm, to the behavior of offspring, and to environmental variation. Many species exhibit two or more of the five main patterns of care (none, mother only, father only, biparental, cooperative), often within the same population (Persson and Öhrström, 1989; Webb et al., 1999). Additionally, offspring vary across species in their ability to influence care delivery, with consequences for investment levels, pre-natal effects and parent-offspring conflict (Hinde et al., 2010). Furthermore, most studies are time- and location-restricted, limiting our understanding of the effects of environmental variation. Integrating these complexities with a more detailed picture of carer behavior remains a key challenge for behavioral ecology.