
==== Front
Ann Biomed Eng
Ann Biomed Eng
Annals of Biomedical Engineering
0090-6964
1573-9686
Springer International Publishing Cham

38955890
3550
10.1007/s10439-024-03550-9
S.I. : Concussions II
Instrumented Mouthguard Decoupling Affects Measured Head Kinematic Accuracy
http://orcid.org/0000-0001-7092-2527
Gellner Ryan A. gryan3@vt.edu

Begonia Mark T.
Wood Matthew
Rockwell Lewis
Geiman Taylor
Jung Caitlyn
Rowson Steve
grid.438526.e 0000 0001 0694 4940 Virginia Tech (Biomedical Engineering and Mechanics), Blacksburg, VA USA
Associate Editor Stefan M. Duma oversaw the review of this article.

2 7 2024
2 7 2024
2024
52 10 28542871
20 2 2024
20 5 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by/4.0/ Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Many recent studies have used boil-and-bite style instrumented mouthguards to measure head kinematics during impact in sports. Instrumented mouthguards promise greater accuracy than their predecessors because of their superior ability to couple directly to the skull. These mouthguards have been validated in the lab and on the field, but little is known about the effects of decoupling during impact. Decoupling can occur for various reasons, such as poor initial fit, wear-and-tear, or excessive impact forces. To understand how decoupling influences measured kinematic error, we fit a boil-and-bite instrumented mouthguard to a 3D-printed dentition mounted to a National Operating Committee on Standards for Athletic Equipment (NOCSAE) headform. We also instrumented the headform with linear accelerometers and angular rate sensors at its center of gravity (CG). We performed a series of pendulum impact tests, varying impactor face and impact direction. We measured linear acceleration and angular velocity, and we calculated angular acceleration from the mouthguard and the headform CG. We created decoupling conditions by varying the gap between the lower jaw and the bottom face of the mouthguard. We tested three gap conditions: 0 mm (control), 1.6 mm, and 4.8 mm. Mouthguard measurements were transformed to the CG and compared to the reference measurements. We found that gap condition, impact duration, and impact direction significantly influenced mouthguard measurement error. Error was higher for larger gaps and in frontal (front and front boss) conditions. Higher errors were also found in padded conditions, but the mouthguards did not collect all rigid impacts due to inherent limitations. We present characteristic decoupling time history curves for each kinematic measurement. Exemplary frequency spectra indicating characteristic decoupling frequencies are also described. Researchers using boil-and-bite instrumented mouthguards should be aware of their limitations when interpreting results and should seek to address decoupling through advanced post-processing techniques when possible.

Keywords

Instrumented mouthguard
Boil-and-bite
Decoupling
Measurement error
Head
Impact
Kinematics
issue-copyright-statement© Biomedical Engineering Society 2024
==== Body
pmcIntroduction

Head impact exposure in living volunteers has been reported in the literature for years [1–7]. Head impact sensors characterize both impact frequency and magnitude and have been used primarily in athlete and military populations [8]. Mouthguard- and retainer-style sensors have become increasingly popular because of their potential to couple directly to the upper dentition and, hence, the skull. Some sensors perform well in the lab under ideal conditions, where a tight fit with the lower dentition clamping the mouthguard in place is often used [9–12]. The use of these sensors by volunteers in the field likely differs from the methods used to validate them in a lab for reasons such as sub-optimal fit and the potential for users to have unclenched jaws during a measured head impact [11, 13, 14].

Mouthpieces with good retention to the upper dentition may remain locked to the teeth during a measured acceleration event. Still, the forces associated with an impact can be large enough to overcome this coupling force, leading to decoupling of the mouthguard from the teeth. In addition, mouthpiece wear and tear may lower protective effects [15] and teeth retention force over time. Because these devices are often used in fast-paced environments, users may communicate with teammates or fail to anticipate impacts, inevitably resulting in some impacts occurring while users’ mouths are open. Impacts occurring concurrently with users’ mouths being open can lead to data artifacts [11, 16] or complete decoupling. The effects of mouthguard decoupling are challenging to predict in advance and detrimental to research studies using the devices due to data loss, poor data quality, and unexplainable data artifacts. Some mouthpiece sensors include algorithm- or hardware-based switches that disable data recording completely if the mouthpiece is not determined to be in a mouth [13, 17]. These switches may cause the mouthguard to forego measurement if enough relative movement occurs before or during an impact, possibly missing critical head impacts associated with injury. If data are captured, noise artifacts may cause erroneously high measurements in linear or angular kinematics as the mouthguard’s small mass and moment of inertia relative to the head encourage higher accelerations when not rigidly attached to the teeth [18, 19].

If researchers can identify mouthguard decoupling in field impacts, the next hurdle is finding the best method of dealing with them. Strict removal of affected impacts is the most straightforward answer, but this reduces sample size, introduces sampling biases, and decreases the chance of catching the most critical impacts—those rare events associated with head injury [18]. Thus, the ability to salvage affected impacts would be the most valuable method of dealing with decoupling events. While this may not be possible for all decoupling events, minimizing error in as many impacts as possible is crucial to understanding head impact biomechanics and injury tolerances in human volunteer studies.

This study aimed to determine the effects of instrumented mouthpiece decoupling on measurement error, as compared to a ground truth measurement placed at a headform center of gravity (CG). We also sought to quantify signal characteristics associated with decoupling. Here, we define decoupling as the mouthpiece sensor moving relative to the teeth. The hypothesis was that the mouthguard would become less accurate with an increased gap from the lower jaw to the mouthguard, which corresponded to increasing movement relative to the headform’s dentition.

Methods

Laboratory head impact tests were conducted using a pendulum impactor and a padded impactor face. A medium National Operating Committee on Standards for Athletic Equipment (NOCSAE) headform housed laboratory-grade instrumentation at its CG and a commercially available instrumented mouthguard on its teeth. We obtained raw mouthguard signals, which we then filtered and transformed to the headform CG for comparison. We calculated error by comparing the peak values of the CG-transformed linear accelerations to those measured directly at the CG. Measurements collected at the CG were also inversely transformed to the teeth location to generate relative acceleration signals to identify decoupling. Decoupling signals were calculated as the difference in kinematic signal traces. We rotated all measurements to match the coordinate system designated by the Society of Automotive Engineers (SAE) with the x-axis in the posterior-anterior direction, y-axis in the medial-lateral direction, and z-axis in the superior-inferior direction [20].

We used a pendulum to impact the bare medium NOCSAE headform mounted on a 50th percentile male Hybrid III neck [21, 22], which was free to slide along a single axis of a 5-degree-of-freedom linear slide table (Fig. 1A). Via a headform modification, 3D-printed dental arches were mounted in an anatomically correct location for mouthguard testing (Fig. 1B) [1]. We instrumented the headform with three linear accelerometers (Endevco 7264b-2000) and a triaxial angular rate sensor (DTS ARS3 Pro 18k) at the headform’s CG. Headform data were collected at 20 kHz for 150 ms, including 50 ms pre-trigger and 100 ms post-trigger. Data collection was triggered when the dominant axis exceeded 5 g (e.g., the x-axis for Front location impacts). The sampling rate of the lab instrumentation was 20 kHz, and the data collection system applied a hardware anti-aliasing filter with a cutoff frequency of 4 kHz to all channels before any other operations.Fig. 1 a Gravity-driven pendulum impactor used to accelerate headform. b NOCSAE headform with 3D-printed dentition and boil and bite mouthguard in place (lower dentition plate not shown)

