Article Text


Evaluation of decompression safety in an occupational diving group using self reported diving exposure and health status
  1. D J Doolette1,
  2. D F Gorman2
  1. 1Anaesthesia & Intensive Care, The University of Adelaide, Australia
  2. 2Department of Medicine, University of Auckland, New Zealand; Anaesthesia & Intensive Care, The University of Adelaide, Australia
  1. Correspondence to:
 Dr D Doolette, Anaesthesia & Intensive Care, The University of Adelaide, Adelaide, Australia 5005; 


Background: Many occupational diving groups have substantially different diving patterns to those for which decompression schedules are validated.

Aims: To evaluate tuna farm occupational diving practice against existing decompression models and describe a method for collecting and modelling self reported field decompression data.

Methods: Machine readable objective depth/time profiles were obtained from depth/time recorders worn by tuna farm occupational divers. Divers’ health status was measured at the end of each working day using a self administered health survey that produces an interval diver health score (DHS) with possible values ranging from 0 to 30. Depth/time profiles were analysed according to existing decompression models. The contribution of diving exposure and between diver variability to DHS was evaluated using linear regression.

Results: The mean risk of decompression sickness was calculated as 0.005 (SD 0.003, n = 383). The mean DHS following diving was 3 (SD 2, n = 383) and following non-diving activities was 1 (SD 1, n = 41). After accounting for between diver variability in intercept, DHS was found to increase one unit for every 1% increase in the risk of decompression sickness.

Conclusions: A method has been established for the collection and analysis of self reported objective decompression data from occupational diving groups that can potentially be used as the basis for development of purpose designed occupational diving decompression schedules.

  • decompression sickness
  • environmental health
  • regression analysis
  • AIC, Akaike Information Criteria
  • DCIEM, Defence and Civil Institute of Environmental Medicine
  • DCS, decompression sickness
  • DHS, diver health score
  • DUR, dive duration
  • MSW, metres sea water
  • NUM, number of divers
  • pDCS, risk of DCS

Statistics from

Decompression sickness (DCS) is a significant health risk for occupational divers. The risk of DCS (pDCS) is minimised by use of a schedule of diving depths and durations and consequent decompression rates. Most decompression schedules are based on rudimentary models of gas kinetics, then tested under the intended operational conditions. pDCS can vary considerably with changes in diving conditions,1 and decompression schedules should be selected on proven reliability for similar exposures in a similar population. Instead, most decompression schedules are used generically such that uncertainty exists in the prediction of DCS. Most generically available decompression schedules originate from the military where it has been possible to conduct extensive human testing. This may be a problem for emerging occupational diving groups without specifically developed decompression schedules.

The tuna farming industry in Australia involves diving practices that are substantially different from those of conventional occupational and military divers.2 The Australian Standard for Occupational Diving3 requires the use of the Defence and Civil Institute of Environmental Medicine (DCIEM) standard air diving schedules4 for all occupational diving. These schedules are well tested and have a good operational record, but have not been evaluated in the unique context of tuna farm diving practices. Tuna farm diving practice involves extensive repetitive diving and multi-day diving (consecutive daily diving exposures); multi-day diving is avoided in military decompression schedule validation data. Consequently, the pDCS for tuna farm diving is unknown.

Widespread use of depth/time recorders by occupational divers makes accurate field data available from operational dives. DCS may not be reliably reported from dives conducted without medical supervision. However, the development of a scoring system for health status self assessment after diving provides a reliable measure of diving outcome.5 With careful audit, these data can be used to evaluate decompression practice against existing decompression models. Traditional deterministic decompression procedures, such as the DCIEM standard air diving, define acceptable diving schedules. More recent, probabilistic decompression models can assign pDCS to any depth/time profile. These latter models are produced using parameter estimation procedures to select decompression models of best fit to observed decompression data.6

The aims of this study are to describe tuna farm occupational diving practice and evaluate this practice against existing deterministic and probabilistic decompression models. A further aim is to describe and evaluate a method for collecting self reported field decompression data with the hypothesis that self reported health status is an indicator of decompression stress.


Data collection

We used a prospective decompression data collection approach (approved by the University of Adelaide Human Ethics Committee and conducted in accordance with the “National statement on ethical conduct in research involving humans, 1999”) that was based on machine readable objective depth/time profiles data and a novel self administered diver health survey. Machine readable depth/time profiles were obtained from diver carried depth/time recorders. These devices record gauge pressure (as depth of water) at fixed, device dependent intervals typically of 5 or 20 seconds.

