Spatially distributed tracer‐aided modelling to explore water and isotope transport,storage and mixing in a pristine,humid tropical catchment

Rapidly transforming headwater catchments in the humid tropics provide important resources for drinking water, irrigation, hydropower, and ecosystem connectivity. However, such resources for downstream use remain unstudied. To improve understanding of the behaviour and influence of pristine rainforests on water and tracer fluxes, we adapted the relatively parsimonious, spatially distributed tracer‐aided rainfall–runoff (STARR) model using event‐based stable isotope data for the 3.2‐km² San Lorencito catchment in Costa Rica. STARR was used to simulate rainforest interception of water and stable isotopes, which showed a significant isotopic enrichment in throughfall compared with gross rainfall. Acceptable concurrent simulations of discharge (Kling–Gupta efficiency [KGE] ~0.8) and stable isotopes in stream water (KGE ~0.6) at high spatial (10 m) and temporal (hourly) resolution indicated a rapidly responding system. Around 90% of average annual streamflow (2,099 mm) was composed of quick, near‐surface runoff components, whereas only ~10% originated from groundwater in deeper layers. Simulated actual evapotranspiration (ET) from interception and soil storage were low (~420 mm/year) due to high relative humidity (average 96%) and cloud cover limiting radiation inputs. Modelling suggested a highly variable groundwater storage (~10 to 500 mm) in this steep, fractured volcanic catchment that sustains dry season baseflows. This groundwater is concentrated in riparian areas as an alluvial–colluvial aquifer connected to the stream. This was supported by rainfall–runoff isotope simulations, showing a “flashy” stream response to rainfall with only a moderate damping effect and a constant isotope signature from deeper groundwater (~400‐mm additional mixing volume) during baseflow. The work serves as a first attempt to apply a spatially distributed tracer‐aided model to a tropical rainforest environment exploring the hydrological functioning of a steep, fractured‐volcanic catchment. We also highlight limitations and propose a roadmap for future data collection and spatially distributed tracer‐aided model development in tropical headwater catchments.     

Late Quaternary (Marine Isotope Stage 3 to Recent) sedimentary evolution of the Basque shelf (southern Bay of Biscay)

Late Quaternary (MIS 3 to Recent) oceanographic evolution of the Basque shelf has been analysed for the first time based on the sedimentological analysis of three cores obtained from the middle and outer shelves. The cores are located in two interfluves separated by the San Sebastian canyon. The variability of the percentage of the planktonic foraminifera species Neogloboquadrina pachyderma sin. and of δ¹⁸O_bull allowed us to identify the influence of colder and warmer waters in the Basque shelf during the late Quaternary. From ～56 cal. ka BP to the end of the Last Glacial Maximum (19 cal. ka BP) the sedimentary record shows a decreasing trend in the mean grain size that correlates with the eustatic sea‐level fall. The last Deglaciation (19–11.5 cal. ka BP) is characterized by a sea‐level rise that produced an important hiatus in the western outer shelf. During the Holocene, the middle and outer shelves present different behaviours. From 11.5 to 6.7 cal. ka BP, in the outer shelf the sea‐level rise that started during the Deglaciation produced a hiatus, whereas in the middle shelf the sedimentary succession records the presence of warm to temperate waters. Between 6.7–4.9 cal. ka BP, the entrance of cold surface water‐masses that only affected the middle shelf has been identified, and temperate to warm waters occurred in the outer shelf. The cold surface water‐masses retreated during 4.9–4.3 cal. ka BP in the middle shelf. Finally, from 4.3 cal. ka BP to Recent, the middle shelf registers a hiatus due to sea‐level stabilization after a generalized transgression, synchronous to a decrease in the energy of the water‐masses in the outer shelf. In conclusion, the environmental changes detected in the Basque shelf are attributed to both regional and eustatic sea‐level changes.     

Holocene paleoclimate inferred from stable isotope (δ18O and δ13C) values in Sphagnum cellulose,Mohos peat bog,Romania

Journal of Paleolimnology - We measured stable isotopes (δ18O and δ13C) in Sphagnum cellulose that was extracted from a long peat core drilled in the ombrotrophic Mohos peat bog, Ciomadul...     

Ecological distribution of Stephanodiscus niagarae Ehrenberg in central Mexico and niche modeling for its last glacial maximum habitat suitability in the Nearctic realm

Journal of Paleolimnology - Stephanodiscus niagarae Ehrenberg is currently restricted to specific regions of central Mexico, however, during the late Pleistocene, it had a wider distribution in the...     

Landscape transform and spatial metrics for mapping spatiotemporal land cover dynamics using Earth Observation data-sets

Analysis of Earth observation (EO) data, often combined with geographical information systems (GIS), allows monitoring of land cover dynamics over different ecosystems, including protected or conservation sites. The aim of this study is to use contemporary technologies such as EO and GIS in synergy with fragmentation analysis, to quantify the changes in the landscape of the Rajaji National Park (RNP) during the period of 19 years (1990–2009). Several statistics such as principal component analysis (PCA) and spatial metrics are used to understand the results. PCA analysis has produced two principal components (PC) and explained 84.1% of the total variance, first component (PC1) accounted for the 57.8% of the total variance while the second component (PC2) has accounted for the 26.3% of the total variance calculated from the core area metrics, distance metrics and shape metrics. Our results suggested that notable changes happened in the RNP landscape, evidencing the requirement of taking appropriate measures to conserve this natural ecosystem.     

