Sublimation’s impact on temporal change of albedo dichotomy on Iapetus

Iapetus, one of the saturnian moons, has an extreme albedo contrast between the leading and trailing hemispheres. The origin of this albedo dichotomy has led to several hypotheses, however it remains controversial. To clarify the origin of the dichotomy, the key approach is to investigate the detailed distribution of the dark material. Recent studies of impact craters and surface temperature from Cassini spacecraft data implied that sublimation of H₂O ice can occur on Iapetus’ surface. This ice sublimation can change the albedo distribution on the moon with time.In this study, we evaluate the effect of ice sublimation and simulate the temporal change of surface albedo. We assume the dark material and the bright ice on the surface to be uniformly mixed with a certain volume fraction, and the initial albedo distribution to incorporate the dark material deposits on the surface. That is, the albedo at the apex is lowest and concentrically increases in a sinusoidal pattern. This situation simulates that dark materials existed around the Iapetus’ orbit billions of years ago, and the synchronously rotating Iapetus swept the material and then deposited it on its surface. The evolution of the surface albedo during 4.0 Gyr is simulated by estimating the surface temperature from the insolation energy on Iapetus including the effect of Saturn’s eccentricity and Iapetus’ obliquity precession, and evaluating the sublimation rate of H₂O ice from the Iapetus’ surface.As a result, we found that the distribution of the surface albedo changed dramatically after 4.0 Gyr of evolution. The sublimation has three important effects on the resultant surface albedo. First, the albedo in the leading hemisphere has significantly decreased to approach the minimum value. Second, the albedo distribution has been elongated along the equator. Third, the edge of the low albedo region has become clear. Considering the effect of ice sublimation, the current albedo distribution can be reconstructed from the sinusoidal albedo distribution, suggesting the apex-antapex cratering asymmetry as a candidate for the origin of the albedo dichotomy. From the model analysis, we obtained an important aspect that the depth of the turn-over layer where the darkening process proceeded for 4 Gyr should be an order of 10 cm, which is consistent with evaluation from the Cassini radar observations.     

Rational approximation formula for Chandrasekhar’s <Emphasis Type="Italic">H</Emphasis>-function for isotropic scattering

In this work, we first establish a simple procedure to obtain with 11-figure accuracy the values of Chandrasekhar’s H-function for isotropic scattering using a closed-form integral representation and the Gauss-Legendre quadrature. Based on the numerical values of the function produced by this method for various combinations of ? ₀, the single scattering albedo, and μ, the cosine of the zenith angle θ of the direction of radiation emergent from or incident upon a semi-infinite scattering-absorbing medium, we propose a rational approximation formula with μ ^1/4 and \(\sqrt{1-\varpi_{0}}\) as the independent variables. This allows us to reproduce the correct values of H(? ₀,μ) within a relative error of 2.1×10^?5 without recourse to any iterative procedure or root-finding process.     

Influence of the Oyashio Current and Tsugaru Warm Current on the circulation and water properties of Otsuchi Bay,Japan

Mooring and hydrographic observations were conducted from September 2012 to May 2014 at the mouth of Otsuchi Bay, a ria along the Pacific coast of Japan. Our observations quantitatively demonstrated that the circulation and the water properties of Otsuchi Bay are strongly influenced by the Tsugaru Warm Current (TWC) and Oyashio Current (OY) at seasonal and subseasonal time scales. Two bottom-mounted velocity profilers and temperature and salinity measurements beneath the near-surface halocline showed a counterclockwise lateral circulation pattern related to the TWC, which was enhanced from summer to autumn. From winter to early spring, the lateral circulation patterns related to the TWC weakened and the influence of the OY occasionally increased. When the OY was weak, surface flows became an overturning structure, with outflows in the upper layer and inflows in the lower layer. When the OY was strong and passed close to the Sanriku coast, the circulation became highly variable and intermittent. Intrusions of the markedly low-salinity OY water were observed on two occasions and persisted for periods of several weeks to several months. Salinity was sometimes less than 33.7, the lower limit of the typical TWC from late summer to autumn even when the TWC dominates. We suggested that this is the seasonal fluctuations of the TWC itself, as the upstream current of the Tsushima Warm Current is freshened in summer as a result of the influence of the Changjiang River. The surface water was generally fresher in the south of the bay than in the north, suggesting the Coriolis deflection of the river plume.     

