Extension and Partitioning in an Oblique Subduction Zone, New Zealand: Constraints from Three-Dimensional Numerical Modeling

Contraction, strike slip, and extension displacements along the Hikurangi margin northeast of the North Island of New Zealand coincide with large lateral gradients in material properties. We use a finite- difference code utilizing elastic and elastic-plastic rheologies to build large- scale, three-dimensional numerical models which investigate the influence of material properties on velocity partitioning within oblique subduction zones. Rheological variation in the oblique models is constrained by seismic velocity and attenuation information available for the Hikurangi margin. We compare the effect of weakly versus strongly coupled subduction interfaces on the development of extension and the partitioning of velocity components for orthogonal and oblique convergence and include the effect of ponded sediments beneath the Raukumara Peninsula. Extension and velocity partitioning occur if the subduction interface is weak, but neither develops if the subduction interface is strong. The simple mechanical model incorporating rheological variation based on seismic observations produces kinematics that closely match those published from the Hikurangi margin. These include extension within the Taupo Volcanic Zone, uplift over ponded sediments, and dextral contraction to the south.

Three-Dimensional Mechanics of Yakutat Convergence in the Southern Alaskan Plate Corner

Three-dimensional numerical models are used to investigate the mechanical evolution of the southern Alaskan plate corner where the Yakutat and the Pacific plates converge on the North American plate. The evolving model plate boundary consists of Convergent, Lateral, and Subduction subboundaries with flow separation of incoming material into upward or downward trajectories forming dual, nonlinear advective thermal/mechanical anomalies that fix the position of major subaerial mountain belts. The model convergent subboundary evolves into two teleconnected orogens: Inlet and Outlet orogens form at locations that correspond with the St. Elias and the Central Alaska Range, respectively, linked to the East by the Lateral boundary. Basins form parallel to the orogens in response to the downward component of velocity associated with subduction. Strain along the Lateral subboundary varies as a function of orogen rheology and magnitude and distribution of erosion. Strain-dependent shear resistance of the plate boundary associated with the shallow subduction zone controls the position of the Inlet orogen. The linkages among these plate boundaries display maximum shear strain rates in the horizontal and vertical planes where the Lateral subboundary joins the Inlet and Outlet orogens. The location of the strain maxima shifts with time as the separation of the Inlet and Outlet orogens increases. The spatiotemporal predictions of the model are consistent with observed exhumation histories deduced from thermochronology, as well as stratigraphic studies of synorogenic deposits. In addition, the complex structural evolution of the St Elias region is broadly consistent with the predicted strain field evolution. Citation: Koons, P. O., B. P. Hooks, T. Pavlis, P. Upton, and A. D. Barker (2010), Three-dimensional mechanics of Yakutat convergence in the southern Alaskan plate corner, Tectonics, 29, TC4008, doi: 10.1029/2009TC002463.

Model-Based Estimates of Calanus Finmarchicus Abundance in the Gulf of Maine

Ocean observing systems and satellites routinely collect a wealth of information on physical conditions in the ocean. With few exceptions, such as chlorophyll concentrations, information on biological properties is harder to measure autonomously. Here, we present a system to produce estimates of the distribution and abundance of the copepod Calanus finmarchicus in the Gulf of Maine. Our system uses satellite-based measurements of sea surface temperature and chlorophyll concentration to determine the developmental and reproductive rates of C. finmarchicus. The rate information then drives a population dynamics model of C. finmarchicus that is embedded in a 2-dimensional circulation field. The first generation of this system produces realistic information on interannual variability in C. finmarchicus distribution and abundance during the winter and spring. The model can also be used to identify key drivers of interannual variability in C. finmarchicus. Experiments with the model suggest that changes in initial conditions are overwhelmed by variability in growth rates after approximately 50 d. Temperature has the largest effect on growth rate. Elevated chlorophyll during the late winter can lead to increased C. finmarchicus abundance during the spring, but the effect of variations in chlorophyll concentrations is secondary to the other inputs. Our system could be used to provide real-time estimates or even forecasts of C. finmarchicus distribution. These estimates could then be used to support management of copepod predators such as herring and right whales.

Influence of microphysical schemes on atmospheric water in the Weather Research and Forecasting model

