Spatial synchrony cascades across ecosystem boundaries and up food webs via resource subsidies.

Cross-ecosystem subsidies are critical to ecosystem structure and function, especially in recipient ecosystems where they are the primary source of organic matter to the food web. Subsidies are indicative of processes connecting ecosystems and can couple ecological dynamics across system boundaries. However, the degree to which such flows can induce cross-ecosystem cascades of spatial synchrony, the tendency for system fluctuations to be correlated across locations, is not well understood. Synchrony has destabilizing effects on ecosystems, adding to the importance of understanding spatiotemporal patterns of synchrony transmission. In order to understand whether and how spatial synchrony cascades across the marine-terrestrial boundary via resource subsidies, we studied the relationship between giant kelp forests on rocky nearshore reefs and sandy beach ecosystems that receive resource subsidies in the form of kelp wrack (detritus). We found that synchrony cascades from rocky reefs to sandy beaches, with spatiotemporal patterns mediated by fluctuations in live kelp biomass, wave action, and beach width. Moreover, wrack deposition synchronized local abundances of shorebirds that move among beaches seeking to forage on wrack-associated invertebrates, demonstrating that synchrony due to subsidies propagates across trophic levels in the recipient ecosystem. Synchronizing resource subsidies likely play an underappreciated role in the spatiotemporal structure, functioning, and stability of ecosystems.

Tail-dependent spatial synchrony arises from nonlinear driver-response relationships.

Spatial synchrony may be tail-dependent, that is, stronger when populations are abundant than scarce, or vice-versa. Here, 'tail-dependent' follows from distributions having a lower tail consisting of relatively low values and an upper tail of relatively high values. We present a general theory of how the distribution and correlation structure of an environmental driver translates into tail-dependent spatial synchrony through a non-linear response, and examine empirical evidence for theoretical predictions in giant kelp along the California coastline. In sheltered areas, kelp declines synchronously (lower-tail dependence) when waves are relatively intense, because waves below a certain height do little damage to kelp. Conversely, in exposed areas, kelp is synchronised primarily by periods of calmness that cause shared recovery (upper-tail dependence). We find evidence for geographies of tail dependence in synchrony, which helps structure regional population resilience: areas where population declines are asynchronous may be more resilient to disturbance because remnant populations facilitate reestablishment.

Disturbance and nutrients synchronise kelp forests across scales through interacting Moran effects.

Spatial synchrony is a ubiquitous and important feature of population dynamics, but many aspects of this phenomenon are not well understood. In particular, it is largely unknown how multiple environmental drivers interact to determine synchrony via Moran effects, and how these impacts vary across spatial and temporal scales. Using new wavelet statistical techniques, we characterised synchrony in populations of giant kelp Macrocystis pyrifera, a widely distributed marine foundation species, and related synchrony to variation in oceanographic conditions across 33 years (1987-2019) and >900 km of coastline in California, USA. We discovered that disturbance (storm-driven waves) and resources (seawater nutrients)-underpinned by climatic variability-act individually and interactively to produce synchrony in giant kelp across geography and timescales. Our findings demonstrate that understanding and predicting synchrony, and thus the regional stability of populations, relies on resolving the synergistic and antagonistic Moran effects of multiple environmental drivers acting on different timescales.

Spectrally simplified approach for leveraging legacy geostationary oceanic observations.

