Millennial-scale oscillations in the Kuroshio–Oyashio boundary during MIS 19 based on the radiolarian record from the Chiba composite section, central Japan

Marine isotope stage (MIS) 19 is considered to be the best orbital analog for the present interglacial. Consequently, clarifying the climatic features of this period can provide us with insights regarding a natural baseline for assessing future climate changes. A high-resolution radiolarian record from 800 to 750 ka (MIS 20 to MIS 18) was examined from the Chiba composite section (CbCS) of the Kokumoto Formation, including the Global Boundary Stratotype Section and Point for the lower–middle Pleistocene boundary on the Boso Peninsula on the Pacific side of central Japan. Millennial-scale oscillations in the Kuroshio warm and Oyashio cold currents were revealed by the Tr index, which is estimated using a simple equation based on radiolarian assemblages. The estimated Tr values ranged between 0.1 and 0.8 for MIS 18 through MIS 19, with minimum and maximum values corresponding to values observed off present day Aomori (41°N) and the Boso Peninsula (35°N), respectively. The observed patterns tended to be synchronous with the total radiolarian abundance associated with their production. Multiple maxima in radiolarian abundance occurred during periods of the Oyashio expanded mode before 785 ka and during periods of Kuroshio extension after 785 ka in MIS 19. Such increases in radiolarian abundance with the Kuroshio extension during MIS 19 are likely related to improvements in nutrient and photic environments with the development of a two-layer structure along the Kuroshio–Oyashio boundary zone. A similar pattern of millennial-scale climatic changes was also recognized in a precipitation record from the Sulmona Basin in central Italy, suggesting a close relationship with the CbCS record as a result of a large-scale climate system similar to the Arctic Oscillation in the northern hemisphere.

Effects of snow manipulation on larch trees in the taiga forest ecosystem in northeastern Siberia

Changes in winter precipitation (snow) may greatly affect vegetation by altering hydrological and biochemical processes. To understand the effects of changing snow cover depth and melt timing on the taiga forest ecosystem, a snow manipulation experiment was conducted in December 2015 at the Spasskaya Pad experimental larch forest in Eastern Siberia, which is characterized by a continental dry climate with extreme cold winters and hot summers. Variables including soil temperature and moisture, oxygen and hydrogen isotope ratios of soil moisture and stem water, foliar nitrogen and carbon contents and their isotopes, phenology, and soil inorganic nitrogen were observed at snow removal (SNOW−), snow addition (SNOW+), and CONTROL plots. After snow manipulation, the soil temperature at the SNOW− plot decreased significantly compared to the CONTROL and SNOW+ plots. At SNOW− plot, snowmelt was earlier and soil temperature was higher than at other plots during spring because of low soil moisture caused by less snowmelt water. Despite the earlier snowmelt and higher soil temperature in the SNOW− plot in the early growing season, needle elongation was delayed. Leaf chemistry also differed between the CONTROL and SNOW− plots. The needle nitrogen content in the SNOW− plot was lower in the middle of July, whereas no difference was observed among the three plots in August. The soil inorganic nitrogen content of each plot corresponded to these results. The amount of soil ammonium was lower in the SNOW− plot than in the other plots at the end of July, however, once production started in August, the amount of soil ammonium in the three plots was comparable. Extremely low soil temperatures in winter and freeze–thaw cycles in spring and dry soil condition in spring and early summer at the SNOW− plot may have influenced the phenology and production of soil inorganic nitrogen.

Detection of Martian dust storms using mask regional convolutional neural networks

Martian dust plays a crucial role in the meteorology and climate of the Martian atmosphere. It heats the atmosphere, enhances the atmospheric general circulation, and affects spacecraft instruments and operations. Compliant with that, studying dust is also essential for future human exploration. In this work, we present a method for the deep-learning-based detection of the areal extent of dust storms in Mars satellite imagery. We use a mask regional convolutional neural network, consisting of a regional-proposal network and a mask network. We apply the detection method to Mars daily global maps of the Mars global surveyor, Mars orbiter camera. We use center coordinates of dust storms from the eight-year Mars dust activity database as ground-truth to train and validate the method. The performance of the regional network is evaluated by the average precision score with $$50\%$$ 50 % overlap ($$mAP_{50}$$ m A P 50 ), which is around $$62.1\%$$ 62.1 % .

An electrodynamics model for Data Interpretation and Numerical Analysis of ionospheric Missions and Observations (DINAMO)

