History of decompression research and development


Decompression in the context of diving derives from the reduction in ambient pressure experienced by the diver during the ascent at the end of a dive or hyperbaric exposure and refers to both the reduction in pressure and the process of allowing dissolved inert gases to be eliminated from the tissues during this reduction in pressure.
When a diver descends in the water column the ambient pressure rises. Breathing gas is supplied at the same pressure as the surrounding water, and some of this gas dissolves into the diver's blood and other tissues. Inert gas continues to be taken up until the gas dissolved in the diver is in a state of equilibrium with the breathing gas in the diver's lungs,, or the diver moves up in the water column and reduces the ambient pressure of the breathing gas until the inert gases dissolved in the tissues are at a higher concentration than the equilibrium state, and start diffusing out again. Dissolved inert gases such as nitrogen or helium can form bubbles in the blood and tissues of the diver if the partial pressures of the dissolved gases in the diver gets too high when compared to the ambient pressure. These bubbles, and products of injury caused by the bubbles, can cause damage to tissues known as decompression sickness or the bends. The immediate goal of controlled decompression is to avoid development of symptoms of bubble formation in the tissues of the diver, and the long-term goal is to also avoid complications due to sub-clinical decompression injury.
The symptoms of decompression sickness are known to be caused by damage resulting from the formation and growth of bubbles of inert gas within the tissues and by blockage of arterial blood supply to tissues by gas bubbles and other emboli consequential to bubble formation and tissue damage. The precise mechanisms of bubble formation and the damage they cause has been the subject of medical research for a considerable time and several hypotheses have been advanced and tested. Tables and algorithms for predicting the outcome of decompression schedules for specified hyperbaric exposures have been proposed, tested, and used, and usually found to be of some use but not entirely reliable. Decompression remains a procedure with some risk, but this has been reduced and is generally considered to be acceptable for dives within the well-tested range of commercial, military and recreational diving.
The first recorded experimental work related to decompression was conducted by Robert Boyle, who subjected experimental animals to reduced ambient pressure by use of a primitive vacuum pump. In the earliest experiments the subjects died from asphyxiation, but in later experiments, signs of what was later to become known as decompression sickness were observed. Later, when technological advances allowed the use of pressurisation of mines and caissons to exclude water ingress, miners were observed to present symptoms of what would become known as caisson disease, the bends, and decompression sickness. Once it was recognized that the symptoms were caused by gas bubbles, and that recompression could relieve the symptoms, further work showed that it was possible to avoid symptoms by slow decompression, and subsequently various theoretical models have been derived to predict low-risk decompression profiles and treatment of decompression sickness.

Timeline

Early decompression theory generally assumed that inert gas bubble formation in the tissues could be avoided during decompression, and the aim of the decompression tables and algorithms was to prevent bubble formation while minimising decompression time. Most dissolved phase models are perfusion limited, and differ mainly by the number of compartments, the range of half-times, and supersaturation tolerances assumed. These models are commonly referred to as Haldanean.

Haldane's theory and tables

was commissioned by the Royal Navy to develop a safe decompression procedure. The current method was a slow linear decompression, and Haldane was concerned that this was ineffective due to additional nitrogen buildup in the slow early stages of the ascent.
Haldane's hypothesis was that a diver could ascend immediately to a depth where the supersaturation reaches but does not exceed the critical supersaturation level, at which depth the pressure gradient for off-gassing is maximized and the decompression is most efficient. The diver would remain at this depth until saturation had reduced sufficiently for him to ascend another 10 feet, to the new depth of critical supersaturation, where the process would be repeated until it was safe for the diver to reach the surface. Haldane assumed a constant critical ratio of dissolved nitrogen pressure to ambient pressure which was invariant with depth.
A large number of decompression experiments were done using goats, which were compressed for three hours to assumed saturation, rapidly decompressed to surface pressure, and examined for symptoms of decompression sickness. Goats which had been compressed to 2.25 bar absolute or less showed no signs of DCS after rapid decompression to the surface. Goats compressed to 6 bar and rapidly decompressed to 2.6 bar also showed no signs of DCS. Haldane and his co-workers concluded that a decompression from saturation with a pressure ratio of 2 to 1 was unlikely to produce symptoms.