This study examines how different microphysical parameterization schemes influence orographically induced precipitation and the distributions of hydrometeors and water vapour for midlatitude summer conditions in the Weather Research and Forecasting (WRF) model. A high-resolution two-dimensional idealized simulation is used to assess the differences between the schemes in which a moist air flow is interacting with a bell-shaped 2 km high mountain. Periodic lateral boundary conditions are chosen to recirculate atmospheric water in the domain. It is found that the 13 selected microphysical schemes conserve the water in the model domain. The gain or loss of water is less than 0.81% over a simulation time interval of 61 days. The differences of the microphysical schemes in terms of the distributions of water vapour, hydrometeors and accumulated precipitation are presented and discussed. The Kessler scheme, the only scheme without ice-phase processes, shows final values of cloud liquid water 14 times greater than the other schemes. The differences among the other schemes are not as extreme, but still they differ up to 79% in water vapour, up to 10 times in hydrometeors and up to 64% in accumulated precipitation at the end of the simulation. The microphysical schemes also differ in the surface evaporation rate. The WRF single-moment 3-class scheme has the highest surface evaporation rate compensated by the highest precipitation rate. The different distributions of hydrometeors and water vapour of the microphysical schemes induce differences up to 49 W m−2 in the downwelling shortwave radiation and up to 33 W m−2 in the downwelling longwave radiation.

Monte Carlo based beam model using a photon MLC for modulated electron radiotherapy

PURPOSE Modulated electron radiotherapy (MERT) promises sparing of organs at risk for certain tumor sites. Any implementation of MERT treatment planning requires an accurate beam model. The aim of this work is the development of a beam model which reconstructs electron fields shaped using the Millennium photon multileaf collimator (MLC) (Varian Medical Systems, Inc., Palo Alto, CA) for a Varian linear accelerator (linac). METHODS This beam model is divided into an analytical part (two photon and two electron sources) and a Monte Carlo (MC) transport through the MLC. For dose calculation purposes the beam model has been coupled with a macro MC dose calculation algorithm. The commissioning process requires a set of measurements and precalculated MC input. The beam model has been commissioned at a source to surface distance of 70 cm for a Clinac 23EX (Varian Medical Systems, Inc., Palo Alto, CA) and a TrueBeam linac (Varian Medical Systems, Inc., Palo Alto, CA). For validation purposes, measured and calculated depth dose curves and dose profiles are compared for four different MLC shaped electron fields and all available energies. Furthermore, a measured two-dimensional dose distribution for patched segments consisting of three 18 MeV segments, three 12 MeV segments, and a 9 MeV segment is compared with corresponding dose calculations. Finally, measured and calculated two-dimensional dose distributions are compared for a circular segment encompassed with a C-shaped segment. RESULTS For 15 × 34, 5 × 5, and 2 × 2 cm(2) fields differences between water phantom measurements and calculations using the beam model coupled with the macro MC dose calculation algorithm are generally within 2% of the maximal dose value or 2 mm distance to agreement (DTA) for all electron beam energies. For a more complex MLC pattern, differences between measurements and calculations are generally within 3% of the maximal dose value or 3 mm DTA for all electron beam energies. For the two-dimensional dose comparisons, the differences between calculations and measurements are generally within 2% of the maximal dose value or 2 mm DTA. CONCLUSIONS The results of the dose comparisons suggest that the developed beam model is suitable to accurately reconstruct photon MLC shaped electron beams for a Clinac 23EX and a TrueBeam linac. Hence, in future work the beam model will be utilized to investigate the possibilities of MERT using the photon MLC to shape electron beams.

Spontaneous supersymmetry breaking in the 2D N=1 Wess-Zumino Model

We study the phase diagram of the two-dimensional N=1 Wess-Zumino model on the lattice using Wilson fermions and the fermion loop formulation. We give a complete nonperturbative determination of the ground state structure in the continuum and infinite volume limit. We also present a determination of the particle spectrum in the supersymmetric phase, in the supersymmetry broken phase and across the supersymmetry breaking phase transition. In the supersymmetry broken phase, we observe the emergence of the Goldstino particle.

Three-dimensional printing for perioperative planning of complex aortic arch surgery

