Advective Flux Research Articles

Does the Operation of a Reservoir Alter Its Interactions with the Atmosphere? Investigating the Role of Advective Fluxes on Energy and Hydrological Balances of the Romaine-2 Subarctic Hydropower Reservoir

Abstract The hydrological processes of cascading hydroelectric reservoirs differ from those of lakes, due to the importance of the inflows and outflows that vary with energy demand. These heat and water advection terms are rarely considered in water body energy balance analyses even though reservoirs are common man-made structures, especially in North America, and thus may affect the regional climate. This study provides a comprehensive assessment of the water and energy balance of the 85-km2 Romaine-2 northern reservoir (50.69°N, 63.24°W), mean depth of 44 m, highlighting the significant contribution of the advection heat fluxes. The water balance input was primarily controlled by upstream (turbine) inflows (77.6%), while lateral (natural) inflows and direct precipitation represented 21.2% and 1.2%, respectively. As for the reservoir’s heat budget, the net advection of heat accounted on average for 25.0% of the input, of which net radiation was the largest component (73.3%). After accounting for the absence of energy balance closure, latent heat and sensible heat fluxes represented 73.2% and 25.1% of total energy output from the reservoir, respectively. The thermal regime was influenced by the hydrological flow conditions, which were regulated by reservoir management. This played a major role in the evolution of the thermocline and the temperature of the epilimnion, and ultimately, in the dynamics of the turbulent heat fluxes. This study suggests that the heat advection term represents a large fraction of the heat budget of northern reservoirs and should be properly considered.

Modelling diagenetic reactions and secondary porosity generation in sandstones controlled by the advection of low‐molecular‐weight organic acids

AbstractHigher secondary porosity was observed in the centre of a sandstone unit in the Eocene Shahejie Formation fan delta front sandstones from the Bozhong Depression, Bohai Bay Basin. This differs from past studies showing secondary porosity mainly in the marginal parts of sandstones adjacent to shales. This study utilized reactive transport models involving low‐molecular‐weight organic acids (LMWOA) to discuss potential processes resulting in the contrary distribution of secondary porosity. An interface model simulating LMWOA diffusion from adjacent shales to the sandstone resulted in secondary porosity in sandstones adjacent to shales. In contrast, an advection model simulating advective transport of LMWOA parallel to the sandstone bedding successfully generated higher secondary porosity in the central part. The central part of the sandstone exhibited better grain sorting (greater depositional porosity) and significantly less early carbonate cements compared to the marginal sandstone parts. Consequently, the central part had greater porosity prior to the dissolution through LMWOA. The initially higher porosity in the central part allowed for a higher advective flux of LMWOA‐rich water and associated lower pH, resulting in decreased oligoclase saturation, higher oligoclase dissolution rates and ultimately higher secondary porosity. This study indicates that grain sorting during sediment deposition, early carbonate cementation, LMWOA production in adjacent shales, and advection processes collectively control the diagenetic reactions and the distribution of secondary porosity in sandstones.

Processes behind the decrease in Congo Basin precipitation during the rainy seasons inferred from ERA‐5 reanalysis

AbstractThe ongoing Congo Basin drying is linked to decreasing atmospheric water vapour. To comprehend the processes reducing precipitation, we analysed water and Moist Static Energy components, associated with heating sources using ERA5 reanalysis data over the period 1981–2022. The findings reveal that the precipitation deficit in the Congo Basin during March–April–May and September–October–November is primarily attributed to a decreasing trend in the vertical moisture advection anomaly induced by vertical velocity anomalies. Following this, a horizontal moisture advection anomaly induced by specific humidity anomalies plays a significant role. However, evaporation, horizontal moisture advection anomalies induced by wind anomalies, and vertical moisture advection anomalies induced by specific humidity anomalies do not appear to contribute significantly to the drying trend in the region. The horizontal propagation of MSE advection anomalies exhibits a similar pattern to that of moisture flux advection anomalies, both at different boundaries and within the basin itself. Conversely, the contribution of vertical advection of MSE anomalies remains weak. By integrating moisture and MSE budgets, we interpret the declining precipitation trend as a consequence of uplift driven by horizontal MSE advection. Furthermore, our analysis indicates a decreasing trend in the atmospheric heating source within the basin in association with reduced precipitation. This reduction could signify a decrease in moisture influx into the basin, thereby providing more explanations for the observed drying.

