Vertical Velocity Field Research Articles

A fully consistent and conservative vertically adaptive coordinate system for SLIM 3D v0.4 with an application to the thermocline oscillations of Lake Tanganyika

Abstract. The discontinuous Galerkin (DG) finite element method is well suited for the modelling, with a relatively small number of elements, of three-dimensional flows exhibiting strong velocity or density gradients. Its performance can be highly enhanced by having recourse to r-adaptivity. Here, a vertical adaptive mesh method is developed for DG finite elements. This method, originally designed for finite difference schemes, is based on the vertical diffusion of the mesh nodes, with the diffusivity controlled by the density jumps at the mesh element interfaces. The mesh vertical movement is determined by means of a conservative arbitrary Lagrangian–Eulerian (ALE) formulation. Though conservativity is naturally achieved, tracer consistency is obtained by a suitable construction of the mesh vertical velocity field, which is defined in such a way that it is fully compatible with the tracer and continuity equations at a discrete level. The vertically adaptive mesh approach is implemented in the three-dimensional version of the geophysical and environmental flow Second-generation Louvain-la-Neuve Ice-ocean Model (SLIM 3D; www.climate.be/slim). Idealised benchmarks, aimed at simulating the oscillations of a sharp thermocline, are dealt with. Then, the relevance of the vertical adaptivity technique is assessed by simulating thermocline oscillations of Lake Tanganyika. The results are compared to measured vertical profiles of temperature, showing similar stratification and outcropping events.

A New Dynamical Index for India-Burma Trough

Based on the vertical velocity field of reanalysis datasets, this study defines a new dynamical index for the India-Burma trough and supports this index’s advantages by analyzing reanalysis and observational datasets. For a convenient understanding, the vertical velocities of 5 levels ranging from 700 hPa to 500 hPa within the area of 15.625°N–24.375°N and 90.625°E–100.625°E are multiplied by −1 and summed up into a time series involving each year from 1979 to 2012. The standardized value of the time series is defined as the index of India-Burma trough (IIBT). IIBT can reflect the characteristics of the annual strength and the interdecadal variation of the India-Burma trough. IIBT can also well reveal the relationship between the India-Burma trough and its upstream teleconnection. What is more, through a correlation analysis on the grid point precipitation field, respectively, with the IIBT and the India-Burma trough indices defined with vorticity and geopotential height, over southern Asia the correlation pattern between the IIBT and the precipitation field is found to nearly be the sum of the correlation patterns of the latter 2 indices with the precipitation. To the south of the TP, the correlation field between the IIBT and the grid point precipitation shows dipolar distribution, which is consistent with the correlation patterns of the IIBT with the vertical velocity, specific humidity, and the mid-level geopotential height in the same spatial location. IIBT is beneficial for more accurate study of the impact of the India-Burma trough on the associated weather and climate.

Temporal-Spatial Surface Seasonal Mass Changes and Vertical Crustal Deformation in South China Block from GPS and GRACE Measurements.

The solid Earth deforms elastically in response to variations of surface atmosphere, hydrology, and ice/glacier mass loads. Continuous geodetic observations by Global Positioning System (CGPS) stations and Gravity Recovery and Climate Experiment (GRACE) record such deformations to estimate seasonal and secular mass changes. In this paper, we present the seasonal variation of the surface mass changes and the crustal vertical deformation in the South China Block (SCB) identified by GPS and GRACE observations with records spanning from 1999 to 2016. We used 33 CGPS stations to construct a time series of coordinate changes, which are decomposed by empirical orthogonal functions (EOFs) in SCB. The average weighted root-mean-square (WRMS) reduction is 38% when we subtract GRACE-modeled vertical displacements from GPS time series. The first common mode shows clear seasonal changes, indicating seasonal surface mass re-distribution in and around the South China Block. The correlation between GRACE and GPS time series is analyzed which provides a reference for further improvement of the seasonal variation of CGPS time series. The results of the GRACE observations inversion are the surface deformations caused by the surface mass change load at a rate of about −0.4 to −0.8 mm/year, which is used to improve the long-term trend of non-tectonic loads of the GPS vertical velocity field to further explain the crustal tectonic movement in the SCB and surroundings.

