Search NASA⌕ Search

SEARCH · Search NASA

Results for “CLIMATE MODELS”

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 1,675 records · Page 93

Natural variation in growth and yield to waterlogging across climate-diverse pennycress accessions

Many winter annual crops, such as pennycress ( Thlaspi arvense L.), are subjected to heavy precipitation events during their growing season. Therefore, it is essential to identify pennycress accessions with natural variation in flooding resilience. We used climate modeling data to assess spring soil moisture levels in the geographic origins of 471 natural pennycress accessions. We selected 34 accessions with variation in predicted soil moisture and tested survivability under prolonged waterlogging at the rosette stage. It took seven weeks for the first accessions to die, indicating that pennycress is hardy to prolonged waterlogging at the vegetative stage. Furthermore, we chose ‘susceptible’ and ‘tolerant’ accessions to waterlog for one week at the reproductive stage, the growth stage aligned with spring rainfall. Six accessions had significantly reduced seed weight at maturity and two had minimal impacts on growth and seed yield after waterlogging and can be further explored for adaptive traits.

59 BASIC BIOLOGICAL SCIENCES↗

Multi-Variate Factorisation of Numerical Simulations

Factorisation (also known as “factor separation”) is widely used in the analysis of numerical simulations. It allows changes in properties of a system to be attributed to changes in multiple variables associated with that system. There are many possible factorisation methods; here we discuss three previously proposed factorisations that have been applied in the field of climate modelling: the linear factorisation, the Stein and Alpert (1993) factorisation, and the Lunt et al. (2012) factorisation. We show that, when more than two variables are being considered, none of these three methods possess all four properties of “uniqueness”, “symmetry”, “completeness”, and “purity”. Here, we extend each of these factorisations so that they do possess these properties for any number of variables, resulting in three factorisations – the “linear-sum” factorisation, the “shared-interaction” factorisation, and the “scaled-residual” factorisation. We show that the linear-sum factorisation and the shared-interaction factorisation reduce to be identical in the case of four or fewer variables, and we conjecture that this holds for any number of variables. We present the results of the factorisations in the context of three past studies that used the previously proposed factorisations.

Numerical simulations↗

Traveling Weather Disturbances in Mars Southern Extratropics: Sway of the Great Impact Basins

As on Earth, between late autumn and early spring on Mars middle and high latitudes within its atmosphere support strong mean thermal contrasts between the equator and poles (i.e. "baroclinicity"). Data collected during the Viking era and observations from both the Mars Global Surveyor (MGS) and Mars Reconnaissance Orbiter (MRO) indicate that this strong baroclinicity supports vigorous, large-scale eastward traveling weather systems (i.e. transient synoptic-period waves). Within a rapidly rotating, differentially heated, shallow atmosphere such as on Earth and Mars, such large-scale, extratropical weather disturbances are critical components of the global circulation. These wave-like disturbances act as agents in the transport of heat and momentum, and moreover generalized tracer quantities (e.g., atmospheric dust, water vapor and water-ice clouds) between low and high latitudes of the planet. The character of large-scale, traveling extratropical synoptic-period disturbances in Mars' southern hemisphere during late winter through early spring is investigated using a high-resolution Mars global climate model (Mars GCM). This global circulation model imposes interactively lifted (and radiatively active) dust based on a threshold value of the instantaneous surface stress. Compared to observations, the model exhibits a reasonable "dust cycle" (i.e. globally averaged, a more dusty atmosphere during southern spring and summer occurs). In contrast to their northern-hemisphere counterparts, southern synoptic-period weather disturbances and accompanying frontal waves have smaller meridional and zonal scales, and are far less intense synoptically. Influences of the zonally asymmetric (i.e. east-west varying) topography on southern large-scale weather disturbances are examined. Simulations that adapt Mars' full topography compared to simulations that utilize synthetic topographies emulating essential large-scale features of the southern middle latitudes indicate that Mars' transient barotropic/baroclinic eddies are significantly influenced by the great impact basins of this hemisphere (e.g., Argyre and Hellas). In addition, the occurrence of a southern storm zone in late winter and early spring is keyed particularly to the western hemisphere via orographic influences arising from the Tharsis highlands, and the Argyre and Hellas impact basins. Geographically localized transient-wave activity diagnostics are constructed that illuminate fundamental differences amongst such simulations and these are described.

