From transient to steady state deformation and grain size: A thermodynamic approach using elasto-visco-plastic numerical modeling

Numerical simulation experiments give insight into the evolving energy partitioning during high-strain torsion experiments of calcite. Our numerical experiments are designed to derive a generic macroscopic grain size sensitive flow law capable of describing the full evolution from the transient regime to steady state. The transient regime is crucial for understanding the importance of micro structural processes that may lead to strain localization phenomena in deforming materials. This is particularly important in geological and geodynamic applications where the phenomenon of strain localization happens outside the time frame that can be observed under controlled laboratory conditions. Ourmethod is based on an extension of the paleowattmeter approach to the transient regime. We add an empirical hardening law using the Ramberg-Osgood approximation and assess the experiments by an evolution test function of stored over dissipated energy (lambda factor). Parameter studies of, strain hardening, dislocation creep parameter, strain rates, temperature, and lambda factor as well asmesh sensitivity are presented to explore the sensitivity of the newly derived transient/steady state flow law. Our analysis can be seen as one of the first steps in a hybrid computational-laboratory-field modeling workflow. The analysis could be improved through independent verifications by thermographic analysis in physical laboratory experiments to independently assess lambda factor evolution under laboratory conditions.

Combined effect of heterogeneity, anisotropy and saturation on steady state flow and transport: Structure recognition and numerical simulation

1 Natural soil profiles may be interpreted as an arrangement of parts which are characterized by properties like hydraulic conductivity and water retention function. These parts form a complicated structure. Characterizing the soil structure is fundamental in subsurface hydrology because it has a crucial influence on flow and transport and defines the patterns of many ecological processes. We applied an image analysis method for recognition and classification of visual soil attributes in order to model flow and transport through a man-made soil profile. Modeled and measured saturation-dependent effective parameters were compared. We found that characterizing and describing conductivity patterns in soils with sharp conductivity contrasts is feasible. Differently, solving flow and transport on the basis of these conductivity maps is difficult and, in general, requires special care for representation of small-scale processes.

Numerical solution of optimal control problems with constant control delays

We investigate a class of optimal control problems that exhibit constant exogenously given delays in the control in the equation of motion of the differential states. Therefore, we formulate an exemplary optimal control problem with one stock and one control variable and review some analytic properties of an optimal solution. However, analytical considerations are quite limited in case of delayed optimal control problems. In order to overcome these limits, we reformulate the problem and apply direct numerical methods to calculate approximate solutions that give a better understanding of this class of optimization problems. In particular, we present two possibilities to reformulate the delayed optimal control problem into an instantaneous optimal control problem and show how these can be solved numerically with a stateof- the-art direct method by applying Bock’s direct multiple shooting algorithm. We further demonstrate the strength of our approach by two economic examples.

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.

Glacier Calving: A Numerical Model of Forces in the Calving-Speed/Water-Depth Relation

Empirical data suggest that the race of calving of grounded glaciers terminating in water is directly proportional to the water depth. Important controls on calving may be the extent to which a calving face tends to become oversteepened by differential flow within the ice and the extent to which bending moments promote extrusion and bottom crevassing at the base of a calving face. Numerical modelling suggests that the tendency to become oversteepened increases roughly linearly with water depth. In addition, extending longitudinal deviatoric stresses at the base of a calving face increase with water depth. These processes provide a possible physical explanation for the observed calving-rate/water-depth relation.

A Numerical Investigation of the Gulf Stream and Its Meanders in Response to Cold Air Outbreaks

The three-dimensional Princeton Ocean Model is used to examine the modification of the Gulf Stream and its meanders by cold air outbreaks. Two types of Gulf Stream meanders are found in the model. Meanders on the shoreward side of the Gulf Stream are baroclinically unstable. They are affected little by the atmospheric forcing because their energy source is stored at the permanent thermocline, well below the influence of the surface forcing. Meanders on the seaward side of the stream are both barotropically and baroclinically unstable. The energy feeding these meanders is stored at the surface front separating the Gulf Stream and the Sargasso Seal which is greatly reduced in case of cold air outbreaks. Thus, meanders there reduce strength and also seem to slow their downstream propagation due to the southward Ekman flow. Heat budget calculations suggest two almost separable processes. The oceanic heal released to the atmosphere during these severe cooling episodes comes almost exclusively from the upper water column. Transport of heat by meanders from the Gulf Stream to the shelf, though it is large, does not disrupt the principal balance. It is balanced nicely with the net heat transport in the downstream direction.