Displacements Of Parcels Research Articles

Stokes drift and particle trajectories induced by surface waves atop a shear flow

Surface waves and currents are crucial to the mass transfer in the air-sea interaction as they can drive a variety of dynamical processes. How mass can be transported by surface waves and current coupling is addressed through a study of their induced motions of fluid parcels. To this end, a weakly nonlinear wavetrain is imposed on the background flow whose direction and magnitude are permitted to vary with water depth and second-order features of this configuration are investigated. A leading-order approximation to the Stokes drift is derived, correct to the second order in wave steepness, and applicable to an arbitrarily depth-dependent background flow. The reduced forms of the approximate Stokes drift are provided in a few limiting cases such as a current with an exponential profile or propagating in an orthogonal direction to the wave propagation. Novel features related to the Stokes drift and particle trajectories have been reported for the first time as a result of the rotation induced by the wave and current coupling. A non-vanishing component of the Stokes drift velocity and net-mean displacement of fluid parcels in the span-wise direction to the wave propagation are observed in the cases where a shear current propagates obliquely to the waves direction. A non-monotonic dependence on water depth of the stream-wise component of the Stokes drift is shown, and thereby the largest mass transport induced no longer occurs on the still water surface but some depth beneath. The non-monotonic behavior occurs beyond the regime of the near-irrotational assumption of wave-induced motions. It can also lead to the change of the signs for the stream-wise Stokes drift throughout the water column, and thus an overall cancellation of the integrated mass transport by waves over the water column, indicating that the depth-integrated models can likely lead to underestimated effects of the mass transport which is non-trivial at a local depth. The results from this study have far-reaching impact. The Stokes drift profile is a direct input to the parametrization of the surface waves forcing in ocean circulations and the obliquely propagating Stokes drift can be plausibly responsible for the formation of oblique Langmuir rolls to wave propagation in the open ocean.

The Influence of Atmospheric Bores on Nocturnal Convection Initiation in Weakly Forced Synoptic Environments

Abstract This study creates a composite sounding for nocturnal convection initiation (CI) events under weakly forced conditions and utilizes an idealized numerical simulation to assess the impact of atmospheric bores on these environments. Thirteen soundings were used to create this composite sounding. Common conditions associated with these weakly forced environments include a nocturnal low-level jet and a Brunt–Väisälä frequency of 0.011 s−1 above 900 hPa. The median lift needed for parcels to realize any convective instability is 490 m, the median convective available potential energy of these convectively unstable parcels is 992 J kg−1, and the median initial pressure of these parcels is 800 hPa. An idealized numerical simulation was utilized to examine the potential influence of bores on CI in an environment based on composite sounding. The characteristics of the simulated bore were representative of observed bores. The vertical velocities associated with this simulated bore were between 1 and 2 m s−1, and the net upward displacement of parcels was between 400 and 650 m. The vertical displacement of air parcels has two notable phases: lift by the bore itself and smaller-scale lift that occurs 100–150 km ahead of the bore passage. The prebore lift is between 50 and 200 m and appears to be related to low-frequency waves ahead of the bores. The lift with these waves was maximized in the low to midtroposphere between 1 and 4 km AGL, and this lift may play a role in assisting CI in these otherwise weakly forced environments.

Reinterpretation of the Thorpe Length Scale

