A review of open software resources in python for electrical resistivity modelling

Geophysical modelling performs to obtain subsurface structures in agreement with measured data. Freeware algorithms for geoelectrical data inversion have not been widely used in geophysical communities; however, different open-source modelling/inversion algorithms were developed in recent years. In this study, we review the structures and applications of openly Python-based inversion packages, such as pyGIMLi (Python Library for Inversion and Modelling in Geophysics), BERT (Boundless Electrical Resistivity Tomography), ResIPy (Resistivity and Induced Polarization with Python), pyres (Python wrapper for electrical resistivity modelling), and SimPEG (Simulation and Parameter Estimation in Geophysics). In addition, we examine the recovering ability of pyGIMLi, BERT, ResIPy, and SimPEG freeware through inversion of the same synthetic model forward responses. A versatile pyGIMLi freeware is highly suitable for various geophysical data inversion. The SimPEG framework is developed to allow the user to explore, experiment with, and iterate over multiple approaches to the inverse problem. In contrast, BERT, pyres, and ResIPy are exclusively designed for geoelectric data inversion. BERT and pyGIMLi codes can be easily modified for the intended applications. Both pyres and ResIPy use the same mesh designs and inversion algorithms, but pyres uses scripting language, while ResIPy uses a graphical user interface (GUI) that removes the need for text inputs. Our numerical modelling shows that all the tested inversion freeware could be effective for relatively larger targets. pyGIMLi and BERT could also obtain reasonable model resolutions and anomaly accuracies for small-sized subsurface structures. Based on the heterogeneous layered model and experimental target scenario results, the geoelectrical data inversion could be more effective in pyGIMLi, BERT, and SimPEG freeware packages. Moreover, this study can provide insight into implementing suitable inversion freeware for reproducible geophysical research, mainly for geoelectrical modelling.

The dependence of upper ocean gyres on wind and buoyancy forcing

A series of numerical simulations with different forcing conditions are carried out, to investigate the roles played by buoyancy and wind forcing on the upper ocean gyres, and to contrast the laminar and eddying regimes. Model experiments show that the buoyancy-driven eastward geostrophic flow tends to suppress the formation of the wind-driven subpolar gyre, but the northward eddy heat transport can homogenize the subpolar water and reduce the meridional temperature gradient by about two-third, thus counteracting the buoyancy effect and saving the subpolar gyre. For the subtropical gyre, its transport is enhanced by eddy mixing, and the role of buoyancy forcing is very sensitive to the choice of diapycnal diffusivity. Our results suggest that eddy effects must be considered in the dynamics of the subpolar gyre, and vertical diffusivity should be selected carefully in simulating the basin-wide circulations.

Upper crustal shear-wave velocity structure Beneath Western Java, Indonesia from seismic ambient noise tomography

This paper presents the depth inversion of Rayleigh wave group velocity to obtain an S-wave velocity model from seismic ambient noise cross-correlation in western Java, Indonesia. This study utilizes the vertical component data of a temporary seismograph network deployed in 2016, which was used in a previous study to estimate fundamental mode Rayleigh wave group velocity maps. In this study, the Neighborhood Algorithm was applied to invert the Rayleigh wave group velocities into 1D shear-wave velocity (Vs) profiles, which were then interpolated to produce a high-resolution, pseudo-3D Vs model. These tomographic images of Vs extend to ~ 20 km depth and show a pronounced NE-SW contrast of low and high Vs in the depth range 1–5 km that correlates well with the Bouguer anomaly map. We interpret the low Vs in the northeastern part of the study area as associated with alluvial and volcanic products from the Sunda Shelf and modern volcanic arc, whereas the high Vs in the southwestern part is associated with volcanic arc products from earlier episodes of subduction. We also obtained the depth of the northern Java Basin, which is in the range of 5–6 km, and the Garut Basin, which extends to 5 km depth. For greater depths, Vs gradually increases throughout western Java, which reflects the crystalline basement. This study provides estimates of the shallow crustal Vs structure underneath West Java with higher resolution than previous tomographic studies, which could be useful for supporting future earthquake studies in the region.

