Climate change during the last 150 million years: reconstruction from a carbon cycle model

Variations of the atmospheric CO₂ level and the global mean surface temperature during the last 150 Ma are reconstructed by using a carbon cycle model with high-resolution input data. In this model, the organic carbon budget and the CO₂ degassing from the mantle, both of which would characterize the carbon cycle during the Cretaceous, are considered, and the silicate weathering process is formulated consistently with an abrupt increase in the marine strontium isotope record for the last 40 Ma. The second-order variations of the atmospheric CO₂ level and the global mean surface temperature in addition to the first-order cooling trend are obtained by using high-resolution data of carbon isotopic composition of marine limestone, seafloor spreading rate, and production rate of oceanic plateau basalt. The results obtained from this model are in good agreement with the previous estimates of palaeo-CO₂ level and palaeoclimate inferred from geological, biogeochemical, and palaeontological models and records. The system analyses of the carbon cycle model to understand the cause of the climate change show that the dominant controlling factors for the first-order cooling trend of climate change during the last 150 Ma are tectonic forcing such as decrease in volcanic activity and the formation and uplift of the Himalayas and the Tibetan Plateau, and, to a lesser extent, biological forcing such as the increase in the soil biological activity. The mid-Cretaceous was very warm because of the high CO₂ level (4–5 PAL) maintained by the enhanced CO₂ degassing rate due to the increased mantle plume activities and seafloor spreading rates at that time, although the enhanced organic carbon burial would have a tendency to decrease the CO₂ level effectively at that period. The variation of organic carbon burial rate may have been responsible for the second-order climate change during the last 150 Ma.     

Incipient subduction and deduction along the eastern margin of the Japan Sea

The eastern margin of the Japan Sea has been under convergent tectonics since sometime after the Pliocene. The suture line is on the ocean-continent boundary between the Japan-Yamato Basins and the Tohoku (Northeast Japan) Arc. Incipient subduction, and obduction, of the oceanic and sub-oceanic crust are observed along the suture line with the occurrence of numerous N-S trending, fault-bounded ridges and troughs. Focal mechanism solutions of several compressional type earthquakes with magnitudes of M = 6.9 to 7.7 which have occurred along the zone are consistent with the observation of geological thrust faults. The convergent stresses in the Japan Sea may be due to India-Eurasia collision and its associated intra-plate or inter-microplate movements in East Asia. The eastward movement of the Amurian Plate causes Baikal extension along its western margin and the Japan Sea convergence along its eastern margin.     

Cohesive zone length of metagabbro at supershear rupture velocity

We investigated the shear strain field ahead of a supershear rupture. The strain array data along the sliding fault surfaces were obtained during the large-scale biaxial friction experiments at the National Research Institute for Earth Science and Disaster Resilience. These friction experiments were done using a pair of meter-scale metagabbro rock specimens whose simulated fault area was 1.5 m?×?0.1 m. A 2.6-MPa normal stress was applied with loading velocity of 0.1 mm/s. Near-fault strain was measured by 32 two-component semiconductor strain gauges installed at an interval of 50 mm and 10 mm off the fault and recorded at an interval of 1 MHz. Many stick-slip events were observed in the experiments. We chose ten unilateral rupture events that propagated with supershear rupture velocity without preceding foreshocks. Focusing on the rupture front, stress concentration was observed and sharp stress drop occurred immediately inside the ruptured area. The temporal variation of strain array data is converted to the spatial variation of strain assuming a constant rupture velocity. We picked up the peak strain and zero-crossing strain locations to measure the cohesive zone length. By compiling the stick-slip event data, the cohesive zone length is about 50 mm although it scattered among the events. We could not see any systematic variation at the location but some dependence on the rupture velocity. The cohesive zone length decreases as the rupture velocity increases, especially larger than \( \sqrt{2} \) times the shear wave velocity. This feature is consistent with the theoretical prediction.     

Age of the Solomon Sea Basin from magnetic lineations

Magnetic anomalies measured in the central to western half of the Solomon Sea, when considered with other magnetic data, reveal the existence of linear patterns. Magnetic lineation anomaly models of the Cenozoic, 65 to 0 Ma, suggest that an age between 34 and 28 Ma and a half-rate spreading speed of 5.8 cm/yr for the northern flank of a former spreading center best fits our present magnetic data in the Solomon Sea Basin. Heat flow and bathymetry data support this preferred model.     

P - V - T equation of state of stishovite to the mantle transition zone conditions