Abstract In 1977, S. A. Thorpe proposed a method to estimate the dissipation rate ε of turbulence kinetic energy (TKE) in an overturning turbulent layer in a lake, by sorting the observed (unstable) density profile to render it stable and thus deriving a length scale LT named after him, from the resulting vertical displacements of water parcels. By further proposing that this purely empirical scale (with no a priori physical basis, unlike many other turbulence length scales) is proportional to the Ozmidov scale LO, definable only for stably (not unstably or neutrally) stratified flows, he was able to extract ε. The simplicity of the approach that requires nothing but CTD (Conductivity, Temperature and Depth) casts in water bodies, including lakes and oceans, made it attractive, until microstructure profilers were developed and perfected in later decades to actually make in-situ measurements of ε. Since equivalent microstructure devices are not available for the atmosphere, Thorpe technique has been resurrected in recent years for application to the atmosphere, using potential temperature profiles obtained from high vertical resolution radiosondes. Its popularity and utility have increased lately, in spite of unresolved issues related to the validity of assuming LT is proportional to LO. In this study, we touch upon these issues and offer an alternative interpretation of the Thorpe length scale as indicative of the turbulence velocity scale σK, which allows Thorpe sorting technique to be applied to all turbulent flows, including those generated by convection.

On the Local Available Potential Energy Perspective of Baroclinic Wave Development

Abstract Extratropical cyclones convert available potential energy (APE) to kinetic energy. However, our current understanding of APE conversion on synoptic scales is limited, as the well-established Lorenz APE framework is only applicable in a global, volume-integrated sense. Here, we employ a recently developed local APE framework to investigate APE and its tendencies in a highly idealized, dispersive baroclinic wave, which leads to the formation of a primary and a downstream cyclone. By utilizing a Lagrangian approach, we demonstrate that locally the downstream cyclone not only consumes APE but also generates it. Initially, APE is transported from both poleward and equatorward reservoirs into the baroclinic zone, where it is then consumed by the vertical displacement of air parcels associated with the developing cyclone. To a lesser extent, APE is also created within the cyclone when air parcels overshoot their reference state; i.e., air colder than its reference state is lifted and air warmer than its reference state is lowered. The volume integral of the APE tendency is dominated by slow vertical displacements of large air masses, whereas the dry intrusion (DI) and warm conveyor belt (WCB) of the cyclone are responsible for the largest local APE tendencies. Diabatic effects within the DI and WCB contribute to the generation of APE in regions where it is consumed adiabatically, thereby enhancing baroclinic conversion in situ. Our findings provide a comprehensive and mechanistic understanding of the local APE tendency on synoptic scales within an idealized setting and complement existing frameworks explaining the energetics of cyclone intensification.

Vertical Motions in Orographic Cloud Systems over the Payette River Basin. Part IV: Controls on Supercooled Liquid Water Content and Cloud Droplet Number Concentrations

Abstract This paper examines the controls on supercooled liquid water content (SLWC) and drop number concentrations (Nt,CDP) over the Payette River basin during the Seeded and Natural Orographic Wintertime Clouds: The Idaho Experiment (SNOWIE) campaign. During SNOWIE, 27.4% of 1-Hz in situ cloud droplet probe samples were in an environment containing supercooled liquid water (SLW). The interquartile range of SLWC, when present, was found to be 0.02–0.18 g m−3 and 13.3–37.2 cm−3 for Nt,CDP, with the most extreme values reaching 0.40–1.75 g m−3 and 150–320 cm−3 in isolated regions of convection and strong shear-induced turbulence. SLWC and Nt,CDP distributions are shown to be directly related to cloud-top temperature and ice particle concentrations, consistent with past research over other mountain ranges. Two classes of vertical motions were analyzed as potential controls on SLWC and Nt,CDP, the first forced by the orography and fixed in space relative to the topography (stationary waves) and the second transient, triggered by vertical shear and instability within passing synoptic-scale cyclones. SLWC occurrence and magnitudes, and Nt,CDP associated with fixed updrafts were found to be normally distributed about ridgelines when SLW was present. SLW was more likely to form at low altitudes near the terrain slope associated with fixed waves due to higher mixing ratios and larger vertical air parcel displacements at low altitudes. When considering transient updrafts, SLWC and Nt,CDP appear more uniformly distributed over the flight track with little discernable terrain dependence as a result of time and spatially varying updrafts associated with passing weather systems. The implications for cloud seeding over the basin are discussed.

