Estimated Injury Risk for Specific Injuries and Body Regions in Frontal Motor Vehicle Crashes

Objective: Injury risk curves estimate motor vehicle crash (MVC) occupant injury risk from vehicle, crash, and/or occupant factors. Many vehicles are equipped with event data recorders (EDRs) that collect data including the crash speed and restraint status during a MVC. This study's goal was to use regulation-required data elements for EDRs to compute occupant injury risk for (1) specific injuries and (2) specific body regions in frontal MVCs from weighted NASS-CDS data.

Methods: Logistic regression analysis of NASS-CDS single-impact frontal MVCs involving front seat occupants with frontal airbag deployment was used to produce 23 risk curves for specific injuries and 17 risk curves for Abbreviated Injury Scale (AIS) 2+ to 5+ body region injuries. Risk curves were produced for the following body regions: head and thorax (AIS 2+, 3+, 4+, 5+), face (AIS 2+), abdomen, spine, upper extremity, and lower extremity (AIS 2+, 3+). Injury risk with 95% confidence intervals was estimated for 15–105 km/h longitudinal delta-Vs and belt status was adjusted for as a covariate.

Results: Overall, belted occupants had lower estimated risks compared to unbelted occupants and the risk of injury increased as longitudinal delta-V increased. Belt status was a significant predictor for 13 specific injuries and all body region injuries with the exception of AIS 2+ and 3+ spine injuries. Specific injuries and body region injuries that occurred more frequently in NASS-CDS also tended to carry higher risks when evaluated at a 56 km/h longitudinal delta-V. In the belted population, injury risks that ranked in the top 33% included 4 upper extremity fractures (ulna, radius, clavicle, carpus/metacarpus), 2 lower extremity fractures (fibula, metatarsal/tarsal), and a knee sprain (2.4–4.6% risk). Unbelted injury risks ranked in the top 33% included 4 lower extremity fractures (femur, fibula, metatarsal/tarsal, patella), 2 head injuries with less than one hour or unspecified prior unconsciousness, and a lung contusion (4.6–9.9% risk). The 6 body region curves with the highest risks were for AIS 2+ lower extremity, upper extremity, thorax, and head injury and AIS 3+ lower extremity and thorax injury (15.9–43.8% risk).

Conclusions: These injury risk curves can be implemented into advanced automatic crash notification (AACN) algorithms that utilize vehicle EDR measurements to predict occupant injury immediately following a MVC. Through integration with AACN, these injury risk curves can provide emergency medical services (EMS) and other patient care providers with information on suspected occupant injuries to improve injury detection and patient triage.     

Development of a Computationally Efficient Full Human Body Finite Element Model

Introduction: A simplified and computationally efficient human body finite element model is presented. The model complements the Global Human Body Models Consortium (GHBMC) detailed 50th percentile occupant (M50-O) by providing kinematic and kinetic data with a significantly reduced run time using the same body habitus.

Methods: The simplified occupant model (M50-OS) was developed using the same source geometry as the M50-O. Though some meshed components were preserved, the total element count was reduced by remeshing, homogenizing, or in some cases omitting structures that are explicitly contained in the M50-O. Bones are included as rigid bodies, with the exception of the ribs, which are deformable but were remeshed to a coarser element density than the M50-O. Material models for all deformable components were drawn from the biomechanics literature. Kinematic joints were implemented at major articulations (shoulder, elbow, wrist, hip, knee, and ankle) with moment vs. angle relationships from the literature included for the knee and ankle. The brain of the detailed model was inserted within the skull of the simplified model, and kinematics and strain patterns are compared.

Results: The M50-OS model has 11 contacts and 354,000 elements; in contrast, the M50-O model has 447 contacts and 2.2 million elements. The model can be repositioned without requiring simulation. Thirteen validation and robustness simulations were completed. This included denuded rib compression at 7 discrete sites, 5 rigid body impacts, and one sled simulation. Denuded tests showed a good match to the experimental data of force vs. deflection slopes. The frontal rigid chest impact simulation produced a peak force and deflection within the corridor of 4.63 kN and 31.2%, respectively. Similar results vs. experimental data (peak forces of 5.19 and 8.71 kN) were found for an abdominal bar impact and lateral sled test, respectively. A lateral plate impact at 12 m/s exhibited a peak of roughly 20 kN (due to stiff foam used around the shoulder) but a more biofidelic response immediately afterward, plateauing at 9 kN at 12 ms. Results from a frontal sled simulation showed that reaction forces and kinematic trends matched experimental results well. The robustness test demonstrated that peak femur loads were nearly identical to the M50-O model. Use of the detailed model brain within the simplified model demonstrated a paradigm for using the M50-OS to leverage aspects of the M50-O. Strain patterns for the 2 models showed consistent patterns but greater strains in the detailed model, with deviations thought to be the result of slightly different kinematics between models. The M50-OS with the deformable skull and brain exhibited a run time 4.75 faster than the M50-O on the same hardware.

