Dynamic Optimization of Nitrogen Use in Agriculture

Agricultural production is highly dependent on inorganic substances including fertilizers. High-yielding crop varieties, such as corn, require large amounts of primary nutrients including nitrogen, phosphorus and potassium. Farmers often add a surplus of nutrients to crops to maximize yields. Utilization of primary nutrients has increased by more than 300% while that of nitrogen alone has increased by more than 600% between 1960 and 2007 (USDA, 2009). From 1964 to 2007, the use of nitrogen in the corn sector alone increased from 1,623,000 to 5,714,000 nutrient tons (USDA, 2009). While increasing production, increased fertilizer use can potentially create negative externalities in the form of nitrate-nitrogen contamination in groundwater. Groundwater is the source of drinking water for about half the total U.S. population and nearly all of the rural population, and it provides over 50 billion gallons per day for agricultural needs (USGS, 2009). In the U.S. the main source of nitrate pollution in the groundwater results from the actions of farmers through the use of fertilizers and other chemicals (Haller, et al. 2009). Nitrogen-nitrate contamination can have adverse human affects including methemoglobinemia or ―blue-baby‖ syndrome (Majumdar, 2003). The potential for nitrate contamination in corn production is especially problematic as corn alone accounts for over 90% of feed grains produced in the U.S. (USDA, 2009). The USDA estimates that approximately 80 million acres of land is planted to corn, with the majority in the Heartland region (the Midwest) of the U.S. (2009). The Heartland region is primarily rural and much of the population there derives its drinking water from groundwater. Therefore, the potential for groundwater contamination is greatly increased in this region.Environmental Economics, Nitrogen/Nitrate Contamination, Dynamic Optimization, Agriculture, Agricultural and Food Policy, Demand and Price Analysis, Environmental Economics and Policy, C61, C63, Q10, Q51, Q53,

Dynamic Optimization of Nitrogen Use in Agriculture

Agricultural and Food Policy, Crop Production/Industries,

Tectonics and volcanism of Eastern Aphrodite Terra: No subduction, no spreading

Eastern Aphrodite Terra is approximately equal in size to the western North American Cordillera, from Mexico to Alaska. Its size and unique landforms make it an important area for understanding the tectonics of Venus, yet models for its formation are diametrically opposed. This region is part of the Equatorial Highlands, which was proposed as a region of lithospheric thinning, isostatic uplift, and attendant volcanism. Eastern Aphrodite Terra is dominated by circular structures within which deformation and volcanism are intimately related. These structures are marked by radial and concentric fractures, and volcanic flows that emanate from a central vent, as well as from concentric fracture sets. Cross-cutting relations between flows and concentric fracture sets indicate that outer concentric fracture sets are younger than inner fracture sets. The circular structures are joined by regional northeast- to east-trending fractures that dominantly postdate formation of the circular structures. We propose that the circular structures 'grow' outward with time. Although these structures probably represent addition of crust to the lithosphere, they do not represent significant lithospheric spreading or convergence, and the region does not mark the boundary between two distinct tectonic plates. This region is not easily explained by analogy with either terrestrial midocean rifts or subduction zones. It is perhaps best explained by upwelling of magma diapirs that blister the surface, but do not cause significant lithospheric spreading. Further study of the structural and volcanic evolution of this region using Magellan altimetry and SAR data should lead to better understanding of the tectonic evolution of this region

Dynamics of rapidly spinning blob-filaments: fluid theory with a parallel kinetic extension

Blob-filaments (or simply 'blobs') are coherent structures formed by turbulence and sustained by nonlinear processes in the edge and scrape-off layer (SOL) of tokamaks and other magnetically confined plasmas. The dynamics of these blob-filaments, in particular their radial motion, can influence the scrape-off layer width and plasma interactions with both the divertor target and with the main chamber walls. Motivated by recent results from the XGC1 gyrokinetic simulation code reported on elsewhere [J. Cheng et al. submitted to Nucl. Fusion and available at arXiv:2302.02877v1], a theory of rapidly spinning blob-filaments has been developed. The theory treats blob filaments in the closed flux surface region or the region that is disconnected from sheaths in the SOL. It extends previous work by treating blob spin, arising from partially or fully adiabatic electrons, as the leading order effect and retaining inertial (ion charge polarization) physics in next order. Spin helps to maintain blob coherency and affects the blob's propagation speed. Dipole charge polarization, treated perturbatively, gives rise to blob-filaments with relatively slow radial velocity, comparable to that observed in the simulations. The theory also treats the interaction of rapidly spinning blob filaments with a zonal flow layer. It is shown analytically that the flow layer can act like a transport barrier for these structures. Finally parallel electron kinetic effects are incorporated into the theory. Various asymptotic parameter regimes are discussed and asymptotic expressions for the radial and poloidal motion of the blob-filaments are obtained.Comment: 31 pages, 2 figures, accepted in the journal Physics of Plasmas 30, 072302 (2023

Hypoxic pre-conditioning increases the infiltration of endothelial cells into scaffolds for dermal regeneration pre-seeded with mesenchymal stem cells.