We modified the NOCSAE headform so that 3D-printed dental arches could be mounted in an anatomical location for this test series [9]. Relative to the headform CG, the mouthguard instrumentation was at + 82 mm along the x-axis, − 9 mm along the y-axis, and + 65 mm along the z-axis. We fitted a boil-and-bite style mouthguard (Prevent Biometrics, Edina, MN) to the dentition. The manufacturer’s instructions required boiling the mouthguard for 20 s, followed by placing the mouthguard on the teeth. The user is then told to bite under suction for 30 s, followed by biting for an additional 30 s. We had to adapt the instructions for our dentition, as suction could not be applied. Instead, we used multiple clamps to provide force around the perimeter of the mouthguard. Five spring clamps and two trigger clamps provided clamping force (Pony Jorgensen, Saddle Brook, NJ—Fig. 2). We used the manufacturer's recommended times for boiling and fitting. We confirmed that fit was acceptable per the manufacturer’s standards by ensuring the mouthguard passed the “open mouth” test—it remained on the teeth under an open-mouth condition by pressing it onto the teeth and then letting go. The mouthguard did not fall off the teeth when no external pressure was applied. During testing, we placed an aluminum plate at varying distances from the bottom surface of the mouthguard to create a “closed mouth” condition and two “open mouth” conditions.Fig. 2 Spring clamps provided “bite” force for fitting the boil-and-bite mouthguard to the 3D-printed dentition

The mouthguard was instrumented with a three-axis accelerometer and a three-axis gyroscope, each sampling at 3.2 kHz. For the given sampling rate, the linear accelerometer had a bandwidth of 1600 Hz, while the gyroscope had a bandwidth of 890 Hz—both above the frequencies of interest [23]. Mouthguard data collection was triggered when any one accelerometer channel exceeded 8 g. The measurement window included 50 ms of data, including 10 ms of pre-trigger data and 40 ms of post-trigger data. The mouthguard output was tailored such that all raw data were exported, enabling the authors to do the signal processing entirely. No data post-processing techniques typically implemented within the mouthguard manufacturer’s data download system were employed for this study.

The pendulum impactor was covered with a vinyl nitrile foam pad (VN600, Dertex Corporation) or a rigid nylon face. The padded impactor face generated a resultant linear acceleration duration of approximately 10–12 ms, while the rigid face generated durations of approximately 3–5 ms. Four locations were impacted on the headform (front, front boss, rear, and rear boss) at four different pendulum angles, as described in detail in previous studies [9]. These pendulum angles corresponded to four target peak resultant linear head acceleration severities (25, 50, 75, and 100 g). Peak linear acceleration targets were said to have been achieved if the impact resulted in a peak resultant linear acceleration at the headform CG of plus or minus 10% of the target linear acceleration. We tested three decoupling conditions for each combination of location and severity: no gap, slight gap, and large gap. The gap conditions were introduced by tightening the bottom aluminum dentition to the lower face of zero, one, or three 1.6 mm shims stacked on top of one another, resulting in a no-gap condition of 0 mm or a gap condition of either 1.6 mm or 4.8 mm (Fig. 3). The shims were removed after setting the gaps. This setup created a preset, fixed opening between the mouthguard and the lower dentition in the gap conditions. The mouthguard stayed attached to the upper dentition until the time of impact in all tests, as verified by high-speed video. Each configuration underwent two trials for a total of 192 tests.Fig. 3 Lower dentition to mouthguard gap conditions from left to right: 0 mm, 1.6 mm, and 4.8 mm

Some rigid impacts were not adequately collected or were missed entirely by the mouthguard sensors. This data loss could have been due to inherent limitations with the mouthguard instrumentation (e.g., sampling rate) or due to data rejection algorithms within the mouthguard firmware (version 2.0.14). Sometimes, the mouthguard missed the primary impact but collected a secondary impact (e.g., when the headform carrier slid into the spring stopper at the end of the track). To ensure these data artifacts did not influence the results, we excluded impacts wherein the mouthguard missed the primary impact (32 impacts, 17%; Table 1).Table 1 Number of trials collected in each impact configuration

Impactor face	Impact direction	Gap condition (mm)	25 g target linear accel	50 g target linear accel	75 g target linear accel	100 g target linear accel	
Pad	Front	0	2	2	2	2	
1.6	2	2	2	2	
4.8	2	2	2	2	
Front boss	0	2	2	2	2	
1.6	2	2	2	2	
4.8	2	2	2	2	
Rear	0	2	2	2	2	
1.6	2	2	2	2	
4.8	2	2	2	2	
Rear boss	0	2	2	2	2	
1.6	2	2	2	2	
4.8	2	2	2	2	
Rigid	Front	0	2	2	1	2	
1.6	2	1	2	1	
4.8	2	1	0	0	
Front boss	0	2	0	0	0	
1.6	2	0	0	0	
4.8	2	0	0	0	
Rear	0	2	2	2	2	
1.6	2	2	2	2	
4.8	2	2	2	2	
Rear boss	0	1	2	2	2	
1.6	1	2	2	0	
4.8	1	2	2	1	
Bold cells denote one of two trials were missed. Bold italics cells denote two of two trials were missed

Raw data output by the mouthguard sensors were filtered using a second-order phaseless Butterworth dual-pass filter [24, 25] and transformed to the headform CG using rigid body dynamics equations (Eq. 1). Angular acceleration was calculated via a five-point stencil from filtered angular velocity measurements [26]. Based on previous analyses, filter cutoff frequencies (− 3 dB point) were chosen for the mouthguard data: in padded impacts, 100 Hz for linear accelerometers and 175 Hz for angular rate sensors; in rigid impacts, 175 Hz for linear accelerometers and 275 Hz for angular rate sensors. These filters best match ground truth CG measurements filtered to SAE J211-1 recommended values [27]. Data collected at the headform CG were filtered using SAE J211 recommended filter cutoff frequencies: 1650 Hz, equivalent to CFC1000, was used for linear accelerometers and 300 Hz, equivalent to CFC180, was used for angular rate sensors [20, 28].1 aCG→=aP→+ω→×ω→×r→+α→×r→

where aCG is linear acceleration at the CG, aP is linear acceleration at the mouthguard, r is the vector distance from point CG to point P, ω is angular velocity, and α is angular acceleration of the rigid body.

Error was computed as the difference between the peak values of the mouthguard signals and the headform CG measurements for both kinematic measurements (linear acceleration and angular velocity) and the calculated value of angular acceleration. We computed percent error by dividing magnitude error by the CG measurement (Eq. 2). We analyzed error with mouthguard peak values reported from the entire time window available from the mouthguard traces (50 ms), as this peak value is what would be reported in the field for a given impact. We noticed that relative differences between the reference and mouthguard signals occurred later in the data collection window; therefore, we repeated our analyses using mouthguard-reported peak values found only within the impact duration. Impact duration was defined as 12 ms for padded impacts and 5.5 ms for rigid impacts, to account for variation in linear acceleration duration by location. We searched the impact from − 1 to 5.5 ms for rigid impacts and − 1 ms to 12 ms for padded impacts, relative to the trigger at time 0 ms, to account for slight variations in trigger and peak timing within impacts.2 ∈%=PKMG-PKCGPKCG×100%

