Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “CESM2”

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 73 records · Page 4

Simulations With the Marine Biogeochemistry Library (MARBL)

Abstract The Marine Biogeochemistry Library (MARBL) is a prognostic ocean biogeochemistry model that simulates marine ecosystem dynamics and the coupled cycles of carbon, nitrogen, phosphorus, iron, silicon, and oxygen. MARBL is a component of the Community Earth System Model (CESM); it supports flexible ecosystem configuration of multiple phytoplankton and zooplankton functional types; it is also portable, designed to interface with multiple ocean circulation models. Here, we present scientific documentation of MARBL, describe its configuration in CESM2 experiments included in the Coupled Model Intercomparison Project version 6 (CMIP6), and evaluate its performance against a number of observational data sets. The model simulates present‐day air‐sea CO 2 flux and many aspects of the carbon cycle in good agreement with observations. However, the simulated integrated uptake of anthropogenic CO 2 is weak, which we link to poor thermocline ventilation, a feature evident in simulated chlorofluorocarbon distributions. This also contributes to larger‐than‐observed oxygen minimum zones. Moreover, radiocarbon distributions show that the simulated circulation in the deep North Pacific is extremely sluggish, yielding extensive oxygen depletion and nutrient trapping at depth. Surface macronutrient biases are generally positive at low latitudes and negative at high latitudes. CESM2 simulates globally integrated net primary production (NPP) of 48 Pg C yr −1 and particulate export flux at 100 m of 7.1 Pg C yr −1 . The impacts of climate change include an increase in globally integrated NPP, but substantial declines in the North Atlantic. Particulate export is projected to decline globally, attributable to decreasing export efficiency associated with changes in phytoplankton community composition.

Meteorology & Atmospheric Sciences↗

Standardized Data to Improve Understanding and Modeling of Soil Nitrogen at Continental Scale

Abstract Nitrogen (N) is a key limiting nutrient in terrestrial ecosystems, but there remain critical gaps in our ability to predict and model controls on soil N cycling. This may be in part due to lack of standardized sampling across broad spatial–temporal scales. Here, we introduce a continentally distributed, publicly available data set collected by the National Ecological Observatory Network (NEON) that can help fill these gaps. First, we detail the sampling design and methods used to collect and analyze soil inorganic N pool and net flux rate data from 47 terrestrial sites. We address methodological challenges in generating a standardized data set, even for a network using uniform protocols. Then, we evaluate sources of variation within the sampling design and compare measured net N mineralization to simulated fluxes from the Community Earth System Model 2 (CESM2). We observed wide spatiotemporal variation in inorganic N pool sizes and net transformation rates. Site explained the most variation in NEON’s stratified sampling design, followed by plots within sites. Organic horizons had larger pools and net N transformation rates than mineral horizons on a sample weight basis. The majority of sites showed some degree of seasonality in N dynamics, but overall these temporal patterns were not matched by CESM2, leading to poor correspondence between observed and modeled data. Looking forward, these data can reveal new insights into controls on soil N cycling, especially in the context of other environmental data sets provided by NEON, and should be leveraged to improve predictive modeling of the soil N cycle.

54 ENVIRONMENTAL SCIENCES↗

Anthropogenic Aerosols Contribute to the Recent Decline in Precipitation Over the U.S. Southwest

Abstract The winter‐spring precipitation over the Southwestern United States (SWUS) decreased since 1980. It is frequently attributed to Pacific internal decadal variability, but recent studies found anthropogenic aerosols (AA) can also induce a transition to a negative Pacific Decadal Variability (PDV) phase. We revisit the attribution of SWUS drying by quantifying the contributions of anthropogenically forced decadal Pacific Sea Surface Temperatures (SSTs). Applying a low‐frequency component analysis to observations, Community Earth System Model version 2 (CESM2) all‐forcings and single‐forcing large ensembles, we find up to 42% of the observed precipitation trend to be related to the AA‐induced negative PDV‐like pattern, which is driven by the emission shift from the Western to the Eastern Hemisphere. In CESM2, other radiative forcings counteract the influence of AA, but it remains unclear whether the model correctly simulates this balance. This implies that the near‐future trajectories of these forcings, in particular Asian aerosols, are important for projections of SWUS precipitation.

Geology↗