Haldane's model

The decompression model formulated from these findings made the following assumptions.
This model was used to compute a set of tables. The method comprises choosing a depth and time exposure, and calculation the nitrogen partial pressure in each of the tissue compartments at the end of that exposure.
Chamber tests and open water dives with two divers were made in 1906. The divers were successfully decompressed from each exposure.
The tables were adopted by the Royal Navy in 1908. The Haldane tables of 1906 are considered to be the first true set of decompression tables, and the basic concept of parallel tissue compartments with half-times and critical supersaturation limits are still in use in several later decompression models, algorithms, tables and decompression computers.

U.S. Navy decompression tables

US Navy decompression tables have gone through a lot of development over the years. They have mostly been based on parallel multi-compartment exponential models. The number of compartments has varied, and the allowable supersaturation in the various compartments during ascent has undergone major development based on experimental work and records of decompression sickness incidents.

C&R tables (1915)

The first decompression tables produced for the U.S. Navy were developed by the Bureau of Construction and Repair in 1915 and were consequently known as the C&R tables. They were derived from a Haldanean model, with oxygen decompression to depths up to 300 ft on air, and were successfully used to depths of slightly over 300 ft

Hawkins Schilling and Hansen (1930s)

Submarine escape training led US Navy personnel to believe that Haldane's allowable supersaturation ratios for fast tissues were unnecessarily conservative, as calculated values indicated that supersaturation in trainees exceeded Haldane's limits, but they did not develop DCS. A large number of experimental dives were conducted over 3 years to derive allowable supersaturation ratios for a Haldanian 5 compartment model with compartment half-times of 5, 10, 20, 40 and 70 minutes. Values for critical supersaturation derived from this experimental work were different for each tissue compartment. Values for slow tissues were close to Haldane's findings, but considerably higher values were found for the fast tissues. These values were so high that the researchers concluded that the 5- and 10-minute tissues were not relevant to the development of DCS. Based on these conclusions, a set of tables was computed which omitted the 5- and 10-minute tissues.

Yarbrough (1937 tables)

Yarbrough's 1937 tables were based on a Haldanean 3 compartment model with compartment half-times of 20, 40 and 70 minutes. Ascent rate was chosen to be 25 ft per minute, which was a convenient rate to pull up a diver in standard dress.

1956 tables

Van der Aue worked on procedures for surface decompression and oxygen use in the early 1950s and during his research found problems with the 1937 tables for long dive times. He also found that the fast tissues which had been dropped in the 1930s would control decompression in some cases, so he re-introduced the fast compartments to the model, and added an extra slower compartment to better model long-duration dives.
Assumptions of the 1956 model:
Ascent rate was chosen at 60 fsw/min as a compromise between the practical requirements for military scuba and surface supplied diving operations.
Repetitive diving was accommodated in the tables using the slowest compartment to control surface off-gassing.
A minimum surface interval of 10 minutes was found necessary to ensure that the 120-minute compartment would have controlling effect for repetitive dives.

U.S. Navy exceptional exposure tables

The U.S. Navy 1956 tables were soon found to be problematic for dives deeper than 100 fsw for longer than 2 to 4 hours.
US Navy exceptional exposure tables use an 8 compartment Haldanean model developed by Workman, with half-times of 5, 10, 20, 40, 80, 120, 160 and 240 minutes, and are not compatible with the rest of the US Navy Air tables for repetitive diving, although for convenience they have been appended to the standard US Navy Air tables. The tables warn that no repetitive diving is permitted following an exceptional exposure dive, and although the 240-minute tissue would only desaturate completely in 24 hours, there is no restriction to assuming an unsaturated diver after 12 hours.

Reformatting of the U.S. Navy 1956 tables by recreational diving community

