Explore the vast range of titles available.

The Canadian Earth System Model version 5 (CanESM5) is a global model developed to simulate historical climate change and variability, to make centennial-scale projections of future climate, and to produce initialized seasonal and decadal predictions. This paper describes the model components and their coupling, as well as various aspects of model development, including tuning, optimization, and a reproducibility strategy. We also document the stability of the model using a long control simulation, quantify the model's ability to reproduce large-scale features of the historical climate, and evaluate the response of the model to external forcing. CanESM5 is comprised of three-dimensional atmosphere (T63 spectral resolution equivalent roughly to 2.8∘) and ocean (nominally 1∘) general circulation models, a sea-ice model, a land surface scheme, and explicit land and ocean carbon cycle models. The model features relatively coarse resolution and high throughput, which facilitates the production of large ensembles. CanESM5 has a notably higher equilibrium climate sensitivity (5.6 K) than its predecessor, CanESM2 (3.7 K), which we briefly discuss, along with simulated changes over the historical period. CanESM5 simulations contribute to the Coupled Model Intercomparison Project phase 6 (CMIP6) and will be employed for climate science and service applications in Canada.

Anthropogenic CO2 ${\\text{CO}}_{2}$ emissions lead to ocean warming, deoxygenation and acidification. Superimposed on long‐term trends are episodic extremes of temperature, oxygen, and acidity. Here we present an innovative method for assessing single and compound extremes using a high‐resolution regional model of the Northeastern Pacific Ocean. We use an unsupervised clustering approach to identify regions with similar habitat characteristics near the seafloor, define extreme thresholds seasonally using a fixed baseline (1996–2020) within each cluster, and quantify the fraction of ocean waters that exceed these thresholds for both single and compound stressors. Compound extremes (most commonly of O2 ${\\mathrm{O}}_{2}$ and acidification) are rare but show an increasing trend in some clusters. Potential predictability of occurrence of extremes is demonstrated by correlation with basin‐scale climate variability.

We document the configuration and emergent simulation features from the Geophysical Fluid Dynamics Laboratory (GFDL) OM4.0 ocean/sea ice model. OM4 serves as the ocean/sea ice component for the GFDL climate and Earth system models. It is also used for climate science research and is contributing to the Coupled Model Intercomparison Project version 6 Ocean Model Intercomparison Project. The ocean component of OM4 uses version 6 of the Modular Ocean Model and the sea ice component uses version 2 of the Sea Ice Simulator, which have identical horizontal grid layouts (Arakawa C‐grid). We follow the Coordinated Ocean‐sea ice Reference Experiments protocol to assess simulation quality across a broad suite of climate‐relevant features. We present results from two versions differing by horizontal grid spacing and physical parameterizations: OM4p5 has nominal 0.5° spacing and includes mesoscale eddy parameterizations and OM4p25 has nominal 0.25° spacing with no mesoscale eddy parameterization. Modular Ocean Model version 6 makes use of a vertical Lagrangian‐remap algorithm that enables general vertical coordinates. We show that use of a hybrid depth‐isopycnal coordinate reduces the middepth ocean warming drift commonly found in pure z* vertical coordinate ocean models. To test the need for the mesoscale eddy parameterization used in OM4p5, we examine the results from a simulation that removes the eddy parameterization. The water mass structure and model drift are physically degraded relative to OM4p5, thus supporting the key role for a mesoscale closure at this resolution. Key Points Documentation is provided for a new generation of NOAA‐GFDL CMIP6/OMIP ocean ice climate models Dynamical core and physical parameterizations are described and key features of interannual CORE simulations are assessed Using hybrid vertical coordinates reduces spurious ocean heat drift