Investigating Future Arctic Sea Ice Loss and Near–Surface Wind Speed Changes Related to Surface Roughness Using the Community Earth System Model

The Arctic is undergoing a pronounced and rapid transformation in response to changing greenhouse gasses, including reduction in sea ice extent and thickness. There are also projected increases in near-surface Arctic wind. This study addresses how the winds trends may be driven by changing surface roughness and/or stability in different Arctic regions and seasons, something that has not yet been thoroughly investigated. We analyze 50 experiments from the Community Earth System Model Version 2 (CESM2) Large Ensemble and five experiments using CESM2 with an artificially decreased sea ice roughness to match that of the open ocean. We find that with a smoother surface there are higher mean wind speeds and slower mean ice speeds in the autumn, winter, and spring. The artificially reduced surface roughness also strongly impacts the wind speed trends in autumn and winter, and we find that atmospheric stability changes are also important contributors to driving wind trends in both experiments. In contrast to the clear impacts on winds, the sea ice mean state and trends are statistically indistinguishable, suggesting that near-surface winds are not major drivers of Arctic sea ice loss. Two major results of this work are: (a) the near-surface wind trends are driven by changes in both surface roughness and near-surface atmospheric stability that are themselves changing from sea ice loss, and (b) the sea ice mean state and trends are driven by the overall warming trend due to increasing greenhouse gas emissions and not significantly impacted by coupled feedbacks with the surface winds.

54 ENVIRONMENTAL SCIENCES↗

Examining Tropical Convection Features at Storm-Resolving Scales Over the Maritime Continent Region

Global Storm Resolving Models (GSRMs) provide a way to understand weather and climate events across scales for better-informed climate impacts. In this work, we apply the recently developed and validated CAM (Community Atmosphere Model)—MPAS (Model for Prediction Across Scales) modeling framework, based on the open-source Community Earth System Model (CESM2), to examine the tropical convection features at the storm resolving scale over the Maritime Continent region at 3 km horizontal spacing. We target two global numerical experiments during the winter season of 2018 for comparison with observation in the region. Here we focus on the investigation of the representations of the convective systems, precipitation statistics, and tropical cyclone behaviors. We found that regional-refined experiments show more accurate precipitation distributions, diurnal cycles, and better agreement with observations for tropical cyclone features in terms of intensity and strength statistics. We expect the exploration of this work will further advance the development and use of the storm-resolving model in precipitation predictions across scales.

54 ENVIRONMENTAL SCIENCES↗

Interannual fires as a source for subarctic summer decadal climate variability mediated by permafrost thawing

Climate model simulations run under the Coupled Model Intercomparison Project Phase 6 (CMIP6) use an inhomogeneous biomass burning aerosol (BBA) emission dataset, which exhibits pronounced interannual variability from 1997–2014 due to the infusion of satellite data. Using the Community Earth System Model version 2 Large Ensemble (CESM2-LE) with original and smoothed CMIP6 BBA forcings, we show that the CMIP6 data inhomogeneity causes spurious decadal subarctic land warming. During years with reduced aerosol concentrations, increased solar radiation can trigger abrupt subarctic permafrost thawing, increased soil water drainage, upper soil drying, and subsequent surface warming. This slow process, which is further amplified by nonlinear cloud-aerosol interactions, cannot be completely offset during years of increased aerosol fluxes, thereby reddening surface temperature spectra in response to large-amplitude interannual aerosol forcing. More generally, our CESM2 experiments identify a pathway for generating decadal variability in high latitudes, involving interannual shortwave forcing and slow nonlinear soil responses.

54 ENVIRONMENTAL SCIENCES↗

A Positive Low Cloud–Sea Surface Temperature Feedback in the East Asian Marginal Seas during El Niño Mature Winters and Their Following Spring

