The setting itself is remarkable. The reserve spans 737,088 hectares of federally protected forest bounded on the east by the Iriri River, a clear-water tributary of the Xingu. No roads or trails reach the study area; access requires roughly two days of navigation along the Xingu and Iriri rivers from Altamira. This remoteness has long kept botanists away, even though recent remote-sensing analyses based on airborne LiDAR data and deep learning had flagged the region as harboring extensive areas of high canopy palm cover. The new field survey confirms that prediction on the ground, and it does so with numbers that would impress any tropical ecologist.

Walking the transect on foot at roughly one kilometer per hour, at least two experienced observers censused a continuous twenty-meter-wide belt, partitioned into one hundred contiguous sampling units of twenty by fifty meters. In total they recorded 1,506 adult palm individuals belonging to eighteen taxa, seventeen identified to species level and one morphospecies to genus. Four species dominated overwhelmingly: Attalea speciosa with 527 individuals across 85 plots, Euterpe longibracteata with 267, Geonoma baculifera with 263, and Euterpe oleracea with 105. Together these four accounted for more than three-quarters of every palm counted. Species accumulation curves and multiple richness estimators, including Chao, Jackknife, and Bootstrap, converged on eighteen species, indicating that the sampling effort captured nearly the entire local palm flora.

Attalea speciosa, the babassu palm, reached an average density of 52.7 individuals per hectare and local peaks of 270 adult individuals per hectare within a single plot. Its hyperdominance is not merely a botanical curiosity. Palms can constitute up to a quarter of stems exceeding ten centimeters in diameter in some tropical forests, and biomass models developed specifically for palms suggest their carbon fixation capacity may be nearly a third higher than estimates based on dicotyledonous trees. In a region flagged by satellite for exceptional palm abundance, documenting which species dominate and where is directly relevant to understanding forest carbon dynamics in the eastern Amazon.

To understand why palms clustered where they did, the researchers extracted environmental data for the centroid of each sampling unit. From the TOPODATA digital elevation model, derived from Shuttle Radar Topography Mission data at thirty-meter resolution, they obtained altitude, slope, relief, exposure, and terrain shape. From HydroSHEDS they calculated the distance to the nearest perennial stream. And from CBERS-4A satellite imagery they derived the normalized difference vegetation index, or NDVI, a proxy for vegetation greenness and forest structure that ranged from 0.597 to 0.662 along the transect. Crucially, the team also recorded in the field the presence or absence of igarapés, the shallow, narrow, slow-moving streams typical of Amazonian environments that are generally invisible in satellite imagery.

The statistical architecture of the study was rigorous. A Redundancy Analysis linked community composition to the continuous environmental variables, with Variance Inflation Factors confirming no problematic multicollinearity. Species richness was modeled with a Gaussian General Linear Model for continuous predictors and analysis of variance for categorical ones, followed by Tukey’s Honestly Significant Difference tests for pairwise comparisons. For the six most abundant species, Bayesian regression models accommodated different response distributions, including zero-inflated count and presence-absence data, quantifying species-specific habitat associations with credible intervals rather than single point estimates.

The results revealed a striking division of labor among environmental factors. Community composition was significantly associated with NDVI, altitude, and distance to streams, which together explained 17.82 percent of the observed variation, with the first constrained axis accounting for 81.84 percent of the constrained variation. Species richness, by contrast, showed no significant association with any continuous variable but responded strongly to the presence of igarapés and to terrain shape. Plots containing igarapés had significantly higher richness than plots without them, with estimated marginal means of 2.07 versus 1.56 species per plot. Terrain shape mattered too: convex, water-divergent areas and concave, water-neutral areas supported higher richness than flat, water-convergent and flat, water-divergent terrains, and the categorical model explained fully half of the observed variation in richness.

The Bayesian models then exposed the ecological personalities of the dominant species. Attalea speciosa showed the strongest associations of all: its expected abundance decreased with increasing NDVI and was lower in plots containing igarapés, indicating a clear preference for relatively open forest patches away from streams. Geonoma baculifera and Euterpe oleracea displayed the opposite pattern, both positively associated with the presence of igarapés and substantially more abundant in riparian plots. Geonoma baculifera formed dense clusters reaching up to 630 individuals per hectare, comparable to densities reported for seasonally flooded sandy soils elsewhere in the eastern Amazon. Euterpe longibracteata responded more to terrain shape than to greenness, favoring flat water-neutral or convex water-convergent topographic positions. Models for the two Astrocaryum species showed lower explanatory power and broad credible intervals, warranting cautious interpretation.