Some of the earliest modifications to the U.S. Navy tables involved changes to their layout by the recreational diving community.
These include:
Decompression theory is not an exact science. Decompression models approximate a physiological process that is incompletely understood, and rather complex, by simple mathematical models, in the hope of producing a useful procedure with acceptably low risk of injury to the user. New information allows theories and models to be modified to provide more reliable results, and the availability of faster and more powerful computer processors at low cost has made more exhaustive numerical methods more practicable, and the computation of relatively far more complex models is now quite possible, even in real time.
Several factors have encouraged researchers to modify existing tables and develop new models:
Jeppesen made the simplest modification to the U.S. Navy tables by drawing a new line to reduce no-stop limits on an otherwise unchanged table. Divers were recommended to remain within the modified no-stop limit. If one of the new time limits was not listed on the U.S. Navy table, the next shorter table entry was to be selected.

Bassett tables

These tables were based on the 1956 US Navy tables and the no-decompression limits recommended by Bruce Bassett.
Changes were also made to the table rules and decompression requirements:
The first NAUI tables were based on reformatted but otherwise unmodified US Navy 1956 tables and issued in the early 1980s.
The next version was a NAUI modification of the US Navy 1956 tables using the following modifications, and released a few years later.
NAUI adapted the 1995 DCIEM Sports Table for use in all NAUI courses and these were used until they were replaced by RGBM based tables in 2002.
NAUI RGBM Trimix and Nitrox tables copyrighted 1999 have also been released.

Pandora tables

These tables were designed for use on the excavation of the wreck of the Pandora
In 1981 Karl Huggins modified the US Navy 6 compartment model using M values derived to follow the Spencer no-decompression limits. The tables are exclusively for no-decompression diving and are presented in the same format as the US Navy tables.
A major difference from the US Navy tables is that the repetitive group designators represent nitrogen levels in all tissues, unlike the USN table which represent only the 120-minute compartment. The Huggins repetitive group indicates a percentage of the M0 for the most saturated tissue, and this is intended to make the tables more applicable to multilevel diving procedures.
The Huggins tables have not been officially tested, but are more conservative than the 1956 US Navy tables. They have been calculated from limits which would theoretically produce venous bubbles 10 to 20% of the time.

Recreational Dive Planner, distributed by PADI

The tables known as the Recreational Dive Planner were developed and tested by Raymond Rogers and DSAT exclusively for no-stop diving. The M values were derived from Spencer's no-stop limits and the repetitive group designators were based on a 60-minute tissue compartment. This combination resulted in more conservative first dives, but less conservative repetitive dives.
The RDP tables were developed for no-stop diving, but recommend a safety stop at 15 fsw for 3 minutes. Emergency decompression for dives which inadvertently exceed the no-stop limit is specified.
The RDP tables are available in two formats:
The RDP was tested for single day multi-level dives and multi-day dives with multiple dives per day. There were no incidences of symptomatic DCS during testing.

Bühlmann tables

Professor A.A. Bühlmann of the Laboratory of Hyperbaric medicine of the Medical Clinic of the University of Zurich developed the Swiss tables, more often referred to as Bühlmann tables, in the early 1960s. The model is Haldanian, with 16 tissue compartments with half-times from 2.65 minutes to 635 minutes, each with linearly varying supersaturation limits depending on the tissue and the ambient pressure, and is based on absolute pressures, which simplifies application to altitude diving.
The full set of Swiss Tables consists of tables for four altitude ranges: 0 to 700 m, 701 to 1500 m, 1501 to 2500 m and 2501 to 3500 m. Ascent rate was chosen as 10 m per minute.
No stop limits and decompression schedules tend to be more conservative than the US Navy air tables.
The Swiss tables use the 80-minute tissue compartment for control of repetitive dive calculations, which tends to be less conservative than the US Navy tables for this application.

Modified Bühlmann tables

Swiss sport diving tables

In 1986 the Bühlmann model was used to generate dive tables for recreational divers. One set was for altitudes from 0 to 700m above sea level and other for altitudes from. The repetitive group designators are based on the 80-minute compartment.

Bühlmann/Hahn tables (German)

The German tables were developed by Dr. Max Hahn using a derivative of the Bühlmann ZH-L16 model using half-times ranging from 2.65 to 635 minutes. Three sets were published for altitude ranges 0–200 m, 201–700 m, and 701–1,200 m. The repetitive group designators are based on the 80-minute compartment.
Safety factors were added to the depths in the table to take into account depth gauge errors. The depths used for calculations were 2.4% greater than the depths listed in the two lower altitude tables, and 3% + 1 msw greater than the depths listed in the highest altitude table.