Conclusions: The simplified GHBMC model is intended to complement rather than replace the detailed M50-O model. It exhibited, on average, a 35-fold reduction in run time for a set of rigid impacts. The model can be used in a modular fashion with the M50-O and more broadly can be used as a platform for parametric studies or studies focused on specific body regions.     

New Methodology for an Expert-Designed Map From International Classification of Diseases (ICD) to Abbreviated Injury Scale (AIS) 3+ Severity Injury

Objective: There has been a longstanding desire for a map to convert International Classification of Diseases (ICD) injury codes to Abbreviated Injury Scale (AIS) codes to reflect the severity of those diagnoses. The Association for the Advancement of Automotive Medicine (AAAM) was tasked by European Union representatives to create a categorical map classifying diagnoses codes as serious injury (Abbreviated Injury Scale [AIS] 3+), minor/moderate injury (AIS 1/2), or indeterminate. This study's objective was to map injury-related ICD-9-CM (clinical modification) and ICD-10-CM codes to these severity categories.

Methods: Approximately 19,000 ICD codes were mapped, including injuries from the following categories: amputations, blood vessel injury, burns, crushing injury, dislocations/sprains/strains, foreign body, fractures, internal organ, nerve/spinal cord injury, intracranial, laceration, open wounds, and superficial injury/contusion. Two parallel activities were completed to create the maps: (1) An in-person expert panel and (2) an electronic survey. The panel consisted of expert users of AIS and ICD from North America, the United Kingdom, and Australia. The panel met in person for 5 days, with follow-up virtual meetings to create and revise the maps. Additional qualitative data were documented to resolve potential discrepancies in mapping. The electronic survey was completed by 95 injury coding professionals from North America, Spain, Australia, and New Zealand over 12 weeks. ICD-to-AIS maps were created for: ICD-9-CM and ICD-10-CM. Both maps indicated whether the corresponding AIS 2005/Update 2008 severity score for each ICD code was AIS 3+, 1/2, or indeterminable. Though some ICD codes could be mapped to multiple AIS codes, the maximum severity of all potentially mapped injuries determined the final severity categorization.

Results: The in-person panel consisted of 13 experts, with 11 Certified AIS specialists (CAISS) with a median of 8 years and an average of 15 years of coding experience. Consensus was reached for AIS severity categorization for all injury-related ICD codes. There were 95 survey respondents, with a median of 8 years of injury coding experience. Approximately 15 survey responses were collected per ICD code. Results from the 2 activities were compared, and any discrepancies were resolved using additional qualitative and quantitative data from the in-person panel and survey results, respectively.

Conclusions: Robust maps of ICD-9-CM and ICD-10-CM injury codes to AIS severity categories (3+ versus <3) were successfully created from an in-person panel discussion and electronic survey. These maps provide a link between the common ICD diagnostic lexicons and the AIS severity coding system and are of value to injury researchers, public health scientists, and epidemiologists using large databases without available AIS coding.     

Driver Injury Risk Variability in Finite Element Reconstructions of Crash Injury Research and Engineering Network (CIREN) Frontal Motor Vehicle Crashes

Objective: A 3-phase real-world motor vehicle crash (MVC) reconstruction method was developed to analyze injury variability as a function of precrash occupant position for 2 full-frontal Crash Injury Research and Engineering Network (CIREN) cases.

Method: Phase I: A finite element (FE) simplified vehicle model (SVM) was developed and tuned to mimic the frontal crash characteristics of the CIREN case vehicle (Camry or Cobalt) using frontal New Car Assessment Program (NCAP) crash test data. Phase II: The Toyota HUman Model for Safety (THUMS) v4.01 was positioned in 120 precrash configurations per case within the SVM. Five occupant positioning variables were varied using a Latin hypercube design of experiments: seat track position, seat back angle, D-ring height, steering column angle, and steering column telescoping position. An additional baseline simulation was performed that aimed to match the precrash occupant position documented in CIREN for each case. Phase III: FE simulations were then performed using kinematic boundary conditions from each vehicle's event data recorder (EDR). HIC15, combined thoracic index (CTI), femur forces, and strain-based injury metrics in the lung and lumbar vertebrae were evaluated to predict injury.

