ABSTRACT

Introduction

Neuromusculoskeletal injuries (NMSKI) are very common in the military, which contribute to short- and long-term disability.

Materials and Methods

Population-level NMSKI, limited duty (LIMDU), and long-term disability episode counts in the U.S. Navy (USN) and U.S. Marine Corps (USMC) from December 2016 to August 2021 were extracted from the Musculoskeletal Naval Epidemiological Surveillance Tool. The incidence of NMSKI, LIMDU, and long-term disability was calculated. A hurdle negative binomial regression evaluated the association of body region, sex, age, rank, age by rank, and service branch on NMSKI, LIMDU, and long-term disability incidence.

Results

From December 2016 to August 2021, there were 2,004,196 NMSKI episodes (USN: 3,285/1,000 Sailors; USMC: 4,418/1,000 Marines), 16,791 LIMDU episodes (USN: 32/1,000 Sailors; USMC: 29/1,000 Marines), and 2,783 long-term disability episodes (USN: 5/1,000 Sailors; USMC: 5/1,000 Marines). There was a large-magnitude protective effect on NMSKI during the pandemic (relative risk, USN: 0.70; USMC: 0.75). Low back and ankle-foot were the most common, primarily affecting female personnel, aged 25-44 years, senior enlisted, in the USMC. Shoulder, arm, pelvis-hip, and knee conditions had the greatest rates of disability, with female sex, enlisted ranks, aged 18-24 years, and service in the USMC having the most salient risk factors.

Conclusion

Body region, sex, age, rank, and branch were the salient factors for NMSKI. The significant protective effect during the pandemic was likely a function of reduced physical exposure and limited access to nonurgent care. Geographically accessible specialized care, aligned with communities with the greatest risk, is needed for timely NMSKI prevention, assessment, and treatment.

INTRODUCTION

Neuromusculoskeletal injuries (NMSKI) are very common in military personnel, which degrade the capability of the armed forces to meet mission objectives that contribute to national security.1 In best cases, these conditions preclude physical activity and limit the ability of a service member to function in their occupational specialties for only a limited period. More frequently, body system impairments and symptoms (such as pain) preclude function for protracted periods, which substantially affects the quality of life, increases health care costs, and degrades mission readiness.2 The association between military activity and injury is likely bidirectional since fluctuations in geopolitics and operational requirements influence the type and volume of exposure, environmental factors, and intrinsic factors associated with NMSKI in military service members.1

There are multiple factors associated with NMSKI, including sex, occupation, age, rank, and the period of study. These factors, in addition to the suspected latent factors such as care-seeking behaviors, access to care, and the mode and duration of exposure, plausibly drive these outcomes in the military.1 This likely includes periods of reduced operational demand,3–6 something that occurred during the 2020 COVID-19 pandemic as the U.S. Military shifted to a health protection posture required to mitigate the spread of this respiratory illness. A decline in NMSKI over time was observed in recent studies of the military, a finding that was attributed to the change in geopolitics, decline in combat operations, return to peacetime operations, and decreased exposure to hazards incurred during training and operations.3–6 This supposition can now be assessed in a more recent study epoch, where training and non-essential operations were halted to mitigate the spread of COVID-19.7

While military personnel frequently experience rapid recovery and return to occupational duties following NMSKI, others have protracted periods of activity limitation that preclude participation in the organizational mission.2 This is especially problematic considering the associated costs and implications for national security.1,2 An assessment of the factors that may be associated with NMSKI outcomes, short-term limited duty (LIMDU), and long-term disability (medical attrition) would allow for greater precision for targeted prevention, optimization of health care delivery models, and identification of factors that may contribute to protracted care and long-term disablement. Therefore, the purpose of this study was to assess the effects of sex, age, rank, service branch, and body region on the incidence of NMSKI, LIMDU, and long-term disability in the Department of the Navy (DoN), specifically members serving in the U.S. Navy (USN) and U.S. Marine Corps (USMC). A secondary aim was to evaluate the impact of the COVID-19 pandemic on the incidence of NMSKI.

METHODS

A population-based epidemiological retrospective cohort study of all military members in the DoN was performed to assess body region, sex, age, rank, and service branch on the incidence of NMSKI, LIMDU, and long-term disability following these injuries from December 2016 to August 2021. The Strengthening the Reporting of Observational Studies in Epidemiology Statement was used to guide reporting.

Neuromusculoskeletal injury episodes, LIMDU, and long-term disability (operationally defined by Disability Evaluation System [DES] referral) counts were extracted from the Musculoskeletal Naval Epidemiological Surveillance Tool (MSK NEST; Navy Bureau of Medicine and Surgery, Falls Church, VA, USA). This surveillance tool leverages existing validated databases, such as the Military Health System Data Repository, LIMDU Sailor and Marine Readiness Tracker System, DES, and Defense Manpower Data Center Reporting System to report injury burden. Neuromusculoskeletal injury episodes within the MSK NEST were derived from health care encounters in the Military Health System Data Repository and defined using the Army Public Health Center’s Taxonomy for Musculoskeletal Injuries based on an ICD-10 classification.8 Individuals with repeat visits for the same diagnosis in a single care episode were only counted once. The MSK NEST provides aggregated count data for NMSKI episodes by body region and de-identified patient characteristics that included sex; age range (18-24, 25-34, 35-44, and ≥45 years); rank (junior enlisted [E1-E5], senior enlisted [E6-E9], junior officers [O1-O3], and senior officers [O4-O10]); and service branch (USN and USMC) for all active and reserve personnel in the DoN. The database does not include any personally identifiable health information.