The hyperdominance of Attalea speciosa carries a deeper historical resonance. In the neighboring Rio Iriri reserve on the opposite bank of the river, the same species was found hyperdominant and associated with historical forest disturbance, where archaeological evidence indicates that this disturbance dates to Indigenous occupation during the eighteenth century. The present-day dominance of babassu in the region may therefore reflect a long-term anthropogenic legacy rather than recent extractive activity, although in the newly surveyed area the relative contributions of natural disturbance regimes and historical human influence remain uncertain, and both may have contributed. This possibility ties the palm community to one of the most compelling narratives in Amazonian science: the recognition that many seemingly pristine forests bear the fingerprints of ancient human management.

The richness documented at RESEX Riozinho, seventeen identified species plus one morphospecies, exceeded that of the two other areas previously surveyed along the eastern margin of the Xingu Basin, where ten species were recorded in the Kayapó Indigenous Territory roughly 400 kilometers to the southwest and eleven near Altamira. This pattern fits the broader westward increase in palm diversity across Amazonia, associated with rising rainfall toward the center and west of the basin. Two widespread palms of the Xingu Basin, Socratea exorrhiza, the famous walking palm adapted to seasonally flooded soils through its stilt roots, and Oenocarpus bataua, long associated with Indigenous management and dispersal, both occurred at comparatively low local abundance, showing that regional hyperdominance is not necessarily mirrored at the plot scale. The occurrence of Iriartea deltoidea, among the most abundant Amazonian tree species despite having been recorded in relatively few floristic inventories, extends its documented distribution within the eastern Amazon.

What does this mean for conservation? Several palm species have been proposed as indicators or surrogates of biodiversity, and describing their local environmental associations may improve their application in conservation planning. The findings may also contribute to future ecological studies, restoration initiatives, and agroforestry projects involving native palms. Yet the researchers are careful about scope. Because these patterns were documented along a single continuous transect, they should be interpreted as spatial associations within the surveyed landscape rather than as general ecological relationships for the eastern Amazon. The relatively low proportion of variation in composition explained by the measured variables suggests that additional factors, such as edaphic properties and disturbance history, also shape the community.

Still, the broader lesson is unmistakable. In one of the least-studied corners of Pará, the architecture of the palm community is written in the language of streams, slopes, and forest greenness. Small watercourses that satellites cannot see turn out to matter more for palm diversity than the perennial rivers that they can, a finding that argues powerfully for the continued value of boots-on-the-ground botany in an age of remote sensing. As the Amazon faces intensifying pressures from deforestation corridors barely a hundred kilometers to the north and east, inventories like this one provide the floristic baselines against which all future change will be measured, and they remind us that the forest’s hidden kingdoms still await their cartographers.

Subject of Research: Palm community structure and its associations with topography, hydrography, and forest structure in the eastern Amazon

Article Title: Palm community structure in relation to topography, hydrography, and forest structure in the Iriri River Valley, eastern Amazon

Article References: Salm, R., von Schmalz, K. E., Ruz, E. H., de Souza, D. A., de Oliveira Jr, A. P., & Pereira, C. C. (2026). Palm community structure in relation to topography, hydrography, and forest structure in the Iriri River Valley, eastern Amazon. Neotropical Biology and Conservation, 21(4), 297-320. https://doi.org/10.3897/neotropical.21.e206551

Image Credits: AI Generated

DOI: 10.3897/neotropical.21.e206551

Keywords: Amazon, palm communities, topography, hydrography, NDVI, species distribution, Palm, community, structure, relation, forest, Iriri

Tracking glacial melt in the ocean is deceptively difficult. When ice shelves and glacier termini discharge freshwater, whether as basal melt from floating ice or as subglacial outflow at grounding lines, that water mixes almost immediately with ambient seawater. Oceanographers traditionally quantify the meltwater fraction using tracer-based calculations that require predefined source water types, known as end members. In the classic approach, an analyst assumes the ocean can be described as a mixture of a small number of pure inputs, for example warm deep water, winter-modified shelf water, and pure glacial melt, each with known temperature and salinity. The meltwater fraction is then inferred from the leftover properties that cannot be explained by the mixing of those assumed sources.

