scholarly journals A Multivariate Empirical Orthogonal Function Method to Construct Nitrate Maps in the Southern Ocean

2018 ◽  
Vol 35 (7) ◽  
pp. 1505-1519 ◽  
Author(s):  
Yu-Chiao Liang ◽  
Matthew R. Mazloff ◽  
Isabella Rosso ◽  
Shih-Wei Fang ◽  
Jin-Yi Yu

AbstractThe ability to construct nitrate maps in the Southern Ocean (SO) from sparse observations is important for marine biogeochemistry research, as it offers a geographical estimate of biological productivity. The goal of this study is to infer the skill of constructed SO nitrate maps using varying data sampling strategies. The mapping method uses multivariate empirical orthogonal functions (MEOFs) constructed from nitrate, salinity, and potential temperature (N-S-T) fields from a biogeochemical general circulation model simulation Synthetic N-S-T datasets are created by sampling modeled N-S-T fields in specific regions, determined either by random selection or by selecting regions over a certain threshold of nitrate temporal variances. The first 500 MEOF modes, determined by their capability to reconstruct the original N-S-T fields, are projected onto these synthetic N-S-T data to construct time-varying nitrate maps. Normalized root-mean-square errors (NRMSEs) are calculated between the constructed nitrate maps and the original modeled fields for different sampling strategies. The sampling strategy according to nitrate variances is shown to yield maps with lower NRMSEs than mapping adopting random sampling. A k-means cluster method that considers the N-S-T combined variances to identify key regions to insert data is most effective in reducing the mapping errors. These findings are further quantified by a series of mapping error analyses that also address the significance of data sampling density. The results provide a sampling framework to prioritize the deployment of biogeochemical Argo floats for constructing nitrate maps.

1997 ◽  
Vol 25 ◽  
pp. 111-115 ◽  
Author(s):  
Achim Stössel

This paper investigates the long-term impact of sea ice on global climate using a global sea-ice–ocean general circulation model (OGCM). The sea-ice component involves state-of-the-art dynamics; the ocean component consists of a 3.5° × 3.5° × 11 layer primitive-equation model. Depending on the physical description of sea ice, significant changes are detected in the convective activity, in the hydrographic properties and in the thermohaline circulation of the ocean model. Most of these changes originate in the Southern Ocean, emphasizing the crucial role of sea ice in this marginally stably stratified region of the world's oceans. Specifically, if the effect of brine release is neglected, the deep layers of the Southern Ocean warm up considerably; this is associated with a weakening of the Southern Hemisphere overturning cell. The removal of the commonly used “salinity enhancement” leads to a similar effect. The deep-ocean salinity is almost unaffected in both experiments. Introducing explicit new-ice thickness growth in partially ice-covered gridcells leads to a substantial increase in convective activity, especially in the Southern Ocean, with a concomitant significant cooling and salinification of the deep ocean. Possible mechanisms for the resulting interactions between sea-ice processes and deep-ocean characteristics are suggested.


2010 ◽  
Vol 40 (5) ◽  
pp. 880-899 ◽  
Author(s):  
Matthew R. Mazloff ◽  
Patrick Heimbach ◽  
Carl Wunsch

Abstract An eddy-permitting general circulation model of the Southern Ocean is fit by constrained least squares to a large observational dataset during 2005–06. Data used include Argo float profiles, CTD synoptic sections, Southern Elephant Seals as Oceanographic Samplers (SEaOS) instrument-mounted seal profiles, XBTs, altimetric observations [Envisat, Geosat, Jason-1, and Ocean Topography Experiment (TOPEX)/Poseidon], and infrared and microwave radiometer observed sea surface temperature. An adjoint model is used to determine descent directions in minimizing a misfit function, each of whose elements has been weighted by an estimate of the observational plus model error. The model is brought into near agreement with the data by adjusting its control vector, here consisting of initial and meteorological boundary conditions. Although total consistency has not yet been achieved, the existing solution is in good agreement with the great majority of the 2005 and 2006 Southern Ocean observations and better represents these data than does the World Ocean Atlas 2001 (WOA01) climatological product. The estimate captures the oceanic temporal variability and in this respect represents a major improvement upon earlier static inverse estimates. During the estimation period, the Drake Passage volume transport is 153 ± 5 Sv (1 Sv ≡ 106 m3 s−1). The Ross and Weddell polar gyre transports are 20 ± 5 Sv and 40 ± 8 Sv, respectively. Across 32°S there is a surface meridional overturning cell of 12 ± 12 Sv, an intermediate cell of 17 ± 12 Sv, and an abyssal cell of 13 ± 6 Sv. The northward heat and freshwater anomaly transports across 30°S are −0.3 PW and 0.7 Sv, with estimated uncertainties of 0.5 PW and 0.2 Sv. The net rate of wind work is 2.1 ± 1.1 TW. Southern Ocean theories involving short temporal- and spatial-scale dynamics may now be tested with a dynamically and thermodynamically realistic general circulation model solution that is known to be compatible with the modern observational datasets.