Results: Tuning the SVM to specific vehicle models resulted in close matches between simulated and test injury metric data, allowing the tuned SVM to be used in each case reconstruction with EDR-derived boundary conditions. Simulations with the most rearward seats and reclined seat backs had the greatest HIC15, head injury risk, CTI, and chest injury risk. Calculated injury risks for the head, chest, and femur closely correlated to the CIREN occupant injury patterns. CTI in the Camry case yielded a 54% probability of Abbreviated Injury Scale (AIS) 2+ chest injury in the baseline case simulation and ranged from 34 to 88% (mean = 61%) risk in the least and most dangerous occupant positions. The greater than 50% probability was consistent with the case occupant's AIS 2 hemomediastinum. Stress-based metrics were used to predict injury to the lower leg of the Camry case occupant. The regional-level injury metrics evaluated for the Cobalt case occupant indicated a low risk of injury; however, strain-based injury metrics better predicted pulmonary contusion. Approximately 49% of the Cobalt occupant's left lung was contused, though the baseline simulation predicted 40.5% of the lung to be injured.

Conclusions: A method to compute injury metrics and risks as functions of precrash occupant position was developed and applied to 2 CIREN MVC FE reconstructions. The reconstruction process allows for quantification of the sensitivity and uncertainty of the injury risk predictions based on occupant position to further understand important factors that lead to more severe MVC injuries.     

Comparing universal kriging and land-use regression for predicting concentrations of gaseous oxides of nitrogen (NOx) for the Multi-Ethnic Study of Atherosclerosis and Air Pollution (MESA Air)

BACKGROUND: Epidemiological studies that assess the health effects of long-term exposure to ambient air pollution are used to inform public policy. These studies rely on exposure models that use data collected from pollution monitoring sites to predict exposures at subject locations. Land use regression (LUR) and universal kriging (UK) have been suggested as potential prediction methods. We evaluate these approaches on a dataset including measurements from three seasons in Los Angeles, CA. METHODS: The measurements of gaseous oxides of nitrogen (NOx) used in this study are from a "snapshot" sampling campaign that is part of the Multi-Ethnic Study of Atherosclerosis and Air Pollution (MESA Air). The measurements in Los Angeles were collected during three two-week periods in the summer, autumn, and winter, each with about 150 sites. The design included clusters of monitors on either side of busy roads to capture near-field gradients of traffic-related pollution. LUR and UK prediction models were created using geographic information system (GIS)-based covariates. Selection of covariates was based on 10-fold cross-validated (CV) R(2) and root mean square error (RMSE). Since UK requires specialized software, a computationally simpler two-step procedure was also employed to approximate fitting the UK model using readily available regression and GIS software. RESULTS: UK models consistently performed as well as or better than the analogous LUR models. The best CV R(2) values for season-specific UK models predicting log(NOx) were 0.75, 0.72, and 0.74 (CV RMSE 0.20, 0.17, and 0.15) for summer, autumn, and winter, respectively. The best CV R(2) values for season-specific LUR models predicting log(NOx) were 0.74, 0.60, and 0.67 (CV RMSE 0.20, 0.20, and 0.17). The two-stage approximation to UK also performed better than LUR and nearly as well as the full UK model with CV R(2) values 0.75, 0.70, and 0.70 (CV RMSE 0.20, 0.17, and 0.17) for summer, autumn, and winter, respectively. CONCLUSION: High quality LUR and UK prediction models for NOx in Los Angeles were developed for the three seasons based on data collected for MESA Air. In our study, UK consistently outperformed LUR. Similarly, the 2-step approach was more effective than the LUR models, with performance equal to or slightly worse than UK.     

Diet quality affects an attractive white plumage pattern in dark-eyed juncos (Junco hyemalis)

Sexually selected traits that act as signals of quality often display some degree of condition dependence. In birds, condition dependence of ornamental plumage is often mediated by production costs related to acquisition or allocation of dietary resources. White plumage ornaments, however, have often been assumed to be inexpensive because their production requires neither pigment nor specialized feather structure. In male dark-eyed juncos (Junco hyemalis), the size of a white patch on the tail contributes to attractiveness and mating success. Using captive males, we examined the effects of diet quality on the size and brightness of the tail-white patch. After removing four tail feathers to induce replacement, we maintained subjects on a subsistence (low-protein) or enriched (high-protein) diet while induced feathers grew. Birds that received an enriched diet grew their feathers more quickly and grew larger, brighter white patches. Feather growth rate was positively correlated with the increase in the size of the tail-white patch, a relationship that was stronger in the subsistence diet group. However, within diet treatments, faster-grown feathers were slightly duller. Taken together, these results suggest that variation in diet quality may lead to condition-dependent expression of tail white and that condition dependence may be stronger in more stressful environments. We suggest a mechanism by which increased feather growth rate may lead to an increase in the size of the tail-white patch and discuss potential trade-offs between signal size and brightness.     