Quantifying the Contribution of Ocean Advection and Surface Flux to the Upper‐Ocean Salinity Variability Resolved by Climate Model Simulations

AbstractThis study examines the impact of ocean advection and surface freshwater flux on the non‐seasonal, upper‐ocean salinity variability in two climate model simulations with eddy‐resolving and eddy‐parameterized ocean components (HR and LR, respectively). We assess the realism of each simulation by comparing their sea surface salinity (SSS) variance with satellite and Argo float estimates. In the extratropics, the HR variance is about five times larger than that in LR and agrees with Argo. In turn, the extratropical satellite SSS variance is smaller than that from HR and Argo by about a factor of two, potentially caused by the insufficient resolution of radiometers to capture mesoscale features and their low sensitivity to SSS in cold waters. Using a simplified salinity conservation equation for the upper‐50‐m ocean, we find that the advection‐driven variance in HR is, on average, 10 times larger than the surface flux‐driven variance, reflecting the action of mesoscale processes.

Predicting the Timing of the Solar Cycle 25 Polar Field Reversal

The process of the Sun’s polar field cancellation reversal commences with the emergence of new cycle Hale’s polarity active regions. Once the Sun undergoes polarity reversal, typically occurring near the peak of solar activity, it begins the process of accumulating the seed field for the forthcoming solar cycle. In recent years, the advective flux transport (AFT) model has proven highly effective in forecasting the progression of polar fields by leveraging observations of surface flows and magnetic flux emergence. In this study, we make use of the predictive capability of the AFT model to simulate the evolution of the polar fields and estimate the timing of the Solar Cycle 25 polarity reversal in both hemispheres of the Sun. We use the statistical properties of active regions along with Solar Cycle 13, which closely resembles the current solar cycle (Cycle 25), to generate synthetic active regions in order to simulate future magnetic flux emergence in AFT to predict the evolution of the polar field. Based on our simulations, we anticipate that the northern hemisphere of the Sun will undergo a polarity reversal between 2024 June and November, with the center of our distribution at 2024 August. In the southern hemisphere, we anticipate a polarity reversal between 2024 November and 2025 July, centered around 2025 February. Additionally, assuming that the reversal of the axial dipole moment coincides with the peak of the solar cycle, our findings indicate that Cycle 25 is expected to peak in 2024 (likely between 2024 April and August).

Lake surface cooling drives littoral-pelagic exchange of dissolved gases.

The extent of littoral influence on lake gas dynamics remains debated in the aquatic science community due to the lack of direct quantification of lateral gas transport. The prevalent assumption of diffusive horizontal transport in gas budgets fails to explain anomalies observed in pelagic gas concentrations. Here, we demonstrate through high-frequency measurements in a eutrophic lake that daily convective horizontal circulation generates littoral-pelagic advective gas fluxes one order of magnitude larger than typical horizontal fluxes used in gas budgets. These lateral fluxes are sufficient to redistribute gases at the basin-scale and generate concentration anomalies reported in other lakes. Our observations also contrast the hypothesis of pure, nocturnal littoral-to-pelagic exchange by showing that convective circulation transports gases such as oxygen and methane toward both the pelagic and littoral zones during the daytime. This study challenges the traditional pelagic-centered models of aquatic systems by showing that convective circulation represents a fundamental lateral transport mechanism to be integrated into gas budgets.

A portable low-cost device to quantify advective gas fluxes from mofettes into the lower atmosphere: First application to Starzach mofettes (Germany)