The problem is that the answers depend heavily on the choices made. Pick a different deep-water definition, adjust the salinity of the meltwater end member, or allow a glacial ice end member in addition to liquid melt, and the estimated freshwater fractions can shift substantially. In regions with complex hydrography, where Antarctic Bottom Water formation, modified Circumpolar Deep Water intrusions, and seasonal sea-ice processes all compete to shape water properties, the ambiguity grows worse. East Antarctica, with its thousands of kilometers of ice front and sparse observations, has been especially vulnerable to these methodological uncertainties, leaving the meltwater budget of the region poorly constrained.

The new study sidesteps the end-member problem entirely. Rather than prescribing source water types and solving for their proportions, the researchers developed a parameterization that identifies glacier-derived freshwater directly from the structure of the hydrographic data itself. The technique exploits the fact that glacial melt alters temperature and salinity along characteristic lines in property space: because meltwater enters the ocean at the freezing point and carries negligible salt, its addition moves water masses in predictable directions in temperature-salinity coordinates. By parameterizing these trajectories without fixing the end points, the method estimates the freshwater contribution at every measured depth, producing not just a surface map but a three-dimensional reconstruction of the meltwater field.

The reconstruction is built from the vast archive of hydrographic observations collected across the East Antarctic shelf and slope, including conductivity-temperature-depth profiles, seal-mounted sensor data, and ship-based measurements gathered over multiple decades. Each profile is processed to separate the meltwater signal from other processes that also modify salinity, such as sea-ice formation and melting, precipitation, and the intrusion of off-shelf water masses. The end-member-independent framework then assembles these individual column estimates into a continuous three-dimensional field, resolved in longitude, latitude, and depth, that captures the horizontal pathways and vertical distribution of glacier-derived freshwater around the continent’s eastern half.

The resulting picture is striking. Meltwater is not distributed uniformly along the coast. Instead, the reconstruction shows concentrated lenses and layers of freshwater that accumulate at intermediate depths, typically well below the surface, where melt-laden water spreads neutrally according to its density. Along several major glacier systems, plumes of meltwater extend tens to hundreds of kilometers offshore, following the contours of shelf banks and canyon systems that steer the flow. In some locations the freshwater signal reaches the upper slope, hinting that glacial melt from East Antarctica may be exported into the broader Southern Ocean circulation rather than being trapped locally over the shelf, as older, two-dimensional assessments often implied.

These vertical details carry significant implications for ocean physics and climate. Freshwater stabilizes the water column by reducing surface density, which suppresses vertical mixing and can alter the formation of dense shelf waters that ultimately feed Antarctic Bottom Water, a key component of the global overturning circulation. By quantifying where melt accumulates at depth, the reconstruction allows scientists to test whether meltwater is interfering with bottom-water formation sites, potentially weakening the engine that ventilates the deep ocean and stores carbon and heat on centennial timescales. The three-dimensional view also provides essential validation data for ocean and coupled climate models, which historically have struggled to represent meltwater pathways realistically and often rely on crude runoff schemes at the ice-ocean boundary.

The methodological advance is as important as the observational findings. Because the parameterization does not require users to specify source water properties, it can be applied consistently across regions and through time, enabling fair comparisons between sectors of Antarctica and between different observational eras. Consistency is precisely what large-scale budget studies need: aggregating meltwater estimates produced with different end-member choices has been a persistent obstacle to constructing a continent-wide picture. An end-member-independent approach also reduces the risk of circular reasoning, in which assumptions about meltwater properties determine the meltwater fraction that is then used to infer melt rates. The authors show that their framework yields robust meltwater distributions under a range of environmental conditions, offering a template that can be transferred to other glacier-influenced seas.

For East Antarctica specifically, the study arrives at a pivotal moment. Long considered more stable than the marine-terminating glaciers of West Antarctica, the eastern ice sheet is increasingly showing signs of change, with warm modified deep water reaching the flanks of some major ice shelves and several basins identified as potential candidates for future accelerated retreat. A reliable reconstruction of where glacier-derived freshwater already enters the ocean provides both a baseline against which future change can be measured and a diagnostic of which systems are presently discharging melt at elevated rates. If meltwater export from the region strengthens, the three-dimensional fields produced by this method will help determine how quickly that signal propagates into the abyssal circulation.

