Vertical heat transfer based on direct microstructure measurements in the ice-free Pacific-side Arctic Ocean: the role and impact of the Pacific water intrusion

This study quantifies diapycnal mixing and vertical heat transfer in the Pacific side of the Arctic Ocean, where sea-ice cover has disappeared between July and September in the last few decades. We conducted microstructure measurements in the open water region around the Canada Basin from late summer to fall in 2009 and 2010 using R/V Mirai. In the study domain, the dissipation rate of turbulent kinetic energy, ε, is typically as low level as O(10^?10) W kg^?1, resulting in vertical heat diffusivity of O(10^?7) m² s^?1, which is close to the molecular diffusivity of heat, suggesting comparatively little predominance of mechanical turbulent mixing. An exception is the case at the Barrow Canyon, where the strong baroclinic throughflow generates substantial vertical mixing, producing ε > O(10^?7) W kg^?1, because of the shear flow instability. Meanwhile, in the confluence region, where the warm/salty Pacific water and the cold/fresh Arctic basin water encounter, the micro-temperature profiles revealed a localized enhancement in vertical diffusivity of heat, reaching O(10^?5) m² s^?1 or greater. In this region, an intrusion of warm Pacific water creates a horizontally interleaved structure, where the double-diffusive mixing facilitates vertical heat transfer between the intruding Pacific water and the surrounding basin waters.     

FDM Simulation of Seismic Waves, Ocean Acoustic Waves, and Tsunamis Based on Tsunami-Coupled Equations of Motion

We have developed a new, unified modeling technique for the total simulation of seismic waves, ocean acoustic waves, and tsunamis resulting from earthquakes, based on a finite difference method simulation of the 3D equations of motion. Using the equilibrium between the pressure gradient and gravity in these equations, tsunami propagation is naturally incorporated in the simulation based on the equations of motion. The performance of the parallel computation for the newly developed tsunami-coupled equations using a domain partitioning procedure shows a high efficiency coefficient with a large number of CPU cores. The simulation results show how the near-field term associated with seismic waves produced by shallow earthquakes leads to a permanent coseismic deformation of the ground surface, which gives rise to the initial tsunami on the sea surface. Propagation of the tsunami along the sea surface as a gravity wave, and ocean acoustic waves in seawater with high-frequency multiple P-wave reflections between the free surface and sea bottom, are also clearly demonstrated by the present simulations. We find a good agreement in the tsunami waveform between our results and those obtained by other simulations based on an analytical model and the Navier–Stokes equations, demonstrating the effectiveness of the tsunami-coupling simulation model. Based on this simulation, we show that the ratio of the amplitude of ocean acoustic waves to the height of the tsunami, both of which are produced by the earthquake, strongly depends on the rise time of the earthquake rupture. This ratio can be used to obtain a more detailed understanding of the source rupture processes of subduction zone earthquakes, and for implementing an improved tsunami alert system for slow tsunami earthquakes.     

Nonlinear sloshing in a floating‐roofed oil storage tank under long‐period seismic ground motion

A hybrid analytical and FEM is proposed to investigate the nonlinear sloshing in a floating‐roofed oil storage tank under long‐period seismic ground motion. The tank is composed of a rigid cylindrical wall and a flat bottom, whereas the floating roof is treated as an elastic plate undergoing large deflection. The contained liquid is assumed to be inviscid and incompressible, and the flow is assumed to be irrotational. The method of analysis is based on representation of the liquid motion by superposing the analytical modes that satisfy the Laplace equation and the rigid wall and bottom boundary conditions. The FEM is then applied to solve the remaining kinematic and dynamic boundary conditions at the moving liquid surface coupled with the nonlinear equation of motion of the floating roof. This requires only the discretization of the liquid surface and the floating roof into finite elements, thus leading to a computationally efficient and accurate method compared with full numerical analysis. As numerical examples to illustrate the applicability of the proposed method, two oil storage tanks with single‐deck type floating roofs damaged during the 2003 Tokachioki earthquake are studied. It is shown that the nonlinear oscillation modes with the circumferential wave numbers 0, 2 and 3 caused by the finite liquid surface elevation as well as the membrane action due to large deflection of the deck produce excessively large stresses in the pontoon, which may cause the catastrophic failure of pontoon followed by the submergence of the roof. Copyright © 2012 John Wiley & Sons, Ltd.     

The distribution of total and methylmercury concentrations in soils near the Idrija mercury mine, Slovenia, and the dependence of the mercury concentrations on the chemical composition and organic carbon levels of the soil