Lagrangian description of the atmospheric flow from Eulerian tracer advection with relaxation

AbstractThe Lagrangian representation of fluid flows offers a natural perspective to study many kinds of physical mechanisms. By contrast, the Eulerian representation is more convenient from a diagnostic point of view. This article attempts to combine elements of both worlds by proposing an Eulerian method that allows one to extract Lagrangian information about the atmospheric flow. The method is based on the offline advection of passive tracer fields and includes a relaxation term. The latter device allows one to run the integration in a continuous fashion without the need for reinitialization. As a result one obtains accumulated Lagrangian information, for example, about the recent parcel displacement or the recent parcel‐based diabatic heating, at each point of an Eulerian grid at any time step. The method is implemented with a pseudospectral algorithm suitable for gridded global atmospheric data and compared with the more traditional trajectory method. The method's utility is demonstrated on the basis of a few examples, which relate to cloud formation and the development of temperature anomalies. The examples highlight that the method provides a convenient diagnostic of parcel‐based changes, paving an intuitive way to explore the physical processes involved. Due to its gridpoint‐based nature, the proposed method can be applied to large data sets in a straightforward and computationally efficient manner, suggesting that the method is particularly useful for climatological analyses.

Three-Dimensional Static Instability of Gravity Waves and a Possible Parameterization of the Associated Wave Breaking

Abstract Parameterizations of subgrid-scale gravity waves (GWs) in atmospheric models commonly involve the description of the dissipation of GWs. Where they dissipate, GWs have an increased effect on the large-scale flow. Instabilities that trigger wave breaking are an important starting point for the route to dissipation. Possible destabilizing mechanisms are numerous, but the classical vertical static instability is still regarded as a key indicator for the disposition to wave breaking. In this work, we investigate how the horizontal variations associated with a GW could alter the criterion for static instability. To this end, we use an extension of the common parcel displacement method. This three-dimensional static stability analysis predicts a significantly larger range of instability than does the vertical static stability analysis. In this case, the Lindzen-type saturation adjustment to a state of marginal stability is perhaps a less suitable ansatz for the parameterization of the GW breaking. To develop a possible ansatz for the GW dissipation due to three-dimensional instability, we apply the methods of irreversible thermodynamics, which are embedded in the Gibbs formalism of dynamics. In this way, the parameterization does not only satisfy the second law of thermodynamics, but it can also be made consistent with the conservation of energy and further (non-)conservation principles. We develop the parameterization for a discrete spectrum of GW packets. Offline computations of GW drag and dissipative heating rates are performed for two vertical profiles of zonal wind and temperature for summer and winter conditions from CIRA data. The results are compared to benchmarks from the literature.

Rapid vertical exchange at fronts in the Northern Gulf of Mexico.

Over the Texas-Louisiana Shelf in the Northern Gulf of Mexico, the eutrophic, fresh Mississippi/Atchafalaya river plume isolates saltier waters below, supporting the formation of bottom hypoxia in summer. The plume also generates strong density fronts, features of the circulation that are known pathways for the exchange of water between the ocean surface and the deep. Using high-resolution ocean observations and numerical simulations, we demonstrate how the summer land-sea breeze generates rapid vertical exchange at the plume fronts. We show that the interaction between the land-sea breeze and the fronts leads to convergence/divergence in the surface mixed layer, which further facilitates a slantwise circulation that subducts surface water along isopycnals into the interior and upwells bottom waters to the surface. This process causes significant vertical displacements of water parcels and creates a ventilation pathway for the bottom water in the northern Gulf. The ventilation of bottom water can bypass the stratification barrier associated with the Mississippi/Atchafalaya river plume and might impact the dynamics of the region’s dead zone.

Stratified flow past a sphere at moderate Reynolds numbers

