Search Results
You are looking at 1 - 10 of 41 items for
- Author or Editor: Frank Bryan x
- Refine by Access: All Content x
Abstract
Experiments with a low resolution, primitive equation ocean general circulation model with idealized basin geometry and surface forcing have been carried out in order to identify the processes controlling the climatically important aspects of the circulation. Emphasis was placed on the sensitivity of the model solutions to the magnitude of the vertical diffusivity. Scaling arguments suggest, and the numerical experiments confirm, that the solutions are most sensitive to the magnitudes of the wind stress curl and the vertical diffusivity. For small vertical diffusivity, wind forcing dominates the solution. The vertical scale of the thermocline is set by the strength of the Ekman pumping, and there is a multiple gyre circulation in the upper layers. For large vertical diffusivity, diabatic surface forcing dominates the solution. Vertical diffusion controls the vertical scale of the thermocline, and there is a single large anticyclonic gyre in the upper layers. Both the meridionally and zonally integrated overturning circulations are sensitive to the vertical diffusivity, though not to the same degree. The poleward heat transport is dominated by the zonally integrated meridional overturning circulation and, hence, also shows a sensitivity to the vertical diffusivity. The maximum poleward heat transport for the model used in this study varies by an order of magnitude as the vertical diffusivity is varied over a range comparable to that of estimates based on observations. The sensitivity to the prescription of the surface diabatic forcing and to other closure parameters is also described.
Abstract
Experiments with a low resolution, primitive equation ocean general circulation model with idealized basin geometry and surface forcing have been carried out in order to identify the processes controlling the climatically important aspects of the circulation. Emphasis was placed on the sensitivity of the model solutions to the magnitude of the vertical diffusivity. Scaling arguments suggest, and the numerical experiments confirm, that the solutions are most sensitive to the magnitudes of the wind stress curl and the vertical diffusivity. For small vertical diffusivity, wind forcing dominates the solution. The vertical scale of the thermocline is set by the strength of the Ekman pumping, and there is a multiple gyre circulation in the upper layers. For large vertical diffusivity, diabatic surface forcing dominates the solution. Vertical diffusion controls the vertical scale of the thermocline, and there is a single large anticyclonic gyre in the upper layers. Both the meridionally and zonally integrated overturning circulations are sensitive to the vertical diffusivity, though not to the same degree. The poleward heat transport is dominated by the zonally integrated meridional overturning circulation and, hence, also shows a sensitivity to the vertical diffusivity. The maximum poleward heat transport for the model used in this study varies by an order of magnitude as the vertical diffusivity is varied over a range comparable to that of estimates based on observations. The sensitivity to the prescription of the surface diabatic forcing and to other closure parameters is also described.
Abstract
In this study, the salinity budget of the North Atlantic subtropical salinity maximum region for control volumes bounded by isohaline surfaces is analyzed. The authors provide closed budgets based on output from a high-resolution numerical simulation and partial budgets based on analyses of observational climatologies of hydrography and surface fluxes. With this choice of control volume, advection is eliminated from the instantaneous volume-integrated salt budget, and time-mean advection is eliminated from the budget evaluated from time-averaged data. In this way, the role of irreversible mixing processes in the maintenance and variability of the salinity maximum are more readily revealed. By carrying out the analysis with both near-instantaneous and time-averaged model output, the role of mesoscale eddies in stirring and mixing for this water mass is determined. This study finds that the small-scale mixing acting on enhanced gradients generated by the mesoscale eddies is approximately equal to that acting on the large-scale gradients estimated from climatological-mean conditions. The isohaline salinity budget can be related to water mass transformation rates associated with surface forcing and mixing processes in a straightforward manner. The authors find that the surface net evaporation in the North Atlantic salinity maximum region accounts for a transformation of 7 Sverdrups (Sv; 1 Sv ≡ 106 m3 s−1) of water across the 37-psu isohaline outcrop into the salinity maximum in the simulation, whereas the estimate based on climatological observations is 9 to 10 Sv.
Abstract
In this study, the salinity budget of the North Atlantic subtropical salinity maximum region for control volumes bounded by isohaline surfaces is analyzed. The authors provide closed budgets based on output from a high-resolution numerical simulation and partial budgets based on analyses of observational climatologies of hydrography and surface fluxes. With this choice of control volume, advection is eliminated from the instantaneous volume-integrated salt budget, and time-mean advection is eliminated from the budget evaluated from time-averaged data. In this way, the role of irreversible mixing processes in the maintenance and variability of the salinity maximum are more readily revealed. By carrying out the analysis with both near-instantaneous and time-averaged model output, the role of mesoscale eddies in stirring and mixing for this water mass is determined. This study finds that the small-scale mixing acting on enhanced gradients generated by the mesoscale eddies is approximately equal to that acting on the large-scale gradients estimated from climatological-mean conditions. The isohaline salinity budget can be related to water mass transformation rates associated with surface forcing and mixing processes in a straightforward manner. The authors find that the surface net evaporation in the North Atlantic salinity maximum region accounts for a transformation of 7 Sverdrups (Sv; 1 Sv ≡ 106 m3 s−1) of water across the 37-psu isohaline outcrop into the salinity maximum in the simulation, whereas the estimate based on climatological observations is 9 to 10 Sv.
Abstract
The National Center for Atmospheric Research’s Climate System Model is a comprehensive model of the physical climate system. A 300-yr integration of the model has been carried out without flux correction. The solution shows very little drift in the surface temperature distribution, sea-ice extent, or atmospheric circulation. The lack of drift in the surface climate is attributed to relatively good agreement in the estimates of meridional heat transport in the uncoupled ocean model and that implied by the uncoupled atmospheric model. On the other hand, there is significant drift in the temperature and salinity distributions of the deep ocean. The ocean loses heat at an area-averaged rate of 0.35 W m−2, the upper ocean becomes fresher, and the deep ocean becomes colder and saltier than in the uncoupled ocean model equilibrium or in observations. The cause of this drift is an unreasonably large meridional transport of freshwater in the sea ice model, resulting in the production of excessively cold and salty Antarctic Bottom Water. There is also significant drift in the Arctic basin, with the complete erosion of the surface halocline early in the coupled integration.
Abstract
The National Center for Atmospheric Research’s Climate System Model is a comprehensive model of the physical climate system. A 300-yr integration of the model has been carried out without flux correction. The solution shows very little drift in the surface temperature distribution, sea-ice extent, or atmospheric circulation. The lack of drift in the surface climate is attributed to relatively good agreement in the estimates of meridional heat transport in the uncoupled ocean model and that implied by the uncoupled atmospheric model. On the other hand, there is significant drift in the temperature and salinity distributions of the deep ocean. The ocean loses heat at an area-averaged rate of 0.35 W m−2, the upper ocean becomes fresher, and the deep ocean becomes colder and saltier than in the uncoupled ocean model equilibrium or in observations. The cause of this drift is an unreasonably large meridional transport of freshwater in the sea ice model, resulting in the production of excessively cold and salty Antarctic Bottom Water. There is also significant drift in the Arctic basin, with the complete erosion of the surface halocline early in the coupled integration.
Abstract
For the first time estimates of divergent eddy heat flux (DEHF) from a high-resolution (0.1°) simulation of the Parallel Ocean Program (POP) are compared with estimates made during the Kuroshio Extension System Study (KESS). The results from POP are in good agreement with KESS observations. POP captures the lateral and vertical structure of mean-to-eddy energy conversion rates, which range from 2 to 10 cm2 s−3. The dynamical mechanism of vertical coupling between the deep and upper ocean is the process responsible for DEHFs in POP and is in accordance with baroclinic instability observed in the Gulf Stream and Kuroshio Extension. Meridional eddy heat transport values are ~14% larger in POP at its maximum value. This is likely due to the more zonal path configuration in POP. The results from this study suggest that HR POP is a useful tool for estimating eddy statistics in the Kuroshio Extension region, and thereby provide guidance in the formulation and testing of eddy mixing parameterization schemes.
Abstract
For the first time estimates of divergent eddy heat flux (DEHF) from a high-resolution (0.1°) simulation of the Parallel Ocean Program (POP) are compared with estimates made during the Kuroshio Extension System Study (KESS). The results from POP are in good agreement with KESS observations. POP captures the lateral and vertical structure of mean-to-eddy energy conversion rates, which range from 2 to 10 cm2 s−3. The dynamical mechanism of vertical coupling between the deep and upper ocean is the process responsible for DEHFs in POP and is in accordance with baroclinic instability observed in the Gulf Stream and Kuroshio Extension. Meridional eddy heat transport values are ~14% larger in POP at its maximum value. This is likely due to the more zonal path configuration in POP. The results from this study suggest that HR POP is a useful tool for estimating eddy statistics in the Kuroshio Extension region, and thereby provide guidance in the formulation and testing of eddy mixing parameterization schemes.
Abstract
We use an interannually forced version of the Parallel Ocean Program, configured to resolve mesoscale eddies, to close the global eddy potential energy (EPE) budget associated with temperature variability. By closing the EPE budget, we are able to properly investigate the role of diabatic processes in modulating mesoscale energetics in the context of other processes driving eddy–mean flow interactions. A Helmholtz decomposition of the eddy heat flux field into divergent and rotational components is applied to estimate the baroclinic conversion from mean to eddy potential energy. In doing so, an approximate two-way balance between the “divergent” baroclinic conversion and upgradient vertical eddy heat fluxes in the ocean interior is revealed, in accordance with baroclinic instability and the relaxation of isopycnal slopes. However, in the mixed layer, the EPE budget is greatly modulated by diabatic mixing, with air–sea interactions and interior diffusion playing comparable roles. Globally, this accounts for ∼60% of EPE converted to EKE (eddy kinetic energy), with the remainder being dissipated by air–sea interactions and interior mixing. A seasonal composite of baroclinic energy conversions shows that the strongest EPE to EKE conversion occurs during the summer in both hemispheres. The seasonally varying diabatic processes in the upper ocean are further shown to be closely linked to this EPE–EKE conversion seasonality, but with a lead. The peak energy dissipation through vertical mixing occurs ahead of the minimum EKE generation by 1–2 months.
Abstract
We use an interannually forced version of the Parallel Ocean Program, configured to resolve mesoscale eddies, to close the global eddy potential energy (EPE) budget associated with temperature variability. By closing the EPE budget, we are able to properly investigate the role of diabatic processes in modulating mesoscale energetics in the context of other processes driving eddy–mean flow interactions. A Helmholtz decomposition of the eddy heat flux field into divergent and rotational components is applied to estimate the baroclinic conversion from mean to eddy potential energy. In doing so, an approximate two-way balance between the “divergent” baroclinic conversion and upgradient vertical eddy heat fluxes in the ocean interior is revealed, in accordance with baroclinic instability and the relaxation of isopycnal slopes. However, in the mixed layer, the EPE budget is greatly modulated by diabatic mixing, with air–sea interactions and interior diffusion playing comparable roles. Globally, this accounts for ∼60% of EPE converted to EKE (eddy kinetic energy), with the remainder being dissipated by air–sea interactions and interior mixing. A seasonal composite of baroclinic energy conversions shows that the strongest EPE to EKE conversion occurs during the summer in both hemispheres. The seasonally varying diabatic processes in the upper ocean are further shown to be closely linked to this EPE–EKE conversion seasonality, but with a lead. The peak energy dissipation through vertical mixing occurs ahead of the minimum EKE generation by 1–2 months.
Abstract
Observational and model evidence has been mounting that mesoscale eddies play an important role in air–sea interaction in the vicinity of western boundary currents and can affect the jet stream storm track. What is less clear is the interplay between oceanic and atmospheric meridional heat transport in the vicinity of western boundary currents. It is first shown that variability in the North Pacific, particularly in the Kuroshio Extension region, simulated by a high-resolution fully coupled version of the Community Earth System Model matches observations with similar mechanisms and phase relationships involved in the variability. The Pacific decadal oscillation (PDO) is correlated with sea surface height anomalies generated in the central Pacific that propagate west preceding Kuroshio Extension variability with a ~3–4-yr lag. It is then shown that there is a near compensation of O(0.1) PW (PW ≡ 1015 W) between wintertime atmospheric and oceanic meridional heat transport on decadal time scales in the North Pacific. This compensation has characteristics of Bjerknes compensation and is tied to the mesoscale eddy activity in the Kuroshio Extension region.
Abstract
Observational and model evidence has been mounting that mesoscale eddies play an important role in air–sea interaction in the vicinity of western boundary currents and can affect the jet stream storm track. What is less clear is the interplay between oceanic and atmospheric meridional heat transport in the vicinity of western boundary currents. It is first shown that variability in the North Pacific, particularly in the Kuroshio Extension region, simulated by a high-resolution fully coupled version of the Community Earth System Model matches observations with similar mechanisms and phase relationships involved in the variability. The Pacific decadal oscillation (PDO) is correlated with sea surface height anomalies generated in the central Pacific that propagate west preceding Kuroshio Extension variability with a ~3–4-yr lag. It is then shown that there is a near compensation of O(0.1) PW (PW ≡ 1015 W) between wintertime atmospheric and oceanic meridional heat transport on decadal time scales in the North Pacific. This compensation has characteristics of Bjerknes compensation and is tied to the mesoscale eddy activity in the Kuroshio Extension region.
Abstract
The water mass transformation (WMT) framework describes how water of one class, such as a discrete interval of density, is converted into another class via air–sea fluxes or interior mixing processes. This paper investigates how this process is modified at the surface when mesoscale ocean eddies are present, using a state-of-the-art high-resolution climate model with reasonable fidelity in the Southern Ocean. The method employed is to coarse-grain the high-resolution model fields to remove eddy signatures, and compare the results with those from the full model fields. This method shows that eddies reduced the WMT by 2–4 Sv (10%–20%; 1 Sv ≡ 106 m3 s−1) over a wide range of densities, from typical values of 20 Sv in the smoothed case. The corresponding water mass formation was reduced by 40% at one particular density increment, namely, between 1026.4 and 1026.5 kg m−3, which corresponds to the lighter end of the range of Indian Ocean Mode Water in the model. The effect of eddies on surface WMT is decomposed into three terms: direct modulation of the density outcrops, then indirectly, by modifying the air–sea density flux, and the combined effect of the two, akin to a covariance. It is found that the first and third terms dominate, i.e., smoothing the outcrops alone has a significant effect, as does the combination of smoothing both outcrops and density flux distributions, but smoothing density flux fields alone has little effect. Results from the coarse-graining method are compared to an alternative approach of temporally averaging the data. Implications for climate model resolution are also discussed.
Abstract
The water mass transformation (WMT) framework describes how water of one class, such as a discrete interval of density, is converted into another class via air–sea fluxes or interior mixing processes. This paper investigates how this process is modified at the surface when mesoscale ocean eddies are present, using a state-of-the-art high-resolution climate model with reasonable fidelity in the Southern Ocean. The method employed is to coarse-grain the high-resolution model fields to remove eddy signatures, and compare the results with those from the full model fields. This method shows that eddies reduced the WMT by 2–4 Sv (10%–20%; 1 Sv ≡ 106 m3 s−1) over a wide range of densities, from typical values of 20 Sv in the smoothed case. The corresponding water mass formation was reduced by 40% at one particular density increment, namely, between 1026.4 and 1026.5 kg m−3, which corresponds to the lighter end of the range of Indian Ocean Mode Water in the model. The effect of eddies on surface WMT is decomposed into three terms: direct modulation of the density outcrops, then indirectly, by modifying the air–sea density flux, and the combined effect of the two, akin to a covariance. It is found that the first and third terms dominate, i.e., smoothing the outcrops alone has a significant effect, as does the combination of smoothing both outcrops and density flux distributions, but smoothing density flux fields alone has little effect. Results from the coarse-graining method are compared to an alternative approach of temporally averaging the data. Implications for climate model resolution are also discussed.
Abstract
Present-day control and 1% yr−1 increasing carbon dioxide runs have been made using two versions of the Community Climate System Model, version 3.5. One uses the standard versions of the ocean and sea ice components where the horizontal resolution is 1° and the effects of mesoscale eddies are parameterized, and the second uses a resolution of 1/10° where the eddies are resolved. This is the first time the parameterization has been tested in a climate change run compared to an eddy-resolving run. The comparison is made not straightforward by the fact that the two control run climates are not the same, especially in their sea ice distributions. The focus is on the Antarctic Circumpolar Current region, where the effects of eddies are of leading order. The conclusions are that many of the differences in the two carbon dioxide transient forcing runs can be explained by the different control run sea ice distributions around Antarctica, but there are some quantitative differences in the meridional overturning circulation, poleward heat transport, and zonally averaged heat uptake when the eddies are parameterized rather than resolved.
Abstract
Present-day control and 1% yr−1 increasing carbon dioxide runs have been made using two versions of the Community Climate System Model, version 3.5. One uses the standard versions of the ocean and sea ice components where the horizontal resolution is 1° and the effects of mesoscale eddies are parameterized, and the second uses a resolution of 1/10° where the eddies are resolved. This is the first time the parameterization has been tested in a climate change run compared to an eddy-resolving run. The comparison is made not straightforward by the fact that the two control run climates are not the same, especially in their sea ice distributions. The focus is on the Antarctic Circumpolar Current region, where the effects of eddies are of leading order. The conclusions are that many of the differences in the two carbon dioxide transient forcing runs can be explained by the different control run sea ice distributions around Antarctica, but there are some quantitative differences in the meridional overturning circulation, poleward heat transport, and zonally averaged heat uptake when the eddies are parameterized rather than resolved.
Abstract
The global distributions of the air–sea fluxes of heat and freshwater and water mass transformation rates from a control integration of the coupled National Center for Atmospheric Research (NCAR) Climate System Model (CSM) are compared with similar fields from an uncoupled ocean model equilibrium spinup and a new surface climatology. The climatology and uncoupled model use the same bulk-flux forcing scheme and are forced with National Centers for Environmental Predicition (formerly the National Meteorological Center) atmospheric reanalysis data and satellite-based cloud cover, solar flux, and precipitation estimates. The climatological fluxes for the open ocean are adjusted to give a global net balance and are in broad general agreement with standard ship-based estimates. An exception is the ice-free Southern Ocean, where the net heat and evaporative fluxes appear to be too weak but where the observational coverage underlying the reanalyis is quite poor. Major differences are observed between the climatology and the NCAR CSM coupled solution, namely, enhanced tropical and subtropic solar insolation, stronger energy and hydrologic cycles, and excessive high-latitude ice formation/melt producing a several-fold increase in Arctic and Antarctic deep water formation through brine rejection. The anomalous fluxes and corresponding water-mass transformations are closely tied to the coupled ocean model drift, characterized by a reorganization of the vertical salinity distribution. Some error features in the heat flux and sea surface temperature fields are common to both the coupled and uncoupled solutions, primarily in the western boundary currents and the Antarctic circumpolar current, and are thus likely due to the poor representation of the circulation field in the coarse-resolution NCAR ocean model. Other problems particular to the uncoupled spinup are related to the bulk-flux forcing scheme, an example being excess freshwater deposition in the western boundary currents arising from the inclusion of a weak open ocean surface salinity restoring term. The effective thermal restoring coefficent, which relates the change in nonsolar surface heat flux to sea surface temperature changes, is on average 14.6 W m−2 K−1 for the coupled solution or about a third of the range from the bulk flux forcing scheme, 40–60 W m−2 K−1.
Abstract
The global distributions of the air–sea fluxes of heat and freshwater and water mass transformation rates from a control integration of the coupled National Center for Atmospheric Research (NCAR) Climate System Model (CSM) are compared with similar fields from an uncoupled ocean model equilibrium spinup and a new surface climatology. The climatology and uncoupled model use the same bulk-flux forcing scheme and are forced with National Centers for Environmental Predicition (formerly the National Meteorological Center) atmospheric reanalysis data and satellite-based cloud cover, solar flux, and precipitation estimates. The climatological fluxes for the open ocean are adjusted to give a global net balance and are in broad general agreement with standard ship-based estimates. An exception is the ice-free Southern Ocean, where the net heat and evaporative fluxes appear to be too weak but where the observational coverage underlying the reanalyis is quite poor. Major differences are observed between the climatology and the NCAR CSM coupled solution, namely, enhanced tropical and subtropic solar insolation, stronger energy and hydrologic cycles, and excessive high-latitude ice formation/melt producing a several-fold increase in Arctic and Antarctic deep water formation through brine rejection. The anomalous fluxes and corresponding water-mass transformations are closely tied to the coupled ocean model drift, characterized by a reorganization of the vertical salinity distribution. Some error features in the heat flux and sea surface temperature fields are common to both the coupled and uncoupled solutions, primarily in the western boundary currents and the Antarctic circumpolar current, and are thus likely due to the poor representation of the circulation field in the coarse-resolution NCAR ocean model. Other problems particular to the uncoupled spinup are related to the bulk-flux forcing scheme, an example being excess freshwater deposition in the western boundary currents arising from the inclusion of a weak open ocean surface salinity restoring term. The effective thermal restoring coefficent, which relates the change in nonsolar surface heat flux to sea surface temperature changes, is on average 14.6 W m−2 K−1 for the coupled solution or about a third of the range from the bulk flux forcing scheme, 40–60 W m−2 K−1.
Abstract
The National Center for Atmospheric Research (NCAR) Ocean Model has been developed for use in NCAR’s Climate System Modeling project, a comprehensive development of a coupled ocean–atmosphere–sea ice–land surface model of the global climate system. As part of this development, new parameterizations of diffusive mixing by unresolved processes have been implemented for the tracer equations in the model. Because the strength of the mixing depends upon the density structure under these parameterizations, it is possible that local explicit mixing may be quite small in selected locations, in contrast to the constant diffusivity model generally used. When a spatially centered advection scheme is used in the standard model configuration, local overshooting of tracer values occurs, leading to unphysical maxima and minima in the fields. While the immediate problem is a local Gibbs’s phenomenon, there is the possibility that such local tracer anomalies might propagate by advection and diffusion far from the source, causing inaccuracies in the tracer fields globally.
Because of these issues, a third-order upwind scheme was implemented for the advection of tracers. Numerical experiments show that this scheme is computationally efficient compared to alternatives (such as the flux-corrected transport scheme) and that it works well with other aspects of the model, such as acceleration (important for spinup efficiency) and the new mixing parameterizations. The scheme mimimizes overshooting effects while keeping the dissipative aspect of the advective operator reasonably small. The net effect is to produce solutions in which the large-scale fields are affected very little while local extrema are nearly (but not completely) removed, leading to physically much more realistic tracer patterns.
Abstract
The National Center for Atmospheric Research (NCAR) Ocean Model has been developed for use in NCAR’s Climate System Modeling project, a comprehensive development of a coupled ocean–atmosphere–sea ice–land surface model of the global climate system. As part of this development, new parameterizations of diffusive mixing by unresolved processes have been implemented for the tracer equations in the model. Because the strength of the mixing depends upon the density structure under these parameterizations, it is possible that local explicit mixing may be quite small in selected locations, in contrast to the constant diffusivity model generally used. When a spatially centered advection scheme is used in the standard model configuration, local overshooting of tracer values occurs, leading to unphysical maxima and minima in the fields. While the immediate problem is a local Gibbs’s phenomenon, there is the possibility that such local tracer anomalies might propagate by advection and diffusion far from the source, causing inaccuracies in the tracer fields globally.
Because of these issues, a third-order upwind scheme was implemented for the advection of tracers. Numerical experiments show that this scheme is computationally efficient compared to alternatives (such as the flux-corrected transport scheme) and that it works well with other aspects of the model, such as acceleration (important for spinup efficiency) and the new mixing parameterizations. The scheme mimimizes overshooting effects while keeping the dissipative aspect of the advective operator reasonably small. The net effect is to produce solutions in which the large-scale fields are affected very little while local extrema are nearly (but not completely) removed, leading to physically much more realistic tracer patterns.