Weather disturbances↗

A Future Mars Environment for Science and Exploration

Today, Mars is arid and cold with a very thin atmosphere that has significant frozen and underground water resources. The thin atmosphere prevents liquid water from residing permanently on its surface and makes it difficult to land missions since it is not thick enough to completely facilitate a soft landing. In its past, under the influence of a significant greenhouse effect, Mars must have had a significant water ocean covering perhaps 30% of the northern hemisphere. Mars lost its protective magnetosphere and therefore much of its atmosphere around 3 Ga ago, due to the solar wind. The atmospheric loss into the solar wind is somewhat balanced by the outgassing of the Mars interior and crust that contributes to the existing atmosphere leading to a global-mean surface atmosphere of ~6 mbar pressure currently. By using our extensive simulation tools and physics capabilities in Space Weather and Mars global climate modeling, we have started to explore the effects on Mars of placing an artificial magnetic dipole field at the Mars L1 Lagrange point putting Mars in a magnetotail. This situation then eliminates many of the solar-wind erosion processes that occur with the planet's ionosphere and upper atmosphere allowing the Martian atmosphere to grow in pressure and bulk temperature over time. Under thicker atmospheres, the global circulation patterns and seasonal changes are much different than at present. An enhanced atmosphere would: allow larger landed mass of equipment to the surface, shield against some cosmic and solar particle radiation, extend the ability for extraction, and provide "open air" greenhouses to exist for plant production, just to name a few. These new conditions on Mars would allow human explorers and researchers to study the planet in much greater detail and enable a truly profound new understanding of the habitability of this planet.

Space Weather and Mars global climate modeling↗

Exploring Causal Relationships and Adjustment Timescales of Aerosol-Cloud Interactions in Geostationary Satellite Observations and CAM6 Using Wavelet Phase Coherence Analysis

We present for the first time within the cloud physics context, the application of wavelet phase coherence analysis to disentangle counteracting physical processes associated with the lead-lag phase difference between cloud-proxy liquid water path (LWP) and aerosol-proxy cloud droplet number concentration ( N d ) in an Eulerian framework using satellite-based observations and climate model outputs. This approach allows us to identify the causality and dominant adjustment timescales governing the correlation between LWP and N d . Satellite observations indicate a more prevalent positive correlation between daytime LWP and N d regardless of whether LWP leads or lags N d . The positive cloud water response, associated with precipitation processes, typically occurs within 1 hr, while the negative response resulting from entrainment drying, usually takes 2–4 hr. CAM6 displays excessively rapid negative responses along with overly strong negative cloud water response and insufficient positive response, leading to a more negative correlation between LWP and N d compared to observations.

54 ENVIRONMENTAL SCIENCES↗

Diagnosing added value of convection-permitting regional models using precipitation event identification and tracking

Dynamical downscaling with high-resolution regional climate models may offer the possibility of realistically reproducing precipitation and weather events in climate simulations. As resolutions fall to order kilometers, the use of explicit rather than parametrized convection may offer even greater fidelity. However, these increased resolutions both allow and require increasingly complex diagnostics for evaluating model fidelity. Here we focus on precipitation evaluation and analyze five 2-month-long dynamically downscaled model runs over the continental United States that employ different convective and microphysics parameterizations, including one high-resolution convection-permitting simulation. All model runs use the Weather Research and Forecasting Model driven by National Center for Environmental Prediction reanalysis data. We show that employing a novel rainstorm identification and tracking algorithm that allocates essentially all rainfall to individual precipitation events allows new insights into model biases. Results include that, at least in these runs, model wet bias is driven by excessive areal extent of individual precipitating events, and that the effect is time-dependent, producing excessive diurnal cycle amplitude. This amplified cycle is driven not by new production of events but by excessive daytime enlargement of long-lived precipitation events. Here, we further show that in the domain average, precipitation biases appear best represented as additive offsets. Of all model configurations evaluated, convection-permitting simulations most consistently reduced biases in precipitation event characteristics.

54 ENVIRONMENTAL SCIENCES↗

Sea Level Rise in the CESM Large Ensemble: The Role of Individual Climate Forcings and Consequences for the Coming Decades