Development and validation of the diver health survey has been described in detail elsewhere.5 It is an inventory of nine standardised items and responses covering five symptoms of decompression illness (paraesthesia, rash, balance, vitality, and pain), five health status indicators (vitality, pain, physical functioning, role limitation, and health perception), and time of onset of symptoms, plus one free response, each item scored from 0 to 3. The resulting summed diver health score (DHS) can be analysed as interval data and is correlated with clinical decompression illness.5 The validated format of the diver health survey and scoring instructions are available from the corresponding author. The DHS was used as the outcome measure without any attempt to categorise dives as resulting in DCS or not.

The symptoms of DCS and of other decompression illnesses that arise from lung barotrauma are often indistinguishable.7 This is a limitation for validation of decompression models, since these can only predict DCS. However, lung barotrauma would not be expected as a common feature of the experienced occupational diver group in the present study.

Decompression data collection was organised through the tuna farm diving supervisors. Generally, for 1–2 month periods, dive team members completed a diver health survey at the end of any workday that included a diving exposure. Health surveys were returned in confidence by reply paid mail to our laboratory. Depth/time profiles were uploaded from recorders by either the diving supervisor or the farm secretarial staff, using the proprietary software for these devices. Depth/time profiles were copied to floppy disk and returned by reply paid mail. For baseline data, some divers completed health surveys following non-diving work activities.

Decompression data were managed using purpose designed, partially automated database and analysis applications programmed in our laboratory (Access 2000, Visual BASIC, Microsoft Corp., Redmond, WA; and Statistica BASIC, Statsoft, Inc., Tulsa, OK). Depth/time profile audit consisted of cross checking between daily log information (diver name, date, and number of dives) supplied separately by individual divers (health survey) and by company (depth/time profiles). Prior to analysis, to reduce computing time, raw depth/time profiles were reduced to a smaller number of depth/time nodes: distance between adjacent nodes was increased by sequential exclusion of intervening data points until reaching a threshold residual sum of squares between the excluded points and a line connecting the nodes. Raw depth/time profile and extracted nodes were visually inspected to exclude data errors. Repetitive dives (surface interval less than 18 hours) and surface intervals were combined into a single daily depth/time profile.

From a larger database of health surveys and depth/time profiles, a diving data subset was identified comprising 383 verified decompression datum pairs, where both the full day’s depth/time profiles and DHS were available (383 health surveys covering 527 dives conducted by 23 divers). Only divers contributing more than one DHS were included. Air was breathed throughout all dives. Each verified decompression datum was considered independent. In addition, nine of these divers contributed health surveys (n = 41) at the end of workdays during which they did not dive (non-diving data subset).

Evaluation of decompression practice

Each daily depth/time profile was checked for compliance with the DCIEM standard air diving schedules. The US Navy linear exponential kinetics probabilistic decompression model (USN93), which is a useful predictor of pDCS in military air diving,8 was used for reporting of pDCS for each daily depth/time profile. An implementation of the USN93 model similar to that previously described,9 using published parameters,8 was written in the BASIC programming language (Visual BASIC, Microsoft Corp., Redmond, WA.) The pDCS was tracked over the daily depth/time profile and subsequent 18 hours.

Modelling of DHS

The DHS has been previously shown to correlate with clinical diagnosis for decompression illness. Here the DHS was used to measure decompression related health status among divers during normal diving operations; therefore, the contribution of diving exposure and interdiver variability to DHS in the absence of diagnosed decompression illness was evaluated by linear regression using the combined diving (n = 383) and non-diving (n = 41) datasets. To accommodate possible interdiver variability we used a linear mixed effect modelling approach. The full model investigated was of the form:


which comprised the dependent variable DHS and fixed explanatory variables, pDCS, maximum daily depth in metres sea water (MSW), total daily dive duration in minutes (DUR), and the number of divers per day (NUM). The 23 divers were considered a random sample from a population where the intercept (β0) of the regression on the explanatory variables depends on the diver. Subscript i denotes divers, subscript j denotes days, and e denotes error.

Parameters of the regression models were estimated by maximising the likelihood. The likelihood is the joint probability density function of the observed values of the dependent variable given the respective regression model. To find the most parsimonious model, explanatory variables were manually removed from the full model and the resulting reduced models also fit to the data. Significant difference between nested models was evaluated by likelihood ratio test, 2(LLp − LLk) ≅ χ2p−k, where LL is the maximised log likelihood of the model and p and k are the number of parameters in the full and reduced model respectively (p > k). Models were selected according to minimum Akaike Information Criteria (AIC) which evaluate model fit to the data (maximised log likelihood, LL) and penalises for model complexity (number of estimated parameters, p), AIC = −2LL + 2p. For the selected models, data were examined for influential values; these were considered outlier values with high leverage. Outliers were data with standardised residual more than 2 standard deviations from the mean. Leverage was taken as the diagonal of the hat matrix and values more than twice the mean were considered high.