Since military end strength fluctuates annually because of attrition and recruitment of replacements,9 the population at risk was a dynamic cohort. The cumulative incidence of NMSKI was calculated and normalized to the subpopulation at risk (with consideration to sex, service branch, and rank) during the 4.75-year study epoch. Rates of LIMDU and long-term disability were calculated per 1,000 episode-years for each subpopulation. In the evaluation of the effects of COVID-19, the mean NMSKI incidence was calculated for the pre-pandemic epoch (December 2016 to March 2020) and during the pandemic (April 2020 to August 2021). Relative risk (RR) and 95% CI were calculated to assess the effects of the pandemic epoch.

A multivariable negative binomial regression was performed to evaluate the association of body region, sex, age, rank, age by rank, and service branch on the incidence of NMSKI episodes, LIMDU, and long-term disability. Because of overdispersion of the data, a negative binomial model was employed over a Poisson model.10 Since there were excess zeros noted in the data, with many demographic categories not having any reported outcomes of interest (especially among the smaller subpopulations), the assessment of convergence between the unadjusted and adjusted (hurdle) negative binomial regression models was performed. In the hurdle negative binomial regression model, a logistic link is employed to remove demographic categories that contributed no episodes of injury, LIMDU, or DES from the final count model. By doing so, excess zeros that contribute to skewed point estimates, standard errors, and overdispersion are parsed in the truncated regression model.10 Results of the hurdle negative binomial regression model are reported using calculations of the predictors regressed on count data (count model assessing the rate of the outcome), as well as a linked logistic regression (zero model that assessed the probability within demographic categories of not having the outcome). Both unadjusted and adjusted models are reported in the Supplementary material. The regression analyses were performed using the “MASS” (version 7.3-54) and “PSCL” (version 1.5) packages on R (version 3.5.1, The R Foundation for Statistical Computing, Vienna, Austria). Ankle-foot injuries served as the reference in the assessment of body region. Female military personnel were selected as the reference group in the assessment of sex. The age range of 18-24 years was selected as the reference for age. Junior enlisted personnel served the reference group for rank because of the greater disease and non-battle injuries in this group compared to commissioned officers.11 Finally, the USMC was the reference group for service branch. The level of significance was P < .05 for all analyses. Statistical significance was evaluated using the convergence of RR CIs that did not cross 1.00 and the P-value reported in the regression analyses.

RESULTS

The incidence of NMSKI, NMSKI-related LIMDU, and long-term disability stratified by body region, sex, age, rank, and service branch in the USN, USMC, and DoN is reported in Supplementary Tables S1-S3. From December 2016 to August 2021, there were 2,004,196 episodes of NMSKI (USN: 3,285 episodes/1,000 Sailors; USMC: 4,418 episodes/1,000 Marines), 16,791 cases of LIMDU (USN: 32/1,000 Sailors; USMC: 29/1,000 Marines), and 2,783 cases of long-term disability (USN: 5/1,000 Sailors; USMC: 5/1,000 Marines) in the DoN. Figure 1 details the incidence of the NMSKI in the USN, USMC, and the total DoN during pre-pandemic and pandemic epochs. The pre-pandemic rate of NMSKI was 761.00 per 1,000 person-years (95% CI: 733.41-788.08) in the USN; 1,004.73 per 1,000 person-years (95% CI: 973.03-1,035.92) in the USMC; and 847.93 per 1,000 person-years (95% CI: 818.81-876.55) in the total DoN. During the pandemic, the rate decreased to 534.65 per 1,000 person-years (95% CI: 511.53-557.27) in the USN; 751.23 per 1,000 person-years (95% CI: 723.82-778.18) in the USMC; and 608.96 per 1,000 person-years (95% CI: 584.28-633.13) in the DoN, an indication of a large-magnitude protective effect (RR, USN: 0.70; USMC: 0.75; DoN: 0.72). A large inflection in the NMSKI rate was also noted in August 2021 during the final month of the study.

Effects of the COVID-19 pandemic on episodic neuromuscular injury incidence in the U.S. Navy and U.S. Marine Corps.
FIGURE 1.

Effects of the COVID-19 pandemic on episodic neuromuscular injury incidence in the U.S. Navy and U.S. Marine Corps.

NMSKI Factors

Although the incidence of ankle-foot injuries was second to only low back conditions across the DoN (Table I, with the full hurdle model detailed in Supplementary Table S4), there were no significant differences between the rates of these conditions when assessed in the individual branches. Although male Sailors had lower risk of NMSKI compared to their female counterparts, there was no significant difference in rates observed between male and female Marines when controlling for all other factors. When compared to the 18-24-year-old group, there were higher incidence of injury observed in Sailors aged 25-44 years and lower rates in Marines aged ≥45 years. Rank was an important factor, with Navy senior enlisted demonstrating significantly higher rates and Marine junior officers with significantly lower rates of NMSKI episodes than junior enlisted. The interactions of rank by body region and rank by age on injury incidence were either non-significant or had observed singularities and were removed from the final model.

TABLE I.

Results of the Hurdle Negative Binomial Regression Assessing Body Region, Sex, Age, Rank, and Service Branch on Risk of Neuromusculoskeletal Injury Episodes (Adjusted for Zero Inflation) in USN and U.S. Marine Corps Personnel

