Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “sublimation”

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.

At least 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Hyperfine spectroscopy and laser cooling of the fermionic isotopes 47 Ti and 49 Ti

Here, we report on magneto-optical trapping of the two fermionic isotopes of atomic titanium, 47 Ti and 49 Ti . Unlike the even mass-number isotopes, which were recently laser cooled, 47 Ti and 49 Ti have nonzero nuclear spins and, consequently, their atomic levels are split by hyperfine structure. Combining and comparing theoretical calculations and atomic beam-spectroscopy measurements, we determine the hyperfine structures and isotope shifts of the 3⁢𝑑 2 ⁢4⁢𝑠 2 𝑎 3 ⁢𝐹 4 → 3⁢𝑑 2 ⁢( 3 𝑃)⁢4⁢𝑠⁢4⁢𝑝⁢( 3 𝑃 𝑜 ) 𝑦 5 ⁢𝐷$^o_4$ optical-pumping transition at optical wavelength 391 nm and the 3⁢𝑑 3 ⁢(4𝐹)⁢4⁢𝑠 𝑎 5 ⁢𝐹 5 → 3⁢𝑑 3 ⁢( 4 𝐹)⁢4⁢𝑝 𝑦 5 ⁢𝐺$^𝑜_6$ laser-cooling transition at wavelength 498 nm. With this information, we produce magneto-optical traps of both 47 Ti and 49 Ti by applying two additional tones of light to repump atoms to the maximum-spin states on the laser-cooling transition. Directly loading from the atomic flux of a titanium sublimation pump, we produce 47 Ti and 49 Ti traps with 731(190) and 1142(240) atoms, and with lifetimes of 330⁢(15) and 310⁢(8)⁢ ms, respectively.

atomic spectra↗
DOE OSTI2024

How Accurate Are Simulations and Experiments for the Lattice Energies of Molecular Crystals?

Molecular crystals play a central role in a wide range of scientific fields, including pharmaceuticals and organic semiconductor devices. However, they are challenging systems to model accurately with computational approaches because of a delicate interplay of intermolecular interactions such as hydrogen bonding and Van der Waals dispersion forces. Here, by exploiting recent algorithmic developments, we report the first set of diffusion Monte Carlo lattice energies for all 23 molecular crystals in the popular and widely used X23 dataset. Comparisons with previous state-of-the-art lattice energy predictions (on a subset of the dataset) and a careful analysis of experimental sublimation enthalpies reveals that high-accuracy computational methods are now at least as reliable as (computationally derived) experiments for the lattice energies of molecular crystals. Overall, this work demonstrates the feasibility of high-level explicitly correlated electronic structure methods for broad benchmarking studies in complex condensed phase systems, and signposts a route towards closer agreement between experiment and simulation. Published by the American Physical Society 2024

Physics↗
DOE OSTI2020

Arsenic Doping of Polycrystalline CdSeTe Devices for Microsecond Life-times with High Carrier Concentrations

We report seminal advances in fabrication and understanding of group V (As) doped thin-film polycrystalline CdTe-based solar cells. The devices are fabricated using a novel approach, by sublimating layers of CdSeTe and CdSeTe:As. This new method allowed us to achieve minority carrier lifetime of over 1 μs, carrier concentration of more than 5×l0 15 cc -1 and external radiative efficiency of over 2 % in a device configuration. We find an increase in open-circuit voltage when comparing As-doped, Cu-doped and undoped devices. The choice of CdSeTe instead of a CdTe-only absorber has been explained using first-principle density functional theory model. A SCAPS device model is used to analyze the potential causes for lower open-circuit voltage.

Munshi, Amit H.↗
DOE OSTI2023

The Surface Atmosphere Integrated Field Laboratory (SAIL) Campaign

