Using multispectral landsat and sentinel-2 satellite data to investigate vegetation change at Mount St. Helens since the great volcanic eruption in 1980

Long-term analyses of vegetation succession after catastrophic events are of high interest for an improved understanding of succession dynamics. However, in many studies such analyses were restricted to plot-based measurements. Contrarily, spatially continuous observations of succession dynamics over extended areas and time-periods are sparse. Here, we applied a change vector analysis (CVA) to investigate vegetation succession dynamics at Mount St. Helens after the great volcanic eruption in 1980 using Landsat. We additionally applied a supervised random forest classification using Sentinel-2 data to map the currently prevailing vegetation types. Change vector analysis was performed with the normalized difference vegetation index (NDVI) and the urban index (UI) for three subsequent decades after the eruption as well as for the whole observation time between 1984 and 2016. The influence of topography on the current vegetation distribution was examined by comparing altitude, slope angles and aspect values of vegetation classes derived by the random forest classification. Wilcox- Rank-Sum test was applied to test for significant differences between topographic properties of the vegetation classes inside and outside of the areas affected by the eruption. For the full time period, a total area of 516 km² was identified as re-vegetated, whereas the area and magnitude of re-growing vegetation decreased during the three decades and migrated closer to the volcanic crater. Vegetation losses were mainly observed in regions unaffected by the eruption and related mostly to timber harvesting. The vegetation type classification reached a high overall accuracy of approximately 90%. 36 years after the eruption, coniferous and deciduous trees have established at formerly devastated areas dominating with a proportion of 66%, whereas shrubs are more abundant in riparian zones. Sparse vegetation dominates at regions very close to the crater. Elevation was found to have a great influence on the reestablishment and distribution of the vegetation classes within the devastated areas showing in almost all cases significant differences in altitude distribution. Slope was less important for the different classes - only representing significantly higher values for meadows, whereas aspect seems to have no notable influence on the reestablishment of vegetation at Mount St. Helens. We conclude that major vegetation succession dynamics after catastrophic events can be assessed and characterized over large areas from freely available remote sensing data and hence contribute to an improved understanding of succession dynamics.     

Modeling regional landslide susceptibility using dynamic soil moisture profiles

A landslide susceptibility mapping study was performed using dynamic hillslope hydrology. The modified infinite slope stability model that directly includes vadose zone soil moisture (SM) was applied at Cleveland Corral, California, US and Krishnabhir, Dhading, Nepal. The variable infiltration capacity (VIC-3L) model simulated vadose zone soil moisture and the wetness index hydrologic model simulated groundwater (GW). The GW model predictions had a 75% NASH-Sutcliffe efficiency when compared to California’s in-situ GW measurements. The model performed best during the wet season. Using predicted GW and VIC-3L vadose zone SM, the developed landslide susceptibility maps showed very good agreement with mapped landslides at each study region. Previous quasi-dynamic model predictions of Nepal’s hazardous areas during extreme rainfall events were enhanced to improve the spatial characterization and provide the timing of hazardous conditions.     

Effects of cold eddy on phytoplankton production and assemblages in Luzon strait bordering the South China Sea

