Body Height Estimation from Femur Measurements in Postmortem Computed Tomography

After attending this presentation, attendees will: (1) understand how body height from computed tomography data can be estimated; and, (2) gain knowledge about the accuracy of estimated body height and limitations. The presentation will impact the forensic science community by providing knowledge and competence which will enable attendees to develop formulas for single bones to reconstruct body height using postmortem Computer Tomography (p-CT) data. The estimation of Body Height (BH) is an important component of the identification of corpses and skeletal remains. Stature can be estimated with relative accuracy via the measurement of long bones, such as the femora. Compared to time-consuming maceration procedures, p-CT allows fast and simple measurements of bones. This study undertook four objectives concerning the accuracy of BH estimation via p-CT: (1) accuracy between measurements on native bone and p-CT imaged bone (F1 according to Martin 1914); (2) intra-observer p-CT measurement precision; (3) accuracy between formula-based estimation of the BH and conventional body length measurement during autopsy; and, (4) accuracy of different estimation formulas available.1 In the first step, the accuracy of measurements in the CT compared to those obtained using an osteometric board was evaluated on the basis of eight defleshed femora. Then the femora of 83 female and 144 male corpses of a Swiss population for which p-CTs had been performed, were measured at the Institute of Forensic Medicine in Bern. After two months, 20 individuals were measured again in order to assess the intraobserver error. The mean age of the men was 53±17 years and that of the women was 61±20 years. Additionally, the body length of the corpses was measured conventionally. The mean body length was 176.6±7.2cm for men and 163.6±7.8cm for women. The images that were obtained using a six-slice CT were reconstructed with a slice thickness of 1.25mm. Analysis and measurements of CT images were performed on a multipurpose workstation. As a forensic standard procedure, stature was estimated by means of the regression equations by Penning & Riepert developed on a Southern German population and for comparison, also those referenced by Trotter & Gleser “American White.”2,3 All statistical tests were performed with a statistical software. No significant differences were found between the CT and osteometric board measurements. The double p-CT measurement of 20 individuals resulted in an absolute intra-observer difference of 0.4±0.3mm. For both sexes, the correlation between the body length and the estimated BH using the F1 measurements was highly significant. The correlation coefficient was slightly higher for women. The differences in accuracy of the different formulas were small. While the errors of BH estimation were generally ±4.5–5.0cm, the consideration of age led to an increase in accuracy of a few millimetres to about 1cm. BH estimations according to Penning & Riepert and Trotter & Gleser were slightly more accurate when age-at-death was taken into account.2,3 That way, stature estimations in the group of individuals older than 60 years were improved by about 2.4cm and 3.1cm.2,3 The error of estimation is therefore about a third of the common ±4.7cm error range. Femur measurements in p-CT allow very accurate BH estimations. Estimations according to Penning led to good results that (barely) come closer to the true value than the frequently used formulas by Trotter & Gleser “American White.”2,3 Therefore, the formulas by Penning & Riepert are also validated for this substantial recent Swiss population.

Body height estimation from post-mortem CT femoral F1 measurements in a contemporary Swiss population.

PURPOSE The present study aimed at the comparison of body height estimations from cadaver length with body height estimations according to Trotter and Gleser (1952) and Penning and Riepert (2003) on the basis of femoral F1 section measurements in post-mortem computed tomography (PMCT) images. METHODS In a post-mortem study in a contemporary Swiss population (226 corpses: 143 males (mean age: 53±17years) and 83 females (mean age: 61±20years)) femoral F1 measurements (403 femora: 199 right and 204 left; 177 pairs) were conducted in PMCT images and F1 was used for body height estimation using the equations after Trotter and Gleser (1952, "American Whites"), and Penning and Riepert (2003). RESULTS The mean observed cadaver length was 176.6cm in males and 163.6cm in females. Mean measured femoral length F1 was 47.5cm (males) and 44.1cm (females) respectively. Comparison of body height estimated from PMCT F1 measurements with body height calculated from cadaver length showed a close congruence (mean difference less than 0.95cm in males and less than 1.99cm in females) for equations both applied after Penning and Riepert and Trotter and Gleser. CONCLUSIONS Femoral F1 measurements in PMCT images are very accurate, reproducible and feasible for body height estimation of a contemporary Swiss population when using the equations after Penning and Riepert (2003) or Trotter and Gleser (1952).

