Pelagic-feeding seabirds deliver nutrient subsidies that enhance the productivity, biodiversity, and resilience of terrestrial and marine ecosystems, particularly in nutrient-poor tropical environments. However, the biogeophysical variables governing the fluxes of these nutrients within and among interconnected ecosystems remain poorly understood. To address this, we examined the spatial distribution of seabird-vectored nutrients in the seascape of Tetiaroa, a semi-enclosed coral atoll in French Polynesia, where seabird populations and associated nutrient cycles are recovering after recent rat eradication. We focus on the nitrogen isotope (δ15N) signatures of a dominant marine alga as evidence of seabird-vectored nutrient uptake. Integrating stable isotope analysis within a seascape ecology framework, we show that breeding seabird biomass, depth, distance to land, geographic location within the atoll, and seafloor curvature drive spatial patterns of nutrient enrichment. Specifically, our models account for up to 88% of the variation in algal δ15N signatures and reveal peak enrichment in shallow, nearshore areas where water flow slows and converges due to localised seafloor curvature. These results extend previous research by highlighting seafloor geomorphology, notably curvature, as a modulator of fine-scale nutrient delivery patterns. Although a complex model incorporating 11 high-resolution biogeophysical variables enhanced spatial predictions by revealing fine-scale variations, a simpler model using only 5 of these variables was comparably effective in capturing overall spatial trends. This study identifies the key seascape configuration and complexity characteristics likely to affect the spatial patterns of recovery potential following the restoration of seabird-driven nutrient cycles, offering valuable guidance for ongoing restoration efforts in this coupled island-reef system. Future investigations could assess how the effects of biogeophysical variables on nutrient delivery vary in magnitude and direction across different geographic, geological, and anthropogenic contexts.

Accurately tracing the sources and fate of excess PO43- in waterways is necessary for sustainable catchment management. The natural abundance isotopic composition of O in PO43- (δ18OP) is a promising tracer of point source pollution, but its ability to track diffuse agricultural pollution is unclear. We tested the hypothesis that δ18OP could distinguish between agricultural PO43- sources by measuring the integrated δ18OP composition and P speciation of contrasting inorganic fertilisers (compound v rock) and soil textures (sand, loam, clay). δ18OP composition differed between the three soil textures sampled across six working livestock farms: sandy soils had lower overall δ18OP values (21 ± 1 ‰) than the loams (23 ± 1 ‰), which corresponded with a smaller, but more readily leachable, PO43- pool. Fertilisers had greater δ18OP variability (~8‰) driven by both fertiliser type and manufacturing year. Upscaling these values showed that ‘agricultural soil leaching’ δ18OP signatures could span from 18 – 25 ‰, and are influenced by both fertiliser type and the time between application and leaching. These findings emphasise the potential of δ18OP to untangle soil-fertiliser P dynamics under controlled conditions, but that its use to trace catchment-scale agricultural PO43- losses is limited by uncertainties in soil biological P cycling and its associated isotopic fractionation.