In-situ synchrotron X-ray diffraction experiments were conducted using the SPEED-1500 multi-anvil press of SPring-8 on stishovite SiO₂ and pressure-volume-temperature data were collected at up to 22.5 GPa and 1,073 K, which corresponds to the pressure conditions of the base of the mantle transition zone. The analysis of room-temperature data yielded V₀=46.56(1) ų, K_{T 0}=296(5) GPa and K _T =4.2(4), and these properties were consistent with the subsequent thermal equation of state (EOS) analyses. A fit of the present data to high-temperature Birch-Murnaghan EOS yielded (K_T /T) _P =–0.046(5) GPa K^–1 and = a + bT with values of a =1.26(11)×10^–5 K^–1 and b =1.29(17)×10^–8 K^–2. A fit to the thermal pressure EOS gives ₀=1.62(9)×10^–5 K^–1, ( K _T / T) _V =–0.027(4) GPa K^–1 and (²P /T ²) _V =27(5)×10^–7 GPa K^–2. The lattice dynamical approach by Mie-Grüneisen-Debye EOS yielded ₀=1.33(6), q =6.1(8) and ₀=1160(120) K. The strong volume dependence of the thermal pressure of stishovite was revealed by the analysis of present data, which was not detectable by the previous high-temperature data at lower pressures, and this yields ( K _T / T) _V 0 and q 1. The analyses for the fictive volume for a and c axes show that relative stiffness of c axis to a axis is similar both on compression and thermal expansion. Present EOS enables the accurate estimate of density of SiO₂ in the deep mantle conditions.     

PSP-toxicification of the carnivorous gastropod Rapana venosa inhabiting the estuary of Nikoh River, Hiroshima Bay, Hiroshima Prefecture, Japan

During surveillance on the toxicity of invertebrates such as bivalves inhabiting the coasts of Hiroshima Bay in 2001 and 2002, the carnivorous gastropod rapa whelk Rapana venosa collected in the estuary of Nikoh River, was found to contain toxins which showed paralytic actions in mice; the maximum toxicities (as paralytic shellfish poison, PSP) were 4.2 MU/g (May 2001) and 11.4 MU/g (April 2002). Their total toxicities were 224 and 206 MU/viscera of one specimen throughout the monitoring period. Attempts were made to identify the toxic principle in the gastropod. The viscera were extracted with 80% ethanol acidified with acetic acid, followed by defatting with dichloromethane. The aqueous layer obtained was treated with activated charcoal and then applied to a Sep-Pak C18 cartridge. The unbound toxic fraction was analyzed by high-performance liquid chromatography techniques. The gastropod toxin was rather unexpectedly identified as PSP. It was comprised of high toxic component (gonyautoxin-3; GTX3, GTX2, saxitoxin; STX) as the major components, which accounted for approximately 91 mol% of all components along with C1 and C2, which are N-sulfocarbamoyl derivatives. Judging from their toxin patterns, it is suggested that the PSP toxification mechanism of the gastropod that PSP toxins produced by phytoplankton such as Alexandrium tamarense, are transferred to and accumulated in plankton feeders such as the short-necked clam, and then transferred to this carnivorous rapa whelk R. venosa through predation.     

Age estimation of the Solomon Sea based on heat flow data

Several heat flow measurements were made during the NAT83 cruise in the central part of the Solomon Sea Basin. The average value of 87 mW/m² (2.08 HFU) calculated from these and other data indicates that the age of the Solomon Sea Basin may range from 24 to 44 Ma. This is supported by the water depth, of approximately 4,500 m, versus age relationship. There is a possibility that the Solomon Sea Basin is not a back-arc basin associated with an arc but was formerly a relatively large oceanic plate. The agreement in age from both heat flow and water depth data favors the latter hypothesis.     

Age of the Solomon Sea Basin from magnetic lineations

Magnetic anomalies measured in the central to western half of the Solomon Sea, when considered with other magnetic data, reveal the existence of linear patterns. Magnetic lineation anomaly models of the Cenozoic, 65 to 0 Ma, suggest that an age between 34 and 28 Ma and a half-rate spreading speed of 5.8 cm/yr for the northern flank of a former spreading center best fits our present magnetic data in the Solomon Sea Basin. Heat flow and bathymetry data support this preferred model.     

Stress fields and fault reactivation angles of the 2000 western Tottori aftershocks and the 2001 northern Hyogo swarm in southwest Japan

We estimated the stress fields of the aftershocks of the 2000 western Tottori earthquake (M_w 6.6) and the northern Hyogo swarm (max M_w 5.2) by a stress tensor inversion of moment tensor solutions reported from the National Research Institute for Earth Science and Disaster Prevention (Japan). The maximum principal stress direction of the western Tottori sequence was estimated as N107°E with a strike–slip regime. In the northern Hyogo swarm, the orientations of the principal stress directions could not be well constrained by the observed data, but after examining the detailed characteristics of the solution, we obtained a most probable solution of N113°E for the σ₁ direction. These solutions are consistent with the maximum horizontal directions roughly estimated from the strike directions of large earthquakes occurring geographically between these two seismic activities. We measured the angle between each fault–slip direction and maximum principal stress direction to investigate the frictional properties of earthquakes. The distribution of the angles was forward modeled to estimate the coefficient of friction and the stress ratio, assuming uniformly distributed fault orientations. For the western Tottori sequence, a homogeneous stress field with a coefficient of friction less than 0.4 was estimated. A high stress level was also suggested because very little change occurred in the stress field during the mainshock. For the northern Hyogo sequence, the coefficient of friction was estimated to be between 0.5 and 1.0.