ORIGINAL RESEARCH article

Front. Psychiatry, 20 May 2021

Sec. Public Mental Health

Volume 12 - 2021 | https://doi.org/10.3389/fpsyt.2021.676256

Mental Health Symptoms Unexpectedly Increased in Students Aged 11–19 Years During the 3.5 Years After the 2016 Fort McMurray Wildfire: Findings From 9,376 Survey Responses

  • 1. Department of Computing Science, University of Alberta, Edmonton, AB, Canada

  • 2. Department of Psychiatry, University of Alberta, Edmonton, AB, Canada

  • 3. Department of Accounting and Business Analytics, Alberta School of Business, University of Alberta, Edmonton, AB, Canada

  • 4. Faculty of Social Work, University of Calgary, Calgary, AB, Canada

  • 5. Department of Sociology and Anthropology, Mount Royal University, Calgary, AB, Canada

  • 6. Faculty of Health Sciences, University of Lethbridge, Lethbridge, AB, Canada

  • 7. Fort McMurray Catholic School District, Fort McMurray, AB, Canada

  • 8. Fort McMurray Public School District, Fort McMurray, AB, Canada

Article metrics

View details

32

Citations

7,5k

Views

2,4k

Downloads

Abstract

In Fort McMurray, Alberta, Canada, the wildfire of May 2016 forced the population of 88,000 to rapidly evacuate in a traumatic and chaotic manner. Ten percentage of the homes in the city were destroyed, and many more structures were damaged. Since youth are particularly vulnerable to negative effects of natural disasters, we examined possible long-term psychological impacts. To assess this, we partnered with Fort McMurray Public and Catholic Schools, who surveyed Grade 7–12 students (aged 11–19) in November 2017, 2018, and 2019—i.e., at 1.5, 2.5, and 3.5 years after the wildfire. The survey included validated measurement scales for post-traumatic stress disorder (PTSD), depression, anxiety, drug use, alcohol use, tobacco use, quality of life, self-esteem, and resilience. Data analysis was done on large-scale anonymous surveys including 3,070 samples in 2017; 3,265 samples in 2018; and 3,041 samples in 2019. The results were unexpected and showed that all mental health symptoms increased from 2017 to 2019, with the exception of tobacco use. Consistent with this pattern, self-esteem and quality of life scores decreased. Resilience scores did not change significantly. Thus, mental health measures worsened, in contrast to our initial hypothesis that they would improve over time. Of note, we observed higher levels of mental health distress among older students, in females compared to male students, and in individuals with a minority gender identity, including transgender and gender-non-conforming individuals. These findings demonstrate that deleterious mental health effects can persist in youth for years following a wildfire disaster. This highlights the need for multi-year mental health support programs for youth in post-disaster situations. The indication that multi-year, post-disaster support is warranted is relatively novel, although not unknown. There is a need to systematically investigate factors associated with youth recovery following a wildfire disaster, as well as efficacy of psychosocial strategies during later phases of disaster recovery relative to early post-disaster interventions.

Introduction

In May 2016, a large wildfire affected Fort McMurray, Alberta, Canada and the surrounding area. Called “The Beast” in the popular media (1), the wildfire necessitated the evacuation of the entire city on May 3, 2016, and over 88,000 residents were displaced for several weeks. In addition to damaging community structures and infrastructure, the fire destroyed 10% of the homes in Fort McMurray, leaving many individuals homeless. The fire burned 590,000 hectares of land before it was brought under control on July 4, 2016. Insurance costs for the damages were estimated at $3.6 billion by the Insurance Bureau of Canada, which made the Fort McMurray wildfire the most expensive insured catastrophe in Canadian history (2). Many individuals were left jobless due to damage and closure of local businesses. Social, emotional, and psychological difficulties also affected the community, as is typical after severe disasters (3, 4).

Disasters tend to impact children and youth in particular (510) given developmental vulnerabilities (cognitive, emotional, social, and physiological) associated with childhood and adolescence, including the need to rely on others for support (11). Our group has previously examined Fort McMurray school mental health surveys completed by Grade 7–12 students in November 2017 (12, 13). Brown et al. (12) reported elevated mental health symptoms for PTSD and depression compared to a control population from the same province which had not experienced a natural disaster. As reported in (13), individuals who were more personally impacted by the 2016 wildfire, such as having their home destroyed, exhibited greater symptoms of PTSD, depression, anxiety, and alcohol and/or substance misuse. These findings were consistent with previous findings of altered mental health in youth following natural disasters (9, 14).

Previous studies have also reported links between wildfires and adverse mental health effects in adults and children, including increased symptoms of PTSD (1524), increased symptoms of depression and anxiety (17, 19, 2328), mental health difficulties (18, 29, 30), reduction in health and well-being (31), increased consumption of anxiolytics-hypnotics (32), and decreased capacity to cope with adversity (33). The combined adverse impacts of wildfire disasters on individuals' mental health and coping abilities, particularly in youth, are important because coping plays a major role in post-disaster recovery, and lack thereof can be detrimental to long-term mental health outcomes [also see (34)]. More generally, non-wildfire disasters, including floods, earthquakes, and tsunamis, have also been linked to adverse effects on mental health in adults and children in the form of increased symptoms or incidences of PTSD, major depressive disorder, suicidality, generalised anxiety disorder, and substance use disorder [see (8, 9, 3539)]. Displacement caused by natural disasters (e.g., short- and long-term evacuation, homelessness) imposes additional physical, emotional, and psychosocial stress on families, which may then adversely affect longer-term development of children and youth (40).

Longitudinal studies, most of which have focused on adults, have reported long-term negative impacts on mental health from natural disasters (17, 25, 4146). The pattern of long-term effects can be complex, with individuals exhibiting better recovery from some symptoms, such as depression or anxiety, but lesser recovery for others such as PTSD symptoms (47). Relatively few longitudinal studies have focused on mental health outcomes among youth following natural disaster, although studies of children and youth following Hurricane Katrina have shown long term impacts on mental health (48), with varying recovery trajectories for different groups (49, 50). We are aware of only three population studies that examined wildfire impacts on youth mental health, including one study focusing on PTSD and depression (19) and our previous Fort McMurray studies (12, 13).

In the present study, our hypothesis was that mental health symptoms among Fort McMurray youth would improve with time following the November 2017 survey. This hypothesis was based on theoretical work on trauma recovery, which consistently emphasises the role of time in the recovery process (4, 51, 52). To test this hypothesis, we examined mental health survey data collected by Fort McMurray school boards in Grades 7–12 in November across three consecutive years, including 2017, 2018, and 2019. This repeated testing was conducted 1.5, 2.5, and 3.5 years after the 2016 wildfire (All data were collected prior to the COVID-19 pandemic). We hypothesised specifically that symptoms of PTSD, depression, anxiety, and alcohol/substance use would steadily improve from 2017 to 2018 to 2019.

Materials and Methods

Overview and Ethical Considerations

Information collected from students included questions on demographics, mental health, resilience, and personal exposure to and direct impacts of the wildfire. Measurement instruments were selected by the school systems, informed by the relevant literature and advice from the University of Alberta research team. Written letters were sent to parents and guardians to inform them of the survey 2 weeks prior. Parents and guardians could have their child(ren) not participate in the survey, and students also independently had the option to participate or not in the survey as explained at the start of each survey session (see details below and Appendix A: Survey Description Script in the Supplementary Material). Survey data collection was intentionally anonymous, and participants were not asked for their names nor any other identifying information. The study design was approved by the University of Alberta's Health Research Ethics Board (ethics protocol number Pro00072669 approved June 26, 2017). This paper reports on findings from the anonymous survey data collected from both school boards.

Mental health surveys were conducted by the two school boards in Fort McMurray—Fort McMurray Public Schools and Fort McMurray Catholic Schools (henceforth, “Schools”). They asked all students in Grades 7–12 to complete the surveys in November 2017 (18 months after the 2016 wildfire), November 2018 (30 months post-wildfire), and November 2019 (42 months post-wildfire) (All data were collected before the COVID-19 pandemic). The survey was administered during regular class time as part of the standard curriculum to evaluate the support programs the Schools had put in place following the wildfire (see Appendix B: Mental Health Support Programs, in the Supplementary Material). The Schools determined that surveys would be done in November 2017, 2018, and 2019 as the month of November worked best given various logistical and staff capacity considerations.

Fort McMurray Public Schools and Fort McMurray Catholic Schools (“Schools”) administered all aspects of survey data collection, including participant consent, in accordance with their standard procedures and policies. The Schools asked researchers from the University of Alberta for assistance in designing the survey and analysing the anonymous dataset.

Survey Questionnaires

The survey included 10 questionnaires (see Table 1 for additional details):

  • Demographics Questionnaire (Demographics, 7 questions)—a custom-designed questionnaire assessing age, gender, and the student's grade and school.

  • Impact of Fire Questionnaire (IOF, 6 questions)—a custom-designed questionnaire to assess the impact of the 2016 wildfire on the student.

  • Child PTSD Symptom Scale (CPSS, 19 questions)—used to assess symptoms of post-traumatic stress disorder (PTSD) (53); total CPSS score ranges from 0 to 51.

  • Patient Health Questionnaire, Adolescent version (PHQ-A, 11 questions)—used to assess symptoms of depression and suicidality (54, 55); total PHQ-A score ranges from 0 to 27.

  • Hospital Anxiety and Depression Scale (HADS, 7 questions, anxiety-related questions only)—used to assess symptoms of anxiety (56); total HADS score ranges from 0 to 21.

  • CRAFFT Questionnaire (CRAFFT, 9 questions)—used to assess symptoms of alcohol and substance misuse (57, 58); total CRAFFT scores ranges from 0 to 6.

  • Tobacco Use Questionnaire (2 questions) – two questions on tobacco use: “During the past month, did you smoke tobacco products?” and “During the past month, did you use smokeless tobacco products?”

  • Rosenberg Self-Esteem Scale (Rosenberg, 10 questions)—used to assess self-esteem (59); total score ranges from 0 to 30 with higher scores indicating higher self-esteem.

  • Kidscreen Questionnaire (Kidscreen-10, 11 questions)—used to assess quality of life (60); total Kidscreen-10 score ranges from 0 to 44 with higher scores indicating better quality of life.

  • Child and Youth Resilience Measure (CYRM-12, 12 questions)—used to assess resilience to adverse experience or trauma (61); total CYRM-12 score ranges from 12 to 60.

Table 1