Extreme IOD induced tropical Indian Ocean warming in 2020

The tropical Indian Ocean (TIO) basin-wide warming occurred in 2020, following an extreme positive Indian Ocean Dipole (IOD) event instead of an El Niño event, which is the first record since the 1960s. The extreme 2019 IOD induced the oceanic downwelling Rossby waves and thermocline warming in the southwest TIO, leading to sea surface warming via thermocline-SST feedback during late 2019 to early 2020. The southwest TIO warming triggered equatorially antisymmetric SST, precipitation, and surface wind patterns from spring to early summer. Subsequently, the cross-equatorial “C-shaped” wind anomaly, with northeasterly–northwesterly wind anomaly north–south of the equator, led to basin-wide warming through wind-evaporation-SST feedback in summer. This study reveals the important role of air–sea coupling processes associated with the independent and extreme IOD in the TIO basin-warming mode, which allows us to rethink the dynamic connections between the Indo-Pacific climate modes.

Earthquake monitoring of the Baribis Fault near Jakarta, Indonesia, using borehole seismometers

The geological setting of Jakarta and its immediate surroundings are poorly understood, yet it is one of the few places in Indonesia that is impacted by earthquakes from both the Java subduction zone and active faults on land. In this study, a borehole seismic experiment with low noise characteristics was deployed to record seismic activity on the ~ E-W oriented Baribis Fault, which is ~ 130 km long, passes to the south of Jakarta, and is only ~ 20 km away at its nearest point. A primary objective of this study is to determine whether this fault is seismically active, and therefore, whether it might pose a threat to nearby population centers, including Jakarta in particular. A total of seven broadband instruments that spanned Jakarta and the surrounding region were installed between the end of July 2019 and August 2020, during which time we were able to detect and locate 91 earthquakes. Two earthquakes were located close to the Baribis Fault line, one of which was felt in Bekasi (southeast of Jakarta) where it registered II-III on the Modified Mercalli Intensity (MMI) scale. The focal mechanism solutions of these events indicate the presence of a thrust fault, which is in good agreement with previous studies, and suggest that the Baribis Fault is active.

Does the wind stress always damp an oceanic eddy?

Up to now, the literature has shown that the relative wind stress does negative work on ocean mesoscale eddies. In other words, the relative wind stress inhibits the development of the eddies. However, based on a newly derived simplified theoretical model, the present study finds that under the action of a steady and uniform wind field, eddies can rapidly obtain kinetic energy from the wind field following several hours of adaption and adjustment, in which the wind stress transitions from doing negative to positive work. The finding is supported by the fact that the relative wind stress work on oceanic eddies over the northeastern tropical Pacific ocean is positive with the nearly constant gap wind. This implies that energy input from the wind is sensitive to eddy velocity structure, and hence, wind stress is not always a killer of eddies.

Role of the eastern boundary-generated waves on the termination of 1997 Indian Ocean Dipole event

The termination of Indian Ocean Dipole (IOD) events is examined in terms of equatorial wave dynamics. In situ and satellite observations combined with an output from a linear wave model are used in this study. Our emphasis is on the 1997 IOD event but our results apply to other positive IOD events as well. We find that the termination of anomalously cold sea surface temperature (SST) in the eastern pole of the dipole is associated with a warming tendency caused by the net surface heat fluxes. However, net surface heat fluxes alone cannot explain the total change in the SST. We show that during the peak phase of an IOD event, the weakening of zonal heat advection caused by eastern boundary-generated Rossby waves combined with the reduction of vertical entrainment and diffusion creates favorable conditions for surface heat fluxes to warm the SST in the eastern basin.

Impact of land-use change on catchment water balance: a case study in the central region of South Africa

