The Prevalence and Factors Associated with Musculoskeletal Pain Among Pilgrims During the Hajj.

BACKGROUND: Musculoskeletal pain is a primary burden on individuals as well as social and health care systems. Annually, 2-3 million pilgrims perform the Hajj in Mecca, Saudi Arabia. The Hajj is highly physically demanding because pilgrims generally move by foot for long distances among a series of religious sites, an effort that may exceed their typical levels of physical activity. To understand the impact of musculoskeletal pain on the completion of the Hajj, it is first necessary to evaluate the extent of the problem. Accordingly, this study aimed to estimate the prevalence of musculoskeletal pain and associated factors among pilgrims during the Hajj.
METHODS: A cross-sectional survey was conducted during the period of the Hajj. The participants were adult pilgrims ≥ 18 years of age. Data regarding demographics, the prevalence of falls and the point prevalence of musculoskeletal pain by anatomical site were recorded. Participants were allowed to report more than one site of pain. Prevalence, crude and adjusted risk ratios were calculated.
RESULTS: A total of 1715 pilgrims were included in the analysis. The prevalence of falls was 13.76%. The prevalence of overall musculoskeletal pain (pain at any site) was 80.46%. Musculoskeletal pain was most commonly reported in the ankle/foot (38.34%), leg (29.89%), lower back (28.47%) and knee (21.84%). In general, musculoskeletal pain at multiple sites was more common in females and in older and obese individuals. However, there were variations in the importance of sex, age and body mass index as associated factors across different pain sites.
CONCLUSION: Musculoskeletal pain is common among pilgrims. Unlike most populations examined in other studies, ankle/foot pain was the most common in pilgrims. These data provide guidance for potential preventative programs and the allocation of resources to optimize pilgrims' experiences and ability to complete the Hajj.

Introduction

Musculoskeletal pain is common and imposes a major burden on individuals as well as social and healthcare systems,1,2 affecting both sexes, all ages, and all socio-cultural groups.1,3 Musculoskeletal pain is considered to be one of the most common causes of pain and physical disability, impacting hundreds of millions of people around the world.4,5 The prevalence of musculoskeletal pain is rising, and it has been described as an epidemic.1,6

“Hajj” is an Arabic word meaning “pilgrimage”,7 and it is the largest annual pilgrimage in the world, undertaken by Muslims at least once in their lifetime as a religious duty8–11 if they are financially and physically capable.12 Every year, around two to three million pilgrims converge simultaneously on the holy city of Mecca in Saudi Arabia (SA) for the rites of the Hajj, where they perform a series of rituals9,10,12,13 that were originally performed by the Prophet Mohammed.9 Based on the lunar Islamic year/calendar, the Hajj begins on the eighth day of Dhu al-Hijjah (the last month of the Islamic year) and ends on the twelfth day of the same month,10,12 although some pilgrims may spend further days to complete the Hajj.11 The date of the Hajj differs with respect to the Gregorian calendar, occurring 11 days earlier each year.8

This mass gathering, which leads to extreme congestion in small areas, may contribute to one of the most important public health problems in the world,9,14 resulting in high environmental and healthcare demands. The Hajj also involves high physical demands because pilgrims move, generally by foot, among a series of religious sites over 5–7 days while following a specific route, with average distances of 5–15 km/day,8,11,15 potentially reaching a total of 63 km during the whole period of the Hajj.16 This likely exceeds the typical physical activity level of most individuals and is further complicated by overcrowding, extreme heat and fatigue.7,9,11,17 Many pilgrims also maximize prayers (Salat) during Hajj,18 which involve repetitive motions between a series of postures, including standing, bowing, prostration, and sitting.19 In addition, individuals’ normal daily routines may be changed during the Hajj period as pilgrims move between places where the geography and climate are different, stay in tents and may neglect their self-health management while they are preoccupied with religious rituals.7,12

Many studies have been conducted to investigate health issues in pilgrims, with the main focus on infectious diseases9–12,14,20,21 and other diseases, including cardiovascular diseases, neurological disorders, trauma, gastrointestinal problems, diabetes, heat exhaustion and dermatological diseases.7,10,14,15,22–25 However, the prevalence of musculoskeletal pain among pilgrims during the Hajj and the potential associated factors have not been investigated. Thus, it is difficult to determine the impact of musculoskeletal pain on individuals who undertake the Hajj. Indeed, it is unclear whether the musculoskeletal pain among pilgrims is more prevalent in some physical sites than in others and whether the pain in these sites is associated with specific factors. To understand the potential impact of musculoskeletal pain on the completion of the Hajj, the development of preparatory advice and the preparation of support services, it is first necessary to evaluate the extent of the problem. Therefore, this study aimed to estimate the point prevalence of musculoskeletal pain by anatomical site among pilgrims during the Hajj and to study potential associations with individual characteristics.