In this study, we introduce a portable low-cost device for in situ gas emission measurement from focused point sources of CO2, such as mofettes. We assess the individual sensors’ precision with calibration experiments and perform an independent verification of the system’s ability to measure gas flow rates in the range of liters per second. The results from one week of continuous CO2 flow observation from a wet mofette at the Starzach site is presented and correlated with the ambient meteorological dynamics. In the observed period, the gas flow rate of the examined mofette exhibits a dominant cycle of around four seconds that is linked to the gas rising upwards through a water column. We find the examined mofette to have a daily emission of 465 kg ±16 %. Furthermore, two events were observed that increased the flow rate abruptly by around 25 % within only a few minutes and a decaying period of 24 hours. These types of events were previously observed by others at the same site but dismissed as measurement errors. We discuss these events as a hydrogeological phenomenon similar to cold-water geyser eruptions. For meteorological events like the passages of high pressure fronts with steep changes in atmospheric pressure, we do not see a significant correlation between atmospheric parameters and the rate of gas exhalation in our one-week time frame, suggesting that on short timescales the atmospheric pumping effect plays a minor role for wet mofettes at the Starzach site.

Experimental investigation of scalar dispersion in indoor spaces

The scalar dispersion from point sources in indoor spaces is experimentally investigated using simultaneous particle–image velocimetry and planar laser-induced fluorescence techniques in a 20:1 and a 60:1 full-to-model scale room model. The ventilation inlets dominate turbulence production, with magnitudes of the velocities and Reynolds stresses observed to increase with air changes per hour (ACH). Mean concentration maps show a dependence on the ACH and source location which is attributed to the flow field at the near-source region. The peak-to-mean concentration shows a weak dependence on the mean concentration and concentration variance maps, indicating risk for toxic chemicals may be underpredicted if based only on these information. The concentration PDFs are generally well-described by exponential distributions with C99/crms′ values never exceeding 5.0. The magnitudes of the advective and turbulent scalar fluxes are strongly dependent on the ACH and source location, neither of which are able to dominate the other by more than an order of magnitude. The eddy diffusivity tensor was measured and a conditional-averaging based method is proposed to approximate it to an isotropic eddy diffusion coefficient, K. For real applications where K is used to estimate magnitudes of the turbulent scalar flux using the gradient transport model, the assumption of isotropic turbulence can introduce an uncertainty of around 17.8%.

Mechanisms of Non-Fresh Groundwater Presence at Water Tables in Highly Permeable Coastal Aquifers.

Coastal aquifers with high hydraulic conductivities on the order of 10-2 m s-1 have unconventional salinity distributions with the presence of non-fresh groundwater at the water table over a wide swath near the coast. This study aims to unravel the mechanisms underlying the phenomenon via numerical simulations for variably saturated, density-driven flow and solute transport in porous media. The simulation results indicate that the existence of non-fresh groundwater at the water table is attributed to the upward mass flux in the saturated zone near the coast, which transports solute from deeper groundwater toward the water table. With high hydraulic conductivity, the upward mass flux becomes prominent at shallower elevations because of the high Darcy flux and the shallow saline groundwater. The upward mass flux has two main drivers, upward advection by the upward flow component and transverse dispersion by the seaward flow component. The advective mass flux dominates over the transverse dispersion in the deep part of the saturated zone where only groundwater with sea water salinity exists. In contrast, the transverse dispersion becomes more pronounced than the upward advection in the shallow saturated zone just beneath the water table and in the unsaturated zone immediately above the water table. Our findings help interpret the unconventional salinity distributions observed and elucidate the unique dynamics of groundwater flow and solute transport in highly permeable coastal aquifers.

Nonlinear diffusion‐limited two‐dimensional colliding‐plume simulations with very high‐order numerical approximations

