Reliability of reanalyses products in simulating precipitation and temperature characteristics over India

Various reanalyses have been utilized in numerous climate related researches around the globe, however, there exists considerable biasedness in these products, especially in precipitation and temperature data. The ability of these reanalysis products to simulate the precipitation and temperature patterns is observed to be satisfactory at global scale, while it differs significantly at regional scale, especially over regions of high spatio-temporal heterogeneity such as India. Therefore, it is essential to evaluate the applicability and robustness of reanalyses in climate related research. The annual and seasonal variability in spatio-temporal patterns and trends of precipitation and temperature data, with respect to the IMD gridded data over 34 yrs, are evaluated for six global reanalyses namely, NCEP/NCAR Reanalysis (NCEP R1), NCEP-DOE AMIP-2 Reanalysis (NCEP R2), Climate Forecast System Reanalysis (CFSR), ECMWF Interim Reanalysis (ERA-Interim), Modern Era Retrospective Analysis for Research and Application Land only model (MERRA-Land) and JMA 55-year Reanalysis (JRA-55). The ability of the reanalyses was tested based on several factors such as statistical and categorical indices, spells and trends, for annual and seasonal daily values. Several regional and seasonal differences were observed, particularly over high rainfall regions such as Western Ghats and northeastern India. MERRA-Land is found to give the best results for precipitation over India, which is attributed to the updated forcing data using gauge-based precipitation observations. Similarly, ERA-Interim and JRA-55 exhibit better performance for temperature than other datasets. All reanalyses failed to correctly reproduce the trends in IMD data, for both precipitation and temperature. These observations will provide a better perception on the reliability and applicability of reanalyses for climate and hydrological studies over India.     

Tracing submarine groundwater discharge flux in Tolo Harbor,Hong Kong (China)

Submarine groundwater discharge (SGD) is an important pathway for groundwater and associated chemicals to discharge to the sea. Groundwater levels monitored along a transect perpendicular to the shoreline are used to calculate SGD flux from the nearshore aquifer to Tolo Harbor, Hong Kong (China). The calculated SGD flux—recharge/discharge measured with Darcy’s Law methods—agrees well with estimates based on geo-tracer techniques and seepage meter in Tolo Harbor during previous studies. The estimated freshwater SGD is 1.69–2.0 m²/d at the study site and 0.3?±?0.04 cm/d for the whole of Tolo Harbor, which is comparable to the river discharge (0.25?±?0.07 cm/d) and precipitation (0.45?±?0.15 cm/d). The tide-driven SGD in the intertidal zone is 13.98–17.59 m²/d at the study site and 2.42?±?0.56 cm/d for the whole of Tolo Harbor. The SGD occurring in the subtidal zone and the bottom of Tolo Harbor is 3.12?±?4.63 cm/d. Fresh SGD accounts for ~5% of the total SGD, while the rest (~95%) is contributed by saline SGD driven by various forces. About 96% of the tide-driven SGD in the intertidal zone occurs in the ebbing tide period because the head difference between the groundwater level and sea level is great during this period. Tide-driven SGD in the spring tide is ~1.2 times that during neap tide. The tidal fluctuation amplitude and tide-driven SGD in the intertidal zone are positively correlated to each other; thus, a spring neap variation of the tide-driven SGD is observed.     

Were merely storm-landslides driven by the 2015-2016 Niño in the Mendoza River valley?

Since Holocene time, above-mean precipitations recorded during the El Niño warm ENSO phase have been linked to the occurrence of severe debris flows in the arid Central Andes. The 2015–2016 El Niño, for its unusual strength, began driving huge and dangerous landslides in the Central Andes (32°) in the recent South Hemisphere summer. The resulting damages negatively impacted the regional economy. Despite this, causes of these dangerous events were ambiguously reported. For this reason, a multidisciplinary study was carried out in the Mendoza River valley. Firstly, a geomorphological analysis of affected basins was conducted, estimating morphometric parameters of recorded events such as velocity, stream flow, and volume. Atmospheric conditions during such events were analyzed, considering precipitations, snow cover, temperature range, and the elevation of the zero isotherm. Based on our findings, the role of El Niño on the slope instability in the Central Andes is more complex in the climate change scenario. Even though some events were effectively triggered by intense summer rainstorm following expectations, the most dangerous events were caused by the progressive uplifting of the zero isotherm in smaller basins where headwaters are occupied by debris rock glaciers. Our research findings give light to the dynamic coupled system ENSO–climate change–landslides (ECCL) at least in this particular case study of the Mendoza River valley. Landslide activity in this Andean region is driven by wetter conditions linked to the ENSO warm phase, but also to progressive warming since the twentieth century in the region. This fact emphasizes the future impact of the natural hazards on Andean mountain communities.     

Ensemble clustering for efficient robust optimization of naturally fractured reservoirs