Purpose In this study, we show the use of three-dimensional printing models for preoperative planning of surgery for patients with complex aortic arch anomalies. Description A 70-year-old man with an extensively arteriosclerotic aneurysm reaching from the ascending aorta to the descending aorta was referred to our center for complete aortic arch replacement. We visualized and reconstructed computed tomography data of the patient and fabricated a flexible three-dimensional model of the aortic arch including the aneurysm. Evaluation This model was very helpful for the preoperative decision making and planning of the frozen elephant trunk procedure owing to the exact and lifelike illustration of the native aortic arch. Conclusions Three-dimensional models are helpful in preoperative planning and postoperative evaluation of frozen elephant trunk procedures in patients with complex aortic anatomy.

Articulated Statistical Shape Model-Based 2D-3D Reconstruction of a Hip Joint

In this paper, reconstruction of three-dimensional (3D) patient-specific models of a hip joint from two-dimensional (2D) calibrated X-ray images is addressed. Existing 2D-3D reconstruction techniques usually reconstruct a patient-specific model of a single anatomical structure without considering the relationship to its neighboring structures. Thus, when those techniques would be applied to reconstruction of patient-specific models of a hip joint, the reconstructed models may penetrate each other due to narrowness of the hip joint space and hence do not represent a true hip joint of the patient. To address this problem we propose a novel 2D-3D reconstruction framework using an articulated statistical shape model (aSSM). Different from previous work on constructing an aSSM, where the joint posture is modeled as articulation in a training set via statistical analysis, here it is modeled as a parametrized rotation of the femur around the joint center. The exact rotation of the hip joint as well as the patient-specific models of the joint structures, i.e., the proximal femur and the pelvis, are then estimated by optimally fitting the aSSM to a limited number of calibrated X-ray images. Taking models segmented from CT data as the ground truth, we conducted validation experiments on both plastic and cadaveric bones. Qualitatively, the experimental results demonstrated that the proposed 2D-3D reconstruction framework preserved the hip joint structure and no model penetration was found. Quantitatively, average reconstruction errors of 1.9 mm and 1.1 mm were found for the pelvis and the proximal femur, respectively.

The Utilization of Classical Spin Monte Carlo Methods to Simulate the Magnetic Behavior of Extended Three-Dimensional Cubic Networks Incorporating M(II) Ions with an S = 5/2 Ground State Spin