AbstractAtmospheric numerical models play a crucial role in operational weather forecasting, as well as in improving our understanding of atmospheric dynamics via research studies. Maximising their accuracy is of paramount importance. Use of advective flux schemes greater than O7 in atmospheric models is largely undocumented, with no studies considering O3–O17 fluxes with formal accuracy‐preserving high‐order interpolation, pressure gradient/divergence, and subgrid‐scale (SGS) turbulent fluxes. Higher‐order numerical approximations can reduce truncation, amplitude and phase errors, and potentially improve model accuracy and effective resolution. Here, simulations are presented using very‐high‐order O3–O17 fluxes, with/without high‐order O2–O18 Lagrangian interpolations, pressure gradient/divergence approximations, and SGS turbulent fluxes for a two‐dimensional, highly‐viscous (Re ∼ 100) diffusion‐limited, nonlinear colliding‐plumes problem using 25–200 m spatial resolutions. The highest‐order flux schemes coupled with higher‐order interpolations, pressure gradient/divergence and SGS flux approximations produced the best solutions, with higher‐order fluxes and interpolations being most important. Overall solution convergence of order about O1–O2 with mode‐split (fast sound/slow advective waves) O3 Runge–Kutta temporal schemes was negatively impacted by order at most O1 temporal convergence with SGS fluxes, divergence damping, and especially spatial filters, compared to about order O3 convergence with these inactivated. While very‐high‐order schemes were shown to improve solution accuracy, few cost‐effective higher‐order highly‐viscous test problem solutions (higher order vs finer resolution) were found using theoretical floating‐point operations (FPO) with Courant‐Friedrichs‐Lewy (CFL)‐limited or constant stable Courant number‐based time steps. However, employing central processing unit (CPU) time, rather than FPOs, demonstrated there was reduced computational burden using higher‐order approximations. We conclude that O9–O17 flux schemes with or without high‐order (≥O4) interpolations, pressure gradient/divergence approximations, and SGS fluxes can improve atmospheric model solution accuracy, without prohibitive computation costs, compared to O3–O7 flux with O2 interpolations, pressure gradient/divergence approximations, and SGS fluxes.

Investigating Atmospheric Responses to and Mechanisms Governing North Atlantic Sea Surface Temperatures over 10-Year Periods

Abstract North Atlantic sea surface temperature (SST) variability plays a critical role in modulating the climate system. However, characterizing patterns of North Atlantic SST variability and diagnosing the associated mechanisms is challenging because they involve coupled atmosphere–ocean interactions with complex spatiotemporal relationships. Here we address these challenges by applying a time-evolving self-organizing map approach to a long preindustrial coupled control simulation and identify a variety of 10-yr spatiotemporal evolutions of winter SST anomalies, including but not limited to those associated with the North Atlantic Oscillation–Atlantic multidecadal variability (NAO–AMV)-like interactions. To assess mechanisms and atmospheric responses associated with various SST spatiotemporal evolutions, composites of atmospheric and oceanic variables associated with these evolutions are investigated. Results show that transient-eddy activities and atmospheric circulation responses exist in almost all the evolutions that are closely correlated to the details of the SST pattern. In terms of the mechanisms responsible for generating various SST evolutions, composites of ocean heat budget terms demonstrate that contributions to upper-ocean temperature tendency from resolved ocean advection and surface heat fluxes rarely oppose each other over 10-yr periods in the subpolar North Atlantic. We further explore the potential for predictability for some of these 10-yr SST evolutions that start with similar states but end with different states. However, we find that these are associated with abrupt changes in atmospheric variability and are unlikely to be predictable. In summary, this study broadly investigates the atmospheric responses to and the mechanisms governing the North Atlantic SST evolutions over 10-yr periods. Significance Statement Climate variability in the North Atlantic Ocean has wide-ranging impacts on global and regional climate. However, the processes involved include interactions between the ocean and atmosphere that vary across both space and time, making it challenging to characterize and predict. Using a novel machine learning approach, this study identifies various time evolutions of North Atlantic sea surface temperature patterns over 10-yr periods. This includes evolutions with similar start states but different trajectories, which have important implications for predictability. Furthermore, we investigate the mechanisms responsible for these evolutions and how different sea surface temperature patterns affect atmospheric circulation through small-scale atmospheric disturbances. These new insights into the complex ocean–atmosphere interactions over time are critical for improving decadal prediction skill.

Nutricline adjustment by internal tidal beam generation enhances primary production in idealized numerical models