The emergence of a spatial pattern in the externally forced response (FR) of dynamic sea level (DSL) during the altimeter era has recently been demonstrated using climate models but our understanding of its initial emergence, drivers, and implications for the future is poor. Here the anthropogenic forcings of the DSL pattern are explored using the Community Earth System Model Large Ensemble (CESM-LE) and Single-Forcing Large Ensemble, a newly available set of simulations where values of individual forcing agents remain fixed at 1920 levels, allowing for an estimation of their effects. Statistically significant contributions to the DSL FR are identified for greenhouse gases (GHGs) and industrial aerosols (AERs), with particularly strong contributions resulting from AERs in the mid-twentieth century and GHGs in the late twentieth and twenty-first century. Secondary, but important, contributions are identified for biomass burning aerosols in the equatorial Atlantic Ocean in the mid-twentieth century, and for stratospheric ozone in the Southern Ocean during the late twentieth century. Key to understanding regional DSL patterns are ocean heat content and salinity anomalies, which are driven by surface heat and freshwater fluxes, ocean dynamics, and the spatial structure of seawater thermal expansivity. Finally, potential implications for the interpretation of DSL during the satellite era and the longer records from tide gauges are suggested as a topic for future research.

54 ENVIRONMENTAL SCIENCES↗

Effect of Recent Sea Surface Temperature Trends on the Arctic Stratospheric Vortex

The springtime Arctic polar vortex has cooled significantly over the satellite era, with consequences for ozone concentrations in the springtime transition season. The causes of this cooling trend are deduced by using comprehensive chemistry-climate model experiments. Approximately half of the satellite era early springtime cooling trend in the Arctic lower stratosphere was caused by changing sea surface temperatures (SSTs). An ensemble of experiments forced only by changing SSTs is compared to an ensemble of experiments in which both the observed SSTs and chemically- and radiatively-active trace species are changing. By comparing the two ensembles, it is shown that warming of Indian Ocean, North Pacific, and North Atlantic SSTs, and cooling of the tropical Pacific, have strongly contributed to recent polar stratospheric cooling in late winter and early spring, and to a weak polar stratospheric warming in early winter. When concentrations of ozone-depleting substances and greenhouse gases are fixed, polar ozone concentrations show a small but robust decline due to changing SSTs. Ozone changes are magnified in the presence of changing gas concentrations. The stratospheric changes can be understood by examining the tropospheric height and heat flux anomalies generated by the anomalous SSTs. Finally, recent SST changes have contributed to a decrease in the frequency of late winter stratospheric sudden warmings.

Arctic↗

Separating the Influences of Low-Latitude Warming and Sea Ice Loss on Northern Hemisphere Climate Change

Abstract Analyzing a multimodel ensemble of coupled climate model simulations forced with Arctic sea ice loss using a two-parameter pattern-scaling technique to remove the cross-coupling between low- and high-latitude responses, the sensitivity to high-latitude sea ice loss is isolated and contrasted to the sensitivity to low-latitude warming. Despite some differences in experimental design, the Northern Hemisphere near-surface atmospheric sensitivity to sea ice loss is found to be robust across models in the cold season; however, a larger intermodel spread is found at the surface in boreal summer, and in the free tropospheric circulation. In contrast, the sensitivity to low-latitude warming is most robust in the free troposphere and in the warm season, with more intermodel spread in the surface ocean and surface heat flux over the Northern Hemisphere. The robust signals associated with sea ice loss include upward turbulent and longwave heat fluxes where sea ice is lost, warming and freshening of the Arctic Ocean, warming of the eastern North Pacific Ocean relative to the western North Pacific with upward turbulent heat fluxes in the Kuroshio Extension, and salinification of the shallow shelf seas of the Arctic Ocean alongside freshening in the subpolar North Atlantic Ocean. In contrast, the robust signals associated with low-latitude warming include intensified ocean warming and upward latent heat fluxes near the western boundary currents, freshening of the Pacific Ocean, salinification of the North Atlantic, and downward sensible and longwave fluxes over the ocean.

54 ENVIRONMENTAL SCIENCES↗

Examining future changes in coastal low-level jet properties offshore California through dynamical downscaling