During the winter and subsequent spring of an El Niño year, the East Asian marginal sea (EAMS) exhibits positive sea surface temperature anomalies (SSTAs) and fewer low clouds, while the western North Pacific experiences negative SSTAs. In this study, we suggest that the positive SSTAs in EAMS are maintained by a positive low cloud–SST feedback. In neutral winters and springs, the EAMS is covered by low clouds, which have a cooling effect on surface temperatures. During an El Niño year, a western North Pacific anomalous anticyclone is established, and along its northwestern flank, there are favorable conditions for convergence of moisture and weaker surface latent heat flux over the EAMS. Once a positive SSTA has been established, a further reduction of turbulent mixing results in less low cloud and enhanced solar heating of the ocean mixed layer; this reinforces and maintains both the positive SSTA and the lack of low cloud via a positive feedback mechanism. Further, the concurrent increase of low cloud–SST feedback and anticyclone circulation strengths is evident in the coupled-model simulations from phase 6 of the Coupled Model Intercomparison Project. Furthermore, sensitivity experiments, performed with the atmospheric components of Community Earth System Model (CESM2), reveal that a positive SSTA helps to maintain the western North Pacific anomalous anticyclone. Four pacemaker-coupled experiments by CESM2, with sea surface temperature in the equatorial Pacific restored to the observational anomalies plus the model climatology and altered low cloud feedback over EAMS, suggest that the low cloud–SST feedback results in more than the maintenance of a positive SSTA over the EAMS: the positive feedback is also a previously overlooked mechanism for the maintenance of the western North Pacific anomalous anticyclone.

54 ENVIRONMENTAL SCIENCES↗

Poleward Migration of the Latitude of Maximum Tropical Cyclone Intensity—Forced or Natural?

Abstract Past studies have shown a significant observed poleward trend in the latitude at which tropical cyclones reach their lifetime maximum intensity (LMI), especially in the northwest Pacific basin. Given the brevity of the historical record, it remains difficult to separate the forced trend from internal variability of the climate system. A recently developed tropical cyclone downscaling model is used to downscale the Community Earth System Model, version 2 (CESM2), preindustrial control simulation. It is found that the observed trend in the latitude at which tropical cyclones reach their LMI in the northwest Pacific is very unlikely to be caused by internal variability. The same downscaling model is then used to downscale CESM2 simulations under historical forcing. The resulting trend distribution shows a significant poleward migration of tropical cyclone LMI even after regressing out both natural variability and the part of the forced warming pattern that projects onto natural variability. The results indicate that the observed poleward migration of the latitude at which tropical cyclones reach their LMI in the northwest Pacific basin is likely to be, at least in part, forced. However, the magnitude of the projected poleward trend in climate models can be significantly modulated by the simulated spatial pattern of ocean warming. This highlights how discrepancies between models and observations, with regard to projected changes to the equatorial zonal sea surface temperature gradient under anthropogenic forcing, can lead to large uncertainties in projected changes to the LMI latitude of tropical cyclones. Significance Statement Observations in the northwest Pacific basin show that the latitude at which tropical cyclones are at their most intense has been trending northward in the recent half century. These changes are important since tropical cyclones could bring hazardous weather to coastal areas that are poorly equipped to handle them. Here, we show that natural variations in Earth’s climate are very unlikely to explain the observed poleward trend in the latitude that tropical cyclone reach their maximum intensity. We find that it is much more likely that the observed trend is forced by human-related emissions, though the spatial pattern of warming in response to greenhouse emissions can have significant impacts on the magnitude of the trend.

Lin, Jonathan↗

How Well Can CMIP6 Models Represent the Observed Influence of the Pacific and Indian Oceans on the Indian Summer Monsoon Rainfall?

This study evaluates the ability of CMIP6 climate models to simulate the observed effects of tropical Pacific and Indian Ocean sea surface temperature anomalies (SSTAs) on Indian summer monsoon rainfall (ISMR) variability. Using observational data and the large ensemble historical simulations of seven CMIP6 models from 1950 to 2014, we applied a cyclostationary linear inverse model (CS-LIM) to isolate the impacts of tropical Pacific SSTAs, Indian Ocean SSTAs and their interaction on the interannual variability of ISMR. Overall, CMIP6 models well reproduced the observed enhanced (reduced) ISMR variability from Pacific SSTAs (Indian Ocean SSTAs and the Indo-Pacific interaction), but with varying spatial patterns and magnitudes. While CESM2 and E3SM-2-0 showed the best agreement with observations for the effects of Pacific SSTAs and the Indo-Pacific interaction, respectively, CMIP6 models showed mixed results for the impacts from Indian Ocean SSTAs. Composite analysis of ISMR anomalies during the developing phases of pure and co-occurring El Niño-Southern Oscillation (ENSO) and Indian Ocean dipole (IOD) events revealed that the impacts from Pacific SSTAs were captured reasonably well by E3SM-2-0, CESM2, MIROC6, and MPI-ESM1-2-LR, while E3SM-2-0 also showed the best agreement with observations for the effects from the Indo-Pacific interaction. However, all models showed substantial biases in simulating the Indian Ocean SSTA impacts on ISMR, especially for pure El Niño events. Overall, this study provides new insights into how individual CMIP6 models simulate the isolated impacts from the tropical Pacific and Indian Oceans, which has important applications for improving ISMR predictions and interpreting ISMR future projections.