2018 ◽  
Vol 116 (37) ◽  
pp. 18251-18256 ◽  
Author(s):  
F. J. Beron-Vera ◽  
A. Hadjighasem ◽  
Q. Xia ◽  
M. J. Olascoaga ◽  
G. Haller

The emergence of coherent Lagrangian swirls (CLSs) among submesoscale motions in the ocean is illustrated. This is done by applying recent nonlinear dynamics tools for Lagrangian coherence detection on a surface flow realization produced by a data-assimilative submesoscale-permitting ocean general circulation model simulation of the Gulf of Mexico. Both mesoscale and submesoscale CLSs are extracted. These extractions prove the relevance of coherent Lagrangian eddies detected in satellite-altimetry–based geostrophic flow data for the arguably more realistic ageostrophic multiscale flow.


2019 ◽  
Vol 632 ◽  
pp. A114 ◽  
Author(s):  
F. Sainsbury-Martinez ◽  
P. Wang ◽  
S. Fromang ◽  
P. Tremblin ◽  
T. Dubos ◽  
...  

Context. The anomalously large radii of hot Jupiters has long been a mystery. However, by combining both theoretical arguments and 2D models, a recent study has suggested that the vertical advection of potential temperature leads to a hotter adiabatic temperature profile in the deep atmosphere than the profile obtained with standard 1D models. Aims. In order to confirm the viability of that scenario, we extend this investigation to 3D, time-dependent models. Methods. We use a 3D general circulation model DYNAMICO to perform a series of calculations designed to explore the formation and structure of the driving atmospheric circulations, and detail how it responds to changes in both the upper and deep atmospheric forcing. Results. In agreement with the previous, 2D study, we find that a hot adiabat is the natural outcome of the long-term evolution of the deep atmosphere. Integration times of the order of 1500 yr are needed for that adiabat to emerge from an isothermal atmosphere, explaining why it has not been found in previous hot Jupiter studies. Models initialised from a hotter deep atmosphere tend to evolve faster toward the same final state. We also find that the deep adiabat is stable against low-levels of deep heating and cooling, as long as the Newtonian cooling timescale is longer than ~3000 yr at 200 bar. Conclusions. We conclude that steady-state vertical advection of potential temperature by deep atmospheric circulations constitutes a robust mechanism to explain the inflated radii of hot Jupiters. We suggest that future models of hot Jupiters be evolved for a longer time than currently done, and when possible that models initialised with a hot deep adiabat be included. We stress that this mechanism stems from the advection of entropy by irradiation-induced mass flows and does not require a (finely tuned) dissipative process, in contrast with most previously suggested scenarios.


2019 ◽  
Author(s):  
Camille Risi ◽  
Joseph Galewsky ◽  
Gilles Reverdin ◽  
Florent Brient

Abstract. Understanding what controls the water vapor isotopic composition of the sub-cloud layer (SCL) over tropical oceans (δD0) is a first step towards understanding the water vapor isotopic composition everywhere in the troposphere. We propose an analytical model to predict δD0 as a function of sea surface conditions, humidity and temperature profiles, and the altitude from which the free tropospheric air originates (zorig). To do so, we extend previous studies by (1) prescribing the shape of δD0 vertical profiles, and (2) linking δD0 to zorig. The model relies on the hypotheses that δD0 profiles are steeper than mixing lines and no clouds are precipitating. We show that δD0 does not depend on the intensity of entrainment, dampening hope that δD0 measurements could help constrain this long-searched quantity. Based on an isotope-enabled general circulation model simulation, we show that δD0 variations are mainly controlled by mid-tropospheric depletion and rain evaporation in ascending regions, and by sea surface temperature and zorig in subsiding regions. When the air mixing into the SCL is lower in altitude, it is moister, and thus it depletes more efficiently the SCL. In turn, could δD0 measurements help estimate zorig and thus discriminate between different mixing processes? Estimates that are accurate enough to be useful would be difficult to achieve in practice, requiring measuring daily δD profiles, and measuring δD0 with an accuracy of 0.1 ‰ and 0.4 ‰ in trade-wind cumulus and strato-cumulus clouds respectively.