When the barotropic tide encounters variable bathymetry, fluctuating flow along a topographic slope generates baroclinic tides, or internal tides. There is growing evidence that these internal tides can affect primary production in the euphotic zone, though the dominant mechanisms are unclear. Internal tides move passive phytoplankton through an exponentially varying light field, enhancing primary production near the base of the euphotic zone. In addition internal tides also increase primary production through vertical nutrient advection into the euphotic zone. Topographically generated internal tides can be separated into two regimes: 1) the often highly nonlinear near-field regime where tidal beams are observed and 2) the more linear far-field regime. This study examines the primary production response to these internal tide processes using the Regional Ocean Modeling System (ROMS) coupled to a simple Nutrient, Phytoplankton, Zooplankton, Detritus (NPZD) model configured for an oligotrophic system with the nutricline positioned below 50 m depth. These idealized simulations generate internal tide beams with an oscillating, horizontal body force at the M2 tidal frequency that is applied to domains with a bathymetric step and uniform stratification. Sensitivity of the primary production response to the energy content of the tidal beam is obtained by adjusting the height and slope of the bathymetric step. Simulation results reveal that primary production intensifies along tidal beams due to the local enhancement of parcel vertical displacement (light effect) and nutrient advective flux divergence (nutrient effect). In the near-field regime across the range of step heights and slopes in this study, the nutrient effect is an order of magnitude larger and explains 92% of the variance in primary production versus only 14% for the light effect. The geometry of the generating feature sets the kinematics of the tidal beam. The light effect is limited in the euphotic zone across our domains because realized changes in light experienced over a tidal cycle are small relative to the amount of light available at a particular depth. In contrast, the magnitude of the nutrient effect increases more substantially with tidal beam energy.

Diffusioosmotic dispersion of solute in a long narrow channel

Solute–surface interactions have garnered considerable interest in recent years as a novel control mechanism for driving unique fluid dynamics and particle transport with potential applications in fields such as biomedicine, the development of microfluidic devices and enhanced oil recovery. In this study, we will discuss dispersion induced by the diffusioosmotic motion near a charged wall in the presence of a solute concentration gradient. Here, we introduce a plug of salt with a Gaussian distribution at the centre of a channel with no background flow. As the solute diffuses, the concentration gradient drives a diffusioosmotic slip flow at the walls, which results in a recirculating flow in the channel; this, in turn, drives an advective flux of the solute concentration. This effect leads to cross-stream diffusion of the solute, altering the effective diffusivity of the solute as it diffuses along the channel. We derive theoretical predictions for the solute dynamics using a multiple-time-scale analysis to quantify the dispersion driven by the solute–surface interactions. Furthermore, we derive a cross-sectionally averaged concentration equation with an effective diffusivity analogous to that from Taylor dispersion. In addition, we use numerical simulations to validate our theoretical predictions.

Water, heat and salt balance of the south-eastern part of the Barents Sea

The south-eastern part of the Barents Sea is located away from the main currents, with a combination of climatic, hydrological and oceanological processes creating conditions that make the region different from the rest of the Barents Sea such that it is seen as a separate region and called sometimes the Pechora Sea. Despite the intensive economic activity in the south-eastern part of the Barents Sea, it is not yet clear to what extent the general Atlantic water transport in the Barents Sea and, consequently, the transport of heat and salt, affects this region. Therefore, the aim of this study was to assess advective flows at open boundaries, as well as other components of the water, heat and salt balances of the south-eastern part of the Barents Sea. Based on monthly average data from the MERCATOR GLORYS12V1 reanalysis for the period 1993–2018, we calculated water transport, heat and salt flows at the boundaries of the south-eastern part of the Barents Sea (at 50° E in the west, at 71° N in the north and in the Kara Gate Strait); to close the balances, an assessment was made of sea-atmosphere interaction characteristics on the sea surface based on ECMWF ERA5 reanalysis data. Water, heat and salt balances were combined with a residual not exceeding 1.6 %. Linear trends for the characteristics obtained were calculated. It is revealed that the average long-term resulting water transport in the south-eastern part of the Barents Sea is directed from the north-west of the region to the Kara Gate Strait (0.40 Sv). This current is associated with the Atlantic waters and also carries heat and salt. The resulting heat input (5.92 TW) creates a heat excess in the water area, which is compensated for by interaction with the atmosphere (1.86 TW). The salt flow through the region is estimated at 13.98 t/s. During the study period, all the main flows have a statistically significant positive trend in the incoming and outgoing parts of the balances: water transport — 0.005 Sv per year; salt flow — 0.18 t/s per year. This indicates an increase in the transit of Atlantic waters through the south-eastern part of the Barents Sea. An increase in the advective heat flux (0.15 TW per year) across the western border is accompanied by an increase in its release into the atmosphere (0.07 TW per year) and an increase in evaporation of 6.9 mm per year. Sea levels are also rising at a rate of 0.27 cm per year. Thus, the increasing dynamics of the processes in the region is a factor to take into account in conducting economic activities.