The coastal low-level jet, or coastal low-level jet (CLLJ), is a synoptically-forced meteorological feature frequently present offshore the western United States (U.S.). Characterized by a wind speed maximum that resides at the top of the marine boundary layer, the CLLJ is largely controlled by the location and strength of the North Pacific High (NPH) as well as the coastal geometry. Considering the rich wind resource available in this offshore region, the Bureau of Ocean Energy Management identified wind energy lease areas offshore California and supported the deployment of two U.S. Department of Energy wind lidar buoys near Morro Bay and Humboldt. Despite our relatively good understanding of the fundamental mechanisms responsible for large-scale CLLJ properties offshore the western U.S., future changes in CLLJ characteristics are less clear. To address this research challenge, and ultimately to better inform future wind turbine deployments, we use simulations driven by three global climate models (GCMs). We apply self-organizing maps to the model outputs for a historical and two future climate periods to show the range of NPH regimes that support CLLJ conditions during the warm seasons, as well as the subtle contribution from land-falling cyclones approaching the mainland during the cold seasons. Compared to the historical period, the three GCM-driven simulations agree that CLLJ conditions will become more (less) prevalent from central California northward (southward). They agree less with respect to future changes in maximum CLLJ wind speeds and CLLJ heights. However, after considering model biases present during the historical period, wind speeds between the models are actually more similar during the 2070–2095 period than during the historical period. The potential combination of more frequent CLLJ conditions characterized by relatively consistent wind speeds occurring at lower heights across northern California suggests that the Humboldt lease area may be ideal for a long-term wind turbine deployment.

54 ENVIRONMENTAL SCIENCES↗

Significant Climate Benefits from Near-Term Climate Forcer Mitigation in Spite of Aerosol Reductions

Near-term climate forcers (NTCFs), including aerosols and chemically reactive gases such as tropospheric ozone and methane, offer a potential way to mitigate climate change and improve air quality--so called "win-win" mitigation policies. Prior studies support improved air quality under NTCF mitigation, but with conflicting climate impacts that range from a significant reduction in the rate of global warming to only a modest impact. Here, we use state-of-the-art chemistry-climate model simulations conducted as part of the Aerosol and Chemistry Model Intercomparison Project (AerChemMIP) to quantify the 21st-century impact of NTCF reductions, using a realistic future emission scenario with a consistent air quality policy. Non-methane NTCF (NMNTCF; aerosols and ozone precursors) mitigation improves air quality, but leads to significant increases in global mean precipitation of 1.3% by mid-century and 1.4% by end-of-the-century, and corresponding surface warming of 0.23 and 0.21 K. NTCF (all-NTCF; including methane) mitigation further improves air quality, with larger reductions of up to 45% for ozone pollution, while offsetting half of the wetting by mid-century (0.7% increase) and all the wetting by end-of-the-century (non-significant 0.1% increase) and leading to surface cooling of -0.15 K by mid-century and -0.50 K by end-of-the-century. This suggests that methane mitigation offsets warming induced from reductions in NMNTCFs, while also leading to net improvements in air quality.

Near-term climate forcers (NTCFs)↗

Investigation of short-term effective radiative forcing of fire aerosols over North America using nudged hindcast ensembles

Aerosols from fire emissions can potentially have large impact on clouds and radiation. However, fire aerosol sources are often intermittent, and their effect on weather and climate is difficult to quantify. Here we investigated the short-term effective radiative forcing of fire aerosols using the global aerosol–climate model Community Atmosphere Model version 5 (CAM5). Different from previous studies, we used nudged hindcast ensembles to quantify the forcing uncertainty due to the chaotic response to small perturbations in the atmosphere state. Daily mean emissions from three fire inventories were used to consider the uncertainty in emission strength and injection heights. The simulated aerosol optical depth (AOD) and mass concentrations were evaluated against in situ measurements and reanalysis data. Overall, the results show the model has reasonably good predicting skills. Short (10-day) nudged ensemble simulations were then performed with and without fire emissions to estimate the effective radiative forcing. Results show fire aerosols have large effects on both liquid and ice clouds over the two selected regions in April 2009. Ensemble mean results show strong negative shortwave cloud radiative effect (SCRE) over almost the entirety of southern Mexico, with a 10-day regional mean value of −3.0 W m −2 . Over the central US, the SCRE is positive in the north but negative in the south, and the regional mean SCRE is small (−0.56 W m −2 ). For the 10-day average, we found a large ensemble spread of regional mean shortwave cloud radiative effect over southern Mexico (15.6 % of the corresponding ensemble mean) and the central US (64.3 %), despite the regional mean AOD time series being almost indistinguishable during the 10-day period. Moreover, the ensemble spread is much larger when using daily averages instead of 10-day averages. This demonstrates the importance of using a large ensemble of simulations to estimate the short-term aerosol effective radiative forcing.