2019 ◽  
Vol 49 (10) ◽  
pp. 2553-2570 ◽  
Author(s):  
Mads B. Poulsen ◽  
Markus Jochum ◽  
James R. Maddison ◽  
David P. Marshall ◽  
Roman Nuterman

AbstractAn interpretation of eddy form stress via the geometry described by the Eliassen–Palm flux tensor is explored. Complimentary to previous works on eddy Reynolds stress geometry, this study shows that eddy form stress is fully described by a vertical ellipse, whose size, shape, and orientation with respect to the mean flow shear determine the strength and direction of vertical momentum transfers. Following a recent proposal, this geometric framework is here used to form a Gent–McWilliams eddy transfer coefficient that depends on eddy energy and a nondimensional geometric parameter α, bounded in magnitude by unity. The parameter α expresses the efficiency by which eddies exchange energy with baroclinic mean flow via along-gradient eddy buoyancy flux—a flux equivalent to eddy form stress along mean buoyancy contours. An eddy-resolving ocean general circulation model is used to estimate the spatial structure of α in the Southern Ocean and assess its potential to form a basis for parameterization. The eddy efficiency α averages to a low but positive value of 0.043 within the Antarctic Circumpolar Current, consistent with an inefficient eddy field extracting energy from the mean flow. It is found that the low eddy efficiency is mainly the result of that eddy buoyancy fluxes are weakly anisotropic on average. The eddy efficiency is subject to pronounced vertical structure and is maximum at ~3-km depth, where eddy buoyancy fluxes tend to be directed most downgradient. Since α partly sets the eddy form stress in the Southern Ocean, a parameterization for α must reproduce its vertical structure to provide a faithful representation of vertical stress divergence and eddy forcing.


1998 ◽  
Vol 11 (8) ◽  
pp. 1883-1905 ◽  
Author(s):  
O. P. Sharma ◽  
H. Le Treut ◽  
G. Sèze ◽  
L. Fairhead ◽  
R. Sadourny

Abstract The sensitivity of the interannual variations of the summer monsoons to imposed cloudiness has been studied with a general circulation model using the initial conditions prepared from the European Centre for Medium-Range Forecasts analyses of 1 May 1987 and 1988. The cloud optical properties in this global model are calculated from prognostically computed cloud liquid water. The model successfully simulates the contrasting behavior of these two successive monsoons. However, when the optical properties of the observed clouds are specified in the model runs, the simulations show some degradation over India and its vicinity. The main cause of this degradation is the reduced land–sea temperature contrast resulting from the radiative effects of the observed clouds imposed in such simulations. It is argued that the high concentration of condensed water content of clouds over the Indian land areas will serve to limit heating of the land, thereby reducing the thermal contrast that gives rise to a weak Somali jet. A countermonsoon circulation is, therefore, simulated in the vector difference field of 850-hPa winds from the model runs with externally specified clouds. This countermonsoon circulation is associated with an equatorial heat source that is the response of the model to the radiative effects of the imposed clouds. Indeed, there are at least two clear points that can be made: 1) the cloud–SST patterns, together, affect the interannual variability; and 2) with both clouds and SST imposed, the model simulation is less sensitive to initial conditions. Additionally, the study emphasizes the importance of dynamically consistent clouds developing in response to the dynamical, thermal, and moist state of the atmosphere during model integrations.


2009 ◽  
Vol 6 (3) ◽  
pp. 333-348 ◽  
Author(s):  
L. Patara ◽  
N. Pinardi ◽  
C. Corselli ◽  
E. Malinverno ◽  
M. Tonani ◽  
...  

Abstract. This paper analyzes the relationship between deep sedimentary fluxes and ocean current vertical velocities in an offshore area of the Ionian Sea, the deepest basin of the Eastern Mediterranean Sea. Sediment trap data are collected at 500 m and 2800 m depth in two successive moorings covering the period September 1999–May 2001. A tight coupling is observed between the upper and deep traps and the estimated particle sinking rates are more than 200 m day−1. The current vertical velocity field is computed from a 1/16°×1/16° Ocean General Circulation Model simulation and from the wind stress curl. Current vertical velocities are larger and more variable than Ekman vertical velocities, yet the general patterns are alike. Current vertical velocities are generally smaller than 1 m day−1: we therefore exclude a direct effect of downward velocities in determining high sedimentation rates. However we find that upward velocities in the subsurface layers of the water column are positively correlated with deep particle fluxes. We thus hypothesize that upwelling would produce an increase in upper ocean nutrient levels – thus stimulating primary production and grazing – a few weeks before an enhanced vertical flux is found in the sediment traps. High particle sedimentation rates may be attained by means of rapidly sinking fecal pellets produced by gelatinous macro-zooplankton. Other sedimentation mechanisms, such as dust deposition, are also considered in explaining large pulses of deep particle fluxes. The fast sinking rates estimated in this study might be an evidence of the efficiency of the biological pump in sequestering organic carbon from the surface layers of the deep Eastern Mediterranean basins.


