A review of the NEMURO and NEMURO.FISH models and their application to marine ecosystem investigations

The evolution of the North Pacific Ecosystem Model for Understanding Regional Oceanography (NEMURO) family of models to study marine ecosystems is reviewed. Applications throughout the North Pacific have shown the models to be robust and to be able to reproduce 1D, 2D and 3D components of nutrient, carbon cycle and biogeochemical cycles as well as aspects of the lower trophic levels ecosystem (phyto- and zooplankton). NEMURO For Including Saury and Herring, an extension that includes higher trophic levels, can be run uncoupled or coupled to NEMURO. In the uncoupled mode, the growth and weight of an individual fish is computed using plankton densities simulated by NEMURO but with no feedback between fish consumption and plankton mortality. In the coupled mode, the feeding, growth and weight of a representative fish are computed, and prey removals due to feeding by fish appear as mortality terms on the prey. The NEMURO family of models continues to evolve, including effects of the microbial loop and iron limitation at lower trophic levels, and full life cycle, multi-species and multi-generational simulations at higher trophic levels. We outline perspectives for future end-to-end modeling efforts that can be used to study marine ecosystems in response to global environmental change.     

Assessing desertification by using soil indices

Desertification generally refers to land degradation in arid, semiarid, and dry semi-humid climatic zones. It involves five principal processes: vegetation degradation, water erosion, wind erosion, salinization and waterlogging, and soil crusting and compaction. The aim of this study is assessing desertification using soil criteria. For this purpose, nine indices including sodium absorption ratio (SAR), soil gypsum percentage, soil texture, the content of HCO₃ ^?, the percentage of the organic matter, electrical conductivity (EC), pH, the content of the soil sodium, and chloride were used. The soil samples were taken in the north of Zayandeh-Rood River in Isfahan province of Iran, using soil data randomly sampled in a depth of 0–20 cm. After assessing the normality of the samples using Kolmogorov-Smirnov test, indices were imported into GIS environment and interpolated with IDW and normal and discrete kriging methods for delineating soil characteristics maps based on MEDALUS model. In this model, the data were firstly changed from 100 to 200. Thus 100 and 200 are estimated as the best and worst quality, respectively. Then the final map of soil criteria has been created by geometric mean of its indicators. The results showed that the maximum area is related to the medium class of desertification and is equal to 44,746 ha. The areas of severe and very severe classes of desertification are equal to 30,949 and 351 ha, respectively. The results also revealed that the indices of the organic matter and soil gypsum percentage are the most influential indices which affect desertification phenomenon.     

Tectonophysics and recent geodynamics

Important problems of tectonophysical analysis of observational results in recent geodynamics are considered. The paradox of low rates of recent horizontal deformations of the Earth’s surface is formulated. This paradox states that, according to GPS measurements, the annual mean rates of relative deformations are 10⁻⁹–10⁻¹⁰ per year, whereas, according to ground-based geodetic and deformographic observations, the rates of recent horizontal deformations of the Earth’s surface are on the order of 10⁻⁵–10⁻⁶ per year. This paradox is shown to appear as a result of the application of the procedure of averaging to the displacements obtained on large measuring bases with a low degree of details of the time of observations. It is established that, according to multiyear (over 40 years) data, the annual mean rate of deformation processes in the Kopet Dagh seismically active region is 4 × 10⁻⁸ per annum along the vertical and 3.2 × 10⁻⁸ per year along the horizontal.     

Time and frequency domains response analyses of April 2015 Greece’s earthquake in the Nile Delta based on GNSS-PPP

Earthquake response characteristic and sensitivity of 1-Hz Global Navigation Satellite System - precise point positioning (GNSS-PPP) for a seismic response in time and frequency domains are the main objectives of this study. The Delta and Cairo areas, Egypt, experiencing the shock of the Greece earthquake occurring on April 16, 2015 (18:50 hours), is presented in this study. The results of this study reveal that the seismic wave effect is very small and can be neglected inside the Delta region and the strong motion occurred in an upward direction for the whole monitoring area. Furthermore, the time-frequency analysis illustrates the ability of the time and frequency domains to use and analyze the surface motion based on a 1-Hz sampling frequency and to detect the small- and high-strength motions of seismic waves.     

Comparing data-driven landslide susceptibility models based on participatory landslide inventory mapping in Purwosari area,Yogyakarta, Java