A numerical study of stably stratified flows past spheres at Reynolds numbers Re=200 and Re=300 is reported. In these flow regimes, a neutrally stratified laminar flow induces distinctly different near-wake features. However, the flow behaviour changes significantly as the stratification increases and suppresses the scale of vertical displacements of fluid parcels. Computations for a range of Froude numbers Fr∈[0.1,∞] show that as Froude number decreases, the flow patterns for both Reynolds numbers become similar. The representative simulations of the lee-wave instability at Fr=0.625 and the two-dimensional vortex shedding at Fr=0.25 regimes are illustrated for flows past single and tandem spheres, thereby providing further insight into the dynamics of stratified flows past bluff bodies. In particular, the reported study examines the relative influence of viscosity and stratification on the dividing streamline elevation, wake structure and flow separation. The solutions of the Navier-Stokes equations in the incompressible Boussinesq limit are obtained on unstructured meshes suitable for simulations involving multiple bodies. Computations are accomplished using the finite volume, non-oscillatory forward-in-time (NFT) Multidimensional Positive Definite Transport Algorithm (MPDATA) based solver. The impact and validity of the numerical approximations, especially for the cases exhibiting strong stratification, are also discussed. Qualitative and quantitative comparisons with available laboratory experiments and prior numerical studies confirm the validity of the numerical approach.

The Role of Wind Speed and Wind Shear for Banner Cloud Formation

Abstract Banner clouds are clouds that appear to be attached to the leeward face of a steep mountain. This paper investigates the role of wind speed and wind shear for the formation of banner clouds. Large-eddy simulations are performed to simulate the flow of dry air past an idealized pyramid-shaped mountain. The potential for cloud formation is diagnosed through the Lagrangian vertical parcel displacement, which in the case of a banner cloud shows a plume of large values in the lee of the mountain. In addition, vortical structures are visualized through streamlines and their curvature. A series of sensitivity experiments indicates that both the flow and the banner cloud occurrence are largely independent of the ambient wind speed U. On the other hand, the shear of the ambient wind has a profound impact on the location of the stagnation point on the windward face as well as on the flow geometry in the lee of the mountain. The relevant measure for shear is H/Hs, where H denotes the height of the mountain and Hs = U/Uz is the scale height of the shear (with Uz denoting the scale of the shear). The simulations are also used to compute the line-of-sight velocity component seen by a hypothetical Doppler wind lidar positioned in the lee of the mountain; the analysis suggests that such sensitivities can potentially be detected using modern wind lidar technology.

A Composite Perspective on Bore Passages during the PECAN Campaign

AbstractAtmospheric bores have been shown to have a role in the initiation and maintenance of elevated convection. Previous observational studies of bores have been case studies of more notable events. However, this creates a selection bias toward extraordinary cases, while discussions of the differences between bores that favor convective initiation and maintenance and bores that do not are lacking from the literature. This study attempts to fill that gap by analyzing a high-temporal-resolution thermodynamic profile composite of eight bores observed by multiple platforms during the Plains Elevated Convection at Night (PECAN) campaign in order to assess the impact of bores on the environment. The time–height cross section of the potential temperature composite displays quasi-permanent parcel displacements up to 900 m with the bore passage. Low-level lifting is shown to weaken the capping inversion and reduce convective inhibition (CIN) and the level of free convection (LFC). Additionally, low-level water vapor increases by about 1 g kg−1 in the composite mean. By assessing variability across the eight cases, it is shown that increases in low-level water vapor result in increases to convective available potential energy (CAPE), while drying results in decreased CAPE. Most cases resulted in decreased CIN and LFC height with the bore passage, but only some cases resulted in increased CAPE. This suggests that bores will increase the potential for convective initiation, but future research should be directed toward better understanding cases that result in increased CAPE as those are the types of bores that will increase severity of convection.

Tidally-generated internal waves in Southeast Hudson Bay

