Low Vertical Resolution Research Articles

A Case Study: Evaluation of PAFOG One‐D Model With Advection in Simulations of Fog/Stratus From C‐FOG Experiment

AbstractFog simulation is a challenge due to the complex microphysical, radiative, and turbulent processes. C‐FOG (Coastal Fog) was a comprehensive experiment which aimed to improve our understanding and forecasting skill of fog in coastal areas of Atlantic Canada. An intermittent fog event, containing two fog periods occurring over the Grand Banks of Newfoundland beginning at 00:00 UTC on September 13, 2018, is investigated in this study. Our approach is to apply the Weather Research and Forecasting (WRF) model, in conjunction with a 1‐D model for the atmospheric boundary layer, PAFOG (Parameterized Fog model). Model results are evaluated against observations collected onboard the research vessel Hugh R. Sharp. We introduced an advection term to PAFOG to investigate the ability of the model to simulate the two‐level fog/stratus system and we developed a set of sensitivity experiments. Results suggest that PAFOG with the advection terms calculated from ERA5 shows a good ability to simulate intermittent fog caused by fog lifting as stratus, and by comparison, outperform results with advection terms derived from WRF simulations. Regarding the latter, in an experiment with high vertical resolution, WRF gives a delayed fog event, whereas with low vertical resolution, WRF gives descending stratus. Our results provide a possible methodology to study the detailed structure of fog/stratus using PAFOG, a sounding profile as initial conditions, and advection profiles derived from ERA5, with low computing cost.

Vertical grid refinement for stratocumulus clouds in the radiation scheme of the global climate model ECHAM6.3-HAM2.3-P3

Abstract. In this study, we implement a vertical grid refinement scheme in the radiation routine of the global aerosol–climate model ECHAM-HAM, aiming to improve the representation of stratocumulus clouds and address the underestimation of their cloud cover. The scheme is based on a reconstruction of the temperature inversion as a physical constraint for the cloud top. On the refined grid, the boundary layer and the free troposphere are separated and the cloud's layer is made thinner. The cloud cover is recalculated either by conserving the cloud volume (SC-VOLUME) or by using the Sundqvist cloud cover routine on the new grid representation (SC-SUND). In global climate simulations, we find that the SC-VOLUME approach is inadequate, as there is a mismatch, in most cases, between the layer of the inversion and the layer of the stratocumulus cloud, which prevents its application and is itself likely caused by an overly low vertical resolution. Additionally, we find that the occurrence frequency of stratocumulus clouds is underestimated in ECHAM-HAM, limiting a priori the potential benefits of a scheme like SC-VOLUME targeting only cloud amount when present. With the SC-SUND approach, the possibility for new clouds to be formed on the refined grid results in a large increase in mean total cloud cover in stratocumulus regions. In both cases, however, the changes exerted in the radiation routine are too weak to produce a significant improvement in the simulated stratocumulus cloud cover. We investigate and discuss the reasons behind this. The grid refinement scheme could be used more effectively for this purpose if implemented directly in the model's cloud microphysics and cloud cover routines, but other possible ways forward are also discussed.

Comprehensive quantification of height dependence of entrainment mixing between stratiform cloud top and environment

Abstract. Different entrainment-mixing processes of turbulence are crucial to processes related to clouds; however, only a few qualitative studies have been concentrated on the vertical distributions of entrainment-mixing mechanisms with low vertical resolutions. To quantitatively study vertical profiles of entrainment-mixing mechanisms with a high resolution, the stratiform clouds observed in the Physics of Stratocumulus Top (POST) project are examined. The unique sawtooth flight pattern allows for an examination of the vertical distributions of entrainment-mixing mechanisms with a 5 m vertical resolution. Relative standard deviation of volume mean radius divided by relative standard deviation of liquid water content is introduced to be a new estimation of microphysical homogeneous mixing degree, to overcome difficulties of determining the adiabatic microphysical properties required in existing measures. The vertical profile of this new measure indicates that entrainment-mixing mechanisms become more homogeneous with decreasing altitudes and are consistent with the dynamical measures of Damköhler number and transition scale number. Further analysis shows that the vertical variation of entrainment-mixing mechanisms with decreasing altitudes is due to the increases of turbulent dissipation rate in cloud and relative humidity in droplet-free air and the decrease of size of droplet-free air. The results offer insights into the theoretical understanding and parameterizations of vertical variation of entrainment-mixing mechanisms.