USNU.S. Marine CorpsTotal
95% CI95% CI95% CI
RRLLULRRLLULRRLLUL
Region
Arm1.100.522.320.090.050.170.200.120.33
Elbow0.150.080.280.160.090.290.160.100.26
Knee0.600.341.060.710.411.230.660.421.04
Leg0.270.150.490.460.260.800.370.230.59
Low back1.520.882.621.630.952.801.571.002.45
Mid back0.460.260.810.490.280.840.470.300.74
Neck0.640.361.110.580.341.000.610.390.96
Pelvis-hip0.270.150.490.490.290.850.370.230.58
Shoulder0.540.300.960.640.371.110.590.380.94
Wrist and hand0.330.180.600.290.170.500.310.200.50
Other0.280.150.500.220.130.390.250.160.40
Sex
Males0.130.100.181.050.811.350.400.320.50
Age range (years)
25-342.011.253.241.380.922.051.571.082.28
35-441.971.183.301.480.942.321.661.112.50
≥451.130.651.940.450.270.750.680.441.06
Rank
Senior enlisted4.662.757.911.530.962.432.641.744.01
Junior officers0.570.321.010.370.230.590.380.250.58
Senior officers1.540.832.871.080.631.880.970.601.58
Branch
USN0.480.390.60
Negative log-likelihood1,011 on 39 df1,141 on 39 df2,284 on 41 df
USNU.S. Marine CorpsTotal
95% CI95% CI95% CI
RRLLULRRLLULRRLLUL
Region
Arm1.100.522.320.090.050.170.200.120.33
Elbow0.150.080.280.160.090.290.160.100.26
Knee0.600.341.060.710.411.230.660.421.04
Leg0.270.150.490.460.260.800.370.230.59
Low back1.520.882.621.630.952.801.571.002.45
Mid back0.460.260.810.490.280.840.470.300.74
Neck0.640.361.110.580.341.000.610.390.96
Pelvis-hip0.270.150.490.490.290.850.370.230.58
Shoulder0.540.300.960.640.371.110.590.380.94
Wrist and hand0.330.180.600.290.170.500.310.200.50
Other0.280.150.500.220.130.390.250.160.40
Sex
Males0.130.100.181.050.811.350.400.320.50
Age range (years)
25-342.011.253.241.380.922.051.571.082.28
35-441.971.183.301.480.942.321.661.112.50
≥451.130.651.940.450.270.750.680.441.06
Rank
Senior enlisted4.662.757.911.530.962.432.641.744.01
Junior officers0.570.321.010.370.230.590.380.250.58
Senior officers1.540.832.871.080.631.880.970.601.58
Branch
USN0.480.390.60
Negative log-likelihood1,011 on 39 df1,141 on 39 df2,284 on 41 df

Abbreviations: RR: relative risk; USN: U.S. Navy. Bolded figures depict statistical significance (P < .05). Reference: body region (ankle-foot), sex (female), age (18-24 years), rank (junior enlisted), and branch (U.S. Marine Corps).

TABLE I.

Results of the Hurdle Negative Binomial Regression Assessing Body Region, Sex, Age, Rank, and Service Branch on Risk of Neuromusculoskeletal Injury Episodes (Adjusted for Zero Inflation) in USN and U.S. Marine Corps Personnel

USNU.S. Marine CorpsTotal
95% CI95% CI95% CI
RRLLULRRLLULRRLLUL
Region
Arm1.100.522.320.090.050.170.200.120.33
Elbow0.150.080.280.160.090.290.160.100.26
Knee0.600.341.060.710.411.230.660.421.04
Leg0.270.150.490.460.260.800.370.230.59
Low back1.520.882.621.630.952.801.571.002.45
Mid back0.460.260.810.490.280.840.470.300.74
Neck0.640.361.110.580.341.000.610.390.96
Pelvis-hip0.270.150.490.490.290.850.370.230.58
Shoulder0.540.300.960.640.371.110.590.380.94
Wrist and hand0.330.180.600.290.170.500.310.200.50
Other0.280.150.500.220.130.390.250.160.40
Sex
Males0.130.100.181.050.811.350.400.320.50
Age range (years)
25-342.011.253.241.380.922.051.571.082.28
35-441.971.183.301.480.942.321.661.112.50
≥451.130.651.940.450.270.750.680.441.06
Rank
Senior enlisted4.662.757.911.530.962.432.641.744.01
Junior officers0.570.321.010.370.230.590.380.250.58
Senior officers1.540.832.871.080.631.880.970.601.58
Branch
USN0.480.390.60
Negative log-likelihood1,011 on 39 df1,141 on 39 df2,284 on 41 df
USNU.S. Marine CorpsTotal
95% CI95% CI95% CI
RRLLULRRLLULRRLLUL
Region
Arm1.100.522.320.090.050.170.200.120.33
Elbow0.150.080.280.160.090.290.160.100.26
Knee0.600.341.060.710.411.230.660.421.04
Leg0.270.150.490.460.260.800.370.230.59
Low back1.520.882.621.630.952.801.571.002.45
Mid back0.460.260.810.490.280.840.470.300.74
Neck0.640.361.110.580.341.000.610.390.96
Pelvis-hip0.270.150.490.490.290.850.370.230.58
Shoulder0.540.300.960.640.371.110.590.380.94
Wrist and hand0.330.180.600.290.170.500.310.200.50
Other0.280.150.500.220.130.390.250.160.40
Sex
Males0.130.100.181.050.811.350.400.320.50
Age range (years)
25-342.011.253.241.380.922.051.571.082.28
35-441.971.183.301.480.942.321.661.112.50
≥451.130.651.940.450.270.750.680.441.06
Rank
Senior enlisted4.662.757.911.530.962.432.641.744.01
Junior officers0.570.321.010.370.230.590.380.250.58
Senior officers1.540.832.871.080.631.880.970.601.58
Branch
USN0.480.390.60
Negative log-likelihood1,011 on 39 df1,141 on 39 df2,284 on 41 df

Abbreviations: RR: relative risk; USN: U.S. Navy. Bolded figures depict statistical significance (P < .05). Reference: body region (ankle-foot), sex (female), age (18-24 years), rank (junior enlisted), and branch (U.S. Marine Corps).

Factors of NMSKI-Related LIMDU and Long-Term Disability