French Navy – Marine Nationale 90 (MN90) decompression tables

The mathematical model used for the development of the MN 90 tables is Haldanian, and was also used for the GERS 1965 table.
Haldane's assumptions about the limiting factors for ascent are:
Specific assumptions and conditions for use of the MN90 model and tables are as follows:

Royal Navy Physiological Laboratory model

In the early 1950s, Hempleman developed a diffusion limited model for gas transfer from the capillaries into the tissues. The basis for this model is radial diffusion from a capillary into the surrounding tissue, but by assuming closely packed capillaries in a plane the model was developed into a "tissue slab" equivalent to one-dimensional linear bulk diffusion in two directions into the tissues from a central surface.
The 1972 RNPL tables were based on a modified Hempleman tissue slab model and are more conservative than the US Navy tables.
A version of the RNPL tables was used by the British Sub-Aqua Club until the production of the BSAC'88 tables in 1988.

DCIEM model and tables

In the mid-1960s, the Canadian Defence and Civil Institute of Environmental Medicine developed the Kidd/Stubbs serial decompression model. This differs from Haldanian models which are parallel models and assume that all compartments are exposed to ambient partial pressures and no gas interchange occurs between compartments. A serial model assumes that the diffusion takes place through a series of compartments, and only one is exposed to the ambient partial pressures and is in effect a compartmentalised version of the Hempelman bulk diffusion slab model.
The Kidd/Stubbs model has four serial compartments, each with a half time of approximately 21 minutes. Allowable surfacing supersaturation ratios for the initial two compartments are taken as 1.92 and 1.73, while the gas concentration in the last two compartments is not considered in the computation.
DCIEM has continuously evaluated and modified the model over the years. A revised set of tables was released in 1984, based on thousands of Doppler evaluated dives.
The DCIEM 1983 decompression model is a decompression calculation model rather than a physiological model. Modifications were made to the model to get it to fit observed data, as the original model had several observed shortcomings, while retaining the basic model structure so that it could be applied to existing hardware with minimal modifications.

Mixed phase models (dissolved and bubble phases)

Thermodynamic model

LeMessurier and Hills published a paper in 1965 on A thermodynamic approach arising from a study on Torres Strait diving techniques which suggests that decompression by conventional models results in bubble formation which is then eliminated by re-dissolving at the decompression stops, which is slower than elimination while still in solution, thus indicating the importance of minimising bubble phase for efficient gas elimination.

Tables du Ministère du Travail

Tables du Ministère du Travail 1974 (MT74)

The first French official air decompression tables were published in 1974 by the Ministère du Travail

Tables du Ministère du Travail 1992 (MT92)

In 1982, the French government funded a research project for the evaluation of the MT74 tables using computer analysis of the dive report database, which indicated that the MT74 tables had limitations for severe exposures. The government then supported a second project to develop and validate new tables. A complete set of air tables, with options of pure oxygen breathing at 6 m, at 12 m, surface decompression, split level diving, repetitive diving, etc. was developed in 1983. This early model already implemented the concept of continuous compartment half-times. For the safe ascent criteria, the Arterial Bubble model was not derived mathematically, but an approximation was defined empirically by fitting mathematical expressions to selected exposures from the Comex database. At the time, the best fit was obtained by the expression now called AB Model-1, which was used to compute a set of decompression tables that was evaluated offshore on selected Comex worksites. In 1986, after some minor adjustments, the tables were included in the Comex diving manuals and used as standard procedures. In 1992, the tables were included in the new French diving regulations under the name of Tables du Ministère du Travail 1992 or MT92 tables

The arterial bubble decompression model