There are different approaches and techniques for landslide susceptibility mapping. However, no agreement has been reached in both the procedure and the use of specific controlling factors employed in the landslide susceptibility mapping. Each model has its own assumption, and the result may differ from place to place. Different landslide controlling factors and the completeness of landslide inventory may also affect the different result. Incomplete landslide inventory may produce significance error in the interpretation of the relationship between landslide and controlling factor. Comparing landslide susceptibility models using complete inventory is essential in order to identify the most realistic landslide susceptibility approach applied typically in the tropical region Indonesia. Purwosari area, Java, which has total 182 landslides occurred from 1979 to 2011, was selected as study area to evaluate three data-driven landslide susceptibility models, i.e., weight of evidence, logistic regression, and artificial neural network. Landslide in the study area is usually affected by rainfall and anthropogenic activities. The landslide typology consists of shallow translational and rotational slide. The elevation, slope, aspect, plan curvature, profile curvature, stream power index, topographic wetness index, distance to river, land use, and distance to road were selected as landslide controlling factors for the analysis. Considering the accuracy and the precision evaluations, the weight of evidence represents considerably the most realistic prediction capacities (79%) when comparing with the logistic regression (72%) and artificial neural network (71%). The linear model shows more powerful result than the nonlinear models because it fits to the area where complete landslide inventory is available, the landscape is not varied, and the occurence of landslide is evenly distributed to the class of controlling factor.     

A fine-granularity scheduling algorithm for parallel XDraw viewshed analysis

Viewshed analysis is widely used in many terrain applications such as siting problem, path planning problem, and etc. But viewshed computation is very time-consuming, in particular for applications with large-scale terrain data. Parallel computing as a mainstream technique with the tremendous potential has been introduced to enhance the computation performance of viewshed analysis. This paper presents a revised parallel viewshed computation approach based on the existing serial XDraw algorithm in a distributed parallel computing environment. A layered data-dependent model for processing data dependency in the XDraw algorithm is built to explore scheduling strategy so that a fine-granularity scheduling strategy on the process-level and thread-level parallel computing model can be accepted to improve the efficiency of the viewshed computation. And a parallel computing algorithm, XDraw-L, is designed and implemented taken into account this scheduling strategy. The experimental results demonstrate a distinct improvement of computation performance of the XDraw-L algorithm in this paper compared with the coarse-partition algorithm, like XDraw-E which is presented by Song et al. (Earth Sci Inf 10(5):511–523, 2016), and XDraw-B that is the basic algorithm of serial XDraw. Our fine-granularity scheduling algorithm can greatly improve the scheduling performance of the grid cells between the layers within a triangle region.     

Stability study on the northern batter of MBC Open Pit using Plaxis 3D

Cracks appeared on the northern batter at Maddingley Brown Coal Open Pit Mine, Victoria, Australia, on 8 November 2013 and a 2-day rainfall event happened 5 days later. This study models the stability of the northern batter considering the effect of the rainfall event and an emergency buttress using finite element method (FEM) encoded in Plaxis 3D. It is found that the batter tended to lead to block sliding after overburden removal. The observed vertical crack would be a combined action of the overburden removal and groundwater flow. The simulated location of cracks agrees well with the actual location, and the simulated heave of the coal seam is in good agreement with the experience in Victoria brown coal open pit mining. The rainfall accelerated the development of the cracks. With the construction of the emergency buttress, the batter became stable that is in good agreement with the monitored data.     

Analysis of radial movement of an unconfined leaky aquifer due to well pumping and injection

Radial movement of an unconfined leaky aquifer was studied with respect to hydraulic forces that are induced by well recharge and discharge. New analytic solutions in the velocity and displacement fields were found and applied to describe transient movement in an unconfined leaky aquifer. Linear momentum and mass balance of saturated porous sediments, the Darcy-Gersevanov law, and the analytic solution of hydraulic drawdown for unsteady flow within the unconfined leaky aquifer were introduced to find the new solutions. Analytic results indicate that the nonlinear relation between the initial hydraulic head (h₀) and the well function has an insignificant effect on the aquifer transient movement when the drawdown s<0.02h ₀. When the well function is simplified with different assumptions and pumping conditions, the new solutions correspondingly reduce to cases that are similar to the Hantush-Jacob, Muskat, and Theis transient movement of a confined leaky aquifer. It was found that large leakance is important in slowing radial movement and reducing aquifer deformation. Flow velocity in the aquifer is more responsive to leakance than to cumulative displacement within the aquifer. The zones and boundary with tensile stress can be located using the same approach applied to a confined aquifer for risk assessment of earth fissuring.

---