A comparative simulation study of coupled THM processes and their effect on fractured rock permeability around nuclear waste repositories

This paper presents an international, multiple-code, simulation study of coupled thermal, hydrological, and mechanical (THM) processes and their effect on permeability and fluid flow in fractured rock around heated underground nuclear waste emplacement drifts. Simulations were conducted considering two types of repository settings (1) open emplacement drifts in relatively shallow unsaturated volcanic rock, and (2) backfilled emplacement drifts in deeper saturated crystalline rock. The results showed that for the two assumed repository settings, the dominant mechanism of changes in rock permeability was thermal–mechanically induced closure (reduced aperture) of vertical fractures, caused by thermal stress resulting from repository-wide heating of the rock mass. The magnitude of thermal–mechanically induced changes in permeability was more substantial in the case of an emplacement drift located in a relatively shallow, low-stress environment where the rock is more compliant, allowing more substantial fracture closure during thermal stressing. However, in both of the assumed repository settings in this study, the thermal–mechanically induced changes in permeability caused relatively small changes in the flow field, with most changes occurring in the vicinity of the emplacement drifts.     

Audio frequency magneto-telluric survey of Norikura Volcano in central Japan

Norikura Volcano has not been active during the last 10,000 years in spite of the activity of the surrounding volcanic mountains. To study past volcanic activities, geological studies were carried out extensively. However, quite a few geophysical investigations were conducted to contribute to volcanology. Our objective is to detect the present subsurface structure of Norikura Volcano and to define volcanic stratifications. In the vicinity of Norikura Volcano, geothermal fields are still active. Subsurface volcanic rocks in this area have been exposed to geothermal activity and altered. To comprehend volcanic stratifications of Norikura and geothermal activity, we conducted audio frequency magneto-telluric (AMT) surveys around Norikura Volcano. AMT survey is useful in clearly defining the resistivity structure related to volcanic regions. The AMT data were acquired over a frequency range 10 Hz–10 kHz. Decomposition analysis was applied to the tensor impedance data. Subsequently, apparent resistivity and phase data were inverted using a two-dimensional magneto-telluric (MT) inversion and a model of Norikura was derived. The final model manifests that the surface resistors are in agreement with andesite lava or dacite lava. As for the deeper structure, a horizontal conductor is situated above resistive basements. The alteration of the conductor was weak, while basement rocks were strongly altered and/or heated through the thermal activity. The existence of these layers seems to indicate the degree of thermal activity of Norikura Volcano.     

Redeposition of volcaniclastic sediments by a tsunami 4600 years ago at Kushima City,south-eastern Kyushu,Japan