The arterial bubble assumption is that the filtering capacity of the lung has a threshold radius of the size of a red blood cell and that sufficiently small decompression bubbles can pass to the arterial side, especially during the initial phase of ascent. Later in the ascent, bubbles grow to a larger size and remain trapped in the lung. This may explain why conventional Doppler measurements have not detected any bubbles in the arterial circulation.
The arterial bubble assumption can introduce variability in the decompression outcome through the lung function. The first variable is individual susceptibility. The filtering capacity of the lung may be assumed to vary between individuals, and for a given individual, from day to day, and may account for the inter-personal and intra-persona variability which have been observed in DCS susceptibility. Basically, a good diver is a good bubble filter. This is a justification for divers who seek top physical fitness for severe decompression exposures.
The second variable is related to dive conditions and speculates an influence of CO2 on the lung filter. Raised levels of CO2 could decrease the lungs' filtration capacity and allow bubbles to pass to the arterial side of the circulation. Thus, diving situations associated with CO2 retention and hypercapnia would be associated with a higher risk of Type II DCS. This could explain why the following situations, which are all related to high levels of CO2, have been identified as contributing factors to DCS:
The arterial bubble assumption is also consistent with the accidental production of arterial bubbles. One scenario considers a shunt at the heart or lung level that passes bubbles from the venous to the arterial side. A patent foramen ovale is thought to only open in certain conditions. A PFO conveniently explains neurological accidents after recreational air diving without any procedure violation, but it does not explain vestibular hits in deep diving. Vestibular symptoms can appear very early in the decompression, long before the massive bubble production required to overload the system.
A second scenario considers pressure increases during decompression that reduce bubble diameters. This can allow bubbles trapped in the lung during a normal decompression to suddenly pass through the capillaries and become responsible for Type II DCS symptoms. This could explain the difference in outcomes of in-water decompression versus surface decompression. Data collected in the North Sea have shown that if the overall incidence rate of the two diving methods is about the same, that surface decompression tends to produce ten times more type II DCS than in-water decompression. It is assumed that when the diver ascends to the surface, bubbles are produced that are trapped by the lung capillaries, and on recompression of the diver in the deck chamber, these bubbles are reduced in diameter and pass to the arterial side, later causing neurological symptoms. The same scenario was proposed for type II DCS recorded after sawtooth diving profiles or multiple repetitive dives.
The arterial bubble assumption also provides an explanation for the criticality of the initial ascent phase. Bubbles associated with symptoms are not necessarily generated on site. There is a growth process at the beginning of the ascent that may last for several cycles until the bubbles have reached a critical size when they are either filtered in the lung or stopped at the tissue level. It is postulated that the production of a shower of small arterial bubbles during the first minutes of the initial ascent is a precursor for DCS symptoms.
An attempt was made to turn this scenario into a decompression model.
The arterial bubble model assumptions
  1. A Diver breathes a compressed gas mixture that contains inert gas which dissolves in the various tissues during the pressure exposure. When the ascent is initiated, the inert gas is off-loaded as soon as a suitable gradient is created.
  2. Bubbles are normally produced in the vascular bed and transported by the venous system to heart, then to the lungs.
  3. The lungs work as a filter and trap the bubbles in the capillaries which have a smaller diameter. Gas transfer into the alveoli eliminates the bubbles.
  4. The critical issue is the filtering capacity of the lung system. Small bubbles may pass through the lungs into the systemic circulation.
  5. At the level of the aortic arch, the distribution of blood likely to carry bubbles to neurological tissue such as the brain or the spinal cord.
  6. The brain is a fast tissue and might be in supersaturated state in the early phase of decompression. It acts as a gas reservoir and feeds any local bubble which will grow. The bubble may just proceed through the capillaries to the venous side for another cycle, but may be trapped and will then grow in place, causing local restriction of the blood supply and finally ischemia. This may develop into central neurological symptoms.
  7. Similarly, arterial bubbles may reach the spinal cord and grow on site from local gas and produce spinal neurological symptoms.
  8. Much later in the decompression, bubbles may reach a significant size and exert a local deformation, particularly in stiffer tissues such as tendons and ligaments, that excites nerve terminations and produces pain.