The science of mountainous hydrology spans the atmosphere through the bedrock and inherently crosses physical and disciplinary boundaries: land-atmosphere interactions in complex terrain enhance clouds and precipitation, while watersheds retain and release water over a large range of spatial and temporal scales. Limited observations in complex terrain challenge efforts to improve predictive models of the hydrology in the face of rapid changes. The Upper Colorado River exemplifies these challenges, especially with ongoing mismatches between precipitation, snowpack, and discharge. Consequently, the U.S. Department of Energy’s (DOE) Atmospheric Radiation Measurement (ARM) user facility has deployed an observatory to the East River Watershed near Crested Butte, Colorado between September 2021 and June 2023 to measure the main atmospheric drivers of water resources, including precipitation, clouds, winds, aerosols, radiation, temperature and humidity. This effort, called the Surface Atmosphere Integrated Field Laboratory (SAIL), is also working in tandem with DOE-sponsored surface and subsurface hydrologists and other federal, state, and local partners. SAIL data can be benchmarks for model development by producing a wide range of observational information on precipitation and its associated processes, including those processes that impact snowpack sublimation and redistribution, aerosol direct radiative effects in the atmosphere and in the snowpack, aerosol impacts on clouds and precipitation, and processes controlling surface fluxes of energy and mass. Preliminary data from SAIL’s first year showcase the rich information content in SAIL’s many data-streams and support testing hypotheses that will ultimately improve scientific understanding and predictability of Upper Colorado River hydrology in 2023 and beyond.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Wet-Bulb Temperature from Pressure, Relative Humidity, and Air Temperature

Abstract The algorithms currently in use for calculating the wet-bulb temperature from pressure, relative humidity, and air temperature have errors of a degree or more. Here, the equation for the thermodynamic wet-bulb temperature is derived using the Rankine–Kirchhoff approximations, which are highly accurate for meteorological and climatological applications. Likewise, equations are derived for the thermodynamic ice-bulb temperature and for the psychrometric wet-bulb and ice-bulb temperatures. The equations are fast to solve, requiring only a microsecond on a current laptop computer, and the results match empirical data to within hundredths of a degree. Furthermore, the equations reveal the existence of air temperatures and humidities that generate bistability: Evaporating liquid water is stable at a liquid temperature above freezing and, simultaneously, sublimating ice is stable at an ice temperature below freezing.

Romps, David M↗
DOE OSTI2020

A New Method for Ice-Ice Aggregation in the Adaptive Habit Model

In this paper, a novel methodology for modeling ice-ice aggregation is presented. This methodology combines a modified hydrodynamic collection algorithm with bulk aggregate characteristic information from an offline simulator that collects ice particles, namely, the Ice Particle and Aggregate Simulator, and has been implemented into the Adaptive Habit Microphysics scheme in the Weather Research and Forecasting Model. Aggregates, or snow, are formed via collection of cloud ice particles, where initial ice characteristics and the resulting geometry determine aggregate characteristics. Upon implementation, idealized squall-line simulations are performed to examine the new methodology in comparison with commonly used bulk microphysics schemes. It is found that the adaptive habit aggregation parameterization develops snow and reduces ice mass and number concentrations compared to other schemes. The development of aggregates through the new methodology cascades into other interesting effects, including enhancements in ice and snow growth, as well as homogeneous freezing. Further microphysical analyses reveal varying sensitivities, where snow processes are most sensitive to the new parameterization, followed by ice, then cloud, rain, and graupel processes. Further, the new scheme results in enhancements in surface precipitation due to the persistence of snow at lower altitudes. This persistence is a result of shape-dependent melting and sublimation, increasing the residence time. Moreover, these low-level enhancements are reflected in increases in radar reflectivity at the surface and its spatial distribution. Finally, the ability to predict snow shape and density allows for the simulation of polarimetric radar quantities, resulting in signature enhancements compared to schemes that do not consider spatial and temporal variations in snow shape and density.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Effects from Time Dependence of Ice Nucleus Activity for Contrasting Cloud Types