The Hyuga-nada Sea, south-eastern Kyushu, Japan, is located between a strong (Nankai Trough) and a weak interplate coupling zone (Ryukyu Trench). Over the past 400 years this area has only experienced Magnitude 7·5 earthquakes or smaller and associated small-scale tsunamis. However, this short historical record most likely does not include the full range of high magnitude, low frequency giant earthquakes that might have occurred in the region. Thus, it is still unclear whether giant earthquakes and their associated tsunamis have occurred in this region. This paper reports on a prehistoric tsunami deposit discovered in a coastal lowland in south-eastern Kyushu facing the Hyuga-nada Sea. There is a reddish-brown pumiceous layer preserved in a non-marine, organic-rich mud sequence obtained from onshore sediment cores. This layer is recognized as the ca 4600 year old Kirishima-Miike tephra (that is now placed around 4500 years ago) sourced from Mount Kirishima, southern Kyushu. Another whitish pumiceous layer is evident below the Kirishima-Miike tephra in almost all of the sediment cores. A relatively high percentage of marine and brackish diatoms is recorded within this lower pumiceous layer (but not in the surrounding muds or in the overlying Kirishima-Miike tephra), indicating a marine or beach sediment source. Plant material obtained from organic-rich mud immediately below the event layer was dated to ca 4430 to 4710 cal yr bp , providing a limiting-maximum age for this marine incursion event. The presence of marine diatoms below the event layer is probably explained by pre-seismic subsidence. An absence of the resting spore of the planktonic brackish diatom Cheatoceros and the appearance of the freshwater diatom Eunotia serra immediately above the event layer probably represents a marked change to a relatively low-salinity environment. Assuming that there were no significant local geomorphological changes, such as drainage obstruction caused by formation of a new barrier spit, it is considered that co-seismic or immediate post-seismic uplift are the most likely explanations for this notable environmental change. Based on the crustal movements noted before and after the marine incursion, this event is interpreted here as an earthquake-generated tsunami. Moreover, because of these notable seismic crustal movements the tsunamigenic earthquake probably occurred immediately offshore of the study site.     

Water masses and currents of deep circulation southwest of the Shatsky Rise in the western North Pacific

We conducted full-depth hydrographic observations in the southwestern region of the Northwest Pacific Basin in September 2004 and November 2005. Deep-circulation currents crossed the observation line between the East Mariana Ridge and the Shatsky Rise, carrying Lower Circumpolar Deep Water westward in the lower deep layer (θ<1.2 °C) and Upper Circumpolar Deep Water (UCDW) and North Pacific Deep Water (NPDW) eastward in the upper deep layer (1.3–2.2 °C). In the lower deep layer at depths greater than approximately 3500 m, the eastern branch current of the deep circulation was located south of the Shatsky Rise at 30°24′–30°59′N with volume transport of 3.9 Sv (1 Sv=10⁶ m³ s⁻¹) in 2004 and at 30°06′–31°15′N with 1.6 Sv in 2005. The western branch current of the deep circulation was located north of the Ogasawara Plateau at 26°27′–27°03′N with almost 2.1 Sv in 2004 and at 26°27′–26°45′N with 2.7 Sv in 2005. Integrating past and present results, volume transport southwest of the Shatsky Rise is concluded to be a little less than 4 Sv for the eastern branch current and a little more than 2 Sv for the western branch current. In the upper deep layer at depths of approximately 2000–3500 m, UCDW and NPDW, characterized by high and low dissolved oxygen, respectively, were carried eastward at the observation line by the return flow of the deep circulation composing meridional overturning circulation. UCDW was confined between the East Mariana Ridge and the Ogasawara Plateau (22°03′–25°33′N) in 2004, whereas it extended to 26°45′N north of the Ogasawara Plateau in 2005. NPDW existed over the foot and slope of the Shatsky Rise from 29°48′N in 2004 and 30°06′N in 2005 to at least 32°30′N at the top of the Shatsky Rise. Volume transport of UCDW was estimated to be 4.6 Sv in 2004, whereas that of NPDW was 1.4 Sv in 2004 and 2.6 Sv in 2005, although the values for NPDW may be slightly underestimated, because they do not include the component north of the top of the Shatsky Rise. Volume transport of UCDW and NPDW southwest of the Shatsky Rise is concluded to be approximately 5 and 3 Sv, respectively. The pathways of UCDW and NPDW are new findings and suggest a correction for the past view of the deep circulation in the Pacific Ocean.     

Variation of the partial pressure of CO2 in surface water from Kuroshio to Oyashio and the relation between environmental factors and the partial pressure at 144°E off Sanriku, northwestern North Pacific in May, 1997