Derivation of the Arterial Bubble Model
A model based on the Arterial Bubble assumption was developed for the calculation of decompression tables.
This gas phase model uses an equation which can be compared to a classic "M-value" associated with a corrective factor that reduces the permitted gradient for small values of the compartment time constant.
The consequence is the introduction of deeper stops than a classic dissolved phase decompression model.
The rationalization of the arterial bubble assumption considers two situations:
The critical volume concept was developed by Hennessy and Hempleman who formulated a simple mathematical condition linking the dissolved gas and the safe ascent pressure:
Where Ptissue represents the dissolved gas tension, Pambient, the ambient pressure and a and b two coefficients. This linear relationship between dissolved gas and ambient pressure has the same mathematical form as an M value, which suggests that all the Haldanean models using M-values, may be considered expressions of the critical volume criterion, though their authors may have argued for other interpretations.

U.S. Navy E-L algorithm and the 2008 tables

The use of simple symmetrical exponential gas kinetics models has shown up the need for a model that would give slower tissue washout. In the early 1980s the US Navy Experimental Diving Unit developed an algorithm using a decompression model with exponential gas absorption as in the usual Haldanian model, but a slower linear release during ascent. The effect of adding linear kinetics to the exponential model is to lengthen the duration of risk accumulation for a given compartment time constant
The model was originally developed for programming decompression computers for constant oxygen partial pressure closed circuit rebreathers. Initial experimental diving using an exponential-exponential algorithm resulted in an unacceptable incidence of DCS, so a change was made to a model using the linear release model, with a reduction in DCS incidence.
The same principles were applied to developing an algorithm and tables for a constant oxygen partial pressure model for heliox diving
The linear component is active when the tissue pressure exceeds ambient pressure by a given amount specific to the tissue compartment. When the tissue pressure drops below this cross-over criterion the tissue is modelled by exponential kinetics. During gas uptake, tissue pressure never exceeds ambient, so it is always modelled by exponential kinetics. This results in a model with the desired asymmetrical characteristics of slower washout than uptake.
The linear/exponential transition is smooth. Choice of cross-over pressure determines the slope of the linear region as equal to the slope of the exponential region at the cross-over point.
During the development of these algorithms and tables, it was recognized that a successful algorithm could be used to replace the existing collection of incompatible tables for various air and Nitrox diving modes currently in the U.S. Navy Diving Manual with a set of mutually compatible decompression tables based on a single model, which was proposed by Gerth and Doolette in 2007. This has been done in Revision 6 of the US Navy Diving Manual published in 2008, though some changes were made.
An independent implementation of the EL-Real Time Algorithm was developed by Cochran Consulting, Inc. for the diver-carried Navy Dive Computer
under the guidance of E. D. Thalmann.

Physiological interpretation

Computer testing of a theoretical bubble growth model reported by Ball, Himm, Homer and Thalmann produced results which led to the interpretation of the three compartments used in the probabilistic LE model, with fast, intermediate and slow time constants, of which only the intermediate compartment uses the linear kinetics modification during decompression, as possibly not representing distinct anatomically identifiable tissues, but three different kinetic processes which relate to different elements of DCS risk.
They conclude that bubble evolution may not be sufficient to explain all aspects of DCS risk, and the relationship between gas phase dynamics and tissue injury requires further investigation.

BSAC '88 Tables

The BSAC '88 Tables are published in the form of a booklet of four table sets giving no calculation repetitive diving solutions from sea level to 3000 metres altitude.
These tables were developed by Tom Hennessy to replace the RNPL/BSAC tables when the Club wanted a set of tables which could approach the versatility of a dive computer.
Very little information on the theoretical model and algorithm for the BSAC 1988 tables appears to be available.
What is known, is that the tables were developed specifically for recreational diving for the British Sub-Aqua Club by Dr Tom Hennessy and were released in 1988.
Also in 1988, a chapter titled Modelling Human Exposure to Altered Pressure Environments, by T.R. Hennessy was published in Environmental Ergonomics, discussing the shortcomings of several decompression models and the associated experimental validation procedures.
In this work Hennessy proposes an alternative combined perfusion/diffusion model. The number of compartments discussed ranges from 4 in model "A", to 2 in model "B".
Hennessy concludes that if the undissolved and dissolved gas content of a tissue cannot be independently measured either directly or indirectly then the safe maximum limits relative to the ambient pressure cannot be accurately determined through decompression trials and it will not be possible to systematically develop a comprehensive biophysical model for gas exchange. He proposes a best fit double compartment model for dissolved gas and a single compartment model for undissolved gas as these are the simplest models consistent with available data.
The parameters used in the design of these tables include:
The initial dive uses table A, and the diver is allocated a Surfacing Code based on depth and time of the dive. After a surface interval of at least 15 minutes the diver can select a new Current Tissue Code which models the residual nitrogen load, and uses this code to select the repetitive dive table.
The BSAC'88 tables are presented in a format which does not require any calculation by the user.