Materials and Methods

Study Design

A cross-sectional survey utilizing convenience sampling was conducted among pilgrims during the period of the Hajj. Ethical approval was obtained for the study from the Physiotherapy Research Committee of the Faculty of Applied Medical Sciences at Umm Al-Qura University in Mecca, SA. This study was registered at in August 2018 (Research Registry Identifying Number: 4352), and it was conducted in accordance with the Declaration of Helsinki.

Participants

The participants were adult pilgrims aged 18 years or older of all nationalities, who had performed the Hajj in August 2018, which corresponded with Dhu al-Hijjah 1439 of the Hijri/Islamic calendar. Any participants who had performed Umrah were excluded because Umrah involves some similar religious tasks and can be performed within a few hours, whereas the Hajj requires over 5–7 days to complete. Furthermore, all individuals who were working to serve and provide services to pilgrims (e.g., transportation, safety, social care, drifters care or healthcare) but did not perform the Hajj themselves were excluded. Informed consent was obtained, and all included participants agreed to participate.

Data Collection

The data were collected after completion of the second day of the Hajj from 21st to 31st August 2018 (10 to 20 Dhu al-Hijjah 1439 of the Hijri/Islamic calendar). The data were collected by qualified healthcare professionals, who approached individuals at different sites in Mecca using an online application or a paper form. Pilgrims were either questioned verbally by healthcare professionals or answered the questionnaire in person while healthcare professionals were nearby to explain any questions if they needed help. All healthcare professionals were trained to perform data collection in a standardized manner. All data were collected from all main sites of Hajj rituals, including Muzdalifah, Mina and the Holy Mosque.

The data from the online survey were collected using a validated website (), and the data collected from the paper form were uploaded to the online application by the authors to obtain one single master file for both methods (online application and paper form). All information collected was anonymous, and only the authors had access to the data.

Survey

The survey data were sequentially collected (one by one) to avoid response bias in the sequence listed below. The survey was provided in two versions: an Arabic version and an English version (). The structure of the survey used in this study was as follows:

Demographics: Participants were asked to answer questions about their nationality/country, sex, age and body mass index (BMI).

Prevalence of falls: Participants were asked if they had lost their balance and fallen during the Hajj to estimate the prevalence of falls.

Point prevalence of musculoskeletal pain: Participants were asked if they had experienced musculoskeletal pain or discomfort (by anatomical site) during the Hajj. All participants were allowed to report more than one site of pain.

Sample Size Calculation

According to a statistical report of the Saudi General Authority for Statistics in 2017, the total number of pilgrims who visited the holy city of Mecca in SA to perform the Hajj in 2017 was 2,352,122. This estimated number includes Saudi (n=209,415) and non-Saudi (n=2,142,707‬) pilgrims.26 Based on this information, it was assumed that around 2–3 million pilgrims would come to SA to perform the Hajj in 2018. Thus, the sample size was calculated by setting the statistical power at a 99% confidence interval, with a population size of 3,000,000 and a margin of error of 5%, indicating that the required sample size for this study was 666 participants. After the completion of the Hajj 2018, the Saudi General Authority for Statistics revealed that the total number of pilgrims who performed the Hajj in 2018 was 2,371,675.13 This estimation was within our assumed range of the total number of pilgrims to perform the Hajj in 2018 (2–3 million).

Statistical Analysis

Incomplete responses were discarded when there were no data available about the falls and musculoskeletal pain together. The findings were considered statistically significant when P<0.05. The data were analyzed using StataIC version 16 (College Station, TX: StataCorp LLC).

Descriptive statistics (frequencies and percentages) were used to analyze participants’ responses. The investigation of sex- and age-based differences has been recommended, which would provide further information of epidemiological data and contribute to a better understanding of the prevalence of musculoskeletal pain.27–29 Therefore, the prevalence and their 95% confidence intervals (CIs) were calculated, and Pearson’s chi-squared two-tailed tests were used to identify sex- and age-based differences in the following measures:

Falls: The prevalence of falls was calculated for each age and sex group.

Musculoskeletal pain: The prevalence of overall musculoskeletal pain (pain at any site), the prevalence of total lower limb pain (pain at hip/pelvis, thigh, knee, leg and ankle/foot) and the prevalence of musculoskeletal pain for each anatomical site (shoulder, upper arm, elbow, forearm, wrist/hand, head, cervical, thoracic, lumbar, hip/pelvis, thigh, knee, leg and ankle/foot) were calculated for each age and sex group. To adequately estimate the prevalence of musculoskeletal pain, only data from pilgrims who did not experience falls were calculated. This is because falls are acute events that may result in musculoskeletal pain, and thus they may be considered a causal pathway, whereas the majority of musculoskeletal pain in the pilgrims was probably due to repetitive motion or overuse.