54 ENVIRONMENTAL SCIENCES↗

A Regime-Oriented Approach to Observationally Constraining Extratropical Shortwave Cloud Feedbacks

The extratropical shortwave (SW) cloud feedback is primarily due to increases in extratropical liquid cloud extent and optical depth. Here, we examine the response of extratropical (35°-75°) marine cloud liquid water path (LWP) to a uniform 4 K increase in sea surface temperature (SST) in global climate models (GCMs) from the fifth Coupled Model Intercomparison Project (CMIP5) and variants of the HadGEM3-GC3.1 GCM. Compositing is used to partition data into periods inside and out of cyclones. The response of extratropical LWP to a uniform SST increase and associated atmospheric response varies substantially among GCMs, but the sensitivity of LWP to cloud controlling factors (CCFs) is qualitatively similar. When all other predictors are held constant, increasing moisture flux drives an increase in LWP. Increasing SST, holding all other predictors fixed, leads to a decrease in LWP. The combinations of these changes lead to LWP, and by extension reflected SW, increasing with warming in both hemispheres. Observations predict an increase in reflected SW over oceans of 0.8 W/sq.m to 1.6 W/sq.m per K SST increase (35°N-75°N) and 1.2 W/sq.m to 1.9 W/sq.m per K SST increase (35°S- 75°S). This increase in reflected SW is mainly due to increased moisture convergence into cyclones because of increasing available moisture. The efficiency at which converging moisture is converted into precipitation determines the amount of liquid cloud. Thus, cyclone precipitation processes are critical to constraining extratropical cloud feedbacks.

extratropical shortwave (SW) cloud feedback↗

CMIP5 Historical Simulations (1850-2012) with GISS ModelE2

Observations of climate change during the CMIP5 extended historical period (1850-2012) are compared to trends simulated by six versions of the NASA Goddard Institute for Space Studies ModelE2 Earth System Model. The six models are constructed from three versions of the ModelE2 atmospheric general circulation model, distinguished by their treatment of atmospheric composition and the aerosol indirect effect, combined with two ocean general circulation models, HYCOM and Russell. Forcings that perturb the model climate during the historical period are described. Five-member ensemble averages from each of the six versions of ModelE2 simulate trends of surface air temperature, atmospheric temperature, sea ice and ocean heat content that are in general agreement with observed trends, although simulated warming is slightly excessive within the past decade. Only simulations that include increasing concentrations of long-lived greenhouse gases match the warming observed during the twentieth century. Differences in twentieth-century warming among the six model versions can be attributed to differences in climate sensitivity, aerosol and ozone forcing, and heat uptake by the deep ocean. Coupled models with HYCOM export less heat to the deep ocean, associated with reduced surface warming in regions of deepwater formation, but greater warming elsewhere at high latitudes along with reduced sea ice. All ensembles show twentieth-century annular trends toward reduced surface pressure at southern high latitudes and a poleward shift of the midlatitude westerlies, consistent with observations.

ocean heat content↗

The Continual Intercomparison of Radiation Codes: Results from Phase I