The numerical simulations of the magnetic properties of extended three-dimensional networks containing M(II) ions with an S = 5/2 ground-state spin have been carried out within the framework of the isotropic Heisenberg model. Analytical expressions fitting the numerical simulations for the primitive cubic, diamond, together with (10−3) cubic networks have all been derived. With these empirical formulas in hands, we can now extract the interaction between the magnetic ions from the experimental data for these networks. In the case of the primitive cubic network, these expressions are directly compared with those from the high-temperature expansions of the partition function. A fit of the experimental data for three complexes, namely [(N(CH3)4][Mn(N3)] 1, [Mn(CN4)]n 2, and [FeII(bipy)3][MnII2(ox)3] 3, has been carried out. The best fits were those obtained using the following parameters, J = −3.5 cm-1, g = 2.01 (1); J = −8.3 cm-1, g = 1.95 (2); and J = −2.0 cm-1, g = 1.95 (3).

Luminescence and Energy Transfer of [Ru(bpy) 3 ] 2+ , [Cr(ox) 3 ] 3- , and [Os(bpy) 3 ] 2+ in Three-Dimensional Oxalato-Networks

Luminescence and energy transfer in [Zn1-xRux(bpy)3][NaAl1-yCry(ox)3] (x ≈ 0.01, y = 0.006 − 0.22; bpy = 2,2‘-bipyridine, ox = C2O42-) and [Zn1-x-yRuxOsy(bpy)3][NaAl(ox)3] (x ≈ 0.01, y = 0.012) are presented and discussed. Surprisingly, the luminescence of the isolated luminophores [Ru(bpy)3]2+ and [Os(bpy)3]2+ in [Zn(bpy)3][NaAl(ox)3] is hardly quenched at room temperature. Steady-state luminescence spectra and decay curves show that energy transfer occurs between [Ru(bpy)3]2+ and [Cr(ox)3]3- and between [Ru(bpy)3]2+ and [Os(bpy)3]2+ in [Zn1-xRux(bpy)3][NaAl1-yCry(ox)3] and [Zn1-x-yRuxOsy(bpy)3] [NaAl(ox)3], respectively. For a quantitative investigation of the energy transfer, a shell type model is developed, using a Monte Carlo procedure and the structural parameters of the systems. A good description of the experimental data is obtained assuming electric dipole−electric dipole interaction between donors and acceptors, with a critical distance Rc for [Ru(bpy)3]2+ to [Cr(ox)3]3- energy transfer of 15 Å and for [Ru(bpy)3]2+ to [Os(bpy)3]2+ energy transfer of 33 Å. These values are in good agreement with those derived using the Förster−Dexter theory.

Fluid-structure interaction of an elastic pipe immersed inside a coaxial rigid pipe: a model for the Tullio phenomenon

An axisymmetric, elastic pipe is filled with an incompressible fluid and is immersed in a second, coaxial rigid pipe which contains the same fluid. A pressure pulse in the outer fluid annulus deforms the elastic pipe which invokes a fluid motion in the fluid core. It is the aim of this study to investigate streaming phenomena in the core which may originate from such a fluid-structure interaction. This work presents a numerical solver for such a configuration. It was developed in the OpenFOAM software environment and is based on the Arbitrary Lagrangian Eulerian (ALE) approach for moving meshes. The solver features a monolithic integration of the one-dimensional, coupled system between the elastic structure and the outer fluid annulus into a dynamic boundary condition for the moving surface of the fluid core. Results indicate that our configuration may serve as a mechanical model of the Tullio Phenomenon (sound-induced vertigo).

Achieving a more realistic assessment of rockfall hazards by coupling three-dimensional process models and field-based tree-ring data

Sound knowledge of the spatial and temporal patterns of rockfalls is fundamental for the management of this very common hazard in mountain environments. Process-based, three-dimensional simulation models are nowadays capable of reproducing the spatial distribution of rockfall occurrences with reasonable accuracy through the simulation of numerous individual trajectories on highly-resolved digital terrain models. At the same time, however, simulation models typically fail to quantify the ‘real’ frequency of rockfalls (in terms of return intervals). The analysis of impact scars on trees, in contrast, yields real rockfall frequencies, but trees may not be present at the location of interest and rare trajectories may not necessarily be captured due to the limited age of forest stands. In this article, we demonstrate that the coupling of modeling with tree-ring techniques may overcome the limitations inherent to both approaches. Based on the analysis of 64 cells (40 m × 40 m) of a rockfall slope located above a 1631-m long road section in the Swiss Alps, we illustrate results from 488 rockfalls detected in 1260 trees. We illustrate that tree impact data cannot only be used (i) to reconstruct the real frequency of rockfalls for individual cells, but that they also serve (ii) the calibration of the rockfall model Rockyfor3D, as well as (iii) the transformation of simulated trajectories into real frequencies. Calibrated simulation results are in good agreement with real rockfall frequencies and exhibit significant differences in rockfall activity between the cells (zones) along the road section. Real frequencies, expressed as rock passages per meter road section, also enable quantification and direct comparison of the hazard potential between the zones. The contribution provides an approach for hazard zoning procedures that complements traditional methods with a quantification of rockfall frequencies in terms of return intervals through a systematic inclusion of impact records in trees.

Soft Tissue Alterations in Esthetic Postextraction Sites: A 3-Dimensional Analysis.

Dimensional alterations of the facial soft and bone tissues following tooth extraction in the esthetic zone play an essential role to achieve successful outcomes in implant therapy. This prospective study is the first to investigate the interplay between the soft tissue dimensions and the underlying bone anatomy during an 8-wk healing period. The analysis is based on sequential 3-dimensional digital surface model superimpositions of the soft and bone tissues using digital impressions and cone beam computed tomography during an 8-wk healing period. Soft tissue thickness in thin and thick bone phenotypes at extraction was similar, averaging 0.7 mm and 0.8 mm, respectively. Interestingly, thin bone phenotypes revealed a 7-fold increase in soft tissue thickness after an 8-wk healing period, whereas in thick bone phenotypes, the soft tissue dimensions remained unchanged. The observed spontaneous soft tissue thickening in thin bone phenotypes resulted in a vertical soft tissue loss of only 1.6 mm, which concealed the underlying vertical bone resorption of 7.5 mm. Because of spontaneous soft tissue thickening, no significant differences were detected in the total tissue loss between thin and thick bone phenotypes at 2, 4, 6, and 8 wk. More than 51% of these dimensional alterations occurred within 2 wk of healing. Even though the observed spontaneous soft tissue thickening in thin bone phenotypes following tooth extraction conceals the pronounced underlying bone resorption pattern by masking the true bone deficiency, spontaneous soft tissue thickening offers advantages for subsequent bone regeneration and implant therapies in sites with high esthetic demand (Clinicaltrials.gov NCT02403700).

Investigation of hemodynamics in an in vitro system simulating left ventricular support through the right subclavian artery using 4-dimensional flow magnetic resonance imaging.

OBJECTIVES Left ventricular assist devices are an important treatment option for patients with heart failure alter the hemodynamics in the heart and great vessels. Because in vivo magnetic resonance studies of patients with ventricular assist devices are not possible, in vitro models represent an important tool to investigate flow alterations caused by these systems. By using an in vitro magnetic resonance-compatible model that mimics physiologic conditions as close as possible, this work investigated the flow characteristics using 4-dimensional flow-sensitive magnetic resonance imaging of a left ventricular assist device with outflow via the right subclavian artery as commonly used in cardiothoracic surgery in the recent past. METHODS An in vitro model was developed consisting of an aorta with its supra-aortic branches connected to a left ventricular assist device simulating the pulsatile flow of the native failing heart. A second left ventricular assist device supplied the aorta with continuous flow via the right subclavian artery. Four-dimensional flow-sensitive magnetic resonance imaging was performed for different flow rates of the left ventricular assist device simulating the native heart and the left ventricular assist device providing the continuous flow. Flow characteristics were qualitatively and quantitatively evaluated in the entire vessel system. RESULTS Flow characteristics inside the aorta and its upper branching vessels revealed that the right subclavian artery and the right carotid artery were solely supported by the continuous-flow left ventricular assist device for all flow rates. The flow rates in the brain-supplying arteries are only marginally affected by different operating conditions. The qualitative analysis revealed only minor effects on the flow characteristics, such as weakly pronounced vortex flow caused by the retrograde flow via the brachiocephalic artery. CONCLUSIONS The results indicate that, despite the massive alterations in natural hemodynamics due to the retrograde flow via the right subclavian and brachiocephalic arteries, there are no drastic consequences on the flow in the brain-feeding arteries and the flow characteristics in the ascending and descending aortas. It may be beneficial to adjust the operating condition of the left ventricular assist device to the residual function of the failing heart.

Effects of Surface Forcing on Interannual Variability of the Fall Phytoplankton Bloom in the Gulf of Maine Revealed Using a Process-Oriented Model

SeaWiFS (Sea-viewing Wide Field-of-view Sensor) chlorophyll data revealed strong interannual variability in fall phytoplankton dynamics in the Gulf of Maine, with 3 general features in any one year: (1) rapid chlorophyll increases in response to storm events in fall; (2) gradual chlorophyll increases in response to seasonal wind-and cooling-induced mixing that gradually deepens the mixed layer; and (3) the absence of any observable fall bloom. We applied a mixed-layer box model and a 1-dimensional physical-biological numerical model to examine the influence of physical forcing (surface wind, heat flux, and freshening) on the mixed-layer dynamics and its impact on the entrainment of deep-water nutrients and thus on the appearance of fall bloom. The model results suggest that during early fall, the surface mixed-layer depth is controlled by both wind-and cooling-induced mixing. Strong interannual variability in mixed-layer depth has a direct impact on short-and long-term vertical nutrient fluxes and thus the fall bloom. Phytoplankton concentrations over time are sensitive to initial pre-bloom profiles of nutrients. The strength of the initial stratification can affect the modeled phytoplankton concentration, while the timing of intermittent freshening events is related to the significant interannual variability of fall blooms.