Characterizing noise and spurious convection in convective data assimilation

The occurrence of noise and spurious convective cells in analyses of convective‐scale data assimilation (DA) is a problem that is not quantified sufficiently to develop reliable set‐ups for the DA of radar observations in an ensemble Kalman filter. This article presents two approaches to quantifying these phenomena using a local ensemble transform Kalman filter testbed, where simulated radar observations are assimilated with varying spatial and temporal DA parameters. Firstly, the dynamical coupling of cold pools and emerging cells of organized deep convection is quantified in terms of spatio‐temporal correlations, so that spurious perturbations can be identified. Secondly, the abundance of gravity‐wave noise is characterized by the variance of the vertical velocity field during DA cycling and the ensemble forecasts. This evaluation is performed separately for subsets of model points inside observed storms, in their vicinity and in the surrounding clear‐air regions. Varying horizontal localization lengths, observation averaging scales and the length of the assimilation cycling interval are compared, with respect to the two diagnostic approaches. It is concluded that the cold‐pool coupling method and the partitioned variance of vertical velocity fields are viable measures to quantify the influence of DA parameters and algorithms on the balance of analysis states that are used to initiate forecasts of convection.

Kinematics of the tectonic blocks and active faults at the post-orogenic stage: Northern Taiwan

Northern Taiwan offers a rare opportunity to quantify the rate and style of deformation in the transition zone, from the collision in the southeast, subduction, rifting, to orogenic collapsing in the northwest. The updated GPS velocity field collected from 2002 to 2014 and leveling vertical velocity field from 2004 to 2010 were inverted by a 3-D block model in this study to infer the tectonic block rotations and fault slip rates in reference to the Eurasian plate. Our optimized modeling results show that the clockwise rotation rate of 20.8°±2.8°/Myr in the Hualien block is caused by the coupling of the waning arc-continent collision between the Luzon volcanic arc and the Eurasian plate, and the northeastward migration of the block toward the Ryukyu trench. In Ilan domain, the rapid clockwise motion was derived at a rate of 34.9°±65.9°/Myr in southern part of the plain. We proposed that the rapid clockwise rotation results from the asymmetry rifting at the southwestward extension of the Okinawa trough. An another new key finding is highlighted in the capability of maximum magnitude ∼6.2 is indicated on the Suao and the Choshui faults in the Ilan plain.

Geodetic constraint on the motion of a slab window: Implication for the Mendocino Crustal Conveyor model

AbstractThe migration of the slab window in the Northern California Coast Ranges provides a unique setting to study the viscous coupling between crust and asthenosphere flow. The mechanisms of these dynamic processes are explained by the Mendocino Crustal Conveyor model, which predicts a 2‐D “double‐humped” surface uplift rate pattern on a 400 km long profile. To evaluate the Mendocino Crustal Conveyor (MCC) model using accurate geodetic measurements, we derive the vertical velocity field from 43 continuous GPS (CGPS) stations in the Coast Ranges region and project it onto the profile along three possible orientations of the slab window. The CGPS measured uplift rates are in good agreement with the MCC prediction, when the slab window orientation is parallel to the symmetry axis of the region of thickened crust. Thus, the CGPS solutions not only provide a complementary means to diagnose the MCC model but also provide an effective way to constrain the orientation of the slab window.

A Physical Interpretation of the Wind-Wave Instability as Interacting Waves

AbstractOne mechanism for the growth of ocean surface waves by wind is through a shear instability that was first described by Miles in 1957. A physical interpretation of this wind-wave instability is provided in terms of the interaction of the surface gravity wave with perturbations of vorticity within the critical layer—a near-singularity in the airflow where the background flow speed matches that of the surface gravity wave. This physical interpretation relies on the fact that the vertical velocity field is slowly varying across the critical layer, whereas both the displacement and vorticity fields vary rapidly. Realizing this allows for the construction of a physically intuitive description of the critical layer vorticity perturbations that may be approximated by a simple vortex sheet model, the essence of the wind-wave instability can then be captured through the interaction of the critical layer vorticity with the surface gravity wave. This simple model is then extended to account for vorticity perturbations in the airflow profile outside of the critical layer and is found to lead to an exact description of the linear stability problem that is also computationally efficient. The interpretation allows, in general, for the incorporation of sheared critical layers into the “wave interaction theory” that is commonly used to provide a physical description and rationalization of results in the stability of stratified shear flows.