The risk ratio (RR) for musculoskeletal pain was calculated to measure the association between an exposure and an outcome and to identify whether a specific exposure is a risk factor for a specific outcome.30 The result of the RR can be interpreted as follows:31 exposure does not affect the risk of the outcome (RR=1), exposure is associated with a higher risk of the outcome (RR>1) and exposure is associated with a lower risk of the outcome (RR<1). Binomial generalized linear models with a log link function were used to obtain RRs.30 The analyses were limited to most common/important sites of musculoskeletal pain. Two models were performed:

Model 1: This model was a univariate regression analysis to estimate crude (unadjusted) RRs and their 95% CIs.

Model 2: This model was a multiple regression analysis to estimate RRs and their 95% CIs with an adjustment for sex, age and BMI.

Results

Characteristics of the Respondents

A total of 2110 responses were received in the study, of which 32 (1.52%) were from pilgrims who refused to participate, and 1715 (81.28%) were from pilgrims who agreed to participate and provided complete information for all the survey sections. Data from 59 nationalities were recorded, the majority from Arabic nationalities. Characteristics of the included respondents for falls and musculoskeletal pain analyses are shown in Table 1.

Table 1

Characteristics of the Respondents

Variables		N (%)
		Male	Female	Total
Respondents Included for Falls Analysis
Sex		920 (53.64)	795 (46.36)	1715 (100)
Age (years)	18–29	307 (52.03)	283 (47.97)	590 (34.40)
	30–39	250 (58.00)	181 (42.00)	431 (25.13)
	40–49	142 (49.65)	144 (50.35)	286 (16.68)
	50–59	133 (53.20)	117 (46.80)	250 (14.58)
	≥60	88 (55.70)	70 (44.30)	158 (9.21)
BMI (kg/cm2)	Underweight (<18.5)	45 (62.50)	27 (37.5‬0)	72 (4.20)
	Normal weight (18.5–24.9)	379 (50.60)	370 (49.40)	749 (43.67)
	Overweight (25–29.9)	210 (51.98)	194 (48.02)	404 (23.56)
	Obese (≥30)	286 (58.37)	204 (41.63)	490 (28.57)
Respondents Included for Musculoskeletal Pain Analysisa
Sex		811 (54.83)	668 (45.17)	1479 (100)
Age (years)	18–29	283 (53.60)	245 (46.40)	528 (35.70)
	30–39	228 (59.22)	157 (40.78)	385 (26.03)
	40–49	123 (50.20)	122 (49.80)	245 (16.57)
	50–59	103 (53.37)	90 (46.63)	193 (13.05)
	≥60	74 (57.81)	54 (42.19)	128 (8.65)
BMI (kg/cm2)	Underweight (<18.5)	39 (61.90)	24 (38.10)	63 (4.26)
	Normal weight (18.5–24.9)	329 (52.14)	302 (47.86)	631 (42.66)
	Overweight (25–29.9)	188 (53.11)	166 (46.89)	354 (23.94)
	Obese (≥30)	255 (59.16)	176 (40.84)	431 (29.14)

Note: aRespondents who experienced falls (n=263) during the Hajj were excluded from the musculoskeletal pain analysis.

Abbreviations: N, number; %, percentage; BMI, body mass index; kg, kilogram; cm, centimeter.

Prevalence of Falls

The prevalence of falls was 13.76%, and the prevalence was higher in females (15.97%) than in males (11.85%). The prevalence of falls was higher in older pilgrims (50–59 and ≥60 age groups) than in younger pilgrims, with the largest difference between the males and the females observed in pilgrims who were older than 60 years of age (Table 2).

Table 2

Prevalence of Falls by Sex and Age

Variables		Male	Female	P-value	Male and Female
Sex		11.85 (9.92–14.10)	15.97 (13.59–18.68)	0.013	13.76 (12.21–15.47)
Age Group	18–29	7.82 (5.31–11.37)	13.43 (9.94–17.89)	0.026	10.51 (8.28–13.24)
	30–39	8.80 (5.88–12.96)	13.26 (9.07–18.97)	0.139	10.67 (8.10–13.94)
	40–49	13.38 (8.74–19.95)	15.28 (10.31–22.05)	0.647	14.34 (10.75–18.87)
	50–59	22.56 (16.28–30.37)	23.08 (16.37–31.49)	0.922	22.80 (18.03–28.39)
	≥60	15.91 (9.72–24.95)	22.86 (14.59–33.95)	0.269	18.99 (13.64–25.81)
	P-value	0.000	0.054		0.000