A categorical data analysis approach to estimation of demographic characteristics for Texas county populations

The need for timely population data for health planning and Indicators of need has Increased the demand for population estimates. The data required to produce estimates is difficult to obtain and the process is time consuming. Estimation methods that require less effort and fewer data are needed. The structure preserving estimator (SPREE) is a promising technique not previously used to estimate county population characteristics. This study first uses traditional regression estimation techniques to produce estimates of county population totals. Then the structure preserving estimator, using the results produced in the first phase as constraints, is evaluated.^ Regression methods are among the most frequently used demographic methods for estimating populations. These methods use symptomatic indicators to predict population change. This research evaluates three regression methods to determine which will produce the best estimates based on the 1970 to 1980 indicators of population change. Strategies for stratifying data to improve the ability of the methods to predict change were tested. Difference-correlation using PMSA strata produced the equation which fit the data the best. Regression diagnostics were used to evaluate the residuals.^ The second phase of this study is to evaluate use of the structure preserving estimator in making estimates of population characteristics. The SPREE estimation approach uses existing data (the association structure) to establish the relationship between the variable of interest and the associated variable(s) at the county level. Marginals at the state level (the allocation structure) supply the current relationship between the variables. The full allocation structure model uses current estimates of county population totals to limit the magnitude of county estimates. The limited full allocation structure model has no constraints on county size. The 1970 county census age - gender population provides the association structure, the allocation structure is the 1980 state age - gender distribution.^ The full allocation model produces good estimates of the 1980 county age - gender populations. An unanticipated finding of this research is that the limited full allocation model produces estimates of county population totals that are superior to those produced by the regression methods. The full allocation model is used to produce estimates of 1986 county population characteristics. ^

AN ESTIMATION OF IN UTERO DRUG EXPOSURE BY DETERMINANTS OF FREQUENCY AND DISTRIBUTION OF MATERNAL DRUG CONSUMPTION

Data derived from 1,194 gravidas presenting at the observation unit of a city/county hospital between October 11, 1979 through December 7, 1979 were evaluated with respect to the proportion ingesting drugs during pregnancy. The mean age of the mother at the time of the interview was 22.0 years; 43.0 percent were Black; 34.0 percent Latin-American, 21.0 percent White and 2.0 percent other; mean gravida was 2.5 pregnancies; mean parity was 1.0; and mean number of previous abortions was 0.34. Completed interview data was available for 1,119 gravida, corresponding urinalyses for 997 subjects. Ninety and one-tenth percent (90.1 percent) of the subjects reported ingestion of one or more drug preparation(s) (prescription, OTC, or substances used for recreational purposes) during pregnancy with a range of 0 to 11 substances and a mean of 2.7. Dietary supplements (vitamins and minerals) were most frequently reported followed by non-narcotic analgesics. Seventy-six and one tenth percent (76.1 percent) of the population reported consumption of prescription medication, 42.5 percent reported consumption of over-the-counter medications, 45.7 percent reported consumption of a substance for recreational purposes and 4.3 percent reported illicit consumption of a substance. For selected substances, no measurable difference was found between obtaining the information from the interview method or from a urinalysis assay. ^

AN EMPIRICAL ESTIMATION OF MEDICAL AND NON-MEDICAL INPUTS IN THE PRODUCTION OF CARDIO-RESPIRATORY PHYSICAL HEALTH STATUS

The purpose of this study was to examine, in the context of an economic model of health production, the relationship between inputs (health influencing activities) and fitness.^ Primary data were collected from 204 employees of a large insurance company at the time of their enrollment in an industrially-based health promotion program. The inputs of production included medical care use, exercise, smoking, drinking, eating, coronary disease history, and obesity. The variables of age, gender and education known to affect the production process were also examined. Two estimates of fitness were used; self-report and a physiologic estimate based on exercise treadmill performance. Ordinary least squares and two-stage least squares regression analyses were used to estimate the fitness production functions.^ In the production of self-reported fitness status the coefficients for the exercise, smoking, eating, and drinking production inputs, and the control variable of gender were statistically significant and possessed theoretically correct signs. In the production of physiologic fitness exercise, smoking and gender were statistically significant. Exercise and gender were theoretically consistent while smoking was not. Results are compared with previous analyses of health production. ^

