Impact of periodic intermediary flows on submarine melting of a Greenland glacier

Author Posting. © American Geophysical Union, 2014. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research: Oceans 119 (2014): 7078–7098, doi:10.1002/2014JC009953.The submarine melting of a vertical glacier front, induced by an intermediary circulation forced by periodic density variations at the mouth of a fjord, is investigated using a nonhydrostatic ocean general circulation model and idealized laboratory experiments. The idealized configurations broadly match that of Sermilik Fjord, southeast Greenland, a largely two layers system characterized by strong seasonal variability of subglacial discharge. Consistent with observations, the numerical results suggest that the intermediary circulation is an effective mechanism for the advection of shelf anomalies inside the fjord. In the numerical simulations, the advection mechanism is a density intrusion with a velocity which is an order of magnitude larger than the velocities associated with a glacier-driven circulation. In summer, submarine melting is mostly influenced by the discharge of surface runoff at the base of the glacier and the intermediary circulation induces small changes in submarine melting. In winter, on the other hand, submarine melting depends only on the water properties and velocity distribution at the glacier front. Hence, the properties of the waters advected by the intermediary circulation to the glacier front are found to be the primary control of the submarine melting. When the density of the intrusion is intermediate between those found in the fjord's two layers, there is a significant reduction in submarine melting. On the other hand, when the density is close to that of the bottom layer, only a slight reduction in submarine melting is observed. The numerical results compare favorably to idealized laboratory experiments with a similar setup.Support to C. Cenedese and F. Straneo was given by the National Science Foundation project OCE-1130008. C. Cenedese received support also from the WHOI Arctic Research Initiative. R. Sciascia and P. Heimbach are supported in part by NSF project OCE-1129746. Additional funding for P. Heimbach comes through NASA's project NNH11ZDA001N-IDS A.28.2015-04-2

Imaging stray magnetic field of individual ferromagnetic nanotubes

We use a scanning nanometer-scale superconducting quantum interference device to map the stray magnetic field produced by individual ferromagnetic nanotubes (FNTs) as a function of applied magnetic field. The images are taken as each FNT is led through magnetic reversal and are compared with micromagnetic simulations, which correspond to specific magnetization configurations. In magnetic fields applied perpendicular to the FNT long axis, their magnetization appears to reverse through vortex states, i.e.\ configurations with vortex end domains or -- in the case of a sufficiently short FNT -- with a single global vortex. Geometrical imperfections in the samples and the resulting distortion of idealized mangetization configurations influence the measured stray-field patterns.Comment: 14 pages, 4 figure

Magnetization reversal of an individual exchange biased permalloy nanotube

We investigate the magnetization reversal mechanism in an individual permalloy (Py) nanotube (NT) using a hybrid magnetometer consisting of a nanometer-scale SQUID (nanoSQUID) and a cantilever torque sensor. The Py NT is affixed to the tip of a Si cantilever and positioned in order to optimally couple its stray flux into a Nb nanoSQUID. We are thus able to measure both the NT's volume magnetization by dynamic cantilever magnetometry and its stray flux using the nanoSQUID. We observe a training effect and temperature dependence in the magnetic hysteresis, suggesting an exchange bias. We find a low blocking temperature $T_B = 18 \pm 2$ K, indicating the presence of a thin antiferromagnetic native oxide, as confirmed by X-ray absorption spectroscopy on similar samples. Furthermore, we measure changes in the shape of the magnetic hysteresis as a function of temperature and increased training. These observations show that the presence of a thin exchange-coupled native oxide modifies the magnetization reversal process at low temperatures. Complementary information obtained via cantilever and nanoSQUID magnetometry allows us to conclude that, in the absence of exchange coupling, this reversal process is nucleated at the NT's ends and propagates along its length as predicted by theory.Comment: 8 pages, 4 figure

Photo- and Electro-Disintegration of 3He at Threshold and pd Radiative Capture

The present work reports results for: pd radiative capture observables measured at center-of-mass (c.m.) energies in the range 0--100 keV and at 2 MeV by the TUNL and Wisconsin groups, respectively; contributions to the Gerasimov-Drell-Hearn (GDH) integral in 3He from the two- up to the three-body breakup thresholds, compared to experimental determinations by the TUNL group in this threshold region; longitudinal, transverse, and interference response functions measured in inclusive polarized electron scattering off polarized 3He at excitation energies below the threshold for breakup into ppn, compared to unpolarized longitudinal and transverse data from the Saskatoon group. The calculations are based on a realistic Hamiltonian with two- and three-nucleon interactions and a realistic current operator, including one- and two-body components. The theoretical predictions obtained by including only one-body currents are in violent disagreement with data. These differences between theory and experiment are, to a large extent, removed when two-body currents are taken into account, although some rather large discrepancies remain in the c.m. energy range 0--100 keV, particularly for the pd differential cross section and tensor analyzing power at small angles, and contributions to the GDH integral. A rather detailed analysis indicates that these discrepancies have, in large part, a common origin, and can be traced back to an excess strength obtained in the theoretical calculation of the E1 reduced matrix element associated with the pd channel having L,S,J=1,1/2,3/2. It is suggested that this lack of E1 strength observed experimentally might have implications for the nuclear interaction at very low energies. Finally, the validity of the long-wavelength approximation for electric dipole transitions is discussed.Comment: 47 pages RevTex file, 10 PostScript figures, submitted to Phys. Rev.

Model Calculations for the Two-Fragment Electro-Disintegration of 4^4He

Differential cross sections for the electro-disintegration process $e + {^4He} \longrightarrow {^3H}+ p + e'$ are calculated, using a model in which the final state interaction is included by means of a nucleon-nucleus (3+1) potential constructed via Marchenko inversion. The required bound-state wave functions are calculated within the integrodifferential equation approach (IDEA). In our model the important condition that the initial bound state and the final scattering state are orthogonal is fulfilled. The sensitivity of the cross section to the input $p{^3H}$ interaction in certain kinematical regions is investigated. The approach adopted could be useful in reactions involving few cluster systems where effective interactions are not well known and exact methods are presently unavailable. Although, our Plane-Wave Impulse Approximation results exhibit, similarly to other calculations, a dip in the five-fold differential cross-section around a missing momentum of $\sim 450 MeV/c$, it is argued that this is an artifact of the omission of re-scattering four-nucleon processes.Comment: 16 pages, 6 figures, accepted for publication by Phys.Rev.

Automated derivation of the adjoint of high-level transient finite element programs

In this paper we demonstrate a new technique for deriving discrete adjoint and tangent linear models of finite element models. The technique is significantly more efficient and automatic than standard algorithmic differentiation techniques. The approach relies on a high-level symbolic representation of the forward problem. In contrast to developing a model directly in Fortran or C++, high-level systems allow the developer to express the variational problems to be solved in near-mathematical notation. As such, these systems have a key advantage: since the mathematical structure of the problem is preserved, they are more amenable to automated analysis and manipulation. The framework introduced here is implemented in a freely available software package named dolfin-adjoint, based on the FEniCS Project. Our approach to automated adjoint derivation relies on run-time annotation of the temporal structure of the model, and employs the FEniCS finite element form compiler to automatically generate the low-level code for the derived models. The approach requires only trivial changes to a large class of forward models, including complicated time-dependent nonlinear models. The adjoint model automatically employs optimal checkpointing schemes to mitigate storage requirements for nonlinear models, without any user management or intervention. Furthermore, both the tangent linear and adjoint models naturally work in parallel, without any need to differentiate through calls to MPI or to parse OpenMP directives. The generality, applicability and efficiency of the approach are demonstrated with examples from a wide range of scientific applications