We introduce a new numerical model developed to assist with Data Interpretation and Numerical Analysis of ionospheric Missions and Observations (DINAMO). DINAMO derives the ionospheric electrostatic potential at low- and mid-latitudes from a two-dimensional dynamo equation and user-specified inputs for the state of the ionosphere and thermosphere (I–T) system. The potential is used to specify the electric fields and associated F-region E × B plasma drifts. Most of the model was written in Python to facilitate the setup of numerical experiments and to engage students in numerical modeling applied to space sciences. Here, we illustrate applications and results of DINAMO in two different analyses. First, DINAMO is used to assess the ability of widely used I–T climatological models (IRI-2016, NRLMSISE-00, and HWM14), when used as drivers, to produce a realistic representation of the low-latitude electrodynamics. In order to evaluate the results, model E × B drifts are compared with observed climatology of the drifts derived from long-term observations made by the Jicamarca incoherent scatter radar. We found that the climatological I–T models are able to drive many of the features of the plasma drifts including the diurnal, seasonal, altitudinal and solar cycle variability. We also identified discrepancies between modeled and observed drifts under certain conditions. This is, in particular, the case of vertical equatorial plasma drifts during low solar flux conditions, which were attributed to a poor specification of the E-region neutral wind dynamo. DINAMO is then used to quantify the impact of meridional currents on the morphology of F-region zonal plasma drifts. Analytic representations of the equatorial drifts are commonly used to interpret observations. These representations, however, commonly ignore contributions from meridional currents. Using DINAMO we show that that these currents can modify zonal plasma drifts by up to ~ 16 m/s in the bottom-side post-sunset F-region, and up to ~ 10 m/s between 0700 and 1000 LT for altitudes above 500 km. Finally, DINAMO results show the relationship between the pre-reversal enhancement (PRE) of the vertical drifts and the vertical shear in the zonal plasma drifts with implications for equatorial spread F.

Time difference between the 1854 CE Ansei–Tokai and Ansei–Nankai earthquakes estimated from distant tsunami waveforms on the west coast of North America

We estimated the time difference between the 1854 CE Ansei–Tokai and Ansei–Nankai earthquakes from tidal records of two tide gauge stations (San Francisco and San Diego) on the west coast of North America. The first signals of the Ansei–Tokai tsunami were apparent, whereas those of the Ansei–Nankai tsunami were obscured by the later waves of the Ansei–Tokai tsunami. Waveforms of the Ansei–Nankai tsunami simulated with nonlinear dispersive wave theory by assuming an origin time of 07:00 GMT on 24 December arrived earlier than in the observations. The normalized root mean square and the misfit between the simulated and observed waveforms of the Ansei–Nankai tsunami showed a time difference between them of approximately 0.4 h. This finding suggests that the actual origin time of the Ansei–Nankai tsunami was approximately 07:24 GMT on 24 December. A previous study estimated the origin time of the Ansei–Tokai tsunami to be about 00:30 GMT on 23 December. Thus, we concluded that the time difference between the 1854 CE Ansei–Tokai and Ansei–Nankai tsunamis was 30.9 h. Despite the significant difference in the time resolution between the seasonal timekeeping system used in Japan in 1854 and waveform digitization, our result is roughly in agreement with historical descriptions of the tsunamis, suggesting that such information can be effectively used to determine the origin times of historical earthquakes.

Formation, preservation and extinction of high-pressure minerals in meteorites: temperature effects in shock metamorphism and shock classification

The goal of classifying shock metamorphic features in meteorites is to estimate the corresponding shock pressure conditions. However, the temperature variability of shock metamorphism is equally important and can result in a diverse and heterogeneous set of shock features in samples with a common overall shock pressure. In particular, high-pressure (HP) minerals, which were previously used as a solid indicator of high shock pressure in meteorites, require complex pressure–temperature–time (P–T–t) histories to form and survive. First, parts of the sample must be heated to melting temperatures, at high pressure, to enable rapid formation of HP minerals before pressure release. Second, the HP minerals must be rapidly cooled to below a critical temperature, before the pressure returns to ambient conditions, to avoid retrograde transformation to their low-pressure polymorphs. These two constraints require the sample to contain large temperature heterogeneities, e.g. melt veins in a cooler groundmass, during shock. In this study, we calculated shock temperatures and possible P–T paths of chondritic and differentiated mafic–ultramafic rocks for various shock pressures. These P–T conditions and paths, combined with observations from shocked meteorites, are used to constrain shock conditions and P–T–t histories of HP-mineral bearing samples. The need for rapid thermal quench of HP phases requires a relatively low bulk-shock temperature and therefore moderate shock pressures below ~ 30 GPa, which matches the stabilities of these HP minerals. The low-temperature moderate-pressure host rock generally shows moderate shock-deformation features consistent with S4 and, less commonly, S5 shock stages. Shock pressures in excess of 50 GPa in meteorites result in melt breccias with high overall post-shock temperatures that anneal out HP-mineral signatures. The presence of ringwoodite, which is commonly considered an indicator of the S6 shock stage, is inconsistent with pressures in excess of 30 GPa and does not represent shock conditions different from S4 shock conditions. Indeed, ringwoodite and coexisting HP minerals should be considered as robust evidence for moderate shock pressures (S4) rather than extreme shock (S6) near whole-rock melting.

Evolution of the geological structure and mechanical properties due to the collision of multiple basement topographic highs in a forearc accretionary wedge: insights from numerical simulations