All statistical calculations were performed using R software base package (version 1.4.1; The R Development Core Team, 2002) and the non-linear mixed effect package (version 3.1–23; Pinheiro J, Bates D, DebRoy S, Sarkar D, 2001).


General diving practice

Table 1 summarises the diving data subset. Repetitive dives were those separated by surface intervals greater than 10 minutes and less than 18 hours; the mean surface interval was 97 minutes (SD 56, n = 144). There were 253 single dives, 116 double repetitive dive series, and 14 triple repetitive dive series. Repetitive diving exposures were generally deeper and longer than single dive exposure. Figure 1 shows typical depth/time profiles. Single dives were often long duration, multi-level dives for inspection and maintenance of the tow cages used to move tuna from the offshore fishing grounds or of the near-shore feeding pens. Repetitive dives were often short duration dives to remove dead tuna from the bottom of the feeding pens. Consecutive daily diving exposure (multi-day diving) was common in tuna farm diving and occurred in 2 day (n = 110), 3 day (n = 75), 4 day (n = 48), and 5 day (n = 10) series. Table 2 gives the contribution of the individual divers to the diving data subset.

Table 1

Diving data subset summary

Table 2

Diving data subset

Figure 1

Typical depth/time profiles. (A) Single net inspection dive. (B) Repetitive “mortality” dives.

Evaluation of decompression practice

The mean interval between diving and health self assessment was 3 (SD 3, n = 382) hours, excluding the one diver health survey completed after a lengthy retrieval (144 hours) from a remote dive site for recompression therapy. This was the only diver who sought medical treatment for decompression illness (DHS = 14). The mean DHS following diving was 3 (SD 2, n = 346) and was not different between single and repetitive dives (table 1). Figure 2 shows the distribution of DHS following all diving activities. Nine of the divers completed health surveys following non-diving activities (table 3); the mean non-diving DHS was 1 (SD 1, n = 41) but there were significant differences between individual divers’ non-diving DHS (one way ANOVA, df = 8, F = 3.63, p = 0.004).

Table 3

Non-diving data subset

Figure 2

Frequencies of DHS measured following 383 diving days (bars) and DHS predicted from pDCS by model 4 (circles).

Tuna farm divers are required to adhere to DCIEM standard air decompression schedules; only one daily exposure violated these schedules and the health outcome was unremarkable (DHS = 2). A pDCS ⩾ 0.0156 is calculated from the USN93 model for the relevant no-decompression stop DCIEM standard air schedules. The mean pDCS for daily tuna farm daily depth/time profiles according to the USN93 model was 0.005 (SD 0.003) (see table 1). The highest risk exposures in the present data were typically short duration repetitive dives to depths in excess of 18 msw (280 kPa) and with rapid decompression to the surface (fig 1B). The instantaneous risk evaluated in the USN93 model always declined to zero before the end of the 18 hour post-dive integration period and typically before 4 hours. This supports treatment of daily depth/time profiles as independent in the subsequent modelling procedures.

Modelling of DHS

The significant differences in non-diving DHS between divers (table 3) suggest there may be considerable inter-diver variability in all DHS. This was supported by factorial ordinary least squares regression with DHS as the dependent variable, pDCS as the independent variable, and diver as the between group factor; this produced highly variable parameters estimates between divers.

This inter-diver variability is accommodated in models compared in table 4; the variance of model intercept between divers (not shown) was 1.4 in each model. In the full model (1), explanatory variable MSW and NUM had non-significant parameter estimates and these were removed to produce a simpler model with similar fit to the data (model 2). Removal from model 2 of either DUR (model 3) or pDCS (model 4) resulted in significantly worse fit to the data. The null model (model 5) that comprised only the random intercept term produced significantly worse fit to the data than all other models.

Table 4

Model comparisons

Although inclusion of explanatory variable DUR in model 2 significantly improved fit to the data compared to model 4, the parameter estimate of β2 indicated only a unit change in DHS unit for every 167 minutes diving duration, unsubstantial in the context of mean dive duration of 32 minutes. The estimates of the common parameters were similar and the difference in fit was small, so the simpler model 4 was chosen for further investigation. In model 4 the estimate of parameter β1 indicated approximately a unit change in DHS for every 1% change in pDCS. The intercept values ranged from 0.1 to 4.7 between divers. The four DHS values >8 were identified as outliers, but these did not have high leverage so were retained in the data. Figure 2 shows model 4 estimates of DHS: for the diving data subset the fitted values of DHS were rounded to the nearest integer and the frequency of each integer value (circles) superimposed on the frequency histogram of observed data (bars).

It was investigated whether there was a random effect of divers on both intercept and slope by allowing parameter β1 to vary between divers, (β1 + ηi)pDCSij. The improved fit to the data (AIC = 1596, LL = −792) was not significantly different from model 4 (likelihood ratio 5.733, df = 5,4, p = 0.0569).