The use of multispectral geostationary satellites to study aquatic ecosystems improves the temporal frequency of observations and mitigates cloud obstruction, but no operational capability presently exists for the coastal and inland waters of the United States. The Advanced Baseline Imager (ABI) on the current iteration of the Geostationary Operational Environmental Satellites, termed the R Series (GOES-R), however, provides sub-hourly imagery and the opportunity to overcome this deficit and to leverage a large repository of existing GOES-R aquatic observations. The fulfillment of this opportunity is assessed herein using a spectrally simplified, two-channel aquatic algorithm consistent with ABI wave bands to estimate the diffuse attenuation coefficient for photosynthetically available radiation, Kd(PAR). First, an in situ ABI dataset was synthesized using a globally representative dataset of above- and in-water radiometric data products. Values of Kd(PAR) were estimated by fitting the ratio of the shortest and longest visible wave bands from the in situ ABI dataset to coincident, in situKd(PAR) data products. The algorithm was evaluated based on an iterative cross-validation analysis in which 80% of the dataset was randomly partitioned for fitting and the remaining 20% was used for validation. The iteration producing the median coefficient of determination (R2) value (0.88) resulted in a root mean square difference of 0.319m-1, or 8.5% of the range in the validation dataset. Second, coincident mid-day images of central and southern California from ABI and from the Moderate Resolution Imaging Spectroradiometer (MODIS) were compared using Google Earth Engine (GEE). GEE default ABI reflectance values were adjusted based on a near infrared signal. Matchups between the ABI and MODIS imagery indicated similar spatial variability (R2=0.60) between ABI adjusted blue-to-red reflectance ratio values and MODIS default diffuse attenuation coefficient for spectral downward irradiance at 490 nm, Kd(490), values. This work demonstrates that if an operational capability to provide ABI aquatic data products was realized, the spectral configuration of ABI would potentially support a sub-hourly, visible aquatic data product that is applicable to water-mass tracing and physical oceanography research.

Climate-driven regime shifts in a mangrove-salt marsh ecotone over the past 250 years.

Climate change is driving the tropicalization of temperate ecosystems by shifting the range edges of numerous species poleward. Over the past few decades, mangroves have rapidly displaced salt marshes near multiple poleward mangrove range limits, including in northeast Florida. It is uncertain whether such mangrove expansions are due to anthropogenic climate change or natural climate variability. We combined historical accounts from books, personal journals, scientific articles, logbooks, photographs, and maps with climate data to show that the current ecotone between mangroves and salt marshes in northeast Florida has shifted between mangrove and salt marsh dominance at least 6 times between the late 1700s and 2017 due to decadal-scale fluctuations in the frequency and intensity of extreme cold events. Model projections of daily minimum temperature from 2000 through 2100 indicate an increase in annual minimum temperature by 0.5 °C/decade. Thus, although recent mangrove range expansion should indeed be placed into a broader historical context of an oscillating system, climate projections suggest that the recent trend may represent a more permanent regime shift due to the effects of climate change.

Evaluating regional resiliency of coastal wetlands to sea level rise through hypsometry-based modeling.

Sea level rise (SLR) threatens coastal wetlands worldwide, yet the fate of individual wetlands will vary based on local topography, wetland morphology, sediment dynamics, hydrologic processes, and plant-mediated feedbacks. Local variability in these factors makes it difficult to predict SLR effects across wetlands or to develop a holistic regional perspective on SLR response for a diversity of wetland types. To improve regional predictions of SLR impacts to coastal wetlands, we developed a model that addresses the scale-dependent factors controlling SLR response and accommodates different levels of data availability. The model quantifies SLR-driven habitat conversion within wetlands across a region by predicting changes in individual wetland hypsometry. This standardized approach can be applied to all wetlands in a region regardless of data availability, making it ideal for modeling SLR response across a range of scales. Our model was applied to 105 wetlands in southern California that spanned a broad range of typology and data availability. Our findings suggest that if wetlands are confined to their current extents, the region will lose 12% of marsh habitats (vegetated marsh and unvegetated flats) with 0.6 m of SLR (projected for 2050) and 48% with 1.7 m of SLR (projected for 2100). Habitat conversion was more drastic in wetlands with larger proportions of marsh habitats relative to subtidal habitats and occurred more rapidly in small lagoons relative to larger sites. Our assessment can inform management of coastal wetland vulnerability, improve understanding of the SLR drivers relevant to individual wetlands, and highlight significant data gaps that impede SLR response modeling across spatial scales. This approach augments regional SLR assessments by considering spatial variability in SLR response drivers, addressing data gaps, and accommodating wetland diversity, which will provide greater insights into regional SLR response that are relevant to coastal management and restoration efforts.

Global patterns of kelp forest change over the past half-century.