Notes: The results of prevalence of falls with 95% CIs are presented as percentages (%). Based on Pearson’s chi-squared tests, the P-values of statistically significant differences (P<0.05) between sex and age groups are shown in bold.

Musculoskeletal Pain

A total of 1479 responses were included in the analysis of the prevalence of musculoskeletal pain, while 236 responses received from participants who experienced falls were excluded from this analysis.

Prevalence of Musculoskeletal Pain

The prevalence of overall musculoskeletal pain was 80.46%, with a prevalence of musculoskeletal pain for the total lower limb, total head/spine and total upper limb of 65.38%, 45.84% and 20.89%, respectively. The most prevalent musculoskeletal pain by anatomical site among pilgrims was ankle/foot pain (38.34%), followed by leg pain (29.89%), low back pain (28.47%) and knee pain (21.84%), as shown in Figure 1. There were sex- and age-based differences in multiple sites of pain (Table 3). The largest difference between males and females was observed for head and leg pain in younger pilgrims and for ankle/foot and knee pain in older pilgrims. The largest difference between younger and older pilgrims was observed for knee pain in both males and females.

Figure 1

The prevalence of musculoskeletal pain by site. This figure shows the point prevalence of musculoskeletal pain among pilgrims. The results are presented as frequencies and percentages for each site.

Table 3

Point Prevalence of Musculoskeletal Pain by Anatomical Site, Sex and Age

Pain Site	Sex	Age Group					P-value
		18–29	30–39	40–49	50–59	≥60
Shoulder	Male	12.37 (9.03–16.71)	13.16 (9.37–18.16)	7.32 (3.90–13.32)	16.50 (10.57–24.85)	14.86 (8.51–24.69)	0.290
	Female	17.14 (12.94–22.36)	17.83 (12.64–24.57)	20.49 (14.28–28.50)	16.67 (10.37–25.69)	11.11 (5.19–22.19)	0.669
Upper Arm	Male	1.41 (0.55–3.58)	3.95 (2.09–7.33)	4.88 (2.25–10.23)	1.94 (0.53–6.81)	2.70 (0.74–9.33)	0.261
	Female	3.27 (1.66–6.31)	3.18 (1.37–7.24)	0.00 (0.00–3.05)	0.00 (0.00–4.09)	1.85 (0.33–9.77)	0.131
Elbow	Male	1.41 (0.55–3.58)	2.19 (0.94–5.03)	1.63 (0.45–5.74)	0.97 (0.17–5.30)	4.05 (1.39–11.25)	0.575
	Female	0.41 (7e-04–2.28)	3.18 (1.37–7.24)	2.46 (0.84–6.98)	4.44 (1.74–10.88)	1.85 (0.33–9.77)	0.146
Forearm	Male	1.77 (0.76–4.07)	2.19 (0.94–5.03)	1.63 (0.45–5.74)	4.85 (2.09–10.86)	4.05 (1.39–11.25)	0.376
	Female	3.27 (1.66–6.31)	1.27 (0.35–4.53)	0.82 (0.14–4.50)	5.56 (2.40–12.35)	1.85 (0.33–9.77)	0.172
Wrist/Hand	Male	3.53 (1.93–6.38)	1.75 (0.68–4.42)	1.63 (0.45–5.74)	5.83 (2.70–12.13)	8.11 (3.77–16.58)	0.045
	Female	3.27 (1.66–6.31)	3.18 (1.37–7.24)	4.10 (1.76–9.24)	10.00 (5.35–17.92)	1.85 (0.33–9.77)	0.056
Head	Male	12.72 (9.33–17.11)	10.09 (6.82–14.68)	12.20 (7.53–19.15)	15.53 (9.79–23.75)	16.22 (9.53–26.24)	0.555
	Female	14.69 (10.81–19.67)	23.57 (17.61–30.79)	19.67 (13.59–27.60)	15.56 (9.50–24.43)	5.56 (1.91–15.11)	0.022
Neck	Male	6.36 (4.06–9.83)	12.28 (8.63–17.18)	11.38 (6.90–18.20)	17.48 (11.35–25.94)	17.57 (10.56–27.77)	0.007
	Female	17.14 (12.94–22.36)	13.38 (8.92–19.58)	13.11 (8.24–20.25)	10.00 (5.35–17.92)	5.56 (1.91–15.11)	0.158
Thoracic	Male	13.07 (9.64–17.50)	10.96 (7.54–15.69)	8.13 (4.48–14.32)	12.62 (7.53–20.40)	13.51 (7.51–23.12)	0.650
	Female	13.88 (10.10–18.77)	12.74 (8.40–18.86)	6.56 (3.36–12.41)	17.78 (11.25–26.94)	11.11 (5.19–22.19)	0.151