Influence of subdaily polar motion model on nutation offsets estimated by very long baseline interferometry

This paper studies the connection between the subdaily model for polar motion used in the processing of very long baseline interferometry (VLBI) observations and the estimated nutation offsets. By convention accepted by the International Earth Rotation Service, the subdaily model for polar motion recommended for routine processing of geodetic observations does not contain any daily retrograde terms due to their one-to-one correlation with the nutation. Nevertheless, for a 24-h VLBI solution a part of the signal contained in the polar motion given by the used subdaily model is numerically mistaken for a retrograde daily sidereal signal. This fictitious retrograde daily signal contributes to the estimated nutation, leading to systematic differences between the nutation offsets from VLBI solutions computed with different subdaily polar motion models. We demonstrate this effect using solutions for all suitable 24-h VLBI sessions over a time span of 11 years (2000–2011). By changing the amplitudes of one tidal term in the underlying subdaily model for polar motion and comparing the estimated parameters to the solutions computed with the unchanged subdaily model, the paper shows and explains theoretically the effects produced by the individual subdaily terms on the VLBI nutation estimates.     

On the consistency of the current conventional EOP series and the celestial and terrestrial reference frames

Precise transformation between the celestial reference frames (CRF) and terrestrial reference frames (TRF) is needed for many purposes in Earth and space sciences. According to the Global Geodetic Observing System (GGOS) recommendations, the accuracy of positions and stability of reference frames should reach 1 mm and 0.1 mm year\(^{-1}\), and thus, the Earth Orientation Parameters (EOP) should be estimated with similar accuracy. Different realizations of TRFs, based on the combination of solutions from four different space geodetic techniques, and CRFs, based on a single technique only (VLBI, Very Long Baseline Interferometry), might cause a slow degradation of the consistency among EOP, CRFs, and TRFs (e.g., because of differences in geometry, orientation and scale) and a misalignment of the current conventional EOP series, IERS 08 C04. We empirically assess the consistency among the conventional reference frames and EOP by analyzing the record of VLBI sessions since 1990 with varied settings to reflect the impact of changing frames or other processing strategies on the EOP estimates. Our tests show that the EOP estimates are insensitive to CRF changes, but sensitive to TRF variations and unmodeled geophysical signals at the GGOS level. The differences between the conventional IERS 08 C04 and other EOP series computed with distinct TRF settings exhibit biases and even non-negligible trends in the cases where no differential rotations should appear, e.g., a drift of about 20 \(\upmu \)as year\(^{-1 }\)in \(y_{\mathrm{pol }}\) when the VLBI-only frame VTRF2008 is used. Likewise, different strategies on station position modeling originate scatters larger than 150 \(\upmu \)as in the terrestrial pole coordinates.     

Spherical integral transforms of second-order gravitational tensor components onto third-order gravitational tensor components

New spherical integral formulas between components of the second- and third-order gravitational tensors are formulated in this article. First, we review the nomenclature and basic properties of the second- and third-order gravitational tensors. Initial points of mathematical derivations, i.e., the second- and third-order differential operators defined in the spherical local North-oriented reference frame and the analytical solutions of the gradiometric boundary-value problem, are also summarized. Secondly, we apply the third-order differential operators to the analytical solutions of the gradiometric boundary-value problem which gives 30 new integral formulas transforming (1) vertical-vertical, (2) vertical-horizontal and (3) horizontal-horizontal second-order gravitational tensor components onto their third-order counterparts. Using spherical polar coordinates related sub-integral kernels can efficiently be decomposed into azimuthal and isotropic parts. Both spectral and closed forms of the isotropic kernels are provided and their limits are investigated. Thirdly, numerical experiments are performed to test the consistency of the new integral transforms and to investigate properties of the sub-integral kernels. The new mathematical apparatus is valid for any harmonic potential field and may be exploited, e.g., when gravitational/magnetic second- and third-order tensor components become available in the future. The new integral formulas also extend the well-known Meissl diagram and enrich the theoretical apparatus of geodesy.     

A spatial error model with continuous random effects and an application to growth convergence

We propose a spatial error model with continuous random effects based on Matérn covariance functions and apply this model for the analysis of income convergence processes (\(\beta \)-convergence). The use of a model with continuous random effects permits a clearer visualization and interpretation of the spatial dependency patterns, avoids the problems of defining neighborhoods in spatial econometrics models, and allows projecting the spatial effects for every possible location in the continuous space, circumventing the existing aggregations in discrete lattice representations. We apply this model approach to analyze the economic growth of Brazilian municipalities between 1991 and 2010 using unconditional and conditional formulations and a spatiotemporal model of convergence. The results indicate that the estimated spatial random effects are consistent with the existence of income convergence clubs for Brazilian municipalities in this period.