Kelp forests (Order Laminariales) form key biogenic habitats in coastal regions of temperate and Arctic seas worldwide, providing ecosystem services valued in the range of billions of dollars annually. Although local evidence suggests that kelp forests are increasingly threatened by a variety of stressors, no comprehensive global analysis of change in kelp abundances currently exists. Here, we build and analyze a global database of kelp time series spanning the past half-century to assess regional and global trends in kelp abundances. We detected a high degree of geographic variation in trends, with regional variability in the direction and magnitude of change far exceeding a small global average decline (instantaneous rate of change = -0.018 y-1). Our analysis identified declines in 38% of ecoregions for which there are data (-0.015 to -0.18 y-1), increases in 27% of ecoregions (0.015 to 0.11 y-1), and no detectable change in 35% of ecoregions. These spatially variable trajectories reflected regional differences in the drivers of change, uncertainty in some regions owing to poor spatial and temporal data coverage, and the dynamic nature of kelp populations. We conclude that although global drivers could be affecting kelp forests at multiple scales, local stressors and regional variation in the effects of these drivers dominate kelp dynamics, in contrast to many other marine and terrestrial foundation species.

Satellite sensor requirements for monitoring essential biodiversity variables of coastal ecosystems.

The biodiversity and high productivity of coastal terrestrial and aquatic habitats are the foundation for important benefits to human societies around the world. These globally distributed habitats need frequent and broad systematic assessments, but field surveys only cover a small fraction of these areas. Satellite-based sensors can repeatedly record the visible and near-infrared reflectance spectra that contain the absorption, scattering, and fluorescence signatures of functional phytoplankton groups, colored dissolved matter, and particulate matter near the surface ocean, and of biologically structured habitats (floating and emergent vegetation, benthic habitats like coral, seagrass, and algae). These measures can be incorporated into Essential Biodiversity Variables (EBVs), including the distribution, abundance, and traits of groups of species populations, and used to evaluate habitat fragmentation. However, current and planned satellites are not designed to observe the EBVs that change rapidly with extreme tides, salinity, temperatures, storms, pollution, or physical habitat destruction over scales relevant to human activity. Making these observations requires a new generation of satellite sensors able to sample with these combined characteristics: (1) spatial resolution on the order of 30 to 100-m pixels or smaller; (2) spectral resolution on the order of 5 nm in the visible and 10 nm in the short-wave infrared spectrum (or at least two or more bands at 1,030, 1,240, 1,630, 2,125, and/or 2,260 nm) for atmospheric correction and aquatic and vegetation assessments; (3) radiometric quality with signal to noise ratios (SNR) above 800 (relative to signal levels typical of the open ocean), 14-bit digitization, absolute radiometric calibration <2%, relative calibration of 0.2%, polarization sensitivity <1%, high radiometric stability and linearity, and operations designed to minimize sunglint; and (4) temporal resolution of hours to days. We refer to these combined specifications as H4 imaging. Enabling H4 imaging is vital for the conservation and management of global biodiversity and ecosystem services, including food provisioning and water security. An agile satellite in a 3-d repeat low-Earth orbit could sample 30-km swath images of several hundred coastal habitats daily. Nine H4 satellites would provide weekly coverage of global coastal zones. Such satellite constellations are now feasible and are used in various applications.

Operationalizing the social-ecological systems framework to assess sustainability.

Environmental governance is more effective when the scales of ecological processes are well matched with the human institutions charged with managing human-environment interactions. The social-ecological systems (SESs) framework provides guidance on how to assess the social and ecological dimensions that contribute to sustainable resource use and management, but rarely if ever has been operationalized for multiple localities in a spatially explicit, quantitative manner. Here, we use the case of small-scale fisheries in Baja California Sur, Mexico, to identify distinct SES regions and test key aspects of coupled SESs theory. Regions that exhibit greater potential for social-ecological sustainability in one dimension do not necessarily exhibit it in others, highlighting the importance of integrative, coupled system analyses when implementing spatial planning and other ecosystem-based strategies.

Fluctuations in population fecundity drive variation in demographic connectivity and metapopulation dynamics.

