Fatigue microdamage in bovine trabecular bone

Microdamage, in the form of small cracks, may accumulate in trabecular bone loaded in fatigue. Specimens of bovine trabecular bone were loaded in compressive fatigue at one of four normalized stresses and loading was stopped after the specimens reached one of six maximum strains. Microdamage was identified using a fluorochrome staining technique, and microdamage parameters, including the number of damaged trabeculae and the damaged area fraction, were measured. No microdamage was observed during loading to strains below the yield strain; at higher strains, all microdamage parameters increased with increasing maximum compressive strain. Few significant differences were observed in the type or amount of microdamage accumulation between specimens loaded to the same maximum strain at different normalized stresses; however, more trabecular fractures were observed at high numbers of cycles, which corresponded to low normalized stresses.     

Trabecular bone microdamage and microstructural stresses under uniaxial compression

The balance between local remodeling and accumulation of trabecular bone microdamage is believed to play an important role in the maintenance of skeletal integrity. However, the local mechanical parameters associated with microdamage initiation are not well understood. Using histological damage labeling, micro-CT imaging, and image-based finite element analysis, regions of trabecular bone microdamage were detected and registered to estimated microstructural von Mises effective stresses and strains, maximum principal stresses and strains, and strain energy density (SED). Bovine tibial trabecular bone cores underwent a stepwise uniaxial compression routine in which specimens were micro-CT imaged following each compression step. The results indicate that the mode of trabecular failure observed by micro-CT imaging agreed well with the polarity and distribution of stresses within an individual trabecula. Analysis of on-axis subsections within specimens provided significant positive relationships between microdamage and each estimated tissue stress, strain and SED parameter. In a more localized analysis, individual microdamaged and undamaged trabeculae were extracted from specimens loaded within the elastic region and to the apparent yield point. As expected, damaged trabeculae in both groups possessed significantly higher local stresses and strains than undamaged trabeculae. The results also indicated that microdamage initiation occurred prior to apparent yield at local principal stresses in the range of 88-121 MPa for compression and 35-43 MPa for tension and local principal strains of 0.46-0.63% in compression and 0.18-0.24% in tension. These data provide an important step towards understanding factors contributing to microdamage initiation and establishing local failure criteria for normal and diseased trabecular bone.     

Detection of fatigue microdamage in whole rat femora using contrast-enhanced micro-computed tomography

Microdamage in bone tissue is typically studied using destructive, two-dimensional histological techniques. Contrast-enhanced micro-computed tomography (micro-CT) was recently demonstrated to enable non-destructive, three-dimensional (3-D) detection of microdamage in machined cortical and trabecular bone specimens in vitro. However, the accumulation of microdamage in whole bones is influenced by variations in the magnitude and mode of loading due to the complex whole bone morphology. Therefore, the objective of this study was to detect the presence, spatial location, and accumulation of fatigue microdamage in whole rat femora in vitro using micro-CT with a BaSO₄ contrast agent. Microdamage was detected and observed to accumulate at specific spatial locations within the cortex of femora loaded in cyclic three-point bending to a 5% or 10% reduction in secant modulus. The ratio of the segmented BaSO₄ stain volume (SV) to the total volume (TV) of cortical bone was adopted as a measure of damage. The amount of microdamage measured by micro-CT (SV/TV) was significantly greater for both loaded groups compared to the control group (p<0.05), but the difference between loaded groups was not statistically significant. At least one distinct region of microdamage, as indicated by the segmented SV, was observed in 85% of loaded specimens. A specimen-specific finite element model confirmed elevated tensile principal strains localized in regions of tissue corresponding to the accumulated microdamage. These regions were not always located where one might expect a priori based upon Euler–Bernoulli beam theory, demonstrating the utility of contrast-enhanced micro-CT for non-destructive, 3-D detection of fatigue microdamage in whole bones in vitro.     

The effects of tensile-compressive loading mode and microarchitecture on microdamage in human vertebral cancellous bone