The biochemical effects of a cold-core eddy that was shed from the Kuroshio Current at the Luzon Strait bordering the South China Sea (SCS) were studied in late spring, a relatively unproductive season in the SCS. The extent of the eddy was determined by time-series images of SeaWiFS ocean color, AVHRR sea surface temperature, and TOPEX/Jason-1 sea surface height anomaly. Nutrient budgets, nitrate-based new production, primary production, and phytoplankton assemblages were compared between the eddy and its surrounding Kuroshio and SCS waters. The enhanced productivity in the eddy was comparable to wintertime productivity in the SCS basin, which is supported by upwelled subsurface nitrate under the prevailing Northeastern Monsoon. There were more Synechococcus, pico-eucaryotes, and diatoms, but less Trichodesmium in the surface water inside the eddy than outside. Prochlorococcus and Richelia intracellularis showed no spatial differences. Water column-integrated primary production (IPP) inside the eddy was 2–3 times that outside the eddy in the SCS (1.09 vs. 0.59 g C m⁻²d⁻¹), as was nitrate-based new production (INP) (0.67 vs. 0.25 g C m⁻²d⁻¹). INP in the eddy was 6 times that in the Kuroshio (0.12 g C m⁻²d⁻¹). IPP and INP in the eddy were higher than the maximum production values ever measured in the SCS basin. Surface chlorophyll a concentration (0.40 mg m⁻³) in the eddy equaled the maximum concentration registered for the SCS basin and was higher than the wintertime average (0.29 ± 0.04 mg m⁻³). INP was 3.5 times as great and IPP was doubled in the eddy compared to the wintertime SCS basin. As cold core eddies form intermittently all year round as the Kuroshio invades the SCS, their effects on phytoplankton productivity and assemblages are likely to have important influences on the biogeochemical cycle of the region.     

Arctic sea-ice motion and its relation to pressure field

Daily Arctic sea-ice motion maps during the winter seasons (December–March) from December 1988 to March 2003 derived from NSCAT, QuikSCAT, SSM/I, and AMSRE data by a wavelet analysis method have been merged with those derived from buoy data. These merged sea-ice motion data have been used to study the circulation regimes and winter-to-winter variability of Arctic sea-ice motion. The relation between sea-ice motion and the pressure field in the Arctic Ocean was also studied by applying Principal Component Analysis (PCA) to the monthly merged sea-ice motion data and the monthly pressure field data from IABP. The mean Arctic sea-ice motion map of the 15 winter seasons has two distinct features: the Beaufort Gyre and a cyclonic circulation system in the Eurasian Basin, which moves ice from the Laptev Sea to Fram Strait. The strengths and sizes of the two features change from one winter season to another. Seasons with a strong or normal Beaufort Gyre alternate with seasons with a weak or no Beaufort Gyre every one to three seasons. The principal components of the first two modes of PCA of the monthly sea-ice motion are closely correlated with their counterparts of the monthly pressure field in the Arctic Ocean. The mode-one components of these two anomalies alternate between anticyclonic and cyclonic circulation systems. The correlation between Arctic Oscillation indexes and the principal components of the first mode of PCA of the monthly Arctic sea-ice motion is low but statistically significant.     

Coastal upwelling along the north coast of Papua New Guinea and SST cooling over the pacific warm pool: A case study for the 2002/03 El Niño event

We investigate an overlooked mechanism—coastal upwelling—for sea surface temperature (SST) cooling in the western side of the mean location of the Pacific warm pool (WSWP: 5°S–5°N, 140°E–150°E) prior to El Niño onset. We analyze various observed data such as the TRIangle Trans-Ocean buoy Network (TRITON) moored buoy data, Conductivity-Temperature-Depth (CTD) data, satellite data and a hindcast experiment output by a high-resolution ocean general circulation model (OGCM). We focus on the precondition of the 2002/03 El Niño event, for which many datasets are available. Relatively cool water upwelled along the north coast of Papua New Guinea (PNG) during December 2001, prior to the onset of the 2002/03 El Niño event, and then spread out over a wider area to the northeast. Simultaneously, strong west-northerly surface winds occur along the north coast. Heat budget analysis of TRITON buoy data in the WSWP reveals that negative zonal heat advection due to eastward current is the main factor for cooling the mixed layer in the WSWP in contrast to the warming effect of the surface heat flux during the period. This cooling requires a source of colder water to the west. Similar analysis of OGCM outputs also suggests that the upwelled relatively cool water along the PNG north coast, and its northeastward extension to the equatorial region, contributes to cooling of the surface water over the WSWP mainly via negative zonal heat advection. Similar mechanisms are confirmed also for the 1982/83 and 1997/98 El Niño events by analyses of OGCM outputs and historical SST data. The low SST in the WSWP generated a positive zonal SST gradient together with high SST east of the WSWP. It may contribute to enhancement of the westerly surface wind in this region, leading to the onset of the 2002/03 El Niño event.     