monsoon↗

Primary and secondary ice production: interactions and their relative importance

A discrepancy of up to 5 orders of magnitude between ice crystal and ice nucleating particle (INP) number concentrations was found in the measurements, indicating the potentially important role of secondary ice production (SIP) in the clouds. However, the interactions between primary and SIP processes and their relative importance remain unexplored. In this study, we implemented five different ice nucleation schemes as well as physical representations of SIP processes (i.e., droplet shattering during rain freezing, ice-ice collisional break-up, and rime splintering) in the Community Earth System Model version 2 (CESM2). We ran CESM2 in the single column mode for model comparisons with the DOE Atmospheric Radiation Measurement (ARM) Mixed-Phase Arctic Cloud Experiment (M-PACE) observations. We found that the model experiments with aerosol-aware ice nucleation schemes and SIP processes yield the best simulation results for the M-PACE single-layer mixed-phase clouds. We further investigated the relative importance of ice nucleation and SIP to ice number and cloud phase as well as interactions between ice nucleation and SIP in the M-PACE single-layer mixed-phase clouds. Our results show that SIP contributes 80 % to the total ice formation and transforms ~30 % of pure liquid-phase clouds simulated in the model experiments without considering SIP into mixed-phase clouds. The SIP is not only a result of ice crystals produced from ice nucleation, but also competes with the ice nucleation by reducing the number concentrations of cloud droplets and cloud-borne dust INPs. Conversely, strong ice nucleation also suppresses SIP by glaciating mixed-phase clouds and thereby reducing the amount of precipitation particles (rain and graupel).

54 ENVIRONMENTAL SCIENCES↗

The Seasonal-to-Multiyear Large Ensemble (SMYLE) prediction system using the Community Earth System Model version 2

Abstract. The potential for multiyear prediction of impactful Earth system change remains relatively underexplored compared to shorter (subseasonal to seasonal) and longer (decadal) timescales. In this study, we introduce a new initialized prediction system using the Community Earth System Model version 2 (CESM2) that is specifically designed to probe potential and actual prediction skill at lead times ranging from 1 month out to 2 years. The Seasonal-to-Multiyear Large Ensemble (SMYLE) consists of a collection of 2-year-long hindcast simulations, with four initializations per year from 1970 to 2019 and an ensemble size of 20. A full suite of output is available for exploring near-term predictability of all Earth system components represented in CESM2. We show that SMYLE skill for El Niño–Southern Oscillation is competitive with other prominent seasonal prediction systems, with correlations exceeding 0.5 beyond a lead time of 12 months. A broad overview of prediction skill reveals varying degrees of potential for useful multiyear predictions of seasonal anomalies in the atmosphere, ocean, land, and sea ice. The SMYLE dataset, experimental design, model, initial conditions, and associated analysis tools are all publicly available, providing a foundation for research on multiyear prediction of environmental change by the wider community.

54 ENVIRONMENTAL SCIENCES↗

A protocol for model intercomparison of impacts of marine cloud brightening climate intervention