The location of the amphidromic point of the M2 tide in Hudson Bay roughly coincides with Belcher Islands, a region where the surface mixed layer stays relatively fresh throughout summer and winter due to significant ice melt and river discharge. High-resolution satellite radar imagery for the ice-free season revealed that the coastal region in the south-east Belcher Islands is a hot spot for short-period internal wave activity. For a first investigation of tidal dynamics in the region, we took advantage of the sea ice platform to deploy an ice-tethered mooring consisting of nine conductivity and temperature sensors and an acoustic Doppler current profiler. The mooring was deployed at 65 m depth in January–March 2014 in a narrow channel between Broomfield and O′Leary islands located in the south-east tip of the Belcher Islands group in Hudson Bay (56°20′N, 79°30′W), northeast Canada. The surface mixed layer under the land-fast ice in this area stays relatively fresh through winter presumably because of significant winter river discharge in nearby James Bay. The mooring recorded oscillations of temperature and salinity throughout the whole water column, which were attributed to vertical displacement caused by internal tidal waves. The tidal harmonic analysis performed for the M2 tidal constituent showed the dominance of the baroclinic tide with maximum velocity amplitudes at the surface and decreasing with depth. Vertical displacements of water parcels derived from both temperature and salinity were statistically similar and displayed the maximum values of 11.9 m at 35 m (instrument depth). The combination of winter hydrographic data and summer satellite observations confirmed that the observed internal waves were generated by the interaction of strong tides, typical for Hudson Bay, with highly variable bottom topography south-east of the Belcher Islands archipelago.

Mountain-Wave Propagation under Transient Tropospheric Forcing: A DEEPWAVE Case Study

The impact of transient tropospheric forcing on the deep vertical mountain-wave propagation is investigated by a unique combination of in situ and remote sensing observations and numerical modeling. The temporal evolution of the upstream low-level wind follows approximately a [Formula: see text] shape and was controlled by a migrating trough and connected fronts. Our case study reveals the importance of the time-varying propagation conditions in the upper troposphere and lower stratosphere (UTLS). Upper-tropospheric stability, the wind profile, and the tropopause strength affected the observed and simulated wave response in the UTLS. Leg-integrated along-track momentum fluxes ([Formula: see text]) and amplitudes of vertical displacements of air parcels in the UTLS reached up to 130 kN m−1 and 1500 m, respectively. Their maxima were phase shifted to the maximum low-level forcing by ≈8 h. Small-scale waves ([Formula: see text] km) were continuously forced, and their flux values depended on wave attenuation by breaking and reflection in the UTLS region. Only maximum flow over the envelope of the mountain range favored the excitation of longer waves that propagated deeply into the mesosphere. Their long propagation time caused a retarded enhancement of observed mesospheric gravity wave activity about 12–15 h after their observation in the UTLS. For the UTLS, we further compared observed and simulated [Formula: see text] with fluxes of 2D quasi-steady runs. UTLS momentum fluxes seem to be reproducible by individual quasi-steady 2D runs, except for the flux enhancement during the early decelerating forcing phase.

Generation of Hydroacoustic Waves by an Oscillating Ice Block in Arctic Zones

The time harmonic problem of propagating hydroacoustic waves generated in the ocean by a vertically oscillating ice block in arctic zones is discussed. The generated acoustic modes can result in orbital displacements of fluid parcels sufficiently high that may contribute to deep ocean currents and circulation. This mechanism adds to current efforts for explaining ocean circulation from a snowball earth Neoproterozoic Era to greenhouse earth arctic conditions and raises a challenge as the extent of ice blocks shrinks towards an ice-free sea. Surprisingly, unlike the free-surface setting, here it is found that the higher acoustic modes exhibit a larger contribution.

Mechanisms for Organization and Echo Training in a Flash-Flood-Producing Mesoscale Convective System