The ocean biogeochemistry components of two new versions of the Canadian Earth System Model (CanESM) are presented and compared to observations and other models. CanESM5 employs the same ocean biology model as CanESM2, whereas CanESM5-CanOE (Canadian Ocean Ecosystem model) is a new, more complex model developed for CMIP6, with multiple food chains, flexible phytoplankton elemental ratios, and a prognostic iron cycle. This new model is described in detail and the outputs (distributions of major tracers such as oxygen, dissolved inorganic carbon, and alkalinity, the iron and nitrogen cycles, plankton biomass, and historical trends in CO2 uptake and export production) compared to CanESM5 and CanESM2, as well as to observations and other CMIP6 models. Both CanESM5 models show gains in skill relative to CanESM2, which are attributed primarily to improvements in ocean circulation. CanESM5-CanOE shows improved skill relative to CanESM5 for most major tracers at most depths. CanESM5-CanOE includes a prognostic iron cycle, and maintains high-nutrient/low-chlorophyll conditions in the expected regions (in CanESM2 and CanESM5, iron limitation is specified as a temporally static “mask”). Surface nitrate concentrations are biased low in the subarctic Pacific and equatorial Pacific, and high in the Southern Ocean, in both CanESM5 and CanESM5-CanOE. Export production in CanESM5-CanOE is among the lowest for CMIP6 models; in CanESM5, it is among the highest, but shows the most rapid decline after about 1980. CanESM5-CanOE shows some ability to simulate aspects of plankton community structure that a single-species model can not (e.g. seasonal dominance of large cells) but is biased towards low concentrations of zooplankton and detritus relative to phytoplankton. Cumulative ocean uptake of anthropogenic carbon dioxide through 2014 is lower in both CanESM5-CanOE (122 PgC) and CanESM5 (132 PgC) than in observation-based estimates (145 PgC) or the model ensemble mean (144 PgC).

Sea levels of different atmosphere–ocean general circulation models (AOGCMs) respond to climate change forcing in different ways, representing a crucial uncertainty in climate change research. We isolate the role of the ocean dynamics in setting the spatial pattern of dynamic sea-level ( ζ ) change by forcing several AOGCMs with prescribed identical heat, momentum (wind) and freshwater flux perturbations. This method produces a ζ projection spread comparable in magnitude to the spread that results from greenhouse gas forcing, indicating that the differences in ocean model formulation are the cause, rather than diversity in surface flux change. The heat flux change drives most of the global pattern of ζ change, while the momentum and water flux changes cause locally confined features. North Atlantic heat uptake causes large temperature and salinity driven density changes, altering local ocean transport and ζ . The spread between AOGCMs here is caused largely by differences in their regional transport adjustment, which redistributes heat that was already in the ocean prior to perturbation. The geographic details of the ζ change in the North Atlantic are diverse across models, but the underlying dynamic change is similar. In contrast, the heat absorbed by the Southern Ocean does not strongly alter the vertically coherent circulation. The Arctic ζ change is dissimilar across models, owing to differences in passive heat uptake and circulation change. Only the Arctic is strongly affected by nonlinear interactions between the three air-sea flux changes, and these are model specific.

Creatine has become one of the most popular dietary supplements in the sports nutrition market. The form of creatine that has been most extensively studied and commonly used in dietary supplements is creatine monohydrate (CM). Studies have consistently indicated that CM supplementation increases muscle creatine and phosphocreatine concentrations by approximately 15-40%, enhances anaerobic exercise capacity, and increases training volume leading to greater gains in strength, power, and muscle mass. A number of potential therapeutic benefits have also been suggested in various clinical populations. Studies have indicated that CM is not degraded during normal digestion and that nearly 99% of orally ingested CM is either taken up by muscle or excreted in urine. Further, no medically significant side effects have been reported in literature. Nevertheless, supplement manufacturers have continually introduced newer forms of creatine into the marketplace. These newer forms have been purported to have better physical and chemical properties, bioavailability, efficacy, and/or safety profiles than CM. However, there is little to no evidence that any of the newer forms of creatine are more effective and/or safer than CM whether ingested alone and/or in combination with other nutrients. In addition, whereas the safety, efficacy, and regulatory status of CM is clearly defined in almost all global markets; the safety, efficacy, and regulatory status of other forms of creatine present in today's marketplace as a dietary or food supplement is less clear.