Varying Permeability Model

This decompression model was developed by D.E. Yount and others at the University of Hawaii to model laboratory observations of bubble formation and growth in both inanimate and in vivo systems exposed to pressure variations. It presumes that microscopic bubble nuclei always exist in aqueous media, including living tissues. These bubble nuclei are spherical gas phases that are small enough to remain in suspension yet strong enough to resist collapse, their stability being provided by an elastic surface layer consisting of surface-active molecules with variable gas permeability. These skins resist the effect of surface tension, as surface tension tends to collapse a small bubble by raising internal pressure above ambient so that the partial pressure gradient favours diffusion out of the bubble in inverse proportion to the radius of the surface.
Any nuclei larger than a specific "critical" size, will grow during decompression. The VPM aims to limit the cumulative volume of these growing bubbles during and after decompression to a tolerable level by limiting the pressure difference between the gas in the bubbles and the ambient pressure. In effect, this is equivalent to limiting the supersaturation, but instead of using an arbitrary linear fit to experimental data, the physics of bubble growth is used to model the acceptable supersaturation for any given pressure exposure history.
Growth in size and number of gas bubbles is computed based on factors representing pressure balances in the bubbles, physical properties of the "skins" and the surrounding environment. If the total volume of gas in the bubbles is predicted to be less than a "critical volume", then the diver is assumed to be within the safe limits of the model.
The bubble model is superposed on a multiple parallel tissue compartment model. Ingassing is assumed to follow the classic Haldanean model.

Bubble population distribution

Bubble size vs number has an exponential distribution

Bubble nucleation

Gas bubbles with a radius greater than 1 micron should float to the surface of a standing liquid, whereas smaller ones should dissolve rapidly due to surface tension. The Tiny Bubble Group has been able to resolve this apparent paradox by developing and experimentally verifying a new model for stable gas nuclei.
According to the varying-permeability model, gas bubble nuclei are simply stable microbubbles. The stability of these microbubbles is due to elastic skins or membranes consisting of surface-active molecules. These skins are normally permeable to gas, and collapse is prevented by their compression strength. These skins can become stiff and effectively impermeable to gas when they are subjected to large compressions, typically exceeding 8 atm, at which stage the pressure inside increases during further compression as predicted by Boyle's law.
Essentially, there are three parameters in the VP model:
the compression strength of the skin; the initial radius; and the onset pressure for impermeability.

Ordering hypothesis

The ordering hypothesis states that nuclei are neither created nor destroyed by the pressure schedule, and initial ordering according to size is preserved.
It follows from the ordering hypothesis that each bubble count is determined by the properties and behavior of that one "critical" nucleus which is right at the bubble formation threshold.
All nuclei that are larger than the critical nucleus will form bubbles, and all nuclei that are smaller will not. Furthermore, a family of pressure schedules which yields the same bubble count N is characterized by the same critical nucleus and hence by the same critical radius, the same crumbling compression, and the same onset of impermeability.

Development of decompression model

The original assumption was that bubble number is directly proportional to decompression stress. This approach worked well for long exposures, but not when the exposure time varied considerably.
A better model was obtained by allowing more bubbles to form on the shorter dives than on the longer dives. The constant bubble number assumption was replaced by a "dynamic-critical-volume hypothesis". As in earlier applications of the critical-volume criterion, it was assumed that whenever the total volume of gas phase accumulated exceeds a critical value, signs or symptoms of DCS will appear. In the special case of long exposures the two models are equivalent.
The "dynamic" aspect of this hypothesis is that gas is continuously entering and leaving the gas phase.
The accumulated volume is calculated as a function of time by integrating over the product of the bubble number and the degree of supersaturation, and subtracting the free gas that is being dissipated continuously by the lung.
Gas uptake and elimination are assumed to be exponential, as in conventional Haldanean models.
As a first approximation only the inert gasses are taken into account. For oxygen partial pressures above 2.4 bar, the quantity of oxygen dissolved in the arterial blood exceeds the amount that the body can use, and the hemoglobin is saturated with oxygen in both the veins and the arteries. If more oxygen is added, the partial pressure of oxygen in the venous blood rises.

