The limit of macroscopic homogeneous ice nucleation at the nanoscale

Nucleation in small volumes of water has garnered renewed interest due to the relevance of pore condensation and freezing under conditions of low partial pressures of water, such as in the upper troposphere. Molecular simulations can in principle provide insight on this process at the molecular scale that is challenging to achieve experimentally. However, there are discrepancies in the literature as to whether the rate in confined systems is enhanced or suppressed relative to bulk water at the same temperature and pressure. In this study, we investigate the extent to which the size of the critical nucleus and the rate at which it grows in thin films of water are affected by the thickness of the film. Our results suggest that nucleation remains bulk-like in films that are barely large enough accommodate a critical nucleus. This conclusion seems robust to the presence of physical confining boundaries. We also discuss the difficulties in unambiguously determining homogeneous nucleation rates in nanoscale systems, owing to the challenges in defining the volume. Our results suggest any impact on a film's thickness on the rate is largely inconsequential for present day experiments.

Does supercooled water retain its universal nucleation behavior under shear at high pressure?

Understanding the nucleation of homogeneous flow systems at high pressures is vital in protein crystallization and cryopreservation, where high pressure prevents the freezing of biological samples. This study examines the behavior of ice nucleation under shear at various pressures and explores the universal nucleation behavior of the sheared systems applied to supercooled water at higher pressures. In this study, the nucleation rates for the TIP4P/Ice model via a seeding method based on extended classical nucleation theory (CNT) are computed at pressures of 1, 100, 500, 700, and 1000 bar and a constant temperature of 240 K. Using extended CNT with explicitly embodying the shear rate, we analyzed the dependence of pressure on the transport and thermodynamic properties. In line with previous studies, we observed that Δμliq-ice and viscosity decrease while diffusivity increases with an increase in pressure. Furthermore, we showed that the dependence of the nucleation rate on shear at higher pressure is non-monotonic, with the maximum at optimal shear rates between 107 and 108 s-1. Our results demonstrate a non-monotonic pressure dependence of the optimal shear rates, which could originate from a violation of the Stokes-Einstein relation.

Anomalous Vapor and Ice Nucleation in Water at Negative Pressures: A Classical Density Functional Theory Study

In contrast to the abundance of work on the anomalous behavior of water, the relationship between the water's thermodynamic anomalies and kinetics of phase transition from metastable water is relatively unexplored. In this work, we have employed classical density functional theory to provide a unified and coherent picture of nucleation (both vapor and ice) from metastable water at negative pressure conditions. Our results suggest a peculiar nonmonotonic temperature dependence of vapor-liquid surface tension at temperatures where vapor-liquid coexistence is metastable with respect to the ice phase. The vapor nucleation barrier on isochoric cooling also shows a nonmonotonic temperature dependence. We further report that, for low density isochores, the temperature of the minimum vapor nucleation barrier (TΔΩv/min*) does not coincide with the temperature of maximum density (TMD) where metastability is maximum. The difference between the TΔΩv/min* and the TMD, however, decreases with increasing the density of the isochore. The vapor nucleation barrier along isobars shows an interesting crossover behavior in the vicinity of the Widom line on lowering the temperature. Our results on the ice nucleation suggest an anomalous retracing behavior of the nucleation barrier along isotherms at negative pressures and theoretically validate the recent findings that the reentrant ice(Ih)-liquid coexistence line can induce a drastic change in the kinetics of ice nucleation. Thus, this study establishes a direct connection between the metastable water's thermodynamic anomalies and the (vapor and ice) nucleation kinetics. In addition, this study provides deeper insights into the origin of the isothermal compressibility maximum on isochoric cooling.

The freezing behavior of aqueous n-alcohol nanodroplets