2019 ◽  
Vol 76 (6) ◽  
pp. 1547-1564 ◽  
Author(s):  
Spencer A. Hill ◽  
Simona Bordoni ◽  
Jonathan L. Mitchell

Abstract We consider the relevance of known constraints from each of Hide’s theorem, the angular momentum–conserving (AMC) model, and the equal-area model on the extent of cross-equatorial Hadley cells. These theories respectively posit that a Hadley circulation must span all latitudes where the radiative–convective equilibrium (RCE) absolute angular momentum satisfies or or where the RCE absolute vorticity satisfies ; all latitudes where the RCE zonal wind exceeds the AMC zonal wind; and over a range such that depth-averaged potential temperature is continuous and that energy is conserved. The AMC model requires knowledge of the ascent latitude , which needs not equal the RCE forcing maximum latitude . Whatever the value of , we demonstrate that an AMC cell must extend at least as far into the winter hemisphere as the summer hemisphere. The equal-area model predicts , always placing it poleward of . As is moved poleward (at a given thermal Rossby number), the equal-area-predicted Hadley circulation becomes implausibly large, while both and become increasingly displaced poleward of the minimal cell extent based on Hide’s theorem (i.e., of supercritical forcing). In an idealized dry general circulation model, cross-equatorial Hadley cells are generated, some spanning nearly pole to pole. All homogenize angular momentum imperfectly, are roughly symmetric in extent about the equator, and appear in extent controlled by the span of supercritical forcing.


2020 ◽  
Vol 20 (11) ◽  
pp. 6607-6630 ◽  
Author(s):  
Peter Kuma ◽  
Adrian J. McDonald ◽  
Olaf Morgenstern ◽  
Simon P. Alexander ◽  
John J. Cassano ◽  
...  

Abstract. Southern Ocean (SO) shortwave (SW) radiation biases are a common problem in contemporary general circulation models (GCMs), with most models exhibiting a tendency to absorb too much incoming SW radiation. These biases have been attributed to deficiencies in the representation of clouds during the austral summer months, either due to cloud cover or cloud albedo being too low. The problem has been the focus of many studies, most of which utilised satellite datasets for model evaluation. We use multi-year ship-based observations and the CERES spaceborne radiation budget measurements to contrast cloud representation and SW radiation in the atmospheric component Global Atmosphere (GA) version 7.1 of the HadGEM3 GCM and the MERRA-2 reanalysis. We find that the prevailing bias is negative in GA7.1 and positive in MERRA-2. GA7.1 performs better than MERRA-2 in terms of absolute SW bias. Significant errors of up to 21 W m−2 (GA7.1) and 39 W m−2 (MERRA-2) are present in both models in the austral summer. Using ship-based ceilometer observations, we find low cloud below 2 km to be predominant in the Ross Sea and the Indian Ocean sectors of the SO. Utilising a novel surface lidar simulator developed for this study, derived from an existing Cloud Feedback Model Intercomparison Project (CFMIP) Observation Simulator Package (COSP) – active remote sensing simulator (ACTSIM) spaceborne lidar simulator, we find that GA7.1 and MERRA-2 both underestimate low cloud and fog occurrence relative to the ship observations on average by 4 %–9 % (GA7.1) and 18 % (MERRA-2). Based on radiosonde observations, we also find the low cloud to be strongly linked to boundary layer atmospheric stability and the sea surface temperature. GA7.1 and MERRA-2 do not represent the observed relationship between boundary layer stability and clouds well. We find that MERRA-2 has a much greater proportion of cloud liquid water in the SO in austral summer than GA7.1, a likely key contributor to the difference in the SW radiation bias. Our results suggest that subgrid-scale processes (cloud and boundary layer parameterisations) are responsible for the bias and that in GA7.1 a major part of the SW radiation bias can be explained by cloud cover underestimation, relative to underestimation of cloud albedo.


Sign in / Sign up

Export Citation Format

Share Document