An assessment study on the impact of land use change through afforestation on catchment water balance was carried out in one of the semi-arid quaternary catchments (C52A) of the Modder River Basin located in the central region of South Africa. The study used ArcGIS and the Soil and Water Assessment Tool (SWAT) to create scenarios of forest cover areas and to simulate the water balance of C52A, respectively. Climatic data from 1993 to 2020 were used to simulate the catchment water balance. The climatic data from 2011 to 2020 were generated using SWAT weather generator while the rest was obtained from South Africa Weather Service at three weather stations located within C52A catchment. In the C52A quaternary catchment it was envisaged to increase forest cover by considering afforestation on pastureland that is found on slopes > 8% and > 3% which created two forest scenarios. The baseline land use [Land use 2000 (LU2000)] which is taken as the base scenario in this study has about 84% of the catchment area covered with pasture. The two forest scenarios created were: forest scenario1 (FRSE1), i.e. conversion of pastureland on slope > 8% to forest, and forest scenario2 (FRSE2), i.e. conversion of pastureland on slope > 3% to forest. The type of forest considered in these scenarios is an evergreen tree, acronym as FRSE by SWAT land use classification. The conversion increased forest cover by 8.3% and 30.5% on FRSE1 and FRSE2, respectively. The result of the water balance of the catchment based on the land use scenarios were compared with the baseline land use scenario (LU2000). The result obtained showed that FRSE1 produced a non-significant change both on the mean monthly surface runoff and water yield compared to LU2000. On the other hand, FRSE2 showed 30% decrease on the mean monthly surface runoff, but increased the mean monthly lateral flow and base flow by 110% and 254%, respectively compared to the LU2000. Thus, in the overall water balance, the mean monthly water yield of the catchment increased by 171% on FRSE2 compared to the LU2000. Although there are considerable number of research reports on the negative effect of forest on catchment water yield, this study showed a significant water yield increase when approximately 30% of the C52A catchment area, which lay on slopes > 3% and covered by grass, was converted to evergreen forest land. The result showed that having a forest cover on a suitable slope range can bring about a positive effect on the total water yield of a catchment. Therefore, this finding is important for catchment management stakeholders and policy-makers when devising land use and water resources management strategies in a catchment.

Local seismotectonic analysis of the July 2019 Molucca Sea earthquake sequence based on moment tensor solutions

We investigate the local seismotectonic of the Molucca Sea area using moment tensor calculations for the earthquakes that occurred in July 2019 at a depth of 10–55 km. The mainshock of Mw 6.8 occurred on July 7, followed by aftershocks until July 18, with magnitudes ranging from Mw 4.6 to Mw 5.8. Moment tensor solutions are calculated by applying Isolated Asperities (ISOLA) software using the full waveform data recorded at regional seismic stations. The analyzed frequency bands used in this study are 0.01–0.03 Hz and 0.04–0.05 Hz for the event with Mw ≥ 5 and Mw < 5, respectively. We provide validations of new moment tensor solutions for Mw < 5 events in the Molucca Sea region for the period during the earthquake sequence. The results show that thrust and oblique faults are dominant during this event, which indicate a compressional stress of divergent double subduction (DDS) of the Sangihe and Halmahera arcs. Only one full moment tensor solution reveals the normal fault mechanism, which may indicate the manifestation of strain release of compressional stress in the surrounding area. Furthermore, these results also support the previous studies suggesting that the Talaud-Mayu Ridge located in the middle of the Molucca Sea has developed as a consequence of the transpressional tectonic activity.

Are only floods with large discharges threatening? Flood characteristics evolution in the Yangtze River Basin

This study focuses on the evolution of flood risk in the Yangtze River Basin under climate change, which is a critical issue for socioeconomic development in future. In this study, we (1) compared the 1998 and 2020 floods and found that the destructiveness of a given discharge is now greater than before; (2) revealed three issues related to the above finding; and (3) prospected the future development of up-to-date technologies to better address the issue that floods with high water levels will frequently threaten us. The outcomes of this study would be of great significance to future flood control operation of large river basins.