The work also demonstrates how reanalysis of existing observations can yield new science without new expeditions. Decades of shipboard hydrography and the growing record of instrumented seals have created an underexploited treasure trove for the Southern Ocean; the challenge has been extracting subtle signals, like glacial freshwater, from noisy, unevenly sampled data. By turning a long-standing methodological weakness, the dependence on assumed source waters, into a solved problem, the researchers have converted scattered profiles into a coherent, multidimensional dataset of one of climate science’s most consequential tracers. As observations accumulate and parameterization techniques mature, the approach promises continuously updated maps of Antarctic meltwater, giving scientists and policymakers a clearer view of how the ice sheet, the ocean, and the global climate system are entangling beneath the surface of the far South.

Subject of Research: Three-dimensional mapping of glacier-derived freshwater in East Antarctic coastal waters using an end-member-independent hydrographic method

Article Title: Three-dimensional reconstruction of glacier-derived freshwater in East Antarctica using an end-member-independent hydrographic parameterization

Article References: Watanabe, Y. W., Hirano, D., Ohashi, Y., Sugita, M., Nakano, Y., Makabe, R., & Mizobata, K. (2026). Three-dimensional reconstruction of glacier-derived freshwater in East Antarctica using an end-member-independent hydrographic parameterization. Nature Communications, 17(1), Article 9498. https://doi.org/10.1038/s41467-026-77441-z

Image Credits: AI Generated

DOI: 10.1038/s41467-026-77441-z

Keywords: East Antarctica, glacial meltwater, hydrography, ice shelf melt, Southern Ocean, Antarctic Bottom Water, temperature-salinity analysis, sea-level rise, ocean circulation, water mass analysis, climate change, oceanography

The study, led by Somayeh Abdollahi of Ghent University together with colleagues from Flanders Hydraulics, the National Geographic Institute, and the Agency for Maritime and Coastal Services, harnessed twenty-three years of high-resolution water-level observations recorded between 2001 and 2023. Four tide gauges anchored the analysis: Nieuwpoort and Oostende on the open coast, the harbour station at Zeebrugge, and the offshore Westhinder platform, designated MP7, which stands sentinel over the shipping approaches far from land. For the first time in Belgian waters, the team systematically compared multiple analysis windows, reconstruction strategies, and tidal-constituent selection procedures, all within the widely used UTide harmonic-analysis framework developed by David Codiga at the University of Rhode Island.

Harmonic analysis is the classical engine of tidal science, tracing its lineage back to Arthur Doodson’s landmark 1921 decomposition of the tide-generating potential. The idea is elegant: the tide at any location can be represented as the sum of dozens of sinusoidal constituents, each with a precisely known astronomical frequency but an amplitude and phase that must be estimated from observations. Once those amplitudes and phases are locked in, the tide can be predicted forward or backward in time indefinitely. LAT is then defined as the lowest water level the astronomical tide alone can produce, excluding the chaotic contributions of storm surges and weather. In practice, however, the answer depends heavily on which constituents are included, how long a record is analysed, and how the reconstruction is performed.

The researchers confronted this methodological fragility head-on. They tested four period-selection approaches, labelled M1 through M4: a single long-record harmonic analysis spanning 2001 to 2019 with annual reconstructions; separate year-by-year analyses; annual analyses used to predict the following nineteen years; and a sophisticated complex averaging of annual harmonic coefficients that preserves the vector nature of tidal amplitudes and phases. Each method was applied at all four stations, and the spread in the resulting LAT estimates turned out to be far from trivial, underscoring that a single ‘correct’ LAT value is an illusion unless the analytical recipe is specified precisely.

Equally consequential was the choice of tidal constituents. The team compared three selection strategies: the default constituent set offered by UTide, a fixed set traditionally used by the Dutch hydrographic service, and a novel iterative procedure guided by the spectrum of the analysis residuals. In the residual-spectrum-guided approach, the researchers analysed the tide, examined the frequency content of what remained unexplained, and added constituents corresponding to prominent residual peaks, repeating the process until the spectrum was clean. This data-driven strategy, implemented entirely with UTide routines, emerged as the preferred configuration, capable of capturing shallow-water constituents that generic or borrowed constituent sets overlook in the complicated hydrodynamics of the southern North Sea.