Abstract The role of time-dependent freezing of ice nucleating particles (INPs) is evaluated with the “Aerosol–Cloud” (AC) model in 1) deep convection observed over Oklahoma during the Midlatitude Continental Convective Cloud Experiment (MC3E), 2) orographic clouds observed over North California during the Atmospheric Radiation Measurement (ARM) Cloud Aerosol Precipitation Experiment (ACAPEX), and 3) supercooled, stratiform clouds over the United Kingdom, observed during the Aerosol Properties, Processes And Influences on the Earth’s climate (APPRAISE) campaign. AC uses the dynamical core of the WRF Model and has hybrid bin–bulk microphysics and a 3D mesoscale domain. AC is validated against coincident aircraft, ground-based, and satellite observations for all three cases. Filtered concentrations of ice (>0.1–0.2 mm) agree with those observed at all sampled levels. AC predicts the INP activity of various types of aerosol particles with an empirical parameterization (EP), which follows a singular approach (no time dependence). Here, the EP is modified to represent time-dependent INP activity by a purely empirical approach, using our published laboratory observations of time-dependent INP activity. In all simulated clouds, the inclusion of time dependence increases the predicted INP activity of mineral dust particles by 0.5–1 order of magnitude. However, there is little impact on the cloud glaciation because the total ice is mostly (80%–90%) from secondary ice production (SIP) at levels warmer than about −36°C. The Hallett–Mossop process and fragmentation in ice–ice collisions together initiate about 70% of the total ice, whereas fragmentation during both raindrop freezing and sublimation contributes <10%. Overall, total ice concentrations and SIP are unaffected by time-dependent INP activity. In the simulated APPRAISE case, the main causes of persistence of long-lived clouds and precipitation are predicted to be SIP in weak embedded convection and reactivation following recirculation of dust particles in supercooled layer cloud.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Characteristics and effects of aerosols during blowing snow events in the central Arctic

Sea salt aerosol (SSaer) significantly impacts aerosol-radiation and aerosol-cloud interactions, and sublimated blowing snow is hypothesized to be an important SSaer source in polar regions. Understanding blowing snow and other wind-sourced aerosols’ climate relevant properties is needed, especially during winter when Arctic amplification is greatest. However, most of our understanding of blowing snow SSaer comes from modeling studies, and direct observations are sparse. Additionally, SSaer can originate from multiple sources, making it difficult to disentangle emission processes. Here, we present comprehensive observations of wind-sourced aerosol during blowing snow events from the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition in the central Arctic. High wind speed strongly enhances total aerosol number, submicron sodium chloride mass, cloud condensation nuclei concentrations, and scattering coefficients. Generally, the relative response of aerosol properties to wind speed enhancement is strongest in fall when Arctic aerosol concentrations are lowest. Blowing snow events showed similar aerosol and environmental properties across events, apart from occasions with high snow age (>6 days since last snowfall). Coarse-mode number concentrations (>1 μm) are better explained by variability in wind speed averaged over 12-h air mass back trajectories arriving at the MOSAiC site compared to local, instantaneous wind speed, suggesting the importance of regional transport and consideration of air mass history for wind-driven aerosol production. These MOSAiC observations provide new insights into wind-driven aerosol in the central Arctic and may help validate modeling studies and improve model parameterizations particularly for aerosol direct and indirect radiative forcing.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Modeled estimates of permafrost hydrology and related fields for pan-Arctic region over the period 1980-2100