Table II details the assessment of LIMDU in the USN, USMC, and total DoN (complete results of the hurdle model are detailed in Supplementary Table S5). The rates of LIMDU were significantly higher in knee, pelvis-hip, and shoulder conditions and lower in the mid back compared to ankle-foot conditions in the USN and USMC. Only Marines had lower rates of LIMDU because of neck, wrist-hand, and other conditions compared to ankle-foot conditions. The rate of NMSKI-related LIMDU was significantly lower in male Marines compared to their female counterparts and in Sailors aged ≥25 years and Marines aged ≥35 years compared to personnel aged 18-24 years. All junior enlisted had significantly higher rates of LIMDU compared to all other ranks.

TABLE II.

Results of the Hurdle Negative Binomial Regression Assessing Body Region, Sex, Age, Rank, and Service Branch on Risk of Neuromusculoskeletal Injury-Related Episodes of Limited Duty (Adjusted for Zero Inflation) in USN and U.S. Marine Corps Personnel

USNU.S. Marine CorpsTotal
95% CI95% CI95% CI
RRLLULRRLLULRRLLUL
Body region
Arm9.614.9918.490.840.551.294.482.896.94
Elbow0.770.451.320.770.531.120.770.541.10
Knee2.711.744.222.201.563.112.461.803.35
Leg1.050.631.720.880.611.280.960.681.35
Low back1.170.761.821.240.891.721.190.881.61
Mid back0.520.310.870.360.240.530.420.290.60
Neck0.930.591.450.670.470.960.760.561.05
Pelvis-hip2.961.894.622.821.993.982.872.103.92
Shoulder2.091.353.262.291.633.202.171.602.95
Wrist and hand0.960.601.550.560.380.820.770.551.08
Other0.910.521.570.350.230.520.600.420.87
Sex
Males0.830.671.020.670.570.780.750.650.87
Age range (years)
25-340.680.510.910.860.711.050.820.681.00
35-440.590.450.780.710.560.900.720.600.88
≥450.600.420.860.740.560.980.830.651.07
Rank
Senior enlisted0.460.350.600.750.610.930.540.470.63
Junior officers0.340.250.460.540.420.700.510.420.61
Senior officers0.330.230.480.460.340.620.380.310.47
Branch
USN0.340.270.44
Negative log-likelihood1,649 on 39 df1,466 on 39 df3,179 on 41 df
USNU.S. Marine CorpsTotal
95% CI95% CI95% CI
RRLLULRRLLULRRLLUL
Body region
Arm9.614.9918.490.840.551.294.482.896.94
Elbow0.770.451.320.770.531.120.770.541.10
Knee2.711.744.222.201.563.112.461.803.35
Leg1.050.631.720.880.611.280.960.681.35
Low back1.170.761.821.240.891.721.190.881.61
Mid back0.520.310.870.360.240.530.420.290.60
Neck0.930.591.450.670.470.960.760.561.05
Pelvis-hip2.961.894.622.821.993.982.872.103.92
Shoulder2.091.353.262.291.633.202.171.602.95
Wrist and hand0.960.601.550.560.380.820.770.551.08
Other0.910.521.570.350.230.520.600.420.87
Sex
Males0.830.671.020.670.570.780.750.650.87
Age range (years)
25-340.680.510.910.860.711.050.820.681.00
35-440.590.450.780.710.560.900.720.600.88
≥450.600.420.860.740.560.980.830.651.07
Rank
Senior enlisted0.460.350.600.750.610.930.540.470.63
Junior officers0.340.250.460.540.420.700.510.420.61
Senior officers0.330.230.480.460.340.620.380.310.47
Branch
USN0.340.270.44
Negative log-likelihood1,649 on 39 df1,466 on 39 df3,179 on 41 df

Abbreviations: RR: relative risk; USN: U.S. Navy. Bolded figures depict statistical significance (P < .05). Reference: body region (ankle-foot), sex (female), age (18-24 years), rank (junior enlisted), and branch (U.S. Marine Corps).

TABLE II.

Results of the Hurdle Negative Binomial Regression Assessing Body Region, Sex, Age, Rank, and Service Branch on Risk of Neuromusculoskeletal Injury-Related Episodes of Limited Duty (Adjusted for Zero Inflation) in USN and U.S. Marine Corps Personnel

USNU.S. Marine CorpsTotal
95% CI95% CI95% CI
RRLLULRRLLULRRLLUL
Body region
Arm9.614.9918.490.840.551.294.482.896.94
Elbow0.770.451.320.770.531.120.770.541.10
Knee2.711.744.222.201.563.112.461.803.35
Leg1.050.631.720.880.611.280.960.681.35
Low back1.170.761.821.240.891.721.190.881.61
Mid back0.520.310.870.360.240.530.420.290.60
Neck0.930.591.450.670.470.960.760.561.05
Pelvis-hip2.961.894.622.821.993.982.872.103.92
Shoulder2.091.353.262.291.633.202.171.602.95
Wrist and hand0.960.601.550.560.380.820.770.551.08
Other0.910.521.570.350.230.520.600.420.87
Sex
Males0.830.671.020.670.570.780.750.650.87
Age range (years)
25-340.680.510.910.860.711.050.820.681.00
35-440.590.450.780.710.560.900.720.600.88
≥450.600.420.860.740.560.980.830.651.07
Rank
Senior enlisted0.460.350.600.750.610.930.540.470.63
Junior officers0.340.250.460.540.420.700.510.420.61
Senior officers0.330.230.480.460.340.620.380.310.47
Branch
USN0.340.270.44
Negative log-likelihood1,649 on 39 df1,466 on 39 df3,179 on 41 df
USNU.S. Marine CorpsTotal
95% CI95% CI95% CI
RRLLULRRLLULRRLLUL
Body region
Arm9.614.9918.490.840.551.294.482.896.94
Elbow0.770.451.320.770.531.120.770.541.10
Knee2.711.744.222.201.563.112.461.803.35
Leg1.050.631.720.880.611.280.960.681.35
Low back1.170.761.821.240.891.721.190.881.61
Mid back0.520.310.870.360.240.530.420.290.60
Neck0.930.591.450.670.470.960.760.561.05
Pelvis-hip2.961.894.622.821.993.982.872.103.92
Shoulder2.091.353.262.291.633.202.171.602.95
Wrist and hand0.960.601.550.560.380.820.770.551.08
Other0.910.521.570.350.230.520.600.420.87
Sex
Males0.830.671.020.670.570.780.750.650.87
Age range (years)
25-340.680.510.910.860.711.050.820.681.00
35-440.590.450.780.710.560.900.720.600.88
≥450.600.420.860.740.560.980.830.651.07
Rank
Senior enlisted0.460.350.600.750.610.930.540.470.63
Junior officers0.340.250.460.540.420.700.510.420.61
Senior officers0.330.230.480.460.340.620.380.310.47
Branch
USN0.340.270.44
Negative log-likelihood1,649 on 39 df1,466 on 39 df3,179 on 41 df

