Abstract
Computational models of complex systems are usually elaborate and sensitive to implementation details, characteristics which often affect their verification and validation. Model replication is a possible solution to this issue. It avoids biases associated with the language or toolkit used to develop the original model, not only promoting its verification and validation, but also fostering the credibility of the underlying conceptual model. However, different model implementations must be compared to assess their equivalence. The problem is, given two or more implementations of a stochastic model, how to prove that they display similar behavior? In this paper, we present a model comparison technique, which uses principal component analysis to convert simulation output into a set of linearly uncorrelated statistical measures, analyzable in a consistent, model-independent fashion. It is appropriate for ascertaining distributional equivalence of a model replication with its original implementation. Besides model-independence, this technique has three other desirable properties: a) it automatically selects output features that best explain implementation differences; b) it does not depend on the distributional properties of simulation output; and, c) it simplifies the modelers' work, as it can be used directly on simulation outputs. The proposed technique is shown to produce similar results to the manual or empirical selection of output features when applied to a well-studied reference model.

Abstract
Agent-based modeling (ABM) is a bottom-up modeling approach, where each entity of the system being modeled is uniquely represented as an independent decision-making agent. Large scale emergent behavior in ABMs is population sensitive. As such, the number of agents in a simulation should be able to reflect the reality of the system being modeled, which can be in the order of millions or billions of individuals in certain domains. A natural solution to reach acceptable scalability in commodity multi-core processors consists of decomposing models such that each component can be independently processed by a different thread in a concurrent manner. In this paper we present a multithreaded Java implementation of the PPHPC ABM, with two goals in mind: (1) compare the performance of this implementation with an existing NetLogo implementation; and, (2) study how different parallelization strategies impact simulation performance on a shared memory architecture. Results show that: (1) model parallelization can yield considerable performance gains; (2) distinct parallelization strategies offer specific trade-offs in terms of performance and simulation reproducibility; and, (3) PPHPC is a valid reference model for comparing distinct implementations or parallelization strategies, from both performance and statistical accuracy perspectives.

Abstract
Optimising the design of water distribution networks (WDNs) is a well-known problem that has been studied by numerous researchers. This work proposes a heuristic based on simulated annealing and improved by using concepts from the cross-entropy method. The proposed optimization approach is presented and used in two case studies of different complexity. The results show not only a fall in the computational effort of the new approach relative to simulated annealing but also include a comparison with other heuristic results from the literature, used to solve the same problems.

Abstract
Knowledge regarding Ecosystem Services (ES) delivery and the socio-ecological factors that influence their proficiency is essential to allow cities to adopt policies that lead to resource-efficient planning and greater resilience. As one of the matrix elements of urban ecological structure, vegetation may play a major role in promoting ES proficiency through planting design. This research addresses the heterogeneity of ES delivered by the urban vegetation of Porto, a Portuguese city. A methodology is proposed to investigate associations between socioeconomic indicators and structural variables of the urban forest, and also which structural variables of the urban forest, if any, differ along a socioeconomic gradient. Our results reveal that before setting planning and management goals, it is crucial to understand local patterns of ES and their relationships with socioeconomic patterns, which can be affected by variables such as building age. This should be followed by the identification of structural variables of the urban forest that better explain the differences, in order to target these through planning and management goals. The conceptual framework adopted in this research can guide adaptation of our methodology to other cities, providing insights for planning and management suitable to site-specific conditions and directly usable by stakeholders.

Abstract
In this study, hyperspectral reflectance (HySR) data derived from a handheld spectroradiometer were used to assess the water status of three grapevine cultivars in two sub-regions of Douro wine region during two consecutive years. A large set of potential predictors derived from the HySR data were considered for modelling/predicting the predawn leaf water potential (Psi(pd)) through different statistical and machine learning techniques. Three HySR vegetation indices were selected as final predictors for the computation of the models and the in-season time trend was removed from data by using a time predictor. The vegetation indices selected were the Normalized Reflectance Index for the wavelengths 554 nm and 561 nm (NRI554;561), the water index (WI) for the wavelengths 900 nm and 970 nm, and the D1 index which is associated with the rate of reflectance increase in the wavelengths of 706 nm and 730 nm. These vegetation indices covered the green, red edge and the near infrared domains of the electromagnetic spectrum. A large set of state-of-the-art analysis and statistical and machine-learning modelling techniques were tested. Predictive Modelling techniques based on generalized boosted model (GBM), bagged multivariate adaptive regression splines (B-MARS), generalized additive model (GAM), and Bayesian regularized neural networks (BRNN) showed the best performance for predicting Psi(pd), with an average determination coefficient (R-2) ranging between 0.78 and 0.80 and RMSE varying between 0.11 and 0.12 MPa. When cultivar Touriga Nacional was used for training the models and the cultivars Touriga Franca and Tinta Barroca for testing (independent validation), the models performance was good, particularly for GBM (R-2 = 0.85; RMSE = 0.09 MPa). Additionally, the comparison of Psi(pd) observed and predicted showed an equitable dispersion of data from the various cultivars. The results achieved show a good potential of these predictive models based on vegetation indices to support irrigation scheduling in vineyard.

Abstract
While an understanding of evolutionary processes in shifting environments is vital in the context of rapid ecological change, one of the most potent selective forces, sexual selection, remains curiously unexplored. Variation in sexual selection across a species range, especially across a gradient of temperature regimes, has the potential to provide a window into the possible impacts of climate change on the evolution of mating patterns. Here, we investigated some of the links between temperature and indicators of sexual selection, using a cold-water pipefish as model. We found that populations differed with respect to body size, length of the breeding season, fecundity, and sexual dimorphism across a wide latitudinal gradient. We encountered two types of latitudinal patterns, either linear, when related to body size, or parabolic in shape when considering variables related to sexual selection intensity, such as sexual dimorphism and reproductive investment. Our results suggest that sexual selection intensity increases toward both edges of the distribution and that the large differences in temperature likely play a significant role. Shorter breeding seasons in the north and reduced periods for gamete production in the south certainly have the potential to alter mating systems, breeding synchrony, and mate monopolization rates. As latitude and water temperature are tightly coupled across the European coasts, the observed patterns in traits related to sexual selection can lead to predictions regarding how sexual selection should change in response to climate change. Based on data from extant populations, we can predict that as the worm pipefish moves northward, a wave of decreasing selection intensity will likely replace the strong sexual selection at the northern range margin. In contrast, the southern populations will be followed by heightened sexual selection, which may exacerbate the problem of local extinction at this retreating boundary.