A modeling protocol (defined by a series of climate model simulations with specified model output) is introduced. Studies using these simulations are designed to improve the understanding of climate impacts using a strategy for climate intervention (CI) known as marine cloud brightening (MCB) in specific regions; therefore, the protocol is called MCB-REG (where REG stands for region). The model simulations are not intended to assess consequences of a realistic MCB deployment intended to achieve specific climate targets but instead to expose responses to interventions in six regions with pervasive cloud systems that are often considered candidates for such a deployment. A calibration step involving simulations with fixed sea surface temperatures (SSTs) is first used to identify a common forcing, and then coupled simulations with forcing in individual regions and combinations of regions are used to examine climate impacts. Synthetic estimates constructed by superposing responses from simulations with forcing in individual regions are considered a means of approximating the climate impacts produced when MCB interventions are introduced in multiple regions. A few results comparing simulations from three modern climate models (CESM2, E3SMv2, and UKESM1) are used to illustrate the similarities and differences between model behavior and the utility of estimates of MCB climate responses that were synthesized by summing responses introduced in individual regions. Cloud responses to aerosol injections differ substantially between models (CESM2 clouds appear much more susceptible to aerosol emissions than the other models), but patterns in precipitation and surface temperature responses were similar when forcing is imposed with similar amplitudes in the same regions. A previously identified La Niña-like response to forcing introduced in the Southeast Pacific is evident in this study, but the amplitude of the response was shown to markedly differ across the three models. Other common response patterns were also found and are discussed. Forcing in the Southeast Atlantic consistently (across all three models) produces weaker global cooling than that in other regions, and the Southeast Pacific and South Pacific show the strongest cooling. This indicates that the efficiency of a given intervention depends on not only the susceptibility of the clouds to aerosol perturbations, but also the strength of the underlying radiative feedbacks and ocean responses operating within each region. These responses were generally robust across models, but more studies and an examination of responses with ensembles would be beneficial.

54 ENVIRONMENTAL SCIENCES↗

The short-term comprehensive impact of the phase-out of global coal combustion on air pollution and climate change

With the continuous intensification of global warming, the reduction and ultimate phase-out of coal combustion is an inevitable trend in the future global energy transformation. This study comprehensively analyzed the impact of phasing out coal combustion on global emissions and concentrations of air pollutants, radiative fluxes, meteorology and climate using Community Earth System Model 2 (CESM2). The results indicate that after the global phase-out of coal combustion, there is a marked decrease in the concentrations of sulfur dioxide (SO 2 ), nitrogen oxides (NO x ) and fine particulate matter (PM 2.5 ), with some regions experiencing a reduction of exceeding 50%. There is no significant change in global ozone (O 3 ) concentration. There are decreasing AOD and positive radiative fluxes globally in the short term, though the cloud contributes minor negative radiative fluxes. The global air temperature may increase by approximately (0.02 ± 0.15) °C on average with regional and seasonal variations, and the precipitation may potentially increase by approximately (2.7 ± 40.6) mm yr −1 globally and over 20% in equatorial regions in the short term. But combined with the decreasing trend of cloud water content in the Northern Hemisphere, it indicates a potential increase in the extremity of precipitation events. In conclusion, this study provides references for global control of air pollution, mitigation strategies of climate change, and transformation of energy structures under the objective of “carbon neutrality”, such as focusing on the negative climate impacts of exacerbating regional warming and increasing extreme precipitation resulting from the rapid reduction of aerosols in the short term.

01 COAL, LIGNITE, AND PEAT↗

Attribution of Snowpack Errors to Simulated Temperature and Precipitation in E3SMv1 Over the Contiguous United States

Abstract Snow water equivalent (SWE), temperature, and precipitation biases and trends are evaluated in the atmosphere‐land simulations of the Energy Exascale Earth System Model version 1 (E3SMv1) in comparison to the Community Earth System Model version 2 (CESM2) and two other models using the ground measurement‐based University of Arizona (UA) snow product. SWE, temperature, and precipitation biases are highest in magnitude in the Western contiguous United States (CONUS). SWE errors are attributed to temperature and precipitation through multiple linear regressions of normalized errors, the coefficients of which represent the sensitivities to temperature and precipitation errors. SWE errors are more sensitive to temperature errors throughout the CONUS. Model SWE and temperature trends are generally opposite from UA product trends in the Western CONUS. SWE trend errors are also attributed to temperature and precipitation trend errors using multiple linear regressions of normalized trend errors. SWE trend errors are more sensitive to those of precipitation at higher elevations (>1,500 m) in the Western CONUS in these simulations. Thus, the sensitivity to temperature and precipitation differ for SWE errors and its trend errors. Furthermore, the SWE trend errors are more sensitive to temperature and precipitation in the atmosphere‐ocean coupled simulations in which the atmosphere‐land is coupled to an active ocean model. These results suggest that both errors in simulated temperature and precipitation contribute to SWE errors.