QuestionsAnswer choices
Demographics questionnaire
1Are you at school right now, while you are taking the survey?Yes, No
2Are you a student?Yes, No
3What gender do you identify with?Female, Male, Other, Prefer not to say
4What is your age?10 years or less, 11 years, 12 years, 13 years, 14 years, 15 years, 16 years, 17 years, 18 years, 19 years, 20 years or more
5What school are you in currently?Select from a list of all Ft McMurray schools with any classes in grades 7–12
6What grade are you in currently?7, 8, 9, 10, 11, 12, other
7What school were you in for grade 6?Select from a list of all Ft McMurray schools with grade 6
Impact of fire questionnaire
1Were in you or near Fort McMurray during any part of the 2016 wildfire?Yes, No
2Did you evacuate because of the fire?Yes, No
3Was your home destroyed by the fire?Yes, No
4Did you see the fire in person?Yes, No
5What school are you in?Select from a list of all Ft McMurray schools with any classes in grades 7–12
6What grade are you in?7, 8, 9, 10, 11, 12, other
PATIENT HEALTH QUESTIONNAIRE (PHQ-A, DEPRESSION SYMPTOMS)
Over the past 2 weeks, how often have you been bothered by any of the following problems?
1Feeling down, depressed, irritable or hopelessNot at all, Several days, More than half the days, Nearly every day
2Little interest or pleasure in doing things?Not at all, Several days, More than half the days, Nearly every day
3Trouble falling or staying asleep, or sleeping too muchNot at all, Several days, More than half the days, Nearly every day
4Poor appetite, weight loss, or overeating?Not at all, Several days, More than half the days, Nearly every day
5Feeling tired, or having little energy?Not at all, Several days, More than half the days, Nearly every day
6Feeling bad about yourself-or that you are a failure or that you have let yourself or your family downNot at all, Several days, More than half the days, Nearly every day
7Trouble concentrating on things, such as school work, reading or watching televisionNot at all, Several days, More than half the days, Nearly every day
8Moving or speaking so slowly that other people could have noticed. Or the opposite-being so figety or restless that you have been moving around a lot more than usualNot at all, Several days, More than half the days, Nearly every day
9Thoughts that you would be better off dead, or of hurting yourself in some wayNot at all, Several days, More than half the days, Nearly every day
Questions 10 and 11 asked only if answer to question 9 is not “Not at all”
10Has there been a time in the past month when you have had serious thoughts about ending your life?Yes, No
11Have you ever, in your WHOLE LIFE, tried to kill yourself or made a suicide attempt?Yes, No
HOSPITAL ANXIETY AND DEPRESSION SCALE (HADS, ANXIETY SYMPTOMS)
Tick the box beside the reply that is closest to how you have been feeling in the past week. Don't take too long over you replies: your immediate is best
1I feel tense or wound up:Most of the time; A lot of the time; From time to time, occasionally; Not at all
2I get a sort of frightened feeling as if something bad is about to happen:Very definitely and quite badly; Yes, but not too badly; A little, but it doesn't worry me; Not at all
3Worrying thoughts go through my mind:A great deal of the time; A lot of the time; From time to time, but not too often; Only occasionally
4I can sit at ease and feel relaxed:Definitely, Usually, Not often, Not at all
5I get a sort of frightened feeling like 'butterflies' in the stomach:Not at all, Occasionally, Quite often, Very often
6I feel restless and have to be on the move:Very much indeed, Quite a lot, Not very much, Not at all
7I get sudden feelings of panic:Very often indeed, Quite often, Not very often, Not at all
CPSS QUESTIONNAIRE (PTSD SYMPTOMS)
Instructions to participant: below is a list of problems that kids sometimes have after experiencing an upsetting event. Read each one carefully
and circle the number (0–3) that best describes how often that problem has bothered you in the last 2 weeks
1Please select your most distressing event:2016 Fort McMurray wildfire, Death of someone close to you, Injury that you suffered, Physical assault against you, Sexual assualt, Other
2How long as it been since the event (in years)?<1 month, 2–5 months, 6–11 months, 1 year, 2 years, 3–5 years, 6–10 years, 11 or more years
Below is a list of problems that kids sometimes have after experiencing an upsetting event. Read each one carefully and circle the number (0–3)
that best describes how often that problem has bothered you in the last 2 weeks
3Having upsetting thoughts or images about the event that came into your head when you didn't want them toNot at all or only at one time, Once a week or less/once in a while, 2 to 4 times a week/half the time, 5 or more times a week/almost always
4Having bad dreams or nightmaresSame as above
5Acting or feeling as if the event was happening again (hearing something or seeing a picture about it and feeling as if I am there again)Same as above
6Feeling upset when you think about it or hear about the event (for e.g., feeling scared, angry, sad, guilty, etc)Same as above
7Having feelings in your body when you think about or hear about the event (for e.g., breaking out into a sweat, heart beating fast)Same as above
8Trying not to think about, talk about, or have feelings about the eventSame as above
9Trying to avoid activities, people, or places that remind you of the traumatic eventSame as above
10Not being able to remember an important part of the upsetting eventSame as above
11Having much less interest or doing things you used to doSame as above
12Not feeling close to people around youSame as above
13Not being able to have strong feelings (for e.g., being unable to cry or unable to feel happy)Same as above
14Feeling as if your future plans or hopes will not come true (for e.g., you will not have a job or getting married or having kids)Same as above
15Having trouble falling or staying asleepSame as above
16Feeling irritable or having fits of angerSame as above
17Having trouble concentrating (for e.g., losing track of a storey on the television, forgetting what you read, not paying attention in class)Same as above
18Being overly careful (for e.g., checking to see who is around you and what is around you)Same as above
19Being jumpy or easily startled (for e.g., when someone walks up behind you)Same as above
CRAFFT QUESTIONNAIRE (DRUGS/ALCOHOL/TABACCO)
During the past 12 months, did you:
1Drink any alcohol (more than a few sips)?Yes, No
2Smoke any marijuana or hashish?Yes, No
3Use anything else to get high?Yes, No
4Have you ever ridden in a CAR driven by someone (including yourself) who was “high” or had been using alcohol or drugs?Yes, No
Questions 5–9 asked only if “yes” to one or more of questions 1–3
5Do you ever use alcohol or drugs to RELAX, feel better about yourself, or fit in?Yes, No
6Do you ever use alcohol or drugs while you are by yourself, or ALONE?Yes, No
7Do you every FORGET things you did while using alcohol or drugs?Yes, No
8Do your FAMILY or FRIENDS ever tell you that you should cut down on your drinking or drug use?Yes, No
9Have you ever gotten into TROUBLE while you were using alcohol or drugs?Yes, No
Tobacco use questionnaire
1During the past month, did you smoke tobacco products?Yes, No
2During the past month, did you use smokeless tobacco products?Yes, No
Rosenberg self-esteem scale
1On the whole, I am satisfied with myselfStrongly agree, Agree, Disagree, Strongly disagree
2At times, I think I am no good at allStrongly agree, Agree, Disagree, Strongly disagree
3I feel that I have a number of good qualitiesStrongly agree, Agree, Disagree, Strongly disagree
4I am able to do things as well as most other peopleStrongly agree, Agree, Disagree, Strongly disagree
5I feel I do not have much to be proud ofStrongly agree, Agree, Disagree, Strongly disagree
6I certainly feel useless at timesStrongly agree, Agree, Disagree, Strongly disagree
7I feel that I'm a person of worth, at least on an equal plane with othersStrongly agree, Agree, Disagree, Strongly disagree
8I wish I could have more respect for myselfStrongly agree, Agree, Disagree, Strongly disagree
9All in all, I am inclined to feel that I am a failureStrongly agree, Agree, Disagree, Strongly disagree
10I take a positive attitude towards myselfStrongly agree, Agree, Disagree, Strongly disagree
KIDSCREEN QUESTIONNAIRE (QUALITY OF LIFE)
Thinking about the last week:
1Have you physically felt fit and well?Not at all, Slightly, Moderately, Very, Extremely
2Have you felt full of energy?Never, Seldom, Quite often, Very often, Always
3Have you felt sad?Never, Seldom, Quite often, Very often, Always
4Have you felt lonely?Never, Seldom, Quite often, Very often, Always
5Have you had enough time for yourself?Never, Seldom, Quite often, Very often, Always
6Have you been able to do the things that you want to do in your free time?Never, Seldom, Quite often, Very often, Always
7Have your parent(s) treated you fairly?Never, Seldom, Quite often, Very often, Always
8Have you had fun with your friends?Never, Seldom, Quite often, Very often, Always
9Have you got on well at school?Not at all, Slightly, Moderately, Very, Extremely
10Have you been able to pay attention?Never, Seldom, Quite often, Very often, Always
11In general, how would you say your health is?Excellent, Very good, Good, Fair, Poor
CHILD AND YOUTH RESILIENCE MEASURE (CYRM-12)
To what extent do the sentences below describe you? Select an answer for each statement
1I am able to solve my problems without harming myself or othersNot at all, A little, Some-what, Quite a bit, A lot
2I know where to go in the community to get helpNot at all, A little, Some-what, Quite a bit, A lot
3Getting an education is important to meNot at all, A little, Some-what, Quite a bit, A lot
4I try to finish what I startNot at all, A little, Some-what, Quite a bit, A lot
5I have people I look up toNot at all, A little, Some-what, Quite a bit, A lot
6My parents/caregivers know a lot about meNot at all, A little, Some-what, Quite a bit, A lot
7My family stands by me during difficult timesNot at all, A little, Some-what, Quite a bit, A lot
8My friends stand by me during difficult timesNot at all, A little, Some-what, Quite a bit, A lot
9I have opportunities to develop skills that will be useful later in lifeNot at all, A little, Some-what, Quite a bit, A lot
10I am treated fairly in my communityNot at all, A little, Some-what, Quite a bit, A lot
11I feel I belong at schoolNot at all, A little, Some-what, Quite a bit, A lot
12I enjoy my cultural and family traditionsNot at all, A little, Some-what, Quite a bit, A lot

Questionnaire details [table reproduced with permission from our previously-published paper Brown et al. (13)].

Survey Administration Procedure