Demographic connectivity is vital to sustaining metapopulations yet often changes dramatically through time due to variation in the production and dispersal of offspring. However, the relative importance of variation in fecundity and dispersal in determining the connectivity and dynamics of metapopulations is poorly understood due to the paucity of comprehensive spatio-temporal data on these processes for most species. We quantified connectivity in metapopulations of a marine foundation species (giant kelp Macrocystis pyrifera) across 11 years and approximately 900 km of coastline by estimating population fecundity with satellite imagery and propagule dispersal using a high-resolution ocean circulation model. By varying the temporal complexity of different connectivity measures and comparing their ability to explain observed extinction-colonization dynamics, we discovered that fluctuations in population fecundity, rather than fluctuations in dispersal, are the dominant driver of variation in connectivity and contribute substantially to metapopulation recovery and persistence. Thus, for species with high variability in reproductive output and modest variability in dispersal (most plants, many animals), connectivity measures ignoring fluctuations in fecundity may overestimate connectivity and likelihoods of persistence, limiting their value for understanding and conserving metapopulations. However, we demonstrate how connectivity measures can be simplified while retaining utility, validating a practical solution for data-limited systems.

Poleward expansion of mangroves is a threshold response to decreased frequency of extreme cold events.

Regional warming associated with climate change is linked with altered range and abundance of species and ecosystems worldwide. However, the ecological impacts of changes in the frequency of extreme events have not been as well documented, especially for coastal and marine environments. We used 28 y of satellite imagery to demonstrate that the area of mangrove forests has doubled at the northern end of their historic range on the east coast of Florida. This expansion is associated with a reduction in the frequency of "extreme" cold events (days colder than -4 °C), but uncorrelated with changes in mean annual temperature, mean annual precipitation, and land use. Our analyses provide evidence for a threshold response, with declining frequency of severe cold winter events allowing for poleward expansion of mangroves. Future warming may result in increases in mangrove cover beyond current latitudinal limits of mangrove forests, thereby altering the structure and function of these important coastal ecosystems.

Connectivity structures local population dynamics: a long-term empirical test in a large metapopulation system.

Ecological theory predicts that demographic connectivity structures the dynamics of local populations within metapopulation systems, but empirical support has been constrained by major limitations in data and methodology. We tested this prediction for giant kelp Macrocystis pyrifera, a key habitat-forming species in temperate coastal ecosystems worldwide, in southern California, USA. We combined a long-term (22 years), large-scale (~500 km coastline), high-resolution census of abundance with novel patch delineation methods and an innovative connectivity measure incorporating oceanographic transport and source fecundity. Connectivity strongly predicted local dynamics (well-connected patches had lower probabilities of extinction and higher probabilities of colonization, leading to greater likelihoods of occupancy) but this relationship was mediated by patch size. Moreover, the relationship between connectivity and local population dynamics varied over time, possibly due to temporal variation in oceanographic transport processes. Surprisingly, connectivity had a smaller influence on colonization relative to extinction, possibly because local ecological factors differ greatly between extinct and extant patches. Our results provide the first comprehensive evidence that southern California giant kelp populations function as a metapopulation system, challenging the view that populations of this important foundation species are governed exclusively by self-replenishment.

Integrating physiological threshold experiments with climate modeling to project mangrove species' range expansion.

Predictions of climate-related shifts in species ranges have largely been based on correlative models. Due to limitations of these models, there is a need for more integration of experimental approaches when studying impacts of climate change on species distributions. Here, we used controlled experiments to identify physiological thresholds that control poleward range limits of three species of mangroves found in North America. We found that all three species exhibited a threshold response to extreme cold, but freeze tolerance thresholds varied among species. From these experiments, we developed a climate metric, freeze degree days (FDD), which incorporates both the intensity and the frequency of freezes. When included in distribution models, FDD accurately predicted mangrove presence/absence. Using 28 years of satellite imagery, we linked FDD to observed changes in mangrove abundance in Florida, further exemplifying the importance of extreme cold. We then used downscaled climate projections of FDD to project that these range limits will move northward by 2.2-3.2 km yr(-1) over the next 50 years.

Patch definition in metapopulation analysis: a graph theory approach to solve the mega-patch problem.