Flow under long two-dimensional dam sluice gate using WSPH

An application of the weakly compressible Smoothed Particle Hydrodynamics (WSPH) numerical method is presented here for the case of two-dimensional flow in a long channel with a partially open sluice gate. The results are compared with an analytical solution provided by shallow water equations (SWE) and available experimental data. Of particular interest is the application of this numerical method to a sluice gate case with a high ratio of channel length to depth, which tends to amplify the effects of the chosen numerical resolution. Good model congruence was observed even for relatively low vertical resolution, and the effects of the equations used to describe the boundary conditions were identified.

Modeling the GABLS4 Strongly‐Stable Boundary Layer With a GCM Turbulence Parameterization: Parametric Sensitivity or Intrinsic Limits?

AbstractThe representation of stable boundary layers (SBLs) still challenges turbulence parameterizations implemented in current weather or climate models. The present work assesses whether these model deficiencies reflect calibration choices or intrinsic limits in currently‐used turbulence parameterization formulations and implementations. This question is addressed for the CNRM atmospheric model ARPEGE‐Climat 6.3 in a single‐column model/large‐eddy simulation (SCM/LES) comparison framework, using the history matching with iterative refocusing statistical approach. The GABLS4 case, which samples a nocturnal strong SBL observed at Dome C, Antarctic Plateau, is used. The standard calibration of the ARPEGE‐Climat 6.3 turbulence parameterization leads to a too deep SBL, a too high low‐level jet and misses the nocturnal wind rotation. This behavior is found for low and high vertical resolution model configurations. The statistical tool then proves that these model deficiencies reflect a poor parameterization calibration rather than intrinsic limits of the parameterization formulation itself. In particular, the role of two lower bounds that were heuristically introduced during the parameterization implementation to increase mixing in the free troposphere and to avoid runaway cooling in snow‐ or ice‐covered region is emphasized. The statistical tool identifies the space of the parameterization free parameters compatible with the LES reference, accounting for the various sources of uncertainty. This space is non‐empty, thus proving that the ARPEGE‐Climat 6.3 turbulence parameterization contains the required physics to capture the GABLS4 SBL. The SCM framework is also used to validate the statistical framework and a few guidelines for its use in parameterization development and calibration are discussed.

The Topographical Optimization of 3D Microgroove Pattern Intervals for Ligamentous Cell Orientations: In Vitro.

Specific orientations of periodontal ligaments (PDLs) to tooth-root surface play an important role in offering positional stabilities of teeth, transmitting and absorbing various stresses under masticatory/occlusal loading conditions, or promoting tissue remodeling by mechanical stimulations to periodontal cells. However, it is still challenging to spatially control PDL orientations and collective PDL cell alignments using 3D scaffold architectures. Here, we investigated the optimization of scaffold topographies in order to control orientations of human PDL cells with predictability in in vitro. The 3D PDL-guiding architectures were designed by computer-aided design (CAD) and microgroove patterns on the scaffold surfaces were created with four different slice intervals such as 25.40 µm (μG-25), 19.05 µm (μG-19), 12.70 µm (μG-12), and 6.35 µm (μG-6) by the digital slicing step. After scaffold design and 3D wax printing, poly-ε-caprolactone (PCL) was casted into 3D printed molds and human PDL cells were cultured for 7 days. In the results, μG-25 with low vertical resolution can angularly organize seeded cells predictably rather than μG-6 created by the highest resolution for high surface quality (or smooth surface). Moreover, nuclear orientations and deformability were quantitatively analyzed and a significant correlation between microgroove pattern intervals and cell alignments was calculated for the topographic optimization. In conclusion, controllable microgroove intervals can specifically organize human PDL cells by 3D printing, which can create various surface topographies with structural consistence. The optimal surface topography (μG-25) can angularly guide human PDL cells, but 6.35 µm-thick patterns (μG-6) showed random organization of cell collectivity.

The impact of vertical resolution in reducing biases in sea surface temperature in a tropical Pacific Ocean model