Abbreviations: RR: relative risk; USN: U.S. Navy. Bolded figures depict statistical significance (P < .05). Reference: body region (ankle-foot), sex (female), age (18-24 years), rank (junior enlisted), and branch (U.S. Marine Corps).

Table III details the assessment of body region, sex, age, rank, and service branch on long-term disability (complete results, including the zero model of the hurdle negative binomial regression, are detailed in Supplementary Table S6). Sailors with arm, elbow, or knee conditions; Navy junior enlisted; and DoN personnel with pelvis-hip and shoulder conditions had the greatest incidence of long-term disability, with Marines having neck conditions and Sailors aged ≥45 years having the lowest rates of long-term disability. When the results of the hurdle model were compared to the standard negative binomial regression in the assessment of adjusted injury risk, the zero-adjusted model outperformed the non-adjusted model in each of the outcomes assessed (Supplementary Tables S7-S9).

TABLE III.

Results of the Hurdle Negative Binomial Regression Assessing Body Region, Sex, Age, Rank, and Service Branch on Risk of Neuromusculoskeletal Injury-Related Episodes of Long-Term Disability (Adjusted for Zero Inflation) in USN and U.S. Marine Corps Personnel

USNU.S. Marine CorpsTotal
95% CI95% CI95% CI
RRLLULRRLLULRRLLUL
Body region
Arm20.047.1656.090.430.190.9614.275.8834.61
Elbow4.561.4414.430.860.401.831.830.883.80
Knee2.571.315.041.720.973.032.101.263.50
Leg1.260.592.690.640.351.180.970.561.71
Low back1.420.742.701.360.782.371.350.832.20
Mid back1.680.694.100.600.331.090.960.531.73
Neck1.550.793.060.680.371.231.040.611.74
Pelvis-hip5.262.6710.382.951.655.263.892.326.53
Shoulder2.241.124.442.751.544.932.671.584.51
Wrist and hand1.970.944.100.660.341.291.470.832.60
Other1.570.693.540.640.321.301.000.541.86
Sex
Males0.730.531.000.390.290.520.570.450.73
Age range (years)
25-340.650.440.971.080.801.460.820.621.10
35-440.890.561.430.530.330.870.890.621.27
≥451.340.622.910.370.150.871.360.672.76
Rank
Senior enlisted0.320.220.450.820.601.110.520.410.67
Junior officers0.410.240.721.020.561.840.430.330.57
Senior officers0.160.050.490.600.321.120.580.390.87
Branch
USN0.350.190.66
Negative log-likelihood857 on 39 Df767 on 39 Df1,679 on 41 Df
USNU.S. Marine CorpsTotal
95% CI95% CI95% CI
RRLLULRRLLULRRLLUL
Body region
Arm20.047.1656.090.430.190.9614.275.8834.61
Elbow4.561.4414.430.860.401.831.830.883.80
Knee2.571.315.041.720.973.032.101.263.50
Leg1.260.592.690.640.351.180.970.561.71
Low back1.420.742.701.360.782.371.350.832.20
Mid back1.680.694.100.600.331.090.960.531.73
Neck1.550.793.060.680.371.231.040.611.74
Pelvis-hip5.262.6710.382.951.655.263.892.326.53
Shoulder2.241.124.442.751.544.932.671.584.51
Wrist and hand1.970.944.100.660.341.291.470.832.60
Other1.570.693.540.640.321.301.000.541.86
Sex
Males0.730.531.000.390.290.520.570.450.73
Age range (years)
25-340.650.440.971.080.801.460.820.621.10
35-440.890.561.430.530.330.870.890.621.27
≥451.340.622.910.370.150.871.360.672.76
Rank
Senior enlisted0.320.220.450.820.601.110.520.410.67
Junior officers0.410.240.721.020.561.840.430.330.57
Senior officers0.160.050.490.600.321.120.580.390.87
Branch
USN0.350.190.66
Negative log-likelihood857 on 39 Df767 on 39 Df1,679 on 41 Df

Abbreviations: RR: relative risk; USN: U.S. Navy. Bolded figures depict statistical significance (P < .05); Reference: body region (ankle-foot), sex (female), age (18-24 years), rank (junior enlisted), and branch (U.S. Marine Corps).

TABLE III.

Results of the Hurdle Negative Binomial Regression Assessing Body Region, Sex, Age, Rank, and Service Branch on Risk of Neuromusculoskeletal Injury-Related Episodes of Long-Term Disability (Adjusted for Zero Inflation) in USN and U.S. Marine Corps Personnel