This dataset contains gridded estimates of total runoff (mm/day); subsurface runoff (mm/day); snow water equivalent (SWE, mm/day); vertical upward water flux of combined surface evaporation, evapotranspiration, and snow sublimation (mm/day); and soil temperature (C) across the pan-Arctic basin from simulations of the Permafrost Water Balance Model (PWBM v4). These data were generated to allow researchers to advance understanding of the changes unfolding in the Arctic's terrestrial water cycle under the influences of climate warming, hydrological cycle intensification, and permafrost thaw. The dataset will enable researchers to address important questions spanning multiple disciplines. For example, how will rising river flows exported to the coast impact critical zone ecosystems, near-shore ice dynamics, and ocean biogeochemical processes. Spatial resolution is 25x25 km with N=39926 grid cells spanning the domain. Three forcing datasets were used to drive simulations at an implicit daily time step over the period 1980-2100. Daily air temperature, precipitation, and wind speed data from the WATCH Forcing Data methodology applied to ERA5 reanalysis (W5E5) were used for meteorological forcings over the period 1980-2019. Simulations for the period 1980-2100 used meteorological data from two global climate models; IPSL-CM6A-LR and MPI-ESM1-2-HR. The PWBM hydrological outputs are organized by grid and month in gzipped ASCII text files within each archived tar file. Soil temperature is organized by grid cell and day of year from DOY 205 to DOY 260. No special software is required to work with these data. The package contains spatial domain and river basin information files. These data are most useful for analyses of the spatial and temporal variations in terrestrial hydrology and soil thermal dynamics. A manuscript describing the model simulations and associated analysis has been accepted for publication in The Cryosphere (Rawlins and Karmalkar, 2024).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Groundwater and Surface Water Flow (GSFLOW) model files to explore bedrock circulation depth and porosity in Copper Creek, Colorado

This data package contains integrated hydrological model input and output files for Copper Creek, Colorado (24 km2), a tributary of the East River located in the headwaters of the Upper Colorado River Basin. The model code is the U.S. Geological Survey (USGS) Groundwater and Surface Water Flow (GSFLOW) model. The model contains a 100-m grid resolution and a daily timestep. The land surface model is dynamically linked to a three-dimensional groundwater flow model that allows for streamflow gaining and losing conditions. The groundwater model contains 12 model layers and extends 400 m below land surface. The original Copper Creek model was modified to contain geologic layers representing saprolite, shallow bedrock, and deep bedrock. Endmember depth versus hydraulic conductivity relationships and porosity values for fractured crystalline rock are simulated. For the shallow case, median flow depths occur in the shallow saprolite at depths <8 m, while the deep case promotes a median groundwater flow depth of 100 m. With this modeling framework we compare streamflow response to a plausible worst-case drought lasting up to five years. Streamflow metrics of analysis include average streamflow, fraction of stream network that is dry, no-flow duration, average groundwater flow to streams and time to recovery following the drought. Results and implications are presented in a paper submitted to Geophysical Research Letters titled, "The role of bedrock circulation depth and porosity in mountain streamflow response to prolonged drought" by Rosemary WH. Carroll, Andrew H. Manning and Kenneth H Williams. A Readme.txt file provides instructions on how to download all model files and execute each model scenario. In addition to the GSFLOW output/prms/copper_drought.csv file containing daily basin water stores and fluxes (refer to GSFLOW manual) and the output/prms/copper_drought_statvar.dat file with output defined in the gsflow3.control file (refer to GSFLOW Manual), output files also include spatially distributed daily values of total evapotranspiration, canopy evaporation, precipitation, snowfall, infiltration, snow water equivalent, potential evapotranspiration, recharge, sublimation, soil moisture, contributing interflow, water table elevations, changes in groundwater storage, groundwater evapotranspiration, interbasin groundwater flow (limited to the alluvium below the stream outlet), and surface-groundwater exchanges within the river system.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

A Bottom-up Approach to Improve the Representation of Deep Convective Clouds in Weather and Climate Models (Final Report)