Bathymetry enhancement by altimetry-derived gravity anomalies in the East Sea (Sea of Japan)

The gravity-geologic method (GGM) was used to enhance the bathymetry of the East Sea (Sea of Japan) with satellite altimetry-derived free-air gravity anomalies and shipborne depth measurements. By comparison with the bathymetry model of Smith and Sandwell’s (SAS) approach (1994), GGM was found to have an advantage with short wavelength (≤12 km) components, while SAS better predicts longer wavelength (≥25 km) components, despite its dependency on density contrast. To mitigate this limitation, a tuning density contrast of 10.25 g/cm³ between seawater and the seafloor was primarily estimated by the downward continuation method and then validated by the check points method with GGM. Similarly, SAS is limited by the “A” value in low-pass part of the Wiener filter, which defines the effective range of the wavelength components on bathymetry. As a final result, we present an enhanced GGM bathymetry model by integrating all available data.     

Trend and dynamic cause of sediment particle size on the adjacent continental shelf of the Yangtze Estuary

Based on the measured data in recent 20 years, the variation trends of the median grain size of the surface sediment, the sand-silt boundary and the mud area on the adjacent continental shelf of the Yangtze Estuary were analyzed in depth, and the effects of natural mechanism and human activities were discussed. The results show that:(1) In recent years(2006-2013), the median grain size of sediment and the distribution pattern of grouped sediments in the adjacent continental shelf area to the Yangtze Estuary have presented no obvious change compared with those before 2006;(2) The median diameter of the surface sediment in the continental shelf area displayed a coarsening trend with the decrease of sediment discharge from the basin and the drop in suspended sediment concentration in the shore area;(3) In 2004-2007, the sand-silt boundary in the north part(31°30′N) of the continental shelf area presented no significant changes, while that in the south part(31°30′S) moved inwards; In 2008-2013, both the sand-silt boundaries in the north and south parts of the continental shelf area moved inwards, mainly due to the fact that in the dry season, a relatively enhanced hydrodynamic force of the tides was generated in the Yangtze River, as well as a decreased suspended sediment concentration and a flow along the banks in North Jiangsu;(4) The mud area where the maximum deposition rate is found in the Yangtze Estuary, tends to shrink due to the drop in sediment discharge from the basin, and the decrease in suspended sediment concentration in the shore area and erosion in the delta. Moreover, it tended to shift to the south at the same time because the implement of the training works on the deep-water channel of the North Passage changed the split ratio between the North and South Passages with an increase in the power of the discharged runoff in the South Passage.     

High-frequency ocean radar derived characteristics of sea surface currents in the Ariake Sea,Japan

Composite analysis was conducted using high-frequency radar data obtained during 2006–2015 in order to gain a better understanding of the current field in the Ariake Sea. The seasonally averaged surface current in the Ariake Sea was directed southward in all seasons, except around river mouths during summer. Heavy rainfall enhanced the outflow along the eastern coast of the Shimabara Peninsula from Isahaya Bay to the southern area 2–5 days after heavy rainfall. Spring–neap differences were clearly seen in the southward current along the Shimabara Peninsula. Interannual variation in the M₂ tidal current amplitude was synchronized with the lunar nodal cycle.     

Free Vibration Analysis of Symmetrically Laminated Composite Plates on Elastic Foundation and Coupled with Stationary Fluid

Free vibration analysis of symmetrically laminated composite plates resting on Pasternak elastic support and coupled with an ideal, incompressible and inviscid fluid is the objective of the present work. The fluid domain is considered to be infinite in the length direction but bounded in the depth and width directions. In order to derive the eigenvalue equation, Rayleigh-Ritz method is applied for the fluid-plate-foundation system. The efficiency of the method is proved by comparison studies with those reported in the open literature. At the end, parametric studies are carried out to examine the impact of different parameters on the natural frequencies.