Lumbar	Male	19.43 (15.25–24.44)	28.07 (22.64–34.23)	28.46 (21.23–36.99)	32.04 (23.81–41.56)	36.49 (26.44–47.87)	0.010
	Female	31.43 (25.94–37.49)	32.48 (25.65–40.15)	32.79 (25.09–41.53)	21.11 (13.95–30.63)	37.04 (25.42–50.37)	0.249
Hip/Pelvis	Male	1.41 (0.55–3.58)	2.19 (0.94–5.03)	1.63 (0.45–5.74)	2.91 (1.00–8.22)	13.51 (7.51–23.12)	0.000
	Female	6.12 (3.75–9.85)	2.55 (1.00–6.37)	6.56 (3.36–12.41)	14.44 (8.64–23.16)	9.26 (4.02–19.91)	0.009
Thigh	Male	11.66 (8.42–15.92)	13.16 (9.37–18.16)	8.94 (5.07–15.31)	10.68 (6.07–18.12)	17.57 (10.56–27.77)	0.445
	Female	19.59 (15.11–25.01)	13.38 (8.92–19.58)	12.30 (7.59–19.30)	11.11 (6.15–19.26)	3.70 (1.02–12.54)	0.019
Knee	Male	12.37 (9.03–16.71)	17.98 (13.54–23.48)	19.51 (13.48–27.39)	30.10 (22.09–39.54)	39.19 (28.86–50.58)	0.000
	Female	17.96 (13.66–23.25)	14.01 (9.44–20.31)	28.69 (21.41–37.27)	46.67 (36.71–56.90)	37.04 (25.42–50.37)	0.000
Leg	Male	24.03 (19.42–29.33)	25.44 (20.22–31.47)	26.02 (19.07–34.41)	24.27 (17.02–33.38)	43.24 (32.57–54.59)	0.018
	Female	36.33 (30.56–42.52)	28.03 (21.59–35.51)	30.33 (22.87–38.98)	38.89 (29.47–49.22)	40.74 (28.68–54.03)	0.209
Ankle/Foot	Male	43.46 (37.81–49.29)	37.72 (31.68–44.17)	35.77 (27.85–44.56)	24.27 (17.02–33.38)	35.14 (25.24–46.50)	0.015
	Female	42.86 (36.82–49.12)	40.76 (33.39–48.58)	34.43 (26.59–43.22)	41.11 (31.51–51.44)	27.78 (17.62–40.89)	0.216
Lower Extremitya	Male	57.60 (51.78–63.22)	60.09 (53.61–66.23)	60.98 (52.15–69.14)	59.22 (49.57–68.22)	77.03 (66.25–85.13)	0.049
	Female	71.84 (65.90–77.10)	66.24 (58.54–73.17)	68.85 (60.17–76.39)	78.89 (69.37–86.05)	72.22 (59.11–82.38)	0.304
Allb	Male	69.96 (64.39–75.01)	77.19 (71.32–82.16)	80.49 (72.61–86.52)	78.64 (69.77–85.45)	83.78 (73.76–90.47)	0.040
	Female	85.71 (80.78–89.55)	82.80 (76.13–87.90)	86.89 (79.75–91.76)	90.00 (82.08–94.65)	87.04 (75.58–93.58)	0.616

Notes: aMusculoskeletal pain at the hip/pelvis, thigh, knee, leg and ankle/foot sites. bMusculoskeletal pain at any site. The results of point prevalence of musculoskeletal pain with 95% CIs are presented as percentages (%). Based on Pearson’s chi-squared tests, the values of statistically significant differences (P<0.05) between sexes (prevalence values) and age groups (P-values) are shown in bold.

Factors Associated with Musculoskeletal Pain

The ARRs results revealed (Table 4) that overall musculoskeletal pain was more common in females, while total lower limb musculoskeletal pain was more common in females and in older (≥60 years) and obese individuals (≥30 kg/cm2). Ankle/foot pain was less common in middle-aged (40–59 years) and older individuals (≥60 years). Musculoskeletal pain at other common sites (lower back, knee and leg) was more common in females and in older and obese individuals (Table 5). There was a significant association between ankle/foot pain and pain in the lower back, knee and leg (Table 6).

Table 4

Risk Ratios for Overall (Pain at Any Site) and Total Lower Limb (Pain in the Hip/Pelvis, Thigh, Knee, Leg and Ankle/Foot) Musculoskeletal Pain in the Hajj Population by Sex, Age and BMI

Variables		Model 1 (CRR)a	P-value	Model 2 (ARR)b	P-value
Overall Musculoskeletal Pain
Sex	Male	1	-	1	-
	Female	1.13 (1.08–1.19)	0.000	1.13 (1.07–1.18)	0.000