The manner in which patches are delineated in spatially realistic metapopulation models will influence the size, connectivity, and extinction and recolonization dynamics of those patches. Most commonly used patch-definition methods focus on identifying discrete, contiguous patches of habitat from a single temporal observation of species occurrence or from a model of habitat suitability. However, these approaches are not suitable for many metapopulation systems where entire patches may not be fully colonized at a given time. For these metapopulation systems, a single large patch of habitat may actually support multiple, interacting subpopulations. The interactions among these subpopulations will be ignored if the patch is treated as a single unit, a situation we term the "mega-patch problem." Mega-patches are characterized by variable intra-patch synchrony, artificially low inter-patch connectivity, and low extinction rates. One way to detect this problem is by using time series data to calculate demographic synchrony within mega-patches. We present a framework for identifying subpopulations in mega-patches using a combination of spatial autocorrelation and graph theory analyses. We apply our approach to southern California giant kelp (Macrocystis pyrifera) forests using a new, long-term (27 years), satellite-based data set of giant kelp canopy biomass. We define metapopulation patches using our method as well as several other commonly used patch delineation methodologies and examine the colonization and extinction dynamics of the metapopulation under each approach. We find that the relationships between patch characteristics such as area and connectivity and the demographic processes of colonizations and extinctions vary among the different patch-definition methods. Our spatial-analysis/graph-theoretic framework produces results that match theoretical expectations better than the other methods. This approach can be used to identify subpopulations in metapopulations where the distributions of organisms do not always reflect the distribution of suitable habitat.

Balancing multiple management objectives as climate change transforms ecosystems.

As climate change facilitates significant and persistent ecological transformations, managing ecosystems according to historical baseline conditions may no longer be feasible. The Resist-Accept-Direct (RAD) framework can guide climate-informed management interventions, but in its current implementations RAD has not yet fully accounted for potential tradeoffs between multiple - sometimes incompatible - ecological and societal goals. Key scientific challenges for informing climate-adapted ecosystem management include (i) advancing our predictive understanding of transformations and their socioecological impacts under novel climate conditions, and (ii) incorporating uncertainty around trajectories of ecological change and the potential success of RAD interventions into management decisions. To promote the implementation of RAD, practitioners can account for diverse objectives within just and equitable participatory decision-making processes.

A marine heatwave changes the stabilizing effects of biodiversity in kelp forests.

Biodiversity can stabilize ecological communities through biological insurance, but climate and other environmental changes may disrupt this process via simultaneous ecosystem destabilization and biodiversity loss. While changes to diversity-stability relationships (DSRs) and the underlying mechanisms have been extensively explored in terrestrial plant communities, this topic remains largely unexplored in benthic marine ecosystems that comprise diverse assemblages of producers and consumers. By analyzing two decades of kelp forest biodiversity survey data, we discovered changes in diversity, stability, and their relationships at multiple scales (biological organizational levels, spatial scales, and functional groups) that were linked with the most severe marine heatwave ever documented in the North Pacific Ocean. Moreover, changes in the strength of DSRs during/after the heatwave were more apparent among functional groups than both biological organizational levels (population vs. ecosystem levels) and spatial scales (local vs. broad scales). Specifically, the strength of DSRs decreased for fishes, increased for mobile invertebrates and understory algae, and were unchanged for sessile invertebrates during/after the heatwave. Our findings suggest that biodiversity plays a key role in stabilizing marine ecosystems, but the resilience of DSRs to adverse climate impacts primarily depends on the functional identities of ecological communities.

Synchrony in dynamics of giant kelp forests is driven by both local recruitment and regional environmental controls.