USNU.S. Marine CorpsTotal
95% CI95% CI95% CI
RRLLULRRLLULRRLLUL
Body region
Arm20.047.1656.090.430.190.9614.275.8834.61
Elbow4.561.4414.430.860.401.831.830.883.80
Knee2.571.315.041.720.973.032.101.263.50
Leg1.260.592.690.640.351.180.970.561.71
Low back1.420.742.701.360.782.371.350.832.20
Mid back1.680.694.100.600.331.090.960.531.73
Neck1.550.793.060.680.371.231.040.611.74
Pelvis-hip5.262.6710.382.951.655.263.892.326.53
Shoulder2.241.124.442.751.544.932.671.584.51
Wrist and hand1.970.944.100.660.341.291.470.832.60
Other1.570.693.540.640.321.301.000.541.86
Sex
Males0.730.531.000.390.290.520.570.450.73
Age range (years)
25-340.650.440.971.080.801.460.820.621.10
35-440.890.561.430.530.330.870.890.621.27
≥451.340.622.910.370.150.871.360.672.76
Rank
Senior enlisted0.320.220.450.820.601.110.520.410.67
Junior officers0.410.240.721.020.561.840.430.330.57
Senior officers0.160.050.490.600.321.120.580.390.87
Branch
USN0.350.190.66
Negative log-likelihood857 on 39 Df767 on 39 Df1,679 on 41 Df
USNU.S. Marine CorpsTotal
95% CI95% CI95% CI
RRLLULRRLLULRRLLUL
Body region
Arm20.047.1656.090.430.190.9614.275.8834.61
Elbow4.561.4414.430.860.401.831.830.883.80
Knee2.571.315.041.720.973.032.101.263.50
Leg1.260.592.690.640.351.180.970.561.71
Low back1.420.742.701.360.782.371.350.832.20
Mid back1.680.694.100.600.331.090.960.531.73
Neck1.550.793.060.680.371.231.040.611.74
Pelvis-hip5.262.6710.382.951.655.263.892.326.53
Shoulder2.241.124.442.751.544.932.671.584.51
Wrist and hand1.970.944.100.660.341.291.470.832.60
Other1.570.693.540.640.321.301.000.541.86
Sex
Males0.730.531.000.390.290.520.570.450.73
Age range (years)
25-340.650.440.971.080.801.460.820.621.10
35-440.890.561.430.530.330.870.890.621.27
≥451.340.622.910.370.150.871.360.672.76
Rank
Senior enlisted0.320.220.450.820.601.110.520.410.67
Junior officers0.410.240.721.020.561.840.430.330.57
Senior officers0.160.050.490.600.321.120.580.390.87
Branch
USN0.350.190.66
Negative log-likelihood857 on 39 Df767 on 39 Df1,679 on 41 Df

Abbreviations: RR: relative risk; USN: U.S. Navy. Bolded figures depict statistical significance (P < .05); Reference: body region (ankle-foot), sex (female), age (18-24 years), rank (junior enlisted), and branch (U.S. Marine Corps).

DISCUSSION

The primary findings of this study were that during the COVID-19 pandemic, there was a significant protective effect for NMSKI that was likely a function of reduced physical exposure and limited access to nonurgent care. In the broader assessment of body region, sex, age, rank, and service branch on burden between 2016 and 2021, low back and ankle-foot were the most common of all NMSKI in the DoN, primarily affecting female personnel, service members aged 25-44 years, senior enlisted, and Marines. Among the cases of LIMDU and long-term disability, female enlisted Marines, service members aged 18-24 years, and those with shoulder, arm, pelvis-hip, and knee conditions demonstrated the greatest risk. This is the first study to the authors’ knowledge to investigate the effects of the pandemic on NMSKI in the military, as well as the burden of NMSKI-related LIMDU and disability across the USN and USMC.

NMSKI Burden

The protective effect observed during the COVID-19 pandemic is likely attributable to a shift in operations to only mission-critical tasks since infectious disease prophylaxis was prioritized for force health protection. Medical operations similarly pivoted to urgent and emergent care, with the capability to provide and document routine care in the electronic health record substantially curtailed. As such, reduced exposure to potentially injurious activities during military operations, combined with reduced routine health care utilization and documentation for NMSKI during the pandemic, is likely responsible for the reduction in reported NMSKI incidence. A similar phenomenon has been observed in NMSKI-related emergency department visits and surgical procedures in civilian hospitals.12 The rapid increase in incidence in the latter months of the study epoch corresponds with a period of resumption of regular operations, including preparation for the fall physical readiness assessment required by the USN and USMC in 2021.

The findings of increased risk of NMSKI in female service members are consistent with recent studies assessing sex within military occupation and across service branches.3–5,13 Occupational exposures differ between enlisted and officer personnel, with the former typically employed in vocational/technical fields and the latter responsible for the administration and supervision of their subordinates. This may explain the lower NMSKI risk in the officer staff. Even within the enlisted ranks, exposures change as more senior personnel assume supervisory and administrative roles. However, the substantially higher injury risk in senior enlisted personnel can be explained. It is highly conceivable that enlisted personnel who did not seek care for NMSKI earlier in their career may experience persistent symptoms, impairments, and activity limitation from chronic injury,1 conditions that require documentation for post-retirement, service-connected disability compensation. It is also plausible that social determinants of care-seeking, such as stigma, social pressures, and access to NMSKI care improved as members advanced in rank. Differences in service cultures between the USN and USMC may influence care-seeking determinants (e.g., mission over self; perceptions of weakness or NMSKI as self-limiting).1 In addition, differences in access to care (geographic and time-based barriers) and exposure based on the specific mission requirements (with vastly different occupational demands and hazards) may also explain why the USMC had a significantly larger risk compared to the USN.1 These suppositions warrant future investigation.

LIMDU and Long-Term Disability

