The response of inorganic carbon distributions and dynamics to upwelling-favorable winds on the northern Gulf of Mexico during summer
Explore the source record for details and available documents.
SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Explore the source record for details and available documents.
Pulse jet mixed process vessels comprise a key element of the U.S. Department of Energy’s strategy to process millions of gallons of legacy nuclear waste slurries. Slurry suctioned into a pulse jet mixer (PJM) tube at the end of one pulse is pneumatically driven from the PJM toward the bottom of the vessel at the beginning of the next pulse, forming a jet. The jet front traverses the distance from nozzle outlet to the bottom of the vessel and spreads out radially. Varying numbers of PJMs are typically arranged in a ring configuration within the vessel at a selected radius and operated concurrently. Centrally directed radial flows from neighboring jets collide to create a central upwell that elevates the solids in the center of the vessel when the PJM tubes expel their contents. An essential goal of PJM operation is to elevate solids to the liquid surface to minimize stratification. Solids stratification may adversely affect throughput of the waste processing plant. Unacceptably high slurry densities at the base of the vessel may plug the pipeline through which the slurry exits the vessel. Additionally, chemical reactions required for processing may not achieve complete conversion. To avoid these conditions, a means of predicting the elevation to which the solids rise in the central upwell that can be used during vessel design remains essential. In this paper we present a simple criterion to evaluate the extent of solids elevation achieved by a turbulent upwell jet. The criterion asserts that at any location in the central upwell the local velocity must be in excess of a cutoff velocity to remain turbulent. We find that local velocities in excess of 0.6 m/s are necessary for turbulent jet flow through both Newtonian and yield stress slurries. By coupling this criterion with the free jet velocity equation relating the local velocity to elevation in the central upwell, we estimate the elevation at which turbulence fails, and consequently the elevation at which the upwell fails to further lift the slurry. Comparing this elevation to the vessel fill level predicts whether the jet flow will achieve the full vertical extent of the vessel at the center. This simple local-velocity criterion determines a minimum PJM nozzle velocity at which the full vertical extent of the central upwell in PJM vessels will be turbulent. The criterion determines a minimum because flow in regions peripheral to the central upwelling jet may not be turbulent, even when the center of the vessel in the upwell is turbulent, if the jet pulse duration is too short. The local-velocity criterion ensures only that there is sufficient wherewithal for the turbulent jet flow to drive solids to the surface in the center of the vessel in the central upwell.
Meteorological data are being collected at several points within four intensive study areas in Barrow. These data assist in the calculation of the energy balance at the land surface and are also useful as inputs into modeling activities. (https://doi.org/10.5440/1164893). Measurements at AREA A (all heights listed are relative to the land surface): Location of the Site A station: N71* 16.888', W156* 36.647' * Air temperature at three heights: 1.4 meters, 2.0 meters, 3.5 meters * Relative Humidity / Dew Point at two heights: 1.4 meters, 3.5 meters * Wind Speed at two heights: 1.3 meters, 3.5 meters * 4 Component Net Radiation @ 1.2 meters ** Upwelling & Downwelling Shortwave Radiation ** Upwelling & Downwelling Longwave Radiation * Summer Precipitation Site A is located on polygonal permafrost and co-located with continuous snow depth measurements and subsurface measurements. Measurements at AREA B (all heights listed are relative to the land surface): Location of the Site B station: N71* 16.785', W156* 36.211' * Air temperature at three heights: 1.4 meters, 2.0 meters, 3.5 meters * Relative Humidity / Dew Point at two heights: 1.4 meters, 3.5 meters * Wind Speed at two heights: 1.3 meters, 3.5 meters * 4 Component Net Radiation @ 1.2 meters ** Upwelling & Downwelling Shortwave Radiation ** Upwelling & Downwelling Longwave Radiation * Summer Precipitation Site B is located on polygonal permafrost and co-located with continuous snow depth measurements and subsurface measurements. ------------------------------------------------- Measurements at AREA C (all heights listed are relative to the land surface): Location of the Site C station: N71* 16.804', W156* 36.072' * Air temperature at three heights: 1.4 meters, 2.0 meters, 3.5 meters * Relative Humidity / Dew Point at two heights: 1.4 meters, 3.5 meters * Wind Speed at two heights: 1.3 meters, 3.5 meters * 4 Component Net Radiation @ 1.2 meters ** Upwelling & Downwelling Shortwave Radiation ** Upwelling & Downwelling Longwave Radiation * Summer Precipitation * Barometric Pressure Site C is located on polygonal permafrost and co-located with continuous snow depth measurements and subsurface measurements. ------------------------------------------------- Measurements at AREA D (all heights listed are relative to the land surface): Location of the Site D station: N71* 16.927', W156* 36.036' * Air temperature at one height: 1.2 meters * Relative Humidity / Dew Point at one height: 1.2 * Wind Speed at one height: 1.2 meters * 4 Component Net Radiation @ 1.0 meters ** Upwelling & Downwelling Shortwave Radiation ** Upwelling & Downwelling Longwave Radiation Site D is located on polygonal permafrost and co-located with continuous snow depth measurements and subsurface measurements. Data are in ASCII files, with columns separated by commas: ','
Modern annual floral and faunal production in the northwest Arabian Sea derives primarily from upwelling induced by strong southwest winds during June, July, and August. Indian Ocean summer monsoon winds are, in turn, driven by differential heating between the Asian continent and the Indian ocean to the south. This differential heating produces a strong pressure gradient resulting in southwest monsoon winds and both coastal and divergent upwelling off the Arabian Peninsula. Over geologic time scales (10(exp 4) to 10(exp 6) years), monsoon wind strength is sensitive to changes in boundary conditions which influence this pressure gradient. Important boundary conditions include the seasonal distribution of solar radiation, global ice volume, Indian Ocean sea surface temperature, and the elevation and albedo of the Asian continent. To the extent that these factors influence monsoon wind strength, they also influence upwelling and productivity. In addition, however, productivity associated with upwelling can be decoupled from the strength of the summer monsoon winds via ocean mechanisms which serve to inhibit or enhance the nutrient supply in the intermediate waters of the Indian Ocean, the source for upwelled waters in the Arabian Sea. To differentiate productivity associated with wind-induced upwelling from that associated with other components of the system such as nutrient sequestering in glacial-age deep waters, we employ a strategy which monitors independent components of the oceanic and atmospheric subsystems. Using sediment records from the Owen Ridge, northwest Arabian Sea, we monitor the strength of upwelling and productivity using two independent indicators, percent G. bulloides and opal accumulation. We monitor the strength of southwest monsoon winds by measuring the grain-size of lithogenic dust particles blown into the Arabian Sea from the surrounding deserts of the Somali and Arabian Peninsulas. Our current hypothesis is that the variability associated with the 41 kyr power in the G. bulloides and opal accumulation records derive from nutrient availability in the intermediate waters which are upwelled via monsoon winds. This hypothesis is testable by comparison with Cd records of intermediate and deep waters of the Atlantic and Indian Ocean.
Our research group comprising personnel from University of Utah/EGI (UU/EGI), University of Nevada at Reno (UNR) and Lawrence Berkeley National Laboratory (LBNL) has been applying magnetotellurics (MT), diagnostic structural affiliations, and soil gas flux and geochemistry to assist in identifying hidden, high-enthalpy geothermal systems in extensional regimes of the U.S. Great Basin. We are specifically looking for upwelling low-resistivity and geologic structures that can carry fluids from magmatic or high-grade metamorphic conditions in the deep crust upward to exploitable depths, and to verify the nature of the deep sources through soil gas geochemical compositions. The project was motivated by prior MT transect coverage of western and central Nevada centered upon the Dixie Valley geothermal system where such favorable indicators were first recognized. The high-angle structures are taken to be fluidized fault zones connecting deep magmatic/metamorphic activity with the geothermal system but the concept requires verification by testing at other systems. The project was set up with a two-phased organization. Phase I was carried out at the McGinness Hills system, central Nevada, where Ormat’s flagship power facility is located and a considerable amount of pre-existing data was available. Transect MT data also showed a strong low-resistivity upwelling originating from interpreted deep crustal magmatic underplating. Controlling structures on production as indicated by Ormat data and new mapping by UNR were favorable to dilatency being an accommodation zone between major normal faults of opposing dip. A 3D MT survey and inversion confirmed the existence of the steep low-resistivity zone dipping ESE toward the deep crust and placed N-S bounds upon the feature. Ormat personnel cooperated with PI Wannamaker in the field in providing fluid samples from the production intervals for He isotope composition. Elevated 3He was verified through analysis at LBNL confirming that there was a magmatic connection into the producing system. High CO2 soil gas flux including possibly metamorphic 13C and 14C component was measured over the area of dilatent structures. Hence, the triad of indicators posed above was confirmed in Phase I. Subsequently, Phase II of the project proceeded in the greenfield Kumiva-Blackrock Desert district of northwestern Nevada to see if a new system could be identified. Transect MT data also showed a strong low-resistivity upwelling originating from interpreted deep crustal magmatic underplating. However, given that there was no exploitable geothermal system recognized at the outset in the region, the collection of new geoscientific data was required. An areal MT survey of 131 sites over the Kumiva-Blackrock Desert district was completed, and a 3D inversion produced using the in-house, DOE-supported finite element inversion algorithm. Particular MT low resistivity upwellings that received followup study occurred under the west flank of the Seven Troughs Range, under the playa immediately west of the Blue Wing Mountains, and under northern Granite Springs Valley. Structural assessment of the project area by Co-I J. Faulds at UNR provided numerous favorable Quaternary fault settings, which were correlated to the MT upwelling structures. Extensive profiling for CO2 soil gas flux over favorable geophysical and geological structure was performed by LBNL personnel. Soil gas flux anomalies generally were not large but did show correlation with MT upwelling structure and favorably dilatent geological structures. Isotope analyses showed presence of possible inorganic/metamorphic 13C but 14C concentrations did not exceed background values. We view the initial concept of a confluence of MT low-resistivity upwelling, favorably dilatent 3D geological structure, and elevated soil gas flux including 13C component to be supported by the further evidence of this project although the indicators in the Phase II study were more diffuse. The northern Granite Springs Valley structure is receiving followup in the UNR Play Fairway Analysis, being in the southern termination of a major normal fault with silicified soil and warm temperature probe results
This study is aimed at modeling mesoscale processes such as upwelling/downwelling and ice edge eddies in the marginal ice zones. A two-dimensional coupled ice-ocean model is used for the study. The ice model is coupled to the reduced gravity ocean model through interfacial stresses. The parameters of the ocean model were chosen so that the dynamics would be nonlinear. The model was tested by studying the dynamics of upwelling. Wings parallel to the ice edge with the ice on the right produce upwelling because the air-ice momentum flux is much greater than air-ocean momentum flux; thus the Ekman transport is greater than the ice than in the open water. The stability of the upwelling and downwelling jets is discussed. The downwelling jet is found to be far more unstable than the upwelling jet because the upwelling jet is stabilized by the divergence. The constant wind field exerted on a varying ice cover will generate vorticity leading to enhanced upwelling/downwelling regions, i.e., wind-forced vortices. Steepening and strengthening of vortices are provided by the nonlinear terms. When forcing is time-varying, the advection terms will also redistribute the vorticity. The wind reversals will separate the vortices from the ice edge, so that the upwelling enhancements are pushed to the open ocean and the downwelling enhancements are pushed underneath the ice.
Upwelling of global deep waters to the sea surface in the Southern Ocean closes the global overturning circulation and is fundamentally important for oceanic uptake of carbon and heat, nutrient resupply for sustaining oceanic biological production, and the melt rate of ice shelves. However, the exact pathways and role of topography in Southern Ocean upwelling remain largely unknown. Here in this paper we show detailed upwelling pathways in three dimensions, using hydrographic observations and particle tracking in high-resolution models. The analysis reveals that the northern-sourced deep waters enter the Antarctic Circumpolar Current via southward flow along the boundaries of the three ocean basins, before spiraling southeastward and upward through the Antarctic Circumpolar Current. Upwelling is greatly enhanced at five major topographic features, associated with vigorous mesoscale eddy activity. Deep water reaches the upper ocean predominantly south of the Antarctic Circumpolar Current, with a spatially nonuniform distribution. The timescale for half of the deep water to upwell from 30° S to the mixed layer is ~60–90 years.
Interannual changes in the strength and seasonal evolution of the 1979 through 1982 surface-level southwest monsoon winds are related to variations in the summer phytoplankton bloom of the northwestern Arabian Sea by synthesis of satellite ocean-color remote sensing with analysis of in-situ hydrographic and meteorological data sets. The 1979-1981 southwest monsoon phytoplankton blooms in the northwest Arabian Sea peaked during August-September, extended from the Omani coast to about 6 E, and appeared to lag the development of open-sea upwelling by at least 1 month. In all 3 years the bloom was driven by spatially distinct upward nutrient fluxes to the euphotic zone forced by the physical processes of coastal upwelling and offshore Ekman pumping. Coastal upwelling was evident from May through September, yielded the most extreme concentrations of phytoplankton biomass, and along the Omani coast was limited in its impact on upper ocean biological variability to the continental shelf. Ekman pumping stimulated the development of a broad open-ocean component of the southwest monsoon phytoplankton bloom oceanward of the Omani shelf. Phytoplankton biomass on the Omani continental shelf was increased during both the early and late phases of the 1980 southwest monsoon due to stronger coastal upwelling under the most intense southwesterly winds of the four summers investigated. Diminished coastal upwelling during the early phase of the weak 1982 southwest monsoon resulted in a coastal bloom that reached a mean phytoplankton-pigment concentration that was 28 percent of that seen in 1980. The lack of a strong regional northwestern Arabian Sea bloom in late summer 1982 is attributed to the development of persistent, shallow temperature stratification that rendered Ekman pumping less effective in driving upward nutrient fluxes.
Eddy-induced upwelling occurs along the western edge of the Gulf Stream between Cape Canaveral, Florida, and Cape Hatteras, North Carolina, in the South Atlantic Bight (SAB). Coastal zone color scanner images of 1-km resolution spanning the period April 13-21, 1979, were processed to examine these eddy features in relation to concurrent shipboard and current/temperature measurements at moored arrays. A quasi-one-dimensional (z), time dependent biological model, using only nitrate as a nutrient source, has been combined with a three-dimensional physical model in an attempt to replicate the observed phytoplankton field at the northward edge of an eddy. The model is applicable only to the SAB south of the Charleston Bump, at approximately 31.5 deg N, since no feature analogous to the bump exists in the model bathymetry. The modeled chlorophyll, nitrate, and primary production fields of the euphotic zone are very similar to those obtained from the satellite and shipboard data at the leading edges of the observed eddies south of the Charleston Bump. The horizontal and vertical simulated fluxes of nitrate and chlorophyll show that only approximately 10% of the upwelled nitrate is utilized by the phytoplankton of the modeled grid box on the northern edge of the cyclone, while approximately 75% is lost horizontally, with the remainder still in the euphotic zone after the 10-day period of the model. Loss of chlorophyll due to sinking is very small in this strong upwelling region of the cyclone. The model is relatively insensitive to variations in the sinking parameterization and the external nitrate and chlorophyll fields but is very sensitive to a reduction of the maximum potential growth rate to half that measured. Given the success of this model in simulating the new production of the selcted upwelling region, other upwelling regions for which measurements or successful models of physical and biological quantities and rates exist could be modeled similarly.
A passive tracer to represent dissolved silicate concentrations, with biologically realistic uptake kinetics, is successfully incorporated into a three-dimensional, eddy-resolving, ocean circulation model of the Indian Ocean. Hypotheses are tested to evaluate physical processes which potentially affect the availability of silicate for diatom production in the Arabian Sea. An alternative mechanism is offered to the idea that open ocean upwelling is primarily responsible for the high, vertical nutrient flux and consequent large-scale phytoplankton bloom in the northwestern Arabian Sea during the southwest monsoon. Model results show that dissolved silicate in surface waters available for uptake by diatoms is primarily influenced by the intensity of nearshore upwelling from soutwest monsoonal wind forcing and by the offshore advective transport of surface waters. The upwelling, which in the model occurs within 200 +/- 50 km of the coast, appears to be a result of a combination of coastal upwelling, Elkman pumping, and divergence of the coastal flow as it turns offshore. Localized intensifications of silicate concentrations appear to be hydrodynamically driven and geographically correlated to coastal topographic features. The absence of diatoms in sediments of the eastern Arabian Basin is consistent with modeled distributional patterns of dissolved silicate resulting from limited westward advection of upwelled coastal waters from the western continental margin of India and rapid uptake of available silicate by diatoms. Concentrations of modeled silicate become sufficiently low to become unavailable for diatom production in the eastern Arabian Sea, a region between 61 deg E and 70 deg E at 8 deg N on the south, with the east and west boundaries converging on the north at approximately 67 deg E, 20 deg N.
Anthropogenic influences have led to a strengthening and poleward shift of westerly winds over the Southern Ocean (SO), especially during austral summer. We use observations, an idealized eddy‐resolving ocean‐sea ice channel model, and a global coupled model to explore the SO response to a step‐change in westerly winds. Previous work hypothesized a two timescale response for sea surface temperature. Initially, Ekman transport cools the surface before sustained upwelling causes warming on decadal timescales. The fast response is robust across our models and the observations: we find Ekman‐driven cooling in the mixed layer, mixing‐driven warming below the mixed layer, and a small upwelling‐driven warming at the temperature inversion. The long‐term response is inaccessible from observations. Neither of our models shows a persistent upwelling anomaly, or long‐term, upwelling‐driven subsurface warming. Mesoscale eddies act to oppose the anomalous wind‐driven upwelling, through a process known as eddy compensation, thereby preventing long‐term warming.
Coastal upwelling regions are one of the most dynamic areas of the world’s oceans. The California and Baja California Coasts are impacted by both coastal upwelling and the California Current, leading to frontal activity that is captured by gradients in both Sea Surface Temperature (SST) and Sea Surface Salinity (SSS). Satellite data are a great source of spatial data to study fronts. However, biases near coastal areas and coarse resolutions can impair its usefulness in upwelling areas. In this work gradients in SST from NASA Multi-Scale Ultra-High Resolution (MUR) and in two SSS products derived from the Soil Moisture Active Passive (SMAP) NASA mission are compared directly with gradients derived from the Saildrone uncrewed vehicles to validate the gradients as well as to assess their ability to detect known frontal features. The three remotely sensed data sets (MURSST/JPL, SMAP SSS/RSS, SMAP SSS) were co-located with the Saildrone data prior to the calculation of the gradients. Wavelet analysis is used to determine how well the satellite derived SST and SSS products are reproducing the Saildrone derived gradients. Overall results indicate the remote sensing products are reproducing features of known areas of coastal upwelling. Differences between the SST and SSS gradients are mainly associated with the limitations of the microwave derived SSS coverage near land and its reduced spatial resolution. The results are promising for using remote sensing data sets to monitor frontal structure along the California Coast and the application to long term changes in coastal upwelling and dynamics.
The Santa Barbara Channel, CA (SBC) is a biodiverse marine ecosystem fueled largely by phytoplankton productivity, and the composition of the phytoplankton community influences the magnitude and fates of this productivity. Here, we create a 22-year monthly time series of phytoplankton biomarker pigment concentrations in the SBC by combining 12 years of high performance liquid chromatography phytoplankton pigment concentrations with bio-optical models and 10 additional years of bio-optical observations. The bio-optical models skillfully predict biomarker pigment concentrations representative of five distinct phytoplankton groups (PGs; diatoms, dinoflagellates, chlorophytes, prymnesiophytes, and picophytoplankton) and resolve seasonal responses to the annual upwelling-relaxation cycle for all PGs except the dinoflagellates. Our observations indicate that nanophytoplankton groups respond most rapidly to seasonal upwelling, followed by diatoms, and then by picophytoplankton as the water column stratifies in the summer. A Regional Ocean Modeling System (ROMS) solution is used to relate advection of different source waters to the observed PG dynamics. The ROMS simulation results suggest that, on seasonal time scales, pronounced cross-SBC differences in PG seasonality are related to cross-SBC differences in source waters. El Niño Southern Oscillation events drive interannual variability in the upwelling response of most PGs. On decadal time scales, dinoflagellate blooms are associated with the warm phase of the North Pacific Gyre Oscillation and anomalous advection of Southern California Bight source waters into the SBC. Taken together, our results provide a novel view of phytoplankton community succession in response to seasonal upwelling by considering the dynamics of pico- and nano-phytoplankton and suggest that regional surface ocean advection plays a substantial role in driving phytoplankton composition in the SBC.
Coastal upwelling of nutrients and metals along eastern boundary currents fuels some of the most biologically productive marine ecosystems. Although iron is a main driver of productivity in many of these regions, iron cycling and acquisition by microbes remain poorly constrained, in part due to the unknown composition of organic ligands that keep bioavailable iron in solution. We investigated organic ligand composition in discrete water samples collected across the highly productive California Coastal upwelling system. Siderophores were observed in distinct nutrient regimes at concentrations ranging from 1 pM to 18 pM. Near the shallow continental shelf, ferrioxamine B was observed in recently upwelled, high chlorophyll surface waters while synechobactins were identified within nepheloid layers at 60–90 m depth. In offshore waters characterized by intermediate chlorophyll, iron, and nitrate concentrations, we found amphibactins and an unknown siderophore with a molecular formula of C 33 H 58 O 8 N 5 Fe. Highest concentrations were measured in the photic zone, however, amphibactins were also found in waters as deep as 1500 m. The distribution of siderophores provides evidence for microbial iron deficiency across a range of nutrient regimes and indicates siderophore production and acquisition is an important strategy for biological iron uptake in iron limited coastal systems. Polydisperse humic ligands were also detected throughout the water column and were particularly abundant near the benthic boundary. Our results highlight the fine-scale spatial heterogeneity of metal ligand composition in an upwelling environment and elucidate distinct sources that include biological production and the degradation of organic matter in suboxic waters.
Accurate multidecadal radiative flux records are vital to understand Arctic amplification and constrain climate model uncertainties. Uncertainty in the NASA Clouds and the Earth’s Radiant Energy System (CERES)-derived irradiances is larger over sea ice than any other surface type and comes from several sources. The year-long Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition in the central Arctic provides a rare opportunity to explore uncertainty in CERES-derived radiative fluxes. First, a systematic and statistically robust assessment of surface shortwave and longwave fluxes was conducted using in situ measurements from MOSAiC flux stations. The CERES Synoptic 1degree (SYN1deg) product overestimates the downwelling shortwave flux by +11.40 Wm –2 and underestimates the upwelling shortwave flux by –15.70 Wm –2 and downwelling longwave fluxes by –12.58 Wm –2 at the surface during summer. In addition, large differences are found in the upwelling longwave flux when the surface approaches the melting point (approximately 0°C). The biases in downwelling shortwave and longwave fluxes suggest that the atmosphere represented in CERES is too optically thin. The large negative bias in upwelling shortwave flux can be attributed in large part to lower surface albedo (–0.15) in satellite footprint relative to surface sensors. Additionally, the results show that the spectral surface albedo used in SYN1deg overestimates albedo in visible and mid-infrared bands. A series of radiative transfer model perturbation experiments are performed to quantify the factors contributing to the differences. The CERES-MOSAiC broadband albedo differences (approximately 20 Wm –2 ) explain a larger portion of the upwelling shortwave flux difference than the spectral albedo shape differences (approximately 3 Wm –2 ). In addition, the differences between perturbation experiments using hourly and monthly MOSAiC surface albedo suggest that approximately 25% of the sea ice surface albedo variability is explained by factors not correlated with daily sea ice concentration variability. Biases in net shortwave and longwave flux can be reduced to less than half by adjusting both albedo and cloud inputs toward observed values. The results indicate that improvements in the surface albedo and cloud data would substantially reduce the uncertainty in the Arctic surface radiation budget derived from CERES data products.
The solar infrared radiation station (SIRS), ground radiometers on stand for upwelling radiation (GNDRAD), and sky radiometers on stand for downwelling radiation (SKYRAD) provide continuous measurements of broadband shortwave (solar) and longwave (atmospheric or infrared) irradiances for downwelling and upwelling components. The following six irradiance measurements are collected from a network of stations to help determine the total radiative flux exchange at each of the primary U.S. Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) field research sites—(Southern Great Plains [SGP], North Slope of Alaska [NSA)]and Eastern North Atlantic [ENA])—as well as the ARM mobile facilities (AMF): Direct normal shortwave (solar beam), Diffuse horizontal shortwave (sky), Global horizontal shortwave (total hemispheric), Upwelling shortwave (reflected), Downwelling longwave (atmospheric infrared), Upwelling longwave (surface infrared).
Hyporheic zones (where surface water and groundwater mix in streambeds) play a critical role in the movement of chemicals within watersheds, particularly in upland headwater streams. Seasonal hydrology controls the expansion and contraction of hyporheic zones, with greater downwelling oxic river water influence during periods of peak river discharge (linked to snowmelt), and greater influence of upwelling anoxic groundwater during periods of low- and base-flow. Through this work, we aimed to quantify upwelling and downwelling fluxes across seasons and assess the impact of these dynamics on microbial community assembly and geochemical gradients at East River, Colorado. Further, we performed high spatial resolution sampling to capture heterogeneous mixing patterns around a characteristic meander on the river and developed reactive transport models to explain seasonal patterns of manganese (Mn) cycling in riverbed sediments. We were able to measure a dominant period of downwelling river water associated with high river discharge in Spring 2017 that led to the mixing of microbial communities across a 60-cm depth profile through the riverbed at three locations around the meander, and higher rates of aerobic respiration via delivery of dissolved oxygen (O2) and dissolved organic carbon (DOC). Conversely, we demonstrated that depth-resolved microbial communities became more distinct and stratified across a depth gradient during periods of low- and base-flow when upwelling groundwater exerted more influence in the riverbed. We also observed an increase in Mn concentrations within deeper streambed sediments during the baseflow season. Reactive transport models were developed to understand seasonal changes in Mn cycling in the streambed because Mn oxides influence the mobility of other metals in streams. In field observations and models, dissolved Mn is flushed from the streambed during spring snowmelt. A shift to upwelling conditions over the subsequent baseflow period allows for groundwater rich in dissolved Mn to mix with oxygenated river water in the shallow subsurface, resulting in net accumulation of Mn-oxides until the bed freezes in winter. We also developed models for hydrograph scenarios with snowmelt events of various size and timing. Our scenarios suggest that in years with less snowpack and a longer baseflow season, Mn oxidation will be favored in the upper riverbed sediments over more of the year, which may increase the sorption capacity of the streambed for other metals. Finally, we assessed fine-scale geochemical and hydrologic heterogeneity in the meandering riverbed by sampling pore water at 20 cm depth across more than 100 locations in August 2018. Sample locations span distinct zones of up-welling groundwater that contrast against dominant down-welling conditions. Clear differences in carbon quality and concentrations of redox sensitive solutes were detected between these zones, highlighting the importance of understanding heterogeneity at small scales across streambeds. Geochemical differences were associated with differences in microbial assemblage, but microbial diversity was generally uniform and high across all locations at the same depth, regardless of the extent of groundwater influence.
An analytical coupled ice-ocean model is considered which is forced by a specified wind stress acting on the open ocean as well as the ice. The analysis supports the conjecture that the upwelling dynamics at ice edges can be understood by means of a simple analytical model. In similarity with coastal problems it is shown that the ice edge upwelling is determined by the net mass flux at the boundaries of the considered region. The model is used to study the sensitivity of the upwelling dynamics in the marginal ice zone to variation in the controlling parameters. These parameters consist of combinations of the drag coefficients used in the parameterization of the stresses on the three interfaces atmosphere-ice, atmosphere-ocean, and ice-ocean. The response is shown to be sensitive to variations in these parameters in that one set of parameters may give upwelling while a slightly different set of parameters may give downwelling.