The Current Feedback on Stress Modifies the Ekman Buoyancy Flux at Fronts

Abstract Ocean surface currents introduce variations into the surface wind stress that can change the component of the stress aligned with the thermal wind shear at fronts. This modifies the Ekman buoyancy flux, such that the current feedback on the stress tends to generate an effective flux of buoyancy and potential vorticity to the mixed layer. Scaling arguments and idealized simulations resolving both mesoscale and submesoscale turbulence suggest that this pathway for air–sea interaction can be important both locally at individual submesoscale fronts with strong surface currents—where it can introduce equivalent advective heat fluxes exceeding several hundred watts per square meter—and in the spatial mean where it reduces the integrated Ekman buoyancy flux by approximately 50%. The accompanying source of surface potential vorticity injection suggests that at some fronts the current feedback modification of the Ekman buoyancy flux may be significant in terms of both submesoscale dynamics and boundary layer energetics, with an implied modification of symmetric instability growth rates and dissipation that scales similarly to the energy lost through the negative wind work generated by the current feedback. This provides an example of how the shift of dynamical regimes into the submesoscale may promote the importance of air–sea interaction mechanisms that differ from those most active at larger scale.

Concurrent assessment of diffusive and advective PAH movement strongly affected by temporal and spatial changes

Chemical movement influences exposure, remediation and interventions. Understanding chemical movement in addition to chemical concentrations at contaminated sites is critical to informed decision making. Using seepage meters and passive sampling devices we assessed both diffusive and advective flux of bioavailable polycyclic aromatic hydrocarbons (PAHs) at three time points, across two seasons, at a former creosote site in St. Helens, Oregon, United States. To our knowledge, this is the first time both diffusive and advective fluxes have been measured simultaneously at a contaminated site. Concentrations of 39 parent PAHs were determined by gas chromatography triple quadrupole mass spectrometry. Across both seasons and all sites, diffusive flux of PAHs was up to three orders of magnitude larger than advective flux. Release of PAHs from sediments and water were identified, likely from legacy contamination, as well as deposition from the air into the site from contemporary and other sources. The majority of PAH movement was comprised of three and four ring PAHs. Chemical movement on the site was found to be spatially and temporally variable. Volatilization decreased and atmospheric deposition increased from summer to fall. At the locations with higher levels of contamination, sum PAH release from sediments decreased by more than two orders of magnitude from summer to late fall. These data reflect the spatial heterogeneity and temporal variability of this site and demonstrate the importance of seasonality in assessing chemical movement at contaminated sites. Results from this study can inform future legacy site assessments to optimize remediation strategies and assess remediation effectiveness.

Helical and non-helical large-scale dynamos in thin accretion discs

ABSTRACT The dynamics of accreting and outgoing flows around compact objects depends crucially on the strengths and configurations of the magnetic fields therein, especially of the large-scale fields that remain coherent beyond turbulence scales. Possible origins of these large-scale magnetic fields include flux advection and disc dynamo actions. However, most numerical simulations have to adopt an initially strong large-scale field rather than allow them to be self-consistently advected or amplified, due to limited computational resources. The situation can be partially cured by using sub-grid models where dynamo actions only reachable at high resolutions are mimicked by artificial terms in low-resolution simulations. In this work, I couple thin-disc models with local shearing-box simulation results to facilitate more realistic sub-grid dynamo implementations. For helical dynamos, detailed spatial profiles of dynamo drivers inferred from local simulations are used, and the non-linear quenching and saturation is constrained by magnetic helicity evolution. In the inner disc region, saturated fields have dipole configurations and the plasma β reaches ≃0.1 to 100, with correlation lengths ≃h in the vertical direction and ≃10 h in the radial direction, where h is the disc scale height. The dynamo cycle period is ≃40 orbital time scale, compatible with previous global simulations. Additionally, I explore two dynamo mechanisms which do not require a net kinetic helicity and have only been studied in shearing-box set-ups. I show that such dynamos are possible in thin accretion discs, but produce field configurations that are incompatible with previous results. I discuss implications for future general-relativistic magnetohydrodynamics simulations.