Résumé Le mouvement radial dans un aquifère semi-perméable libre a été étudié en considérant les forces hydrauliques induites par la recharge et le prélèvement au niveau d’un puits. De nouvelles solutions analytiques des champs de vitesse et de déplacement ont été trouvées et appliquées pour décrire le mouvement transitoire dans un aquifère semi-perméable libre. Pour trouver les nouvelles solutions, le moment linéaire et le bilan de masse dans les sédiments poreux saturés, la loi de Darcy-Gersevanov, et la solution analytique du rabattement hydraulique pour un écoulement transitoire dans un aquifère semi-perméable libre sont pris en compte. Les résultats analytiques montrent que la relation non-linéaire entre le potentiel hydraulique initial (h₀) et la fonction de puits a un effet insignifiant sur le mouvement transitoire de l’aquifère quand le rabattements s <0.02h₀. Si la fonction de puits est simplifiée en tenant compte de différentes hypothèses et conditions de pompage, les nouvelles solutions se réduisent en conséquence à des cas similaires au mouvement transitoire dans un aquifère semi-perméable captif selon Hantush-Jacob, Muskat, et Theis. Il s’est avéré qu’une drainance élevée est importante pour le ralentissement du mouvement radial et la réduction de la déformation de l’aquifère. La vitesse d’écoulement dans l’aquifère est plus sensible à la drainance qu’au déplacement cumulé dans l’aquifère. Les zones et la frontière présentant des contraintes de traction peuvent être localisées en utilisant la même approche appliquée à un aquifère captif pour l’évaluation des risques de fissuration du sol.

---

Resumen Se ha estudiado el movimiento radial en un acuífero no confinado con fuga en relación a las fuerzas hidráulicas que son inducidas por la descarga y recarga del pozo. Se encontraron nuevas soluciones analíticas en los campos de desplazamiento y velocidad, las cuales se aplicaron para describir el movimiento transitorio en un acuífero no confinado con fuga. Se han introducido cálculos de balance de masa y momentum lineal de los sedimentos porosos saturados, la Ley Darcy-Gersevanov, y la solución analítica de descenso hidráulico para flujo sin régimen permanente dentro del acuífero no confinado con goteo para encontrar las nuevas soluciones. Los resultados analíticos indican que la relación no lineal entre la presión hidráulica inicial (h₀) y la función del pozo tiene efecto insignificante en el movimiento transitorio del acuífero cuando el descenso s < 0.02h₀. Cuando se simplifica la función del pozo en base a diferentes supuestos y condiciones de bombeo se reducen las nuevas soluciones a casos que son similares a el movimiento transitorio Hantush-Jacob, Muskat, y Theis de un acuífero confinado con goteo. Se encontró que las fugas grandes son importantes en la disminución del movimiento radial y en la reducción de la deformación del acuífero. La velocidad de flujo en el acuífero responde más a las fugas que al desplazamiento acumulativo dentro del acuífero. Pueden localizarse zonas y límites con esfuerzos tensionales usando el mismo enfoque aplicado a acuíferos confinados para evaluación de riesgo de generación de fracturas en el terreno.

---

     

Spatial and Temporal Variations in <Emphasis Type="Italic">Pyrodinium bahamense</Emphasis> Cyst Concentrations in the Sediments of Bioluminescent Mangrove Lagoon,St. Croix,USVI

Mangrove Lagoon, located on the island of St. Croix, US Virgin Islands (USVI), is one of few actively bioluminescent lagoons in a location experiencing significant anthropogenic impacts. The bioluminescence is due to an abundance of the dinoflagellate Pyrodinium bahamense in the water column. We recovered surface sediments and sediment cores from Mangrove Lagoon to analyze the spatial distribution and temporal variability of P. bahamense cysts in this system. Surface sediment P. bahamense cyst concentrations ranged from 0 to 466 cysts g^?1 dry sediment, with higher abundances associated with elevated surface water nutrient concentrations and a mixed terrestrial–marine organic matter source regime. In combination with available bioassay data, we hypothesize that phytoplankton utilize nutrients rapidly and subsequent decay of organic matter makes nutrients available for dinoflagellates at the sediment–water interface in the eastern and northern quadrants of the lagoon. However, the nutrients are rapidly exhausted during counterclockwise lagoon circulation resulting in the decline of primary productivity and dinoflagellate abundance in the western quadrants. Downcore profiles suggest that P. bahamense blooms have been occurring for decades, declining in recent years. No cysts were present in sediments predating dredging activities of the 1960s that created Mangrove Lagoon. Recent reductions in cyst abundance may be the result of limited primary productivity caused by restricted water exchange with Salt River Bay due to shallowing of a sill at the mouth of the lagoon. This research highlights the need for more comprehensive geochemical and fossil analyses to better understand long-term ecological variability and inform conservation efforts of these unique habitats.