where ∈% is percent error, PKMG is the peak of the kinematic signal from the mouthguard measurement, and PKCG is the peak of the kinematic signal from the measurement at the CG.

To test our hypothesis that the mouthguard would become less accurate with an increasing gap, we used a linear model, performed an ANOVA on the model, and also tested contrasts between factor variables using least squares means. We included gap condition, target linear acceleration (impact severity), impactor face (impact duration), and general location (frontal: front and front boss locations, or rearward: rear and rear boss locations) as independent variables in our linear model. The response variables were the percent errors in peak kinematics from each of the three kinematic variables. We repeated this using the duration-windowed dataset, for a total of six linear models.

Additionally, we used weighted least product regression (WLP) analysis with bootstrapped 95% confidence intervals to determine whether fixed bias, proportional bias, or both were present with each mouthguard condition relative to CG measurements [29, 30]. Least product regression considers error in both x- and y-values when calculating the line of best fit. We used least product regression because although the laboratory measurement was considered the gold standard, some variance exists in the laboratory measurements. We used weighted regression (rather than ordinary) to account for the observed heteroscedasticity. We bootstrapped the data to 10,000 replicates to generate 95% confidence intervals. Fixed bias occurs when the 95% confidence interval for the WLP regression line’s intercept does not include zero: this means the measurements are consistently higher or lower than the reference value. Proportional bias measures for the inclusion of one in the 95% confidence interval for the WLP line’s slope. If one is not included, the measurements are proportionally higher or lower than the reference value. Both can be present in the same dataset.

Headform data collected at the CG were also inversely transformed to the headform dentition location using rigid body dynamics. This enabled us to calculate relative acceleration of the mouthguard with respect to the teeth. We filtered raw laboratory data using CFC180 for both linear acceleration and angular velocity for maximum accuracy in transformation-based 6DOF measurements [25]. We downsampled laboratory data (MATLAB: resample) to match the mouthguard sampling rate, and the mouthguard and laboratory signals were time-aligned at the point of maximum correlation using cross-correlation (MATLAB: xcorr). We visually confirmed alignment for all signal pairs. Laboratory and mouthguard data were filtered to CFC 180 (− 3 dB point of 300 Hz). We chose this cutoff frequency because it was the optimal for transforming lab measurements to the teeth and we wanted frequency attenuation to be equal across measurement devices [25]. We calculated relative linear acceleration and relative angular velocity as the differences between mouthguard and headform measurements at the teeth location. We calculated relative angular acceleration as the difference between the derived angular acceleration signals from each sensor package. Individual axis measurement and the resultant signals were time-aligned and subtracted (mouthguard minus CG) to generate relative motion along each axis.

We transformed relative motion time history curves to the frequency domain using a Fast Fourier Transform (FFT, MATLAB: fft). Relative motion curves were zero-padded to five times their original length to increase visual frequency resolution. An FFT was run on the zero-padded curves to generate frequency domain information for each curve. A single-sided FFT, up to the Nyquist frequency (1600 Hz), was used. We calculated amplitude by dividing the FFT output by the length of the input signal, taking the absolute value, and multiplying by two.

We completed all data processing in MATLAB R2023a and R2023b (Mathworks—Natwick, MA). We completed the statistical analysis (including bootstrapping) and visualization using RStudio 2023.06.0 (R Foundation for Statistical Computing, Vienna, Austria). We used a predetermined level of p = 0.05 for significance in our linear models.

Results

The major factors influencing mouthguard measurement error were gap size, impact direction, and impactor face. Larger gaps in frontal impacts (front and front boss) generally produced greater errors (Fig. 4). Rigid impacts had less error than padded impacts, but this was likely influenced by the missing rigid impacts, which were generally frontal impacts of higher severity. Fixed and proportional bias were found in mouthguard measurements; these generally decreased in severity or were eliminated with less decoupling (smaller gap) and when the analysis was limited to the impact duration only. Lastly, decoupling signals were characterized by their time- and frequency-signals. Specific axes were shown to have higher frequency amplitudes in relative motion signals when decoupling was present.Fig. 4 Probability distributions of percent error for each kinematic by impact location and duration. Differences in error by impact location were observed regardless of gap size, with differences becoming more pronounced with increases in gap size. Front and front boss impacts tended to exhibit wider error distributions than rear and rear boss impacts. Windowing the signals to the impact duration resulted in tighter error distributions

Average Error in Full Measurement Duration Analysis

We first quantified measurement error and bias for peak kinematics by using the entire duration of measurement reported by the mouthguard (10 ms pre-trigger to 40 ms post-trigger) because in typical use, instrumented mouthguards report peak kinematics over the entire collection window. We combined the front and front boss locations (Frontal) and rear and rear boss locations (Rearward) because these impact types had similar error trends. In general, mean error and variance increased with gap size for each kinematic measure (Table 2).Table 2 Differences in percent error in peak resultant kinematics for full measurement window

Impactor face	Impact direction	Gap condition (mm)	Peak linear accel	Peak angular velocity	Peak angular accel	
Percent error (%)* †‡	Percent error (%)*†‡	Percent error (%)*†‡◊	
Pad	Frontal	0	1.9 ± 4.8	− 5.0 ± 4.7	− 8.1 ± 10.2	
1.6	54.8 ± 29.9	11.2 ± 17.7	116.0 ± 95.9	
4.8	140.0 ± 105.0	21.3 ± 21.8	332.0 ± 213.0	
Rearward	0	0.4 ± 4.0	− 1.7 ± 4.4	− 5.3 ± 3.5	
1.6	− 0.1 ± 4.7	− 0.8 ± 4.2	10.6 ± 30.7	
4.8	23.3 ± 41.4	− 0.8 ± 2.9	82.5 ± 164.0	
Rigid	Frontal	0	6.2 ± 9.5	− 0.4 ± 3.2	24.6 ± 8.8	
1.6	28.6 ± 52.1	43.5 ± 64.9	124.0 ± 141.0	
4.8	26.6 ± 48.9	28.1 ± 39.7	91.5 ± 119.4	
Rearward	0	− 2.5 ± 7.1	1.7 ± 3.7	4.2 ± 4.3	
1.6	4.9 ± 4.4	9.1 ± 15.5	17.4 ± 13.3	
4.8	− 4.9 ± 9.8	2.1 ± 7.5	7.05 ± 38.6	
ANOVA results: *significantly different by gap size

†Significantly different by impact direction

‡Significantly different by impactor face

◊Significantly different by impact severity

Percent error in peak linear acceleration was significantly influenced by the gap size (p < 0.001), impactor face (p = 0.001), and impact direction (p < 0.001). Each 1 mm increase in gap size was associated with 11% greater error. Padded impacts had an average of 25% greater error than rigid. Rearward impacts generated 46% less error on average.

The model for peak angular velocity percent error showed significant influence from gap size (p = 0.008), impactor face (p = 0.019), and impact direction (p < 0.001). A gap size increase of 1 mm corresponded to a 2% increase in error for peak angular velocity. Padded impacts generated 7.9% lower error than rigid impacts, on average. Impacts to the front portion of the headform (front and front boss) had 13% higher error than rearward directed (rear and rear boss) impacts.