We present a neutral diffusion operator appropriate for an ocean model making use of general vertical coordinates. The diffusion scheme uses polynomial reconstructions in the vertical, along with a horizontally local but vertically nonlocal stencil for estimates of tracer fluxes. These fluxes are calculated on a vertical grid that is the superset of model columns in a neutral density space. Using flux‐limiters, the algorithm dissipates tracer extrema locally, and no new extrema are created. A demonstration using a linear equation of state in an idealized configuration shows that the algorithm is perfectly neutral. When using the nonlinear TEOS‐10 equation of state with a constant reference pressure, the algorithm compares nearly exactly to a case discretized onto isopycnal surfaces and using along‐layer diffusion. The algorithm's cost is comparable to that of tracer advection and can be readily implemented into ocean general circulation models. Plain Language Summary One prominent effect of ocean turbulence is the mixing of quantities such as temperature, salinity, and other tracers carried by the seawater. For turbulence associated with geostrophic motions, this mixing occurs along neutral directions, which are directions where buoyancy does not change. In computer models that do not explicitly represent such geostrophic eddy motions, this type of turbulent mixing is parameterized by diffusion oriented along the neutral directions. A common means to represent neutral diffusion is via rotating the diffusion fluxes to align with neutral directions (which are typically not aligned with the model grid lines). Unfortunately, this rotation can result in spurious numerical artifacts such as the creation of new minimum or maximum values of tracer concentration (referred to as extrema). In this paper, we propose a new algorithm for calculating and applying neutral diffusion in sublayers of a numerical model. We show that this new method does not create extrema, and it provides a physically accurate representation of turbulent diffusive mixing along neutral directions. Key Points A new discretization of a neutral diffusion operator is detailed, with the new method not based on a rotated tensor The new algorithm overcomes many numerical deficiencies of previous approaches Idealized test cases show minimal spurious dianeutral diffusion for linear and nonlinear equations of state

The climatological mean barotropic vorticity budget is analyzed to investigate the relative importance of surface wind stress, topography, planetary vorticity advection, and nonlinear advection in dynamical balances in a global ocean simulation. In addition to a pronounced regional variability in vorticity balances, the relative magnitudes of vorticity budget terms strongly depend on the length‐scale of interest. To carry out a length‐scale dependent vorticity analysis in different ocean basins, vorticity budget terms are spatially coarse‐grained. At length‐scales greater than 1,000 km, the dynamics closely follow the Topographic‐Sverdrup balance in which bottom pressure torque, surface wind stress curl and planetary vorticity advection terms are in balance. In contrast, when including all length‐scales resolved by the model, bottom pressure torque and nonlinear advection terms dominate the vorticity budget (Topographic‐Nonlinear balance), which suggests a prominent role of oceanic eddies, which are of O(10–100)$\\mathcal{O}(10\\mbox{--}100)$km in size, and the associated bottom pressure anomalies in local vorticity balances at length‐scales smaller than 1,000 km. Overall, there is a transition from the Topographic‐Nonlinear regime at scales smaller than 1,000 km to the Topographic‐Sverdrup regime at length‐scales greater than 1,000 km. These dynamical balances hold across all ocean basins; however, interpretations of the dominant vorticity balances depend on the level of spatial filtering or the effective model resolution. On the other hand, the contribution of bottom and lateral friction terms in the barotropic vorticity budget remains small and is significant only near sea‐land boundaries, where bottom stress and horizontal viscous friction generally peak. Plain Language Summary Vorticity provides a measure of the local rotation and shear of fluid flow. The analysis of physical processes contributing to ocean vorticity has proven fundamental to our understanding of how those processes drive ocean flows, ranging from large‐scale ocean gyres to boundary currents such as the Gulf Stream, which is tens of km in width. Furthermore, a vorticity analysis can inform us about the relative importance of different physical processes in generating flow structures having different length scales. In the present work, we perform a length‐scale dependent vorticity budget analysis using a coarse‐graining method to remove signals finer than a fixed length scale. We coarse‐grain the climatological mean vorticity budget terms over a range of length scales, and then compare the relative magnitudes to identify the dominant vorticity balances as a function of length scale. We find that the spatial structure of the meridional transport is mainly controlled by atmospheric winds, variations in ocean depth and the momentum transport by ocean currents. However, the relative magnitudes of these factors change drastically at different length scales. We conclude that physical interpretations of the primary vorticity balances are fundamentally dependent on the chosen length scale of the analysis. Key Points Relative magnitudes of barotropic vorticity budget terms display significant length‐scale dependence Bottom pressure torque and wind stress curl control the depth‐integrated meridional flow at length scales larger than 1,000 km Nonlinear advection and bottom pressure torque dominate the barotropic vorticity budget at smaller length scales