Partial pressure of CO₂ in surface sea water (pCO₂) was measured continuously off Sanriku in May, 1997 by a new pCO₂ measurement system. We have examined the relation of pCO₂ to physical factors such as temperature, salinity and density, chemical and biological factors such as nutrients and carbonate system and chlorophylla. In the Kuroshio region pCO₂ was not correlated to physical, chemical and biological factors in the range of 260 to 290 μatom. In transition water (Tr1) between Kuroshio and the Oyashio second branch, pCO₂ was weakly correlated to physical factors and strongly correlated to nutrients. In transition water (Tr2) between the Oyashio first and second branches, pCO₂ was highly correlated to temperature (SD: 10.9 μatom) and salinity (SD: 8.6 μatom) and also to nutrients. In transition water (Tr1+Tr2), pCO₂ was highly multivariately correlated to temperature (T), salinity (S), chlorophylla (CH) (or nitrate+nitrite (N)) as follows, pCO₂(μatom)= 10.8×T(°C)+27.7×S+2.57CH(μg/1) −769, R²= 0.86, SD = 20.9, or pCO₂(μatom)= 3.9×T(°C)+25.5×S+16.0NO₃(μM) −686, R²= 0.99, SD = 6.4. Moreover, pCO₂ was predicted by only two factors, one physical (S) and the other chemical/biological (N) as follows: pCO₂ (μatom)=32.8×S+19.4N−908, R²=0.97, SD=8.4. The pH measured at 25°C was well correlated with normalized pCO₂ at a fixed temperature. In the Oyashio region pCO₂ was decreased to 160 μatom, probably because of spring bloom, but was not correlated linearly to chlorophylla. The results obtained showed the possibility of estimating pCO₂ of the Oyashio and transition regions in May by satellite remote sensing of SST, but the problem of estimation of pCO₂ in Kuroshio water remains to be solved.     

A computerized geologic mapping system based on logical models of geologic structures

In this article we describe the basic framework of the computerized geologic mapping system cigma. The system, whic is based on a mathematical formulation of geologic concepts, consists of the following six subsystems: (1) input of geologic data set; (2) inference of stratigraphic sequence; (3) construction of logical models of geologic structures; (4) determination of three-dimensional geologic boundary surfaces; (5) construction of three-dimensional solid model of geologic structures; and (6) graphical presentation. Geologic structures are summarized in several tables called logical models of geologic structures. Each model is constructed automatically from input data on structural relations between geologic bodies. The model interprets the data automatically to create data files necessary to determine the shapes of geologic boundaries; it also provides a threedimensional solid model of geologic structures referring to the shapes of boundaries. As a prototype, we introduce two types of contacts corresponding to conformity and unconformity into the logical model and show that it is possible to draw a geologic map automatically. More complex geologic structures can be introduced into the geologic mapping system through further formulation of geologic structures.     

Experimental investigation of the α-β transformation of San Carlos olivine single crystal

Single crystalline San Carlos olivine (1 mm cube) was transformed to (Mg,Fe)₂SiO₄β-phase at 13.5–15 GPa, 1030–1330 °C for 0–600 min using a multi-anvil high pressure apparatus. The α-β transformation occurred by incoherent surface nucleation and interface-controlled growth and recovered partially transformed samples showed sharply defined reaction rim. The growth rate of the β-phase rim significantly decreased with time and the growth eventually ceased. TEM observations revealed that many dislocations were created in both the relict olivine just near the α-β interface and the β-phase in the rim, which show evidence for deformation caused by interfacial stresses associated with the misfit elastic strain of the transformation. The observed tangled dislocation texture in β-phase suggested that the β-phase rim was hardened and relaxation of the interfacial stress was retarded. This probably caused a localized pressure drop in the relict olivine and decreased the growth rate. Time-dependent growth rates of β-phase is possibly controlled by the rheology of β-phase, which must be considered for the prediction of the olivine metastability in the subducting slabs. Received: 24 January 1997 / Revised, accepted: 24 July 1998