Finally, the linear model with angular acceleration percent error as its response variable showed significance for gap size (p < 0.001), impactor face (p = 0.001), impact direction (p < 0.001), and impact severity (p = 0.005). An increase in gap size by 1 mm generated an average of 29% greater error in angular acceleration. Padded impacts generated 37% greater error than rigid impacts, on average. Impact direction had the greatest effect on error, with impacts to the front or front boss locations generating 107% greater error than those to the rearward locations on average. Every 1 g increase in impact severity corresponded to a 0.9% increase in peak angular acceleration error.

Average Error in Impact Duration-Windowed Analysis

Measurement error and bias were also quantified for peak kinematics during the duration of impact only. Impact duration was defined as − 1 to 5.5 ms for rigid impacts and − 1 ms to 12 ms for padded impacts, relative to trigger. We chose these times so that impact-relevant peaks were captured (see example time history plots, Fig. C1). Impact-relevant peaks were defined based on the lab instrumentation measurements of linear acceleration for each impact because the lab instrumentation was not subject to decoupling conditions. Linear acceleration was chosen because it was the best indicator of the duration of force input to the head.

After windowing the signals to include only the impact duration, mean error and variance did not always increase with gap size (Table 3). In general, the magnitudes of errors were much lower after windowing in the full measurement duration analysis. Resultant peak linear acceleration percent error was significantly associated with impactor face (p = 0.003) and the impact direction (p < 0.001), but not gap size (p = 0.619) or impact severity (p = 0.91). Padded impacts had, on average, 9% greater error in peak resultant linear acceleration. Frontal impacts generated 17% higher error per impact than rearward.Table 3 Differences in percent error in peak resultant kinematics for impact duration window

Impactor face	Impact direction	Gap condition (mm)	Peak linear accel
Percent error (%)†‡	Peak angular velocity
Percent error (%)†	Peak angular accel
Percent error (%)
*†◊	
Pad	Frontal	0	1.9 ± 4.8	2.6 ± 4.0	− 3.9 ± 6.2	
1.6	44.8 ± 25.1	11.4 ± 9.7	87.6 ± 58.8	
4.8	14.7 ± 20.7	10.8 ± 18.5	63.6 ± 87.2	
Rearward	0	0.4 ± 3.7	6.4 ± 4.4	− 5.3 ± 3.5	
1.6	− 0.1 ± 4.7	5.9 ± 4.4	− 1.2 ± 19.4	
4.8	− 1.3 ± 4.5	5.1 ± 3.7	− 8.1 ± 3.9	
Rigid	Frontal	0	6.2 ± 9.5	5.7 ± 5.0	24.6 ± 8.8	
1.6	2.6 ± 22.7	19.7 ± 10.6	86.4 ± 81.3	
4.8	26.6 ± 48.9	43.9 ± 60.0	90.1 ± 120.6	
Rearward	0	− 2.5 ± 7.1	3.7 ± 2.3	4.2 ± 4.3	
1.6	4.9 ± 4.4	8.0 ± 8.3	17.4 ± 13.3	
4.8	− 7.1 ± 9.4	− 1.5 ± 4.1	− 6.8 ± 7.8	
ANOVA results: *significantly different by gap size

†Significantly different by impact direction

‡Significantly different by impactor face

◊Significantly different by impact severity

In the impact duration-windowed model, peak angular velocity percent error had even fewer significant factors than linear acceleration: only impact location (p < 0.001) was significantly associated with error in angular velocity. Frontal impacts generated 7.6% greater error than rearward impacts for peak angular velocity.

The final linear model revealed significant association between percent error in peak angular acceleration and gap size (p = 0.028), impact direction (p < 0.001), and impact severity (p = 0.001). Every 1 mm increment in gap size corresponded to a 4% increase in peak angular acceleration error. Frontal impacts generated 56% higher errors, on average. And each 1 g increase in impact severity was associated with 0.43% greater error in peak angular acceleration.

Weighted Least Products Bias

Weighted least product (WLP) regression revealed fixed and proportional bias among mouthguard measurements relative to headform measurements for the same impacts, as determined by the 95% confidence intervals (Fig. 5). Interestingly, both proportional bias (slope) and fixed bias (intercept) severity decreased in most cases when changing from the full impact duration to the impact-windowed analysis, even if the WLP confidence intervals did not show that the bias disappeared with this change. Severity of bias was also generally greater for the 1.6 mm and 4.8 mm gap conditions than for the no gap condition.Fig. 5 Weighted least product slope and intercept estimates and bootstrapped confidence intervals. Gap conditions generally had slope values further from one and intercept values further from zero. Limiting the analysis window to the impact duration (orange dots) generally produced slopes closer to one and intercepts closer to zero, especially in the 1.6 mm and 4.8 mm gap conditions. Filled points indicate statistically significant proportional (slope) or fixed (intercept) bias

Over the full measurement duration window, padded impacts displayed proportional bias (slope ≠ 1) for all gap sizes except for peak linear acceleration (PLA) at 0 mm gap, peak angular velocity (PAV) at 1.6 mm gap, and peak angular acceleration (PAA) at 0 mm gap. Fixed bias existed (intercept ≠ 0) in the largest gap condition in PLA and PAA, but in 0 mm and 1.6 mm gap condition for PAV. Rigid impacts displayed proportional bias only in angular measures: 0 mm and 1.6 mm in PAV, and 0 mm and 1.6 mm in PAA. Fixed bias existed for only 0 mm and 1.6 mm conditions in PAV for rigid impacts.

In the impact duration window, padded impacts displayed proportional bias for all conditions except for PLA 0 mm gap and PAA 0 mm gap. Fixed bias existed in the 1.6 mm and 4.8 mm conditions for PLA, but only in the 0 mm condition for PAA; PAV displayed no fixed bias for padded impacts in the impact duration windowed data. Rigid impacts displayed proportional bias for only angular measures: 0 mm and 1.6 mm in PAV and PAA. Fixed bias existed for no conditions when limiting the analysis to the duration window in rigid impacts. Full scatter plots with WLP lines and confidence intervals can be seen in Appendix 1.

Decoupling Signatures

Next, we analyzed decoupling signal characteristics by subtracting optimally aligned headform signals from mouthguard signals. Differences between measurement signals revealed relative motion artifacts in the mouthguard signals for gap conditions. These relative motion artifacts were increasingly recognizable in the mouthguard measurements as impact severity and lower dentition gap increased (Fig. 6) and were observed in both linear and angular measurements.Fig. 6 Exemplar average relative linear acceleration time traces at the teeth for each gap condition in padded impacts. Most of the relative motion occurred after the impact event was over (indicated by the gray box). Front and front boss impacts had considerably more relative motion, on average, than rear and rear boss impacts