Many therapies using mesenchymal stem cells (MSC) rely on their ability to produce and release paracrine signals with chemotactic and pro-angiogenic activity. These characteristics, however, are mostly studied under standard in vitro culture conditions. In contrast, various novel cell-based therapies imply pre-seeding MSC into bio-artificial scaffolds. Here we describe human bone marrow-derived MSC seeded in Integra matrices, a common type of scaffold for dermal regeneration (SDR). We show and measured the distribution of MSC within the SDR, where cells clearly establish physical interactions with the scaffold, exhibiting constant metabolic activity for at least 15 days. In the SDR, MSC secrete VEGF and SDF-1α and induce transwell migration of CD34(+) hematopoietic/endothelial progenitor cells, which is inhibited in the presence of a CXCR4/SDF-1α antagonist. MSC in SDR respond to hypoxia by altering levels of angiogenic signals such as Angiogenin, Serpin-1, uPA, and IL-8. Finally, we show that MSC-containing SDR that have been pre-incubated in hypoxia show higher infiltration of endothelial cells after implantation into immune deficient mice. Our data show that MSC are fully functional ex vivo when implanted into SDR. In addition, our results strongly support the notion of hypoxic pre-conditioning MSC-containing SDR, in order to promote angiogenesis in the wounds

Radio-frequency wave interactions with a plasma sheath in oblique-angle magnetic fields using a sheath impedance model

The physics of interactions between waves in plasmas and sheaths for background magnetic fields which make oblique angles with sheath surfaces is studied with the use of the self-consistent finite element code rfSOL incorporating the recently developed sheath impedance model [J. R. Myra, Phys. Plasmas 24, 072507 (2017)]. The calculation based on this model employs the generalized sheath boundary condition (sheath BC), which surpasses the previously used capacitive sheath BC in reliability by taking into account the contributions of the ion and electron currents in the sheath and the displacement current. A series of numerical simulations is carried out in two-dimensional slab geometry with a flat or curved sheath surface as part of the boundary. It is shown that the sheath–plasma wave appears when the equilibrium magnetic field line angle with respect to the sheath surface is small, the absolute value of the radio-frequency (RF) sheath voltage is large, and the plasma density is slightly higher than the lower hybrid resonance density (LHR density), all of which bring the sheath property closer to being capacitive. It is also shown that the sharp variation of the magnetic field line angle along the sheath surface can sensitively affect the maximum absolute value of the RF sheath voltage at a plasma density slightly lower than the LHR density

Calculation of RF sheath properties from surface wave-fields: a post-processing method

In ion cyclotron range of frequency (ICRF) experiments in fusion research devices, radio frequency (RF) sheaths form where plasma, strong RF wave fields and material surfaces coexist. These RF sheaths affect plasma material interactions such as sputtering and localized power deposition, as well as the global RF wave fields themselves. RF sheaths may be modeled by employing a sheath boundary condition (BC) in place of the more customary conducting wall (CW) BC; however, there are still many ICRF computer codes that do not implement the sheath BC. In this paper we present a method for post-processing results obtained with the CW-BC. The post-processing method produces results that are equivalent to those that would have been obtained with the RF sheath BC, under certain assumptions. The post-processing method is also useful for verification of sheath BC implementations and as a guide to interpretation and understanding of the role of RF sheaths and their interactions with the waves that drive them

Linear Analysis and Verification Suite for Edge Turbulence

The edge and scrape-off-layer region of a tokamak plasma is subject to well known resistive and ideal instabilities that are driven by various curvature- and sheath-related mechanisms. While the boundary plasma is typically strongly turbulent in experiments, it is useful to have computational tools that can analyze the linear eigenmode structure, predict quantitative trends in growth rates and elucidate and the underlying drive mechanisms. Furthermore, measurement of the linear growth rate of unstable modes emerging from a known, established equilibrium configuration provides one of the few quantitative ways of rigorously benchmarking large-scale plasma turbulence codes with each other and with a universal standard. In this report, a suite of codes that can describe linearized, nonlocal (e.g. separatrix-spanning) modes in axisymmetric (realistic divertor), toroidal geometry is discussed. Examples of several benchmark comparisons are given, and future development plans for a new eigenvalue edge code are presented

Radio frequency wave interactions with a plasma sheath: The role of wave and plasma sheath impedances

Radio frequency (RF) sheaths form near surfaces where plasma and strong RF fields coexist. The effect of these RF sheaths on wave propagation near the boundary can be characterized by an effective sheath impedance that includes both resistive and capacitive contributions describing RF sheath rectification and RF power absorption in the sheath [J. R. Myra and D. A. D\u27Ippolito, Phys. Plasmas 22, 062507 (2015)]. Here, we define a dimensionless parameter, the ratio of incoming wave impedance to the sheath impedance, which determines the characteristics of the interaction, ranging from quasi-conducting to quasi-insulating, or in the case of matched impedances, to either perfect absorption or a sheath-plasma resonance. A semi-analytical analysis is carried out for electrostatic slow waves in the ion cyclotron range of frequencies. For the propagating slow wave case, where the incident wave is partially reflected, the fraction of power dissipated in the sheath is calculated. For the evanescent slow wave case, which admits a sheath-plasma resonance, an amplification factor is calculated. Using the impedance ratio approach, RF sheath interactions are characterized for a range of RF wave and plasma parameters including plasma density, magnetic field angle with respect to the surface, wave frequency, and wave-vector components tangent to the surface. For a particularly interesting example case, results are compared with the rfSOL code [H. Kohno and J. R. Myra, Comput. Phys. Commun. 220, 129 (2017)]. Finally, electromagnetic effects, absent from the semi-analytical analysis, are assessed