Magnesium/Calcium ratios and sea surface temperature estimation for sediments of the Timor Sea

We use a multiproxy approach to monitor changes in the vertical profile of the Indonesian Throughflow as well as monsoonal wind and precipitation patterns in the Timor Sea on glacial-interglacial, precessional, and suborbital timescales. We focus on an interval of extreme climate change and sea level variation: marine isotope (MIS) 6 to MIS 5e. Paleoproductivity fluctuations in the Timor Sea follow a precessional beat related to the intensity of the Australian (NW) monsoon. Paired Mg/Ca and d18O measurements of surface- and thermocline-dwelling planktonic foraminifers (G. ruber and P. obliquiloculata) indicate an increase of >4°C in both surface and thermocline water temperatures during Termination II. Tropical sea surface temperature changed synchronously with ice volume (benthic d18O) during deglaciation, implying a direct coupling of high- and low-latitude climate via atmospheric and/or upper ocean circulation. Substantial cooling and freshening of thermocline waters occurred toward the end of Termination II and during MIS 5e, indicating a change in the vertical profile of the Indonesian Throughflow from surface- to thermocline-dominated flow.

Western Arabian Sea SST estimation

Millennial-scale records of planktonic foraminiferal Mg/Ca, bulk sediment UK37', and planktonic foraminiferal d18O are presented across the last two deglaciations in sediment core NIOP929 from the Arabian Sea. Mg/Ca-derived temperature variability during the penultimate and last deglacial periods falls within the range of modern day Arabian Sea temperatures, which are influenced by monsoon-driven upwelling. The UK37'-derived temperatures in MIS 5e are similar to modern intermonsoon values and are on average 3.5°C higher than the Mg/Ca temperatures in the same period. MIS 5e UK37' and Mg/Ca temperatures are 1.5°C warmer than during the Holocene, while the UK37'-Mg/Ca temperature difference was about twice as large during MIS 5e. This is surprising as, nowadays, both proxy carriers have a very similar seasonal and depth distribution. Partial explanations for the MIS 5e UK37'-Mg/Ca temperature offset include carbonate dissolution, the change in dominant alkenone-producing species, and possibly lateral advection of alkenone-bearing material and a change in seasonal or depth distribution of proxy carriers. Our findings suggest that (1) Mg/Ca of G. ruber documents seawater temperature in the same way during both studied deglaciations as in the present, with respect to, e.g., season and depth, and (2) UK37'-based temperatures from MIS 5 (or older) represent neither upwelling SST nor annual average SST (as it does in the present and the Holocene) but a higher temperature, despite alkenone production mainly occurring in the upwelling season. Further we report that at the onset of the deglacial warming, the Mg/Ca record leads the UK37' record by 4 ka, of which a maximum of 2 ka may be explained by postdepositional processes. Deglacial warming in both temperature records leads the deglacial decrease in the d18O profile, and Mg/Ca-based temperature returns to lower values before d18O has reached minimum interglacial values. This indicates a substantial lead in Arabian Sea warming relative to global ice melting.

Cretaceous stable isotopic record and sea surface temperture estimation for ODP Site 207-1259

The Turonian (93.5 to 89.3 million years ago) was one of the warmest periods of the Phanerozoic eon, with tropical sea surface temperatures over 35°C. High-amplitude sea-level changes and positive d18O excursions in marine limestones suggest that glaciation events may have punctuated this episode of extreme warmth. New d18O data from the tropical Atlantic show synchronous shifts ~91.2 million years ago for both the surface and deep ocean that are consistent with an approximately 200,000-year period of glaciation, with ice sheets of about half the size of the modern Antarctic ice cap. Even the prevailing supergreenhouse climate was not a barrier to the formation of large ice sheets, calling into question the common assumption that the poles were always ice-free during past periods of intense global warming.

Stable isotope ratios of foraminifera and sea surface temperature estimation of sediment core W2804A-14

