SEARCH · Engineering Papers
Results for “Logs”
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Litter production and foliar nutrient resorption in fast- and slow-growing tree species in the Central Amazon
Litterfall is crucial for forest maintenance, serving as a primary mechanism for nutrient return to the nutrient-poor soils of tropical forests. Foliar nutrient resorption likewise represents an important nutrient-conservation mechanism. Yet, little is known about how these processes vary between fast- and slow-growing species in post-logging areas of the Amazon forest. Here, the objective of this study was to quantify litterfall production and the resorption of foliar nutrients in fast-growing and slow-growing tree species of the Central Amazon, in a forest that was experimentally logged in 1987. The study was conducted from May 2022 to April 2023. Litterfall was collected biweekly using four collectors that were systematically distributed beneath the canopy of each monitored tree, totaling 72 collectors. Three fast-growing and three slow-growing species were selected, each with three replicates, totaling 18 monitored individuals. Species-specific samples of fresh (green) and senesced (litter) leaves were collected and analyzed for their nutrient content and resorption efficiency. Fast-growing species had a monthly leaf litter deposition of 13.53 ± 1.6 g m −2 month −1 , compared to 2.59 ± 0.4 g m −2 month −1 for slow-growing species. The average annual litter production across both functional types was 8.6 ± 2.6 Mg ha −1 year −1 . Nutrient inputs through litterfall were higher in fast-growing species for all elements, particularly nitrogen (N), with 21.92 ± 4.9 kg ha −1 year −1 . Phosphorus (P) and potassium (K) exhibited the highest foliar resorption. P resorption efficiency was 68.3 % in fast-growing species and 57.8 % in slow-growing species. For K, efficiencies were 59.0 % and 41.7 %, respectively. These results highlight the substantial role that fast-growing species play in restoring forest productivity in managed Amazon forests, both through higher litter deposition and nutrient fluxes, and through nutrient conserving-mechanisms such as foliar nutrient resorption.
Reductive Dissolution of Goethite and Hematite by Reduced Flavins
The abiotic reductive dissolution of goethite and hematite by the reduced forms of flavin mononucleotide (FMNH2) and riboflavin (RBFH2), electron transfer mediators (ETM) secreted by the dissimilatory iron-reducing bacterium Shewanella, was investigated under stringent anaerobic conditions. In contrast to the rapid redox reaction rate observed for ferrihydrite and lepidocrocite (Shi et al., 2012), the reductive dissolution of crystalline goethite and hematite was slower, with the extent of reaction limited by the thermodynamic driving force at circumneutral pH. Both the initial reaction rate and reaction extent increased with decreasing pH. On a unit surface area basis, goethite was less reactive than hematite between pH 4.0 and 7.0. AH2DS, the reduced form of the well-studied synthetic ETM anthraquinone-2,6-disulfonate (AQDS), yielded higher rates than FMNH2 under most reaction conditions, despite the fact that FMNH2 was a more effective reductant than AH2DS for ferryhydrite and lepidocrocite. Two additional model compounds, methyl viologen and benzyl viologen, were investigated under similar reaction conditions to explore the relationship between reaction rate and thermodynamic properties. Relevant kinetic data from the literature were also included in the analysis to span a broad range of half-cell potentials. Other conditions being equal, the surface area normalized initial reaction rate (ra) increased as the redox potential of the reductant became more negative. A non-linear, parabolic relationship was observed between log ra and the redox potential for eight reducants at pH 7.0, as predicted by Marcus theory for electron transfer. When pH and reductant concentration were fixed, log ra was positively correlated to the redox potential of four Fe(III) oxides over a wide pH range, following a non-linear parabolic relationship as well.
Local structural variation with oxygen fugacity in Fe 2 SiO 4+x fayalitic iron silicate melts
Here, the structure of molten Fe 2 SiO 4+x has been studied using both high-energy X-ray diffraction and Fe K-edge X-ray absorption near-edge structure (XANES) spectroscopy, combined with aerodynamic levitation and laser beam heating. A wide range of Fe 3+ contents were accessed by varying the levitation and atmospheric gas composition. Diffraction measurements were made in the temperature ( T ) and oxygen partial pressure ranges 1624(21) < T < 2183(94) K (uncertainties in parentheses) and –5.6(3) < ΔFMQ < +2.8(5) log units (relative to the Fayalite-Magnetite-Quartz buffer). Iron K-edge XANES measurements covered the ranges 1557(33) < T < 1994(36) K and –2.1(3) < ΔFMQ < +4.4(3) log units. Fe 3+ contents, x = Fe 3+ /ΣFe, estimated directly from the pre-edge peaks of the XANES spectra varied between 0.15(1) and 0.40(2). While these agree in some cases with semi-empirical models, notable discrepancies are discussed in the context of the redox kinetics and the limitations in both the models and in the calibrations used to derive oxidation state from XANES spectra. XANES pre-edge peak areas imply average Fe–O coordination numbers, n FeO , close to 5 for all Fe 3+ /ΣFe. Diffraction measurements yielded values of 4.4(2) < n FeO < 4.7(1). There is limited evidence for a linear trend n FeO (x) = 4.46(3) + 0.4(1)x. Asymmetric Fe–O bond length distributions peak at around 1.96 Å and have a shoulder arising from longer interatomic distances. Mean r FeO lie close to 2.06 Å, consistent with n FeO close to 5. These observations suggest that Fe 2+ is less efficient at stabilizing tetrahedral Fe 3+ compared to large monovalent alkali cations. Comparison of in-situ XANES estimates of Fe 3+ /ΣFe in the melts to those of the quenched solids obtained from XANES as well as Mössbauer spectroscopy indicate rapid oxidation during cooling, enabled by stirring of the melt by the levitation gas flow. As such, the oxidation state of hot komatiitic and other highly fluid melts may not be retained, even during rapid cooling, as it is for cooler basaltic and more silicic magmas.
Rates and mechanisms of uranyl oxyhydroxide mineral dissolution
Uranyl oxyhydroxide minerals are important weathering products in uranium-contaminated surface and subsurface environments that regulate dissolved uranium concentrations. However, dissolution rates for this class of minerals and associated dissolution mechanisms have not been previously reported for circumneutral pH conditions, particularly for the case of flow through porous media. In this paper, the dissolution rates of K- and Na-compreignacite (K 2 (UO 2 ) 6 O 4 (OH) 6 ·8H 2 O and Na 2 (UO 2 ) 6 O 4 (OH) 6 ·8H 2 O respectively) were measured using flow-through columns reacted with two simulated background porewater (BPW) solutions of low and high dissolved total carbonate content (ca. 0.2 and 2.8 mmol L -1 ). Column materials were characterized before and after reaction with electron microscopy, bulk chemistry, and EXAFS to identify structural and chemical changes during dissolution and to obtain insight into molecular-scale processes. The reactive transport code CrunchFlow was used to calculate overall dissolution rates while accounting for fluid transport and changes in mineral volume and reactive surface area and results were compared to steady-state dissolution rate calculations. In low carbonate BPW systems, interlayer K and Na were initially leached from both minerals, and in Na-compreignacite, K and minor divalent cations from the input solution were incorporated into the mineral structure. Results of characterization analyses suggested that after reaction both K- and Na-compreignacite resembled a disordered K-compreignacite with altered surfaces. A 10-fold increase in dissolved carbonate concentration and corresponding increase in pH (from 6.65 to 8.40) resulted in a net removal of 58-87% of total uranium mass from the columns, compared to <1% net loss in low carbonate BPW systems. Steady-state release of dissolved uranium was not observed with high carbonate solutions and post-reaction characterizations indicated a lack of development of leached or altered surfaces. Dissolution rates (normalized to specific surface area) were about 2.5-3 orders-of-magnitude faster in high versus low carbonate BPW systems, with Na-compreignacite dissolving more rapidly than K-compreignacite under both BPW conditions, possibly due to greater ion exchange (1.57·10 -10 vs. 1.28·10 -13 mol m -2 s -1 [log R = -9.81 and -12.89] and 5.79·10 -10 vs. 3.71·10 -13 mol m -2 s -1 [log R = -9.24 and -12.43] for K- and Na-compreignacite respectively). Experimental and spectroscopic results suggest that the dissolution rate is controlled by bond breaking of a uranyl group and detachment from polyhedral layers of the mineral structure. With higher dissolved carbonate concentrations, this rate-determining step is accelerated by the formation of Ca-uranyl carbonate complexes (dominant species under these conditions), which resulted in an increase of the dissolution rates. Optimization of both dissolution rate and mineral volume fraction in the reactive transport model to account for uranium mass removal during dissolution more accurately reproduced effluent data in high carbonate systems, and resulted in faster overall rates compared with a steady-state dissolution assumption. Finally, this study highlights the importance of coupling reaction and transport processes during the quantification of mineral dissolution rates to accurately predict the fate of contaminants such as uranium in porous geomedia.
The solubility of thorium in carbonate-bearing solutions at hydrothermal conditions
Thorium mineralization is frequently hosted in carbonate-bearing rocks, and thorium commonly substitutes into the structures of carbonate-bearing minerals that have precipitated from or been modified by hydrothermal fluids. Given this common association, it is reasonable to consider the hypothesis that the presence of carbonate ligands in hydrothermal solutions promotes the transport of Th through the formation of stable aqueous complexes. Our ability to evaluate this hypothesis, however, is hindered by the lack of experimental data for Th-carbonate species at conditions beyond ambient. The low-temperature data indicate that carbonate is a strong complexing agent for Th. In this contribution, we investigate the solubility of Th in carbonate-bearing fluids relevant to natural systems (0.05–0.5 m NaHCO 3 /Na 2 CO 3 ; pH T ~ 7.8–9.8) at elevated temperature (175–250 °C). We demonstrate that, in contrast to the behavior of Th at low temperature, the stability of Th-carbonate complexes is not sufficient for them to predominate at these conditions. Instead, the solubility of Th is governed by hydrolysis reactions. Under the experimental conditions investigated, the predominant hydroxyl complexes are Th(OH) 4 0 and Th(OH) 5 -. Furthermore, thermodynamic formation constants were derived for these species at the temperatures considered in our experiments (log β 4 = 43.34 and 44.31 at 175 and 200 °C, respectively, and log β 5 = 46.15 and 47.9 at 225 and 250 °C, respectively) to permit forward modeling of Th mobility in natural systems. Our study indicates that carbonate ions are unlikely to play a role in transporting Th in hydrothermal fluids. Summarizing the results of this study and our previous studies of the solubility of Th in hydrothermal fluids, we conclude that SO 4 2- is the primary ligand responsible for the hydrothermal transport of Th.
New insights into Mn 2+ and Mg 2+ inhibition of calcite growth
Impurity ion and isotope partitioning into carbonate minerals provide a window into the molecular processes occurring at the fluid-mineral interface during crystal growth. Here, we employ calcium isotope fractionation together with process-based modeling to elucidate the mechanisms by which two divalent cations with starkly contrasting compatibility, magnesium and manganese, inhibit calcite growth and incorporate into the mineral lattice. Calcite growth inhibition by Mg 2+ is log-linear and K Mg is on the order of 0.02–0.03 throughout the range of {Mg 2+ }/{Ca 2+ } studied here (0.01–2.6). Mn 2+ exhibits much stronger log-linear growth rate inhibition at low Mn 2+ concentrations (fluid {Mn 2+ }/{Ca 2+ } = 0.001–0.02). Mn 2+ is readily incorporated into the calcite lattice to form a calcite-rhodochrosite solid solution, with large partition coefficients (K Mn 4.6–15.6) inversely correlated to growth rate. For both Mn 2+ and Mg 2+ , calcium isotope fractionation is found to be invariant with {Me 2+ }/{Ca 2+ } despite more than an order of magnitude decline in growth rate. This invariant Δ 44/40 Ca suggests that the presence of Mn 2 + or Mg 2+ does not significantly change the relative rates of Ca 2+ attachment and detachment at kink sites during growth, indicative of a dominantly kink blocking inhibition mechanism. Because the partitioning behavior dictates that Mn 2+ must attach to the surface significantly faster than Ca 2+ , attachment of Mn 2+ is likely to be as a non-monomer species such as an ion pair or possibly a larger polynuclear cluster. We propose that calcite growth rate inhibition by Mn is determined by the kinetics of carbonate attachment at Mn-occupied kink sites, potentially due to slow re-orientation kinetics of carbonate ions that have formed an inner-sphere complex with Mn 2+ at the surface but must reorient to incorporate into the lattice. We demonstrate that patterns in Mg 2+ partitioning and inhibition behavior are broadly consistent with growth inhibition driven by slow Mg 2+ -aquo complex dehydration relative to Ca 2+ but argue that this mechanism likely represents one endmember scenario, seen in Mg-calcite growth at low supersaturations and net precipitation rates. During growth at faster net precipitation rates, some portion of Mg 2+ is likely incorporated as a partially hydrated or otherwise complexed species, but calcite growth remains significantly inhibited by the kinetics of CO 3 2- attachment at Mg 2+ kink sites. These findings suggest a hybrid classical/nonclassical growth mechanism whereby Ca 2+ incorporates largely as a free ion at kink sites while Mn 2+ and some portion of Mg 2+ are incorporated via non-monomer attachment. This pattern may be generalizable; trace constituent cations with aquo-complex desolvation rates significantly slower than the mineral growth rate preferentially incorporate as a non-monomer species during otherwise classical crystal growth.
Experimental study of the partitioning of some platinum group elements (Pd and Ir) between orthopyroxene and silicate melt
Past experiments and observations on natural samples have largely focused on the roles of olivine and chromite in controlling the behaviour of the platinum-group elements (PGE) during melting and solidification, whereas other phases, such as pyroxene, have gone largely uncharacterized. Here, to address this, experiments have been done to measure the partitioning of Pd (with a subset of results for Ir), between orthopyroxene and silicate melt at 1340 °C, 0.1 MPa and log fO 2 of FMQ - 1 to FMQ + 6 (FMQ = Fayalite-Magnetite-Quartz). The X-ray Absorption Near-Edge Structure (XANES) was measured in a subset of experiment glasses. Glass concentrations of Pd (corrected to unit Pd activity) increase from ~6 to ~650 ug/g over the fO 2 range of experiments. The slope of the solubility-fO 2 relation is consistent with Pd 1+ as the dominant oxidation state, with evidence for Pd 0 and Pd 2+ at the lowest and highest experiment fO 2 , respectively. Consistent with this result, the XANES reveal spectral features similar to Pd 0 and Pd 2+ spectral reference materials (specRM) at the most reduced and oxidized synthesis conditions, respectively. Other lines of evidence require the presence of a third melt species, here interpreted to be Pd 1+ . Values of orthopyroxene/melt partition coefficients for Pd (D Pd Opx/melt ) are 0.0051 (+/-0.006) at log fO 2 < ΔFMQ + 3, increasing with fO 2 to a maximum of 0.013 at ~FMQ + 6. Sodium partition coefficients, expected to be similar to Pd, range from 0.0061 (+/-0.00061) at FMQ + 3, increase to 0.007–0.009 at higher fO 2 , but with no clear systematic trend. A value for D Ir Opx-melt of ~0.6 was measured at ~FMQ + 4, indicating significantly more compatible behaviour for Ir relative to Pd. Partitioning results are interpreted in the context of the Blundy-Wood elastic strain model in which the variation in partitioning is related to ionic radius mismatch to an optimal crystallographic site size. Based on the trend in ionic radius with oxidation state, the estimated ionic radius of Pd 1+ in octahedral coordination is similar to Na 1+ , and comparison to previous orthopyroxene-melt partitioning experiments suggests D Pd1+ opx/melt and DNaopx/melt should be nearly identical, consistent with the results of this study. The ionic radius of VI-fold Pd 2+ is close to the optimal M2 site size, so an increased proportion of this species with fO 2 accounts for the larger values of D Pd opx/melt at the highest fO 2 investigated. The much larger partition coefficient for Ir is consistent with the presence of Ir 2+ , whose estimated ionic radius is close to Fe 2+ and Mg 2+ , as well as predictions for the optimal M1 site size. With the assumption that D Pd opx/melt = D Na opx/melt , combined with a revised value for the Pd content of the primitive mantle, a melting model is presented that better reproduces the Pd concentration of high degree melts from sulfide-free mantle sources.
The effect of composition, temperature and pressure on the oxidation state and coordination environment of copper in silicate melts
Copper is a redox variable element that may occur as Cu 0 , Cu + , and Cu 2+ in the Earth’s crust. The oxidation state will affect its partitioning between coexisting minerals, melts and fluids and hence its behaviour in magmatic processes. Copper bearing silicate glasses were quenched from melts with 23 synthetic compositions (19 CaO-Na 2 O-MgO-Al 2 O 3 -SiO 2 (CNMAS), two “granites” containing K 2 O ± H 2 O, and Fe-bearing “MORB” and “andesite”) equilibrated at oxygen fugacities (fO 2 ), expressed in log units relative to the fayalite-magnetite-quartz (FMQ) buffer, ranging from −0.7 to 14, temperatures from 900 to 1500 °C and pressures from 0 to 2.5 GPa. Cu K-edge X-ray absorption near edge structure (XANES) spectra were recorded from the glasses and a pre-edge feature in the XANES spectra was found to scale with the proportion of Cu + . Cu + /ΣCu (where ΣCu = Cu + + Cu 2+ ) was quantified by fitting the intensity of the pre-edge feature as a function of fO 2 to the thermodynamically expected relationship. Cu + /ΣCu was found to only weakly depend on melt composition, with more basic melts (e.g., basalts rather than granites) preferentially stabilising Cu 2+ . Increasing temperature stabilises Cu + , while increasing pressure had little effect on Cu + /ΣCu in CNMAS melts but preferentially stabilised Cu 2+ in granite melts. Cu + /ΣCu can be predicted in silicate melts by the empirical equation: log(Cu 2+ /Cu + ) = 0.25(ΔFMQ + 8.58 − 25050/T + 940P/T – 0.02P) − 4.73 + 5400/T + 1.99Λ + (280P – 90P2) where T is temperature in K, P is pressure in GPa and is the optical basicity of the composition. The effects of fO 2 , melt composition, temperature and pressure on Cu + /ΣCu indicate that Cu + will be the dominant oxidation state in terrestrial silicate melts (e.g., Cu + /ΣCu = 99% in an andesitic melt at 900 °C, 1 GPa and ΔFMQ = 1). The electron exchange reaction Cu 2+ + Fe 2+ → Cu + + Fe 3+ occurs on cooling and, given the abundance of Fe in natural melts, the oxidation state of Cu in natural glasses is unlikely to correspond to that in the original melt.
The solubility of ErPO 4 and Er speciation in hydrothermal fluids at varying pH and salinity between 350 and 450 °C
The rare earth elements (REE) are important for the green-energy transition and can be incorporated into the REE phosphates, such as xenotime-(Y), which also hosts heavy REE (Tb– Lu). Xenotime-(Y) is a common accessory mineral in metamorphic rocks and a range of mineral deposits where it controls the mobility of heavy REE, however, the impact of high temperature aqueous fluids on the behavior of heavy REE is largely unknown. Thermodynamic modeling can be utilized as a tool to predict the mobility of REE in hydrothermal aqueous fluids, but must be supported by accurate experimental data. Here, we measured the solubility of endmember synthetic xenotime-structured ErPO4 in NaCl-HCl-NaOH-bearing aqueous solutions at 350 °C and water vapor saturation pressure, at 400 and at 450 °C and 500 bar using batch-type Inconel reactors. Erbium speciation was investigated as a function of pH from 2.8 to 8, where Er chloride species are predominant at acidic conditions (pH <3) and Er hydroxyl complexes are predominant at near- neutral to alkaline conditions (pH >3). At pH 7–9, the measured ErPO 4 solubility (-9.8 to -7.5 log m Er ) is up to 2.5 orders of magnitude lower than thermodynamic predictions (-9.4 to -6.7 log m Er ) using existing thermodynamic databases. At pH 2–3, the predicted ErPO 4 solubility is ~0.5 orders of magnitude higher at 350 °C and ~1 order of magnitude lower at 450 °C compared to experimentally measured Er concentrations. The thermodynamic properties of aqueous Er species were therefore revised in this study. The partial molal Gibbs energy of formation (Δ f G 0 T,P ) for aqueous Er hydroxyl and chloride species are optimized using GEMSFITS and the logarithmic formation constants (logβ n ( Cl,OH) ) were derived at each experimental temperature and pressure. The updated thermodynamic properties for Er hydroxyl species (Er(OH) +2 , Er(OH) 2 + , and Er(OH) 3 0 ) show that their stability shifts to more acidic conditions at and below 400 °C. The Er chloride species (ErCl +2 and ErCl 2 + ) show increased stability compared to Er hydroxyl species at temperatures of 450 °C and 0.01 mol/kg NaCl. The updated thermodynamic properties are implemented into the GEM-Selektor modeling package to investigate the mobility of Er in saline hydrothermal fluids in equilibrium with alkaline rocks. Importantly, the updated properties for Er hydroxyl species result in low Er solubility at rock equilibrated pH conditions due to an expanded hydroxyl predominance zone, but lower aqueous complex stability overall, whereas previous models suggest greater stability for aqueous Er species. Furthermore, ErPO 4 solubility increases with decreasing temperature due to the deprotonation of HCl, which increases the acidity of hydrothermal fluids and the availability of Cl - to complex with the REE. These simulations highlight how fluid-rock reaction and temperature affect the mobility of REE in hydrothermal ore-forming systems.
A comprehensive analysis of transient pressure and rate data from CO 2 storage projects in a depleted pinnacle reef oil field complex, Michigan, USA
Pressure and rate data are commonly recorded as part of a basic monitoring program in CCS projects. This paper discusses the application of multiple analytical techniques to interpret pressure and rate transient data from CO 2 injection and storage operations. The techniques of interest, i.e., injection-falloff analysis, injectivity/productivity index analysis and pressure pulse arrival time analysis, are commonly used in the oil and gas industry to assess reservoir properties, but not well known in the CCS literature (especially the last two). Injection-falloff analysis involves log-log pressure derivative plotting for the falloff data and history-matching of the entire injection-falloff sequence to determine permeability. In the injectivity/productivity index analysis, rate-normalized pressure buildup is plotted against material balance time or ratio of cumulative injection to injection rate to determine the injectivity index (ratio of injection rate to stabilized pressure buildup) which can be related to the permeability-thickness product. The arrival time analysis identifies the arrival of a pressure disturbance (~0.1 psi change from ambient) to determine the hydraulic diffusivity from which permeability can be estimated. The applicability of these techniques is demonstrated via illustrative examples from multiple wells in different pinnacle carbonate reefs undergoing CO 2 -EOR in Northern Michigan. The paper ends with a discussion of the relative merits of each interpretive technique, as well as recommendations that could be useful for other field projects.
A direct numerical simulation study for confined non-isothermal jet impingement at moderate nozzle-to-plate distances: Capturing jet-to-ambient density effects
A direct numerical simulation (DNS) campaign is deployed for a series of confined downward oriented, non-isothermal turbulent impinging jet configurations. A baseline Reynolds number of 9960 is obtained through a precursor DNS pipe flow simulation (Re τ = 505). Three jet temperature configurations (confinement height to nozzle diameter of three) enter a cylindrical domain that share ambient and impingement plate temperatures (298.15 K). The range of jet temperatures are crafted such that the ratio of inlet to ambient density varies from unity to 0.52, showcasing the effect of density disparity on flow characteristics such as core collapse, radial mixing of momentum and energy, near-wall stagnation behavior, wall-jet profiles, and large-scale vortical structures. Surface quantities provided include mean radial heat flux and wall-shear stress profiles, and heat flux histograms at select radial stations. Results showcase increased radial normal stresses for higher temperature jets that support increased mixing, resulting in large-scale recirculation structures that are smaller, while retaining similar normalized radial wall profiles for shear stress, heat flux and pressure. Furthermore, radial plots for wall shear stress and Nusselt number showcase strong radial decay as compared to previous configurations that share similar jet and ambient temperatures. For the 373.15 K case, a Gaussian-like histogram for heat fluxes at the impingement plate transitions to a log-normal profile as radial distances increase. In contrast, the 573.15 K configuration displays a bi-modal heat flux characteristic at the impingement plate, and in similar manner to the moderate temperature counterpart, transitions to a log-normal profile at larger radial distances.
Combined Recipe for Clinical Target Volume and Planning Target Volume Margins
Purpose: To develop a combined recipe for clinical target volume (CTV) and planning target volume (PTV) margins. Methods and Materials: A widely accepted PTV margin recipe is M{sub geo} = aΣ{sub geo} + bσ{sub geo}, with Σ{sub geo} and σ{sub geo} standard deviations (SDs) representing systematic and random geometric uncertainties, respectively. On the basis of histopathology data of breast and lung tumors, we suggest describing the distribution of microscopic islets around the gross tumor volume (GTV) by a half-Gaussian with SD Σ{sub micro}, yielding as possible CTV margin recipe: M{sub micro} = ƒ(N{sub i}) × Σ{sub micro}, with N{sub i} the average number of microscopic islets per patient. To determine ƒ(N{sub i}), a computer model was developed that simulated radiation therapy of a spherical GTV with isotropic distribution of microscopic disease in a large group of virtual patients. The minimal margin that yielded D{sub min} <95% in maximally 10% of patients was calculated for various Σ{sub micro} and N{sub i}. Because Σ{sub micro} is independent of Σ{sub geo}, we propose they should be added quadratically, yielding for a combined GTV-to-PTV margin recipe: M{sub GTV-PTV} = √([aΣ{sub geo}]{sup 2} + [ƒ(N{sub i})Σ{sub micro}]{sup 2}) + bσ{sub geo}. This was validated by the computer model through numerous simultaneous simulations of microscopic and geometric uncertainties. Results: The margin factor ƒ(N{sub i}) in a relevant range of Σ{sub micro} and N{sub i} can be given by: ƒ(N{sub i}) = 1.4 + 0.8log(N{sub i}). Filling in the other factors found in our simulations (a = 2.1 and b = 0.8) yields for the combined recipe: M{sub GTV-PTV} = √((2.1Σ{sub geo}){sup 2} + ([1.4 + 0.8log(N{sub i})] × Σ{sub micro}){sup 2}) + 0.8σ{sub geo}. The average margin difference between the simultaneous simulations and the above recipe was 0.2 ± 0.8 mm (1 SD). Calculating M{sub geo} and M{sub micro} separately and adding them linearly overestimated PTVs by on average 5 mm. Margin recipes based on tumor control probability (TCP) instead of D{sub min} criteria yielded similar results. Conclusions: A general recipe for GTV-to-PTV margins is proposed, which shows that CTV and PTV margins should be added in quadrature instead of linearly.
Fat Composition Changes in Bone Marrow During Chemotherapy and Radiation Therapy
Purpose: To quantify changes in bone marrow fat fraction and determine associations with peripheral blood cell counts. Methods and Materials: In this prospective study, 19 patients received either highly myelotoxic treatment (radiation therapy plus cisplatin, 5-fluorouracil mitomycin C [FU/MMC], or cisplatin/5-FU/cetuximab) or less myelotoxic treatment (capecitabine-radiation therapy or no concurrent chemotherapy). Patients underwent MR imaging and venipuncture at baseline, midtreatment, and posttreatment visits. We performed mixed effects modeling of the mean proton density fat fraction (PDFF[%]) by linear time, treatment, and vertebral column region (lumbar [L]4-sacral [S]2 vs thoracic [T]10-L3 vs cervical[C]3-T9), while controlling for cumulative mean dose and other confounders. Spearman rank correlations were performed by white blood cell (WBC) counts versus the differences in PDFF(%) before and after treatment. Results: Cumulative mean dose was associated with a 0.43% per Gy (P=.004) increase in PDFF(%). In the highly myelotoxic group, we observed significant changes in PDFF(%) per visit within L4-S2 (10.1%, P<.001) and within T10-L3 (3.93%, P=.01), relative to the reference C3-T9. In the less myelotoxic group, we did not observe significant changes in PDFF(%) per visit according to region. Within L4-S2, we observed a significant difference between treatment groups in the change in PDFF(%) per visit (5.36%, P=.04). Rank correlations of the inverse log differences in WBC versus the differences in PDFF(%) overall and within T10-S2 ranged from 0.69 to 0.78 (P<.05). Rank correlations of the inverse log differences in absolute neutrophil counts versus the differences in PDFF(%) overall and within L4-S2 ranged from 0.79 to 0.81 (P<.05). Conclusions: Magnetic resonance imaging fat quantification is sensitive to marrow composition changes that result from chemoradiation therapy. These changes are associated with peripheral blood cell counts. This study supports a rationale for bone marrow-sparing treatment planning to reduce the risk of hematologic toxicity.
Prospective Randomized Trial of Prone Accelerated Intensity Modulated Breast Radiation Therapy With a Daily Versus Weekly Boost to the Tumor Bed
Purpose: To report the results of a prospective randomized trial comparing a daily versus weekly boost to the tumor cavity during the course of accelerated radiation to the breast with patients in the prone position. Methods and Materials: From 2009 to 2012, 400 patients with stage 0 to II breast cancer who had undergone segmental mastectomy participated in an institutional review board–approved trial testing prone breast radiation therapy to 40.5 Gy in 15 fractions 5 d/wk to the whole breast, after randomization to a concomitant daily boost to the tumor bed of 0.5 Gy, or a weekly boost of 2 Gy, on Friday. The present noninferiority trial tested the primary hypothesis that a weekly boost produced no more acute toxicity than did a daily boost. The recurrence-free survival was estimated for both treatment arms using the Kaplan-Meier method; the relative risk of recurrence or death was estimated, and the 2 arms were compared using the log-rank test. Results: At a median follow-up period of 45 months, no deaths related to breast cancer had occurred. The weekly boost regimen produced no more grade ≥2 acute toxicity than did the daily boost regimen (8.1% vs 10.4%; noninferiority Z = −2.52; P=.006). No statistically significant difference was found in the cumulative incidence of long-term fibrosis or telangiectasia of grade ≥2 between the 2 arms (log-rank P=.923). Two local and two distant recurrences developed in the daily treatment arm and three local and one distant developed in the weekly arm. The 4-year recurrence-free survival rate was not different between the 2 treatment arms (98% for both arms). Conclusions: A tumor bed boost delivered either daily or weekly was tolerated similarly during accelerated prone breast radiation therapy, with excellent control of disease and comparable cosmetic results.
Prognostic Role of Pre–Radiation Therapy {sup 18}F-Fluorodeoxyglucose Positron Emission Tomography for Primary Mediastinal B-Cell Lymphomas Treated with R-CHOP or R-CHOP-Like Chemotherapy Plus Radiation
Purpose: To validate, in a monoinstitutional cohort with extended follow-up, that post–rituximab chemotherapy (R-CT) {sup 18}F-fluorodeoxyglucose positron emission tomography ({sup 18}FDG-PET) is a prognostic factor allowing discrimination of primary mediastinal B-cell lymphoma (PMBCL) patients at higher risk for progression after radiation therapy. Methods and Materials: We analyzed 51 patients, and {sup 18}FDG-PET scans were re-examined evaluating both the Deauville 5-point scale (D5PS) score and the standardized uptake value (SUV) of residual activity, if present. These parameters were then tested by univariate analysis for a potential correlation with progression-free survival (PFS) as the primary study endpoint. Results: Median follow-up time was 51 months (range, 9-153 months). After R-CT, D5PS score was 1 in 10 (19.6%), 2 in 11 (21.6%), 3 in 7 (13.8%), 4 in 17 (33.3%), and 5 in 6 patients (11.7%). Forty-three out of 51 patients (84.3%) had an SUV{sub max} ≤5, and 8 out of 51 (15.7%) had an SUV{sub max} ≥5. Overall, 6 patients experienced progression or relapse: 1 had a D5PS score 2 (with SUV{sub max} ≤5), and 5 had a D5PS score 5 (and SUV{sub max} ≥5). Patients with a D5PS score 5 showed significantly lower PFS rates versus all other scores (log-rank P<.001), as did patients with SUV{sub max} ≥5 when compared with those with SUV{sub max} ≤5 (log-rank P<.001). Conclusions: The present study confirmed the prognostic role of {sup 18}FDG-PET after R-CT, with patients with a D5PS score of 5 and/or an SUV{sub max} ≥5 being at high risk of progression/relapse after RT.
The Effect of Biologically Effective Dose and Radiation Treatment Schedule on Overall Survival in Stage I Non-Small Cell Lung Cancer Patients Treated With Stereotactic Body Radiation Therapy
Purpose: To determine the effect of biologically effective dose (BED{sub 10}) and radiation treatment schedule on overall survival (OS) in patients with early-stage non-small cell lung cancer (NSCLC) undergoing stereotactic body radiation therapy (SBRT). Methods and Materials: Using data from 65 treatment centers in the United States, we retrospectively reviewed the records of T1-2 N0 NSCLC patients undergoing SBRT alone from 2006 to 2014. Biologically relevant covariates, including dose per fraction, number of fractions, and time between fractions, were used to quantify BED{sub 10} and radiation treatment schedule. The linear-quadratic equation was used to calculate BED{sub 10} and to generate a dichotomous dose variable of <105 Gy versus ≥105 Gy BED{sub 10}. The primary outcome was OS. We used the Kaplan-Meier method, the log–rank test, and Cox proportional hazards regression with propensity score matching to determine whether prescription BED{sub 10} was associated with OS. Results: We identified 747 patients who met inclusion criteria. The median BED{sub 10} was 132 Gy, and 59 (7.7%) had consecutive-day fractions. Median follow-up was 41 months, and 452 patients (60.5%) had died by the conclusion of the study. The 581 patients receiving ≥105 Gy BED{sub 10} had a median survival of 28 months, whereas the 166 patients receiving <105 Gy BED{sub 10} had a median survival of 22 months (log–rank, P=.01). Radiation treatment schedule was not a significant predictor of OS on univariable analysis. After adjusting for T stage, sex, tumor histology, and Eastern Cooperative Oncology Group performance status, BED{sub 10} ≥105 Gy versus <105 Gy remained significantly associated with improved OS (hazard ratio 0.78, 95% confidence interval 0.62-0.98, P=.03). Propensity score matching on imbalanced variables within high- and low-dose cohorts confirmed a survival benefit with higher prescription dose. Conclusions: We found that dose escalation to 105 Gy BED{sub 10} and beyond may improve survival in NSCLC patients treated with SBRT.