A fundamental yet outstanding question about deep cumulus convection regards how convective updrafts, downdrafts, and surface-based cold pools are inter-related within organized storms, and how these three convective components are modulated by external and internal factors. We considered this question through numerical model simulations as well as data collected during the MC3E field campaign. One of our basic hypotheses was that deep cold pools are forced by wide downdrafts which arise from wide updrafts. This “bottom-up” perspective required consideration of the external controls on updraft width, and then the internal controls on cold-pool characteristics including depth. Model experiments showed that environmental vertical wind shear (an external factor) exerts a particularly strong control on updraft width; analyses of linear and nonlinear dynamic pressure forcing readily explain this control. The experiments also showed that increases in updraft width are accompanied by corresponding increases in downdraft width and cold-pool depth. The updraft–downdraft connection is through precipitation: the wide area of precipitation that will be initiated in a wide updraft will in turn lead to a wide downdraft. The downdraft–cold pool connection is through mass continuity: at the downdraft terminus, the large downdraft mass flux through a wide downdraft must lead to a deep cold pool. Model experiments also showed that the cold-pool depth and other characteristics can be significantly modulated by microphysical processes (internal factors). Specifically, the cold-pool expansion rate, depth, and intensity was found to depend most strongly on the amount of rain evaporation near the ground. However, the sustenance of the cold pool itself depended upon the downdrafts strengthened by sublimation of graupel higher in the storm. Finally, field campaign data revealed compelling inter-relationships between the cold-pool characteristics, and downdraft transience in terms of the microphysics. For example, pulses of hail and graupel (also internal factors) within downdrafts resulted in short-lived increases in the cold-pool depth. However, the cold-pool speed appeared to be relatively independent of the cold-pool depth and intensity, suggesting limitations of density current theory to convective cold pools in nature.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Kinetics of TATB and LX-17 Decomposition by Thermal Analysis

The kinetics of TATB decomposition have been investigated by thermal analysis (thermogravimetric analysis, TGA, and differential scanning calorimetry, DSC) at both LLNL and the University of Alabama Birmingham. This report summarizes progress in understanding from those results within the context of identifying the global reaction mechanism, i.e., sequential, autocatalytic, competitive, or some mixture thereof. A particularly important issue is how pressure affects both the reaction rate and enthalpy. The effect of pressure helps identify dominant mechanisms and construct and constrain a global mechanism that will work over a wide range of conditions. This memorandum is simply a snapshot of current progress and not a definitive interpretation. I have written it to organize work of others and my interpretations in a systematic manner with the intention of stimulating further advances in understanding. SDT data comes from Jen Montgomery at LLNL, and high-pressure DSC comes from Victoria Stanford and Sergey Vyazovkin at University of Alabama at Birmingham. Major observations and conclusions, some new and some not, are (1) outgassing below 150 °C prior to decomposition releases adsorbed water vapor, (2) outgassing between 150 and 250 °C may come primarily from intercalated water (3) initial thermal decomposition above 250 °C or so involves dehydration of TATB to create furazans, (4) quickly heating to higher temperatures in open pans with flowing gas leads primarily to TATB evaporation between 270 and 390 °C, (5) using a pinhole crucible at atmospheric pressure largely eliminates TATB sublimation, (6) at typical thermal analysis heat rates, the main stage of decomposition is described by a single mass loss process, and decomposition of residual material is spread over a wide range of thermal severity, (7) at typical thermal analysis heat rates and constant temperatures above 300 °C, heat release occurs in two overlapping peaks, with the first one a little faster than mass loss and the second on a little slower than mass loss, (8) total heat release increases significantly with heating rate and substantially with increased pressure, and (9) isothermal experiments at 340 °C indicate that the first exotherm increases by more with LX-17 than TATB, which suggests a possible chemical role of Kel-F decomposition.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Photovoltaics Research and Development: Device Architecture for Next-Generation CdTe PV. Cooperative Research and Development Final Report, CRADA Number CRD-17-00662

Despite remarkable efficiency improvements at the cell and module level, thin-film CdTe PV suffers from a significant voltage deficit compared with the ideal. The proposed solution is to build cells with fully depleted absorbers and with energy barriers front and back to direct electrons and holes in appropriate directions. Combined with minimal interfacial recombination, this approach, which is similar to the successful silicon HIT cells, circumvents the need for higher lifetimes and carrier concentration to achieve higher voltage, Colorado State is particularly well-positioned to the address the depleted absorbed architecture. It has robust single-vacuum closed-spaced-sublimation deposition chamber, currently being enhanced with new sources for CdTe alloys and grading. It has made reasonable-efficiency cells with absorbers as thin as 0.3 um, it has demonstrated electron reflection at the rear of the absorber with a continuous CdMgTe layer, it has replaced CdS with MgZnO, which has several advantages including tunable band offsets, and it has achieved 16.8% efficiency. It also has a mature suite of measurement equipment for cell evaluation and the experience to effectively analyze its results for refinement in cell fabrication. Several healthy collaborations will continue within the CdTe community.