The amount of microdamage in bone tissue impairs mechanical performance and may act as a stimulus for bone remodeling. Here we determine how loading mode (tension vs. compression) and microstructure (trabecular microarchitecture, local trabecular thickness, and presence of resorption cavities) influence the number and volume of microdamage sites generated in cancellous bone following a single overload. Twenty paired cylindrical specimens of human vertebral cancellous bone from 10 donors (47–78 years) were mechanically loaded to apparent yield in either compression or tension, and imaged in three dimensions for microarchitecture and microdamage (voxel size 0.7×0.7×5.0 μm³). We found that the overall proportion of damaged tissue was greater (p=0.01) for apparent tension loading (3.9±2.4%, mean±SD) than for apparent compression loading (1.9±1.3%). Individual microdamage sites generated in tension were larger in volume (p<0.001) but not more numerous (p=0.64) than sites in compression. For both loading modes, the proportion of damaged tissue varied more across donors than with bone volume fraction, traditional measures of microarchitecture (trabecular thickness, trabecular separation, etc.), apparent Young?s modulus, or strength. Microdamage tended to occur in regions of greater trabecular thickness but not near observable resorption cavities. Taken together, these findings indicate that, regardless of loading mode, accumulation of microdamage in cancellous bone after monotonic loading to yield is influenced by donor characteristics other than traditional measures of microarchitecture, suggesting a possible role for tissue material properties.     

Spatial correlations of trabecular bone microdamage with local stresses and strains using rigid image registration

Although microdamage is known to accumulate in trabecular bone with overloading and aging, the tissue-level stresses and strains associated with local bone failure are not well known. Local correlation of microdamage with microstructural stresses and strains requires methods to accurately register histological sections with micro-computed tomography (micro-CT) based finite element models. In addition, the resolution of correlation (i.e., grid size) selected for analysis may affect the observed results. Therefore, an automated, repeatable, and accurate image registration algorithm was developed to determine the range of local stresses and strains associated with microdamage initiation. Using a two-dimensional rigid registration algorithm, bone structures from histology and micro-CT imaging were aligned. Once aligned, microdamaged regions were spatially correlated with local stresses and strains obtained from micro-CT based finite element analysis. Using this more sophisticated registration technique, we were able to analyze the effects of varying spatial grid resolution on local stresses and strains initiating microdamage. The results indicated that grid refinement to the individual pixel level (pixel-by-pixel method) more precisely defined the range of microdamage initiation compared to manually selected individual damaged and undamaged trabeculae. Using the pixel-by-pixel method, we confirmed that trabecular bone from younger cows sustained higher local strains prior to microdamage initiation compared to older bone.     

Synchrotron radiation micro-CT at the micrometer scale for the analysis of the three-dimensional morphology of microcracks in human trabecular bone

Bone quality is an important concept to explain bone fragility in addition to bone mass. Among bone quality factors, microdamage which appears in daily life is thought to have a marked impact on bone strength and plays a major role in the repair process. The starting point for all studies designed to further our understanding of how bone microdamage initiate or dissipate energy, or to investigate the impact of age, gender or disease, remains reliable observation and measurement of microdamage. In this study, 3D Synchrotron Radiation (SR) micro-CT at the micrometric scale was coupled to image analysis for the three-dimensional characterization of bone microdamage in human trabecular bone specimens taken from femoral heads. Specimens were imaged by 3D SR micro-CT with a voxel size of 1.4 μm. A new tailored 3D image analysis technique was developed to segment and quantify microcracks. Microcracks from human trabecular bone were observed in different tomographic sections as well as from 3D renderings. New 3D quantitative measurements on the microcrack density and morphology are reported on five specimens. The 3D microcrack density was found between 3.1 and 9.4/mm3 corresponding to a 2D density between 0.55 and 0.76 /mm2. The microcrack length and width measured in 3D on five selected microcrack ranged respectively from 164 μm to 209 μm and 100 μm to 120 μm. This is the first time that various microcracks in unloaded human trabecular bone--from the simplest linear crack to more complex cross-hatch cracks--have been examined and quantified by 3D imaging at this scale. The suspected complex morphology of microcracks is here considerably more evident than in the 2D observations. In conclusion, this technique opens new perspective for the 3D investigation of microcracks and the impact of age, disease or treatment.     

Heterogeneous glycation of cancellous bone and its association with bone quality and fragility

Non-enzymatic glycation (NEG) and enzymatic biochemical processes create crosslinks that modify the extracellular matrix (ECM) and affect the turnover of bone tissue. Because NEG affects turnover and turnover at the local level affects microarchitecture and formation and removal of microdamage, we hypothesized that NEG in cancellous bone is heterogeneous and accounts partly for the contribution of microarchitecture and microdamage on bone fragility. Human trabecular bone cores from 23 donors were subjected to compression tests. Mechanically tested cores as well as an additional 19 cores were stained with lead-uranyl acetate and imaged to determine microarchitecture and measure microdamage. Post-yield mechanical properties were measured and damaged trabeculae were extracted from a subset of specimens and characterized for the morphology of induced microdamage. Tested specimens and extracted trabeculae were quantified for enzymatic and non-enzymatic crosslink content using a colorimetric assay and Ultra-high Performance Liquid Chromatography (UPLC). Results show that an increase in enzymatic crosslinks was beneficial for bone where they were associated with increased toughness and decreased microdamage. Conversely, bone with increased NEG required less strain to reach failure and were less tough. NEG heterogeneously modified trabecular microarchitecture where high amounts of NEG crosslinks were found in trabecular rods and with the mechanically deleterious form of microdamage (linear microcracks). The extent of NEG in tibial cancellous bone was the dominant predictor of bone fragility and was associated with changes in microarchitecture and microdamage.     

Trabecular bone strain changes associated with subchondral stiffening of the proximal tibia

Subchondral stiffening is a hallmark pathologic feature of osteoarthritis but its mechanical and temporal relationship to the initiation or the progression of osteoarthritis is not established. The mechanical effect of subchondral stiffening on the surrounding trabecular bone is poorly understood. This study employs a relatively new application of digital image correlation to measure strain in the trabecular region of the proximal medial tibia in normal specimens and in specimens with simulated subchondral bone stiffening. Coronal sections from eight normal human cadaveric proximal tibiae were loaded in static compression and high resolution contact radiographs were made. Repeat contact radiographs were collected after the subchondral bone near the jointline was stiffened using polymethylmethacrylate. Digital images, made from loaded and unloaded contact radiographs, were compared using in-house software to measure trabecular displacement and calculate trabecular bone strain. Overall strain was higher in the stiffened specimens suggesting experimental artifiact significantly affected our results. Consistent increases in median maximum shear strain, median maximum principal strain, median minimum principal strain, and peak shear strain were measured near the inner and outer edges of the stiffened segment. Our experiment provides direct experimental measurement of increases in trabecular bone strain caused by subchondral stiffening, however, the clinical and biologic importance of strain increases is unknown.     

Development of a fluorescent light technique for evaluating microdamage in bone subjected to fatigue loading.

A new method using fluorescent light microscopy has been developed to visualize and evaluate bone microdamage. We report the findings of two different experiments with a common aim of comparing the fluorescent light technique to the brightfield method for quantifying microdamage in bone. In Experiment 1, 36 canine femurs were tested in four-point cyclic bending until they had lost between 5 and 43% of their stiffness. The loaded portion of the bone was stained en bloc with basic fuchsin for the presence of damage. Standard point counting techniques were used to calculate fractional damaged area (Dm.Ar = Cr.Ar/B.Ar, mm2/mm2) under brightfield and fluorescent microscopy. In Experiment 2, bone microdamage adjacent to endosseous implants, subjected to fatigue loading (150,000 cycles, 2 Hz and 37 degrees C) ex vivo was examined. The bone around the implant was either allowed to heal (adapted specimen) for 12 weeks after placement in dog mid-femoral diaphyses prior to testing or was loaded immediately to simulate non-healed bone surrounding endosseous implants (non-adapted). Crack numerical density (Cr.Dn = Cr.N/B.Ar, #/mm2), crack surface density (Cr.S.Dn = Tt.Cr.Le/B.Ar, mm/mm2) and fractional damaged area were calculated separately by both techniques in the adapted and non-adapted specimens. In both Experiments 1 and 2, significantly more microdamage was detected by the fluorescent technique than by the brightfield method. Also, there was a trend towards higher intraobserver repeatability when using the fluorescent method. These results suggest that the brightfield technique underestimates microdamage accumulation and that the fluorescent technique better represents the actual amounts of microdamage present. The results demonstrate that the fluorescent method provides an accurate and precise approach for bone microdamage evaluation, and that it improves the prediction of stiffness loss from damage accumulation.     

A Comparison of Micro-CT and Dental CT in Assessing Cortical Bone Morphology and Trabecular Bone Microarchitecture

Objective

The objective of this study was to evaluate the relationship between the trabecular bone microarchitecture and cortical bone morphology by using micro-computed tomography (micro-CT) and dental cone-beam computed tomography (dental CT).

Materials and Methods

Sixteen femurs and eight fifth lumbar vertebrae were collected from eight male Sprague Dawley rats. Four trabecular bone microarchitecture parameters related to the fifth lumbar vertebral body (percent bone volume [BV/TV], trabecular thickness [TbTh], trabecular separation [TbSp], and trabecular number [TbN]) were calculated using micro-CT. In addition, the volumetric cancellous bone grayscale value (vCanGrayscale) of the fifth lumbar vertebral body was measured using dental CT. Furthermore, four cortical bone morphology parameters of the femoral diaphysis (total cross-sectional area [TtAr], cortical area [CtAr], cortical bone area fraction [CtAr/TtAr], and cortical thickness [CtTh]) were calculated using both micro-CT and dental CT. Pearson analysis was conducted to calculate the correlation coefficients (r) of the micro-CT and dental CT measurements. Paired-sample t tests were used to compare the differences between the measurements of the four cortical bone morphology parameters obtained using micro-CT and dental CT.

Results

High correlations between the vCanGrayscale measured using dental CT and the trabecular bone microarchitecture parameters (BV/TV [r = 0.84] and TbTh [r = 0.84]) measured using micro-CT were observed. The absolute value of the four cortical bone morphology parameters may be different between the dental CT and micro-CT approaches. However, high correlations (r ranged from 0.71 to 0.90) among these four cortical bone morphology parameters measured using the two approaches were obtained.

Conclusion

We observed high correlations between the vCanGrayscale measured using dental CT and the trabecular bone microarchitecture parameters (BV/TV and TbTh) measured using micro-CT, in addition to high correlations between the cortical bone morphology measured using micro-CT and dental CT. Further experiments are necessary to validate the use of dental CT on human bone.     

Effects of thresholding techniques on microCT-based finite element models of trabecular bone

Microimaging based finite element analysis is widely used to predict the mechanical properties of trabecular bone. The choice of thresholding technique, a necessary step in converting grayscale images to finite element models, can significantly influence the predicted bone volume fraction and mechanical properties. Therefore, we investigated the effects of thresholding techniques on microcomputed tomography (micro-CT) based finite element models of trabecular bone. Three types of thresholding techniques were applied to 16-bit micro-CT images of trabecular bone to create three different models per specimen. Bone volume fractions and apparent moduli were predicted and compared to experimental results. In addition, trabecular tissue mechanical parameters and morphological parameters were compared among different models. Our findings suggest that predictions of apparent mechanical properties and structural properties agree well with experimental measurements regardless of the choice of thresholding methods or the format of micro-CT images.     

Fatigue microcracks that initiate fracture are located near elevated intracortical porosity but not elevated mineralization

In vivo microcracks in cortical bone are typically observed within more highly mineralized interstitial tissue, but postmortem investigations are inherently limited to cracks that did not lead to fracture which may be misleading with respect to understanding fracture mechanisms. We hypothesized that the one fatigue microcrack which initiates fracture is located spatially adjacent to elevated intracortical porosity but not elevated mineralization. Therefore, the spatial correlation between intracortical porosity, elevated mineralization, and fatigue microdamage was investigated by combining, for the first time, sequential, nondestructive, three-dimensional micro-computed tomography (micro-CT) measurements of each in cortical bone specimens subjected to compressive fatigue loading followed by a tensile overload to fracture. Fatigue loading resulted in significant microdamage accumulation and compromised mechanical properties upon tensile overload compared to control specimens. The microdamage that initiated fracture upon tensile overload was able to be identified in all fatigue-loaded specimens using contrast-enhanced micro-CT and registered images. Two-point (or pair) correlation functions revealed a spatial correlation between microdamage at the fracture initiation site and intracortical porosity, but not highly mineralized tissue, confirming the hypothesis. This difference was unique to the fracture initiation site. Intracortical porosity and highly mineralized tissue exhibited a significantly lower and higher probability, respectively, of being located spatially adjacent to all sites of microdamage compared to the fracture initiation site. Therefore, the results of this study suggest that human cortical bone is tolerant of most microcracks, which are generally compartmentalized within the more highly mineralized interstitial tissue, but a single microcrack of sufficient size located in spatial proximity to intracortical porosity can compromise fracture resistance.     

Shear strength and toughness of trabecular bone are more sensitive to density than damage

Microdamage occurs in trabecular bone under normal loading, which impairs the mechanical properties. Architectural degradation associated with osteoporosis increases damage susceptibility, resulting in a cumulative negative effect on the mechanical properties. Treatments for osteoporosis could be targeted toward increased bone mineral density, improved architecture, or repair and prevention of microdamage. Delineating the relative roles of damage and architectural degradation on trabecular bone strength will provide insight into the most beneficial targets. In this study, damage was induced in bovine trabecular bone samples by axial compression, and the effects on the mechanical properties in shear were assessed. The damaged shear modulus, shear yield stress, ultimate shear stress, and energy to failure all depended on induced damage and decreased as the architecture became more rod-like. The changes in ultimate shear strength and toughness were proportional to the decrease in shear modulus, consistent with an effective decrease in the cross-section of trabeculae based on cellular solid analysis. For typical ranges of bone volume fraction in human bone, the strength and toughness were much more sensitive to decreased volume fraction than to induced mechanical damage. While ultimately repairing or avoiding damage to the bone structure and increasing bone density both improve mechanical properties, increasing bone density is the more important contributor to bone strength.     

Age-related differences in the morphology of microdamage propagation in trabecular bone

Microdamage density has been shown to increase with age in trabecular bone and is associated with decreased fracture toughness. Numerous studies of crack propagation in cortical bone have been conducted, but data in trabecular bone is lacking. In this study, propagation of severe, linear, and diffuse damage was examined in trabecular bone cores from the femoral head of younger (61.3±3.1 years) and older (75.0±3.9 years) men and women. Using a two-step mechanical testing protocol, damage was first initiated with static uniaxial compression to 0.8% strain then propagated at a normalized stress level of 0.005 to a strain endpoint of 0.8%. Coupling mechanical testing with a dual-fluorescent staining technique, the number and length/area of propagating cracks were quantified. It was found that the number of cycles to the test endpoint was substantially decreased in older compared to younger samples (younger: 77,372±15,984 cycles; older: 34,944±11,964 cycles, p=0.06). This corresponded with a greater number of severely damaged trabeculae expanding in area during the fatigue test in the older group. In the younger group, diffusely damaged trabeculae had a greater damage area, which illustrates an efficient energy dissipation mechanism. These results suggest that age-related differences in fatigue life of human trabecular bone may be due to differences in propagated microdamage morphology.     

Theoretical analysis of alendronate and risedronate effects on canine vertebral remodeling and microdamage

Bisphosphonates suppress bone remodeling activity, increase bone volume, and significantly reduce fracture risk in individuals with osteoporosis and other metabolic bone diseases. The objectives of the current study were to develop a mathematical model that simulates control and 1 year experimental results following bisphosphonate treatment (alendronate or risedronate) in the canine fourth lumbar vertebral body, validate the model by comparing simulation predictions to 3 year experimental results, and then use the model to predict potential long term effects of bisphosphonates on remodeling and microdamage accumulation. To investigate the effects of bisphosphonates on bone volume and microdamage, a mechanistic biological model was modified from previous versions to simulate remodeling in a representative volume of vertebral trabecular bone in dogs treated with various doses of alendronate or risedronate, including doses equivalent to those used for treatment of post-menopausal osteoporosis in humans. Bisphosphonates were assumed to affect remodeling by suppressing basic multicellular unit activation and reducing resorption area. Model simulation results for trabecular bone volume fraction, microdamage, and activation frequency following 1 year of bisphosphonate treatment are consistent with experimental measurements. The model predicts that trabecular bone volume initially increases rapidly with 1 year of bisphosphonate treatment, and continues to slowly rise between 1 and 3 years of treatment. The model also predicts that microdamage initially increases rapidly, 0.5–1.5-fold for alendronate or risedronate during the first year of treatment, and reaches its maximum value by 2.5 years before trending downward for all dosages. The model developed in this study suggests that increasing bone volume fraction with long term bisphosphonate treatment may sufficiently reduce strain and damage formation rate so that microdamage does not accumulate above that which is initiated in the first two years of treatment.     

Microdamage propagation in trabecular bone due to changes in loading mode

Microdamage induced by falls or other abnormal loads that cause shear stress in trabecular bone could impair the mechanical properties of the proximal femur or spine. Existing microdamage may also increase the initiation and propagation of further microdamage during subsequent normal, on-axis, loading conditions, resulting in atraumatic or "spontaneous" fractures. Microdamage formation due to shear and compressive strains was studied in 14 on-axis cylindrical bovine tibial trabecular bone specimens. Microdamage was induced by a torsional overload followed by an on-axis compressive overload and quantified microscopically. Fluorescent agents were used to label microdamage and differentiate damage due to the two loading modes. Both the microcrack density and diffuse damage area caused by the torsional overload increased with increasing shear strain from the center to the edge of the specimen. However, the mean microcrack length was uniform across the specimen, suggesting that microcrack length is limited by microstructural features. The mean density of microcracks caused by compressive overloading was slightly higher near the center of the specimen, and the diffuse damage area was uniform across the specimen. Over 20% of the microcracks formed in the initial torsional overloading propagated during compression. Moreover the propagating microcracks were, on average, longer than microcracks formed by a single overload. As such, changes in loading mode can cause propagation of microcracks beyond the microstructural barriers that normally limit the length. Damage induced by in vivo off-axis loads such as falls may similarly propagate during subsequent normal loading, which could affect both remodeling activity and fracture susceptibility.     

Modeling modulus reduction in bovine trabecular bone damaged in compression.

Loading bone beyond its yield point creates microdamage, leading to reduction in stiffness. Previously, we related microdamage accumulation to changes in mechanical properties. Here, we develop a model that predicts stiffness loss based on the presence of microdamage. Modeling is done at three levels: (1) a single trabecula, (2) a cellular solid consisting of intact, damaged, and fractured trabeculae, and (3) a specimen with a localized damage band. Predictions of a reduced modulus agree well with experimental measured modulus reductions of post-yield compression of bovine trabecular bone. The predicted reduced modulus is relatively insensitive to changes in the input parameters.     

Effects of damage on the orthotropic material symmetry of bovine tibial trabecular bone

The macroscopic mechanical properties of trabecular bone can be predicted by its architecture using theoretical relationships between the elastic and architectural properties. Microdamage caused by overloading or fatigue decreases the apparent elastic moduli of trabecular bone requiring these relationships to be modified to predict the damaged elastic properties. In the case of isotropic damage, the apparent level elastic properties could be determined by multiplying all of the elastic constants by a single scalar factor. If the damage is anisotropic, the elastic constants may change by differing factors and the material coordinate system could become misaligned with the fabric coordinate system. High-resolution finite element models were used to simulate damage overloading on seven trabecular bone specimens subjected to pure shear strain in two planes. Comparison of the apparent elastic moduli of the specimens before and after damage showed that the reduction of the elastic moduli was anisotropic. This suggests that the microdamage within the specimens was inhomogeneous. However, after damage the specimens exhibited nearly orthotropic material symmetry as they did before damage. Changes in the orientation of the orthotropic material coordinate system were also small and occurred primarily in the transverse plane. Thus, while damage in trabecular bone is anisotropic, the material coordinate system remains aligned with the fabric tensor.     

Preparation of on-axis cylindrical trabecular bone specimens using micro-CT imaging

The Orientation of trabecular bone specimens for mechanical testing must be carefully controlled. A method for accurately preparing on-axis cylindrical specimens using high-resolution micro-CT imaging was developed. Sixteen cylindrical specimens were prepared from eight bovine tibiae. High-resolution finite element models were generated from micro-CT images of parallelepipeds and used to determine the principal material coordinate system of each parallelepiped. A cylindrical specimen was then machined with a diamond coring bit. The resulting specimens were scanned again to evaluate the orientation. The average deviation between the principal fabric orientation and the longitudinal axis of the cylindrical specimen was only 4.70 +/- 3.11 degrees.     

Simulation of vertebral trabecular bone loss using voxel finite element analysis

Trabecular bone loss in human vertebral bone is characterised by thinning and eventual perforation of the horizontal trabeculae. Concurrently, vertical trabeculae are completely lost with no histological evidence of significant thinning. Such bone loss results in deterioration in apparent modulus and strength of the trabecular core. In this study, a voxel-based finite element program was used to model bone loss in three specimens of human vertebral trabecular bone. Three sets of analyses were completed. In Set 1, strain adaptive resorption was modelled, whereby elements which were subject to the lowest mechanical stimulus (principal strain) were removed. In Set 2, both strain adaptive and microdamage mechanisms of bone resorption were included. Perforation of vertical trabeculae occurred due to microdamage resorption of elements with strains that exceeded a damage threshold. This resulted in collapse of the trabecular network under compression loading for two of the specimens tested. In Set 3, the damage threshold strain was gradually increased as bone loss progressed, resulting in reduced levels of microdamage resorption. This mechanism resulted in trabecular architectures in which vertical trabeculae had been perforated and which exhibited similar apparent modulus properties compared to experimental values reported in the literature. Our results indicate that strain adaptive remodelling alone does not explain the deterioration in mechanical properties that have been observed experimentally. Our results also support the hypothesis that horizontal trabeculae are lost principally by strain adaptive resorption, while vertical trabeculae may be lost due to perforation from microdamage resorption followed by rapid strain adaptive resorption of the remaining unloaded trabeculae.