An Examination of Upper Ocean Processes in the Greenland Sea Utilizing Archived Fluxes, Satellite-Derived Ice Data, and a 1-D Upper Ocean Model
Explore the source record for details and available documents.
SEARCH · Search NASA
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.
This is a one-year cost extension of previous grant but carrying a new award number for the administrative purpose. Supported by this one-year extension, the following research has continued and obtained significant results. 20 papers have been published (9) or submitted (11) to scientific journals in this one-year period. A brief summary of scientific results on: 1. A new method for estimation of the sensible heat flux using satellite vector winds, 2. Pacific warm pool excitation, earth rotation and El Nino Southern Oscillations, 3. A new study of the Mediterranean outflow and Meddies at 400-meter isopycnal surface using multi-sensor data, 4. Response of the coastal ocean to extremely high wind, and 5. Role of wind on the estimation of heat flux using satellite data, are provided below as examples of our many research results conducted in the last year,
The relative importance of ocean and atmospheric dynamics in generating Atlantic Multidecadal Variability (AMV) remains an open question. Comparisons between climate models with SLAB and fully-dynamic (FULL) ocean components are often used to explore this question, but cannot reveal how individual ocean processes generate these differences. We build a hierarchy of physically interpretable stochastic models to investigate the contribution of two upper-ocean processes to AMV: the role of seasonal variation and mixed-layer entrainment. This interpretability arises from the stochastic model’s simplified representation of sea surface temperature (SST), considering only the local upper ocean response to white-noise atmospheric forcing and its impact on surface heat exchange. We focus on understanding differences between SLAB and FULL non-eddy resolving pre-industrial control simulations of the Community Earth System Model 1 (CESM), and estimate the stochastic model parameters from each respective simulation. Despite its simplicity, the stochastic model reproduces temporal characteristics of SST variability in the SPG, including reemergence, seasonal-to-interannual persistence and power spectra. Furthermore, unrealistically persistent SST of the CESM-SLAB ocean simulation is reproduced in the equivalent stochastic model configuration where the mixed-layer depth (MLD) is constant. The stochastic model also reveals that vertical entrainment primarily damps SST variability, thus explaining why SLAB exhibits larger SST variance than FULL. Here, the stochastic model driven by temporally stochastic, spatially coherent forcing patterns reproduces the canonical AMV pattern. However, the amplitude of low-frequency variability remains underestimated, suggesting a role for ocean dynamics beyond entrainment.
We calculate rates of ice thickness change and bottom melt for ice shelves in West Antarctica and the Antarctic Peninsula from a combination of elevation measurements from NASA–CECS Antarctic ice mapping campaigns and NASA Operation IceBridge corrected for oceanic processes from measurements and models, surface velocity measurements from synthetic aperture radar, and high-resolution outputs from regional climate models. The ice thickness change rates are calculated in a Lagrangian reference frame to reduce the effects from advection of sharp vertical features, such as cracks and crevasses, that can saturate Eulerian-derived estimates. We use our method over different ice shelves in Antarctica, which vary in terms of size, repeat coverage from airborne altimetry, and dominant processes governing their recent changes. We find that the Larsen-C Ice Shelf is close to steady state over our observation period with spatial variations in ice thickness largely due to the flux divergence of the shelf. Firn and surface processes are responsible for some short-term variability in ice thickness of the Larsen-C Ice Shelf over the time period. The Wilkins Ice Shelf is sensitive to short-timescale coastal and upper-ocean processes, and basal melt is the dominant contributor to the ice thickness change over the period. At the Pine Island Ice Shelf in the critical region near the grounding zone, we find that ice shelf thickness change rates exceed 40m/yr, with the change dominated by strong submarine melting. Regions near the grounding zones of the Dotson and Crosson ice shelves are decreasing in thickness at rates greater than 40m/yr, also due to intense basal melt. NASA–CECS Antarctic ice mapping and NASA Operation IceBridge campaigns provide validation datasets for floating ice shelves at moderately high resolution when coregistered using Lagrangian methods.
Air-sea turbulent fluxes determine the exchange of momentum, heat, freshwater, and gas between the atmosphere and ocean. These exchange processes are critical to a broad range of research questions spanning length scales from meters to thousands of kilometers and time scales from hours to decades. Examples are discussed (section 2). The estimation of surface turbulent fluxes from satellite is challenging and fraught with considerable errors (section 3); however, recent developments in retrievals (section 3) will greatly reduce these errors. Goals for the future observing system are summarized in section 4. Surface fluxes are defined as the rate per unit area at which something (e.g., momentum, energy, moisture, or CO Z ) is transferred across the air/sea interface. Wind- and buoyancy-driven surface fluxes are called surface turbulent fluxes because the mixing and transport are due to turbulence. Examples of nonturbulent processes are radiative fluxes (e.g., solar radiation) and precipitation (Schmitt et al., 2010). Turbulent fluxes are strongly dependent on wind speed; therefore, observations of wind speed are critical for the calculation of all turbulent surface fluxes. Wind stress, the vertical transport of horizontal momentum, also depends on wind direction. Stress is very important for many ocean processes, including upper ocean currents (Dohan and Maximenko, 2010) and deep ocean currents (Lee et al., 2010). On short time scales, this horizontal transport is usually small compared to surface fluxes. For long-term processes, transport can be very important but again is usually small compared to surface fluxes.
Previous case studies have linked cyclone-induced atmospheric forcing and/or upper-ocean processes to notable Arctic sea ice loss events in the summers of 2012 and 2016. This study examines a more recent and noteworthy case in late summer 2021 in which substantial sea ice loss followed a period of surface meteorological extremes in the Beaufort Sea region of the Arctic. We focus on the period from mid-August to mid-September 2021 that coincided with the Office of Naval Research THINICE Pilot Field Campaign and investigate stormy and windy conditions with respect to air-sea processes impacting sea ice conditions. We find that during the stormy first half of the campaign, cyclone-induced energy fluxes into the marginal ice zone and surrounding waters preconditioned the ice pack for more rapid melt later in the campaign. The second half of the campaign, in contrast, was marked by non-cyclone wind events that enhanced turbulent (namely sensible) heat fluxes into the ice and upper ocean that increased melt. Moreover, this latter period had enhanced advection of the Beaufort Sea ice pack into above-freezing waters, increasing bottom melt to >1 cm d −1 over the remainder of the campaign. While findings are shown to vary by surface type and at relatively small (i.e. ice-floe) scales, insights are offered on the roles of late summer coupled processes on rapid ice loss events in today’s Arctic environment.
Abstract We develop a hierarchy of simplified ocean models for coupled ocean, atmosphere, and sea ice climate simulations using the Community Earth System Model version 1 (CESM1). The hierarchy has four members: a slab ocean model, a mixed‐layer model (MLM) with entrainment and detrainment, an Ekman MLM, and an ocean general circulation model (OGCM). Flux corrections of heat and salt are applied to the simplified models ensuring that all hierarchy members have the same climatology. We diagnose the needed flux corrections from auxiliary simulations in which we restore the temperature and salinity to the daily climatology obtained from a target CESM1 simulation. The resulting three‐dimensional corrections contain the interannual variability fluxes that maintain the correct vertical gradients of temperature and salinity in the tropics. We find that the inclusion of mixed‐layer entrainment and Ekman flow produces sea surface temperature and surface air temperature fields whose means and variances are progressively more similar to those produced by the target CESM1 simulation. We illustrate the application of the hierarchy to the problem of understanding the response of the climate system to the loss of Arctic sea ice. We find that the shifts in the positions of the mid‐latitude westerly jet and of the Inter‐tropical Convergence Zone (ITCZ) in response to sea‐ice loss depend critically on upper ocean processes. Specifically, heat uptake associated with the mixed‐layer entrainment influences the shift in the westerly jet and ITCZ. Moreover, the shift of ITCZ is sensitive to the form of Ekman flow parameterization.
We develop a hierarchy of simplified ocean models for coupled ocean, atmosphere, and sea ice climate simulations using the Community Earth System Model version 1 (CESM1). The hierarchy has four members: a slab ocean model, a mixed-layer model (MLM) with entrainment and detrainment, an Ekman MLM, and an ocean general circulation model (OGCM). Flux corrections of heat and salt are applied to the simplified models ensuring that all hierarchy members have the same climatology. We diagnose the needed flux corrections from auxiliary simulations in which we restore the temperature and salinity to the daily climatology obtained from a target CESM1 simulation. The resulting three-dimensional corrections contain the interannual variability fluxes that maintain the correct vertical gradients of temperature and salinity in the tropics. We find that the inclusion of mixed-layer entrainment and Ekman flow produces sea surface temperature and surface air temperature fields whose means and variances are progressively more similar to those produced by the target CESM1 simulation. We illustrate the application of the hierarchy to the problem of understanding the response of the climate system to the loss of Arctic sea ice. We find that the shifts in the positions of the mid-latitude westerly jet and of the Inter-tropical Convergence Zone (ITCZ) in response to sea-ice loss depend critically on upper ocean processes. Specifically, heat uptake associated with the mixed-layer entrainment influences the shift in the westerly jet and ITCZ. Moreover, the shift of ITCZ is sensitive to the form of Ekman flow parameterization.
The Seasat satellite launched in Summer 1978 carried a synthetic aperture radar (SAR). Although Seasat failed after 105 days in orbit, it provided observations that demonstrate the potential to examine and monitor upper oceanic processes. Seasat made five passes over lower Cook Inlet and Shelikof Strait, Alaska, during Summer 1978. SAR images from the passes show oceanographic features, including a meander in a front, a pair of mesoscale eddies, and internal waves. These features are compared with contemporary and representative images from a satellite-borne Advanced Very High Resolution Radiometer (AVHRR) and Coastal Zone Color Scanner (CZCS), with water property data, and with current observations from moored instruments. The results indicate that SAR data can be used to monitor mesoscale oceanographic features.
The Harmony satellite mission was recently approved as the next European Space Agency (ESA) Earth Explorer 10. The mission science objectives cover several applications related to solid earth, the cryosphere, upper-ocean dynamics and air–sea interactions. The mission consists of a constellation of two satellites, flying with the Copernicus Sentinel 1 (C or D) spacecraft, each hosting a C-band receive-only radar and a thermal infrared (TIR) payload. From an ocean dynamics/air–sea interaction perspective, the mission will provide the unique opportunity to observe simultaneously the signature of submesoscale upper-ocean processes via synthetic aperture radar and TIR imagery. The TIR imager is based on microbolometer technology and its acquisitions will rely on four channels: three narrow-band channels yielding observations at a ≃1 km spatial sampling distance (SSD) and a panchromatic (PAN, 8–12 μm) channel characterized by a ≃300 m SSD. Our study investigates the potential of Harmony in retrieving spatial features related to sea surface temperature (SST) gradients from the high-resolution PAN channel, relying on top-of-atmosphere (TOA) observations. Compared to a standard SST gradient retrieval, our approach does not require atmospheric correction, thus avoiding uncertainties due to inter-channel co-registration and radiometric consistency, with the possibility of exploiting the higher resolution of the PAN channel. The investigations were carried out simulating the future Harmony TOA radiances (TARs), as well as relying on existing state-of-the-art level 1 satellite products. Our approach enables the correct description of SST features at the sea surface avoiding the generation of spurious features due to atmospheric correction and/or instrumental issues. In addition, analyses based on existing satellite products suggest that the clear-sky TOA observations, in a typical mid-latitude scene, allow the reconstruction of up to 85% of the gradient magnitudes found at the sea-surface level. The methodology is less efficient in tropical areas, suffering from smoothing effects due to the high concentrations of water vapor.
Considerable uncertainty remains as to the importance of mechanisms governing decadal and longer variability of the Walker Circulation, its connection to the tropical climate system, and prospects for tropical climate change in the face of anthropogenic forcing. Most contemporary climate models suggest that in response to elevated CO2 and a warmer but more stratified atmosphere, the required upward mass flux in tropical convection will diminish along with the Walker component of the tropical mean circulation as well. Alternatively, there is also evidence to suggest that the shoaling and increased vertical stratification of the thermocline in the eastern Pacific will enable a muted SST increase there-- preserving or even enhancing some of the dynamical forcing for the Walker cell flow. Over the past decade there have been observational indications of an acceleration in near-surface easterlies, a strengthened Pacific zonal SST gradient, and globally-teleconnected dislocations in precipitation. But is this evidence in support of an ocean dynamical thermostat process posited to accompany anthropogenic forcing, or just residual decadal fluctuations associated with variations in warm and cold ENSO events and other stochastic forcing? From a modeling perspective we try to make headway on this question by examining zonal variations in surface energy fluxes and dynamics governing tropical upper ocean heat content evolution in the WCRP CMIP5 model projections. There is some diversity among model simulations; for example, the CCSM4 indicates net ocean warming over the IndoPacific region while the CSIRO model concentrates separate warming responses over the central Pacific and Indian Ocean regions. The models, as with observations, demonstrate strong local coupling between variations in column water vapor, downward surface longwave radiation and SST; but the spatial patterns of changes in the sign of this relationship differ among models and, for models as a whole, with observations. Our analysis focuses initially on probing the inter-model differences in energy fluxes / transports and Walker Circulation response to forcing. We then attempt to identify statistically the El Nino- / La Nina-related ocean heat content variability unique to each model and regress out the associated energy flux, ocean heat transport and Walker response on these shorter time scales for comparison to that of the anthropogenic signals.
In this talk, I will present a brief overview of some recent results and current research priorities regarding the significant participation of upper ocean / submesoscale processes in regional to global Earth system dynamics in a warming climate, such as ocean overturning, biogeochemistry, and heat and carbon uptake as well as atmospheric, coastal, and cryospheric dynamics. In addition to reviewing the evidence for these links by topic, I will highlight how advances in both understanding and prediction are derived from the integrated triad of process studies, modeling, and regional and global observing systems. Finally, I will discuss how triadic investments focused on these upper ocean processes might lead to better-informed decision-making for climate change adaptation and mitigation strategies.
Coupling a cloud resolving model (CRM) with an ocean mixed layer (OML) model can provide a powerful tool for better understanding impacts of atmospheric precipitation on sea surface temperature (SST) and salinity. The objective of this study is twofold. First, by using the three dimensional (3-D) CRM-simulated (the Goddard Cumulus Ensemble model, GCE) diabatic source terms, radiation (longwave and shortwave), surface fluxes (sensible and latent heat, and wind stress), and precipitation as input for the OML model, the respective impact of individual component on upper ocean heat and salt budgets are investigated. Secondly, a two-way air-sea interaction between tropical atmospheric climates (involving atmospheric radiative-convective processes) and upper ocean boundary layer is also examined using a coupled two dimensional (2-D) GCE and OML model. Results presented here, however, only involve the first aspect. Complete results will be presented at the conference.
Earth system models parameterize ocean surface fluxes of heat, moisture, and momentum with empirical bulk flux algorithms, which introduce biases and uncertainties into simulations. We investigate the atmosphere and ocean model sensitivity to algorithm choice in the Energy Exascale Earth System Model (E3SM). Flux differences between algorithms are larger in atmosphere simulations (where wind speeds can vary) than ocean simulations (where wind speeds are fixed by forcing data). Surface flux changes lead to global scale changes in the energy and water cycles, notably including ocean heat uptake and global mean precipitation rates. Compared to the control algorithm, both COARE and University of Arizona (UA) algorithms reduce global mean precipitation and top of atmosphere radiative biases. Further, UA may slightly reduce biases in ocean meridional heat transport. We speculate that changes seen here, especially in the ocean, could be even larger in coupled simulations.
The abrupt warming of the north central Pacific Ocean from 1997 to 1999 is studied using an ocean data assimilation product. During this period, the average mixed-layer temperature in the region of 170-210(deg)E, 25-40(deg)N rises by 1.8 K. The major contributors to the warming are surface heat flux (1.3 K), geostrophic advection (0.7 K), and entrainment (0.7 K). For the geostrophic advection, the contributions by the zonal, meridional, and vertical components are 0.4, -0.1 and 0.3 K, respectively. Mixing and meridional Ekman advection have cooling effect. The significance of the geostrophic advection indicates the importance of ocean dynamics in controlling the abrupt warming tendency during the 1997-99 period and the inadequacy of a slab-mixed-layer model in simulating such warming tendency.
The cosmogenic radionuclide 7[superscript]Be has been applied as a tracer of dynamical processes in the upper ocean and of atmospheric deposition of trace elements at the sea surface. These applications usually assume that 7[superscript]Be is entirely in the dissolved form, and that scavenging and downward export of 7[superscript]Be by settling particles can be neglected. In this work, we explore these assumptions and more generally assess the significance of the 7[superscript]Be activity in the particulate fraction, through the generation of vertical profiles of particulate 7[superscript]Be in the open ocean. From detailed measurements obtained from low-background gamma spectrometers placed in underground facilities, we report vertical profiles of 7[superscript]Be activity in suspended particles (7[superscript]Bep[subscript]) collected in various oceanic regions: the Mediterranean Sea (DYFAMED station), the Indian Sector of the Southern Ocean (station A3-2 from the KEOPS2 cruise), the Sargasso Sea (OFP station), and the subpolar North Atlantic Ocean (GEOVIDE cruise). We find that, in each oceanic region, 7[superscript]Bep[subscript] activities are generally higher in the mixed layer than in the thermocline. They vary in the mixed layer from 3.0 dpm/m3 at DYFAMED to 33.1dpm/m3 at GEOVIDE cross-over station 51/60, i.e., within a range consistent with previous 7[superscript]Bep[subscript] measurements for the open ocean. For the GEOVIDE cruise, the 7[superscript]Bep[subscript] activities measured on different filter types at different depths are corrected for filter offsets derived from multiple 7[superscript]Bep[subscript] measurements at a near-coastal station in the western Mediterranean Sea. We then combine measurements of total 7[superscript]Be activity (Shelley et al., 2017) with our measurements of 7[superscript]Bep[subscript] activity to estimate the solid/solution partitioning of 7[superscript]Be. On average, the particulate fraction would represent 6% of total 7[superscript]Be activity at 5-m water depth (n = 6), 22% at 20 m (n = 2) and 9% at 70 m (n = 3). At GEOVIDE stations, 7[superscript]Bep[subscript] inventories range from 5% to 19% of the total 7[superscript]Be inventories. In the Labrador Sea, the measured 7[superscript]Bep[subscript] inventories are lower than the dry 7[superscript]Be deposition fluxes estimated from aerosol samples collected during GEOVIDE, suggesting that a significant portion of 7[superscript]Bep[subscript] may be removed by sinking particles. The distribution coefficient Kd for 7Be increases with depth, with log10Kd averaging 5.2 ± 0.1 at 5 m to 6.1 ± 0.1 between 70 and 150 m, suggesting that scavenging of dissolved 7[superscript]Be by particles is more pronounced in the thermocline than in the mixed layer when differences in particle concentrations are taken into account. Overall, our study suggests that, at least in some oceanic regions, the removal of 7[superscript]Be by marine particles may be significant and that it may need to be considered in applications of 7[superscript]Be as a tracer of oceanic processes and atmospheric deposition.
A method has been developed to determine the oceanic daily mean mixed layer depth from satellite observations and a mixed layer thermal inertia (MLTI) model. The algorithms were developed to use remotely-sensed values of sea surface temperature, albedo, and surface wind speeds to calculate the thermal inertia and to predict changes in subsurface diurnal mixed layer depth. The MLTI model, based on a mixed layer model of the upper ocean, has been used to simulate the diurnal mixing process and thermal inertia distribution in the Sargasso Sea around 34 deg N, 70 deg W. Sea surface temperature and albedo have been obtained from the NOAA7-AVHRR images. Surface wind speeds have been derived from the Scanning Multichannel Microwave Radiometer (SMMR) aboard Nimbus 7. Image processing was performed for images gathered between June and July 1982. The daily mean mixed layer depths predicted by the MLTI model agree well with data gathered at the LOTUS mooring located in the Sargasso Sea. This suggests that vertical mixing is the dominant physical process that controls the thermal inertia distribution in the midocean, far from major current systems, and that remote sensing is a promising tool to study such upper ocean processes.
The climatological seasonal cycle of sea surface temperature (SST) in the tropical Pacific is simulated using a newly developed upper ocean model. The roles of vertical mixing, solar radiation, and wind stress are investigated in a hierarchy of numerical experiments with various combinations of vertical mixing algorithms and surface-forcing products. It is found that the large SST annual cycle in the eastern equatorial Pacific is, to a large extent, controlled by the annually varying mixed layer depth which, in turn, is mainly determined by the competing effects of solar radiation and wind forcing. With the application of our hybrid vertical mixing scheme the model-simulated SST annual cycle is much improved in both amplitude and phase as compared to the case of a constant mixed layer depth. Beside the strong effects on vertical mixing, solar radiation is the primary heating term in the surface layer heat budget, and wind forcing influences SST by driving oceanic advective processes that redistribute heat in the upper ocean. For example, the SST seasonal cycle in the western Pacific basically follows the semiannual variation of solar heating, and the cycle in the central equatorial region is significantly affected by the zonal advective heat flux associated with the seasonally reversing South Equatorial Current. It has been shown in our experiments that the amount of heat flux modification needed to eliminate the annual mean SST errors in the model is, on average, no larger than the annual mean uncertainties among the various surface flux products used in this study. Whereas a bias correction is needed to account for remaining uncertainties in the annual mean heat flux, this study demonstrates that with proper treatment of mixed layer physics and realistic forcing functions the seasonal variability of SST is capable of being simulated successfully in response to external forcing without relying on a relaxation or damping formulation for the dominant surface heat flux contributions.