Age (years)	18–29	1	-	1	-
	30–39	1.03 (0.96–1.10)	0.421	1.02 (0.95–1.09)	0.632
	40–49	1.08 (1.01–1.16)	0.031	1.06 (0.99–1.14)	0.121
	50–59	1.09 (1.01–1.17)	0.036	1.07 (0.99–1.15)	0.092
	≥60	1.10 (1.01–1.20)	0.027	1.07 (0.98–1.16)	0.140
BMI (kg/cm2)	<18.5	1.00 (0.88–1.15)	0.958	1.03 (0.91–1.17)	0.652
	18.5–24.9	1	-	1	-
	25–29.9	1.03 (0.96–1.10)	0.358	1.02 (0.96–1.09)	0.469
	≥30	1.04 (0.98–1.10)	0.251	1.03 (0.97–1.09)	0.338
Total Lower Limb Musculoskeletal Pain
Sexc	Male	1	-	1	-
	Female	1.17 (1.08–1.26)	0.000	1.16 (1.08–1.25)	0.000
Age (years)c	18–29	1	-	1	-
	30–39	0.97 (0.88–1.08)	0.620	0.98 (0.88–1.08)	0.618
	40–49	1.01 (0.90–1.13)	0.851	0.99 (0.89–1.12)	0.983
	50–59	1.07 (0.95–1.20)	0.282	1.07 (0.96–1.20)	0.228
	≥60	1.17 (1.04–1.32)	0.010	1.14 (1.02–1.28)	0.024
BMI (kg/cm2)d	<18.5	0.99 (0.82–1.21)	0.933	1.03 (0.85–1.26)	0.731
	18.5–24.9	1	-	1	-
	25–29.9	1.01 (0.92–1.11)	0.834	1.03 (0.93–1.13)	0.610
	≥30	1.07 (0.98–1.16)	0.154	1.09 (1.01–1.19)	0.039

Notes: aUnadjusted risk ratio. bRisk ratio adjusted for sex, age and BMI. cRisk ratio adjusted for only sex and age as the convergence was not achieved when a third factor (BMI) was included in the model 2. dRisk ratio adjusted for only sex and BMI as the convergence was not achieved when a third factor (age) was included in the model 2. The results are presented as risk ratios with 95% CIs. Based on log-binomial models, the P-values of statistically significant associations (P<0.05) are shown in bold.

Abbreviations: CRR, crude risk ratio; ARR, adjusted risk ratio; CI, confidence interval; BMI, body mass index; kg, kilogram; cm, centimeter.

Table 5

Risk Ratios for the Most Common Sites of Musculoskeletal Pain in the Hajj Population by Sex, Age and BMI

Variables		Model 1 (CRR)a	P-value	Model 2 (ARR)b	P-value	Model 1 (CRR)a	P-value	Model 2 (ARR)b	P-value
		Ankle/Foot Pain				Leg Pain
Sex	Male	1	-	1	-	1	-	1	-
	Female	1.05 (0.92–1.20)	0.457	1.05 (0.92–1.20)	0.444	1.28 (1.10–1.50)	0.002	1.28 (1.09–1.49)	0.002
Age (years)	18–29	1	-	1	-	1	-	1	-
	30–39	0.90 (0.77–1.06)	0.204	0.90 (0.77–1.06)	0.202	0.89 (0.72–1.10)	0.286	0.88 (0.71–1.10)	0.260
	40–49	0.81 (0.67–0.99)	0.039	0.81 (0.66–0.99)	0.045	0.95 (0.75–1.20)	0.656	0.92 (0.72–1.17)	0.494
	50–59	0.74 (0.59–0.93)	0.011	0.74 (0.59–0.93)	0.011	1.05 (0.82–1.34)	0.725	1.03 (0.80–1.33)	0.821
	≥60	0.74 (0.57–0.97)	0.031	0.74 (0.56–0.97)	0.032	1.42 (1.11–1.81)	0.005	1.39 (1.08–1.78)	0.010
BMI (kg/cm2)	<18.5	0.99 (0.72–1.37)	0.969	0.97 (0.70–1.33)	0.842	0.86 (0.55–1.33)	0.492	0.85 (0.55–1.32)	0.463
	18.5–24.9	1	-	1	-	1	-	1	-
	25–29.9	0.85 (0.71–1.01)	0.065	0.91 (0.76–1.08)	0.278	1.03 (0.84–1.25)	0.774	1.01 (0.83–1.24)	0.914
	≥30	0.99 (0.85–1.15)	0.872	1.07 (0.91–1.25)	0.408	1.03 (0.85–1.24)	0.791	1.01 (0.84–1.23)	0.889
		Low Back Pain				Knee Pain