Each impact severity and gap condition combination had unique relative motion signal characteristics. Most of the relative linear acceleration was isolated to after the impact had concluded (approximately 10 ms post-trigger in all impacts). Angular velocity continued to increase after this time point because, once up to speed, the head continued to rotate even after the applied force to the head (i.e., the pendulum impactor) was no longer acting on the head. We generated characteristic decoupling signature corridors by using the average and standard deviation of all impacts with non-trivial relative motion (greater than or equal to 50 g target linear acceleration) for the two gap conditions. These corridors were generated for linear and angular measures (Fig. 7).Fig. 7 Relative linear acceleration (top row) and angular velocity (bottom row) corridors at the teeth for padded (left column) and rigid (right column) impacts, averaged across both gap conditions and all severities ≥ 50 g resultant linear acceleration (non-trivial motion). Each padded corridor includes data from 12 impacts. The rigid front corridor includes 5 impacts, rigid rear includes 12 impacts, and rigid rear boss includes 9 impacts. The mouthguard did not collect any front boss impacts at ≥ 50 g target severity

Decoupling Spectra

We noted specific frequency bands associated with decoupling that were not present in tight condition impacts from relative motion data. Linear acceleration along the z-axis had consistently higher frequency content above 100 Hz for gap condition tests. Additionally, X and Y channels from angular acceleration showed higher frequency content above 100 Hz when a gap condition was tested (Fig. 8).Fig. 8 Frequency spectra of padded (left column) and rigid (right column) relative motion impact kinematics. Linear acceleration (top row) showed the highest frequency content along the Z-axis in decoupling relative to the no-gap condition. Angular acceleration (bottom row) X- and Y-axes had the highest frequency content in decoupling, compared to the 0 mm gap condition. F front, FB front boss, R rear, RB rear boss

Discussion

This study sought to determine if lower dentition gap, which elicited mouthguard decoupling, affected measurements from instrumented mouthpieces. In general, increasing gap size caused an increase in relative motion between the mouthguard and the headform teeth, which led to a corresponding increase in error for reported peak kinematics from the mouthguard. When decoupling occurred, impacts to the front and front boss locations were characterized by higher relative motion and errors than impacts to the rear and rear boss locations. Specific signal characteristics exist for the decoupling of instrumented mouthpieces, which are generally recognizable by impact location and contaminate the signal with higher frequency content as gap size and relative motion increase.

Recent studies have reported finding similar frequency content caused by decoupling events in instrumented mouthguards used in the field. Luke et al. described ice hockey and rugby field events categorized into one of four categories describing the pre- and post-impact coupling status [13]. Luke combined three categories (decoupled, decoupling, and recoupling) into a poorly coupled group and compared them to the fourth category (coupled). Poorly coupled impacts had significantly higher content in medium (50–200 Hz) and high (> 200 Hz) frequency ranges than coupled impacts for linear and angular acceleration. Our study agrees with this finding: decoupled impacts in a laboratory setting also exhibit higher frequency magnitudes in the frequency ranges described by Luke et al. Additionally, Wu et al. described instrumented mouthguard frequency response in a volunteer soccer heading experiment [14]. These impacts were under confirmed coupled conditions, as the volunteer had clenched teeth. In their study, instrumented mouthguard linear acceleration traces for these coupled impacts produced most of their frequency content below 100 Hz. Our study further corroborates Wu’s results for coupled impacts.

Kuo et al. observed the effects of mandible motion on instrumented mouthguard error [11]. Their study included a clenched jaw, an unconstrained articulating jaw, and a no-jaw condition. They used custom-fit mouthguards, which produced little or no relative motion in the no-jaw condition. The difference in relative motion findings suggests that custom-fit mouthguards (used in Kuo et al.’s study) may have better coupling and less error caused by relative motion than boil-and-bite mouthguards (used in our study). Kuo et al. also found that impact location influenced error, with impacts to the front and top of the head having greater errors than other locations. Our study shows similar location dependence. Finally, Kuo et al. identified characteristic frequency content in error signals from gyroscopes in instrumented mouthguards. Thus, their study presented characteristic frequency content from an unconstrained lower jaw moving into an instrumented mouthguard during impact, while our study presents the characteristic frequency content caused by decoupling and relative motion. Together, our studies point to some of the largest potential sources for measurement error when deploying instrumented mouthguards, namely lower jaw interference and decoupling.

Decoupling error is related to impact location, possibly because the headform is forced into the mouthguard in mostly rearward (rear and rear boss) impacts but away from the teeth in mostly frontal (front and front boss) impacts. This is especially noteworthy because a majority of impacts collected with mouthpieces in some studies are to the front of the head [31, 32]. The results of the present study agree with previous analyses that error depends on impact location and frontal impacts are generally characterized by greater error than rear impacts [12]. Interestingly, even when severe decoupling occurs, head kinematic data with errors similar to the no-gap condition can be collected in some impacts for a short period of the collection window. Average error decreases across all kinematic measures when the measurement window is limited to the impact duration in decoupling events. The windows we used (5.5 ms for rigid impacts and 12 ms for padded impacts) were not optimized, and may not represent the ideal signal window in decoupling events. Additionally, decoupling events collected in the field using similar devices may have different time windows during which the kinematics are valid.

Some level of significant proportional bias existed for most mouthpiece measurements in this study. Bias was of greater magnitude in open-mouth cases than in tight conditions, evidenced by WLP slope estimates drifting further from a value of one. Meaningful proportional bias was seen most often in open-mouth conditions. Fixed bias was of higher magnitude and more often significant in gap conditions. Padded impacts showed trends more clearly, likely because of the missing high-severity rigid impacts not in the dataset due to mouthguard recording errors, primarily frontal high-magnitude impacts (Table 1). Point estimates for WLP slope and intercept in nearly all tight conditions were closer to the target values (1 for slope or 0 for intercept) than corresponding point estimates for gap conditions. Most estimates moved closer to the target value when windowing the signals to the impact duration, suggesting that there is a window during which peak kinematics are valid, even for some severe decoupling events.

Decoupling signals were quantified as relative linear acceleration or angular velocity in the headform frame of reference at the teeth. Characteristic signals, differing by impact locations, were observed when relative motion occurred. Examining frequency-domain transformations of these signals also revealed a range of frequencies (approximately 100–500 Hz) excited excessively when relative motion occurred. Kuo et al. found that jaw motion created error frequencies with an average of 87 Hz in angular gyroscopes when compared to ATD measurements. His value is slightly lower than the range we identified as having excessive frequency content in acceleration measurements from decoupling, suggesting that jaw motion and decoupling may need to be dealt with via different data post-processing methods.

Our identified decoupling frequency range is greater than the greatest magnitude frequency observed via lab instrumentation in padded impacts from this study (approximately 50 Hz) but does overlap with peak frequency content from rigid impacts (approximately 150 Hz). In the future, these decoupling frequencies could be eliminated with more aggressive low-pass filtering, but an understanding of the effect of this filtering on reported peak kinematics, especially in high-frequency rigid impacts, must accompany this data-cleaning step to ensure peak kinematic reporting accuracy is not sacrificed in an attempt to salvage contaminated signals [27].

