Sensitivity of resolved and parameterized surface drag to changes in resolution and parameterization

The relative contribution of resolved and parameterized surface drag towards balancing the atmospheric angular momentum flux convergence (AMFC), and their sensitivity to horizontal resolution and parameterization, are investigated in an atmospheric model. This sensitivity can be difficult to elucidate in free-running climate models, in which the AMFC varies with changing climatologies and, as a result, the relative contributions of surface terms balancing the AMFC also vary. While the sensitivity question has previously been addressed using short-range forecasts, we demonstrate that a nudging framework is an effective method for constraining the AMFC. The Met Office Unified Model is integrated at three horizontal resolutions ranging from 130 km (N96) to 25 km (N512) while relaxing the model’s wind and temperature fields towards the ERAinterim reanalysis within the altitude regions of maximum AMFC. This method is validated against short range forecasts and good agreement is found. These experiments are then used to assess the fidelity of the exchange between parameterized and resolved orographic torques with changes in horizontal resolution. Although the parameterized orographic torque reduces substantially with increasing horizontal resolution, there is little change in resolved orographic torque over 20N to 50N. The tendencies produced by the nudging routine indicate that the additional drag at lower horizontal resolution is excessive. When parameterized orographic blocking is removed at the coarsest of these resolutions, there is a lack of compensation, and even compensation of the opposite sense, by the boundary layer and resolved torques which is particularly pronounced over 20N to 50N. This study demonstrates that there is strong sensitivity in the behaviour of the resolved and parameterized surface drag over this region.

Using laboratory and field measurements to constrain a single habit shortwave optical parameterization for cirrus

A single habit parameterization for the shortwave optical properties of cirrus is presented. The parameterization utilizes a hollow particle geometry, with stepped internal cavities as identified in laboratory and field studies. This particular habit was chosen as both experimental and theoretical results show that the particle exhibits lower asymmetry parameters when compared to solid crystals of the same aspect ratio. The aspect ratio of the particle was varied as a function of maximum dimension, D, in order to adhere to the same physical relationships assumed in the microphysical scheme in a configuration of the Met Office atmosphere-only global model, concerning particle mass, size and effective density. Single scattering properties were then computed using T-Matrix, Ray Tracing with Diffraction on Facets (RTDF) and Ray Tracing (RT) for small, medium, and large size parameters respectively. The scattering properties were integrated over 28 particle size distributions as used in the microphysical scheme. The fits were then parameterized as simple functions of Ice Water Content (IWC) for 6 shortwave bands. The parameterization was implemented into the GA6 configuration of the Met Office Unified Model along with the current operational long-wave parameterization. The GA6 configuration is used to simulate the annual twenty-year short-wave (SW) fluxes at top-of-atmosphere (TOA) and also the temperature and humidity structure of the atmosphere. The parameterization presented here is compared against the current operational model and a more recent habit mixture model.

The role of convective parameterization in the simulation of a cyclone over the South Atlantic

Numerical simulations are carried out to examine the role of the Kuo and Kain-Fritsch (KF) cumulus parameterization schemes and dry dynamics on a cyclone development, in a weak baroclinic atmosphere, over subtropical South Atlantic Ocean. The initial phase of the cyclone development is investigated with a coarse horizontal mesh (75 km) and when the cyclone reaches the mature stage two different horizontal resolutions are used (75 and 25 km). The best performance simulation for the cyclone initial phase occurs when the Kuo convective scheme is applied, and this may be attributed to a greater diabatic warming in the troposphere. On the other hand, the dry simulation is not capable of simulating the correct location and intensity of the cyclone in its initial phase. During the mature phase, a cyclone over deepening occurs in the Kuo scheme experiment associated with larger latent heat release in a deep vertical column. The presence of downdraft currents in the KF scheme, which acts to cool and dry the lower levels, is essential to stabilize the atmosphere and to reproduce the nearest observation cyclone deepening rate. The largest cyclone deepening is found in the Kuo scheme high resolution experiment. This suggests that the KF convective scheme is less sensitive to the horizontal grid resolution. It was also revealed that the diabatic processes are crucial to simulate the observed features of this marine cyclone over subtropical region.

Bayesian inference for a skew-normal IRT model under the centred parameterization