Short-term LIMDU and long-term disability were most common in conditions of the shoulder, arm, pelvis-hip, and knee. Impairments in these regions can limit an individual’s ability to push, pull, lift, kneel, squat, and negotiate terrain, ladders, and stairs during occupational tasks. They also affect a service member’s ability to perform the annual physical fitness assessments used in the USN and USMC.14,15 These assessments encompass a timed running event (USN: 2.4 km, USMC: 4.8 km), a timed plank event (USN and USMC) or crunches (USMC), and an upper extremity endurance event (maximum push-ups [USN and USMC] or pull-ups [USMC] in 2 minutes).14,15 Since the performance of the biannual physical fitness assessment is a requirement in both the USN and USMC,14,15 it is not uncommon for clinicians to use performance on these tasks as criteria when making return to duty determination. Interestingly, service members frequently have low confidence in their ability to perform occupational and physical readiness requirements following NMSKI, with persistent pain and reinjury common in personnel returned to full duty.16 The associations between neuromusculoskeletal impairment, perceptual and psychological factors, pain, and physical activity in personnel with injuries warrant further exploration.

Female sex, junior enlisted, service members aged 18-24 years, and Marines demonstrated the greatest risk of short-term LIMDU and long-term disability. Female sex has been found to be a salient factor for activity limitations and disability in veterans, but not active duty personnel, in a study of data derived from the U.S. Census.17 This divergence may be attributed to changing roles of female personnel, with the current study epoch occurring following a policy change that opened all military occupations to women.18 Higher disability in younger, junior enlisted members is likely a function of occupation-related physical demand. Senior enlisted and officer staff members are often able to continue performing administrative duty requirements while resting, rehabilitating, and recovering from NMSKI, whereas more junior personnel with highly physical occupations often cannot. The lower rates of disability in older, more senior personnel are also likely attributed to the retirement benefit. In these groups, personnel are either emotionally or financially vested in the system (especially when they have more than 10 years of service toward the 20-year retirement) but are too young to retire. It is also plausible that the probability of experiencing a disabling condition following NMSKI that would preclude continued military service is substantially reduced because of lower exposure in the more senior personnel.

Limitations

There are limitations to this study. First, this study relied on deidentified, aggregated data provided in the NEST database. Although this allowed us to ascertain the population-level burden of NMSKI incidence, short-term LIMDU, and long-term disability, the lack of individual-level data precluded the assessment of comorbidities and other intrinsic factors that may have contributed to the outcomes under study. The NEST database only provided data from the USN and USMC, which precluded assessment of the other branches of the U.S. military and a comparison of service-related differences. Since this study relied on encounters in the electronic health record, the true burden of NMSKI (to include injuries that were self-managed) is likely much greater than that reported in this study. Lastly, the ability to evaluate the post-pandemic NMSKI burden was limited based on data availability and timing of this study. Empirical evaluation of these outcomes during the post-pandemic period is certainly warranted in future studies.

Clinical Practice and Policy Implications

From a medical planning perspective, the allocation of geographically accessible specialized care and resources needed for the timely prevention, assessment, and treatment of NMSKI should be aligned with communities that have the greatest risk.1 This evidence helps to prioritize the allocation of staffing and supplies needed to address the populations with the greatest risk of injury, short-term LIMDU, and long-term disability in order to optimize readiness. Although global pandemics are infrequent and only occur once every few generations, the release of biological agents by nations and non-state actors is a persistent threat to the readiness of the armed forces. It is paramount that military leadership considers the physical and operational readiness of personnel during periods of restricted movement and reduced social interaction (e.g., pandemics).

CONCLUSION

There was a significant protective effect during the COVID-19 pandemic that was likely a function of reduced physical exposure and access to nonurgent care. In the broader assessment, injuries to the low back and ankle-foot were the most common of all NMSKI in the DoN, affecting a higher proportion of female personnel, those aged 25-44 years, senior enlisted, and members of the USMC. Among the cases of LIMDU and long-term disability, conditions of the shoulder and arm in the upper quarter and pelvis-hip and knee in the lower quarter had the greatest burden, with personnel who were female, enlisted, aged 18-24 years, and serving in the USMC having the greatest risk.

Patient and Public Involvement

The results of this study were provided to the Office of the Medical Officer of the Marine Corps and to the U.S. Navy Bureau of Medicine and Surgery Neuromusculoskeletal Clinical Community for dissemination and translation.

ACKNOWLEDGMENT

None declared.

SUPPLEMENTARY MATERIAL

Supplementary Material is available at Military Medicine online.

FUNDING

None declared.

CONFLICT OF INTEREST STATEMENT

CDR Fraser reports grants from Defense Health Agency, Congressionally Directed Medical Research Programs, and the Office of Naval Research, outside of the submitted work. In addition, CDR Fraser has a patent pending for an Adaptive and Variable Stiffness Ankle Brace, U.S. Provisional Patent Application No. 63254,474.

DATA AVAILABILITY

Data will be available on request.

CLINICAL TRIAL REGISTRATION

Not applicable.

INSTITUTIONAL REVIEW BOARD (HUMAN SUBJECTS)

The study protocol was approved by the Naval Health Research Center’s Institutional Review Board in compliance with all applicable federal regulations governing the protection of human subjects. Research data were derived from an approved Naval Health Research Center, Institutional Review Board protocol number NHRC.2022.0201.NHSR.

INSTITUTIONAL ANIMAL CARE AND USE COMMITTEE (IACUC)

Not applicable.

INSTITUTIONAL CLEARANCE

Institutional clearance approved.

INDIVIDUAL AUTHOR CONTRIBUTION STATEMENT

J.J.F. contributed to the design, data acquisition, analysis, and drafting of the manuscript. All authors contributed to the interpretation of results, critical revision of the manuscript, and approval of the final product.

REFERENCES

1.

