Large-eddy simulations were employed over half-ice and half-water surfaces, with varying surface temperatures, wind speeds, directions, as to test if the atmospheric interaction with the heterogeneous surface can be predicted via a heterogeneity Richardson number. This dataset was used to determine that surface heat fluxes over ice, water, and the aggregate surface seem to be captured reasonably well by the wind direction and the heterogeneity Richardson number, but the mean wind and turbulent kinetic energy (TKE) profiles were not, suggesting that not only the difference in stability between the two surface, but also the individual stabilities over each surface influence the dynamics.
Large-eddy simulations were employed over five different sea ice patterns, with a constant ice fraction, to test if the overlying atmospheric boundary layer (ABL) dynamics and thermodynamics differs. The results of these simulations were used to determine that there were differences in vertical heat flux, momentum flux, and horizontal wind speed, and that more surface information is needed to predict the ABL over the sea ice surface. To see what other surface information is needed, twenty-two landscape metrics were calculated over forty-four different maps at differing resolutions, using the FRAGSTATs program. The results of that analysis are available in a .csv file in this dataset.
Myers, Clayton; Yamada, Masaaki; Ji, Hantao; Yoo, Jongsoo; Jara-Almonte, Jonathan; Fox, William
The loss-of-equilibrium is a solar eruption mechanism whereby a sudden breakdown of the magnetohydrodynamic force balance in the Sun's corona ejects a massive burst of particles and energy into the heliosphere. Predicting a loss-of-equilibrium, which has more recently been formulated as the torus instability, relies on a detailed understanding of the various forces that hold the pre-eruption magnetic flux rope in equilibrium. Traditionally, idealized analytical force expressions are used to derive simplified eruption criteria that can be compared to solar observations and modeling. What is missing, however, is a validation that these idealized analytical force expressions can be applied to the line-tied, low-aspect-ratio conditions of the corona. In this paper, we address this shortcoming by using a laboratory experiment to study the forces that act on long-lived, arched, line-tied magnetic flux ropes. Three key force terms are evaluated over a wide range of experimental conditions: (1) the upward hoop force; (2) the downward strapping force; and (3) the downward toroidal field tension force. First, the laboratory force measurements show that, on average, the three aforementioned force terms cancel to produce a balanced line-tied equilibrium. This finding validates the laboratory force measurement techniques developed here, which were recently used to identify a dynamic toroidal field tension force that can prevent flux rope eruptions [Myers et al., Nature 528, 526 (2015)]. The verification of magnetic force balance also confirms the low-beta assumption that the plasma thermal pressure is negligible in these experiments. Next, the measured force terms are directly compared to their corresponding analytical expressions. While the measured and analytical forces are found to be well correlated, the low-aspect-ratio, line-tied conditions in the experiment are found to both reduce the measured hoop force and increase the measured tension force with respect to analytical expectations. These two co-directed effects combine to generate laboratory flux rope equilibria at lower altitudes than are predicted analytically. Such considerations are expected to modify the loss-of-equilibrium eruption criteria for analogous flux ropes in the solar corona.
Myers, Clayton; Yamada, Masaaki; Ji, Hantao; Yoo, Jongsoo; Jara-Almonte, Jonathan; Fox, William
Solar eruptions are often driven by magnetohydrodynamic instabilities such as the torus and kink instabilities that act on line-tied magnetic flux ropes. Recent laboratory experiments designed to study these eruptive instabilities have demonstrated the key role of both dynamic (Myers et al 2015 Nature 528, 526) and quasi-static (Myers et al 2016 Phys. Plasmas, in press) magnetic tension forces in contributing to the equilibrium and stability of line-tied magnetic flux ropes. In this paper, we synthesize these laboratory results and explore the relationship between the dynamic and quasi-static tension forces. While the quasi-static tension force is found to contribute to the flux rope equilibrium in a number of regimes, the dynamic tension force is substantial mostly in the so-called failed torus regime where magnetic self-organization events prevent the flux rope from erupting.
One aspect of the interaction between fast ions and magnetohydrodynamic (MHD) instabilities is the fast ion transport. Coupled kink and tearing MHD instabilities have also been reported to cause fast ion transport. Recently, the ''kick" model has been developed to compute the evolution of the fast ion distribution from the neutral beam injection using instabilities as phase space resonance sources. The goal of this paper is to utilize the kick model to understand the physics of fast ion transport caused by the coupled kink and tearing modes. Soft X-ray diagnostics are used to identify the mode parameters in NSTX. The comparison of neutron rates measured and computed from time-dependent TRANSP simulation with the kick model shows the coupling of kink and tearing mode is important in determination of the fast ion transport. The numerical scan of the mode parameters shows that the relative phase of the kink and tearing modes and the overlapping of kink and tearing mode resonances in the phase space can affect the fast ion transport, suggesting that the synergy of the coupled modes may be causing the fast ion transpor
A compact and multi-view Solid State Neutral Particle Analyzer (SSNPA) diagnostic based on silicon photodiode arrays has been successfully tested on the National Spherical Torus Experiment-Upgrade (NSTX-U). The SSNPA diagnostic provides spatially, temporally, and pitch-angle resolved measurements of fast-ion distribution by detecting fast neutral flux resulting from charge exchange (CX) reactions. The system consists of three 16-channel subsystems: t-SSNPA viewing the plasma mid-radius and neutral beam (NB) line #2 tangentially, r-SSNPA viewing the plasma core and NB line #1 radially and p-SSNPA with no intersection with any NB lines. Due to the setup geometry, the active CX signals of t-SSNPA and r-SSNPA are mainly sensitive to passing and trapped particles respectively. In addition, both t-SSNPA and r-SSNPA utilize three vertically stacked arrays with different filter thickness to obtain coarse energy information. The experimental data show that all channels are operational. The signal to noise ratio is typically larger than 10 and the main noise is x-ray induced signal. The active and passive CX signals are clearly observed on t-SSNPA and r-SSNPA during NB modulation. The SSNPA data also indicate significant losses of passing particles during sawteeth, while trapped particles are weakly affected. Fluctuations up to 120 kHz, have been observed on SSNPA, and they are strongly correlated with magnetohydrodynamics (MHD) instabilities.
The growth of magnetic islands in NSTX is modeled successfully, with the consideration of passing fast ions. It is shown that a good quantitative agreement between simulation and experimental measurement can be achieved when the uncompensated cross-field current induced by passing fast ions is included in the island growth model. The fast ion parameters,
along with other equilibrium parameters, are obtained self-consistently using the TRANSP code with the assumptions of the ‘kick’ model (Podestà et al 2017 Plasma Phys. Control. Fusion 59 095008). The results show that fast ions can contribute to overcoming the stabilizing effect of polarization current for magnetic island growth.
Particle distribution functions evolving under the Lorentz operator can be simulated with the Langevin equation for pitch angle scattering. This approach is frequently used in particle based Monte-Carlo simulations of plasma collisions, among others. However, most numerical treatments do not guarantee energy conservation, which may lead to unphysical artifacts such as numerical heating and spectra distortions. We present a novel structure-preserving numerical algorithm for the Langevin equation for pitch angle scattering. Similar to the well-known Boris algorithm, the proposed numerical scheme takes advantage of the structure-preserving properties of the Cayley transform when calculating the velocity-space rotations. The resulting algorithm is explicitly solvable, while preserving the norm of velocities down to machine precision. We demonstrate that the method has the same order of numerical convergence as the traditional stochastic Euler-Maruyama method.
Explosive volcanic eruptions have large climate impacts, and can serve as observable tests of the climatic response to radiative forcing. Using a high resolution climate model, we contrast the climate responses to Pinatubo, with symmetric forcing, and those to Santa Maria and Agung, which had meridionally asymmetric forcing. Although Pinatubo had larger global-mean forcing, asymmetric forcing strongly shifts the latitude of tropical rainfall features, leading to larger local precipitation/TC changes. For example, North Atlantic TC activity over is enhanced/reduced by SH-forcing (Agung)/NH-forcing (Santa Maria), but changes little in response to the Pinatubo forcing. Moreover, the transient climate sensitivity estimated from the response to Santa Maria is 20% larger than that from Pinatubo or Agung. This spread in climatic impacts of volcanoes needs to be considered when evaluating the role of volcanoes in global and regional climate, and serves to contextualize the well-observed response to Pinatubo.
Khanna, Jaya; Medvigy, David; Fueglistaler, Stephan; Walko, Robert
More than 20% Amazon rainforest has been cleared in the past three decades triggering important hydroclimatic changes. Small-scale (~few kilometers) deforestation in the 1980s has caused thermally-triggered atmospheric circulations that increase regional cloudiness and precipitation frequency. However, these circulations are predicted to diminish as deforestation increases. Here we use multi-decadal satellite records and numerical model simulations to show a regime shift in the regional hydroclimate accompanying increasing deforestation in Rondônia, Brazil. Compared to the 1980s, present-day deforested areas in downwind western Rondônia are found to be wetter than upwind eastern deforested areas during the local dry season. The resultant precipitation change in the two regions is approximately ±25% of the deforested area mean. Meso-resolution simulations robustly reproduce this transition when forced with increasing deforestation alone, showing a negligible role of large-scale climate variability. Furthermore, deforestation-induced surface roughness reduction is found to play an essential role in the present-day dry season hydroclimate. Our study illustrates the strong scale-sensitivity of the climatic response to Amazonian deforestation and suggests that deforestation is sufficiently advanced to have caused a shift from a thermally- to a dynamically-driven hydroclimatic regime.
Li, Zhongshu; Gallagher, Kevin P.; Mauzerall, Denise L.
The dataset include a list of power projects outside of China that receive Chinese foreign direct investment from 2000 to 2018. Detailed information including project capacity, location, share of Chinese ownership, type of power generating technologies are collected for each power project.
Muniz, Maria Carolina; Gartner III, Thomas E.; Riera, Marc; Knight, Christopher; Yue, Shuwen; Paesani, Francesco; Panagiotopoulos, Athanassios Z.
This dataset contains all data (including input files, simulation trajectories as well as other data files and analysis scripts) related to the publication "Vapor-liquid equilibrium of water with the MB-pol many-body potential" by Muniz et al. in preparation (2021). In this work, we assessed the performance of the MB-pol many-body potential with respect to water's vapor-liquid equilibrium properties. Through the use of direct coexistence molecular dynamics, we calculated properties such as coexistence densities, surface tension, vapor pressures and enthalpy of vaporization. We found that MB-pol is able to predict these properties in good agreement with experimental data. The results attest to the chemical accuracy of MB-pol and its large range of application across water's phase diagram.
This dataset contains all data relevant to a forthcoming publication in which we used molecular simulation methods to study the phase behavior of supercooled water. The dataset contains simulation input and output files, processed data files, and image files used to create all plots in the manuscript. Python analysis scripts are also included, including instructions for how to re-generate all plots in the manuscript.