Populations of many species display spatially synchronous fluctuations in abundance. Synchrony is most commonly attributed to three processes: factors that influence recruitment (e.g., dispersal, early survival), large-scale environmental variability, and spatially autocorrelated trophic interactions. However it is often difficult to link population synchrony to a specific dominant process, particularly when multiple synchronizing forces are operating. We utilized a new satellite-based data set of giant kelp (Macrocystis pyrifera) canopy biomass to examine population synchrony in southern California kelp forests on spatial scales ranging from 50 m to 300 km and temporal scales ranging from 1 to 11 years. We examined the relationship between synchrony and distance for adult kelp populations, kelp recruits, sea urchin abundance (a major grazer of kelp), and environmental variables known to influence kelp population dynamics. Population synchrony in giant kelp decreased with distance between populations: an initial rapid exponential decrease between 50 m and 1.3 km was followed by a second, large-scale decrease between distances of 1.3 km and 172 km. The 50-m to 1.3-km spatial scale corresponded to the scales of synchrony in the abundance of sea urchins and young kelp recruits, suggesting that local drivers of predation and recruitment influence small-scale synchrony in kelp populations. The spatial correlation patterns of environmental variables, particularly wave height, were similar to the synchrony-distance relationship of kelp populations from 1.3 km to 172 km, suggesting that regional environmental variability, i.e., the Moran effect, was the dominant process affecting synchrony at larger spatial scales. This two-step pattern in the relationship between kelp biomass synchrony and distance was apparent in each of the 11 years of our study. Our results highlight the potential for synthesizing approaches from both landscape and population ecology in order to identify the multiple processes that generate synchrony in population dynamics.

A comparison of global mangrove maps: Assessing spatial and bioclimatic discrepancies at poleward range limits.

Mangrove distribution maps are used for a variety of applications, ranging from estimates of mangrove extent, deforestation rates, quantify carbon stocks, to modelling response to climate change. There are multiple mangrove distribution datasets, which were derived from different remote sensing data and classification methods, and so there are some discrepancies among these datasets, especially with respect to the locations of their range limits. We investigate the latitudinal discrepancies in poleward mangrove range limits represented by these datasets and how these differences translate climatologically considering factors known to control mangrove distributions. We compare four widely used global mangrove distribution maps - the World Atlas of Mangroves, the World Atlas of Mangroves 2, the Global Distribution of Mangroves, the Global Mangrove Watch. We examine differences in climate among 21 range limit positions by analysing a set of bioclimatic variables that have been commonly related to the distribution of mangroves. Global mangrove maps show important discrepancies in the position of poleward range limits. Latitudinal differences between mangrove range limits in the datasets exceed 5°, 7° and 10° in western North America, western Australia and northern West Africa, respectively. In some range limit areas, such as Japan, discrepancies in the position of mangrove range limits in different datasets correspond to differences exceeding 600 mm in annual precipitation and > 10 °C in the minimum temperature of the coldest month. We conclude that dissimilarities in mapping mangrove range limits in different parts of the world can jeopardise inferences of climatic thresholds. We expect that global mapping efforts should prioritise the position of range limits with greater accuracy, ideally combining data from field-based surveys and very high-resolution remote sensing data. An accurate representation of range limits will contribute to better predicting mangrove range dynamics and shifts in response to climate change.

Kelpwatch: A new visualization and analysis tool to explore kelp canopy dynamics reveals variable response to and recovery from marine heatwaves.

Giant kelp and bull kelp forests are increasingly at risk from marine heatwave events, herbivore outbreaks, and the loss or alterations in the behavior of key herbivore predators. The dynamic floating canopy of these kelps is well-suited to study via satellite imagery, which provides high temporal and spatial resolution data of floating kelp canopy across the western United States and Mexico. However, the size and complexity of the satellite image dataset has made ecological analysis difficult for scientists and managers. To increase accessibility of this rich dataset, we created Kelpwatch, a web-based visualization and analysis tool. This tool allows researchers and managers to quantify kelp forest change in response to disturbances, assess historical trends, and allow for effective and actionable kelp forest management. Here, we demonstrate how Kelpwatch can be used to analyze long-term trends in kelp canopy across regions, quantify spatial variability in the response to and recovery from the 2014 to 2016 marine heatwave events, and provide a local analysis of kelp canopy status around the Monterey Peninsula, California. We found that 18.6% of regional sites displayed a significant trend in kelp canopy area over the past 38 years and that there was a latitudinal response to heatwave events for each kelp species. The recovery from heatwave events was more variable across space, with some local areas like Bahía Tortugas in Baja California Sur showing high recovery while kelp canopies around the Monterey Peninsula continued a slow decline and patchy recovery compared to the rest of the Central California region. Kelpwatch provides near real time spatial data and analysis support and makes complex earth observation data actionable for scientists and managers, which can help identify areas for research, monitoring, and management efforts.