The computer codes that calculate the energy budget of solar and thermal radiation in Global Climate Models (GCMs), our most advanced tools for predicting climate change, have to be computationally efficient in order to not impose undue computational burden to climate simulations. By using approximations to gain execution speed, these codes sacrifice accuracy compared to more accurate, but also much slower, alternatives. International efforts to evaluate the approximate schemes have taken place in the past, but they have suffered from the drawback that the accurate standards were not validated themselves for performance. The manuscript summarizes the main results of the first phase of an effort called "Continual Intercomparison of Radiation Codes" (CIRC) where the cases chosen to evaluate the approximate models are based on observations and where we have ensured that the accurate models perform well when compared to solar and thermal radiation measurements. The effort is endorsed by international organizations such as the GEWEX Radiation Panel and the International Radiation Commission and has a dedicated website (i.e., http://circ.gsfc.nasa.gov) where interested scientists can freely download data and obtain more information about the effort's modus operandi and objectives. In a paper published in the March 2010 issue of the Bulletin of the American Meteorological Society only a brief overview of CIRC was provided with some sample results. In this paper the analysis of submissions of 11 solar and 13 thermal infrared codes relative to accurate reference calculations obtained by so-called "line-by-line" radiation codes is much more detailed. We demonstrate that, while performance of the approximate codes continues to improve, significant issues still remain to be addressed for satisfactory performance within GCMs. We hope that by identifying and quantifying shortcomings, the paper will help establish performance standards to objectively assess radiation code quality, and will guide the development of future phases of CIRC

Oreopoulos, Lazaros↗

The implications of uncertain renewable resource potentials for global wind and solar electricity projections

Studies exploring long-term energy system transitions rely on resource cost-supply curves derived from estimates of renewable energy (RE) potentials to generate wind and solar power projections. However, estimates of RE potentials are characterized by large uncertainties stemming from methodological assumptions that vary across studies, including factors such as the suitability of land and the performance and configuration of technology. Based on a synthesis of modeling approaches and parameter values used in prior studies, we explore the implications of these uncertain assumptions for onshore wind and solar photovoltaic electricity generation projections globally using the Global Change Analysis Model. We show that variability in parametric assumptions related to land use (e.g. land suitability) are responsible for the most substantial uncertainty in both wind and solar generation projections. Additionally, assumptions about the average turbine installation density and turbine technology are responsible for substantial uncertainty in wind generation projections. Under scenarios that account for climate impacts on wind and solar energy, we find that these parametric uncertainties are far more significant than those emerging from differences in climate models and scenarios in a global assessment, but uncertainty surrounding climate impacts (across models and scenarios) have significant effects regionally, especially for wind. Our analysis suggests the need for studies focusing on long-term energy system transitions to account for this uncertainty.

14 SOLAR ENERGY↗

CICE, The Los Alamos Sea Ice Model

The Los Alamos sea ice model (CICE) is the result of an effort to develop a computationally efficient sea ice component for a fully coupled atmosphere–land–ocean–ice global climate model. It was originally designed to be compatible with the Parallel Ocean Program (POP), an ocean circulation model developed at Los Alamos National Laboratory for use on massively parallel computers. CICE has several interacting components: a vertical thermodynamic model that computes local growth rates of snow and ice due to vertical conductive, radiative and turbulent fluxes, along with snowfall; an elastic-viscous-plastic model of ice dynamics, which predicts the velocity field of the ice pack based on a model of the material strength of the ice; an incremental remapping transport model that describes horizontal advection of the areal concentration, ice and snow volume and other state variables; and a ridging parameterization that transfers ice among thickness categories based on energetic balances and rates of strain. It also includes a biogeochemical model that describes evolution of the ice ecosystem. The CICE sea ice model is used for climate research as one component of complex global earth system models that include atmosphere, land, ocean and biogeochemistry components. It is also used for operational sea ice forecasting in the polar regions and in numerical weather prediction models.

Hunke, Elizabeth↗

CICE, The Los Alamos Sea Ice Model

The Los Alamos sea ice model (CICE) is the result of an effort to develop a computationally efficient sea ice component for a fully coupled atmosphere–land–ocean–ice global climate model. It was originally designed to be compatible with the Parallel Ocean Program (POP), an ocean circulation model developed at Los Alamos National Laboratory for use on massively parallel computers. CICE has several interacting components: a vertical thermodynamic model that computes local growth rates of snow and ice due to vertical conductive, radiative and turbulent fluxes, along with snowfall; an elastic-viscous-plastic model of ice dynamics, which predicts the velocity field of the ice pack based on a model of the material strength of the ice; an incremental remapping transport model that describes horizontal advection of the areal concentration, ice and snow volume and other state variables; and a ridging parameterization that transfers ice among thickness categories based on energetic balances and rates of strain. It also includes a biogeochemical model that describes evolution of the ice ecosystem. The CICE sea ice model is used for climate research as one component of complex global earth system models that include atmosphere, land, ocean and biogeochemistry components. It is also used for operational sea ice forecasting in the polar regions and in numerical weather prediction models.

Hunke, Elizabeth↗