ROMSPath v1.0: Offline Particle Tracking for the Regional Ocean Modeling System (ROMS)

Offline particle tracking (OPT) is a widely used tool for the analysis of data in oceanographic research. Given the output of a hydrodynamic model, OPT can provide answers to a wide variety of research questions involving fluid kinematics, zooplankton transport, the dispersion of pollutants, and the fate of chemical tracers, among others. In this paper, we introduce ROMSPath, an OPT model designed to complement the Regional Ocean Modelling System (ROMS). Based on the Lagrangian TRANSport (LTRANS) model (North et al., 2008), ROMSPath is written in Fortran 90 and provides advancements in functionality and efficiency compared to LTRANS. First, ROMSPath now calculates particle trajectories using the ROMS native grid, which provides advantages in interpolation, masking, and boundary interaction, while improving accuracy. Second, ROMSPath enables simulated particles to pass between nested ROMS grids, which are an increasingly popular tool to simulate the ocean over multiple scales. Third, the ROMSPath vertical turbulence module enables the turbulent (diffusion) time step and advection time step to be specified separately, adding flexibility and improving computational efficiency. Lastly, ROMSPath includes new infrastructure enabling input of auxiliary parameters for added functionality. In particular, Stokes drift can be input and added to particle advection. Here we describe the details of these updates and improvements.

Study on Full-scale Measurements of Wind Field Characteristics of Coastal Complex Mountainous Terrains

Based on full-scale wind field measurements of coastal complex mountainous terrains, data of fluctuating wind velocities at the height of 30m for four sites, including mountaintop and hillsides, are obtained. The wind load characteristics of mean wind velocities, wind directions, turbulence intensities and speed-up ratios of wind velocities are comprehensively examined. Results show that the maximum mean wind velocity at the mountaintop site is 12.4 m/s. The probability density distributions of mean wind velocities for the four measurement sites can be well represented by the Weibull probability model. The predominant wind directions are around the northeast and southwest. The longitudinal, lateral and vertical turbulence intensities decrease with the increase of mean wind velocities. The turbulence intensities for the mountaintop site are as many as 0.13,0.12 and 0.089 under maximum wind velocities, respectively for the previously mentioned three directions. The speed-up ratios of wind velocities between mountaintop and hillside sites are reduced, as the wind velocities increase. However, in cases of intensive wind with mean wind velocities larger than 8 m/s, the speed-up effects of wind velocities can also appear. The maximum speed-up ratio can reach 1.17.

Observed features of the water masses in the Halmahera Sea in November 2016

Halmahera Sea is one of the locations in the eastern route of Indonesian Throughflow (ITF), where high salinity water is mainly transported by the ITF. The description of water mass in the Halmahera Sea from the perspective of water mass, and related mixing is important. It is not only useful for understanding water mass features, but it can also be used to determine the strength of the turbulent mixing, and so allow how it relates to the water transformation. Here, we report the water mass properties and estimation of mixing quantities in the Halmahera Sea from the CTD profiles based on recent onboard observations during the IOCAS cruise in November 2016. The water mass analysis was done by examining the characteristics of water types in the Temperature-Salinity (T-S) diagram. The mixing estimation uses the density profile derived from temperature and salinity profiles and the quantification of vertical turbulence observed by density overturn. Halmahera Sea is to be found as the location where the thermocline salinity changes abruptly, it is shown from the erosion of salinity maximum in the density of 22-26σθ decreased from the north to the south of the basin. It is associated with strong mixing with spots of higher vertical diffusivity in the thermocline and intermediate layer. In the upper layer, the mixed layer depth in the Halmahera Sea is relatively shallow with an average of about 16.95 m and it is associated with weak wind stress during this month.

Acceleration-Based In Situ Eddy Dissipation Rate Estimation with Flight Data

Inducing civil aviation aircraft to bumpiness, atmospheric turbulence is a typical risk that seriously threatens flight safety. The Eddy Dissipation Rate (EDR) value, as an aircraft-independent turbulence severity indicator, is estimated by a vertical wind-based or aircraft vertical acceleration-based algorithm. Based on the flight data of civil aviation aircraft, the vertical turbulence component is obtained as the input of both algorithms. A new method of computing vertical acceleration response in turbulence is put forward through the Unsteady Vortex Lattice Method (UVLM). The lifting surface of the target aircraft is assumed to be a combination of wing and horizontal tail in a turbulent flight scenario. Vortex rings are assigned on the mean camber surface, forming a non-planar UVLM, to further improve the accuracy. Moreover, the neighboring vortex lattices are placed as close as possible to the structural edge of control surfaces. Thereby, a complete algorithm for estimating vertical acceleration and in situ EDR value from Quick Access Recorder (QAR) flight data is proposed. Experiments show that the aerodynamic performance is computed accurately by non-planar UVLM. The acceleration response by non-planar UVLM is able to track the recorded acceleration data with higher accuracy than that of the linear model. Different acceleration responses at different locations are also obtained. Furthermore, because the adverse effects of aircraft maneuvers are separated from turbulence-induced aircraft bumpiness, the new acceleration-based EDR algorithm shows better accuracy and stability.

Sensitivity of the ocean circulation model to the k–omega vertical turbulence parametrization

Ocean general circulation model (OGCM) of the INM RAS with embedded k turbulent model is developed. The solution of the k model equations depends on the frequencies of buoyancy and velocity shift which are generated by the OGCM. The coefficients of vertical turbulence in OGCM depend on k and omega. The numerical algorithms of both models are based on the splitting method for physical processes. The model equations are split into two stages, describing the three-dimensional transport-diffusion of the kinetic energy of turbulence and frequency and their local generation-dissipation. The system of ordinary differential equations arising at the second stage is solved analytically, which ensures the efficiency of the algorithm. Analytical solution also written for the vertical turbulence coefficient equation. The model is used to study the sensitivity of the model circulation of the North AtlanticArctic Ocean to variations in the parameters of vertical turbulence. Experiments show that varying the coefficients of the analytical model solution can improve the adequacy of the simulation.