Palaeointensities of the Auckland geomagnetic excursions by the LTD-DHT Shaw method

We report new palaeointensity results concerning the Auckland geomagnetic excursions using the double heating technique of the Shaw method with low temperature demagnetisation (LTD-DHT Shaw method). The excursional palaeodirections recorded in six volcanoes of the Auckland volcanic field, New Zealand, have been classified into three groups: north-down (ND), west-up (WU) and south-up (SU) directions. In the present study, five to six consistent palaeointensities have been obtained from each of five volcanoes recording the Auckland geomagnetic excursions. The Wiri (27 ka), Crater Hill and Puketutu volcanoes (ND group) yielded mean palaeointensities of 10.6 ± 1.2 (1σ), 11.8 ± 2.8 and 11.1 ± 0.4 μT, respectively. The Hampton Park volcano (55 ka; WU group) gave 9.5 ± 1.2 μT while the McLennan Hills volcano (SU group) gave 2.5 ± 0.5 μT. It is notable that consistent palaeointensities have been obtained from the three different volcanoes which have almost the same palaeodirections (ND group), possibly supporting the reliability of the palaeointensity data. These five palaeointensities for the Auckland geomagnetic excursions correspond to virtual dipole moments (VDMs) of 0.6-2.1 × 10²² A m², whereas three mean palaeointensities obtained from the Auckland volcanoes having non-excursional palaeodirections are 13.1-40.0 μT giving stronger VDMs of 2.1-6.9 × 10²² A m². These results suggest that the dipole component of the geomagnetic field reduced to about 2 × 10²² A m² or less during the Auckland geomagnetic excursions.     

Paleotectonic structures and their influence on recent seismo-tectonics in the south Kuril subduction zone

Abstract Bathymetric data from south of Hokkaido obtained during a cruise of R/V Hakuho-Maru are summarized, and their correlation with earthquake occurrence is discussed. There are structural lineations on the seaward slope of the Kuril Trench, oblique to the Kuril Trench axis and parallel to the magnetic lineations in the Pacific plate. The structural lineations comprise horst-grabens generated by normal faulting. This suggests that Cretaceous tectonic structures originating at the spreading centre affect present seismotectonics around the trench axis. The structural-magnetic relation is compared to the case of the Japan Trench. North-east of the surveyed area, there are two major fracture zones (Nosappu Fracture Zone and Iturup Fracture Zone) that divide the oceanic plate into three segments. If the fracture zones (FZ) and the zone of paleo-mechanical weakness, represented by magnetic lineations, can control the direction of normal faults at a trench, the extent of the resulting topographic roughness on the seaward slope of the trench would be different across an FZ because of the differences in ages. By studying recent large earthquakes occurring in the south Kuril region, it is shown that several main-aftershock distributions for large earthquakes in this region are bounded by the Nosappu FZ and the Iturup FZ. Two models (Barrier model and Rebound model) are presented to interpret earthquake occurrence near the south Kuril Islands. The Barrier model explains seismic boundaries seen in several examples for earthquake occurrence in the south Kuril regions. The fracture zone forming the boundary of two segments with different magnetic lineations is also the boundary of two different normal fault systems on their ocean bottom, and the difference in sea-bottom roughness between two normal fault systems should affect the seismic coupling at a plate interface. Due to the difference of seismic coupling, earthquake occurrence is controlled by an FZ and then the FZ acts as a seismic boundary (Barrier model). Existing normal faults created by plate bending of subducting oceanic plate should rebound after its subduction (Rebound model). This rebound of normal faults may cause intraplate earthquakes with a high-angle reverse-fault mechanism such as the 1994 Shikotan Earthquake. The energy released by an intraplate earthquake generated by normal-fault rebounding is not directly related to that of interplate earthquakes such as low-angle thrust earthquakes. It is a reason why large earthquakes occurred in the same region during a relatively short period.     

Assessing the suitability of data from Sentinel-1A and 2A for crop classification

Sentinel-1A C-SAR and Sentinel-2A MultiSpectral Instrument (MSI) provide data applicable to the remote identification of crop type. In this study, six crop types (beans, beetroot, grass, maize, potato, and winter wheat) were identified using five C-SAR images and one MSI image acquired during the 2016 growing season. To assess the potential for accurate crop classification with existing supervised learning models, the four different approaches namely kernel-based extreme learning machine (KELM), multilayer feedforward neural networks, random forests, and support vector machine were compared. Algorithm hyperparameters were tuned using Bayesian optimization. Overall, KELM yielded the highest performance, achieving an overall classification accuracy of 96.8%. Evaluation of the sensitivity of classification models and relative importance of data types using data-based sensitivity analysis showed that the set of VV polarization data acquired on 24 July (Sentinel-1A) and band 4 data (Sentinel-2A) had the greatest potential for use in crop classification.     

Responses of diatoms to iron-enrichment (SEEDS) in the western subarctic Pacific, temporal and spatial comparisons