Sex	Male	1	-	1	-	1	-	1	-
	Female	1.17 (0.99–1.38)	0.051	1.18 (1.01–1.39)	0.045	1.24 (1.02–1.50)	0.031	1.23 (1.02–1.48)	0.029
Age (years)	18–29	1	-	1	-	1	-	1	-
	30–39	1.19 (0.97–1.48)	0.101	1.13 (0.91–1.40)	0.257	1.09 (0.81–1.48)	0.564	1.06 (0.78–1.44)	0.724
	40–49	1.22 (0.96–1.56)	0.097	1.13 (0.89–1.45)	0.323	1.61 (1.19–2.18)	0.002	1.50 (1.10–2.05)	0.010
	50–59	1.07 (0.82–1.41)	0.594	0.98 (0.73–1.29)	0.863	2.53 (1.93–3.32)	0.000	2.33 (1.75–3.10)	0.000
	≥60	1.47 (1.12–1.93)	0.005	1.36 (1.03–1.80)	0.028	2.56 (1.90–3.45)	0.000	2.37 (1.74–3.22)	0.000
BMI (kg/cm2)	<18.5	0.60 (0.33–1.07)	0.086	0.62 (0.35–1.12)	0.114	0.99 (0.57–1.74)	0.979	1.06 (0.61–1.84)	0.836
	18.5–24.9	1	-	1	-	1	-	1	-
	25–29.9	1.11 (0.90–1.37)	0.324	1.08 (0.87–1.33)	0.494	1.28 (0.99–1.66)	0.055	1.10 (0.85–1.42)	0.458
	≥30	1.23 (1.02–1.48)	0.033	1.20 (0.99–1.46)	0.062	1.60 (1.27–2.00)	0.000	1.30 (1.03–1.64)	0.025

Notes: aUnadjusted risk ratio. bRisk ratio adjusted for sex, age and BMI. The results are presented as risk ratios with 95% CIs. Based on log-binomial models, the P-values of statistically significant associations (P<0.05) are shown in bold.

Abbreviations: BMI, body mass index; CRR, crude risk ratio; ARR, adjusted risk ratio; CI, confidence interval; kg, kilogram; cm, centimeter.

Table 6

Association Between Ankle/Foot Pain and Pain in the Leg, Knee and Lower Back

Variables		Model 1 (CRR)a	P-value	Model 2 (ARR)b	P-value
Leg	No pain	1	-	1	-
	Pain	1.44 (1.27–1.64)	0.000	1.47 (1.29–1.67)	0.000
Knee	No pain	1	-	1	-
	Pain	1.35 (1.17–1.54)	0.000	1.43 (1.25–1.65)	0.000
Lower back	No pain	1	-	1	-
	Pain	1.21 (1.06–1.39)	0.004	1.24 (1.08–1.42)	0.002

Notes: aUnadjusted risk ratio. bRisk ratio adjusted for sex, age and BMI. Based on log-binomial models, the P-values of statistically significant associations (P<0.05) are shown in bold.

Abbreviations: CRR, crude risk ratio; ARR, adjusted risk ratio; CI, confidence interval; BMI, body mass index; kg, kilogram; cm, centimeter.

Discussion

This study found that most pilgrims were at high risk of experiencing musculoskeletal pain during the Hajj, with the most prevalent forms being ankle/foot pain, leg pain, low back pain and knee pain. The study findings also confirmed significant associations between the prevalence of musculoskeletal pain and factors related to sex, age and BMI.

Prevalence of Musculoskeletal Pain

Most pilgrims (80.46%) had at least one musculoskeletal complaint during the Hajj. This high prevalence rate is consistent with the findings of previous studies where the prevalence of musculoskeletal pain was 70% or higher for different specific populations including occupational drivers, agricultural farmers, marathon runners, nursing assistants, dentists and physiotherapists.32–38 However, several studies reported lower rates of musculoskeletal pain (<25%), particularly among the general population in different countries, such as Sweden, Brazil and Japan,39–41 indicating that the prevalence of musculoskeletal pain among pilgrims is higher than in the general population.

Although a few studies are available in the literature that estimated or discussed the prevalence of musculoskeletal pain during religious mass gatherings such as the Hajj or Arbaeenia, these studies have not comprehensively investigated the prevalence of musculoskeletal pain. For example, one study was conducted in 2016 in Iraq during Arbaeenia in Karbala City and found that 28.2% of individuals had joint pain related to walking long distances.42 Another study was conducted during the Hajj in 2010 and reported that the prevalence of limb pain among pilgrims was 7%.43 The findings of both studies were inconsistent with the findings of the present study as we found a higher prevalence of musculoskeletal pain. This difference in prevalence rates could be explained by the fact that the two studies collected data only from healthcare facilities/systems, whereas the majority of data in our study was from outside of healthcare facilities/systems since many pilgrims may not seek treatment or health services for mild musculoskeletal pain during the Hajj.