In this form of modelling the data are censored at the time the diver health survey is completed as symptoms of DCS may subsequently arise. However, this censoring is probably not severe as symptom onset occurs by three hours (mean interval between diving and health self assessment) in approximately 80% of cases of DCS.10 The contribution of this interval was examined by adding the fixed explanatory variable and parameter β5SINCEij to model 4. The resulting model could only be fit to the diving data subset as there was no SINCE variable in the non-diving data subset. In addition the highly influential outlier value of DHS = 14, SINCE = 144 was excluded. Model 4 (AIC = 1417, LL = −705) and the expanded model (AIC = 1419, LL = −704) were fit to this data subset and were not significantly different (likelihood ratio 0.921, df = 6,4, p = 0.337), and the estimate of parameter β5 was not significant (p = 0.336), indicating that the interval was not important in this dataset.


The tuna farm diving practice reported here is relatively conservative as indicated by DHS skewed towards asymptomatic values and the general adherence to DCIEM standard air diving decompression schedules. The pDCS for tuna farm diving calculated according to the USN93 model was low compared to that calculated for no-decompression stop DCIEM standard air schedules; however, these latter figures assume that the full allowed time is spent at the maximum depth, a practice that is uncommon in the field. In addition, repetitive and multi-day diving was common and both may alter risk factors for DCS. The present conservative diving practice contrasts with anecdotal reports of high risk diving practice in the early days of the tuna farming industry, which are supported by the high incidence of decompression illness from that period,11 and is indicative of a recent positive reform of diving practice.

The USN93 model appears to be a reasonable predictor of tuna farm diving practice. The mean calculated pDCS of 0.005 equates to 1.9 incidents of DCS in 383 daily diving exposures and compares favourably with the one incident of treated DCS and the three additional DHS greater than 8 that may be untreated DCS.5,11 The military decompression data used to calibrate the USN93 model includes some diving practices relevant to tuna farming, but there are differences; in particular, multi-day diving is not included in the USN93 calibration data. For such reasons, occupational diving groups would benefit from specific decompression model development.

After accounting for inter-diver variability, the DHS is a useful measure of decompression related health outcome for field data collection and is positively correlated with decompression stress arising from normal diving operations in the absence of diagnosis of decompression illness. The methods for field decompression data collection and modelling reported here are promising, and the longer term prospect is to use such data to develop decompression schedules specific for tuna farm divers. This will require a more sophisticated modelling effort, possibly in combination with data from other sources, and continued field data collection.

In this regard, the self reporting methods described here are useful both in decompression data collection and subsequent analysis. First, field data can be easily collected, although careful procedures are necessary to ensure adequate data quality. Particular attention must be paid to the method of data audit to ensure full and accurate exposure data and to minimise selective reporting. In this context, crosschecking of data between divers and employers eliminated incomplete datasets and showed no indication of selective reporting from either source. In the absence of medical supervision, confidential self assessment of health status is probably the only valid method of measuring diving outcome. Divers may not report symptoms of DCS for a variety of reasons.12 Routine, confidential self reporting of health status can alleviate peer or economic pressures against reporting DCS and eliminates the need to evaluate the significance of symptoms in the decision to report. The DHS has been shown to be a valid and reliable measure of decompression related diving outcome.5

Secondly, interval measures such as the DHS have some statistical advantages compared to categorical data, particularly for small datasets. To avoid overfitting a binary dataset, the number of independent variables and fitted parameters in a model should be constrained to a fraction (for example, 0.1) of the number of the least frequent outcome13; thus to fit complex decompression models to health outcome categorised as the presence or absence of DCS requires a dataset with many incidents of DCS. An enormous data collection programme would be required to amass such a dataset from the field, where DCS is rare. On the other hand, use of an interval data DHS for each dive, as in the present study, has advantages in the modelling process because the number of fitted parameters is constrained by total sample size. Using these methods, a useful dataset can be acquired relatively easily. The present methods cannot replace well controlled, medically supervised dives, which produce the primary data for decompression model development, but can complement this latter approach by making available field decompression data describing a variety of different diving practices.

Main messages

  • Self assessment of health status can be used as an indicator of decompression outcome and correlates with decompression stress following normal diving operations.

  • Recent tuna farm diving practice is conservative; nevertheless decompression sickness still occurs in spite of conservative diving practice.

Policy implications

  • A relatively simple method exists for evaluating decompression practice in emerging diving industries.


This work was supported by a grant from WorkCover Corporation. We would like to thank the Tuna Boat Owners Association of Australia and the diving supervisors, divers, and secretarial staff from the tuna farming industry who participated in data collection.


View Abstract

Request permissions

If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.