Item response theory (IRT) comprises a set of statistical models which are useful in many fields, especially when there is interest in studying latent variables. These latent variables are directly considered in the Item Response Models (IRM) and they are usually called latent traits. A usual assumption for parameter estimation of the IRM, considering one group of examinees, is to assume that the latent traits are random variables which follow a standard normal distribution. However, many works suggest that this assumption does not apply in many cases. Furthermore, when this assumption does not hold, the parameter estimates tend to be biased and misleading inference can be obtained. Therefore, it is important to model the distribution of the latent traits properly. In this paper we present an alternative latent traits modeling based on the so-called skew-normal distribution; see Genton (2004). We used the centred parameterization, which was proposed by Azzalini (1985). This approach ensures the model identifiability as pointed out by Azevedo et al. (2009b). Also, a Metropolis Hastings within Gibbs sampling (MHWGS) algorithm was built for parameter estimation by using an augmented data approach. A simulation study was performed in order to assess the parameter recovery in the proposed model and the estimation method, and the effect of the asymmetry level of the latent traits distribution on the parameter estimation. Also, a comparison of our approach with other estimation methods (which consider the assumption of symmetric normality for the latent traits distribution) was considered. The results indicated that our proposed algorithm recovers properly all parameters. Specifically, the greater the asymmetry level, the better the performance of our approach compared with other approaches, mainly in the presence of small sample sizes (number of examinees). Furthermore, we analyzed a real data set which presents indication of asymmetry concerning the latent traits distribution. The results obtained by using our approach confirmed the presence of strong negative asymmetry of the latent traits distribution. (C) 2010 Elsevier B.V. All rights reserved.

Bias correction in a multivariate normal regression model with general parameterization

This paper derives the second-order biases Of maximum likelihood estimates from a multivariate normal model where the mean vector and the covariance matrix have parameters in common. We show that the second order bias can always be obtained by means of ordinary weighted least-squares regressions. We conduct simulation studies which indicate that the bias correction scheme yields nearly unbiased estimators. (C) 2009 Elsevier B.V. All rights reserved.

Comparison Of Demand Pattern Calibration In Water Distribution Networks With Geographic And Non-Geographic Parameterization

Demands are one of the most uncertain parameters in a water distribution network model. A good calibration of the model demands leads to better solutions when using the model for any purpose. A demand pattern calibration methodology that uses a priori information has been developed for calibrating the behaviour of demand groups. Generally, the behaviours of demands in cities are mixed all over the network, contrary to smaller villages where demands are clearly sectorised in residential neighbourhoods, commercial zones and industrial sectors. Demand pattern calibration has a final use for leakage detection and isolation. Detecting a leakage in a pattern that covers nodes spread all over the network makes the isolation unfeasible. Besides, demands in the same zone may be more similar due to the common pressure of the area rather than for the type of contract. For this reason, the demand pattern calibration methodology is applied to a real network with synthetic non-geographic demands for calibrating geographic demand patterns. The results are compared with a previous work where the calibrated patterns were also non-geographic.

Parameterization Of Turbulence Models Using 3DVAR Data Assimilation

In this research the 3DVAR data assimilation scheme is implemented in the numerical model DIVAST in order to optimize the performance of the numerical model by selecting an appropriate turbulence scheme and tuning its parameters. Two turbulence closure schemes: the Prandtl mixing length model and the two-equation k-ε model were incorporated into DIVAST and examined with respect to their universality of application, complexity of solutions, computational efficiency and numerical stability. A square harbour with one symmetrical entrance subject to tide-induced flows was selected to investigate the structure of turbulent flows. The experimental part of the research was conducted in a tidal basin. A significant advantage of such laboratory experiment is a fully controlled environment where domain setup and forcing are user-defined. The research shows that the Prandtl mixing length model and the two-equation k-ε model, with default parameterization predefined according to literature recommendations, overestimate eddy viscosity which in turn results in a significant underestimation of velocity magnitudes in the harbour. The data assimilation of the model-predicted velocity and laboratory observations significantly improves model predictions for both turbulence models by adjusting modelled flows in the harbour to match de-errored observations. 3DVAR allows also to identify and quantify shortcomings of the numerical model. Such comprehensive analysis gives an optimal solution based on which numerical model parameters can be estimated. The process of turbulence model optimization by reparameterization and tuning towards optimal state led to new constants that may be potentially applied to complex turbulent flows, such as rapidly developing flows or recirculating flows.

An Improved Parameterization Technique for the Continuation Power Flow