Despite elevated amplitudes in decoupling impacts, there was substantial overlap between decoupling and non-decoupling events in frequency content. As a supplementary analysis, we decided to generate short-time Fourier Transforms (STFT, MATLAB: stft) from the relative motion signals to explore if analyzing frequency with respect to time could better separate the decoupling signals from the tight condition. We applied a Hamming window to the zero-padded signals and allowed for 90% overlap in the time domain, then plotted signals up to the Nyquist frequency. As an exemplar impact, the most severe padded front boss impact (100 g target linear acceleration) is shown (Fig. 9). Time-frequency analysis revealed an increase in frequency magnitude toward the latter half of the data collection window in gap conditions after the linear acceleration had dropped to zero via the reference measurement at the headform CG. This frequency spike was most pronounced along the Z axis in linear acceleration signals and X and Y axes in angular signals (Figs. 12, 13, Appendix 2), which matched the pattern observed in the FFT plots. This pattern is physically consistent with what we would expect during a decoupling event: as the mouthguard is forced off the teeth it moves downward (Z axis, linear) relative to the teeth. In addition, some parts of the mouthguard would likely separate from the teeth before other parts (e.g., the mouthguard slips off the front of the teeth before the rear, or the right molars before the left), inducing a rotation about the X or Y axes. Thus, the frequency content of the relative motion signals describes decoupling distinctly, and the frequency content can be more clearly separated when considering frequency with respect to time.Fig. 9 Short-time Fourier Transform of the relative linear acceleration signals (mouthguard minus CG transformed to teeth) along each axis from the 100 g padded front boss impact. The post-impact frequency content was progressively more noticeable with increasing gap size, especially along the z-axis for linear acceleration

The presented analyses may help to explain on-field artifacts that continue to plague instrumented mouthpiece studies [8]. Any studies using mouthpieces where relative motion is possible, suspected, or not described should be interpreted with the knowledge that measurement error may affect these data. Falsely high data may exist within the dataset, and efforts should be made to remove or salvage impacts with motion artifacts. In addition, falsely low data may exist within the dataset if over-filtering is employed without first understanding the effects of filtering on accuracy. Improvements to current systems are warranted, as even one of the best performing instrumented mouthguards on the market [8, 9] reported most impacts in this study as valid measurements despite the large amounts of relative motion observed.

Our study has several limitations. In our analyses, we included only impacts that were captured and reported by the mouthguard. Therefore, this study portrays the mouthguard’s capability when it performs as intended. Some rigid impacts were missed, resulting in data missing from our intended test matrix. These false negatives could have been included as zeros in our peak kinematics estimates, resulting in even greater errors in the mouthguard estimates for rigid impacts. However, because our rigid lab impacts represent edge cases that occur infrequently in the field, we chose to drop the missed impact conditions to show the capabilities of the mouthguard in the impacts it did capture. Thresholding effects have been shown to introduce bias to the triggering of instrumented mouthguards [33, 34] and the debate is ongoing to find the most optimal triggering solution for these devices. However, we do not necessarily believe that triggering bias was the primary cause of the mouthguard missing impacts in our study because the same acceleration peaks had different outcomes in rigid impacts than in padded. Location effects may play a role, as more front impacts were missed than rear.

Additionally, we used a lower dentition-to-mouthguard gap to evoke decoupling in this study. This successfully generated a range of mouthguard motion relative to the teeth, but decoupling events in the field could differ in their genesis and affects. For example, we did not specifically test how wear-and-tear may degrade mouthguard fit and create decoupling over time. Our tested conditions likely represent an extreme case for mouthguard decoupling. We saw several mouthguards leave the mouth entirely after the data collection window was over. While we confirmed with high-speed video that these occurrences did not affect our measurements, they do point to the severity of our testing and the marginal nature of our mouthguard’s fit. We used a single, commercially available instrumented mouthguard brand. While the sensors in this mouthguard represent other commonly available brands, results could differ slightly with different brands. Finally, the loss of several rigid impacts, especially higher severity and larger gap conditions, makes patterns less recognizable in characteristic time history curves and introduces bias by location and/or severity in the rigid dataset. Still, rigid impacts exhibited relative motion errors dependent upon severity and location.

Conclusions

This study builds upon previous work to better understand mouthguard interactions that could occur in the field. We have shown an in-depth analysis of decoupling in a boil-and-bite style instrumented mouthguard and provided characteristic signatures and frequency content of relative motion for the potential identification of this phenomenon in the field. Open mouth conditions increased errors in mouthguard kinematic measurements, even in a properly fit boil-and-bite mouthguard. Impact duration and location effects were observed, pointing to conditional decoupling characteristics. All of these factors should be understood as potential limitations of field data from instrumented mouthguards and dealt with appropriately via advanced post-processing techniques when possible. Relative motion and the measurement errors that result should be avoided whenever possible by researchers by using the best obtainable coupling. When better coupling is not possible or when examining historical data, attempts should be made to further process or remove impacts in which relative motion is suspected. Both researchers and manufacturers can use these data to better deploy and interpret data from instrumented mouthguards.

Appendix 1: Weighted Least Products Scatter Plots

Weighted least products scatter plots are presented here, with lines and confidence intervals displayed with respect to a one-to-one line. Padded impacts show a clear trend of worsening bias with increasing gap size, while rigid impacts, which are sparser, show that gap conditions are generally worse than no gap (Fig. A1). When limiting the duration in which we searched for the peak (impact-windowed), we saw the fixed and proportional bias severity drop or disappear compared to the full duration analysis (Fig. A2). Differences are less obvious between gap conditions in the impact-windowed analysis, though angular acceleration still shows an obvious trend of worse bias with larger gaps.

Full Measurement Duration

See Fig. 10.Fig. 10 Mouthguard vs headform measurements for padded (left column) and rigid (right column) peak linear acceleration (top row), peak angular velocity (center row), and peak angular acceleration (bottom row), reported for the full measurement window. Solid lines are WLP estimates for slope. Dash-dotted lines are WLP 95% confidence intervals using a bootstrapping method. Gray solid lines represent a one-to-one relationship

Impact-Windowed Duration

See Fig. 11.Fig. 11 Mouthguard vs headform measurements for padded (left column) and rigid (right column) peak linear acceleration (top row), peak angular velocity (center row), and peak angular acceleration (bottom row), reported for the impact duration window. Impact durations: − 1 to 5 ms for rigid, − 1 to 12 ms for padded impacts. Solid lines are WLP estimates for slope. Dash-dotted lines are WLP 95% confidence intervals using a bootstrapping method. Gray solid lines represent a one-to-one relationship

Appendix 2: Short Time Fourier Transform of Angular Acceleration

Frequency content with respect to time may be a very helpful analysis tool in analyzing mouthguard impacts where decoupling is suspected. Our study showed post-impact frequency spikes in rotation about the X and Y axes for relative angular signals made by the mouthguard (Fig. 12). These translated into similar spikes about the same axes for angular acceleration (Fig. 13). These areas of high frequency noise are most likely caused by the mouthguard separating from the front teeth more quickly than the molars (or vice versa) and from the impact side more quickly than the non-impact side, leading to rotation about the X and Y axes that worsens with increasing gap condition.Fig. 12 Short-time Fourier Transform of the relative angular velocity signals (mouthguard minus CG transformed to teeth) along each axis from the 100 g padded front boss impact. The post-impact frequency content in the x- and y-axes was progressively of higher amplitude with increasing gap size

Fig. 13 Short-time Fourier Transform of the relative angular acceleration signals (mouthguard minus CG transformed to teeth) along each axis from the 100 g padded front boss impact. Like angular velocity, the post-impact frequency content in the x- and y-axes was progressively of higher amplitude with increasing gap size

Appendix 3: Example Time History Plots