The mixing of tracers by mesoscale eddies, parameterized in many ocean general circulation models (OGCMs) as a diffusive‐advective process, contributes significantly to the distribution of tracers in the ocean. In the ocean interior, diffusive contribution occurs mostly along the direction parallel to local neutral density surfaces. However, near the surface of the ocean, small‐scale turbulence and the presence of the boundary itself break this constraint and the mesoscale transport occurs mostly along a plane parallel to the ocean surface (horizontal). Although this process is easily represented in OGCMs with geopotential vertical coordinates, the representation is more challenging in OGCMs that use a general vertical coordinate, where surfaces can be tilted with respect to the horizontal. We propose a method for representing the diffusive horizontal mesoscale fluxes within the surface boundary layer of general vertical coordinate OGCMs. The method relies on regridding/remapping techniques to represent tracers in a geopotential grid. Horizontal fluxes are calculated on this grid and then remapped back to the native grid, where fluxes are applied. The algorithm is implemented in an ocean model and tested in idealized and realistic settings. Horizontal diffusion can account for up to 10% of the total northward heat transport in the Southern Ocean and Western boundary current regions of the Northern Hemisphere. It also reduces the vertical stratification of the upper ocean, which results in an overall deepening of the surface boundary layer depth. Finally, enabling horizontal diffusion leads to meaningful reductions in the near‐surface global bias of potential temperature and salinity. Plain Language Summary Mesoscale ocean eddies, which are analogous to the weather systems in the atmosphere, are crucial to the distribution of heat, salt, carbon, and nutrients throughout the global ocean. Most of the ocean models used in climate simulations do not have enough horizontal resolution to resolve these eddies and, therefore, their effects must be parameterized. In the ocean interior, where no appreciable heat or mass is exchanged across density surfaces, the mixing of tracers due to mesoscale eddies occurs along surfaces of constant density. However, as ocean boundaries are approached, mixing occurs on a plane parallel to the boundary. For example, near the surface of the ocean, which is where the scheme presented here is designed to work, this plane is mostly horizontal. There is a class of ocean models that use a vertical coordinate system whose layer thicknesses can vary horizontally, thus complicating the implementation of horizontal diffusive parameterizations. This paper presents and evaluates a method that allows horizontal fluxes to be calculated and applied within the surface layer of this class of ocean models. Key Points General vertical coordinate ocean models lack the option to shift diffusive fluxes from neutral surfaces to a purely horizontal direction Regridding/remapping techniques are used to represent tracers in z*‐coordinates, where horizontal fluxes are easily applied The method reduces global tracer biases in forced simulations, regardless of the coordinate system employed

Background Numerous surveys show that dietary supplements are used by a large proportion of the general public, but there have been relatively few surveys on the prevalence of dietary supplement use among health professionals, including physicians and nurses. Even less information is available regarding the extent to which physicians and nurses recommend dietary supplements to their patients. Methods An online survey was administered in October 2007 to 900 physicians and 277 nurses by Ipsos Public Affairs for the Council for Responsible Nutrition (CRN), a trade association representing the dietary supplement industry. The health professionals were asked whether they used dietary supplements and their reasons for doing so, and whether they recommend dietary supplements to their patients. Results The \"Life...supplemented\" Healthcare Professionals Impact Study (HCP Impact Study) found that 72% of physicians and 89% of nurses in this sample used dietary supplements regularly, occasionally, or seasonally. Regular use of dietary supplements was reported by 51% of physicians and 59% of nurses. The most common reason given for using dietary supplements was for overall health and wellness (40% of physicians and 48% of nurses), but more than two-thirds cited more than one reason for using the products. When asked whether they \"ever recommend dietary supplements\" to their patients, 79% of physicians and 82% of nurses said they did. Conclusion Physicians and nurses are as likely as members of the general public to use dietary supplements, as shown by comparing the results of this survey with data from national health and nutrition surveys. Also, most physicians and nurses recommend supplements to their patients, whether or not the clinicians use dietary supplements themselves.