We propose a conceptual geological model for the collision of multiple basement topographic highs (BTHs; e.g., seamounts, ridges, and horsts) with a forearc accretionary wedge. Even though there are many BTHs on an oceanic plate, there are few examples of modeling the collision of multiple BTHs. We conducted numerical simulations using the discrete element method to examine the effects of three BTH collisions with forearcs. The typical geological structure associated with a BTH collision was reproduced during the collision of the first BTH, and multiple BTH collisions create a cycle of formation of BTH collisional structures. Each BTH forces the basal décollement to move up to the roof décollement, and the roof décollement becomes inactive after the passage of the BTH, and then the décollement moves down to the base. As the active décollement position changes, the sequences of underthrust sediments and uplifted imbricate thrusts are sandwiched between the décollements and incorporated into the wedge. At a low horizontal compressive stress, a “shadow zone” is formed behind (i.e., seaward of) the BTH. When the next BTH collides, the horizontal compressive stress increases and tectonic compaction progresses, which reduce the porosity in the underthrust sediments. Heterogeneous evolution of the geological and porosity structure can generate a distinctive pore pressure pattern. The underthrust sediments retain fluid in the “shadow” of the BTH. Under the strong horizontal compressive stresses associated with the next BTH collision, pore pressure increases along with a rapid reduction of porosity in the underthrust sediments. The distinctive structural features observed in our model are comparable to the large faults in the Kumano transect of the Nankai Trough, Japan, where a splay fault branches from the plate boundary and there are old and active décollements. A low-velocity and high-pore-pressure zone is located at the bottom of the accretionary wedge and in front (i.e., landward) of the subducting ridge in the Kumano transect. This suggests that strong horizontal compressive stresses associated with the current BTH collision has increased the pore pressure within the underthrust sediments associated with previous BTHs.

Frictional heating in a thick fault zone with empirical slip-weakening friction: implications for slip parameter estimation from temperature observations in deep fault drilling

The strain energy released during an earthquake is consumed by processes related to seismic radiation or dissipation. Deep fault drilling and subsequent temperature measurements in a thick fault zone immediately after an event have provided important insights into this dissipation process. By employing an analytical solution to the heat conduction problem, which involves the sudden injection of an infinitesimally thin heat source into an infinite medium, previous drilling projects have estimated the strength of the heat source and the level of shear stress from observed temperature anomalies. However, it is unclear under what conditions this analytical source solution can be regarded as a good approximation for the thick fault problem, a situation which has led to uncertainty of the approximation error in these previous studies. In this study, I first derived an analytical solution for the thick fault problem that accounted for experimentally derived slip-weakening friction. I then validated the derived solution both analytically and numerically. Using the derived thick solution, I next demonstrated that the thick, planar, and source solutions can be considered equivalent under the typical conditions of the previous drilling projects. Therefore, the slip parameters estimated by using the source solution obtained by these studies are appropriate. These results suggest that coseismic information with spatio-temporal extent, such as shear stress and friction coefficient, are lost due to heat diffusion when the temperature observations are conducted; thus, they cannot be inferred directly from observed temperature anomalies. These results also suggest that for most drilling projects, including future ones, the observed temperature distribution can be well explained by using the source solution instead of the thick solution as long as coseismic slip is not markedly delocalized and the spatial extent of the temperature measurements is not significantly larger than the diffusion length.

Classification of imbalanced cloud image data using deep neural networks: performance improvement through a data science competition

Image data classification using machine learning is an effective method for detecting atmospheric phenomena. However, extreme weather events with a small number of cases cause a decrease in classification prediction accuracy owing to the imbalance in data between the target class and the other classes. To build a highly accurate classification model, I held a data analysis competition to determine the best classification performance for two classes of cloud image data, specifically tropical cyclones including precursors and other classes. For the top models in the competition, minority data oversampling, majority data undersampling, ensemble learning, deep layer neural networks, and cost-effective loss functions were used to improve the classification performance of the imbalanced data. In particular, the best model of 209 submissions succeeded in improving the classification capability by 65.4% over similar conventional methods in a measure of the low false alarm ratio.

Changes in surface water masses in the northern East China Sea since the Last Glacial Maximum based on diatom assemblages

We investigated diatom assemblages for the last 19 ka in a piston core KY07-04 PC01 collected from the northern Okinawa Trough, East China Sea. Most of the diatom taxa in the Core KY07-04 PC01 were subtropical to temperate marine species throughout the core. Pronounced turnover of major diatom groups from periphytic to planktonic taxa occurred during deglaciation between 17 and 14 ka ago. During the Last Glacial Maximum (LGM) to early deglaciation, the representative taxa were Paralia spp., a periphytic diatom abundant in continental shelves of the modern East China Sea covering the Chinese Coastal Water. Relative abundances of Paralia spp. rapidly decreased during 15–14 ka ago, corresponding to sea-level rise, displaced by planktonic taxa such as Thalassionema. Thalassionema nitzschioides varieties living in subtropical water increased since 12 ka and have reached the present level at 8 ka ago, suggesting strengthening the influences of Kuroshio water in the northern East China Sea during the early to middle Holocene. Warm planktonic diatom taxa such as Nitzschia bicapitata also gradually increased from 12 to 8 ka ago. Near the core top, a sudden increase in N. bicapitata might relate to the intensifying Kuroshio during the latest Holocene.