Meteorology & Atmospheric Sciences↗

Can Biomass Burning Aerosol Induced Surface Cooling Be Amplified Through Sea Surface Temperature‐Cloud Feedback Over the Southeast Atlantic?

Biomass burning (BB) aerosols exert a strong surface cooling effect over the southeast Atlantic (SEA) via aerosol-radiation and aerosol-cloud interactions. The reduction of the sea surface temperature (SST) can trigger the SST-low cloud feedback. Whether this feedback can amplify the surface cooling effect is examined. The modeling results from the Community Earth System Model version 2 (CESM2) demonstrate that counterintuitively the cloud radiative effect (CRE) caused by the BB aerosols is weaker if SST-low cloud feedback is considered compared to fixed-SST simulation (-2.99 W m -2 vs. -4.79 W m -2 ). This is caused by (a) stronger sea breeze due to larger sea-land temperature contrast causing less smoke transport over SEA and (b) less moisture supply from surface due to colder SST. Changes in SST also lead to counterclockwise rotation of ocean circulation anomalies. Consequently, the excess heat transport from the equator reverses the direction of SST-cloud feedback in this region.

0305↗

Potential Non-Linearities in the High Latitude Circulation and Ozone Response to Stratospheric Aerosol Injection

The impacts of Stratospheric Aerosol Injection (SAI) on the atmosphere and surface climate depend on when and where the sulfate aerosol precursors are injected, as well as on how much surface cooling is to be achieved. We use a set of CESM2(WACCM6) SAI simulations achieving three different levels of global mean surface cooling and demonstrate that unlike some direct surface climate impacts driven by the reflection of solar radiation by sulfate aerosols, the SAI-induced changes in the high latitude circulation and ozone are more complex and could be non-linear. This manifests in our simulations by disproportionally larger Antarctic springtime ozone loss, significantly larger intra-ensemble spread of the Arctic stratospheric jet and ozone responses, and non-linear impacts on the extratropical modes of surface climate variability under the strongest-cooling SAI scenario compared to the weakest one. These potential non-linearities may add to uncertainties in projections of regional surface impacts under SAI.

54 ENVIRONMENTAL SCIENCES↗

Future Anthropogenic Land Use Change Impacts on Carbonaceous Aerosol and Implications for Climate and Air Quality

Future anthropogenic land use change (LUC) may alter atmospheric carbonaceous aerosol (black carbon and organic aerosol) burden by perturbing biogenic and fire emissions. However, there has been little investigation of this effect. We examine the global evolution of future carbonaceous aerosol under the Shared Socioeconomic Pathways projected reforestation and deforestation scenarios using the CESM2 model from present‐day to 2100. Compared to present‐day, the change in future biogenic volatile organic compounds emission follows changes in forest coverage, while fire emissions decrease in both projections, driven by trends in deforestation fires. The associated carbonaceous aerosol burden change produces moderate aerosol direct radiative forcing (−0.021 to +0.034 W/m 2 ) and modest mean reduction in PM 2.5 exposure (−0.11 μg/m 3 to −0.23 μg/m 3 ) in both scenarios. We find that future anthropogenic LUC may be more important in determining atmospheric carbonaceous aerosol burden than direct anthropogenic emissions, highlighting the importance of further constraining the impact of LUC.>

anthropogenic land use change↗

Mechanisms of Projected Changes in Thunderstorm Downburst Environments Across the United States

Responses of downdraft convective available potential energy (DCAPE) to global warming were investigated using the Community Earth System Model (CESM2) under a high‐emission scenario through the year 2100. DCAPE is projected to increase by 5%–12% on average in most areas, independently of wind shear. A diagnostic of downdraft buoyancy is introduced to understand the mechanisms of DCAPE responses. Much of the increase in mean DCAPE is temperature‐driven, with additional contributions from changes in relative humidity and downdraft origin heights. However, extreme values increase at much faster rates than can be explained by local warming. In winter, the latitude of significant DCAPE and CAPE shifts poleward by more than 5° due to larger changes in downburst environments within midlatitude cyclones. The projected increase in cold‐season extremes indicates an interaction between weather events and warming trends that increases the potential for downbursts and straight‐line winds in winter.

Williams, Ian N. [Iowa State Univ., Ames, IA (Unit↗