Below, we provide example time history plots for select impact conditions (Figs. 14, 15). These plots demonstrate the relative accuracy of the mouthguard during the initial impact duration (gray boxes) and the artifacts present in the data thereafter for gap conditions. Impact duration was defined as − 1 ms to 12 ms for padded impacts and − 1 to 5.5 ms for rigid impacts.Fig. 14 Exemplar time history plots of the kinematics from the 100 g padded front boss impacts. Reference curves were recorded at the CG and transformed to the teeth. Mouthguard curves were recorded at the teeth. We observed increasing oscillatory artifacts with greater gap conditions, but most of these artifacts occurred after the impact (gray boxes)

Fig. 15 Exemplar time history plots of the kinematics from the 100 g rigid rear impacts. Reference curves were recorded at the CG and transformed to the teeth. Mouthguard curves were recorded at the teeth. We observed increasing oscillatory artifacts with greater gap conditions, but most of these artifacts occurred after the impact (gray boxes)

Declarations

Conflict of interest

The authors report no conflicts of interest.

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

1. Crisco JJ Wilcox BJ Magnitude of head impact exposures in individual collegiate football players J. Appl. Biomech. 2012 28 2 174 183 10.1123/jab.28.2.174 21911854
Crisco, J. J., B. J. Wilcox, et al. Magnitude of head impact exposures in individual collegiate football players. J. Appl. Biomech. 28(2):174–183, 2012.21911854 10.1123/jab.28.2.174
2. Mihalik JP Guskiewicz KM Head impact biomechanics in youth hockey: comparisons across playing position, event types, and impact locations Ann. Biomed. Eng. 2012 40 1 141 149 10.1007/s10439-011-0405-3 22012083
Mihalik, J. P., K. M. Guskiewicz, et al. Head impact biomechanics in youth hockey: comparisons across playing position, event types, and impact locations. Ann. Biomed. Eng. 40(1):141–149, 2012.22012083 10.1007/s10439-011-0405-3
3. Crisco JJ Wilcox BJ Head impact exposure in collegiate football players J. Biomech. 2011 44 15 2673 2678 10.1016/j.jbiomech.2011.08.003 21872862
Crisco, J. J., B. J. Wilcox, et al. Head impact exposure in collegiate football players. J. Biomech. 44(15):2673–2678, 2011.21872862 10.1016/j.jbiomech.2011.08.003
4. Rowson S Bland ML Biomechanical perspectives on concussion in sport Sports Med. Arthrosc. Rev. 2016 24 3 100 107 10.1097/JSA.0000000000000121 27482775
Rowson, S., M. L. Bland, et al. Biomechanical perspectives on concussion in sport. Sports Med. Arthrosc. Rev. 24(3):100–107, 2016.27482775 10.1097/JSA.0000000000000121
5. Campolettano ET Gellner RA Rowson S High-magnitude head impact exposure in youth football J. Neurosurg. 2017 20 6 604 612
Campolettano, E. T., R. A. Gellner, and S. Rowson. High-magnitude head impact exposure in youth football. J. Neurosurg. 20(6):604–612, 2017.
6. Press JN Rowson S Quantifying head impact exposure in collegiate women's soccer Clin. J. Sport Med. 2016 27 2 104 110 10.1097/JSM.0000000000000313
Press, J. N., and S. Rowson. Quantifying head impact exposure in collegiate women’s soccer. Clin. J. Sport Med. 27(2):104–110, 2016.10.1097/JSM.0000000000000313
7. Rowson S Duma S Six Degree of Freedom Head Acceleration Measurements in Collegiate Football 2008 Blacksburg Virginia Tech - Wake Forest, Center for Injury Biomechanics
Rowson, S., and S. Duma. Six Degree of Freedom Head Acceleration Measurements in Collegiate Football. Blacksburg: Virginia Tech - Wake Forest, Center for Injury Biomechanics, 2008.
8. Jones B Tooby J Ready for impact? A validity and feasibility study of instrumented mouthguards (iMGs) Br. J. Sports Med. 2022 56 20 1171 1179 10.1136/bjsports-2022-105523
Jones, B., J. Tooby, et al. Ready for impact? A validity and feasibility study of instrumented mouthguards (iMGs). Br. J. Sports Med. 56(20):1171–1179, 2022.10.1136/bjsports-2022-105523
9. Kieffer EE Begonia MT Tyson AM Rowson S A two-phased approach to quantifying head impact sensor accuracy: In-laboratory and on-field assessments Ann. Biomed. Eng. 2020 48 2613 2625 10.1007/s10439-020-02647-1 33051745
Kieffer, E. E., M. T. Begonia, A. M. Tyson, and S. Rowson. A two-phased approach to quantifying head impact sensor accuracy: In-laboratory and on-field assessments. Ann. Biomed. Eng. 48:2613–2625, 2020.33051745 10.1007/s10439-020-02647-1
10. Tyson AM Duma SM Rowson S Laboratory evaluation of low-cost wearable sensors for measuring head impacts in sports J. Appl. Biomech. 2018 34 4 320 326 10.1123/jab.2017-0256 29613824
Tyson, A. M., S. M. Duma, and S. Rowson. Laboratory evaluation of low-cost wearable sensors for measuring head impacts in sports. J. Appl. Biomech. 34(4):320–326, 2018.29613824 10.1123/jab.2017-0256
11. Kuo C Wu LC Effect of the mandible on mouthguard measurements of head kinematics J. Biomech. 2016 49 9 1845 1853 10.1016/j.jbiomech.2016.04.017 27155744
Kuo, C., L. C. Wu, et al. Effect of the mandible on mouthguard measurements of head kinematics. J. Biomech. 49(9):1845–1853, 2016.27155744 10.1016/j.jbiomech.2016.04.017
12. Liu Y Domel AG Validation and comparison of instrumented mouthguards for measuring head kinematics and assessing brain deformation in football impacts Ann. Biomed. Eng. 2020 48 2580 2598 10.1007/s10439-020-02629-3 32989591
Liu, Y., A. G. Domel, et al. Validation and comparison of instrumented mouthguards for measuring head kinematics and assessing brain deformation in football impacts. Ann. Biomed. Eng. 48:2580–2598, 2020.32989591 10.1007/s10439-020-02629-3
13. Luke D Kenny R Bondi D Clansey A Wu LC On-field instrumented mouthguard coupling J. Biomech. 2023 162 111889 10.1016/j.jbiomech.2023.111889 38071791
Luke, D., R. Kenny, D. Bondi, A. Clansey, and L. C. Wu. On-field instrumented mouthguard coupling. J. Biomech. 162:111889, 2023.38071791 10.1016/j.jbiomech.2023.111889
14. Wu LC Nangia V In vivo evaluation of wearable head impact sensors Ann. Biomed. Eng. 2016 44 4 1234 1245 10.1007/s10439-015-1423-3 26289941
Wu, L. C., V. Nangia, et al. In vivo evaluation of wearable head impact sensors. Ann. Biomed. Eng. 44(4):1234–1245, 2016.26289941 10.1007/s10439-015-1423-3
15. Sousa AM Pinho AC Messias A Piedade AP Present status in polymeric mouthguards. A future area for additive manufacturing? Polymers 2020 12 7 1490 10.3390/polym12071490 32635307
Sousa, A. M., A. C. Pinho, A. Messias, and A. P. Piedade. Present status in polymeric mouthguards. A future area for additive manufacturing? Polymers. 12(7):1490, 2020.32635307 10.3390/polym12071490
16. Tierney, G., S. Rowson, et al. Frequency content and filtering of head sensor kinematics: A method to enable field-based inter-study comparisons, 2023. https://arxiv.org/abs/2303.03043.
17. Wu LC Zarnescu L Nangia V Cam B Camarillo DB A head impact detection system using SVM classification and proximity sensing in an instrumented mouthguard IEEE Trans. Bio-Med. Eng. 2014 61 11 2659 2668 10.1109/TBME.2014.2320153
Wu, L. C., L. Zarnescu, V. Nangia, B. Cam, and D. B. Camarillo. A head impact detection system using SVM classification and proximity sensing in an instrumented mouthguard. IEEE Trans. Bio-Med. Eng. 61(11):2659–2668, 2014.10.1109/TBME.2014.2320153
18. Rowson S Mihalik J Consensus head acceleration measurement practices (CHAMP): Study design and statistical analysis Ann. Biomed. Eng. 2022 50 11 1346 1355 10.1007/s10439-022-03101-0 36253602
Rowson, S., J. Mihalik, et al. Consensus head acceleration measurement practices (CHAMP): Study design and statistical analysis. Ann. Biomed. Eng. 50(11):1346–1355, 2022.36253602 10.1007/s10439-022-03101-0
19. Campolettano ET Gellner RA Development of a concussion risk function for a youth population using head linear and rotational acceleration Ann. Biomed. Eng. 2020 48 92 103 10.1007/s10439-019-02382-2 31659605
Campolettano, E. T., R. A. Gellner, et al. Development of a concussion risk function for a youth population using head linear and rotational acceleration. Ann. Biomed. Eng. 48:92–103, 2020.31659605 10.1007/s10439-019-02382-2
20. SAE, S. SURFACE VEHICLE RECOMMENDED PRACTICE, in Instrumentation for Impact Test-Part 1-Electronic Instrumentation. SAE International, 2014.
21. Cobb BR Zadnik AM Rowson S Comparative analysis of helmeted impact response of Hybrid III and National Operating Committee on Standards for Athletic Equipment headforms Proc. Inst. Mech. Eng. P 2016 230 1 50 60
Cobb, B. R., A. M. Zadnik, and S. Rowson. Comparative analysis of helmeted impact response of Hybrid III and National Operating Committee on Standards for Athletic Equipment headforms. Proc. Inst. Mech. Eng. P. 230(1):50–60, 2016.
22. Rowson B Rowson S Duma SM Hockey STAR: A methodology for assessing the biomechanical performance of hockey helmets Ann. Biomed. Eng. 2015 43 10 2429 2443 10.1007/s10439-015-1278-7 25822907
Rowson, B., S. Rowson, and S. M. Duma. Hockey STAR: A methodology for assessing the biomechanical performance of hockey helmets. Ann. Biomed. Eng. 43(10):2429–2443, 2015.25822907 10.1007/s10439-015-1278-7
23. Wu LC Laksari K Bandwidth and sample rate requirements for wearable head impact sensors J. Biomech. 2016 49 13 2918 2924 10.1016/j.jbiomech.2016.07.004 27497499
Wu, L. C., K. Laksari, et al. Bandwidth and sample rate requirements for wearable head impact sensors. J. Biomech. 49(13):2918–2924, 2016.27497499 10.1016/j.jbiomech.2016.07.004
24. Winter DA Biomechanics and Motor Control of Human Movement 2009 Hoboken Wiley
Winter, D. A. Biomechanics and Motor Control of Human Movement. Hoboken: Wiley, 2009.
25. Gellner, R., M. Begonia, S. Gagliardi, G. Tierney, and S. Rowson, Optimal Cutoff Frequencies for Filtering Linear Acceleration and Angular Velocity Signals Associated with Laboratory Head Impacts Measured with Externally Mounted Sensors. Proceedings of International research conference on the biomechanics of impact, 2023.
26. SAE. Calculation Guidelines for Impact Testing J1727_201502, 2015.
27. Gellner RA Ch 2—Optimal Head Impact Signal Processing and the Description and Perception of Impact Exposure in Field Adolescent Ice Hockey Players 2023 Blacksburg Virginia Tech
Gellner, R. A. Ch 2—Optimal Head Impact Signal Processing and the Description and Perception of Impact Exposure in Field Adolescent Ice Hockey Players. Blacksburg: Virginia Tech, 2023.
28. Cobb BR Tyson AM Rowson S Head acceleration measurement techniques: Reliability of angular rate sensor data in helmeted impact testing Proc. Inst. Mech. Eng. P 2018 232 2 176 181
Cobb, B. R., A. M. Tyson, and S. Rowson. Head acceleration measurement techniques: Reliability of angular rate sensor data in helmeted impact testing. Proc. Inst. Mech. Eng. P. 232(2):176–181, 2018.
29. Ludbrook J Comparing methods of measurement Clin. Exp. Pharmacol. Physiol. 1997 24 2 198 203 10.1111/j.1440-1681.1997.tb01807.x
Ludbrook, J. Comparing methods of measurement. Clin. Exp. Pharmacol. Physiol. 24(2):198–203, 1997.10.1111/j.1440-1681.1997.tb01807.x
30. Ludbrook J Linear regression analysis for comparing two measurers or methods of measurement: but which regression? Clin. Exp. Pharmacol. Physiol. 2010 37 7 692 699 10.1111/j.1440-1681.2010.05376.x 20337658
Ludbrook, J. Linear regression analysis for comparing two measurers or methods of measurement: but which regression? Clin. Exp. Pharmacol. Physiol. 37(7):692–699, 2010.20337658 10.1111/j.1440-1681.2010.05376.x
31. Bartsch AJ Hedin D High energy side and rear American football head impacts cause obvious performance decrement on video Ann. Biomed. Eng. 2020 48 2667 2677 10.1007/s10439-020-02640-8 33111969
Bartsch, A. J., D. Hedin, et al. High energy side and rear American football head impacts cause obvious performance decrement on video. Ann. Biomed. Eng. 48:2667–2677, 2020.33111969 10.1007/s10439-020-02640-8
32. Bartsch, A. J., D. S. Hedin, et al. Laboratory and on-field data collected by a head impact monitoring mouthguard. Proceedings of 2019 41st Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC), 2019.
33. Wang T Kenny R Wu LC Head impact sensor triggering bias introduced by linear acceleration thresholding Ann. Biomed. Eng. 2021 49 3189 3199 10.1007/s10439-021-02868-y 34622314
Wang, T., R. Kenny, and L. C. Wu. Head impact sensor triggering bias introduced by linear acceleration thresholding. Ann. Biomed. Eng. 49:3189–3199, 2021.34622314 10.1007/s10439-021-02868-y
34. Tooby J Till K When to pull the trigger: conceptual considerations for approximating head acceleration events using instrumented mouthguards Sports Med. 2024 10.1007/s40279-024-02012-5 38922555
Tooby, J., K. Till, et al. When to pull the trigger: conceptual considerations for approximating head acceleration events using instrumented mouthguards. Sports Med. 2024. 10.1007/s40279-024-02012-5.38922555 10.1007/s40279-024-02012-5