14 SOLAR ENERGY↗
DOE OSTI2021

Surface Atmosphere Integrated Field Laboratory (SAIL) Backgrounder

Half of the world relies on water from the mountains, but in a warming world, it is difficult to predict what these resources will look like. The ability to predict future water requires a deep scientific understanding not only of snow and rain, but also of solar and thermal energy, evaporation and transpiration of water, sublimation of snow, impacts of aerosols, and subsurface hydrological processes that control the ultimate discharge from mountainous watersheds.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Developing a Low Cost, High Volume and Scalable Manufacturing Technology for Undoped and Heavily P-Type Doped CdTe Feedstock Materials (Final Report)

The goal of this project was to establish a cost-effective and scalable production approach of feedstock (CdTe and Cd-Se-Te) for the solar industry, with the material properties necessary to improve performance and reduce costs. The primary focus was to develop a CdTe synthesis and growth process that pushed the p-type doping of CdTe to previously unattainable levels, while creating a path to scalable production. In our approach, the first step was bulk CdTe growth doped with Group V dopants (P, As, and Sb) using the Modified Vertical Bridgman (MVB) and high Pressure Bridgman (HPB) techniques, and the second step was examination of films deposited by Close-Spaced Sublimation (CSS) and Vapor Transport Deposition (VTD) using this specialized feedstock.

14 SOLAR ENERGY↗
DOE OSTI2022

Developing CdTe homojunctions applying high throughput deposition

In this project, we attempted to make CdTe homojunctions using single crystal p-type doped substrates and closed-space sublimation epitaxy (CSSE) n-type doped films. The project included modeling which determined the ideal thickness for the n-type layer as being <200 nm. Boules of high concentration indium-doped CdTe and CdSe 0.4 Te 0.6 (CST) were grown using modified vertical Bridgman (MVB) methods. Similarly, iodine-doped CdTe crystals were grown for the first time. Washington State University (WSU) stock phosphorus-doped CdTe was used as the p-type substrate layer. The crystals were characterized by Hall effect and time-resolved photoluminescence (TRPL, for electrical properties), photoluminescence microscopy (for uniformity), X-ray diffraction (for crystal structure), and glow discharge mass spectrometry (GDMS, for dopant and impurity concentration). CdTe:I crystals were also characterized by visible and infrared transmission measurements, and various Cd or Te heat treatments were performed to assess changes in optical and electrical properties. The grown n-type materials – CdTe:In, CdSe0.4Te0.6:In, and CdTe:I – were provided to the National Renewable Energy Laboratory (NREL) for growth of CSSE thick films for characterization by two photon TRPL (for carrier lifetime), electron back-scatter diffraction (EBSD, to assess epitaxy), and Hall effect. Several measurements of secondary ion mass spectroscopy (SIMS) were performed. Nearly 100% of the indium from the crystal was incorporated into the measured thick films, while only ~2-22% of the iodine from the crystal was incorporated. The net result of the diffusion issue is that homojunction devices created using CSSE have a buried homojunction, as indicated by the near infrared peak in the external quantum efficiency (EQE). Various parameterization of front and back contacts suggested that the poor device performance was primarily a result of this buried junction and not due to other effects. There may also be an issue with the CSSE film lifetime in addition to the dopant profile.

14 SOLAR ENERGY↗
DOE OSTI2022

High concentrations of ice: Investigations using polarimetric radar observations combined with in situ measurements and cloud modeling (Final Report)