Abstract In this research, a numerical simulation of an observed training line/adjoining stratiform (TL/AS)-type mesoscale convective system (MCS) was used to investigate processes leading to upwind propagation of convection and quasi-stationary behavior. The studied event produced damaging flash flooding near Dubuque, Iowa, on the morning of 28 July 2011. The simulated convective system well emulated characteristics of the observed system and produced comparable rainfall totals. In the simulation, there were two cold pool–driven convective surges that exited the region where heavy rainfall was produced. Low-level unstable flow, which was initially convectively inhibited, overrode the surface cold pool subsequent to these convective surges, was gradually lifted to the point of saturation, and reignited deep convection. Mechanisms for upstream lifting included persistent large-scale warm air advection, displacement of parcels over the surface cold pool, and an upstream mesolow that formed between 0500 and 1000 UTC. Convection tended to propagate with the movement of the southeast portion of the outflow boundary, but did not propagate with the southwest outflow boundary. This was explained by the vertical wind shear profile over the depth of the cold pool being favorable (unfavorable) for initiation of new convection along the southeast (southwest) cold pool flank. A combination of a southward-oriented pressure gradient force in the cold pool and upward transport of opposing southerly flow away from the boundary layer moved the outflow boundary southward. Upward transport of southerly momentum by convection along the southward-moving outflow boundary, along with convectively induced southward pressure gradient forces cut off southerly inflow to the MCS, which temporarily disrupted backbuilding.

On the validity of single-parcel energetics to assess the importance of internal energy and compressibility effects in stratified fluids

AbstractIt is often assumed on the basis of single-parcel energetics that compressible effects and conversions with internal energy are negligible whenever typical displacements of fluid parcels are small relative to the scale height of the fluid (defined as the ratio of the squared speed of sound to the gravitational acceleration). This paper shows that the above approach is flawed, however, and that a correct assessment of compressible effects and internal energy conversions requires the consideration of the energetics of at least two parcels or, more generally, of mass-conserving parcel rearrangements. As a consequence, it is shown that it is the adiabatic lapse rate and its derivative with respect to pressure, rather than the scale height, that controls the relative importance of compressible effects and internal energy conversions when considering the global energy budget of a stratified fluid. Only when mass conservation is properly accounted for is it possible to explain why the available internal energy can account for up to 40 % of the total available potential energy in the oceans. This is considerably larger than the prediction of single-parcel energetics, according to which this number should be no more than approximately 2 %.

A Semihydrostatic Theory of Gravity-Dominated Compressible Flow

Abstract From Hamilton’s least-action principle, compressible equations of motion with density diagnosed from potential temperature through hydrostatic balance are derived. Slaving density to potential temperature suppresses the degrees of freedom supporting the propagation of acoustic waves and results in a soundproof system. The linear normal modes and dispersion relationship for an isothermal state of rest on f and β planes are accurate from hydrostatic to nonhydrostatic scales, except for deep internal gravity waves. Specifically, the Lamb wave and long Rossby waves are not distorted, unlike with anelastic or pseudoincompressible systems. Compared to similar equations derived by A. Arakawa and C. S. Konor, the semihydrostatic system derived here possesses an additional term in the horizontal momentum budget. This term is an apparent force resulting from the vertical coordinate not being the actual height of an air parcel but its hydrostatic height (the hypothetical height it would have after the atmospheric column it belongs to has reached hydrostatic balance through adiabatic vertical displacements of air parcels). The Lagrange multiplier λ introduced in Hamilton’s principle to slave density to potential temperature is identified as the nonhydrostatic vertical displacement (i.e., the difference between the actual and hydrostatic heights of an air parcel). The expression of nonhydrostatic pressure and apparent force from λ allow the derivation of a well-defined linear symmetric positive definite problem for λ. As with hydrostatic equations, vertical velocity is diagnosed through Richardson’s equation. The semihydrostatic system has therefore precisely the same degrees of freedom as the hydrostatic primitive equations, while retaining much of the accuracy of the fully compressible Euler equations.

Equations of Atmospheric Motion in Non-Eulerian Vertical Coordinates: Vector-Invariant Form and Quasi-Hamiltonian Formulation