Although the mining activity of the Idrija mine in Slovenia ceased in 1995, a large amount of mining dregs containing high concentrations of mercury remains in the area. The mining dregs were transported with river flow and deposition along the Idrija River. To estimate the dispersion and change in the chemical form of mercury, a total of 28 soil core samples were taken around the river. The individual core samples were separated into layers for the analysis of their chemical composition, carbon contents, total mercury (T-Hg) and methylmercury (MeHg) concentrations. The chemical composition measured by X-ray fluorescence spectrometry was useful to estimate the dispersion of tailings: the fluvial terrace soil had a chemical composition similar to that of the tailings and could be distinguished clearly from the forest soil. The highest T-Hg concentration, 1,100 mg kg⁻¹, was observed in the fluvial terrace soil near the mine. Although the concentration decreased gradually along with distance from the mine, concentrations higher than 200 mg kg⁻¹ of T-Hg were still observed in the fluvial terrace soil approximately 20 km downstream from the mine. In the vertical distribution of T-Hg in the hillslope soil, a higher value was observed in the upper layers, which suggests the recent atmospheric deposition of mercury. The concentration of MeHg was the lowest at the riverside and higher in the hillslope soil, which was the opposite of the T-Hg distribution. The total organic carbon content tracked similarly with the distribution of MeHg and a linear relation with a significantly high correlation coefficient was obtained. The distinction may be related to the different dispersion process of mercury, and the organic carbon contents may be an important factor for MeHg formation.     

Spinel inclusions in olivine and plagioclase crystals in a layered gabbro: a marker and a tracer for primary phenocrysts in a differentiating magma reservoir

The problem of whether cumulate rocks were formed by crystal settling or by in situ crystallization after magma emplacement is an important issue concerning the mechanisms of magmatic differentiation. However, it is hard to distinguish these two processes for plutonic rocks because the primary texture and chemical composition have generally been modified by postcumulus processes. To contribute this problem, we studied the distribution and compositions of Cr-spinel inclusions hosted in olivine and plagioclase in the Murotomisaki Gabbroic Intrusion (MGI), SW Japan. It is shown that the olivine-hosted inclusions are restricted to specific horizons where accumulation of olivine phenocrysts is thought to have occurred and that the compositional variations of the Cr-spinel are explained by a secondary compositional modification that probably took place after the magma emplacement. It is also shown that the Cr-spinel inclusions in a chilled margin have suffered the least compositional modification and nearly retains the primary composition. Those in the interior of the intrusion, on the contrary, have been significantly modified by re-equilibration with residual melt driven by cation diffusions through the host phases. Those in plagioclase have been less modified. It is shown that all the spinel inclusions had primarily the same and common composition at the time of magma emplacement. This implies that all the inclusion-bearing crystals, olivine and plagioclase, represent primary phenocrysts that had already existed in the emplaced magma. In this way, spinel inclusion in the MGI may be regarded to be a useful petrographic “marker” for identifying intratelluric phenocrysts and also as a “tracer” to trace the motion of the primary phenocrysts after the magma emplacement.     

Evidence for ocean water invasion into the Chicxulub crater at the Cretaceous/Tertiary boundary

Abstract— The possibility of ocean water invasion into the Chicxulub crater following the impact at the Cretaceous/Tertiary boundary was investigated based on examination of an impactite between approximately 794.63 and 894.94 m in the Yaxcopoil‐1 (Yax‐1) core. The presence of cross lamination in the uppermost part of the impactite suggests the influence of an ocean current at least during the sedimentation of this interval. Abundant occurrence of nannofossils of late Campanian to early Maastrichtian age in the matrices of samples from the upper part of the impactite suggests that the carbonate sediments deposited on the inner rim margin and outside the crater were eroded and transported into the crater most likely by ocean water that invaded the crater after its formation. The maximum grain size of limestone lithics and vesicular melt fragments, and grain and bulk chemical compositions show a cyclic variation in the upper part of the impactite. The upward fining grain size and the absence of erosional contact at the base of each cycle suggest that the sediments were derived from resuspension of units elsewhere in the crater, most likely by high energy currents association with ocean water invasion.     

Magnetic detection of heated soils at paleolithic sites in Japan

We present a methodological approach to detect heated soil on ancient sites, using magnetic measurements. The method is based on changes in magnetic signals of soil by heating. The following three types of soil were used for testing the method: silty soil (SS), weathered volcanic ash (WVA, = loam) and fairly fresh volcanic ash (VA) called Odori tephra. On heating above 250–600°C, the magnetic susceptibility and remanent magnetization intensity increased for the SS and WVA samples, reflecting chemical alteration of magnetic minerals (from goethites to magnetites through hematites). The VA sample showed no susceptibility change suggesting the absence of goethites within it. On heating below 250°C, only the intensities of all the samples increased. This is possibly due to acquisition of thermal remanent magnetization. The largest change of the magnetic signals was identified for the SS sample and the smallest one was seen for the VA sample. Therefore, the in situ susceptibility measurement, which is the nondestructive and indirect method, seems to be effective to detect heated soil for sites of aqueous deposits as the SS. On the other hand, for sites of aeolian deposits as the WVA (loam) and VA, the intensity measurement of collected soils seems to be the most reliable method to detect evidence of heating. The degree of the magnetic stability (coercivity) against progressive alternating-field demagnetization was also an important parameter, indicating whether the investigated soils were heated or unheated. © 1999 John Wiley & Sons, Inc.