Phytoplankton species composition was analyzed inside and outside of the iron-enriched patch during the SEEDS experiment. Before the iron-enrichment, the phytoplankton community consisted of similar proportions of pico-, nano- and micro-sized phytoplankton, and the micro-phytoplankton was dominated by the pennate diatom Pseudo-nitzschia turgidula. Although all the diatoms, except the nano-sized Fragilariopsis sp., increased during the two weeks of the observation period, the flora in the patch dramatically changed with the increase of phytoplankton biomass to a centric diatom-dominated community. Neritic diatoms, especially Chaetoceros debilis, showed higher growth rates than other diatoms, without any delay in the initiation of growth after the enrichment, and accounted for 90% of the micro-phytoplankton after day 9. In contrast, the oceanic diatoms showed distinct delays in the initiation of growth. We conclude that the responses of the diatoms to the manipulation of iron concentration were different by species, and the fast and intensive response of the phytoplankton to iron-enrichment resulted from the presence of a small amount of neritic diatoms at the study site. The important factors that determine the dominant species in the bloom are the potential growth rates under an iron-replete condition and the growth lag. Abundant species in the patch are widely distributed in the North Pacific and their relative contributions in the Oyashio area and at Stn KNOT are high from spring to summer. However, a characteristic difference of species composition between the SEEDS bloom and natural blooms was the lack of Thalassiosira and Coscinodiscus species in the patch, which usually account for a major part of the phytoplankton community under blooming conditions in the western North Pacific.     

Anisotropic tomography of P-wave traveltimes

The ray path of a P -wave is specified in terms of the ray parameter and three Euler angles. the P -wave traveltime depends only on the ray parameter for a spherically symmetric earth. If we introduce an aspherical perturbation, including general ani-sotropy, the dependence on Euler angles can be expanded in terms of the rotation matrix for a fixed ray parameter. If the perturbation is isotropic, the expansion coefficients satisfy certain relations which may be used to obtain definite evidence for anisotropy rather than isotropic lateral heterogeneity.     

A possible role for unstable coupled waves affected by resonance between Kelvin waves and seasonal warming in the development of the strong 1997–1998 El Niño

Time-varying air–sea coupled processes in the central to eastern equatorial Pacific associated with strong El Niño development during the 1997–1998 period are examined using a newly developed reanalysis dataset obtained from four-dimensional variational ocean–atmosphere coupled data assimilation experiments. The time series of this data field exhibits realistic features of El Niño evolution. Our analysis indicates that resonance between eastward-propagating oceanic downwelling Kelvin waves and the seasonal rise of sea-surface temperature (SST) in the central to eastern equatorial Pacific generates relatively persistent high SST conditions accompanied by a deeper thermocline and more relaxed easterly winds than usual. The surface condition resulting from the wave-seasonal SST resonance represents a preconditioned state that leads to an enhancement in incident downwelling Kelvin waves to levels sufficient to induce large-amplitude unstable coupled waves in the central to eastern equatorial region. Heat balance estimates using our reanalysis dataset suggest that the unstable coupled waves are categorized within the intermediate regime of coupled Kelvin and Rossby waves and have the potential to grow rapidly. We argue that the seasonal resonance and the unstable coupled waves should play crucial roles in the development of the largest historical El Niño event, which was recorded between late 1997 and early 1998.     

Evaluating metrics derived from Landsat 8 OLI imagery to map crop cover

Developing techniques are required to generate agricultural land cover maps to monitor agricultural fields. Landsat 8 Operational Land Imager (OLI) offers reflectance data over the visible to shortwave-infrared range. OLI offers several advantages, such as adequate spatial and spectral resolution, and 16 day repeat coverage, furthermore, spectral indices derived from Landsat 8 OLI possess great potential for evaluating the status of vegetation. Additionally, classification algorithms are essential for generating accurate maps. Recently, multi-Grained Cascade Forest, which is also called deep forest, was proposed, and it was shown to give highly competitive performance for classification. However, the ability of this algorithm to generate crop maps with satellite data had not yet been evaluated. In this study, the reflectance at 7 bands and 57 spectral indices calculated from Landsat 8 OLI data were evaluated for its potential for crop type identification.     

Timeline analysis and wavelet multiscale analysis of the AKARI All-Sky Survey at 90 μm

We present a careful analysis of the point-source detection limit of the AKARI All-Sky Survey in the WIDE-S 90-μm band near the North Ecliptic Pole (NEP). Timeline analysis is used to detect IRAS ( Infrared Astronomy Satellite ) sources and then a conversion factor is derived to transform the peak timeline signal to the interpolated 90-μm flux of a source. Combined with a robust noise measurement, the point-source flux detection limit at signal-to-noise ratio  (S/N) > 5  for a single detector row is  1.1 ± 0.1 Jy  which corresponds to a point-source detection limit of the survey of ∼0.4 Jy.
Wavelet transform offers a multiscale representation of the Time Series Data ( tsd ). We calculate the continuous wavelet transform of the tsd and then search for significant wavelet coefficients considered as potential source detections. To discriminate real sources from spurious or moving objects, only sources with confirmation are selected. In our multiscale analysis, IRAS sources selected above 4σ can be identified as the only real sources at the Point Source Scales. We also investigate the correlation between the non- IRAS sources detected in timeline analysis and cirrus emission using wavelet transform and contour plots of wavelet power spectrum. It is shown that the non- IRAS sources are most likely to be caused by excessive noise over a large range of spatial scales rather than real extended structures such as cirrus clouds.