Motivated by observations of fine scale vertical shear and its contribution to mixing in the tropical ocean, this study explores the impact of vertical resolution in an ocean model on sea surface temperature in the tropical Pacific Ocean. We conduct two model experiments that differ in the vertical discretization only, with the grid spacing in one being significantly smaller than the other in the upper ocean. We examine the temperature difference between the high and low vertical resolution experiments. We find that the difference in the upper-most layer is positive in the equatorial cold tongue and negative along the South American coast, thus reducing the commonly seen cool and warm biases in the two regions, respectively. The change in the structure of the vertical diffusivity, as determined by the K-profile parameterization, is identified as the primary cause in reducing the biases. In the central equatorial Pacific, the change in the vertical diffusivity from low to high vertical resolution in the upper pycnocline results in a positive temperature difference that propagates eastward as an equatorial Kelvin wave, rising to the sea surface in the central and eastern regions to increase the sea surface temperature there. In the far eastern equatorial Pacific, the change in the vertical diffusivity in the lower pycnocline produces a negative temperature difference that propagates poleward as coastal Kelvin waves along the west coast of the American continent, outcropping along the South American coast to reduce the sea surface temperature there. The high vertical resolution experiment captures much of the small-scale vertical velocity shear and resolves the fine details of the stratification in the upper ocean. Our analysis suggests that the shear-generated turbulence is the primary contributor to the change in the vertical diffusivity in the central region whereas stratification is the dominant factor in the far eastern region.

Resolving the thickness of peat deposits with contact-less electromagnetic methods: A case study in the Venice coastland

Peat soils are typical deposits characterizing wetlands and reclaimed farmlands. They are important carbon reservoirs and when degraded (e.g., erosive processes, fires, draining and plowing) massive carbon dioxide volumes are released. This leads to increase greenhouse effect and induce serious land subsidence. Thus, mapping the volume of peat deposits is crucial in order to estimate the carbon mass and the potential release of carbon dioxide and consequent loss in soil elevation. Despite the importance of such estimations, forecasting and quantifying the peat thickness is still a challenge. Direct sediment coring provides local information that is difficult to extend to large territories. Indirect geophysical methods are unable to resolve lithological contrasts in the presence of saltwater contamination in coastal areas. In this work, we show the results obtained using two contact-less electromagnetic methods for the characterization of peat deposits in a peatland site of the Venice coastland, Italy. Specifically, a multi-frequency portable instrument (FDEM) and an airborne time-domain electromagnetic one (AEM), known for their very high and relatively low vertical resolution respectively, were used to collect data over a former wetland then reclaimed for agricultural purposes. Additional electrical resistivity tomography (ERT) data are used together with sediment core data to assess the effectiveness and accuracy of the contact-less methods. Results show that both FDEM and AEM are very effective in detecting the presence of the peat layer, despite its low thickness (<2 m) and the high electro-conductive subsoil because of saltwater contamination. However, the AEM method overestimated the peat thickness while the FDEM could accurately resolve the peat thickness even where the layer was thinner than 1 m. When compared to the electrical features extracted from the ERT, discrepancies are on average lower than 30%; when compared to the borehole data, discrepancies are on average slightly higher than 6%.

An iterative inversion of Dual Induction Tool logs from thin-bedded sandy\u2013shaly formations of the Carpathian Foredeep using a modified simulated annealing method

In thin-bedded sandy–shaly Miocene formations of the Carpathian Foredeep, the main source of errors in gas saturation evaluation is the underestimation of resistivity of thin, hydrocarbon-bearing beds, which is the result of the low vertical resolution of induction logging tools. This problem is especially visible in older boreholes drilled in times where the Dual Induction Tool (DIT) was the primary induction tool used for determining the formation resistivity, and in shallowest depth intervals of newer boreholes where the DIT was used instead of newer array tools for cost-saving reasons. In this paper, we show how a global inversion algorithm was used to improve the vertical resolution of DIT logs. Our implementation of an iterative inversion utilizes a one-dimensional formation model, vertical response functions of the DIT, and a modified simulated annealing algorithm to determine the true vertical distribution of the formation resistivity. The algorithm was tested on resistivity logs recorded in a borehole drilled in the Carpathian Foredeep in Poland, where the DIT and the High-Resolution Array Induction (HRAI) tool were run in the same depth interval.