We generate water-rich aerosols containing 1-propanol and 1-pentanol in a supersonic nozzle to study the effects of these solutes on the freezing behavior of water. Condensation and freezing are characterized by two complementary techniques, pressure trace measurements and Fourier Transform Infrared spectroscopy. When 1-pentanol and 1-propanol are present, condensation occurs at higher temperatures because particle formation from the vapor phase is enhanced by the decrease in interfacial free energy of mixed aqueous-alcohol critical clusters relative to those of pure water. FTIR results suggest that when ∼6 nm radius droplets freeze, the tetrahedral structure of the ice is well preserved up to an overall alcohol mole fraction of 0.031 for 1-propanol and 0.043 for 1-pentanol. In this concentration range, the ice nucleation temperature decreases continuously with increasing 1-propanol concentration, whereas the onset of freezing is not significantly perturbed by 1-pentanol up to a mole fraction of 0.03. Furthermore, once freezing starts the ice nucleation rates in the aqueous-alcohol droplets are very close to those for pure water. In contrast, at the highest mole fractions of either alcohol it is not clear whether droplets freeze to form crystalline ice since the final state of the particles cannot be adequately characterized with the available experimental techniques.

Anomalous Behavior in the Nucleation of Ice at Negative Pressures

Ice nucleation is a phenomenon that, despite the relevant implications for life, atmospheric sciences, and technological applications, is far from being completely understood, especially under extreme thermodynamic conditions. In this work we present a computational investigation of the homogeneous ice nucleation at negative pressures. By means of the seeding technique we estimate the size of the ice critical nucleus N_{c} for the TIP4P/Ice water model. This is done along the isotherms 230, 240, and 250 K, from positive to negative pressures until reaching the liquid-gas kinetic stability limit (where cavitation cannot be avoided). We find that N_{c} is nonmonotonic upon depressurization, reaching a minimum at negative pressures in the doubly metastable region of water. According to classical nucleation theory we establish the nucleation rate J and the surface tension γ, revealing a retracing behavior of both when the liquid-gas kinetic stability limit is approached. We also predict a reentrant behavior of the homogeneous nucleation line. The reentrance of these properties is related to the reentrance of the coexistence line at negative pressure, revealing new anomalies of water. The results of this work suggest the possibility of having metastable samples of liquid water for long times at negative pressure provided that heterogeneous nucleation is suppressed.

Uncertainty quantification based cloud parameterization sensitivity analysis in the NCAR community atmosphere model