This DOE-funded joint project had the over-arching goal of understanding the reasons for observed size distributions of ice particles in cold clouds. Its approach involves the use of cloud models and field observations by radar and aircraft of real storms. The project addresses three major research objectives: (1) Utilize a novel polarimetric radar technique to retrieve size distributions and amounts of ice from radar data collected during the previous DOE ARM field campaigns; (2) Evaluate the results of ice microphysical retrievals using available in situ aircraft measurements and other remote sensors (e.g., vertically pointing cloud radars and wind profilers); (3) Improve the treatment of microphysical processes leading to cloud glaciation and ice multiplication in numerical cloud models. The study was performed by three research organizations: University of Oklahoma, The Hebrew University of Jerusalem, Israel, and Lund University, Sweden. The Israel and Swedish partners were funded by the sub-awards from the University of Oklahoma. Cloud modeling studies have been performed by the research teams at The Hebrew University of Jerusalem (HUJ) and Lund University to identify the origins of high concentrations of cloud ice in areas of high ice water content (HIWC). The Hebrew University Cloud Model (HUCM) with full spectral bin microphysics and the Lund University aerosol-cloud (AC) model with a hybrid bin / bulk microphysics scheme complementing HUCM were utilized for simulations. Both research teams had particular focus on secondary ice production (SIP) as one of the possible sources of enhanced ice concentration. The HUJ group suggested a novel concept of ice multiplication during droplet freezing. It is assumed that splintering and droplet fragmentation during droplet freezing takes place because of dendritic growth within a supercooled drop. The resulting simulations of SIP generated small ice in concentrations exceeding hundreds per liter similar to what was observed in the HIWC regions of the tropical storms. The Lund team explored the SIP mechanisms such as breakup of ice particles due to ice-ice collisions and ice sublimation that are expected to dominate the continental storms. They also quantified the impact of homogeneous nucleation of cloud droplets on the total number concentration of ice at very low temperatures near the tops of the clouds. Additionally, the Lund AC model is able to simulate the effect of aerosols of various types (including biological) on the cloud life cycle and the corresponding ice production. The HUCM / AC model was used to simulate one of the “golden” cases of the DOE MC3E campaign on 20 May 2011. The model output was converted into the fields of polarimetric radar variables using the polarimetric radar forward operator developed at the University of Oklahoma and compared with radar observations and in situ microphysical measurements onboard research aircraft. It was demonstrated that specific differential phase KDP is the best radar parameter to identify the HIWC areas and quantify the corresponding ice parameters. For the first time, the shape of the vertical profile of KDP was realistically reproduced by the cloud model with the KDP maximum in the dendritic growth layer (DGL) centered at the -15°C isotherm. The University of Oklahoma team has developed a methodology for polarimetric radar retrievals of such microphysical parameters of ice as ice water content (IWC), mean volume diameter (Dm), and total number concentration (Nt) of ice particles. These retrievals have been validated using in situ aircraft measurements during 6 field campaigns and proved to be quite robust and reliable. This allowed to build the first climatology of the vertical profiles of polarimetric radar variables and retrieved microphysical parameters for the three types of weather systems: continental MCSs, maritime MCSs, and tropical cyclones / hurricanes (Hu and Ryzhkov 2022). The data were collected by a multitude of the WSR-88D radars in 13 continental and 10 maritime MCSs and 11 landfalling hurricanes. The HIWC areas were identified within the examined storms and the corresponding “HIWC statistics” was compared with the “background” one without HIWC. An overarching conclusion of the study is that maritime tropical storms (MCSs and hurricanes) are characterized by smaller size ice in higher concentration compared to the continental MCSs. High ice water content in the HIWC areas is primarily caused by a strong jump in a number concentration of ice particles rather than the increase of their size compared to the “background” environment. This may point to the homogeneous nucleation of excessive amounts of supercooled droplets and / or secondary ice production as the possible origins of HIWC. Such a climatology provides a good observational reference for the modelers to evaluate the performance of their models. As an example, the in-depth analysis of the 20 May 2011 MC3E case shows that the advanced cloud models developed in the course of this study still tend to underestimate the number concentration of ice in the HIWC areas although they succeed in reproducing realistically looking vertical profiles of IWC and Nt. The results of the project research are summarized in 13 journal papers.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Health Hazards of Exposures to Radioiodine

