Abstract
Adolescent depression is a major concern for public health and is associated with negative consequences and outcomes. Identifying adolescent characteristics that might relate to the risk for developing depression is crucial. This study investigated bidirectional associations between coping strategies and depressive symptoms over time. The participants were 1,341 secondary school students from the Netherlands (
Adolescent depression is a major concern for public health, as clinical and subclinical depression rates show a sharp increase in this critical developmental period (Kessler, Petukhova, Sampson, Zaslavsky, & Wittchen, 2012; Reef, Diamantopoulou, Van Meurs, Verhulst, & Van Der Ende, 2009). Most depressed adults experienced their first depressive disorder in adolescence, with the period between 13 and 18 years of age being the most critical time for the onset of depression (Kessler et al., 2005; Kim-Cohen et al., 2003). International epidemiological studies have shown that 2.7% of the 8- to 15-year-old adolescents and 7.5% of the 13- to18-year-old adolescents experience depression (Avenevoli, Swendsen, He, Burstein, & Merikangas, 2015; Merikangas et al., 2010). Although the depression rates are already concerning, a depressed mood and elevated symptoms of depression are not included in these rates and even more common in adolescence, with prevalence rates of 20% in the adolescent population (Meijer, Smit, Schoemaker, & Cuijpers, 2006).
There are several risk factors that are known to be associated with depression in adolescence, such as individual characteristics (e.g., perfectionism), social and family processes (e.g., abuse or bullying), or the presence of a psychiatric disorder. Gender has also proved to be an evident risk factor for the development of depression. Depressive symptoms among boys are relatively stable across adolescence, whereas girls report a significant increase in depressive symptoms from the age of 13 (Avenevoli et al., 2015; D. M. Costello, Swendsen, Rose, & Dierker, 2008).
Besides these, mostly stable risk factors, another factor that is related to the development of depressive symptoms is coping. Substantial individual differences exist in the ways in which adolescents handle events and stressors, and it has been suggested that the use of certain coping strategies is related to greater vulnerability of developing depressive symptoms in response to stressors (Garnefski, Legerstee, Kraaij, van Den Kommer, & Teerds, 2002). Therefore, most depression prevention programs include techniques to enhance coping skills and help adolescents cope with negative thoughts and feelings (Gillham et al., 2007). However, these programs show small effect sizes (Merry et al., 2012; Stice, Shaw, Bohon, Marti, & Rohde, 2009), and there is limited evidence as to whether the techniques included are effective for the prevention of depressive symptoms. To further our understanding of the association between coping and depressive symptoms, the present study investigated whether specific coping strategies are associated with levels of depressive symptoms over time among adolescents.
Coping Strategies
Coping can be defined as “conscious volitional efforts to regulate emotion, cognition, behavior, physiology, and the environment in response to stressful events or circumstances” (Compas, Connor-Smith, Saltzman, Thomsen, & Wadsworth, 2001). Several theorists have studied coping strategies of adolescents, arguing that adolescents who have difficulties with their emotional response to daily stressors experience more severe stress, which may ultimately lead to depression (Aldao, Nolen-Hoeksema, & Schweizer, 2010; Compas, Orosan, & Grant, 1993; Nolen-Hoeksema, Wisco, & Lyubomirsky, 2008). Active coping strategies (e.g., problem-solving) seem to provide a buffering effect, while strategies that involve disengagement from the source of stress put adolescents at risk for developing depressive symptoms (Aldao et al., 2010; Compas et al., 2001; Garnefski, Boon, & Kraaij, 2003; Wadsworth & Compas, 2002).
This was also confirmed in a recent meta-analysis among adolescents, conducted by Schäfer, Naumann, Holmes, Tuschen-Caffier, and Samson (2017), which indicated that avoidance and distraction were related to more depressive symptoms, whereas problem focusing and positive cognitive reframing were related to less depressive symptoms. Nevertheless, there are coping strategies that show mixed findings in relation to psychopathology. For example, social support was found to be unrelated to psychopathology according to a meta-analysis by Compas et al. (2017) but was associated with depressive symptoms in a meta-analysis by Rueger, Malecki, Pyun, Aycock, and Coyle (2016). This might be an example of a strategy that can be effective depending on an individual’s age. Social support might be adaptive for an early adolescent in the regulation of emotions, but it might be maladaptive for the mid-adolescent when it leads to co-rumination (Compas et al., 2017; Stone, Hankin, Gibb, & Abela, 2011). Therefore, researchers should allow for developmental effects, for example, by using longitudinal designs.
Longitudinal designs are especially important in this age group as strategies to cope with stress are not yet stable, and change might occur. Earlier research in the coping area showed that coping strategies increase with age in both number and variety. For example, older adolescents (14- to 18-year-olds) tend to use more emotion regulation strategies, like relaxation or distraction compared with younger adolescents (Donaldson, Prinstein, Danovsky, & Spirito, 2000). Moreover, developmental shifts ensure that adolescents start to use social partners to cope with stress, experience an increased ability to use cognitively complex processes, and change from more overt behavior (e.g., crying or yelling) to less overt strategies like distraction (Thompson & Goodman, 2010; Zimmer-Gembeck & Skinner, 2010).
Nonetheless, according to recent meta-analytic reviews most of these findings have come from cross-sectional studies (Aldao et al., 2010; Compas et al., 2017; Schäfer et al., 2017) and the existing literature has several limitations. First, most of the research between coping and depressive symptoms included adults and had small sample sizes. Second, there is limited research about the reversed effect of depressive symptoms on coping. This might be interesting as it is also possible that the level of symptoms is associated with an increased or decreased use of coping strategies, as this may interfere with the ability to use strategies when facing stress (Compas et al., 2017).
Furthermore, gender differences are often not included in analyses. The risk of depression differs when girls and boys reach adolescence (Avenevoli et al., 2015), and therefore the associations of coping strategies with depressive symptoms may also differ, depending on the moment in development at which coping strategies and depressive symptoms are measured (Carlson & Grant, 2008). Preliminary empirical evidence also suggests that boys and girls differ in coping strategies. According to some studies, girls are more likely to seek support and use more active approach coping strategies compared with boys (Eschenbeck, Kohlmann, & Lohaus, 2007; Frydenberg & Lewis, 1993). However, the findings are inconsistent and there are also studies that have shown opposite effects (Griffith, Dubow, & Ippolito, 2000; Hampel & Petermann, 2005). Allowing for gender differences is important as prevention and intervention programs can adapt to these differences and increase the effectivity.
To address these limitations, the present study provided a longitudinal and bidirectional perspective on the association between coping strategies and depressive symptoms using a sample of 1,341 adolescents who participated in a randomized controlled trial (RCT) examining the effectiveness of a universal depression prevention program. In addition, we allowed for gender differences in the associations between coping strategies and depressive symptoms. We also controlled for condition (intervention versus control), educational level, and ethnicity as they may interfere with depressive symptoms (E. J. Costello & Maughan, 2015). Although coping strategies can be classified in domains, such as adaptive and maladaptive, or primary control (acting directly on the source of stress), secondary control (strategies to adapt to the source of stress), and disengagement coping (strategies to take distance to the source of stress; Compas et al., 2017), the present study will examine single coping strategies in relation to depressive symptoms. This might yield a more nuanced picture of this relationship, which can be translated into practical guidelines that can be incorporated into prevention programs. Therefore, the present study investigates the following strategies: problem focusing, cognitive restructuring, avoidance, distraction, and seeking support in relation to depressive symptoms over time. Although the list of coping strategies is extensive, these are the most used strategies in literature and in coping questionnaires (Skinner, Edge, Altman, & Sherwood, 2003).
Present Study
We analyzed the associations between five single coping strategies and depressive symptoms using bidirectional cross-lagged models. Based on the findings of Garnefski et al. (2002; Garnefski et al., 2003) and Schäfer et al. (2017), we expected to find a negative association over time between problem focusing and positive cognitive reframing and depressive symptoms and we proposed positive associations over time between distraction and avoidance and depressive symptoms. While we expected to find negative associations over time between seeking support and depressive symptoms in this early adolescent phase, we were aware that this could change as adolescents reach the mid-adolescent phase (Compas et al., 2017; Rueger et al., 2016). In addition, we explored whether the reverse direction, depressive symptoms predicting an increased or greater use of specific coping strategies, was true. Given the limited literature on the association between depressive symptoms and coping, these analyses were explorative. Finally, due to inconsistencies regarding the role of gender in coping strategies, we performed additional analyses to explore the role of gender in the relationship between coping strategies and depressive symptoms.
Method
Participants and Procedure
The sample comprised 1,341 participants from 54 classes of nine secondary schools in the southern and middle part of the Netherlands who participated in the RCT of the universal depression prevention program “Op Volle Kracht” (OVK), which translates to “On Full Power” (see Tak, Lichtwarck-Aschoff, Gillham, Van Zundert, & Engels, 2016). This study was approved by the ethical committee of the research institute. There were no differences between the intervention and control group in levels of depressive symptoms and coping strategies on baseline and follow-ups. In addition, there was no effect of OVK on the levels of depressive symptoms or coping strategies at any time points. The participants’ mean age was 13.91 (SD = .55), and 47.5% were girls. Most participants were of Dutch origin (83.1%), and 7% were involved in prevocational education, 51.2% in higher general education, and 41.8% in pre university education.
Schools decided whether to participate in the study. All adolescents in the eighth grade of participating schools were eligible to participate (age range of 11–14). However, the pupils and their parents were informed about the study and were allowed to withdraw from the study at any point. Schools were randomly assigned to the intervention or control condition. Schools in the intervention condition integrated the OVK intervention in the school curriculum during mentor lessons. Schools in the control condition followed the schedule as usual.
Assessments were conducted by administering questionnaires at six time points, baseline (T1), and at 6-month (T2), 12-month (T3), 18-month (T4), 24-month (T5), and 30-month (T6) after baseline. Coping strategies were measured at five time points, and the questionnaire was not administered at T6 due to the length of the questionnaire and for practical reasons. Adolescents received incentives in the form of a gift voucher to complete each assessment. Retention rates were high across all assessments: 96.5% of the participants completed pre-intervention, 89.4% completed post-intervention, 89.3% completed the 6-month follow-up, 83.7% completed the 12-month follow-up, 77.4% completed the 18-month follow-up, and 84.5% completed the 24-month follow-up.
Measures
Coping strategies
We measured coping strategies with the Dutch version of the 54-itemed Children Coping Strategies Checklist–Revised (CCSC-R; de Boo & Wicherts, 2009). Adolescents had to rate whether they use a specific coping style when they face a problem on a 4-point scale ranging from 1 (almost never) to 4 (almost always). The questionnaire comprised five scales and 13 subscales. The five subscales were problem focusing that assessed cognitive decision-making, direct problem-solving, and seeking understanding (12 items; e.g., “Do something to make things better”); positive cognitive reframing that assessed positive and optimistic thinking and control (12 items; e.g., “Tell yourself that you can handle the problem”); distraction strategies that assessed the physical release of emotions and distracting actions (9 items; e.g., “Listen to music”); avoidance strategies that assessed avoidant actions such as repressing and wishful thinking (12 items; e.g., “Just forget about it”); and seeking support, which assessed seeking support to cope with actions and feelings (9 items; e.g., “Tell others how you feel about the problem”). Cronbach’s alpha ranged from .88 to .91 for problem focusing, from .87 to .91 for positive cognitive reframing, from .74 to .79 for distraction, from .73 to .85 for avoidance, and from .91 to .93 for seeking support across the different time points.
Depressive symptoms
We measured depressive symptoms with the Children’s Depression Inventory (CDI; Kovacs, 1985; Timbremont, Braet, & Roelofs, 2008). The CDI is a self-report questionnaire comprising 27 items, each consisting of three statements rated in severity from 0 to 2 (e.g., I don’t feel alone = 0, I often feel alone = 1, I always feel alone = 2). The sum of the scores of depressive symptoms ranged from 0 to 54. Cronbach’s alpha ranged from .84 to .91, indicating a high reliability of all assessments. Item 9 measured suicidal ideation; it was excluded from the questionnaire as this was beyond the scope of this research project.
Demographic variables
This study included the following demographic variables: age, gender, condition, educational level, and ethnicity. We included these variables as covariates as these factors are known to be associated with depressive symptoms and therefore might influence the results (E. J. Costello & Maughan, 2015). Ethnicity was measured by asking adolescents in which country they and their parents were born. When the adolescent or one of the parents was not born in the Netherlands, the adolescent was labeled as having a migration background. Although Tak et al. (2016) found no differences in results between adolescents in the experimental condition versus the control condition, we included this variable to control for possible effects.
Strategy of Analyses
Means, standard deviations, and correlations were computed for all outcome variables included in this study. The associations between coping strategies and depressive symptoms over time (five time points) were tested with five cross-lagged models for each of the five coping strategies. Because the number of parameters to be estimated in the model would increase rapidly by using items as indicators of the latent variables, with the consequence that power to detect important parameters will decrease (Yang & Dunson, 2010) and estimation problems will increase (Sass & Smith, 2006), we decided to use four parcels as indicators of the latent variable depressive symptoms and three parcels for each of the coping strategies. The items of each construct at T1 were allocated to three or four equivalent parts (parcels) according to the item-to-construct balance method (T. D. Little, Cunningham, Shahar, & Widaman, 2002). Parcels for the latent variables at T2 to T5 had identical indicators as the parcels at T1.
Adolescents were nested within 54 classes and classes were nested within nine schools. Correcting for clustering effects with three-level multilevel analysis was not possible because the number of schools was too low. To correct for school effects, eight dummy variables representing nine schools (see Cohen et al., 2003, pp. 303-307) were regressed on intercept and slope as the first step in our analyses. To correct for clustering effects within classes, the TYPE=COMPLEX procedure in Mplus version 7.2 was used (Muthén & Muthén, 1998-2015). The full information robust maximum likelihood (FIML) estimator was used to account for missing values. This estimator requires that missing values are missing at random (MAR), but there are no statistical MAR tests (Nakagawa, 2015). Therefore, we tested for the data mechanism missing completely at random (MCAR) with all outcome variables. Little’s MCAR test was significant, χ2(791) = 865.22, p = .034, indicating that missing values were not MCAR (R. J. Little, 1995). This result does not provide a definite answer for MAR so we used the FIML estimator under the assumption of MAR. Besides χ2(df) and p, two model fit measures were used: (a) the root mean square error of approximation (RMSEA; Byrne, 2005; Kaplan, 2000) and (b) the Bentler comparative fit index (CFI; Kaplan, 2000; Kline, 1998). RMSEA values lower than or equal to .05 are preferred, but values under .08 are acceptable while CFI values above .95 (.90) are indicative of a fair (acceptable) fit. Occasionally, unacceptable or untrusted parameter estimates are found during structural equation modeling (SEM) analyses. In that case, causes and consequences are examined and solutions are found. No unacceptable or untrusted parameter estimates were found in this study.
Cross-lagged analyses were performed for each coping strategy in combination with depressive symptoms. An overall model including depressive symptoms and the five coping strategies would have created a very complex model with too many parameters to be estimated (1,054) in relation to the sample size (1,341). For this reason, we decided to analyze five cross-lagged models separately. All latent variables of the cross-lagged model were regressed on the eight dummy variables for school effects and on the control variables age, condition, education level, and ethnicity. The parameters of the cross-lagged model were estimated freely (in the baseline model). To test equality of regression weights of cross paths, we first tested a model by constraining the cross paths between every consecutive time points to be equal. If this constrained model showed a significant difference with the baseline model (a significant increase of chi-square), post hoc χ2 difference tests were used to test which pair(s) of cross paths were significantly different. Because a robust ML estimator was used (due to the COMPLEX procedure), the χ2 values were robust and had to be rescaled to unbiased χ2 values before using the χ2 difference tests (Satorra & Bentler, 2010). This test statistic is known as the Sattora-Bentler (SB)-scaled χ2 difference test. The robust ML estimator is also robust against non-normality of the variables (see depressive symptoms in Table A1 in the Appendix); parameters in the cross-lagged models are estimated with robust standard errors.
To explore gender differences in regression weights of cross-lagged relations, we used multiple group analysis. First, a baseline model was estimated without constraints between boys and girls. The fit of this model was the baseline χ2. In the second model, the regression weights of cross paths were constrained to be equal across gender, and the χ2 of this model was compared with the baseline. A significant difference in χ2 is an indication that one or more paths are different for boys and girls.
Results
Descriptive statistics (mean, SD, and correlations for each time point) of depressive symptoms and coping strategies are presented in Tables A1 and A2 in the Appendix. The mean for depressive symptoms was 7.55, which is far below the cutoff of 13 that is used to screen for depressive symptoms (Timbremont, Braet, & Dreessen, 2004), indicating that this is a relative healthy sample that is comparable with other universal prevention studies in non-high-risk population samples (e.g., Challen, Machin, & Gillham, 2014; Horowitz & Garber, 2006). Girls reported a greater use of seeking support on all time points compared with boys. In addition, they reported more problem focusing at T3, T4, and T5, positive cognitive reframing at T5, and avoidance at T1 and T5. Regarding depressive symptoms, girls reported higher levels of depressive symptoms at T1 but lower levels of symptoms at T4 compared with boys. No baseline differences in study variables were found for age, condition, educational level, or ethnicity.
Prior to the final cross-lagged analyses, we tested measurement equivalence for each of the six latent variables over time and across gender. To compare unstandardized regression coefficients of the cross-lagged models, configural equivalence and weak metric equivalence are required (Guenole & Brown, 2014; Keith, 2014). Configural equivalence means that the number of factors and pattern of loadings (but not the strengths of the loadings) are the same over time or across gender. Weak metric equivalence is an additional requirement, and means equal metrics (= equal scale intervals) of a latent variable over time or across gender. Factor loadings define the metric of measurement and represent the strength of the relationship between a factor (latent variable) and indicators (parcels in our case). When factor loadings of a latent variable are equal over time, or across gender, the unit of measurement is also equal and relationships with other (latent) variables can be studied (Wang, Chen, Dai, & Richardson, 2018). This allows the comparisons of unstandardized regression coefficients and covariances over time and across groups (Steenkamp & Baumgartner, 1998).
If the fit measures of the configural model are acceptable, configural equivalence is supported. For weak metric equivalence, the factor loadings of the identical parcels over time are constrained to be equal. The fit of this constrained model is compared with the baseline or configural model. Strong metric equivalence was tested by additionally constraining the corresponding intercepts to be equal over time. The fit of this model was compared with the factor loadings constrained model. A decrease of CFI < .01 and an increase of RMSEA < .015 supported the weak or strong metric equivalence of the longitudinal factor model (Chen, Curran, Bollen, Kirby, & Paxton, 2008; Cheung & Rensvold, 2002). The results are presented in Table A3 in the Appendix. For weak metric equivalence, the decrease of CFI for all models was < .01 and the increase of RMSEA was just a little above the threshold of .015, only for gender in relation to seeking support. Weak measurement invariance is supported for all six latent variables and (with one small exception) for gender. For strong measurement equivalence (allowing comparisons of latent means over time and across groups), a decrease of CFI by more than .01 was found for gender in relation to distraction, avoidance, and seeking support. An increase of RMSEA by more than .015 was found for positive cognitive reframing and gender of positive cognitive reframing, for gender of distraction, for avoidance and gender of avoidance, and for seeking support. This means that strong metric equivalence was not always met for some latent variables in combination with gender. As already noted, for testing paths of the cross-lagged models, weak measurement invariance is sufficient. In the last four columns of Table A3, the range of the factor loadings across five time points are presented, including their means and standard deviations. As can be seen in Table A3, all factor loadings are high indicating that the parcels represent their factor very well for each time point.
Cross-Lagged Models of Coping Strategies With Depressive Symptoms
Problem focusing
The fit of the cross-lagged model was good; χ2(769) = 1426.38, p < .001, CFI = .978, RMSEA = .025, 90% confidence interval (CI): [.023, .027], p(RMSEA) = 1.000. The model chi-square test was significant, indicating a poor model fit. However, in large samples chi-square is almost always significant and global fit measures can be used as an alternative. CFI was > .95 and RMSEA was < .05 and, according to the norms used for global fit measures, these values indicate a good model fit. Problem focusing did not predict depressive symptoms over time. Depressive symptoms had a significantly negative association with problem focusing over time, from T1 to T2, T2 to T3, and T3 to T4, indicating that elevated depressive symptoms were associated with lower levels of problem focusing 6 months later. Figure 1 presents the completely standardized parameter estimates. Using the SB-scaled χ2 difference test, we first tested whether the four pairs of cross paths were significantly different by constraining the cross paths between T1 and T2 to be equal, as well as between T2 and T3, between T3 and T4, and between T4 and T5. This test was significant, χ2(4) = 26.08, p < .001. Post hoc difference tests showed that the cross paths were different between T1 and T2, χ2(1) = 8.03, p = .005, between T2 and T3, χ2(1) = 27.05, p < .001, and between T3 and T4, χ2(1) = 3.99, p = .046. This means that the negative associations between depressive symptoms and problem focusing were significantly greater compared with the associations between problem focusing and depressive symptoms at the same time points. The SB-scaled χ2 difference test showed no significant difference in cross paths between boys and girls, χ2(8) = 24.26, p = .075.

Cross-lagged model of depressive symptoms and problem focusing.
Positive cognitive reframing
The fit of the cross-lagged model was χ2(769) = 1521.13, p < .001, CFI = .975, RMSEA = .027, 90% CI: [.025, .029], p(RMSEA) = 1.000. The global fit measures CFI and RMSEA indicated a good fit. Positive cognitive reframing did not predict depressive symptoms over time. Depressive symptoms were significantly negatively associated with positive cognitive reframing over time, from T1 to T2 and from T2 to T3, indicating that higher depressive symptoms were associated with lower levels of positive cognitive reframing 6 months later (see Figure 2). With the SB-scaled χ2 difference test, we found that one or more of the four pairs of cross paths were significantly different, χ2(4) = 11.01, p = .026. Post hoc difference tests showed that the cross paths were different between T1 and T2, χ2(1) = 3.95, p = .047, and between T2 and T3, χ2(1) = 4.41, p = .036. This indicates that the negative associations between depressive symptoms and positive cognitive reframing were significantly greater compared with the associations between positive cognitive reframing and depressive symptoms at the same time points. The SB-scaled χ2 difference test showed no significant difference in cross paths between boys and girls, χ2(8) = 7.31, p = .504.

Cross-lagged model of depressive symptoms and positive cognitive reframing.
Distraction
The fit of the cross-lagged model was χ2(769) = 1588.13, p < .001, CFI = .965, RMSEA = .028, 90% CI: [.026, .030], p(RMSEA) = 1.000. The global fit measures indicated a good model fit. Distraction did not predict depressive symptoms over time. Depressive symptoms had a significant negative association with distraction over time, from T3 to T4 (see Figure 3). Using the SB-scaled χ2 difference test, we found that cross paths were not significantly different, χ2(4) = 5.26, p = .262. The SB-scaled χ2 difference test showed a significant difference in cross paths between boys and girls, χ2(8) = 16.61, p = .034. Post hoc testing showed that cross path Dep T2 to DS T3 was significantly different for boys and girls, χ2(1) = 3.91, p = .048, with B = −.01 and p = .955 for boys and B = .19 and p = .004 for girls. Girls’ elevated depressive symptoms at T2 predicted an increase in distraction strategies at T3, but this association was not found for boys.

Cross-lagged model of depressive symptoms and distraction.
Avoidance
The fit of the cross-lagged model was χ2(769) = 1491.348, p < .001, CFI = .971, RMSEA = .026, 90% CI: [.024, .028], p(RMSEA) = 1.000. The fit measures CFI and RMSEA indicated a good model fit. No significant cross-associations were found between avoidance and depressive symptoms (see Figure 4). The SB-scaled χ2 difference test showed a significant difference in cross paths between boys and girls, χ2(8) = 25.55, p = .001. Post hoc testing revealed that cross path Dep T2 to AV T3 was significantly different for boys and girls, χ2(1) = 10.31, p = .001, with B = −.02 and p = .838 for boys and B = .28 and p = .011 for girls; cross path Dep T3 to AV T4 was significantly different for boys and girls, χ2(1) = 7.92, p = .005, with B = −.29 and p = .005 for boys and B = .09 and p = .394 for girls; and cross path Dep T4 to AV T5 was significantly different for boys and girls, χ2(1) = 6.17, p = .013, with B = .02 and p = .824 for boys and B = .24 and p = .001 for girls. The results suggest that elevated depressive symptoms in boys at T3 were associated with lower levels of avoidance strategies 6 months later. Regarding girls, elevated depressive symptoms at T2 and T4 were associated with an increase in avoidance strategies 6 months later.

Cross-lagged model of depressive symptoms and avoidance.
Seeking support
The fit of the cross-lagged model was χ2(769) = 1566.92, p < .001, CFI = .975, RMSEA = .028, 90% CI: [.026, .030], p(RMSEA) = 1.000. The fit measures CFI and RMSEA indicated a good model fit. Seeking support at T2 showed significantly negative association with depressive symptoms at T3 and vice versa. Depressive symptoms showed a significantly negative associations with seeking support over time, from T1 to T2 and from T2 to T3, indicating that elevated depressive symptoms at T1 and T2 were associated with lower levels of seeking support 6 months later (see Figure 5). With the SB-scaled χ2 difference test, we found that cross paths were significantly different, χ2(4) = 15.17, p = .004. Post hoc difference tests showed that the cross paths between T2 and T3, χ2(1) = 9.44, p = .002, were significantly different. This indicates that the negative associations between depressive symptoms and seeking support were significantly greater at 6 months compared with the associations between seeking support and depressive symptoms on the same time points. The SB-scaled χ2 difference test showed no significant difference in cross paths between boys and girls, χ2(8) = 10.91, p = .207.

Cross-lagged model of depressive symptoms and seeking support.
Discussion
The aim of the present study was to investigate the bidirectional associations between coping strategies and depressive symptoms over time. It was expected that adolescents high in problem focusing, positive cognitive reframing, and seeking support would experience fewer depressive symptoms. Likewise, we expected that adolescents high in distraction and avoidance would experience more depressive symptoms. The analyses regarding the reversed effect of depressive symptoms on coping strategies were exploratory. In addition, we explored the role of gender differences in these relations.
In contrast to the hypotheses, no robust significant associations between coping strategies and depressive symptoms emerged. Of the 40 associations tested, we found that coping affected depressive symptoms over time only in one association. However, there was an indication for the reverse relationship. When adolescents experienced elevated depressive symptoms, they used fewer adaptive techniques to target stress in practical ways, using, for example, problem-solving and cognitive decision-making. In addition, depressive symptoms had a negative association with positive cognitive reframing and seeking support but only at the first two time points. One cross path showed a significant negative association between depressive symptoms and distraction. Although this effect was not consistent over time, it was contrary to expectations and in contrast with research that classified distraction as maladaptive (Aldao et al., 2010).
Significant gender differences were found in depressive symptoms and avoidance across three time points. For boys, a higher score of depressive symptoms was related with lower levels of avoidance strategies. For girls, a higher score on depressive symptoms was related to an increase in avoidance and distraction strategies. However, significant findings among gender differences were neither consistent nor robust. Future studies should clarify these findings to draw further conclusions.
The prospective one-way association between depressive symptoms and the use of problem focusing and positive cognitive reframing was inconsistent with our hypotheses, previous studies (e.g., Aldao & Nolen-Hoeksema, 2010; Garnefski et al., 2003), and models in which bidirectional associations between cognitive variables and depressive symptoms are integrated, such as the transactional framework (Hankin & Abramson, 2001). In these models, maladaptive cognitive schemas predicted depressive symptoms and vice versa. Inconsistencies with previous studies could be caused by differences in the developmental period that were captured or the timing and spacing of measurement intervals (Collins & Graham, 2002); however, there could be other reasons for this inconsistency.
One explanation could be that adolescents with elevated depressive symptoms are less able to use active problem-solving and they tend to revert to other strategies, which suggests the presence of a third factor. This argument is also supported by the meta-analysis by Compas et al. (2017). The researchers stated that due to the lack of longitudinal designs, little is known about the actual direction of coping with psychopathology symptoms. While researchers in general might be more interested in the effect of coping strategies on symptoms for practical reasons, Compas et al. (2017) argued that high levels of symptoms may impede the development of these skills or the ability to use them effectively. It would be interesting to investigate the initial level of symptoms or specific characteristics of depressive symptoms in relation to coping strategies.
Following this argument, another explanation might concern methodological issues. Most studies included cross-sectional designs; hence, the stability over time might have been missed. Multiple subsequent time points should be included in studies to generate substantial conclusions regarding the associations between coping strategies and depressive symptoms (Sameroff & Mackenzie, 2003). The current finding regarding significant associations between depressive symptoms and problem focusing was consistent across four measurements, which increases the reliability of our conclusion. Moreover, the associations between depressive symptoms and coping strategies were analyzed by including the within time correlations and stability in depressive symptoms and coping strategies subscales over time. Inclusion of these factors is important as this provides a more accurate image of the complex associations between depressive symptoms and coping strategies. Nevertheless, current findings need to be replicated before firm conclusions can be made.
Although the findings in the coping literature provide evidence for gender-specific coping preferences, little is known about the direction and effect of coping strategies on psychopathology symptoms. Adult studies have shown that females often report using seeking support, rumination, positive self-talk, and problem-solving while males tend to use avoidant strategies or passivity (Eschenbeck et al., 2007; Tamres, Janicki, & Helgeson, 2002). The differences between men and women regarding the efficacy of specific strategies have not been clarified yet. This study showed that for boys, elevated depressive symptoms were associated with lower levels of avoidance strategies, whereas in girls, elevated depressive symptoms were associated with higher levels of avoidance strategies. This is not in accordance with the findings in a recent meta-analysis by Schäfer et al. (2017), in which it was stated that maladaptive strategies, such as avoidance, were more strongly associated with depressive symptoms among men compared with women. Despite the fact that literature on this subject in adolescence is scarce, adult literature suggests that there are both differences and similarities in the relationship between coping and depressive symptoms in men and women. Holahan, Moos, Holahan, Brennan, and Schutte (2005) found that for women, avoidant coping was linked to depressive symptoms both directly as well as indirectly through life stressors. For men, life stressors completely explained the relationship between avoidant coping and depressive symptoms. Future research should examine the effect of gender on the relation between depressive symptoms and avoidance more closely to replicate our findings and to examine the consequences of this difference.
Strengths and Limitations
This study has some strengths and limitations that provide direction for future research. We used longitudinal data from a large general community sample of adolescents, which facilitates the generalization of the results. In addition, a sample drawn from the general population reduces biases known in clinical research, such as comorbidity, severity of symptoms, and treatment seeking (Goodman et al., 1997). However, most adolescents experienced no depressive symptoms and a small group experienced low to mild depressive symptoms which might have limited the results because of the lack of variance and attenuation of the associations that were examined. In addition, since the study included relatively healthy adolescents, the associations between coping strategies and depressive symptoms might be different for adolescents who experience depression or a higher level of depressive symptoms. Furthermore, the study was based on self-reported data and could therefore be vulnerable to biases. Also, we did not include an assessment regarding sources of the stress that the adolescents experienced, which is an important limitation because some strategies might be adaptive in one situation but not adaptive in another (e.g., problem focusing is less adaptive when one needs to cope with the loss of a loved one). Therefore, we recommend that future studies should include more information about the chronicity or controllability of the stressors.
Regarding the statistical/methodological strengths and limitations, we made use of item parceling which has several advantages and disadvantages (see T. D. Little et al., 2002). Regarding this discussion, researchers concluded that parceling is effective when items within a parcel measure the same construct and this was the case in the current study (Bandalos & Finney, 2001). A statistical limitation is that the large autocorrelations of the latent variables over time could decrease the chances for finding cross-lagged effects. Controlling for past levels of predictors will dramatically reduce the magnitude of cross-lagged effects, especially if autoregressive paths (stability paths) are strong (Adachi & Willoughby, 2015). However, small effect sizes for cross-lagged effects (standardized regression weights) in longitudinal models are the rule rather than the exception. Current guidelines to judge effect sizes are mostly based on cross-sectional research (were stability over time is not controlled). The small effect sizes in longitudinal studies must be assessed according to standards based on longitudinal research. In models with high stability paths, as in this study, small cross-lagged significant effects are not trivial (Adachi & Willoughby, 2015).
Implications for Practice and Future Research
This study has implications for practice and provides interesting new directions for future research. First, the finding that depressive symptoms only seem to interfere with the ability to use problem-focused coping strategies rather than the other way around when facing stress implies that altering coping strategies through depression prevention would not necessarily decrease depressive symptoms. This is in line with a review of Stice et al. (2009) who stated that the content of depression prevention programs was weakly or not associated with a decrease in depressive symptoms. Although more research is necessary to draw conclusions about the effectiveness or ineffectiveness of targeting coping strategies in prevention or intervention, practitioners and researchers should investigate other concepts in the prevention of depression, such as examining techniques that are more behavioral rather than cognitive.
Behavioral activation technique, for instance, focuses on mood monitoring and daily activities to increase positive activities and positive interactions with the environment. Behavioral activation has already proven to be effective in depression treatments and would be a cost-effective solution for prevention programs targeting adolescents with elevated depressive symptoms (Cuijpers et al., 2007). To investigate the effect of activation on adolescents’ mood, the experience sample method (ESM) is an appropriate methodology. ESM is a daily diary method that requires the participants are to report their thoughts and feelings in the moment at multiple time intervals, using, for instance, smart phone mobile apps. In addition to investigating fluctuations over short-time rather than long-term developments (de Haan-Rietdijk, Voelkle, Keijsers, & Hamaker, 2017), technological developments have made it possible to give participants personalized feedback and thus give insight into their mood and activities (Myin-Germeys, Klippel, Steinhart, & Reininghaus, 2016). Research has shown that this is a promising and feasible solution to increase positive affect in people with depressive symptoms (Van Roekel et al., 2017).
Second, future research should include longitudinal and dynamic methodology to investigate the concept of coping in relation to depressive symptoms to draw more substantial conclusions. Facing these demands, a promising solution would be to utilize person-centered and data-driven methods to identify profiles of adolescents who use similar patterns of coping strategies (Aldao, 2013). The few studies that have used a person-centered approach have shown that this approach can be successful in extracting coping profiles (Dixon-Gordon, Aldao, & De Los Reyes, 2015; Lougheed & Hollenstein, 2012; Van den Heuvel, Stikkelbroek, Bodden, & van Baar, 2020). Moreover, this provides the opportunity to investigate divergence in developmental trajectories of coping strategies as well as the development of depressive symptoms. Regarding coping, it would be interesting to investigate whether strategies are stable across adolescence or transform with development. Prevention programs could benefit from information on the course of depressive symptoms by adapting programs to the needs of subgroups or to establish the ideal starting point for prevention. Innovative approaches in future research are required to improve our prevention programs focused on reducing depression among adolescents.
Footnotes
Appendix
Unstandardized Results of Cross-Lagged Analyses Between Coping Strategies and Depressive Symptoms.
| Model | Path | B | SE | p |
|---|---|---|---|---|
| 1. Problem Focusing | PF1–Dep2 | −.01 | −.01 | .40 |
| Dep1–PF2 | −.23 | .08 | .00 | |
| PF2–Dep3 | −.01 | .01 | .25 | |
| Dep2–PF3 | −.29 | .08 | .00 | |
| PF3–Dep4 | −.01 | .01 | .38 | |
| Dep3–PF4 | −.24 | .11 | .03 | |
| PF4–Dep5 | −.00 | .01 | .92 | |
| Dep4–PF5 | −.14 | .10 | .17 | |
| 2. Positive Cognitive Reframing | PCF1–Dep2 | −.01 | .01 | .22 |
| Dep1–PCF2 | −.20 | .09 | .02 | |
| PCF2–Dep3 | −.02 | .01 | .08 | |
| Dep2–PCF3 | −.19 | .09 | .03 | |
| PCF3–Dep4 | −.01 | .02 | .58 | |
| Dep3–PCF4 | −.16 | .09 | .08 | |
| PCF4–Dep5 | −.00 | .01 | .75 | |
| Dep4–PCF5 | −.04 | .08 | .61 | |
| 3. Distraction | DS1–Dep2 | .00 | .01 | .92 |
| Dep1–DS2 | −.03 | .06 | .60 | |
| DS2–Dep3 | .01 | .02 | .47 | |
| Dep2–DS3 | .07 | .06 | .27 | |
| DS3–Dep4 | .01 | .02 | .59 | |
| Dep3–DS4 | −.17 | .08 | .04 | |
| DS4–Dep5 | .03 | .02 | .18 | |
| Dep4–DS5 | .05 | .07 | .51 | |
| 4. Avoidance | AV1–Dep2 | .00 | .01 | .90 |
| Dep1–AV2 | −.07 | .09 | .40 | |
| AV2–Dep3 | .02 | .02 | .43 | |
| Dep2–AV3 | .11 | .07 | .09 | |
| AV3–Dep4 | .02 | .02 | .33 | |
| Dep3–AV4 | −.12 | .08 | .10 | |
| AV4–Dep5 | .02 | .02 | .22 | |
| Dep4–AV5 | .12 | .08 | .11 | |
| 5. Seeking support | SS1–Dep2 | −.01 | .01 | .49 |
| Dep1–SS2 | −.28 | .11 | .01 | |
| SS2–Dep3 | −.03 | .01 | .00 | |
| Dep2–SS3 | −.24 | .09 | .01 | |
| SS3–Dep4 | .00 | .01 | .80 | |
| Dep3–SS4 | −.17 | .10 | .09 | |
| SS4–Dep5 | −.01 | .01 | .16 | |
| Dep4–SS5 | −.13 | .11 | .22 |
Note. Dep = depressive symptoms, PF = problem focusing, PCF = positive cognitive reframing, AV = avoidance, DS = distraction, SS = seeking support.
Declaration of Conflicting Interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: The municipality of Oss, The Netherlands, provided funding for this study. ZonMw, The Netherlands Organization for Health Research and Development, funded the original study of
.
Ethical Approval
The current submission does not overlap with any other published, in press, or in preparation articles, books, or proceedings and has not been posted on a website. Our research is not under consideration elsewhere and has been conducted in accordance with ethical standards in the field.
Informed Consent
Passive parental consent was obtained, and all students included in the study participated voluntarily.
Data
The dataset generated analyzed during the current study are available from the corresponding author on reasonable request.