Diffusive and advective cross-frontal fluxes of inorganic nutrients and dissolved inorganic carbon in the Barents Sea in autumn

The Atlantic Water, entering the Arctic through the Barents Sea and Fram Strait, is the main source of nutrients in the Arctic Ocean. The Barents Sea is divided by the Polar Front into an Atlantic-dominated domain in the south, and an Arctic-dominated domain in the north. The Polar Front is a thermohaline structure, which is topographically-steered at sub-surface, and influenced by the seasonal sea ice edge near the surface. Exchanges of nutrients between the inflowing Atlantic Water and the surrounding waters are key for the primary production in the Barents Sea. In October 2020, we measured nutrients (nitrate, phosphate and silicic acid), dissolved inorganic carbon (DIC), ocean stratification, currents and turbulence in the vicinity of the Polar Front in the Barents Sea within the framework of the Nansen Legacy project, allowing estimates of horizontal and vertical advective fluxes and turbulent fluxes of nitrate and DIC. We studied the autumn situation when primary production was declining. We found a substantial transfer of nitrate and DIC across the Polar Front from the Atlantic domain to the Arctic domain. Up to one quarter of the replenishment of the nitrate in the mixed layer during winter could be attributed to vertical mixing during wind events, shared approximately equally between advective and turbulent fluxes. The vertical turbulent fluxes bring nutrients from the subsurface Atlantic Water to the surface. We also identified an export of nitrate and DIC from the Barents Sea to the Nordic Seas occurring along the eastern shelf of Svalbard. Our study shows the role of vertical fluxes in fall and winter to precondition for the following spring bloom.

Iron Limitation and Biogeochemical Effects in Southern California Current Coastal Upwelling Filaments

AbstractIn the spring and summer, high rates of primary production occur in the California Current Ecosystem (CCE) when nutrients are supplied to the euphotic zone. During periods of intense coastal upwelling, a flux of the micronutrient iron comes from nearshore sedimentary sources. In this upwelling region, mesoscale filament features distribute iron laterally, leading to distinct iron‐influenced ecological zones. This study is the first to focus on the biogeochemical links between iron, the macronutrients, and particulates in coastal upwelling filaments. Broad spatial patterns of iron and biogenic silica concentrations, and proxies of iron‐stress of diatoms, support results from microcosm amendment studies conducted during CCE Long Term Ecological Research process cruises in the summers of 2017 and 2019. We found that the benthic boundary layer and shoreward filament endmember supply dissolved and total dissolvable iron to this feature, but rapid assimilation and sinking by biogenic particles (e.g., diatoms) depletes the surface concentrations. Subsequently, diatom blooms which form in recently upwelled water masses become iron limited over time, thereby affecting the ratios of surface macronutrient reservoirs and biogeochemical advective fluxes. The development of Fe‐limitation during lateral advection may lead to efficient carbon export downstream and offshore of the region with the highest phytoplankton growth rates and productivity.

Avoiding methane emission rate underestimates when using the divergence method

Methane is a powerful greenhouse gas, and a primary target for mitigating climate change in the short-term future due to its relatively short atmospheric lifetime and greater ability to trap heat in Earth’s atmosphere compared to carbon dioxide. Top-down observations of atmospheric methane are possible via drone and aircraft surveys as well as satellites such as the TROPOspheric Monitoring Instrument (TROPOMI). Recent work has begun to apply the divergence method to produce regional methane emission rate estimates. Here we show that when the divergence method is applied to spatially incomplete observations of methane, it can result in negatively biased time-averaged regional emission rates. We show that this effect can be counteracted by adopting a procedure in which daily advective fluxes of methane are time-averaged before the divergence method is applied. Using such a procedure with TROPOMI methane observations, we calculate yearly Permian emission rates of 3.1, 2.4 and 2.7 million tones per year for the years 2019 through 2021. We also show that highly-resolved plumes of methane can have negatively biased estimated emission rates by the divergence method due to the presence of turbulent diffusion in the plume, but this is unlikely to affect regional methane emission budgets constructed from TROPOMI observations of methane. The results from this work are expected to provide useful guidance for future implementations of the divergence method for emission rate estimation from satellite data—be it for methane or other gaseous species in the atmosphere.