Fraser
JJ
,
Schmied
E
,
Rosenthal
MD
,
Davenport
TE
:
Physical therapy as a force multiplier: population health perspectives to address short-term readiness and long-term health of military service members
.
Cardiopulm Phys Ther J
2020
;
31
(
1
):
22
8
.doi: .

2.

Grimm
PD
,
Mauntel
TC
,
Potter
BK
:
Combat and noncombat musculoskeletal injuries in the US Military
.
Sports Med Arthrosc Rev
2019
;
27
(
3
):
84
91
.doi: .

3.

Fraser
JJ
,
MacGregor
AJ
,
Ryans
CP
,
Dreyer
MA
,
Gibboney
MD
,
Rhon
DI
:
Sex and occupation are salient factors associated with lateral ankle sprain risk in military tactical athletes
.
J Sci Med Sport
2021
;
24
(
7
):
677
82
.doi: .

4.

MacGregor
AJ
,
Fogleman
SA
,
Dougherty
AL
,
Ryans
CP
,
Janney
CF
,
Fraser
JJ
:
Sex differences in the incidence and risk of ankle–foot complex stress fractures among U.S. military personnel
.
J Womens Health
2021
;
31
(4):
586
92
.doi: .

5.

Glaviano
NR
,
Boling
MC
,
Fraser
JJ
:
Anterior knee pain risk in male and female military tactical athletes
.
J Athl Train
2021
;
56
(
11
):
1180
7
.doi: .

6.

Aguero
AD
,
Irrgang
JJ
,
MacGregor
AJ
,
Rothenberger
SD
,
Hart
JM
,
Fraser
JJ
:
Sex, military occupation and rank are associated with risk of anterior cruciate ligament injury in tactical-athletes
.
BMJ Mil Health
2022
.doi: .

7.

Cacian
MF
,
Saxton
A
:
COVID-19 and the military: maintaining operations while supporting civil society
. Available at https://www.csis.org/analysis/covid-19-and-military-maintaining-operations-while-supporting-civil-society; accessed
December 6, 2021
.

8.

Hauschild
V
,
Hauret
K
,
Richardson
M
,
Jones
BH
,
Lee
T
:
A Taxonomy of Injuries for Public Health Monitoring and Reporting. Addendum 1, Body Regions and Injury Types. Addendum 2, Fiscal Year 2018 Update
.
Army Public Health Center Aberdeen Proving Ground United States
;
2017
.

9.

Marrone
JB
:
Predicting 36-month attrition in the U.S. military: a comparison across service branches
.
2020
. Available at https://www.rand.org/content/dam/rand/pubs/research_reports/RR4200/RR4258/RAND_RR4258.pdf; accessed
February 11, 2021
.

10.

Feng
CX
:
A comparison of zero-inflated and hurdle models for modeling zero-inflated count data
.
J Stat Distrib Appl
2021
;
8
(
1
):
1
19
.doi: .

11.

Belmont
LPJ
,
Schoenfeld
MAJ
:
Epidemiology of combat wounds in Operation Iraqi Freedom and Operation Enduring Freedom: orthopaedic burden of disease
.
J Surg Orthop Adv
2010
;
19
(
1
):
2
7
.

12.

Crenn
V
,
El Kinani
M
,
Pietu
G
, et al. :
Impact of the COVID-19 lockdown period on adult musculoskeletal injuries and surgical management: a retrospective monocentric study
.
Sci Rep
2020
;
10
(
1
): 22442.doi: .

13.

Rhon
DI
,
Teyhen
DS
,
Kiesel
K
, et al. :
Recovery, rehabilitation, and return to full duty in a military population after a recent injury: differences between lower-extremity and spine injuries
.
Arthrosc Sports Med Rehabil
2022
;
4
(
1
):
e17
27
.doi: .

14.

Office of the Chief of Naval Operations
:
Physical readiness program (OPNAVINST 6110.1J)
.
2019
. Available at https://www.secnav.navy.mil/doni/secnav.aspx; accessed
April 14, 2022
.

15.

Office of the Commandant of the Marine Corps
:
Marine Corps physical fitness and combat fitness tests (PFT/CFT) (Marine Corps Order 6100.13A CH-2)
.
2019
, 76. Available at https://www.fitness.marines.mil/; accessed
April 14, 2022
.

16.

Rhon
DI
,
Molloy
JM
,
Monnier
A
,
Hando
BR
,
Newman
PM
:
Much work remains to reach consensus on musculoskeletal injury risk in military service members: a systematic review with meta-analysis
.
Eur J Sport Sci
2021
;
22
(1):
16
34
.doi: .

17.

Wilmoth
JM
,
London
AS
,
Parker
WM
:
Sex differences in the relationship between military service status and functional limitations and disabilities
.
Popul Res Policy Rev
2011
;
30
(
3
):
333
54
.doi: .

18.

Dempsey
ME
,
Panetta
LE
:
Elimination of the 1994 Direct Ground Combat Definition and Assignment Rule. Memorandum for Secretaries of the Military Departments, Acting Under Secretary of defense for Personnel and Readiness, and Chiefs of the Military Services
.
Joint Chiefs of Staff
;
2013
.

Author notes

The study was archived as a preprint on medRxiv and is available at doi: 10.1101/2022.06.09.22276233

The authors are military service members or employees of the U.S. Government. This work was prepared as part of our official duties. Title 17 of U.S.C. §105 provides that copyright protection under this title is not available for any work of the U.S. Government. Title 17 of U.S.C. §101 defines a U.S. Government work as work prepared by a military service member or employee of the U.S. Government as part of that person’s official duties. The views expressed in this article are those of the authors and do not necessarily reflect the official policy or position of the Department of the Navy, DoD, or the U.S. Government.

This work is written by (a) US Government employee(s) and is in the public domain in the US.

Supplementary data