Comparison of VPM profiles with other models

Comparisons of VPM profiles with USN decompression schedules for extreme exposure dives consistently produce similar total ascent times, but significantly deeper first decompression stops.

Reduced Gradient Bubble Model

The RGBM developed by Dr Bruce Wienke at Los Alamos National Laboratory is a hybrid model which modifies a Haldanian model with factors to take some account of bubble mechanics to model gas phase production during decompression. The bubble factor modifies the M-values of the Haldanian model, making it more conservative.
Features of the modifying factor ξ include:
The effect is to reduce no-stop dive time or increase decompression requirements for repetitive dive in the following categories:
The model has been used to some extent in some Suunto dive computers, and in the HydroSpace Explorer computer, where it is a user selected option for computation formula, with a choice of additional conservatism factors.
The complete RGBM treats coupled perfusion-diffusion transport as a two-stage process, with perfusion providing a boundary condition for gas penetration of the tissues by diffusion. Either process can dominate the exchange depending on time and rate coefficients.
Simplified implementations which require less computational power are available for use in personal decompression computers. These are dominated by perfusion. The inherent biological unsaturation of tissues is considered in the calculations.
The model assumes that bubble nuclei are always present in a specific size distribution, and that a certain number are induced to grow by compression and decompression. An iterative computation is used to model ascent to limit the combined volume of the gas phase. Gas mixtures of helium, nitrogen, and oxygen contain bubble distributions of different sizes, but the same phase volume limit is used.
The model postulates bubble nuclei with aqueous and/or lipid skin structure, in a number and size distribution quantified by an equation-of-state. Like the VPM, RGBM assumes the size distribution is exponentially decreasing in size. Unlike the varying permeability model, bubble seeds are assumed permeable to gas transfer across skin boundaries under all pressures.
The size of nuclei which will grow during decompression is inversely proportional to the supersaturation gradient.
At higher pressures, skin tension of the bubble nuclei reduces gas diffusion to a slower rate. The model assumes that bubble skins are stabilized by surfactants over calculable times scales, which results in variable persistence of the bubble nuclei in the tissues.

Modifications to models and algorithms for diluent gases other than nitrogen

Decompression models and algorithms developed for binary mixtures of nitrogen and oxygen can not be used for gases containing significant amounts of other diluent gases without modification to take into account the different solubilities and diffusion constants of the alternative or added diluents. It is also highly desirable to test any such modifications, to make sure the schedules produced by them are acceptably safe.

Alternative diluent gases

To a large extent commercial offshore diving uses heliox tables that have been developed by the major commercial diving enterprises such as Comex, Oceaneering International Alpha tables, American Oilfield Diving Company gas tables, though modifications of the US Navy Partial pressure tables are also used. In 2006 the unmodified US Navy tables were considered to result in an unacceptably high rate of decompression sickness for commercial applications.
"Cx70" heliox tables were developed and used by Comex between 1970 and 1982. The tables were available in two versions. One was designed for surface-supplied diving and limited to 75 m. The diver breathed heliox as the bottom mix and 100% oxygen at the 6 m stop. The other was designed for closed bell bounce diving and allowed for exposures up to 120 minutes, and depths to 120 m. The diver breathed heliox in the water and in the bell, air after transfer into the deck decompression chamber, and finally oxygen on built in breathing system from 12 m to the surface. These tables produced a relatively high incidence of decompression sickness.
The French Tables du Ministère du Travail 1974 and Tables du Ministère du Travail 1992 were developed specifically for commercial diving.
Norwegian Diving and Treatment Tables,, referenced in NORSOK Standard U100 2.24 for manned underwater operations, are available in Norwegian, Danish and English text and are approved for commercial diving.

Other reading