Continuation methods have been long used in P-V curve tracing due to their efficiency in the resolution of ill-conditioned cases, with close to singular Jacobian matrices, such as the maximum loading point of power systems. Several parameterization techniques have been proposed to avoid matrix singularity and successfully solve those cases. This paper presents a simple geometric parameterization technique to overcome the singularity of the Jacobian matrix by the addition of a line equations located at the plane determined by a bus voltage magnitude and the loading factor. This technique enlarges the set of voltage variables that can be used to whole P-V curve tracing, without ill-conditioning problems and no need of parameter changes. Simulation results, obtained for large realistic Brazilian and American power systems, show that the robustness and efficiency of the conventional power flow are not only preserved but also improved.

A Parameterization Technique for the Continuation Power Flow Developed from the Analysis of Power Flow Curves

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Study of alternative schemes for the parameterization step of the continuation power flow method based on physical parameters, part I: Mathematical modeling

The conventional power flow method is considered to be inadequate to obtain the maximum loading point because of the singularity of Jacobian matrix. Continuation methods are efficient tools for solving this kind of problem since different parameterization schemes can be used to avoid such ill-conditioning problems. This paper presents the details of new schemes for the parameterization step of the continuation power flow method. The new parameterization options are based on physical parameters, namely, the total power losses (real and reactive), the power at the slack bus (real or reactive), the reactive power at generation buses, and transmission line power losses (real and reactive). The simulation results obtained with the new approach for the IEEE test systems (14, 30, 57, and 118 buses) are presented and discussed in the companion paper. The results show that the characteristics of the conventional method are not only preserved but also improved.

Study of alternative schemes for the parameterization step of the continuation power flow method based on physical parameters, part II: Performance evaluation

New parameterization schemes have been proposed by the authors in Part I of this paper. In this part these new options for the parameterization of power flow equations are tested, namely, the total power losses (real and reactive), the power at the slack bus (real or reactive), the reactive power at generation buses, and the transmission line power losses (real and reactive). These different parameterization schemes can be used to obtain the maximum loading point without ill-conditioning problems, once the singularity of Jacobian matrix is avoided. The results obtained with the new approach for the IEEE test systems (14, 30, 57, and 118 buses) show that the characteristics of the conventional method are not only preserved but also improved. In addition, it is shown that the proposed method and the conventional one can be switched during the tracing of PV curves to determine, with few iterations, all points of the PV curve. Several tests were also carried out to compare the performance of the proposed parameterization schemes for the continuation power flow method with the use of both the secant and tangent predictors.

An efficient geometric parameterization technique for the continuation power flow

Continuation methods have been shown as efficient tools for solving ill-conditioned cases, with close to singular Jacobian matrices, such as the maximum loading point of power systems. Some parameterization techniques have been proposed to avoid matrix singularity and successfully solve those cases. This paper presents a new geometric parameterization scheme that allows the complete tracing of the P-V curves without ill-conditioning problems. The proposed technique associates robustness to simplicity and, it is of easy understanding. The Jacobian matrix singularity is avoided by the addition of a line equation, which passes through a point in the plane determined by the total real power losses and loading factor. These two parameters have clear physical meaning. The application of this new technique to the IEEE systems (14, 30, 57, 118 and 300 buses) shows that the best characteristics of the conventional Newton's method are not only preserved but also improved. (C) 2006 Elsevier B.V. All rights reserved.

Parameterization technique for the continuation power flow using the trivial and tangent predictor

This paper presents efficient geometric parameterization techniques using the tangent and the trivial predictors for the continuation power flow, developed from observation of the trajectories of the load flow solution. The parameterization technique eliminates the Jacobian matrix singularity of load flow, and therefore all the consequent problems of ill-conditioning, by the addition of the line equations which pass through the points in the plane determined by the variables loading factor and the real power generated by the slack bus, two parameters with clear physical meaning. This paper also provides an automatic step size control around the maximum loading point. Thus, the resulting method enables not only the calculation of the maximum loading point, but also the complete tracing of P-V curves of electric power systems. The technique combines robustness with ease of understanding. The results to the IEEE 300-bus system and of large real systems show the effectiveness of the proposed method. © 2012 IEEE.

Aprendizagem Baseada em Problemas (PBL – Problem-Based Learning): estratégia para o ensino e aprendizagem de algoritmos e conteúdos computacionais

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Parameterization of bio-optical models for estimating chlorophyll-a concentration in a tropical eutrophic reservoir

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)