Abstract The curl form of equations of inviscid atmospheric motion in general non-Eulerian coordinates is obtained. Narrowing down to a general vertical coordinate, a quasi-Hamiltonian form is then obtained in a Lagrangian, isentropic, mass-based or z-based vertical coordinate. In non-Lagrangian vertical coordinates, the conservation of energy by the vertical transport terms results from the invariance of energy under the vertical relabeling of fluid parcels. A complete or partial separation between the horizontal and vertical dynamics is achieved, except in the Eulerian case. The horizontal–vertical separation is especially helpful for (quasi-)hydrostatic systems characterized by vanishing vertical momentum. Indeed for such systems vertical momentum balance reduces to a simple statement: total energy is stationary with respect to adiabatic vertical displacements of fluid parcels. From this point of view the purpose of (quasi-)hydrostatic balance is to determine the vertical positions of fluid parcels, for which no evolution equation is readily available. This physically appealing formulation significantly extends previous work. The general formalism is exemplified for the fully compressible Euler equations in a Lagrangian vertical coordinate and a Cartesian (x, z) slice geometry, and the deep-atmosphere quasi-hydrostatic equations in latitude–longitude horizontal coordinates. The latter case, in particular, illuminates how the apparent intricacy of the time-dependent metric terms and of the additional forces can be absorbed into a proper choice of prognostic variables. In both cases it is shown how the quasi-Hamiltonian form leads straightforwardly to the conservation of energy using only integration by parts. Relationships with previous work and implications for stability analysis and the derivation of approximate sets of equations and energy-conserving numerical schemes are discussed.

The influence of winter convection on primary production: A parameterisation using a hydrostatic three-dimensional biogeochemical model

In the recent past observational and modelling studies have shown that the vertical displacement of water parcels, and therefore, phytoplankton particles in regions of deep-reaching convection plays a key role in late winter/early spring primary production. The underlying mechanism describes how convection cells capture living phytoplankton cells and recurrently expose them to sunlight.This study presents a parameterisation called ‘phytoconvection’ which focusses on the influence of convection on primary production. This parameterisation was implemented into a three-dimensional physical–biogeochemical model and applied to the Northwestern European Continental Shelf and areas of the adjacent Northeast Atlantic. The simulation was compared to a ‘conventional’ parameterisation with respect to its influence on phytoplankton concentrations during the annual cycle and its effect on the carbon cycle.The simulation using the new parameterisation showed good agreement with observation data recorded during winter, whereas the reference simulation did not capture the observed phytoplankton concentrations. The new parameterisation had a strong influence on the carbon export through the sinking of particulate organic carbon. The carbon export during late winter/early spring significantly exceeded the export of the reference run.Furthermore, a non-hydrostatic convection model was used to evaluate the major assumption of the presented parameterisation which implies the matching of the mixed layer depth with the convective mixing depth. The applied mixed layer depth criterion principally overestimates the actual convective mixing depth. However, the results showed that this assumption is reasonable during late winter, while indicating a mismatch during spring.

Numerical Modeling of Atmospheric Water Content and Probability Evaluation. Part II

Starting from a numerical simulation of tropical atmosphere in radiative convective equilibrium a perturbation to the surface condition has been applied and the subsequent adjustment of the atmosphere has been analyzed. 3D data have been averaged in space and time using as energy variable to sort them the moist static energy, that is conserved for a moist atmosphere during adiabatic fluid parcel displacements. Since the energy balance of the system corresponds to vertically averaged terms of the MSE differential budget, a entropy derived variable- corresponding to a total thermodynamic energy variable- the θe equivalent potential temperature, has been analyzed using the vertically block averaged column data sorted by MSE. A close correlation between the energy distribution analyzed as quartile and mean for the θe and RH vertical profiles is found, highlighting the dynamical link existing between energy and moisture/heat statistical distribution of the system.