Vertical Structure and Seasonal Features of the Heat Island and Humidity Distribution over Moscow Derived from Satellite Data

The vertical structure and annual variations in the urban heat island (UHI) and dry/wet islands over Moscow are investigated using MetOp satellite data for 2015–2017. It is shown that Moscow is characterized by the dominance of the evening maximum of UHI intensity; however, the UHI intensification during the cold season mainly occurs in the daytime. According to satellite data, the maximum development of UHI over Moscow is observed in the forenoon in January. The urban effect propagates to the height of ~2000 m during the warm season and to the height of 1250 m in the forenoon and to 800 m in the evening during the cold season. During the warm season, UHI is accompanied by the dry island in the forenoon and by the wet island at the upper levels in the evening. The existence of the wet island in the lower air layers was not corroborated. Due to the current accuracy of temperature and humidity retrieval (0.5–1 K and 10%, respectively) as well as due to the low vertical resolution of satellite data, the present research is mainly methodological, and the results are preliminary and open to discussion.

Sediment classification in a Brazilian reservoir: Pros and cons of parametric low frequencies

Sediment is the main factor that limits the reservoir lifetime. Therefore, sediment classification is an essential tool for planning and operating reservoir management measures. There has been important development in the hydroacoustic classification of lakebed, especially with linear systems. The main restrictions while using linear hydroacoustic systems for lakebed classification are the shallow penetration in high-frequency applications or the low vertical and horizontal resolution when using low frequencies. With the new developments in the area of echo sounders, parametric systems can achieve high penetration while preserving the high vertical and lateral resolution. To investigate the performance of parametric systems, a new lakebed classification approach was implemented by using a SES2000 Compact. The area studied was the Passauna reservoir in Parana State, Brazil. We used the first echo division method for processing the acoustic data combined with sediment core and grab sampling. The two physical parameters investigated, were the share of the finest fraction (&lt;63 µm) and wet bulk density (WBD). The results showed a high correlation between the primary frequency of 100 kHz (166 µs pulse length) and the physical parameters. Additionally, a significant correlation was observed with the acoustic parameters at 10 kHz frequency. The best correlating acoustic parameter was Attack/Decay (E1´/E1). The gas presence was found to be an important factor determining the penetration depth of the parametric system and the performance of the classification. The advantages of parametric systems, such small directivity and layering effect, represent the major restrictions in sediment classification applications.

A Decision Tree Based Road Recognition Approach Using Roadside Fixed 3D LiDAR Sensors

As one of the most important elements in the intelligent transportation system (ITS), the road traffic monitoring system (RTMS) needs to be functioned with a road recognition mechanism. Current works on road recognition mainly target at the field of automatic driving and cannot be directly used in the RTMS. In this paper, we propose a decision tree-based road recognition algorithm using roadside fixed light detection and ranging (LiDAR) sensors in the RTMS. These LiDAR sensors have a low vertical resolution, which implies that we cannot get a clear far boundary and obvious features of roads from the point cloud data. Point cloud data obtained by the roadside LiDAR sensors are projected onto a plane rasterized to grids of points. Using a decision tree, these grids are first classified into background grids and road grids. For reducing misclassification, these grids are further reclassified using a mean filtering algorithm. Finally, a minimum circumscribed rectangle algorithm is employed to obtain accurate road boundaries. The experiment results show that compared to existing road recognition algorithms, the proposed approach has advantages of being completely automatic, requiring shorter recognition time and having a wider detection range.

Focus on software, data acquisition, and instrumentation

Focus on software, data acquisition, and instrumentation

The Simulation of Stratospheric Water Vapor over the Asian Summer Monsoon Region in CESM1(WACCM) Models.