The ability of a benthic elasmobranch to discriminate between biological and artificial electric fields

To investigate the ability of elasmobranchs to distinguish between differing prey-type electric fields we examined the electroreceptive foraging behaviour of a model species, Scyliorhinus canicula (small-spotted catshark). Catshark preferences were studied by behaviourally conditioning them to swim through narrow tunnels, and on exit presenting them simultaneously with two different electric fields. Their subsequent choices of the following paired options were recorded; (i) Two artificial electric fields (dipole electrodes) with different magnitude direct current (D.C.), (ii) Two artificial electric fields, one D.C. and the other alternating current (A.C.), of the same magnitude, and (iii) similar magnitude, natural and artificial D.C. electric fields associated with shore crabs and dipole electrodes respectively. We found a highly significant preference for the stronger D.C. electric field and a less pronounced, but still significant, preference for the A.C. electric field rather than the D.C. electric field. No preference was demonstrated between the artificial and natural D.C. electric fields. The findings are discussed in relation to the animal’s diet and ecology and with regard to anthropogenic sources of electric fields within their habitat.     

Modeling traffic air pollution in street canyons in New York City for intra-urban exposure assessment in the US Multi-Ethnic Study of atherosclerosis and air pollution

We evaluated the Danish AirGIS air quality and exposure model system using air quality measurement data from New York City in the Multi-Ethnic Study of Atherosclerosis and Air Pollution (MESA Air). Measurements were used from three US EPA Air Quality System (AQS) monitoring stations and a comprehensive MESA Air measurement campaign including about 150 different locations and about 650 samples of about 2 week measurements of NO_x, NO₂ and PM_2.5. AirGIS is a deterministic exposure model system based on the dispersion models Operational Street Pollution Model (OSPM) and the Urban Background Model (UBM). The UBM model reproduced the annual levels within 1–26% depending on station and pollutant at the three urban background EPA monitor stations, and generally reproduced well the seasonal and diurnal variation. The full model with OSPM and UBM reproduced the MESA Air measurements with a correlation coefficient of r² = 0.51 for NO_x, r² = 0.28 for NO₂ and r² = 0.73 for PM_2.5.     

Some Lessons Learned from Public Health on the Process of Adaptation

Lessons learned from more than 150 years of public health research and intervention can provide insights to guide public health professionals and institutions as they design and implement specific strategies, policies, and measures to increase resilience to climate variability and change. This paper identifies both some modifications to public health systems that may enhance adaptive capacity, and lessons drawn from the history of managing environmental and other threats in the public health sector that may have relevance for other sectors as they design approaches to increase their adaptive capacity to more effectively cope with climate variability and change. The views expressed are the author’s own and do not represent official US EPA policy.     

Biological and environmental influences on the trophic ecology of leatherback turtles in the northwest Atlantic Ocean

Understanding the causes and consequences of variability in trophic status is important for interpreting population dynamics and for identifying important habitats for protected species like marine turtles. In the northwest Atlantic Ocean, many leatherback turtles (Dermochelys coriacea) from distinct breeding stocks throughout the Wider Caribbean region migrate to Canadian waters seasonally to feed, but their trophic status during the migratory and breeding cycle and its implications have not yet been described. In this study, we used stable carbon and nitrogen isotope analyses of bulk skin to characterize the trophic status of leatherbacks in Atlantic Canadian waters by identifying trophic patterns among turtles and the factors influencing those patterns. δ¹⁵N values of adult males and females were significantly higher than those of turtles of unknown gender (i.e., presumed to be subadults), and δ¹⁵N increased significantly with body size. We found no significant differences among average stable isotope values of turtles according to breeding stock origin. Significant inter-annual variation in δ¹⁵N among cohorts probably reflects broad-scale oceanographic variability that drives fluctuations in stable isotope values of nutrient sources transferred through several trophic positions to leatherbacks, variation in baseline isotope values among different overwintering habitats used by leatherbacks, or a combination of both. Our results demonstrate that understanding effects of demographic and physiological factors, as well as oceanographic conditions, on trophic status is key to explaining observed patterns in population dynamics and for identifying important habitats for widely distributed, long-lived species like leatherbacks.