The vast majority (>98%) of students who completed the survey did so during regular school hours. A few students citing special circumstances completed the survey from home on their own computers. Depending on their school, students either completed the survey using a desktop computer in a computer laboratory or used laptops brought to their classroom. The survey website was based on an HTML/CSS front end and a back end server written in the Clojure programming language (http://clojure.org). A survey description script was read to each class at the beginning of the survey session (reproduced in the Supplementary Material). The script explained the purpose of the survey and provided instructions for completing the survey. It also explained that the survey was anonymous (participants were not asked for their names nor date of birth) and that participation was voluntary. Students had an opportunity to ask questions before participating. The survey battery included 96 questions in total. Participation required <20 min for most students, but a small number of students took up to 50 min. Participants were able to skip questions, but the survey description script and the survey website did encourage them to answer all questions.

Repeated Measurements in 2017, 2018, and 2019

The first mental health survey of Fort McMurray Grade 7–12 students was completed in November 2017 and was then repeated in November 2018 and November 2019. Many students were present for the survey across all 3 years. Other students were present for only 1 or 2 years. For example, students entering Grade 7 in September 2018 or September 2019 were not surveyed in earlier years (i.e., November 2017 and November 2018, respectively). Similarly, students graduating from Grade 12 in June 2018 or June 2019 were not re-surveyed the following November. Some students also moved to, or away from, Fort McMurray between 2017 and 2019. In summary, some of the surveys collected in 2017, 2018, and 2019 were completed by unique students (i.e., no repetition), although the majority were completed by the same students (i.e., repeated measures over 2 or 3 years). Because survey responses were anonymous, we were unable to link surveys completed by specific students from 1 year to the next. We therefore employed a between-subject analysis (treating all survey samples as independent) rather than using a repeated-measures, within-subject analysis over successive years. It is important to note that adopting a between-subject approach is statistically conservative as there is an expectation of increased overall error variance compared to the likely advantage of within-subject analysis.

Cut-Off Scores and Probable Diagnoses

Probable diagnoses of four different psychiatric conditions were established by thresholding each participant's scores on specific scales. Threshold values for probable diagnoses were derived from the relevant literature for each scale, as described below. Specific probable diagnoses included PTSD (based on the CPSS scale), depression (from the PHQ-A), anxiety (from the HADS), and alcohol/substance use disorder (based on the CRAFFT).

The term “probable diagnosis” is used here, as opposed to “clinical diagnosis,” because the scores were based on self-report scales rather than psychiatric clinical interviews. The literature reports good agreement between psychiatric clinical diagnoses of PTSD, depression, anxiety, and alcohol/substance use disorder with probable diagnoses derived from widely-published threshold scores for the above four questionnaires (57, 58, 6265), and we have previously used this approach (66).

Probable PTSD was determined based on a CPSS score of 15 or more (65, 67). Probable depression was determined based on a PHQ-A score of 11 or more (63). Probable moderately severe depression was determined based on a PHQ-A score of 15 or more (62). Suicidal thinking was determined from responses to two questions from the PHQ-A: question 9 “Over the past 2 weeks, how often have you been bothered by any of the following problems: Thoughts that you would be better off dead, or of hurting yourself in some way?” and question 10 “Has there been a time in the past month when you have had serious thoughts about ending your life?” Participants were assessed as exhibiting suicidal thinking if they answered “Several days,” “More than half the days,” or “Nearly every day” to PHQ-A question 9 and “Yes” to question 10. Participants answering “Not at all” to question 9 skipped (were not shown) question 10, and they were assessed as not exhibiting suicidal thinking. In addition, participants answering “Several days,” “More than half the days,” or “Nearly every day” to PHQ-A question 9 and “No” to question 10 were assessed as not exhibiting suicidal thinking (as distinct from thinking about self-harm) (PHQ-A question 11 was not considered in the definition of suicidal thinking). Probable anxiety was determined based on a HADS score of 11 or more (64). Probable alcohol/substance use disorder was determined based on a CRAFFT score of 2 or more (57, 58). Tobacco use was determined as answering “yes” to either of the two questions on the Tobacco Use Questionnaire. Finally, an “Any of 4 probable diagnoses” criterion was defined as being positive for one or more of the four probable diagnoses: PTSD, depression, anxiety, or alcohol/substance use disorder.

Dependent Variables and Statistical Effects Tested

We examined the following 15 dependent measures: (1) CPSS PTSD score, (2) PHQ-A depression score, (3) HADS anxiety score, (4) CRAFFT alcohol/substance use score, (5) Rosenberg self-esteem score, (6) Kidscreen quality of life score, (7) CYRM-12 resilience score, (8) percent probable PTSD, (9) percent probable depression, (10) percent probable moderately severe depression, (11) percent suicidal thinking, (12) percent probable anxiety, (13) percent probable alcohol/substance use disorder, (14) percent tobacco use, and (15) percent any of four probable diagnoses (See “Cut-off Scores and Probable Diagnoses” section above for details of scoring cut-offs used to determine probable diagnoses). In analysing data for a given measure (e.g., mean CPSS) or probable diagnosis (e.g., probable PTSD), we included only those participants who provided answers for all questions in the relevant questionnaire.

For each of the 15 dependent measures, we tested five statistical effects: (1) linear effect of time (2017 vs. 2018 vs. 2019), (2) linear effect of age (11 to 19 years old), (3) effect of female vs. male gender identity, (4) effect of other vs. female/male gender identity, and (5) effect of preferred not to say vs. female/male gender identity.

Participant age was determined by Demographics Questionnaire question 4 “What is your age?,” with answer choices “10 years or less,” “11 years,” “12 years,” “13 years,” “14 years,” “15 years,” “16 years,” “17 years,” “18 years,” “19 years,” and “20 years or more.” Note that participants answering “10 years or less” or “20 years or more” were excluded to enable modelling of age as a linear variable (see “Data Exclusion” section in the Results).

Participant gender identity was determined based on Demographics Questionnaire question 3 “What gender do you identify with?,” with answer choices “female,” “male,” “other,” and “prefer not to say.” Statistical test 3 (female vs. male) compared participants answering “female” vs. those answering “male.” Test 4 (other vs. female/male) compared participants answering “other” vs. those answering either “female” or “male.” Test 5 (preferred not to say vs. female/male) compared participants answering “prefer not to say” vs. those answering either “female” or “male.” Each test of a gender effect included only those participants with the relevant gender identities (test 3: female and male; test 4: female, male, and other; test 5: female, male, and preferred not to say).

Details of Statistical Analysis

All statistical comparisons were done using permutation testing on the slope parameter from a fitted linear model, with a null hypothesis of zero slope. The linear model included a slope parameter for the effect variable (time, age, female vs. male, other vs. female/male, preferred not to say vs. female/male) as well as parameters for “nuisance variables” as described below. Mathematical details of the linear modelling procedure are included in the “Details of Linear Modelling” section below.

Permutation testing is a non-parametric method and was chosen for its robustness against non-normality. The number of iterations was 105 for all permutation tests. All tests were two-tailed (That is, to compute the p-value for a given test, the absolute value of the slope parameter fitted to the real data was compared against the absolute values of the 105 simulated slope parameters fitted to permuted data).

In total, our analysis of the five statistical effects for each of 15 dependent variables included 75 individual statistical tests. We addressed multiple comparisons using the Benjamini-Hochberg method for false discovery rate (FDR) correction. This method computed a threshold of p = 0.025 for FDR correction across all 75 tests.

Distributions of gender identities and ages were similar across time, and the distribution of gender identities was similar across different ages (see “Demographics” in the Results section). Nonetheless, to address the possibility that results for one effect might have been driven in part by some small difference in one or more of the other effect variables, we included “nuisance variables” in the linear models to which permutation testing was applied. To test effects of time (2017 vs. 2018 vs. 2019), we used a linear model with a term for time as well as nuisance variables including a covariate for age and four indicator variables for gender identity: female, male, other, and preferred not to say. For analyses on the effects of time, only the fitted slope parameter for time was used to generate p-values. Including the other nuisance variables allowed the model to separate out effects of age and gender from effects of time. Similarly, analyses of age used five nuisance variables, including a time covariate and four indicator variables for gender identities. Analyses of gender effects (female vs. male, other vs. female/male, preferred not to say vs. female/male) included time and age as nuisance variables.

As expected from a sample of students in Grade 7–12, there were substantially fewer participants who were 11 or 18–19 years old, compared to those who were 12–17 years old, at the time of data collection (see “Demographics” in the Results section). To address the possibility that results for effects of age might be driven by leverage effects from smaller sample numbers in the extremes (aged 11 or 18–19 years old), we ran a follow-up analysis which included only participants aged 12–17.

We performed all analyses using in-house computer code written in the Clojure programming language (http://clojure.org). The code for statistical testing and FDR correction is available at http://github.com/mbrown/mrgbstats.

Details of Linear Modelling

For a given analysis, we defined an effect variable xi and a dependent variable yi, as well as J covariate nuisance variables vj, i and K indicator nuisance variables wk, i, where i∈[1, N] with N being the number of surveys used in the analysis.

For analyses on effects of time or effects of age, the effect variable xi was mean-centred to make it “independent” (orthogonal) to the intercept (i.e., constant offset). For analyses on effects of gender (female vs. male, other vs. female/male, preferred not to say vs. female/male), the effect variable was categorical and therefore not mean-centred. Covariate nuisance variables vj, i were mean-centred, and indicator nuisance variables wk, i were not mean-centred.

For analyses on effects of time or effects of age, we did not include an intercept term in the model (i.e., a constant offset column containing all ones in the model matrix). Because analyses of time and age included four nuisance indicator variables for gender identities, an intercept column would have been a linear combination of those four indicator variables, rendering the model matrix degenerate. The effect variable for time or for age was mean-centred and therefore orthogonal to the constant offset, in any case. For analyses of effects of gender, which did not have the four nuisance variables for gender, we did include a constant offset term.

For each analysis, we created a model matrix X from the effect variable xi, J covariate nuisance variables vj, i, and K indicator nuisance variables wk, i, as well as a constant offset term for analyses of gender effects. Each survey i contributed one effect variable value, one dependent variable value, and one value for each of (J + K) nuisance variables. For example, for analyses on effects of time, xi was time; v1, i was age; w1, i, w2, i, w3, i, and w4, i were indicator variables for gender identity, including female, male, other, and preferred not to say; and there was no constant offset column. To take a second example, for analyses of other vs. female/male gender identity, xi was 0 for participants identifying as female or male and 1 for those with other gender identity; v1, i and v2, i were time and age, respectively; there were no indicator nuisance variables; and there was a constant offset column. We created a model matrix X with size N by (1 + J + K) or else N by (2 + J + K), for analyses without and with a constant offset column, respectively. The first column of X consisted of the effect variable values xi. If a constant offset column was included, it was the second column. The next J columns were comprised of the mean-centred covariate nuisance variables vj, i, and the last K columns were comprised of the indicator nuisance variables wk, i. We defined the vector where was the vector of dependent variable values yi. The slope parameter used for permutation testing was , the first element of , representing the fitted scaling parameter for the effect variable.

Results

The survey was administered to all Grade 7–12 students in Fort McMurray, Alberta, Canada who were attending either the Public or Catholic schools on the days the survey was conducted during the month of November in 2017, 2018, and 2019. Five Public schools and two Catholic schools participated in the survey. In total, 9,920 surveys were collected during the period of 2017 to 2019. Forty five percentage were collected from Public schools and 55% from Catholic schools. As all surveys were anonymous, there was no way to identify which students did or did not repeat the survey over successive years.

Data Exclusion

Data from 544 surveys were excluded based the following exclusion criteria:

  • Participant answered “10 years or less” to the age question (Demographics Questionnaire question 4).

  • Participant answered “20 years or more” to the age question (Demographics Questionnaire question 4).

  • Participant gave inconsistent answers for positive and negative questions included in the Rosenberg questionnaire (for details, see Appendix C: Exclusion Criteria in the Supplementary Material).

  • Participant gave inconsistent answers for positive questions from the Rosenberg questionnaire and positive questions from the Kidscreen questionnaire (see Appendix C: Exclusion Criteria the Supplementary Material, for details).

  • Participant gave inconsistent answers among the non-reversed and reversed questions from the HADS questionnaire (For two of the HADS questions, the answer order is reversed to test for consistency; see details in Supplementary Material).

Criteria 1 and 2 above excluded participants with ambiguous age. The remaining participants had non-ambiguous ages in the range 11–19 years, allowing us to model effects of age as a linear variable. Criteria 3 to 5 above excluded participants who gave inconsistent answers, possibly because they were not paying attention to the survey or did not understand the questions. After exclusions, the final dataset included 9,376 surveys.

Demographics

Demographics for the 9,376 surveys were as follows. Self-reported gender identity was 47.3% female, 48.6% male, 1.7% other, and 2.4% preferred not to say. Age ranged from 11 to 19, and the mean age of participants was 14.3 years ± 1.8 (standard deviation). For additional demographic details, see Table 2. Distributions of gender identities and ages were similar across the 3 years of data collection (Table 2). The distribution of gender identities was similar across different ages as well (Table 2).

Table 2

Total collectedAfter exclusionsTotal students enrolledRecruitment rate
Number surveys collected and recruitment rates
All 3 years9,9209,37613,84867.7%
20173,2523,0704,40769.7%
20183,4513,2654,59271.1%
20193,2173,0414,84962.7%
FemaleMaleOtherDid not say
Gender distributions
All 3 years47.3%48.6%1.7%2.4%
201748.3%47.8%1.6%2.4%
201846.6%49.9%1.6%2.0%
201947.0%48.0%1.9%3.0%
11 yrs12 yrs13 yrs14 yrs15 yrs16 yrs17 yrs18 yrs19 yrs
Age distributions
All 3 years2.7%18.1%17.7%17.2%16.4%14.5%11.0%1.9%0.6%
20172.4%17.9%17.2%17.7%15.5%14.4%12.2%2.0%0.7%
20182.6%17.7%17.2%16.3%17.9%15.4%10.2%2.2%0.6%
20193.0%18.7%18.8%17.7%15.7%13.7%10.6%1.5%0.4%
11 yrs12 yrs13 yrs14 yrs15 yrs16 yrs17 yrs18 yrs19 yrs
Age distributions by gender
All2.7%18.1%17.7%17.2%16.4%14.5%11.0%1.9%0.6%
Female2.9%17.9%17.2%17.4%16.7%14.4%11.2%1.7%0.6%
Male2.3%18.5%17.7%17.2%16.2%14.6%10.8%2.1%0.5%
Other4.5%7.1%23.7%16.7%18.6%14.7%11.5%3.2%0.0%
Did Not Say3.1%19.6%24.0%16.0%13.3%13.3%8.9%0.9%0.9%

Demographics.

Did Not Say” refers to participants who answered the question on gender identity with “Prefer not to say.”

Changes Over Time (2017, 2018, 2019)

Table 3 shows results for analysis of changes over time. Permutation testing of the linear effect of time, with age and gender partialled out, indicated that 13 of the 15 dependent measures changed over time during the period from 2017 to 2019, with statistical significance surviving FDR multiple comparison correction. Scores for PTSD (CPSS), depression (PHQ-A), anxiety (HADS), and alcohol/substance use (CRAFFT) increased from 2017 to 2019. Self-esteem (Rosenberg) and quality of life (Kidscreen) scores decreased. Resilience (CYRM-12) scores did not change significantly. Rates of probable diagnoses of PTSD, depression, moderately severe depression, suicidal thinking, anxiety, alcohol/substance use disorder, and the “Any of 4 probable diagnoses” category all increased. Rates of tobacco use did not change significantly.

Table 3

MeasureN 2017N 2018N 2019Score 2017Score 2018Score 2019SlopeP-valueFDR
CPSS score (PTSD)2,8693,0352,82912.83 ± 11.4713.03 ± 11.3313.75 ± 11.870.470.0025*
PHQ-A score (depression)2,9643,1492,9268.02 ± 6.478.23 ± 6.488.69 ± 6.700.360.00001*
HADS score (anxiety)2,9833,1702,9397.74 ± 4.727.72 ± 4.688.01 ± 4.790.140.019*
CRAFFT score (alcohol/drugs)2,9943,1902,9470.55 ± 1.260.63 ± 1.300.61 ± 1.320.040.014*
Rosenberg score (self-esteem)2,9673,1382,92118.22 ± 6.5718.01 ± 6.5117.58 ± 6.68−0.330.00011*
Kidscreen score (quality of life)2,9773,1712,92626.97 ± 8.1926.75 ± 8.3726.40 ± 8.45−0.330.0028*
CYRM-12 score (resilience)2,9293,0962,85746.48 ± 9.1746.97 ± 9.1246.69 ± 9.310.120.32
MeasureN2017N2018N2019Rate 2017Rate 2018Rate 2019SlopeP-valueFDR
Probable PTSD2,8693,0352,82937%39%41%2.180.00070*
Probable depression2,9643,1492,92631%31%35%2.200.00033*
Probable moderately severe depression2,9643,1492,92617%18%20%1.660.0011*
Suicidal thinking2,9903,1752,94816%16%18%1.150.018*
Probable anxiety2,9833,1702,93927%28%31%2.160.00020*
Probable alcohol/substance use disorder2,9943,1902,94715%17%16%1.210.012*
Tobacco use3,0043,2012,95613%15%12%−0.160.71
Any of 4 probable diagnoses3,0523,2403,02346%49%51%2.650.00004*

Changes over time (2017, 2018, 2019).

Table shows participant numbers, summary statistics, and results for analysis of effects of time. N X columns indicate number of participants from year X. Score X columns show mean ± standard deviation (not standard error) for year X. Rate X columns indicate rates of probable diagnosis for year X. Slope is the fitted slope parameter from a linear model with age and gender partialled out (see discussion of nuisance variables in the Materials and Methods). The slope parameter indicates the linear relationship between time and the dependent variable, independent of age and gender. P-values were derived from permutation testing on the slope parameters. Comparisons with an asterisk under the FDR column survived FDR multiple comparison correction. P-threshold for FDR correction was 0.025.

Effects of Age (11–19 Years Old)

Table 4 shows results for the analysis of age effects. Permutation testing of the linear effect of age over the range 11–19 years old, with time and gender partialled out, revealed that all 15 dependent measures changed with age, with statistical significance surviving FDR multiple comparison correction. Scores for PTSD (CPSS), depression (PHQ-A), anxiety (HADS), and alcohol/substance use (CRAFFT) increased with age. Self-esteem (Rosenberg), quality of life (Kidscreen), and resilience (CYRM-12) scores decreased with age. Rates of probable diagnoses of PTSD, depression, moderately severe depression, suicidal thinking, anxiety, alcohol/substance use disorder, tobacco use, and the “Any of 4 probable diagnoses” category all increased with age.

Table 4

MeasureN 11N 12N 13N 14N 15N 16N 17N 18N 19Score 11Score 12Score 13Score 14Score 15Score 16Score 17Score 18Score 19SlopeP-valueFDR
CPSS score (PTSD)2161,5401,5251,5211,4441,2949771684813.09 ± 11.1911.97 ± 10.8312.99 ± 11.4113.44 ± 11.6812.61 ± 11.5113.85 ± 11.7314.39 ± 12.2316.20 ± 12.2116.44 ± 11.830.420.00001*
PHQ-A score (depression)2251,6111,5921,5721,4861,3261,002175507.23 ± 5.786.75 ± 5.867.62 ± 6.348.38 ± 6.698.43 ± 6.519.33 ± 6.8510.03 ± 6.7710.47 ± 6.4711.28 ± 6.870.600.00001*
HADS score (anxiety)2361,6321,5971,5731,4981,3281,005173507.40 ± 4.357.13 ± 4.467.58 ± 4.657.72 ± 4.787.82 ± 4.738.33 ± 4.838.71 ± 4.928.91 ± 4.738.48 ± 4.270.270.00001*
CRAFFT score (alcohol/drugs)2361,6471,6171,5781,4971,3271,007172500.04 ± 0.320.12 ± 0.600.30 ± 0.920.49 ± 1.130.73 ± 1.411.01 ± 1.561.26 ± 1.681.11 ± 1.651.92 ± 2.330.220.00001*
Rosenberg score (self-esteem)2331,6181,5901,5621,4811,3239981715018.63 ± 6.3619.08 ± 6.3518.13 ± 6.7617.57 ± 6.7617.83 ± 6.5517.60 ± 6.6217.17 ± 6.3516.32 ± 5.8715.84 ± 5.96−0.330.00001*
Kidscreen score (quality of life)2371,6251,6001,5721,4911,3269981755028.93 ± 8.2029.22 ± 8.0927.92 ± 8.1326.72 ± 8.4526.17 ± 7.9424.92 ± 8.2924.14 ± 8.1123.42 ± 7.5622.42 ± 8.02−0.970.00001*
CYRM-12 score (resilience)2251,5851,5501,5371,4701,3089881714847.48 ± 9.3647.98 ± 8.8746.46 ± 9.4946.13 ± 9.4646.78 ± 9.0846.64 ± 8.8946.32 ± 9.0144.50 ± 9.7645.73 ± 10.65−0.250.00002*
MeasureN11N12N13N14N15N16N17N18N19Rate 11Rate 12Rate 13Rate 14Rate 15Rate 16Rate 17Rate 18Rate 19SlopeP-valueFDR
Probable PTSD2161,5401,5251,5211,4441,2949771684838%34%38%40%36%41%43%53%48%1.650.00001*
Probable depression2251,6111,5921,5721,4861,3261,0021755026%22%28%33%33%38%44%43%46%3.680.00001*
Probable moderately severe depression2251,6111,5921,5721,4861,3261,0021755012%11%16%20%17%23%27%26%34%2.580.00001*
Suicidal thinking2301,6381,6081,5781,4961,3281,0091765017%14%18%17%15%17%19%23%28%0.570.0088*
Probable anxiety2361,6321,5971,5731,4981,3281,0051735026%22%27%28%28%33%36%40%28%2.350.00001*
Probable alcohol/substance use disorder2361,6471,6171,5781,4971,3271,007172501%3%8%13%20%28%35%30%42%5.930.00001*
Tobacco use2351,6481,6201,5821,5061,3351,010175503%4%7%11%17%21%24%22%44%4.010.00001*
Any of 4 probable diagnoses2471,6831,6491,6041,5251,3511,0231805339%37%43%49%49%57%65%63%62%4.640.00001*

Effects of age (11–19 years old).

Table shows participant numbers, summary statistics, and results of analysis of effects of age 11 to 19 years old. N X columns indicate number of participants from age group X. Score X columns show mean ± standard deviation (not standard error) for age group X. Rate X columns indicate rates of probable diagnosis for age group X. Slope is the fitted slope parameter from a linear model with time and gender partialled out (see discussion of nuisance variables in the Materials and Methods). The slope parameter indicates the linear relationship between age and the dependent variable, independent of time and gender. P-values were derived from permutation testing on the slope parameters. Comparisons with an asterisk under the FDR column survived FDR multiple comparison correction. P-threshold for FDR correction was 0.025.

We did follow-up analyses of age restricted to participants aged 12–17 years old, with time and gender partialled out. Suicidal thinking did not change significantly over time with participants aged 12–17 years (p = 0.11), in contrast to analysis of age 11–19 years, which found a statistically significant increase in the rate of suicidal thinking with age (p = 0.0088). The other 14 dependent measures exhibited the same pattern of statistically significant changes in analyses with age 12–17 years (CYRM-12 p = 0.00038, CPSS p = 0.00002, all other tests p = 0.00001) as compared to analyses with age 11–19 years.

Effects of Gender Identity

Tables 57 present results of three analyses of gender identity. Gender identity was determined based on Demographics Questionnaire question 3 “What gender do you identify with?,” with answer choices “female,” “male,” “other,” and “prefer not to say.” Analyses of gender identity compared female vs. male (Table 5), other vs. female/male (Table 6), and preferred not to say vs. female/male (Table 7) (See “Dependent Variables and Statistical Effects Tested” in the Materials and Methods section for additional details). Analyses of gender effects were done using permutation testing, with time and age partialled. All 15 dependent measures showed significant differences for all three gender identity comparisons, with statistical significance surviving FDR multiple comparison correction for all tests. Scores for PTSD (CPSS), depression (PHQ-A), anxiety (HADS), and alcohol/substance use (CRAFFT) were higher in females vs. males, in those with other gender identity vs. females/males, and in those who preferred not to say vs. females/males. Self-esteem (Rosenberg), quality of life (Kidscreen), and resilience (CYRM-12) scores were lower in females, in participants with other gender identity, and in participants who preferred not to say. Rates of probable diagnoses of PTSD, depression, moderately severe depression, suicidal thinking, anxiety, alcohol/ substance use disorder, tobacco use, and the “Any of 4 probable diagnoses” category were higher in females, in those with other gender identity, and in those who preferred not to say.

Table 5

MeasureN maleN femaleScore maleScore femaleSlopeP-valueFDR
CPSS score (PTSD)4,2014,19010.61 ± 10.3015.39 ± 11.964.790.00001*
PHQ-A score (depression)4,3914,2976.76 ± 5.889.66 ± 6.752.900.00001*
HADS score (anxiety)4,4054,3256.52 ± 4.279.01 ± 4.772.500.00001*
CRAFFT score (alcohol/drugs)4,4384,3300.52 ± 1.200.64 ± 1.330.120.00002*
Rosenberg score (self-esteem)4,3754,28719.55 ± 6.1516.53 ± 6.53−3.020.00001*
Kidscreen score (quality of life)4,3974,31428.91 ± 7.7224.72 ± 8.29−4.200.00001*
CYRM-12 score (resilience)4,2924,24847.25 ± 9.0546.64 ± 9.05−0.620.0014*
MeasureNmaleNfemaleRate maleRate femaleSlopeP-valueFDR
Probable PTSD4,2014,19029%47%17.600.00001*
Probable depression4,3914,29723%40%17.560.00001*
Probable moderately severe depression4,3914,29712%24%12.540.00001*
Suicidal thinking4,4264,32612%21%8.480.00001*
Probable anxiety4,4054,32517%38%21.030.00001*
Probable alcohol/substance use disorder4,4384,33014%17%3.460.00001*
Tobacco use4,4514,34112%14%2.050.0040*
Any of 4 probable diagnoses4,5324,40239%58%19.020.00001*

Effects female vs. male gender identity.

Table shows participant numbers, summary statistics, and results of analysis of effects of female vs. male gender identity. N X columns indicate number of participants identifying as gender X (female or male). Score X columns show mean ± standard deviation (not standard error) for gender X. Rate X columns indicate rates of probable diagnosis for gender X. Slope is the fitted slope parameter from a linear model with time and age partialled out (see discussion of nuisance variables in the Materials and Methods). The slope parameter indicates the linear relationship between the female vs. male effect variable and the dependent variable, independent of time and age. P-values were derived from permutation testing on the slope parameters. Comparisons with an asterisk under the FDR column survived FDR multiple comparison correction. P-threshold for FDR correction was 0.025.

Table 6

MeasureN female/maleN otherScore female/maleScore otherSlopeP-valueFDR
CPSS score (PTSD)8,39114213.00 ± 11.4219.38 ± 14.026.300.00001*
PHQ-A score (depression)8,6881458.20 ± 6.4912.52 ± 7.674.180.00001*
HADS score (anxiety)8,7301507.75 ± 4.699.93 ± 5.112.130.00001*
CRAFFT score (alcohol/drugs)8,7681490.58 ± 1.271.13 ± 1.770.520.00001*
Rosenberg score (self-esteem)8,66214818.06 ± 6.5214.25 ± 8.25−3.740.00001*
Kidscreen score (quality of life)8,71114826.83 ± 8.2822.98 ± 9.74−3.640.00001*
CYRM-12 score (resilience)8,54013946.95 ± 9.0540.28 ± 11.33−6.640.00001*
MeasureNfemale/maleNotherRate female/maleRate otherSlopeP-valueFDR
Probable PTSD8,39114238%56%17.780.00003*
Probable depression8,68814532%61%28.190.00001*
Probable moderately severe depression8,68814518%42%23.710.00001*
Suicidal thinking8,75215016%38%21.540.00001*
Probable anxiety8,73015028%43%15.030.00006*
Probable alcohol/substance use disorder8,76814916%28%11.510.00024*
Tobacco use8,79215213%24%10.130.00027*
Any of 4 probable diagnoses8,93415648%69%20.030.00002*

Effects of other vs. female/male gender identity.

Table shows participant numbers, summary statistics, and results of analysis of effects of other vs. female/male gender identity. N X columns indicate number of participants identifying as gender X (other or female/male). Score X columns show mean ± standard deviation (not standard error) for gender X. Rate X columns indicate rates of probable diagnosis for gender X. Slope is the fitted slope parameter from a linear model with time and age partialled out (see discussion of nuisance variables in the Materials and Methods). The slope parameter indicates the linear relationship between the other vs. female/male effect variable and the dependent variable, independent of time and age. P-values were derived from permutation testing on the slope parameters. Comparisons with an asterisk under the FDR column survived FDR multiple comparison correction. P-threshold for FDR correction was 0.025.

Table 7

MeasureN female/maleN preferred not to sayScore female/maleScore preferred not to saySlopeP-valueFDR
CPSS score (PTSD)8,39120013.00 ± 11.4217.09 ± 13.664.110.00001*
PHQ-A score (depression)8,6882068.20 ± 6.4910.25 ± 7.212.140.00001*
HADS score (anxiety)8,7302127.75 ± 4.699.11 ± 5.501.400.00004*
CRAFFT score (alcohol/drugs)8,7682140.58 ± 1.270.93 ± 1.620.410.00001*
Rosenberg score (self-esteem)8,66221618.06 ± 6.5215.54 ± 6.87−2.560.00001*
Kidscreen score (quality of life)8,71121526.83 ± 8.2824.33 ± 8.97−2.700.00001*
CYRM-12 score (resilience)8,54020346.95 ± 9.0541.62 ± 10.59−5.400.00001*
MeasureNfemale/maleNpreferred not to sayRate female/male %Rate preferred not to say %SlopeP-valueFDR
Probable PTSD8,391200385112.880.00022*
Probable depression8,688206324312.190.00024*
Probable moderately severe depression8,688206182911.800.00003*
Suicidal thinking8,75221116269.810.00018*
Probable anxiety8,730212284315.350.00001*
Probable alcohol/substance use disorder8,768214162510.660.00006*
Tobacco use8,79221713197.440.0012*
Any of 4 probable diagnoses8,934225486214.500.00003*

Effects of “prefer not to say” vs. female/male gender identity.

Table shows participant numbers, summary statistics, and results of analysis of effects of preferred not to say vs. female/male gender identity. N X columns indicate number of participants identifying as gender X (preferred not to say or female/male). Score X columns show mean ± standard deviation (not standard error) for gender X. Rate X columns indicate rates of probable diagnosis for gender X. Slope is the fitted slope parameter from a linear model with time and age partialled out (see discussion of nuisance variables in the Materials and Methods). The slope parameter indicates the linear relationship between the preferred not to say vs. female/male effect variable and the dependent variable, independent of time and age. P-values were derived from permutation testing on the slope parameters. Comparisons with an asterisk under the FDR column survived FDR multiple comparison correction. P-threshold for FDR correction was 0.025.

Discussion

This study investigated the multi-year impacts of wildfires on youth mental health. We examined 9,376 mental health survey samples from students in Grades 7–12 in Fort McMurray, Alberta. Surveys were conducted in 2017, 2018, and 2019, in the aftermath of the 2016 Fort McMurray wildfire.

As we have previously reported, the Fort McMurray student population exhibited elevated rates of probable depression, suicidal thinking, and tobacco use; elevated symptoms of anxiety; and reduced scores for quality of life and self-esteem 18 months after the 2016 wildfire, as compared to a control population that had not recently experienced a natural disaster (12). At that time, we observed similar mental health patterns for youth who were not actually present in Fort McMurray during the 2016 wildfire, although youth with greater personal exposure to impacts of the fire (e.g., home destroyed) exhibited worse symptoms of PTSD, depression, anxiety, and alcohol/substance use and lower scores for self-esteem and quality of life (13). The current study provides evidence of ongoing long-term mental health impacts on youth 3.5 years (42 months) post-wildfire.

Longer Term Mental Health Impacts

The results from the present study indicate a slow but statistically significant trend of worsening mental health from 2017 to 2019, including increased symptom scores and increased rates of probable diagnoses of PTSD, depression, anxiety, drug use, and alcohol use. Quality of life and self-esteem scores also decreased from 2017 to 2019. Tobacco use and resilience scores did not change significantly. These findings are consistent with other studies reporting long-term negative impacts on mental health from natural disasters (4145), although these studies are not directly comparable because they focused on mental health in adulthood and/or used different outcome measures.

Our results do not support our hypothesis that mental health would improve with time following the 2016 wildfire. One possibility is that 3.5 years may not be sufficient time for this population to recover from the adverse mental health effects of the wildfire, though that possibility seems unlikely in light of reports that many individuals affected by disaster do recover within 1–2 years (4). Theories of recovery from trauma identify various factors important to the recovery process including a sense of safety and stability, self- and community efficacy, hope, support from family and friends, social support, and social connectedness (4, 51, 52). It is possible that one or more factors important to recovery may be at issue. For example, Fort McMurray has experienced an economic downturn related to the reduction in oil prices starting in 2014, and this may negatively affect the community's sense of stability and hope. Future studies would be needed to test this suggestion.

It is noteworthy that CYRM-12 resilience measures did not change significantly with time (p = 0.27), indicating that the above-discussed changes in mental health measures may not be attributable to a change in resilience.

There was a statistically significant increase in suicidal thinking with age in the analysis of participants aged 11–19 (p = 0.0088) but not in the analysis of participants aged 12–17 (p = 0.11). The difference in results seems to be driven by the larger proportion of students aged 18 and 19 exhibiting suicidal thinking (18 years: 23%, 19 years: 28%) compared to younger students (14–19%) (see Table 4).

Age-Related Differences in Mental Health Impact

We observed worse average scores on all 15 dependent measures in older vs. younger students. Specifically, older students exhibited higher mental health symptom scores, higher rates of probable diagnoses, and lower scores for self-esteem, quality of life, and resilience. These results are consistent with previous reports of increased mental health impairment among older youth compared to younger youth post-disaster (68), as well as higher rates of mental health symptoms in older adolescents more generally (69). As has been suggested by others (68), one possible interpretation of this finding, that would bear future exploration, is that greater awareness among older youth regarding challenges facing their families and the larger community (rebuilding, economic implications, etc.), as well as concerns they may have regarding their future, may negatively influence their well-being compared to younger youth. That older students also exhibited higher mental health symptom scores and lower resilience scores is consistent previous research showing an association between lesser resilience and worse mental health outcomes following disaster (28, 35, 70, 71) and with theoretical conceptions of resilience and its role in buffering the individual's mental health from harm due to adverse experiences (7274). Reduced resilience may have played a role in making older students more vulnerable to developing negative mental health symptoms.

Differences in Mental Health Related to Gender Identity

Our analyses revealed worse average scores on all 15 dependent measures in students identifying as female vs. male, in students with other gender identity vs. females/males, and students who preferred not to say their gender identity vs. females/males. These specific groups of students exhibited higher mental health symptom scores, higher rates of probable diagnoses, and lower scores for self-esteem, quality of life, and resilience.

We interpret participants answering “other” to the question of their gender identity as belonging to a gender minority, including transgender or gender non-conforming. For participants who preferred not to answer the demographics question on gender identity, some presumably identified as female or male but did not want to say so, while others presumably identified as a gender minority, including transgender or gender-non-conforming. Given that the group who preferred not to answer exhibited significantly worse mental health results than those identifying as female or male, we suspect that a majority of the group who preferred not to answer in fact belonged to a gender minority, including transgender or gender-non-conforming.

The finding of worse mental health scores in females compared to males is consistent with a previous report of higher rates of mental health symptoms post-natural disaster in females vs. males (75). We are not aware of any previous studies examining the impact of natural disaster on the mental health of gender minorities, including transgender and gender non-conforming individuals. It is a benefit of our population study approach, with its large sample size, that we are able to do so for the first time. More generally, our results are consistent with previous reports indicating higher rates of mental health symptoms in gender minorities, including transgender and gender-non-conforming individuals, both adults and youth (7679). Resilience scores were lower for females vs. males and for gender minorities vs. females/males, which is consistent with previous studies showing an association between lesser resilience and worse mental health outcomes following disaster (28, 35, 70, 71). This suggests that reduced resilience may have been a factor in specific groups' developing more negative mental health symptoms.

Implications

The destructive nature of disasters tends to attract funding from government and charities to address the immediate aftermath of a disaster. Our results provide an example of worsening mental health impacts in youth during the period 1.5–3.5 years following a wildfire disaster. This has occurred in the context of an ongoing, whole-of-community “build back better” approach to post-wildfire recovery in the area (80), as well as multiple challenges that have faced the community since the 2016 wildfire, including a downturn in the economy (Additional challenges include disastrous flooding following the 2020 spring ice breakup and the outbreak of the COVID-19 pandemic in 2020, though these occurred after the last survey data collection was completed in 2019). Our results underscore the need for multi-year funding, interventions, and policies to address not only the short-term physical damage but also the long-term negative mental health effects of natural disasters.

In addition to focusing on symptomatology, there is a need to investigate factors associated with post-disaster recovery processes, as well as efficacy of psychosocial strategies during later phases of recovery relative to early interventions, while also recognising complex and evolving, social-ecological post-disaster contexts. As an example, given that specific groups exhibited greater negative mental health effects, namely older youth, females, and gender minorities, interventions and policies may be more effective if they take into account developmental stage and gender identity with respect to mitigating these effects.

Limitations

Our analysis included a large dataset of 9,376 survey responses. Conducting full clinical interviews with this large number of participants is not feasible, and so we used clinical measures based on self-report questionnaires, which is a limitation of the study. As noted above, in addition to the 2016 wildfire, Fort McMurray has experienced an economic downturn since 2014. The resulting job losses and financial impacts on families also likely had an effect on youth mental health in addition to the 2016 wildfire. We are not aware of any study specifically on the effects of Alberta's economic downturn on youth mental health, though one report found a negative mental health impact in adults (81). Additionally, given that this study utilised anonymous data, we were not able to identify longitudinal trends in individuals, only in groups.

Conclusion

This study presents a cross-sectional statistical analysis of longitudinal mental health measurements in a population of Fort McMurray youth in Grades 7–12 during the period of 1.5–3.5 years following the 2016 Alberta wildfire. Findings indicate that there was a long-term trend of worsening mental health during that period. These observations support previous reports that youth and communities experience long-term mental health impacts following major natural disasters, such as wildfires. Our findings emphasise the need for multi-year funding and programs to support child and youth mental health in communities that have experienced such disasters.

Statements

Author’s note

Some portions of the manuscript, in the Materials and Methods section, are based on similar material from two previous papers (12, 13) presenting findings from separate analyses of the data from Fort McMurray collected in 2017 only.

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation. Requests for data access should be made by email to the corresponding author.

Ethics statement

This study involving human participants was reviewed and approved by the University of Alberta's Health Research Ethics Board (ethics protocol number Pro00072669 approved June 26, 2017). Data were collected as part of the Fort McMurray Public and Catholic Schools' standard curriculum to evaluate their support programs. Written informed consent from the participants' legal guardian/next of kin was not required to participate in this study, in accordance with the national legislation and the institutional requirements. Parents/guardians were notified of the study 2 weeks prior and were given the opportunity to opt their child(ren) out of the study. Participants themselves were given the opportunity to opt out of the study, as was explained to them at the beginning of each study survey session.

Author contributions

MB, VA, AG, PB-M, JD, CM-H, JO, MM, SN, DK, and PS: study design. MM, SN, and DK: data collection. MB, HP, IC, and PS: analysis. MB, HP, AG, IC, PB-M, JD, CM-H, BL, and PS: manuscript preparation. All authors contributed to the article and approved the submitted version.

Funding

Funding to support the Fort McMurray schools in developing their survey and to conduct statistical analyses and reporting was provided by a collaborative grant from the Canadian Institutes of Health Research, Canadian Red Cross, and Alberta Innovates Health Solutions (Grant No. 201600546).

Acknowledgments

We thank Fort McMurray Public and Catholic School Boards for allowing us to analyze an anonymous dataset derived from their survey. We also thank all individuals involved in the survey data collection including Fort McMurray teachers, counsellors, and other staff members.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpsyt.2021.676256/full#supplementary-material

    Abbreviations

  • CPSS

    Child PTSD Symptom Scale

  • CRAFFT

    CRAFFT Questionnaire (proper name of the questionnaire is CRAFFT)

  • CYRM-12

    Child and Youth Resilience Measure

  • EMPATHY

    Empowering a Multimodal Pathway Towards Healthy Youth project

  • FDR

    false discovery rate

  • HADS

    Hospital Anxiety and Depression Scale

  • Kidscreen-10

    Kidscreen Questionnaire

  • PHQ-A

    The Patient Health Questionnaire, Adolescent version

  • PTSD

    post-traumatic stress disorder

  • Rosenberg

    Rosenberg Self-Esteem Scale.

References

  • 1.

    CBC. Devastating Fort McMurray Wildfire Declared Out 15 Months Later. CBC (2017). Available online at: http://www.cbc.ca/news/canada/edmonton/fort-mcmurray-fire-beast-extinguished-out-1.4271604 (accessed February 16, 2021).

  • 2.

    CrydermanK. Fort McMurray Wildfires to Cost Insurers $3.6-Billion - The Globe and Mail. (2016). Available online at: https://www.theglobeandmail.com/report-on-business/fort-mcmurray-wildfire-damage-to-cost-36-billion-insurance-bureau/article30788517/ (accessed February 16, 2021).

  • 3.

    CohanCLColeSW. Life course transitions and natural disaster: marriage, birth, and divorce following hurricane hugo. J Fam Psychol. (2002) 16:1425. 10.1037/0893-3200.16.1.14

  • 4.

    BonannoGABrewinCRKaniastyKGrecaAM. Weighing the costs of disaster: consequences, risks, and resilience in individuals, families, and communities. Psychol Sci Public Interest. (2010) 11:149. 10.1177/1529100610387086

  • 5.

    SpeierAH. Psychosocial Issues for Children and Adolescents in Disasters. 2nd ed. Rockville, MD. (2000). Available online at: http://cretscmhd.psych.ucla.edu/nola/Video/MHR/Governmentreports/Psychosocial%20Issues%20for%20Children%20and%20Adolescents%20in%20Disasters.pdf (accessed February 16, 2021).

  • 6.

    PeekL. Children and disasters: understanding vulnerability, developing capacities, and promoting resilience - an introduction. Child Youth Environ. (2008) 18:129.

  • 7.

    Substance Abuse and Mental Health Services Administration (SAMHSA). Disaster Technical Assistance Center Supplemental Research Bulletin: Behavioral Health Conditions in Children and Youth Exposed to Natural Disasters. Rockville, MD (2018). Available online at: https://www.samhsa.gov/sites/default/files/srb-childrenyouth-8-22-18.pdf (accessed February 16, 2021).

  • 8.

    NorrisFHFriedmanMJWatsonPJByrneCMDiazEKaniastyK. 60,000 disaster victims speak: part I. An empirical review of the empirical literature, 1981-2001. Psychiatry. (2002) 65:20739. 10.1521/psyc.65.3.207.20173

  • 9.

    NorrisFHFriedmanMJWatsonPJ. 60,000 disaster victims speak: part II. Summary and implications of the disaster mental health research. Psychiatry. (2002) 65:24060. 10.1521/psyc.65.3.240.20169

  • 10.

    PeekLFothergillA. Post-Disaster decline: understanding children's vulnerability before, during, and after katrina. Dialogue Can Philos Assoc. (2014) 10:113. Available online at: https://www.samhsa.gov/sites/default/files/dialogue-vol10-issue3.pdf

  • 11.

    World Health Organization. Adolescent Mental Health. WHO (2019). Available online at: https://www.who.int/news-room/fact-sheets/detail/adolescent-mental-health (accessed February 16, 2021).

  • 12.

    BrownMRGAgyapongVGreenshawAJCribbenIBrett-MacLeanPDroletJet al. After the Fort McMurray wildfire there are significant increases in mental health symptoms in grade 7–12 students compared to controls. BMC Psychiatry. (2019) 19:18. 10.1186/s12888-019-2074-y

  • 13.

    BrownMRGAgyapongVGreenshawAJCribbenIBrett-MacLeanPDroletJet al. Significant PTSD and other mental health effects present 18 months after the Fort McMurray wildfire: findings from 3,070 grades 7–12 students. Front Psychiatry. (2019) 10:623. 10.3389/fpsyt.2019.00623

  • 14.

    RubensSLFelixEDHambrickEP. A meta-analysis of the impact of natural disasters on internalizing and externalizing problems in youth. J Trauma Stress. (2018) 31:33241. 10.1002/jts.22292

  • 15.

    JonesRTRibbeDPCunninghamPBWeddleJDLangleyAK. Psychological impact of fire disaster on children and their parents. Behav Modif. (2002) 26:16386. 10.1177/0145445502026002003

  • 16.

    PsarrosCTheleritisCEconomouMTzavaraCKioulosKTMantonakisLet al. Insomnia and PTSD one month after wildfires: evidence for an independent role of the “fear of imminent death.” Int J Psychiatry Clin Pract. (2017) 21:13741. 10.1080/13651501.2016.1276192

  • 17.

    BryantRAGibbsLGallagherHCPattisonPLusherDMacDougallCet al. Longitudinal study of changing psychological outcomes following the victorian black saturday bushfires. Aust New Zeal J Psychiatry. (2017) 52:54251. 10.1177/0004867417714337

  • 18.

    McDermottBMLeeEMJuddMGibbonP. Posttraumatic stress disorder and general psychopathology in children and adolescents following a wildfire disaster. Can J Psychiatry. (2005) 50:13743. 10.1177/070674370505000302

  • 19.

    PapadatouDGiannopoulouIBitsakouPBellaliTTaliasMATselepiK. Adolescents' reactions after a wildfire disaster in Greece. J Trauma Stress. (2012) 25:5763. 10.1002/jts.21656

  • 20.

    KuligJCTownshendIBoteyAPShepardB. ‘Hope is in our hands:' impacts of the slave lake wildfires in Alberta, Canada on children. In: SzenteJ, editor. Assisting Young Children Caught in Disasters. Educating the Young Child (Advances in Theory and Research, Implications for Practice), Vol. 13. Cham: Springer International Publishing (2018). p. 14356. 10.1007/978-3-319-62887-5_14

  • 21.

    KuligJCDabravolskajJ. The psychosocial impacts of wildland fires on children, adolescents and family functioning: a scoping review. Int J Wildl Fire. (2020) 29:93103. 10.1071/WF18063

  • 22.

    TownshendIAwosogaOKuligJPujadasBAShepardBMcFarlaneB. Impacts of wildfires on school children: a case study of Slave Lake. Int J Mass Emerg Disasters. (2015) 33:148187.

  • 23.

    AgyapongVIORitchieABrownMRGNobleSMankowsiMDengaEet al. Long-term mental health effects of a devastating wildfire are amplified by socio-demographic and clinical antecedents in elementary and high school staff. Front Psychiatry. (2020) 11:448. 10.3389/fpsyt.2020.00448

  • 24.

    RitchieASautnerBOmegeJDengaENwakaBAkinjiseIet al. Long-term mental health effects of a devastating wildfire are amplified by sociodemographic and clinical antecedents in college students. Disaster Med Public Health Prep. (2020) 111. 10.1017/dmp.2020.87

  • 25.

    KirschKRFeldtBAZaneDFHaywoodTJonesRWHorneyJA. Longitudinal community assessment for public health emergency response to wildfire, Bastrop County, Texas. Heal Secur. (2016) 14:93104. 10.1089/hs.2015.0060

  • 26.

    PapanikolaouVAdamisDMellonRCProdromitisG. Psychological distress following wildfires disaster in a rural part of Greece: a case-control population-based study. Int J Emerg Ment Health. (2011) 13:1126.

  • 27.

    FreedyJRSaladinMEKilpatrickDGResnickHSSaundersBE. Understanding acute psychological distress following natural disaster. J Trauma Stress. (1994) 7:25773. 10.1002/jts.2490070207

  • 28.

    AfifiWAFelixEDAfifiTD. The impact of uncertainty and communal coping on mental health following natural disasters. Anxiety, Stress Coping. (2012) 25:32947. 10.1080/10615806.2011.603048

  • 29.

    FreedyJResnickHKilpatrickD. Conceptual framework for evaluating disaster impact: implications for clinical intervention. In: AustinLS, editor. Responding to Disaster: A Guide for Mental Health ProfessionalsWashington, DC: American Psychiatry Press (1992). p. 323.

  • 30.

    MarshallGNSchellTLElliottMNRayburnNRJaycoxLH. Psychiatric disorders among adults seeking emergency disaster assistance after a wildland-urban interface fire. Psychiatr Serv. (2007) 58:50914. 10.1176/ps.2007.58.4.509

  • 31.

    PaveglioTBKooistraCHallTPickeringM. Understanding the effect of large wildfires on residents' well-being: what factors influence wildfire impact?For Sci. (2016) 62:5969. 10.5849/forsci.15-021

  • 32.

    Caamano-IsornaFFigueirasASastreIMontes-MartínezATaracidoMet al. Respiratory and mental health effects of wildfires: an ecological study in galician municipalities (north-west Spain). Environ Health. (2011) 10:48. 10.1186/1476-069X-10-48

  • 33.

    LangleyAK. Coping Efforts and Efficacy, Acculturation, and Post-Traumatic Symptomatology in Adolescents Following Wildfire: A Latent Variable Path Analysis. (Unpublished Ph.D., dissertation) (2000). Available online at: https://vtechworks.lib.vt.edu/handle/10919/26473

  • 34.

    WeemsCFGrahamRA. Resilience and trajectories of posttraumatic stress among youth exposed to disaster. J Child Adolesc Psychopharmacol. (2014) 24:28. 10.1089/cap.2013.0042

  • 35.

    GoldmannEGaleaS. Mental health consequences of disasters. Annu Rev Public Health. (2014) 35:16983. 10.1146/annurev-publhealth-032013-182435

  • 36.

    KarN. Psychological impact of disasters on children: review of assessment and interventions. World J Pediatr. (2009) 5:511. 10.1007/s12519-009-0001-x

  • 37.

    TangBLiuXLiuYXueCZhangL. A meta-analysis of risk factors for depression in adults and children after natural disasters. BMC Public Health. (2014) 14:623. 10.1186/1471-2458-14-623

  • 38.

    DaiWChenLLaiZLiYWangJLiuA. The incidence of post-traumatic stress disorder among survivors after earthquakes: a systematic review and meta-analysis. BMC Psychiatry. (2016) 16:188. 10.1186/s12888-016-0891-9

  • 39.

    KõlvesKKõlvesKEDe LeoD. Natural disasters and suicidal behaviours: a systematic literature review. J Affect Disord. (2013) 146:114. 10.1016/j.jad.2012.07.037

  • 40.

    MastenASOsofskyJD. Disasters and their impact on child development: introduction to the special section. Child Dev. (2010) 81:102939. 10.1111/j.1467-8624.2010.01452.x

  • 41.

    ForresiBSonciniFBottossoEDi PietroEScarpiniGScainiSet al. Post-traumatic stress disorder, emotional and behavioral difficulties in children and adolescents 2 years after the 2012 earthquake in Italy: an epidemiological cross-sectional study. Eur Child Adolesc Psychiatry. (2020) 29:22738. 10.1007/s00787-019-01370-0

  • 42.

    McFarlaneACVan HooffM. Impact of childhood exposure to a natural disaster on adult mental health: 20-Year longitudinal follow-up study. Br J Psychiatry. (2009) 195:1428. 10.1192/bjp.bp.108.054270

  • 43.

    FergussonDMHorwoodLJBodenJMMulderRT. Impact of a major disaster on the mental health of a well-studied cohort. JAMA Psychiatry. (2014) 71:102531. 10.1001/jamapsychiatry.2014.652

  • 44.

    GaoXLengYGuoYYangJCuiQGengBet al. Association between earthquake experience and depression 37 years after the tangshan earthquake: a cross-sectional study. BMJ Open. (2019) 9:e026110. 10.1136/bmjopen-2018-026110

  • 45.

    IsaranuwatchaiWCoytePCMcKenzieKNohS. The 2004 tsunami and mental health in Thailand: a longitudinal analysis of one-and two-year post-disaster data. Disasters. (2017) 41:15070. 10.1111/disa.12188

  • 46.

    BryantRAWatersEGibbsLGallagherHCPattisonPLusherDet al. Psychological outcomes following the victorian black saturday bushfires. Aust N Z J Psychiatry. (2014) 48:63443. 10.1177/0004867414534476

  • 47.

    SchwartzRMGillezeauCNLiuBLieberman-CribbinWTaioliE. Longitudinal impact of hurricane sandy exposure on mental health symptoms. Int J Environ Res Public Health. (2017) 14:957. 10.3390/ijerph14090957

  • 48.

    RakerEJLoweSRArcayaMCJohnsonSTRhodesJWatersMC. Twelve years later: the long-term mental health consequences of hurricane katrina. Soc Sci Med. (2019) 242:112610. 10.1016/j.socscimed.2019.112610

  • 49.

    Self-BrownSLaiBSThompsonJEMcGillTKelleyM Lou. Posttraumatic stress disorder symptom trajectories in hurricane katrina affected youth. J Affect Disord. (2013) 147:198204. 10.1016/j.jad.2012.11.002

  • 50.

    KronenbergMEHanselTCBrennanAMOsofskyHJOsofskyJDLawrasonB. Children of katrina: lessons learned about postdisaster symptoms and recovery patterns. Child Dev. (2010) 81:124159. 10.1111/j.1467-8624.2010.01465.x

  • 51.

    KleberRJ. Trauma and public mental health: a focused review. Front Psychiatry. (2019) 10:451. 10.3389/fpsyt.2019.00451

  • 52.

    DückersMLA. Five essential principles of post-disaster psychosocial care: looking back and forward with Stevan Hobfoll. Eur J Psychotraumatol. (2013) 4:21914. 10.3402/ejpt.v4i0.21914

  • 53.

    FoaEBJohnsonKMFeenyNCTreadwellKR. The child PTSD symptom scale: a preliminary examination of its psychometric properties. J Clin Child Psychol. (2001) 30:37684. 10.1207/S15374424JCCP3003_9

  • 54.

    JohnsonJGHarrisESSpitzerRLWilliamsJBW. The patient health questionnaire for adolescents: validation of an instrument for the assessment of mental disorders among adolescent primary care patients. J Adolesc Health. (2002) 30:196204. 10.1016/S1054-139X(01)00333-0

  • 55.

    SpitzerRLKroenkeKWilliamsJBWGroup the PHQPCS. Validation and utility of a self-report version of PRIME-MD: the PHQ primary care study. JAMA. (1999) 282:1737. 10.1001/jama.282.18.1737

  • 56.

    ZigmondASSnaithRP. The hospital anxiety and depression scale. Acta Psychiatr Scand. (1983) 67:36170. 10.1111/j.1600-0447.1983.tb09716.x

  • 57.

    KnightJRShrierLABravenderTDFarrellMVander BiltJShafferHJ. A new brief screen for adolescent substance abuse. Arch Pediatr Adolesc Med. (1999) 153:5916. 10.1001/archpedi.153.6.591

  • 58.

    KnightJRSherrittLHarrisSKGatesECChangG. Validity of brief alcohol screening tests among adolescents: a comparison of the AUDIT, POSIT, CAGE, and CRAFFT. Alcohol Clin Exp Res. (2003) 27:6773. 10.1111/j.1530-0277.2003.tb02723.x

  • 59.

    RosenbergM. Society and the Adolescent Self-Image. Princeton, NJ: Princeton University Press (1965).

  • 60.

    Ravens-SiebererUErhartMRajmilLHerdmanMAuquierPBruilJet al. Reliability, construct and criterion validity of the KIDSCREEN-10 score: a short measure for children and adolescents' well-being and health-related quality of life. Qual Life Res. (2010) 19:1487500. 10.1007/s11136-010-9706-5

  • 61.

    LiebenbergLUngarMLeBlancJC. The CYRM-12: a brief measure of resilience. Can J Public Health. (2013) 104:e1315. 10.1007/BF03405676

  • 62.

    KroenkeKSpitzerRLWilliamsJB. The PHQ-9: validity of a brief depression severity measure. J Gen Intern Med. (2001) 16:60613. 10.1046/j.1525-1497.2001.016009606.x

  • 63.

    RichardsonLPMcCauleyEGrossmanDCMcCartyCARichardsJRussoJEet al. Evaluation of the patient health questionnaire-9 item for detecting major depression among adolescents. Pediatrics. (2010) 126:111723. 10.1542/peds.2010-0852

  • 64.

    SnaithRP. The hospital anxiety and depression scale. Health Qual Life Outcomes. (2003) 1:29. 10.1186/1477-7525-1-29

  • 65.

    StewartRWEbesutaniCDrescherCFYoungJ. The child PTSD symptom scale: an investigation of its psychometric properties. J Interpers Violence. (2015) 32:223756. 10.1177/0886260515596536

  • 66.

    MrklasKShalabyRHrabokMGusnowskiAVuongWSuroodSet al. Prevalence of perceived stress, anxiety, depression, and obsessive-compulsive symptoms in health care workers and other workers in alberta during the COVID-19 pandemic: cross-sectional survey. JMIR Ment Health. (2020) 7:e22408. 10.2196/22408

  • 67.

    ISSTS. Child PTSD Symptom Scale. Available online at: https://istss.org/clinical-resources/assessing-trauma/child-ptsd-symptom-scale-for-dsm-5-(cpss-5) (accessed February 16, 2021).

  • 68.

    WeemsCFOsofskyJDOsofskyHJKingLSHanselTCRussellJD. Three-year longitudinal study of perceptions of competence and well-being among youth exposed to disasters. Appl Dev Sci. (2018) 22:2942. 10.1080/10888691.2016.1219229

  • 69.

    FordTGoodmanRMeltzerH. The british child and adolescent mental health survey 1999: the prevalence of DSM-IV disorders. J Am Acad Child Adolesc Psychiatry. (2003) 42:120311. 10.1097/00004583-200310000-00011

  • 70.

    BryantRAGallagherHCGibbsLPattisonPMacDougallCHarmsLet al. Mental health and social networks after disaster. Am J Psychiatry. (2017) 174:27785. 10.1176/appi.ajp.2016.15111403

  • 71.

    FelixEAfifiTKia-KeatingMBrownLAfifiWReyesG. Family functioning and posttraumatic growth among parents and youth following wildfire disasters. Am J Orthopsychiatry. (2015) 85:191200. 10.1037/ort0000054

  • 72.

    DavydovDMStewartRRitchieKChaudieuI. Resilience and mental health. Clin Psychol Rev. (2010) 30:47995. 10.1016/j.cpr.2010.03.003

  • 73.

    HerrmanHStewartDEDiaz-GranadosNBergerELJacksonBYuenT. What is resilience?Can J Psychiatry. (2011) 56:25865. 10.1177/070674371105600504

  • 74.

    ZolkoskiSMBullockLM. Resilience in children and youth: a review. Child Youth Serv Rev. (2012) 34:2295303. 10.1016/j.childyouth.2012.08.009

  • 75.

    LoweSRBonumweziJLValdespino-HaydenZGaleaS. Posttraumatic stress and depression in the aftermath of environmental disasters: a review of quantitative studies published in 2018. Curr Environ Heal reports. (2019) 6:34460. 10.1007/s40572-019-00245-5

  • 76.

    SuraceTFusar-PoliLVozzaLCavoneVArcidiaconoCMammanoRet al. Lifetime prevalence of suicidal ideation and suicidal behaviors in gender non-conforming youths: a meta-analysis. Eur Child Adolesc Psychiatry. (2020) 1:3. 10.1007/s00787-020-01508-5

  • 77.

    ValentineSEShipherdJC. A systematic review of social stress and mental health among transgender and gender non-conforming people in the United States. Clin Psychol Rev. (2018) 66:2438. 10.1016/j.cpr.2018.03.003

  • 78.

    ReisnerSLKatz-WiseSLGordonARCorlissHLAustinSB. Social epidemiology of depression and anxiety by gender identity. J Adolesc Health. (2016) 59:2038. 10.1016/j.jadohealth.2016.04.006

  • 79.

    StreedCGMcCarthyEPHaasJS. Association between gender minority status and self-reported physical and mental health in the United States. JAMA Intern Med. (2017) 177:12102. 10.1001/jamainternmed.2017.1460

  • 80.

    Regional Municipality of Wood Buffalo. RMWB 2016 Wildfire Recovery Plan. (2018) Available online at: https://fusesocial.ca/wp-content/uploads/2018/05/RMWB-recovery-plan2106.pdf (accessed February 19, 2021).

  • 81.

    CheungSMarriottB. Impact of an Economic Downturn on Addiction and Mental Health Service Utilization: A Review of the Literature. (2015). Available online at: https://www.albertahealthservices.ca/assets/info/res/mhr/if-res-mhr-kn-12-economic-downturn.pdf (accessed February 22, 2021).

Summary

Keywords

youth, population mental health, wildfire natural disaster, post-traumatic stress disorder (PTSD), depression, anxiety, alcohol and substance use, resilience

Citation

Brown MRG, Pazderka H, Agyapong VIO, Greenshaw AJ, Cribben I, Brett-MacLean P, Drolet J, McDonald-Harker CB, Omeje J, Lee B, Mankowsi M, Noble S, Kitching DT and Silverstone PH (2021) Mental Health Symptoms Unexpectedly Increased in Students Aged 11–19 Years During the 3.5 Years After the 2016 Fort McMurray Wildfire: Findings From 9,376 Survey Responses. Front. Psychiatry 12:676256. doi: 10.3389/fpsyt.2021.676256

Received

04 March 2021

Accepted

15 April 2021

Published

20 May 2021

Volume

12 - 2021

Edited by

Lawrence T. Lam, University of Technology Sydney, Australia

Reviewed by

Je-Yeon Yun, Seoul National University Hospital, South Korea; Yongxin Li, Henan University, China

Updates

Copyright

*Correspondence: Peter H. Silverstone

This article was submitted to Public Mental Health, a section of the journal Frontiers in Psychiatry

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics