Development of skeletal system for mesh-type ICRP reference adult phantoms

Development of Respiratory Tract Organs for ICRP Pediatric Mesh-type Reference Computational Phantoms

ABSTRACT

As part of the activities of the International Commission on Radiological Protection (ICRP) Task Group 103, the present study developed a new set of respiratory tract organs consisting of the extrathoracic, bronchial, bronchiolar, and alveolar-interstitial regions for newborn, 1-, 5-, 10-, and 15-y-old males and females for use in pediatric mesh-type reference computational phantoms. The developed respiratory tract organs, while preserving the original topologies of those of the pediatric voxel-type reference computational phantoms of ICRP Publication 143, have improved anatomy and detailed structure and also include μm-thick target and source regions prescribed in ICRP Publication 66. The dosimetric impact of the developed respiratory tract organs was investigated by calculating the specific absorbed fraction for internal electron exposures, which were then compared with the ICRP Task Group 96 values. The results showed that except for the alveolar-interstitial region as a source region, the pediatric mesh phantoms showed larger specific absorbed fractions than the Task Group 96 values. The maximum difference was a factor of ~3.5 for the extrathoracic-2 basal cell and surface as target and source regions, respectively. These results reflect the differences in the target masses and geometry caused by the anatomical enhancement of the pediatric mesh phantoms. For the alveolar-interstitial region as a source region, the pediatric mesh phantoms showed larger values for low energy ranges and lower values with increasing energies, owing to the differences in the size and shape of the alveolar-interstitial region.

Incorporation of micro-CT-based detailed bone models into ICRP adult mesh-type reference computational phantoms

Objective. The red bone marrow (RBM) and bone endosteum (BE), which are required for effective dose calculation, are macroscopically modeled in the reference phantoms of the international commission on radiological protection (ICRP) due to their microscopic and complex histology. In the present study, the detailed bone models were developed to simplify the dose calculation process for skeletal dosimetry. Approach. The detailed bone models were developed based on the bone models developed at the University of Florida. A new method was used to update the definition of BE region by storing the BE location indices using virtual sub-voxels. The detailed bone models were then installed in the spongiosa regions of the ICRP mesh-type reference computational phantoms (MRCPs) via the parallel geometry feature of the Geant4 code. Main results. Comparing the results between the detailed-bone-installed MRCPs and the original MRCPs with the absorbed dose to spongiosa and fluence-to-dose response function (DRF)-based methods, the DRF-based method showed much smaller but still significant differences. Compared with the values given in ICRP Publications 116 and 133, the differences were very large (i.e. several orders of magnitudes), due mainly to the anatomical improvement of the skeletal system in the MRCPs; that is, spongiosa and medullary cavity are fully enclosed by cortical bone in the MRCPs but not in the ICRP-110 phantoms. Significance. The detailed bone models enable the direct calculation of the absorbed doses to the RBM and BE, simplifying the dose calculation process and potentially improving the consistency and accuracy of skeletal dosimetry.

Development of alimentary tract organs for ICRP pediatric mesh-type reference computational phantoms

In line with the activities of Task Group 103 under the International Commission on Radiological Protection (ICRP), the present study was conducted to develop a new set of alimentary tract organs consisting of the oral cavity, oesophagus, stomach, small intestine, and colon for the newborn, 1 year-old, 5 year-old, 10 year-old, and 15 year-old males and females for use in the pediatric mesh-type reference computational phantoms (MRCPs). The developed alimentary tract organs of the pediatric MRCPs, while nearly preserving the original topology and shape of those of the pediatric voxel-type reference computational phantoms (VRCPs) of ICRPPublication 143, present considerable anatomical improvement and include all micrometre-scale target and source regions as prescribed in ICRPPublication 100. To investigate the dosimetric impact of the developed alimentary tract organs, organ doses and specific absorbed fractions were computed for certain external exposures to photons and electrons and internal exposures to electrons, respectively, which were then compared with the values computed using the current ICRP models (i.e. pediatric VRCPs and ICRP-100 stylised models). The results showed that for external exposures to penetrating radiations (i.e. photons >0.04 MeV), there was generally good agreement between the compared values, within a 10% difference, except for the oral mucosa. For external exposures to weakly penetrating radiations (i.e. low-energy photons and electrons), there were significant differences, up to a factor of ∼8300, owing to the geometric difference caused by the anatomical enhancement in the MRCPs. For internal exposures of electrons, there were significant differences, the maximum of which reached a factor of ∼73 000. This was attributed not only to the geometric difference but also to the target mass difference caused by the different luminal content mass and organ shape.

Development of skeletal systems for ICRP pediatric mesh-type reference computational phantoms

In 2016, the International Commission on Radiological Protection (ICRP) launched Task Group 103 (TG 103) for the explicit purpose of developing a new generation of adult and pediatric reference computational phantoms, named 'mesh-type reference computational phantoms (MRCPs)', that can overcome the limitations of voxel-type reference computational phantoms (VRCPs) of ICRPPublications 110and143due to their finite voxel resolutions and the nature of voxel geometry. After completing the development of the adult MRCPs, TG 103 has started the development of pediatric MRCPs comprising 10 phantoms (male and female versions of the reference newborn, 1-year-old, 5-year-old, 10-year-old, and 15-year-old). As part of the TG 103 project, within the present study, the skeletal systems, one of the most important and complex organ systems of the body, were developed for each phantom age and sex. The developed skeletal systems, while closely preserving the original bone topology of the pediatric VRCPs, present substantial improvements in the anatomy of complex and/or small bones. In order to investigate the dosimetric impact of the developed skeletons, the average absorbed doses and the specific absorbed fractions for radiosensitive skeletal tissues (i.e. active marrow and bone endosteum) were computed for some selected external and internal exposure cases, which were then compared with those calculated with the skeletons of pediatric VRCPs. The comparison result showed that the dose values of the pediatric MRCPs were generally similar to those of the pediatric VRCPs for highly penetrating radiations (e.g. photons >200 keV); however, for weakly penetrating radiations (e.g. photons ⩽200 keV and electrons), significant differences up to a factor of 140 were observed.

Effect of microdistribution of alpha and beta-emitters in targeted radionuclide therapies on delivered absorbed dose in a GATE model of bone marrow

Acute hematologic toxicity is a frequent adverse effect of beta-emitter targeted radionuclide therapies (TRTs). Alpha emitters have the potential of delivering high linear energy transfer (LET) radiation to the tumor attributed to its shorter range. Antibody-based TRTs have increased blood-pool half-lives, and therefore increased marrow toxicity, which is a particular concern with alpha emitters. Accurate 3D absorbed dose calculations focusing on the interface region of blood vessels and bone can elucidate energy deposition patterns. Firstly, a cylindrical geometry model with a central blood vessel embedded in the trabecular tissue was modeled. Monte Carlo simulations in GATE were performed considering beta (177Lu, 90Y) and alpha emitters (211At, 225Ac) as sources restricted to the blood pool. Subsequently, the radioactive sources were added in the trabecular bone compartment in order to model bone marrow metastases infiltration (BMMI). Radial profiles, dose-volume histograms and voxel relative differences were used to evaluate the absorbed dose results. We demonstrated that alpha emitters have a higher localized energy deposition compared to beta emitters. In the cylindrical geometry model, when the sources are confined to the blood pool, the dose to the trabecular bone is greater for beta emitting radionuclides, as alpha emitters deposit the majority of their energy within 70 μm of the vessel wall. In the BMMI model, alpha emitters have a lower dose to untargeted trabecular bone. Our results suggest that when alpha emitters are restricted to the blood pool, as when labeled to antibodies, hematologic toxicities may be lower than expected due to differences in the microdistribution of delivered absorbed dose.

Dosimetric considerations of 99mTc-MDP uptake within the epiphyseal plates of the long bones of pediatric patients

Skeletal scintigraphy is most performed in pediatric patients using the radiopharmaceutical ^99mTc labelled methylene diphosphonate (^99mTc-MDP). Reference biokinetic models for ^99mTc-MDP indicate 50% of the administered activity is uniformly localized to the interior bone surfaces (trabecular and cortical regions), yet imaging data clearly show some preferential uptake to the epiphyseal growth plates of the long bones. To explore the dosimetric consequences of these regional activity concentrations, we have modified mesh-type computational phantoms of the International Commission on Radiological Protection (ICRP) reference pediatric series to explicitly include geometric models of the epiphyseal growth plates (2 mm in thickness) within the left/right, distal/proximal ends of the humeri, radii, ulnae, femora, tibia, and fibulae. Bone mineral activity from the ICRP Publication 128 biokinetic model for ^99mTc-MDP (ICRP 2015) was then partitioned to the growth plates at values of 0.5%, 4.4%, 8.3%, 12.2%, 16.1%, and 20%. Radiation transport simulations were performed to compute ^99mTc S-values and organ dose coefficients to the soft tissues and to bone site-specific regions of spongiosa. As the percentage of bone activity assigned to the growth plates was increased (from 0.5% to 20%), absorbed doses to the soft tissue organs, active bone marrow, bone endosteum (BE), as well as the detriment-weighted dose, were shown to decrease from their nominal values (no substantial growth plate activity), while epiphyseal plate self-doses increased. In the 15 year old male phantom, moving from 0.5% to 20% relative bone activity within the epiphyseal plates resulted in a 15% reduction in active marrow (AM) and BE dose, a 10% reduction in mean soft tissue and detriment-weighted dose, and a 6.3-fold increase in epiphyseal plate self-dose. In the newborn female phantom, we observed a 18% decrease in AM and BE dose, a 10% decrease in mean soft tissue dose, a 15% decrease in detriment-weighted dose, and 12.8-fold increase in epiphyseal plate self-dose. Increases (to 3 mm) and decreases (to 1 mm) in the assumed growth plate thickness of our models were shown to impact only the growth plate self-dose. Future work in differential quantification of ^99mTc-MDP activity-growth plates versus other bone surfaces-is required to provide clinically realistic data on activity partitioning as a function of patient age, and perhaps skeletal site. The phantom series presented here may be used to develop more optimized age-related guidance on ^99mTc-MDP administered activities to children.

INVESTIGATION OF THE INFLUENCE OF THYROID LOCATION ON IODINE-131 S VALUES

The use of iodine-131 S values based on reference computational phantoms with fixed thyroid model may lead to significant dosimetric errors in patients who may have different thyroid location from the reference phantoms. In the present study, we investigated individual thyroid location variation by examining the computed tomography image sets of 40 adult male and female patients. Subsequently, the thyroid location of the adult male and female mesh-type reference phantoms of the International Commission on Radiological Protection (ICRP) was adjusted to match each the highest, mean and the lowest locations of the thyroid observed in this dataset. The thyroid-adjusted phantoms were implemented into the Geant4 Monte Carlo code to calculate thyroid location-dependent iodine-131 S values (rT ← thyroid) for a total of 30 target regions. The maximum variation among the observed thyroid locations was 39 mm and 36 mm for male and female patients, respectively. The mean thyroid locations of both male and female patients showed a good agreement with the ICRP reference phantoms. The thyroid location-dependent Iodine-131 S values were significantly different from the reference phantoms for most target regions by up to a factor of 3. The use of thyroid location-dependent S values in dose reconstructions should help quantify the dosimetric uncertainty in epidemiologic investigations of patients receiving iodine-131 therapy for hyperthyroidism and thyroid cancer.

DOSE COEFFICIENTS OF MESH-TYPE ICRP REFERENCE COMPUTATIONAL PHANTOMS FOR EXTERNAL EXPOSURES OF NEUTRONS, PROTONS, AND HELIUM IONS

To overcome inherent limitations of the Voxel-type Reference Computational Phantoms (VRCPs) due to the limited voxel resolutions and the nature of voxel geometry, the International Commission on Radiological Protection (ICRP) has developed the adult male and female Mesh-type Reference Computational Phantoms (MRCPs). We previously used the MRCPs to calculate a complete set of dose coefficients (DCs) for idealized external exposures of photons and electrons (Yeom et al. NET in press). In the present study, we extended the previous study to include additional radiation particles (neutrons, protons, and helium ions) into the DC library by conducing Monte Carlo radiation transport simulations with the Geant4 code. The MRPC-based DCs were compared with the existing reference DCs of ICRP Publication 116 which are based on the ICRP VRCPs to investigate impact of the new mesh-type reference phantoms on the DC values. We found that the MRCPs generally provide DCs of organ/tissue doses and effective doses similar to those from the VRCPs for penetrating radiations (uncharged particles), whereas significant DC differences were observed for weakly penetrating radiations (charged particles) mainly due to the improved representation of the detailed anatomical structures in the MRCPs over the VRCPs.

Body-size-dependent phantom library constructed from ICRP mesh-type reference computational phantoms

Recently, ICRP Task Group 103 developed new mesh-type reference computational phantoms (MRCPs) for the adult male and female by converting the current voxel-type reference computational phantoms (VRCPs) of ICRP Publication 110 into a high-quality/fidelity mesh format. Utilizing the great deformability/flexibility of the MRCPs compared with the VRCPs, in the present study, we established a body-size-dependent phantom library by modifying the MRCPs. The established library includes 108 adult male and 104 adult female phantoms in different standing heights and body weights, covering most body sizes representative of Caucasian and Asian populations. Ten secondary anthropometric parameters with respect to standing height and body weight were derived from various anthropometric databases and applied in the construction of the phantom library. An in-house program for automatic phantom adjustment was developed and applied for practical construction of such a large number of phantoms in the library with minimized human intervention. Organ/tissue doses calculated with three male phantoms in different standing heights (165, 175, and 190 cm) with a fixed body weight of 80 kg for external exposures to broad parallel photon beams from 0.01 to 10⁴ MeV were compared, observing there are significant dose differences particularly for the photon energies <0.1 MeV in which the organ/tissue doses tended to increase with increasing standing height. In addition, the organ/tissue doses of three female phantoms in different body weights (45, 65, and 140 kg) with a fixed standing height of 165 cm were compared, showing a significant decreasing tendency with increasing body weight for the photon energies <10 MeV. For the higher energies, the opposite trend, interestingly, was observed; that is, the organ/tissue doses tended to increase with increasing body weight. The results, despite the limited number of exposure cases, suggest that the use of the body-size-dependent phantom library can improve the accuracy of individual dose estimates for many retrospective dosimetry studies by taking the body size of individuals into account.

Dose coefficients of percentile-specific computational phantoms for photon external exposures

The use of dose coefficients (DCs) based on the reference phantoms recommended by the International Commission on Radiological Protection (ICRP) with a fixed body size may produce errors to the estimated organ/tissue doses to be used, for example, for epidemiologic studies depending on the body size of cohort members. A set of percentile-specific computational phantoms that represent 10th, 50th, and 90th percentile standing heights and body masses in adult male and female Caucasian populations were recently developed by modifying the mesh-type ICRP reference computational phantoms (MRCPs). In the present study, these percentile-specific phantoms were used to calculate a comprehensive dataset of body-size-dependent DCs for photon external exposures by performing Monte Carlo dose calculations with the Geant4 code. The dataset includes the DCs of absorbed doses for 29 individual organs/tissues from 0.01 to 104 MeV photon energy, in the antero-posterior, postero-anterior, right lateral, left lateral, rotational, and isotropic geometries. The body-size-dependent DCs were compared with the DCs of the MRCPs in the reference body size, showing that the DCs of the MRCPs are generally similar to those of the 50th percentile standing height and body mass phantoms over the entire photon energy region except for low energies (≤ 0.03 MeV); the differences are mostly less than 10%. In contrast, there are significant differences in the DCs between the MRCPs and the 10th and 90th percentile standing height and body mass phantoms (i.e., H10M10 and H90M90). At energies of less than about 10 MeV, the MRCPs tended to under- and over-estimate the organ/tissue doses of the H10M10 and H90M90 phantoms, respectively. This tendency was revised at higher energies. The DCs of the percentile-specific phantoms were also compared with the previously published values of another phantom sets with similar body sizes, showing significant differences particularly at energies below about 0.1 MeV, which is mainly due to the different locations and depths of organs/tissues between the different phantom libraries. The DCs established in the present study should be useful to improve the dosimetric accuracy in the reconstructions of organ/tissue doses for individuals in risk assessment for epidemiologic investigations taking body sizes into account.

Dosimetric impact of voxel resolutions of computational human phantoms for external photon exposure

Several research teams have developed computational phantoms in polygonal-mesh (PM) and/or Non-Uniform Rational B-Spline format, but it has not been systematically evaluated if the existing voxel phantoms are still dosimetrically valid. We created three voxel phantoms with the resolutions of 1,000, 125, and 1 mm3 and simulated the irradiation in antero-posterior geometry with photons of 0.1, 1, and 10 MeV using voxel Monte Carlo codes, and compared the energy deposition to their organs/tissues with the values from the original PM phantom using mesh Monte Carlo codes. The coefficient of variation in energy deposition overall showed about five-fold decrease as the voxel resolution increased but differences were mostly less than 5% for any voxel resolution. We conclude that PM phantoms and mesh Monte Carlo techniques may not be necessary for external photon exposure (0.1 - 10 MeV) and the existing voxel phantoms can provide enough dosimetric accuracy in those exposure conditions.

Mesh-type reference Korean phantoms (MRKPs) for adult male and female for use in radiation protection dosimetry

In the present study, to overcome the dosimetric limitations of the previous voxel-type reference Korean computational phantoms due to their limited voxel resolutions (i.e. on the order of millimeters) and the nature of voxel geometry, a pair of new reference Korean phantoms, called mesh-type reference Korean phantoms (MRKPs), were developed for the adult male and female in a high-quality/fidelity mesh format. The developed phantoms include all target and source regions required for effective dose calculation, even micrometer-scale target and source regions of the respiratory and alimentary tract organs, skin, urinary bladder, and eye lens. The developed phantoms, which are in either the polygon-mesh (PM) format or the tetrahedral-mesh (TM) format as necessary, can be directly used in several general-purpose Monte Carlo codes (e.g. Geant4, MCNP6, and PHITS) without voxelization. In order to understand the dosimetric impact of the new phantoms, the dose coefficients (=fluence-to-effective dose conversion coefficients) were calculated for photons and electrons with energies ranging from 10 keV to 10 GeV for the anterior-posterior (AP) irradiation geometry and compared with those of the previous voxel-type reference Korean phantoms. The results demonstrate that the effective dose coefficients of the MRKPs were generally similar to those of the previous voxel-type reference phantoms for photons; however, for electrons, significant differences were observed at energies lower than 1 MeV that were mainly due to the explicit definition of the 50 µm-thick radiosensitive target layer in the skin of the new mesh phantoms.

New mesh-type phantoms and their dosimetric applications, including emergencies

Committee 2 of the International Commission on Radiological Protection (ICRP) has constructed mesh-type adult reference computational phantoms by converting the voxel-type ICRP Publication 110 adult reference computational phantoms to a high-quality mesh format, and adding those tissues that were below the image resolution of the voxel phantoms and therefore not included in the Publication 110 phantoms. The new mesh phantoms include all the necessary source and target tissues for effective dose calculations, including the 8-40-µm-thick target layers of the alimentary and respiratory tract organs, thereby obviating the need for supplemental organ-specific stylised models (e.g. respiratory airways, alimentary tract organ walls and stem cell layers, lens of the eye, and skin basal layer). To see the impact of the new mesh-type reference phantoms, dose coefficients for some selected external and internal exposures were calculated and compared with the current reference values in ICRP Publications 116 and 133, which were calculated by employing the Publication 110 phantoms and the supplemental stylised models. The new mesh phantoms were also used to calculate dose coefficients for industrial radiography sources near the body, which can be used to estimate the organ doses of the worker who is accidentally exposed by an industrial radiography source; in these calculations, the mesh phantoms were deformed to reflect the size of the worker, and also to evaluate the effect of posture on dose coefficients.

Computational phantoms, ICRP/ICRU, and further developments

Phantoms simulating the human body play a central role in radiation dosimetry. The first computational body phantoms were based upon mathematical expressions describing idealised body organs. With the advent of more powerful computers in the 1980s, voxel phantoms have been developed. Being based on three-dimensional images of individuals, they offer a more realistic anatomy. Hence, the International Commission on Radiological Protection (ICRP) decided to construct voxel phantoms representative of the adult Reference Male and Reference Female for the update of organ dose coefficients. Further work on phantom development has focused on phantoms that combine the realism of patient-based voxel phantoms with the flexibility of mathematical phantoms, so-called 'boundary representation' (BREP) phantoms. This phantom type has been chosen for the ICRP family of paediatric reference phantoms. Due to the limited voxel resolution of the adult reference computational phantoms, smaller tissues, such as the lens of the eye, skin, and micron-thick target tissues in the respiratory and alimentary tract regions, could not be segmented properly. In this context, ICRP Committee 2 initiated a research project with the goal of producing replicas of the ICRP Publication 110 phantoms in polygon mesh format, including all source and target regions, even those with micron resolution. BREP phantoms of the fetus and the pregnant female at various stages of gestation complete the phantoms available for radiation protection computations.

Individualized adjustments to reference phantom internal organ dosimetry-scaling factors given knowledge of patient internal anatomy

Various computational tools are currently available that facilitate patient organ dosimetry in diagnostic nuclear medicine, yet they are typically restricted to reporting organ doses to ICRP-defined reference phantoms. The present study, while remaining computational phantom based, provides straightforward tools to adjust reference phantom organ dose for both internal photon and electron sources. A wide variety of monoenergetic specific absorbed fractions were computed using radiation transport simulations for tissue spheres of varying size and separation distance. Scaling methods were then constructed for both photon and electron self-dose and cross-dose, with data validation provided from patient-specific voxel phantom simulations, as well as via comparison to the scaling methodology given in MIRD Pamphlet No. 11. Photon and electron self-dose was found to be dependent on both radiation energy and sphere size. Photon cross-dose was found to be mostly independent of sphere size. Electron cross-dose was found to be dependent on sphere size when the spheres were in close proximity, owing to differences in electron range. The validation studies showed that this dataset was more effective than the MIRD 11 method at predicting patient-specific photon doses for at both high and low energies, but gave similar results at photon energies between 100 keV and 1 MeV. The MIRD 11 method for electron self-dose scaling was accurate for lower energies but began to break down at higher energies. The photon cross-dose scaling methodology developed in this study showed gains in accuracy of up to 9% for actual patient studies, and the electron cross-dose scaling methodology showed gains in accuracy up to 9% as well when only the bremsstrahlung component of the cross-dose was scaled. These dose scaling methods are readily available for incorporation into internal dosimetry software for diagnostic phantom-based organ dosimetry.