In the development of naturally fractured reservoirs (NFRs), the existence of natural fractures induces severe fingering and breakthrough. To manage the flooding process and improve the ultimate recovery, we propose a numerical workflow to generate optimal production schedules for smart wells, in which the inflow control valve (ICV) settings can be controlled individually. To properly consider the uncertainty introduced by randomly distributed natural fractures, the robust optimization would require a large ensemble size and it would be computationally demanding. In this work, a hierarchical clustering method is proposed to select representative models for the robust optimization in order to avoid redundant simulation runs and improve the efficiency of the robust optimization. By reducing the full ensemble of models into a small subset ensemble, the efficiency of the robust optimization algorithm is significantly improved. The robust optimization is performed using the StoSAG scheme to find the optimal well controls that maximize the net-present-value (NPV) of the NFR’s development. Due to the discrete property of a natural fracture field, traditional feature extraction methods such as model-parameter-based clustering may not be directly applicable. Therefore, two different kinds of clustering-based optimization methods, a state-based (e.g., s _w profiles) clustering and a response-based (e.g., production rates) clustering, are proposed and compared. The computational results show that the robust clustering optimization could increase the computational efficiency significantly without sacrificing much expected NPV of the robust optimization. Moreover, the performance of different clustering algorithms varies widely in correspondence to different selections of clustering features. By properly extracting model features, the clustered subset could adequately represent the uncertainty of the full ensemble.     

A method of alternating characteristics with application to advection-dominated environmental systems

We present a numerical method for solving a class of systems of partial differential equations (PDEs) that arises in modeling environmental processes undergoing advection and biogeochemical reactions. The salient feature of these PDEs is that all partial derivatives appear in linear expressions. As a result, the system can be viewed as a set of ordinary differential equations (ODEs), albeit each one along a different characteristic. The method then consists of alternating between equations and integrating each one step-wise along its own characteristic, thus creating a customized grid on which solutions are computed. Since the solutions of such PDEs are generally smoother along their characteristics, the method offers the potential of using larger time steps while maintaining accuracy and reducing numerical dispersion. The advantages in efficiency and accuracy of the proposed method are demonstrated in two illustrative examples that simulate depth-resolved reactive transport and soil carbon cycling.     

The geological sources of Hg contamination in Kejimkujik National Park,Nova Scotia,Canada: a GIS approach

Recent studies have shown that mercury (Hg) levels in many fish from remote lakes exceed the recommended guidelines for human consumption. Most of these studies conclude that the source of contamination lies in the atmosphere. Kejimkujik National Park (KNP), Nova Scotia, Canada, is considered to be a pristine ecosystem in which fish and loon Hg levels are anomalously high. Studies in the park have shown that atmospheric Hg concentrations may not be high enough to account for the Hg levels in the biota, indicating that the park may be an unusual system in terms of Hg distribution and migration. In an attempt to summarise and synthesise the numerous Hg data sets which have been produced in the park over the last 5-10 years, a geographic information systems (GIS) approach was used to create a common database using the watersheds in the park as the common parameter. By using a GIS database, new relationships and correlations are established and the spatial distribution of Hg levels is more readily evaluated and quantified. The results indicate that geological sources of Hg, arising from biotite-rich granite rocks, may play a larger role in the contamination of the park than previously thought.     

Tsunami hazard in the Eastern Mediterranean: strong earthquakes and tsunamis in Cyprus and the Levantine Sea

A modern tsunami catalogue has been compiled for the region of Cyprus-Levantine Sea in which 24 certain or possible local tsunamis are listed from antiquity up to the present time, while six regional tsunamis, generated in the Hellenic arc, are documented which affected the region. Another set of 13 doubtful events not included in the catalogue are discussed. Tsunami intensities k and K were re-evaluated using the classic 6-grade and the new 12-grade intensity scales, respectively. The strongest tsunamis reported in the region of interest are those of 551 AD, 749, 1068, 1201, 1222, 1546 and 1759, all occurring along the Levantine coast from Gaza northward, with the exception of the 1222 wave which occurred in the Cyprean arc. The causative earthquakes, however, occur on land and are associated with the left-lateral strike-slip Levantine rift and, as such, remain unexplained. In this paper we speculate on the mechanism of these events. A second tsunami zone follows the Cyprean arc, where the situation of subaqueous seismogenic sources favours the generation of tsunamis by co-seismic fault displacements. Submarine or coastal earth slumping, however, may be an additional tsunamigenic component. Based on historical data, the average tsunami recurrence in the Cyprus-Levantine Sea region is roughly estimated to be around 30 years, 120 years and 375 years for moderate (k/K ≥ 2/3), strong (k/K ≥ 3/5) and very strong (k/K ≥ 5/8) events, respectively. The rate of tsunami occurrence equals 0.033, 8.3 × 10⁻³ and 2.7 × 10⁻³ events/year for intensity k/K ≥ 2/3, 3/5 and 5/8, respectively. For a Poissonian (random) process the probabilities of observing at least one moderate, strong or very strong tsunami are 0.28, 0.01 and 3 × 10⁻³ within 1 year, 0.81, 0.34 and 0.13 within 50 years and 0.96, 0.56 and 0.24 within 100 years, respectively. The tsunami potential in the Cyprus-Levantine Sea area is low relative to other Mediterranean tsunamigenic regions. However, the destructiveness of some historical events indicates the need to evaluate tsunami hazard by all available means. In addition, remote tsunamigenic sources, such as those of 1303 and 1481 in the eastern Hellenic arc, are able to threaten the coasts of the Cyprus-Levantine region and, therefore, such regional tsunamis should be taken into account in the evaluation of the tsunami risk of the region.