Carbon isotopically based estimates of CO2 levels have been generated from a record of the photosynthetic fractionation of 13C (epsilon p) in a central equatorial Pacific sediment core that spans the last ~255 ka. Contents of 13C in phytoplanktonic biomass were determined by analysis of C37 alkadienones. These compounds are exclusive products of Prymnesiophyte algae which at present grow most abundantly at depths of 70-90 m in the central equatorial Pacific. A record of the isotopic compostion of dissolved CO2 was constructed from isotopic analyses of the planktonic foraminifera Neogloboquadrina dutertrei, which calcifies at 70-90 m in the same region. Values of epsilon p, derived by comparison of the organic and inorganic delta values, were transformed to yield concentrations of dissolved CO2 (c e) based on a new, site-specific calibration of the relationship between epsilon p and c e. The calibration was based on reassessment of existing epsilon p versus c e data, which support a physiologically based model in which epsilon p is inversely related to c e. Values of PCO2, the partial pressure of CO2 that would be in equilibrium with the estimated concentrations of dissolved CO2, were calculated using Henry's law and the temperature determined from the alkenone-unsaturation index UK 37. Uncertainties in these values arise mainly from uncertainties about the appropriateness (particularly over time) of the site-specific relationship between epsilon p and 1/c e. These are discussed in detail and it is concluded that the observed record of epsilon p most probably reflects significant variations in Delta pCO2, the ocean-atmosphere disequilibrium, which appears to have ranged from ~110 µatm during glacial intervals (ocean > atmosphere) to ~60 µatm during interglacials. Fluxes of CO2 to the atmosphere would thus have been significantly larger during glacial intervals. If this were characteristic of large areas of the equatorial Pacific, then greater glacial sinks for the equatorially evaded CO2 must have existed elsewhere. Statistical analysis of air-sea pCO2 differences and other parameters revealed significant (p < 0.01) inverse correlations of Delta pCO2 with sea surface temperature and with the mass accumulation rate of opal. The former suggests response to the strength of upwelling, the latter may indicate either drawdown of CO2 by siliceous phytoplankton or variation of [CO2]/[Si(OH)4] ratios in upwelling waters.

Salinity estimation from Gulf of Guinea, sediment core MD03-2707

The isotopic composition of surface seawater is widely used to infer past changes in sea surface salinity using paired foraminiferal Mg/Ca and d18O from marine sediments. At low latitudes, paleosalinity reconstructions using this method have largely been used to document changes in the hydrological cycle. This method usually assumes that the modern seawater d18O (d18Osw)/salinity relationship remained constant through time. Modelling studies have shown that such assumptions may not be valid because large-scale atmospheric circulation patterns linked to global climate changes can alter the seawater d18Osw/salinity relationship locally. Such processes have not been evidenced by paleo-data so far because there is presently no way to reconstruct past changes in the seawater d18Osw/salinity relationship. We have addressed this issue by applying a multi-proxy salinity reconstruction from a marine sediment core collected in the Gulf of Guinea. We measured hydrogen isotopes in C37:2 alkenones (dDa) to estimate changes in seawater dD. We find a smooth, long-term increase of ~10 per mil in dDa between 10 and 3 kyr BP, followed by a rapid decrease of ~10 per mil in dDa between 3 kyr BP and core top to values slightly lighter than during the early Holocene. Those features are inconsistent with published salinity estimations based on d18Osw and foraminiferal Ba/Ca, as well as nearby continental rainfall history derived from pollen analysis. We combined dDa and d18Osw values to reconstruct a Holocene record of salinity and compared it to a Ba/Ca-derived salinity record from the same sedimentary sequence. This combined method provides salinity trends that are in better agreement with both the Ba/Ca-derived salinity and the regional precipitation changes as inferred from pollen records. Our results illustrate that changes in atmospheric circulation can trigger changes in precipitation isotopes in a counter-intuitive manner that ultimately impacts surface salinity estimates based on seawater isotopic values. Our data suggest that the trends in Holocene rainfall isotopic values at low latitudes may not uniquely result from changes in local precipitation associated with the amount effect.

Sea surface temperature estimation for the last 70 ka from IODP Site 306-U1313