The preferred configuration, based on the 2001–2019 input window, yielded deterministic LAT values of −0.648 metres at Nieuwpoort, −0.499 metres at Oostende, −0.251 metres at Zeebrugge, and −0.300 metres at the offshore Westhinder station, all expressed relative to the Belgian national height datum TAW/DNG. The pronounced spatial variation, with LAT nearly half a metre deeper at Nieuwpoort than at Zeebrugge, reflects the complex interplay of coastal geometry, bathymetry, and shallow-water tidal dynamics along the Belgian coast. No single conversion grid, the researchers caution, can faithfully represent such station-specific behaviour without careful local calibration.

To quantify how confident one can be in these numbers, the team propagated uncertainty in the harmonic parameters using 1000 Gaussian Monte Carlo realisations, perturbing the estimated amplitudes and phases according to their statistical uncertainties and recomputing the LAT minimum for each realisation. The resulting conditional standard deviations were remarkably tight, ranging from 0.0031 to 0.0034 metres, suggesting that once the harmonic model and constituent set are fixed, the LAT estimate is numerically robust. The authors are careful, however, to frame this spread as conditional on the chosen model rather than a complete uncertainty budget; it captures parameter uncertainty but not the full spectrum of errors that could affect an operational datum.

The Monte Carlo experiments also delivered a reassuring message about computational practicality. Comparing ensembles of 500 and 1000 realisations across the four period-selection methods produced maximum absolute differences of only 0.007 metres in the Monte Carlo mean, 0.009 metres in the conditional standard deviation, and 0.020 metres in one percentile-interval bound. For the preferred residual-spectrum-guided configuration, those differences shrank further to 0.003, 0.004, and 0.006 metres respectively. In other words, halving the ensemble size barely moves the answer, a finding that will ease the computational burden for hydrographic offices considering similar analyses elsewhere.

Throughout the study, the researchers treated existing operational products, including station values published by the Agency for Maritime and Coastal Services and the official LAT conversion grid, as consistency benchmarks rather than independent ground truth. This distinction matters: the operational values were not used to validate the new estimates in an absolute sense, but to check whether the updated analysis remains compatible with the datums that Belgian mariners and coastal engineers rely on today. The deterministic LAT, defined as the minimum of the unperturbed harmonic reconstruction, and the Monte Carlo mean, the average of minima across perturbed realisations, were deliberately reported separately, since the two quantities answer subtly different questions.

The broader significance of the work extends beyond Belgium’s short coastline. The 18.61-year nodal cycle, in which the lunar orbit’s precession modulates tidal amplitudes by up to tens of centimetres, means that short analysis windows can bias LAT estimates substantially, a phenomenon documented in tide-gauge records worldwide. By demonstrating the value of a nineteen-year input window, roughly one full nodal cycle, and by providing a reproducible framework for constituent selection and uncertainty propagation, the Belgian team has offered a template that other hydrographic services can adapt. As sea levels rise and extreme water levels evolve, the low end of the tidal envelope deserves the same scrutiny as the high end, and this study shows exactly how that scrutiny should be applied. The tide-gauge data underpinning the analysis remain freely available through the Meetnet Vlaamse Banken portal, inviting researchers everywhere to test, refine, and extend the approach.

Subject of Research: Assessment of Lowest Astronomical Tide at Belgian North Sea tide gauges using long-term harmonic analysis and residual-spectrum-guided constituent selection

Article Title: Assessment of lowest astronomical tide at belgian north sea tide gauges using long-term harmonic analysis and residual-spectrum-guided constituent selection

Article References: Abdollahi, S., Vanlede, J., Verbeurgt, J., Verstraeten, J., De Kuyper, A., Gurdebeke, P., & De Wulf, A. (2026). Assessment of lowest astronomical tide at belgian north sea tide gauges using long-term harmonic analysis and residual-spectrum-guided constituent selection. Ocean Dynamics, 76(9), Article 96. https://doi.org/10.1007/s10236-026-01850-3

Image Credits: AI Generated

DOI: 10.1007/s10236-026-01850-3

Keywords: Lowest Astronomical Tide, Belgian North Sea, harmonic analysis, tide gauges, tidal constituents, UTide, Monte Carlo uncertainty, vertical datum, ocean dynamics, hydrography, coastal sciences, residual spectrum