Iodine is a chemical element with atomic number 53. Iodine-127 is stable (non-radioactive) and commonly found in nature. Elemental iodine is a purple-colored solid at room temperature and pressure, but spontaneously sublimates (turns into vapor). Iodine is an essential element for life, and is required for proper functioning of the thyroid. Iodine is present in many foods, and is readily absorbed by the body and concentrated in the thyroid gland. A fraction of iodine ingested or inhaled is rapidly removed by the kidneys. The rest of the inhaled or ingested iodine is absorbed the by thyroid and retained for many months. Iodine has a biological half-life of approximately 120 days in health individuals. The biological half-life can be shorter in individuals with hyperthyroidism, and longer in individuals with hypothyroidism. Iodine has a number of radioactive isotopes, most of which have relatively short half-lives (days or weeks). Short half-life iodine isotopes are useful for a variety of medical applications, including imaging and cancer therapy. For example, Iodine-123 (half-life 13 hours) is commonly used for medical imaging of the thyroid, while iodine-131 (half-life 8 days) is used for suppressing thyroid function in individuals with hyperthyroidism or ablating (killing) thyroid cells to treat thyroid cancer. Iodine-125 (half-life 59 days) is produced in nuclear reactors, and has medical uses. Although iodine-125 can be used for thyroid imaging, Iodine-123 is more commonly used for that purpose because of its shorter half-life and higher-energy emissions. Iodine-125 is more commonly used for cancer treatment, and can be processed into small metal pellets (seeds) inserted directly into a tumor. Iodine-125 emits low-energy x-rays which can kill tumor cells and generally cannot escape the tumor, sparing other tissues. Medical iodine for imaging or treatment is typically administered orally in the form a pill or liquid solution. A typical adult thyroid scan using iodine-123 involves having the patient swallow between one and four 0.1 millicuries pills, with the exact dose dependent on the patient’s weight. This results in a whole-body committed effective dose of 80 – 320 mrem, and a thyroid equivalent dose of 1443 – 5772 mrem. Note that the whole-body effective dose relates to the overall cancer risk, while the larger equivalent dose to the thyroid only indicates that most of this risk is the result of exposure to the thyroid. These doses are considered safe, although the procedure is not recommended for pregnant or breastfeeding women. In contrast, the quantities of iodine-131 used for treatment of hyperthyroidism and thyroid cancers are much higher. For treatment of hyperthyroidism, 4 – 10 millicuries are administered, while for thyroid cancer the administration can range from 50 – 150 millicuries of I-131. In addition to medical exposures, large populations were exposed to radioiodine as result of the atomic bombings of Hiroshima and Nagasaki in Japan, and the Chernobyl nuclear accident. These populations have been carefully followed for many years to assess the effect of their radiation exposures on cancer risk. As a result, a great deal is known about the cancer risks associated with radioiodine exposure. Because iodine is concentrated in the thyroid, the principal risk of exposure to radioiodine is thyroid cancer. Children have the highest risk of thyroid cancer after exposure to radioactive iodine. According to a large study of Japanese atomic bomb survivors, an effective dose of one Sievert (100,000 mrem) has been observed to increase the risk of thyroid cancer by a factor of 9.5 in children aged zero to nine years old, by a factor of 3 in children aged 10 to 19 years old, and by barely detectable amounts in adults. Another way of quantifying the risk from radioiodine exposure is from risk coefficients, which provide the risk per unit intake of radionuclides in terms of both morbidity (any cancer) and mortality (death). Both morbidity (risk of cancer) and mortality (death) risks are shown in the table below. Note that because thyroid cancer is almost never fatal, the morbidity coefficients are much larger than the mortality coefficients.

61 RADIATION PROTECTION AND DOSIMETRY↗