In the present study, the most common site of musculoskeletal pain was the ankle/foot (38.34%), with a prevalence higher than that identified in most other studies. While some studies reported higher rates of ankle/foot pain among some specific populations, such as nurses (43.8%)44 and housekeeping staff (58.3–60.4%),45 most studies reported lower rates of ankle/foot pain, with a prevalence lower than 20% in the general population or in specific populations such as nursing assistants, occupational drivers, construction workers, agricultural farmers and marathon runners.32,33,35,39,46–50 There are possible explanations why our study showed a higher prevalence of ankle/foot pain among pilgrims. One is that the Hajj is highly physically demanding. Pilgrims move long distances, generally by foot, which may exceed their typical physical activity levels.7,9,11,17 Several studies have documented an association between jobs or activities/tasks with high physical demands and musculoskeletal pain in different populations.34,39,42 For example, individuals who walk for long distances42 or primarily work while standing and for longer periods of time51,52 as well as those who are exposed to high physical workloads46 have an increased risk of musculoskeletal pain, particularly in the lower limbs.46,52 Another possible explanation could be related to the use of poorly fitting footwear/shoes among pilgrims, which may not be suitable to wear during the Hajj. Multiple studies have shown that poorly and incorrectly fitting shoes are significantly associated with ankle and foot pain,53,54 implying the importance of using properly fitting footwear. Therefore, pilgrims are at a higher risk of developing musculoskeletal pain, especially pain in the ankle/foot.

Factors Associated with Musculoskeletal Pain

The present study showed that the prevalence of overall musculoskeletal pain was greater in females and older individuals, in line with previous studies showing that sex and age are a risk factor for musculoskeletal pain. These studies reported a higher prevalence of musculoskeletal pain in females than in males39–41, 55–57 and in older individuals than in younger individuals.39,40,57 BMI was not associated with overall musculoskeletal pain in our study. In contrast, numerous studies reported that musculoskeletal pain was more common in individuals with a BMI exceeding that of normal-weight individuals.40,51 However, we found that a BMI of 30 kg/cm2 or greater was associated with a higher prevalence of total lower limb musculoskeletal pain.

In terms of the most prevalent site of musculoskeletal pain among the pilgrims in our study, sex and BMI were not associated with ankle/foot pain. In contrast, several studies found that females39,47,49–51,58 and obese individuals44,49,58 were at higher risk of ankle/foot pain. In our study, age was associated with ankle/foot pain in pilgrims, as older pilgrims reported less ankle/foot pain than younger pilgrims. This finding is inconsistent with several previous studies reporting that older pilgrims were at higher risk of ankle/foot pain than younger individuals.33,47,49,50,59 This difference could be explained by the fact that the majority of our participants were young or middle-aged adults (≤59 years, 91.35%). Additionally, older pilgrims may seek help to move between rituals sites during the Hajj, using wheelchairs or transportation services. However, we found that older pilgrims were at increased risk of musculoskeletal pain in other common sites, such as the leg, knee and lower back. In our study, we also found that pilgrims who reported pain in the leg, knee and lower back were more likely to report ankle/foot pain. This finding is similar to those of previous studies that found an association between ankle/foot pain and pain in other anatomical sites, including the leg,58 knee49,50,58 and back.49,50,58

Study Strengths and Limitations

To the best of our knowledge, this is the first study to comprehensively investigate the prevalence of musculoskeletal pain among pilgrims. One strength of this study is that a sample size calculation was performed, and the required sample size was achieved. The data were collected by qualified healthcare professionals, who were trained to perform the data collection in a standardized manner. Another strength is the investigation of the association between musculoskeletal pain and potential associated factors (sex, age and BMI). Nonetheless, there were some limitations to the current study. There was a potential sampling bias due to the use of convenience sampling, so the sample may not be representative of the entire pilgrim population that visited Mecca to perform the Hajj. However, an attempt was made to reduce the sampling bias by distributing the survey at different major sites of Hajj rituals. This study also used a cross-sectional design, which is limited to demonstrating a causal association.

Conclusion

The findings of this study indicate that musculoskeletal pain is common among pilgrims, and the reported prevalence of musculoskeletal pain is greater among pilgrims than that reported for the general population. This signals a significant public health issue that must be addressed by the Saudi Ministry of Health. Unlike most populations, the ankle/foot was the most common site of musculoskeletal pain among pilgrims. The study also demonstrates that sex, age and BMI were associated with the prevalence of musculoskeletal pain, and the importance of these factors varied across different pain sites.