During the six Heinrich Events of the last 70 ka episodic calving from the circum-Atlantic ice sheets released large numbers of icebergs into the North Atlantic. These icebergs and associated melt-water flux are hypothesized to have led to a shutdown of Atlantic Meridional Overturning Circulation (AMOC) and severe cooling in large parts of the Northern Hemisphere. However, due to the limited availability of high-resolution records the magnitude sea surface temperature (SST) changes related to the impact of Heinrich Events on the mid-latitude North Atlantic is poorly constrained. Here we present a record of UK37'-based SSTs derived from sediments of Integrated Ocean Drilling Project (IODP) Site U1313, located at the southern end of the ice-rafted debris (IRD)-belt in the mid-latitude North Atlantic (41°N). We demonstrate that all six Heinrich Events are associated with a rapid warming of surface waters by 2 to 4°C in a few thousand years. The presence of IRD leaves no doubt about the simultaneous timing and correlation between rapid surface water warming and Heinrich Events. We argue that this warming in the mid-latitude North Atlantic is related to a northward expansion of the subtropical gyre during Heinrich Events. As a wide-range of studies demonstrated that in the central IRD-belt Heinrich Events are associated with low SSTs, these results thus identify an anti-phased (seesaw) pattern in SSTs during Heinrich Events between the mid-latitude (warm) and northern North Atlantic (cold). This highlights the complex response of surface water characteristics in the North Atlantic to Heinrich Events that is poorly reproduced by fresh water hosing experiments and challenges the widely accepted view that within the IRD-belt of the North Atlantic Heinrich Events coincide with periods of low SSTs.

Sea surface temperature estimation for MIS 16 of IODP Site 306-U1313

Hudson Strait (HS) Heinrich Events, ice-rafting events in the North Atlantic originating from the Laurentide ice sheet (LIS), are among the most dramatic examples of millennial-scale climate variability and have a large influence on global climate. However, it is debated as to whether the occurrence of HS Heinrich Events in the (eastern) North Atlantic in the geological record depends on greater ice discharge, or simply from the longer survival of icebergs in cold waters. Using sediments from Integrated Ocean Drilling Program (IODP) Site U1313 in the North Atlantic spanning the period between 960 and 320 ka, we show that sea surface temperatures (SSTs) did not control the first occurrence of HS Heinrich(-like) Events in the sedimentary record. Using mineralogy and organic geochemistry to determine the characteristics of ice-rafting debris (IRD), we detect the first HS Heinrich(-like) Event in our record around 643 ka (Marine Isotope Stage (MIS) 16), which is similar as previously reported for Site U1308. However, the accompanying high-resolution alkenone-based SST record demonstrates that the first HS Heinrich(-like) Event did not coincide with low SSTs. Thus, the HS Heinrich(-like) Events do indicate enhanced ice discharge from the LIS at the end of the Mid-Pleistocene Transition, not simply the survivability of icebergs due to cold conditions in the North Atlantic.

Sea surface temperature estimation of the South Tasman Rise

The Subtropical Front (STF) marking the northern boundary of the Southern Ocean has a steep gradient in sea surface temperature (SST) of approximately 4°C over 0.5° of latitude. Presently, in the region south of Tasmania, the STF lies nominally at 47°S in the summer and 45°S in the winter. We present here SST reconstructions in a latitudinal transect of cores across the South Tasman Rise, southeast of Australia, during the late Quaternary. SST reconstructions are based on two paleotemperature proxies, alkenones and faunal assemblages, which are used to assess past changes in SST in spring and summer. The north-south alignment in core locations allows reconstruction of movement of the STF over the last 100 ka. Surface water temperatures during the last glaciation in this region were ~4°C colder than today. Additional temperature changes greater in magnitude than 4°C seen in individual cores can be attributed to changes in the water mass overlying the core site caused by the movement of the front across that location. During the penultimate interglacial, SST was ~2°C warmer and the STF was largely positioned south of 47°S. Movement of the STF to the north occurred during cool climate periods such as the last marine isotope stages 3 and 4. In the last glaciation, the front was at its farthest north position, becoming pinned against the Tasmanian landmass. It moved south by 4° latitude to 47°S in summer during the deglaciation but remained north of 45°S in spring throughout the early deglaciation. After 11 ka B.P. inferred invigoration of the East Australia Current appears to have pushed the STF seasonally south of the East Tasman Plateau, until after 6 ka B.P. when it achieved its present configuration.