Previous observational studies have found a persistent maximum in stratospheric water vapor (SWV) in the upper troposphere lower stratosphere (UTLS) confined by the upper-level anticyclone over the Asian summer monsoon region. This study investigates the simulation of SWV in the Whole Atmosphere Community Climate Model (WACCM). WACCM generally tends to simulate a SWV maximum over the central Pacific Ocean, but this bias is largely improved in the high vertical resolution version. The high vertical resolution model with increased vertical layers in the UTLS is found to have a less stratified UTLS over the central Pacific Ocean compared with the low vertical resolution model. It therefore simulates a steepened PV gradient over the central Pacific Ocean that better closes the upper-level anticyclone and confines the SWV within the enhanced transport barrier.

Regularized inversion of amplitude-versus-incidence angle (AVA) based on a piecewise-smooth model

Different from the stacked seismic data, pre-stack data includes abundant information about shear wave and density. Through inversing the shear wave and density information from the pre-stack data, we can determine oil-bearing properties from different incident angles. The state-of-the-art inversion methods obtain either low vertical resolution or lateral discontinuities. However, the practical reservoir generally has sharp discontinuities between different layers in vertically direction and is horizontally smooth. Towards obtaining the practical model, we present an inversion method based on the regularized amplitude-versus-incidence angle (AVA) data to estimate the piecewise-smooth model from pre-stack seismic data. This method considers subsurface stratum as a combination of two parts: a piecewise smooth part and a constant part. To fix the ill-posedness in the inversion, we adopt four terms to define the AVA inversion misfit function: the data misfit itself, a total variation regularization term acting as a sparsing operator for the piecewise constant part, a Tikhonov regularization term acting as a smoothing operator for the smooth part, and the last term to smoothly incorporate a priori information for constraining the magnitude of the estimated model. The proposed method not only can incorporate structure information and a priori model constraint, but also is able to derive into a convex objective function that can be easily minimized using iterative approach. Compared with inversion results of TV and Tikhonov regularization methods, the inverted P-wave velocity, S-wave velocity and density of the proposed method can better delineate the piecewise-smooth characteristic of strata.

Obstacle Prediction for Automated Guided Vehicles Based on Point Clouds Measured by a Tilted LIDAR Sensor

Environment analysis of automatic vehicles needs the detection from 3-D point cloud information. This paper addresses this task when only partial scanning data are available. Our method develops the detection capabilities of autonomous vehicles equipped with 3-D range sensors for navigation purposes. In industrial practice, the safety scanners of automated guided vehicles (AGVs) and a localization technology provide an additional possibility to gain 3-D point clouds from planar contour points or low vertical resolution. Based on this data and a suitable evaluation algorithm, intelligence of vehicles can be significantly increased without the need for installation of additional sensors. In this paper, we propose a solution for an obstacle categorization problem for partial point clouds without shape modeling. The approach is tested for a known database, as well as for real-life scenarios. In case of AGVs, real-time run is provided by on-board computers of usual complexity.

Correction of seismic time-lapse data with different geometries

ABSTRACTIn the interpretation of time-lapse seismic data, we first assume that basal and monitor seismic data have been gathered and processed under equal condition. Basically having two sets of seismic data (basal and monitor data) does not confirm the possibility of 4D seismic data analysis and repeatability term commands that data should have equal data gathering and process condition. This is due to the low vertical resolution of seismic data which makes the changes in amplitude of seismic in basal and monitoring operations to occur finitely. Therefore, differences in amplitudes in two seismic surveys that have a source other than variations in reserval fluid saturation cause errors in 4D interpretations. In this paper, we represent a procedure to unify different geometries and also study post-stack processing methods for seismic time-lapse data and execute them on existing data at the studied field. In the, first step geometrical conception of basal and monitoring data should be equalized. This section has been implemented using a new methodology, developed in this paper and is called “similarity table”. Then, data will be analyzed with spectra point of view and by executing shaping filters and by using cross equalization approach, basal and monitor seismic information will be closer as far as possible. Finally, analyzes were performed to validate the obtained results and evaluate the quality of time -lapse seismic analyzes.

Lagrangian condensation microphysics with Twomey CCN activation

Abstract. We report the development of a novel Lagrangian microphysics methodology for simulations of warm ice-free clouds. The approach applies the traditional Eulerian method for the momentum and continuous thermodynamic fields such as the temperature and water vapor mixing ratio, and uses Lagrangian super-droplets to represent condensed phase such as cloud droplets and drizzle or rain drops. In other applications of the Lagrangian warm-rain microphysics, the super-droplets outside clouds represent unactivated cloud condensation nuclei (CCN) that become activated upon entering a cloud and can further grow through diffusional and collisional processes. The original methodology allows for the detailed study of not only effects of CCN on cloud microphysics and dynamics, but also CCN processing by a cloud. However, when cloud processing is not of interest, a simpler and computationally more efficient approach can be used with super-droplets forming only when CCN is activated and no super-droplet existing outside a cloud. This is possible by applying the Twomey activation scheme where the local supersaturation dictates the concentration of cloud droplets that need to be present inside a cloudy volume, as typically used in Eulerian bin microphysics schemes. Since a cloud volume is a small fraction of the computational domain volume, the Twomey super-droplets provide significant computational advantage when compared to the original super-droplet methodology. Additional advantage comes from significantly longer time steps that can be used when modeling of CCN deliquescence is avoided. Moreover, other formulation of the droplet activation can be applied in case of low vertical resolution of the host model, for instance, linking the concentration of activated cloud droplets to the local updraft speed. This paper discusses the development and testing of the Twomey super-droplet methodology, focusing on the activation and diffusional growth. Details of the activation implementation, transport of super-droplets in the physical space, and the coupling between super-droplets and the Eulerian temperature and water vapor field are discussed in detail. Some of these are relevant to the original super-droplet methodology as well and to the ice phase modeling using the Lagrangian approach. As a computational example, the scheme is applied to an idealized moist thermal rising in a stratified environment, with the original super-droplet methodology providing a benchmark to which the new scheme is compared.

Changing transport processes in the stratosphere by radiative heating of sulfate aerosols

Abstract. The injection of sulfur dioxide (SO2) into the stratosphere to form an artificial stratospheric aerosol layer is discussed as an option for solar radiation management. Sulfate aerosol scatters solar radiation and absorbs infrared radiation, which warms the stratospheric sulfur layer. Simulations with the general circulation model ECHAM5-HAM, including aerosol microphysics, show consequences of this warming, including changes of the quasi-biennial oscillation (QBO) in the tropics. The QBO slows down after an injection of 4 Tg(S) yr−1 and completely shuts down after an injection of 8 Tg(S) yr−1. Transport of species in the tropics and sub-tropics depends on the phase of the QBO. Consequently, the heated aerosol layer not only impacts the oscillation of the QBO but also the meridional transport of the sulfate aerosols. The stronger the injection, the stronger the heating and the simulated impact on the QBO and equatorial wind systems. With increasing injection rate the velocity of the equatorial jet streams increases, and the less sulfate is transported out of the tropics. This reduces the global distribution of sulfate and decreases the radiative forcing efficiency of the aerosol layer by 10 to 14 % compared to simulations with low vertical resolution and without generated QBO. Increasing the height of the injection increases the radiative forcing only for injection rates below 10 Tg(S) yr−1 (8–18 %), a much smaller value than the 50 % calculated previously. Stronger injection rates at higher levels even result in smaller forcing than the injections at lower levels.

Single-shot grating-based x-ray differential phase contrast imaging with a modified analyzer grating**Project supported by the National Natural Science Foundation of China (Grant Nos. 11275204, 11475175, and 11405175), the China Postdoctoral Science Foundation (Grant No. 2017M612097), and the Fundamental Research Funds for the Central Universities (Grant No. WK2310000065)

X-ray grating interferometer has attracted widely attention in the past years due to its capability in achieving x-ray phase contrast imaging with low brilliance source. However, the widely used phase stepping information extraction method reduces system stability and prolongs data acquisition time by several times compared with conventional x-ray absorption-based imaging. The mechanical stepping can be avoided by using a staggered grating, but at the cost of low vertical spatial resolution. In this paper, employing a modified staggered grating and the angular signal radiography, we proposed a single-shot grating-based x-ray differential phase contrast imaging with decent vertical spatial resolution. The theoretical framework was deduced and proved by numerical experiments. Absorption, phase, and scattering computed tomography can be performed without phase stepping. Therefore, we believe this fast and highly stable imaging method with decent resolution would be widely applied in x-ray grating-based phase contrast imaging.