Using uncertainty quantification techniques, we carry out a sensitivity analysis of a large number (17) of parameters used in the NCAR CAM5 cloud parameterization schemes. The LLNL PSUADE software is used to identify the most sensitive parameters by performing sensitivity analysis. Using Morris One-At-a-Time (MOAT) method, we find that the simulations of global annual mean total precipitation, convective, large-scale precipitation, cloud fractions (total, low, mid, and high), shortwave cloud forcing, longwave cloud forcing, sensible heat flux, and latent heat flux are very sensitive to the threshold-relative-humidity-for-stratiform-low-clouds ([Formula: see text] and the auto-conversion-size-threshold-for-ice-to-snow [Formula: see text] The seasonal and regime specific dependence of some parameters in the simulation of precipitation is also found for the global monsoons and storm track regions. Through sensitivity analysis, we find that the Somali jet strength and the tropical easterly jet associated with the south Asian summer monsoon (SASM) show a systematic dependence on [Formula: see text] and [Formula: see text]. The timing of the withdrawal of SASM over India shows a monotonic increase (delayed withdrawal) with an increase in [Formula: see text]. Overall, we find that [Formula: see text], [Formula: see text], [Formula: see text] and [Formula: see text] are the most sensitive cloud parameters and thus are of high priority in the model tuning process, in order to reduce uncertainty in the simulation of past, present, and future climate.

Spaceborne Cloud and Precipitation Radars: Status, Challenges, and Ways Forward

Spaceborne radars offer a unique three-dimensional view of the atmospheric components of the Earth's hydrological cycle. Existing and planned spaceborne radar missions provide cloud and precipitation information over the oceans and land difficult to access in remote areas. A careful look into their measurement capabilities indicates considerable gaps that hinder our ability to detect and probe key cloud and precipitation processes. The international community is currently debating how the next generation of spaceborne radars shall enhance current capabilities and address remaining gaps. Part of the discussion is focused on how to best take advantage of recent advancements in radar and space platform technologies while addressing outstanding limitations. First, the observing capabilities and measurement highlights of existing and planned spaceborne radar missions including TRMM, CloudSat, GPM, RainCube, and EarthCARE are reviewed. Then, the limitations of current spaceborne observing systems, with respect to observations of low-level clouds, midlatitude and high-latitude precipitation, and convective motions, are thoroughly analyzed. Finally, the review proposes potential solutions and future research avenues to be explored. Promising paths forward include collecting observations across a gamut of frequency bands tailored to specific scientific objectives, collecting observations using mixtures of pulse lengths to overcome trade-offs in sensitivity and resolution, and flying constellations of miniaturized radars to capture rapidly evolving weather phenomena. This work aims to increase the awareness about existing limitations and gaps in spaceborne radar measurements and to increase the level of engagement of the international community in the discussions for the next generation of spaceborne radar systems.

Interfacial Concentration Effect Facilitates Heterogeneous Nucleation from Solution

Crystal nucleation from solution plays an important role in environmental, biological, and industrial processes and mainly occurs at interfaces, although the mechanisms are not well understood. We performed nucleation experiments on glycine aqueous solutions and found that an oil-solution interface dramatically accelerates glycine nucleation compared to an air-solution interface. This is surprising given that nonpolar, hydrophobic oil (tridecane) would not be expected to favor heterogeneous nucleation of highly polar, hydrophilic glycine. Molecular dynamics simulations found significantly enhanced vs depleted glycine concentrations at the oil-solution vs air-solution interfaces, respectively. We propose that this interfacial concentration effect facilitates heterogeneous nucleation, and that it is due to dispersion interactions. This interface effect is distinct from previously described mechanisms, including surface functionalization, templating, and confinement and is expected to be present in a wide range of solution systems. This work provides new insight that is essential for understanding and controlling heterogeneous nucleation.

Anthropogenic combustion iron as a complex climate forcer

Atmospheric iron affects the global carbon cycle by modulating ocean biogeochemistry through the deposition of soluble iron to the ocean. Iron emitted by anthropogenic (fossil fuel) combustion is a source of soluble iron that is currently considered less important than other soluble iron sources, such as mineral dust and biomass burning. Here we show that the atmospheric burden of anthropogenic combustion iron is 8 times greater than previous estimates by incorporating recent measurements of anthropogenic magnetite into a global aerosol model. This new estimation increases the total deposition flux of soluble iron to southern oceans (30-90 °S) by 52%, with a larger contribution of anthropogenic combustion iron than dust and biomass burning sources. The direct radiative forcing of anthropogenic magnetite is estimated to be 0.021 W m^-2 globally and 0.22 W m^-2 over East Asia. Our results demonstrate that anthropogenic combustion iron is a larger and more complex climate forcer than previously thought, and therefore plays a key role in the Earth system.

Anthropogenic combustion iron as a complex climate forcer

Atmospheric iron affects the global carbon cycle by modulating ocean biogeochemistry through the deposition of soluble iron to the ocean. Iron emitted by anthropogenic (fossil fuel) combustion is a source of soluble iron that is currently considered less important than other soluble iron sources, such as mineral dust and biomass burning. Here we show that the atmospheric burden of anthropogenic combustion iron is 8 times greater than previous estimates by incorporating recent measurements of anthropogenic magnetite into a global aerosol model. This new estimation increases the total deposition flux of soluble iron to southern oceans (30–90 °S) by 52%, with a larger contribution of anthropogenic combustion iron than dust and biomass burning sources. The direct radiative forcing of anthropogenic magnetite is estimated to be 0.021 W m⁻² globally and 0.22 W m⁻² over East Asia. Our results demonstrate that anthropogenic combustion iron is a larger and more complex climate forcer than previously thought, and therefore plays a key role in the Earth system.

As a source of soluble iron, anthropogenic combustion iron is considered less important than natural sources. Here, the authors combine new measurements with a global aerosol model and show the atmospheric burden of anthropogenic combustion iron to be 8 times greater than previous estimates.

Scenario dependence of future changes in climate extremes under 1.5 °C and 2 °C global warming

The 2015 Paris Agreement aims to limit global warming below 2 °C and pursue efforts to even limit it to 1.5 °C relative to pre-industrial levels. Decision makers need reliable information on the impacts caused by these warming levels for climate mitigation and adaptation measures. We explore the changes in climate extremes, which are closely tied to economic losses and casualties, under 1.5 °C and 2 °C global warming and their scenario dependence using three sets of ensemble global climate model simulations. A warming of 0.5 °C (from 1.5 °C to 2 °C) leads to significant increases in temperature and precipitation extremes in most regions. However, the projected changes in climate extremes under both warming levels highly depend on the pathways of emissions scenarios, with different greenhouse gas (GHG)/aerosol forcing ratio and GHG levels. Moreover, there are multifold differences in several heavily polluted regions, among the scenarios, in the changes in precipitation extremes due to an additional 0.5 °C warming from 1.5 °C to 2 °C. Our results demonstrate that the chemical compositions of emissions scenarios, not just the total radiative forcing and resultant warming level, must be considered when assessing the impacts of global 1.5/2 °C warming.

Challenges in constraining anthropogenic aerosol effects on cloud radiative forcing using present-day spatiotemporal variability

A large number of processes are involved in the chain from emissions of aerosol precursor gases and primary particles to impacts on cloud radiative forcing. Those processes are manifest in a number of relationships that can be expressed as factors dlnX/dlnY driving aerosol effects on cloud radiative forcing. These factors include the relationships between cloud condensation nuclei (CCN) concentration and emissions, droplet number and CCN concentration, cloud fraction and droplet number, cloud optical depth and droplet number, and cloud radiative forcing and cloud optical depth. The relationship between cloud optical depth and droplet number can be further decomposed into the sum of two terms involving the relationship of droplet effective radius and cloud liquid water path with droplet number. These relationships can be constrained using observations of recent spatial and temporal variability of these quantities. However, we are most interested in the radiative forcing since the preindustrial era. Because few relevant measurements are available from that era, relationships from recent variability have been assumed to be applicable to the preindustrial to present-day change. Our analysis of Aerosol Comparisons between Observations and Models (AeroCom) model simulations suggests that estimates of relationships from recent variability are poor constraints on relationships from anthropogenic change for some terms, with even the sign of some relationships differing in many regions. Proxies connecting recent spatial/temporal variability to anthropogenic change, or sustained measurements in regions where emissions have changed, are needed to constrain estimates of anthropogenic aerosol impacts on cloud radiative forcing.

Enhanced marine sulphur emissions offset global warming and impact rainfall

Artificial fertilisation of the ocean has been proposed as a possible geoengineering method for removing carbon dioxide from the atmosphere. The associated increase in marine primary productivity may lead to an increase in emissions of dimethyl sulphide (DMS), the primary source of sulphate aerosol over remote ocean regions, potentially causing direct and cloud-related indirect aerosol effects on climate. This pathway from ocean fertilisation to aerosol induced cooling of the climate may provide a basis for solar radiation management (SRM) geoengineering. In this study, we investigate the transient climate impacts of two emissions scenarios: an RCP4.5 (Representative Concentration Pathway 4.5) control; and an idealised scenario, based on RCP4.5, in which DMS emissions are substantially enhanced over ocean areas. We use mini-ensembles of a coupled atmosphere-ocean configuration of CESM1(CAM5) (Community Earth System Model version 1, with the Community Atmosphere Model version 5). We find that the cooling effect associated with enhanced DMS emissions beneficially offsets greenhouse gas induced warming across most of the world. However, the rainfall response may adversely affect water resources, potentially impacting human livelihoods. These results demonstrate that changes in marine phytoplankton activity may lead to a mixture of positive and negative impacts on the climate.

Evaluations of tropospheric aerosol properties simulated by the community earth system model with a sectional aerosol microphysics scheme

A sectional aerosol model (CARMA) has been developed and coupled with the Community Earth System Model (CESM1). Aerosol microphysics, radiative properties, and interactions with clouds are simulated in the size-resolving model. The model described here uses 20 particle size bins for each aerosol component including freshly nucleated sulfate particles, as well as mixed particles containing sulfate, primary organics, black carbon, dust, and sea salt. The model also includes five types of bulk secondary organic aerosols with four volatility bins. The overall cost of CESM1-CARMA is approximately ∼2.6 times as much computer time as the standard three-mode aerosol model in CESM1 (CESM1-MAM3) and twice as much computer time as the seven-mode aerosol model in CESM1 (CESM1-MAM7) using similar gas phase chemistry codes. Aerosol spatial-temporal distributions are simulated and compared with a large set of observations from satellites, ground-based measurements, and airborne field campaigns. Simulated annual average aerosol optical depths are lower than MODIS/MISR satellite observations and AERONET observations by ∼32%. This difference is within the uncertainty of the satellite observations. CESM1/CARMA reproduces sulfate aerosol mass within 8%, organic aerosol mass within 20%, and black carbon aerosol mass within 50% compared with a multiyear average of the IMPROVE/EPA data over United States, but differences vary considerably at individual locations. Other data sets show similar levels of comparison with model simulations. The model suggests that in addition to sulfate, organic aerosols also significantly contribute to aerosol mass in the tropical UTLS, which is consistent with limited data.

Short-term time step convergence in a climate model

This paper evaluates the numerical convergence of very short (1 h) simulations carried out with a spectral-element (SE) configuration of the Community Atmosphere Model version 5 (CAM5). While the horizontal grid spacing is fixed at approximately 110 km, the process-coupling time step is varied between 1800 and 1 s to reveal the convergence rate with respect to the temporal resolution. Special attention is paid to the behavior of the parameterized subgrid-scale physics. First, a dynamical core test with reduced dynamics time steps is presented. The results demonstrate that the experimental setup is able to correctly assess the convergence rate of the discrete solutions to the adiabatic equations of atmospheric motion. Second, results from full-physics CAM5 simulations with reduced physics and dynamics time steps are discussed. It is shown that the convergence rate is 0.4-considerably slower than the expected rate of 1.0. Sensitivity experiments indicate that, among the various subgrid-scale physical parameterizations, the stratiform cloud schemes are associated with the largest time-stepping errors, and are the primary cause of slow time step convergence. While the details of our findings are model specific, the general test procedure is applicable to any atmospheric general circulation model. The need for more accurate numerical treatments of physical parameterizations, especially the representation of stratiform clouds, is likely common in many models. The suggested test technique can help quantify the time-stepping errors and identify the related model sensitivities.

Cirrus cloud seeding: a climate engineering mechanism with reduced side effects?

Climate engineering, the intentional alteration of Earth's climate, is a multifaceted and controversial topic. Numerous climate engineering mechanisms (CEMs) have been proposed, and the efficacies and potential undesired consequences of some of them have been studied in the safe environments of numerical models. Here, we present a global modelling study of a so far understudied CEM, namely the seeding of cirrus clouds to reduce their lifetimes in the upper troposphere, and hence their greenhouse effect. Different from most CEMs, the intention of cirrus seeding is not to reduce the amount of solar radiation reaching Earth's surface. This particular CEM rather targets the greenhouse effect, by reducing the trapping of infrared radiation by high clouds. This avoids some of the caveats that have been identified for solar radiation management, for example, the delayed recovery of stratospheric ozone or drastic changes to Earth's hydrological cycle. We find that seeding of mid- and high-latitude cirrus clouds has the potential to cool the planet by about 1.4 K, and that this cooling is accompanied by only a modest reduction in rainfall. Intriguingly, seeding of the 15% of the globe with the highest solar noon zenith angles at any given time yields the same global mean cooling as a seeding strategy that involves 45% of the globe. In either case, the cooling is strongest at high latitudes, and could therefore serve to prevent Arctic sea ice loss. With the caveat that there are still significant uncertainties associated with ice nucleation in cirrus clouds and its representation in climate models, cirrus seeding appears to represent a powerful CEM with reduced side effects.

Impact of Antarctic mixed-phase clouds on climate

Precious little is known about the composition of low-level clouds over the Antarctic Plateau and their effect on climate. In situ measurements at the South Pole using a unique tethered balloon system and ground-based lidar reveal a much higher than anticipated incidence of low-level, mixed-phase clouds (i.e., consisting of supercooled liquid water drops and ice crystals). The high incidence of mixed-phase clouds is currently poorly represented in global climate models (GCMs). As a result, the effects that mixed-phase clouds have on climate predictions are highly uncertain. We modify the National Center for Atmospheric Research (NCAR) Community Earth System Model (CESM) GCM to align with the new observations and evaluate the radiative effects on a continental scale. The net cloud radiative effects (CREs) over Antarctica are increased by +7.4 Wm(-2), and although this is a significant change, a much larger effect occurs when the modified model physics are extended beyond the Antarctic continent. The simulations show significant net CRE over the Southern Ocean storm tracks, where recent measurements also indicate substantial regions of supercooled liquid. These sensitivity tests confirm that Southern Ocean CREs are strongly sensitive to mixed-phase clouds colder than -20 °C.