MHD simulations of oscillating cusp-filling tori around neutron stars – missing upper kHz QPO

Abstract We performed axisymmetric, grid-based, ideal magnetohydrodynamic simulations of oscillating cusp-filling tori orbiting a non-rotating neutron star. A pseudo-Newtonian potential was used to construct the constant angular momentum tori in equilibrium. The inner edge of the torus is terminated by a ‘cusp’ in the effective potential. The initial motion of the model tori was perturbed with uniform sub-sonic vertical and diagonal velocity fields. As the configuration evolved in time, we measured the mass accretion rate on the neutron star surface and obtained the power spectrum. The prominent mode of oscillation in the cusp torus is the radial epicyclic mode. It would appear that vertical oscillations are suppressed by the presence of the cusp. From our analysis, it follows that the mass accretion rate carries a modulation imprint of the oscillating torus, and hence so does the boundary layer luminosity.

Ageostrophic Secondary Circulation in a Subtropical Intrathermocline Eddy

AbstractVertical motions play a key role in the enhancement of primary production within mesoscale eddies through the introduction of nutrients into the euphotic layer. However, the details of the vertical velocity field w driving these enhancements remain under discussion. For the first time the mesoscale w associated with an intrathermocline eddy is computed and analyzed using in situ high-resolution three-dimensional (3D) fields of density and horizontal velocity by resolving a generalized omega equation valid for high Rossby numbers. In the seasonal pycnocline the diagnosed w reveals a multipolar structure with upwelling and downwelling cells located at the eddy periphery. In the main pycnocline w is characterized by a dipolar structure with downwelling velocities upstream of the propagation path and upwelling velocities downstream. Maximum values of w reach 6.4 m day−1. An observed enhancement of chlorophyll-a at the eddy periphery coincides with the location of the upwelling and downwelling cells. Analysis of the forcing terms of the generalized omega equation indicates that the mechanisms behind the dipolar structure of the w field are a combination of horizontal deformation and advection of vertical relative vorticity by ageostrophic vertical shear. The wind during the eddy sampling was rather constant and uniform with a speed of 5 m s−1. Diagnosed nonlinear Ekman pumping leads to a dipolar pattern that mirrors the inferred w. Horizontal ageostrophic secondary circulation is dominated by centripetal acceleration and closes the dipole w structure. Vertical fluxes act to maintain the intrathermocline eddy structure.

Convoluted Intake Distortion Measurements Using Stereo Particle Image Velocimetry

The unsteady distorted flowfields generated within convoluted aeroengine intakes can compromise the engine performance and operability. Therefore, there is a need for a better understanding of the complex characteristics of the distorted flow at the exit of S-shaped intakes. This work presents a detailed analysis of the unsteady swirl distortion based on synchronous, high-spatial-resolution measurements using stereoscopic particle image velocimetry. Two S-duct configurations with different centerline offsets are investigated. The high-offset duct shows greater levels of dynamic and steady swirl distortion and a notably greater tendency toward bulk swirl patterns associated with high swirl distortion. More discrete distortion patterns with locally high swirl levels and the potential to impact the engine operability are identified. The most energetic coherent structures of the flowfield are observed using proper orthogonal decomposition. A switching mode is identified that promotes the alternating swirl switching mechanism and is mostly associated with the occurrence of potent bulk swirl events. A vertical mode that characterizes a perturbation of the vertical velocity field promotes most of the twin swirl flow distortion topologies. It is postulated that it is associated with the unsteadiness of the centerline shear layer.

Three-dimensional study of parallel shear flow around an obstacle in water channel and air tunnel

The purpose of this study is to show the effect of some physical parameters on a parallel shear flow. The variation of the velocity inlet and the influence of an obstacle placed at the bottom of a channel were analyzed. The governing equations based on K –e model in a line source downstream around a three-dimensional obstacle are determined by the finite volume method with SIMPLEC algorithm. The obtained results have allowed to establish the dynamic characteristics of this kind of flow. Horizontal and vertical velocity fields, Reynolds shear stress field, turbulence intensity and scalar concentration in the water channel and a wind tunnel were exposed. These results permit us to localize the maximum and minimum zones of turbulence, which help us to exploit those fields for different purposes such as energy extraction and port infrastructure.

水平来流对扰动成长和对流周期性的影响

Numerical simulation of the 2D fully hydrodynamic equations for the pure fluid in a rectangular channel with horizontal flow for Prandtl number Pr=0.0272 was conducted. Growth and spatiotemporal evolution of the 1D traveling wave patterns in the RayleighBenard convection of the pure fluid were investigated. It is found that the convective growth process can be divided into 3 stages: the development stage, the exponential growth stage and the periodic variation stage. Through analysis on the variation of the maximum vertical velocity field with time for different relative Rayleigh numbers Ra_r in the exponential growth stage, a formula of variation of linear growth rate γ_m was obtained with respect to Ra_r. Furthermore, the traveling wave convection periodicity and its dependence on the horizontal flow Reynolds number were revealed.

Experimental study of flow field characteristics on bed configurations in the pebble bed reactor

The flow field characteristics are of fundamental importance in the design work of the pebble bed high temperature gas cooled reactor (HTGR). The different effects of bed configurations on the flow characteristics of pebble bed are studied through the PTV (Particle Tracking Velocimetry) experiment. Some criteria, e.g. flow uniformity (σ) and mass flow level (α), are proposed to estimate vertical velocity field and compare the bed configurations. The distribution of the Δθ (angle difference between the individual particle velocity and the velocity vector sum of all particles) is also used to estimate the resultant motion consistency level. Moreover, for each bed configuration, the thickness of displacement is analyzed to measure the effect of the funnel flow zone based on the boundary layer theory. Detailed information shows the quantified characteristics of bed configuration effects on flow uniformity and other characteristics; and the sequence of levels of each estimation criterion is obtained for all bed configurations. In addition, a good design of the pebble bed configuration is suggested and these estimation criteria can be also applied and adopted in testing other geometry designs of pebble bed.

POLARIZED LINE FORMATION IN NON-MONOTONIC VELOCITY FIELDS

ABSTRACT For a correct interpretation of the observed spectro-polarimetric data from astrophysical objects such as the Sun, it is necessary to solve the polarized line transfer problems taking into account a realistic temperature structure, the dynamical state of the atmosphere, a realistic scattering mechanism (namely, the partial frequency redistribution—PRD), and the magnetic fields. In a recent paper, we studied the effects of monotonic vertical velocity fields on linearly polarized line profiles formed in isothermal atmospheres with and without magnetic fields. However, in general the velocity fields that prevail in dynamical atmospheres of astrophysical objects are non-monotonic. Stellar atmospheres with shocks, multi-component supernova atmospheres, and various kinds of wave motions in solar and stellar atmospheres are examples of non-monotonic velocity fields. Here we present studies on the effect of non-relativistic non-monotonic vertical velocity fields on the linearly polarized line profiles formed in semi-empirical atmospheres. We consider a two-level atom model and PRD scattering mechanism. We solve the polarized transfer equation in the comoving frame (CMF) of the fluid using a polarized accelerated lambda iteration method that has been appropriately modified for the problem at hand. We present numerical tests to validate the CMF method and also discuss the accuracy and numerical instabilities associated with it.

On the Nature of the Transition Between Roll and Cellular Organization in the Convective Boundary Layer

Both observational and numerical studies of the convective boundary layer (CBL) have demonstrated that when surface heat fluxes are small and mean wind shear is strong, convective updrafts tend to organize into horizontal rolls aligned within 10–20\(^\circ \) of the geostrophic wind direction. However, under large surface heat fluxes and weak to negligible shear, convection tends to organize into open cells, similar to turbulent Rayleigh-Benard convection. Using a suite of 14 large-eddy simulations (LES) spanning a range of \(-z_i/L\) between zero (neutral) and 1041 (highly convective), where \(z_i\) is the CBL depth and L is the Obukhov length, the transition between roll- and cellular-type convection is investigated systematically for the first time using LES. Mean vertical profiles including velocity variances and turbulent transport efficiencies, as well the “roll factor,” which characterizes the rotational symmetry of the vertical velocity field, indicate the transition occurs gradually over a range of \(-z_i/L\); however, the most significant changes in vertical profiles and CBL organization occur from near-neutral conditions up to about \(-z_i/L \approx \) 15–20. Turbulent transport efficiencies and quadrant analysis are used to characterize the turbulent transport of momentum and heat with increasing \(-z_i/L\). It is found that turbulence transports heat efficiently from weakly to highly convective conditions; however, turbulent momentum transport becomes increasingly inefficient as \(-z_i/L\) increases.

Evaluating a campaign GNSS velocity field derived from an online precise point positioning service

Traditional processing of Global Navigation Satellite System (GNSS) data using dedicated scientific software has provided the highest levels of positional accuracy, and has been used extensively in geophysical deformation studies. To achieve these accuracies a significant level of understanding and training is required, limiting their availability to the general scientific community. Various online GNSS processing services, now freely available, address some of these difficulties and allow users to easily process their own GNSS data and potentially obtain high quality results. Previous research into these services has focused on Continually Operating Reference Station (CORS) GNSS data. Less research exists on the results achievable with these services using large campaign GNSS data sets, which are inherently noisier than CORS data. Even less research exists on the quality of velocity fields derived from campaign GNSS data processed through online precise point positioning services. Particularly, whether they are suitable for geodynamic and deformation studies where precise and reliable velocities are needed. In this research, we process a very large campaign GPS data set (spanning 10 yr) with the online Jet Propulsion Laboratory Automated Precise Positioning Service. This data set is taken from a GNSS network specifically designed and surveyed to measure deformation through the central North Island of New Zealand. This includes regional CORS stations. We then use these coordinates to derive a horizontal and vertical velocity field. This is the first time that a large campaign GPS data set has been processed solely using an online service and the solutions used to determine a horizontal and vertical velocity field. We compared this velocity field to that of another well utilized GNSS scientific software package. The results show a good agreement between the CORS positions and campaign station velocities obtained from the two approaches. We discuss the implications of these results for how future GNSS campaign field surveys might be conducted and how their data might be processed

GPS Imaging of vertical land motion in California and Nevada: Implications for Sierra Nevada uplift.

We introduce Global Positioning System (GPS) Imaging, a new technique for robust estimation of the vertical velocity field of the Earth's surface, and apply it to the Sierra Nevada Mountain range in the western United States. Starting with vertical position time series from Global Positioning System (GPS) stations, we first estimate vertical velocities using the MIDAS robust trend estimator, which is insensitive to undocumented steps, outliers, seasonality, and heteroscedasticity. Using the Delaunay triangulation of station locations, we then apply a weighted median spatial filter to remove velocity outliers and enhance signals common to multiple stations. Finally, we interpolate the data using weighted median estimation on a grid. The resulting velocity field is temporally and spatially robust and edges in the field remain sharp. Results from data spanning 5–20 years show that the Sierra Nevada is the most rapid and extensive uplift feature in the western United States, rising up to 2 mm/yr along most of the range. The uplift is juxtaposed against domains of subsidence attributable to groundwater withdrawal in California's Central Valley. The uplift boundary is consistently stationary, although uplift is faster over the 2011–2016 period of drought. Uplift patterns are consistent with groundwater extraction and concomitant elastic bedrock uplift, plus slower background tectonic uplift. A discontinuity in the velocity field across the southeastern edge of the Sierra Nevada reveals a contrast in lithospheric strength, suggesting a relationship between late Cenozoic uplift of the southern Sierra Nevada and evolution of the southern Walker Lane.

ENSO regulation of far‐ and mid‐infrared contributions to clear‐sky OLR

AbstractNASA Aqua‐derived thermodynamic profiles, calculated spectral clear‐sky outgoing longwave radiation (OLR), and vertical velocity fields from meteorological reanalyses are combined to determine the relative proportion of the far‐infrared (FIR) and mid‐infrared (MIR) spectral contributions to the total clear‐sky OLR during different phases of El Niño–Southern Oscillation (ENSO). In the ascending branch of the tropical circulation, the spatial variance of upper tropospheric water vapor is shown to be larger during La Niña than El Niño and is consistent with zonal symmetry changes in the tropical waveguide and associated tropical‐extratropical mixing. In the descending branch, upper tropospheric water vapor shows weaker coupling to lower layers that is evidenced by changes in the ratio of FIR to MIR in the clear‐sky OLR. Diagnostics from the Geophysical Fluid Dynamics Laboratory AM3 model simulation are generally similar to satellite data, but the ratio of FIR to MIR is 5–10% larger with respect to dynamic regime.

Erratum to: Inhibition of mixed-layer deepening during winter in the northeastern Arabian Sea by the West India Coastal Current

Though the deep mixed layers (MLs) that form in the northeastern Arabian Sea (NEAS) during the winter monsoon (November–February) have been attributed to convective mixing driven by dry, cool northeasterly winds from the Indian subcontinent, data show that the deepest MLs occur in the northern NEAS and the maxima of latent-heat and net heat fluxes in the southern NEAS. We use an oceanic general circulation model to show that the deep MLs in the NEAS extend up to ~20°N till the end of December, but are restricted poleward of ~22°N (~23°N) in January (February). This progressive restriction of the deep mixed layers within the NEAS is due to poleward advection of water of lower salinity by the West India Coastal Current (WICC). The deep MLs are sustained till February in the northern NEAS because convective mixing deepens the ML before the waters of lower salinity reach this region and the wind stirring and convective overturning generate sufficient turbulent energy for the ML to maintain the depth attained in January. Though the atmospheric fluxes tend to cool the ML in the southern NEAS, this cooling is countered by the warming due to horizontal advection. Likewise, the cooling due to entrainment, which continues in the southern NEAS even as the ML shallows during January–February, is almost cancelled by the warming caused by a downwelling vertical velocity field. Therefore, the SST changes very little during December–February even as the ML shallows dramatically in the southern NEAS. These deep MLs of the NEAS also preclude a strong intraseasonal response to the intraseasonal variability in the fluxes. This role of horizontal advection implies that the ML depth in the NEAS is determined by an interplay of physical processes that are forced differently. The convective mixing depends on processes that are local to the region, but the advection is due to the WICC, whose seasonal cycle is primarily forced by remote winds. By inhibiting the formation of deep MLs in the southern NEAS, the WICC limits the region of formation of the high-salinity water masses of this region. Since the deep MLs in the NEAS have been linked to the high chlorophyll concentration there, our results imply that the conventional approach of averaging over boxes for studying the impact of physics on biogeochemistry can mask important details that are due to advection because it is the advective component of any budget that is most affected by the averaging process.

Dust devil dynamics

AbstractA self‐consistent hydrodynamic model for the solar heating‐driven onset of a dust devil vortex is derived and analyzed. The toroidal flows and vertical velocity fields are driven by an instability that arises from the inversion of the mass density stratification produced by solar heating of the sandy surface soil. The nonlinear dynamics in the primary temperature gradient‐driven vertical airflows drives a secondary toroidal vortex flow through a parametric interaction in the nonlinear structures. While an external tangential shear flow may initiate energy transfer to the toroidal vortex flow, the nonlinear interactions dominate the transfer of vertical‐radial flows into a fast toroidal flow. This secondary flow has a vertical vorticity, while the primary thermal gradient‐driven flow produces the toroidal vorticity. Simulations for the complex nonlinear structure are carried out with the passive convection of sand as test particles. Triboelectric charging modeling of the dust is used to estimate the charging of the sand particles. Parameters for a Dust Devil laboratory experiment are proposed considering various working gases and dust particle parameters. The nonlinear dynamics of the toroidal flow driven by the temperature gradient is of generic interest for both neutral gases and plasmas.