Six bias correction (BC) methods; delta-method (DT), linear scaling (LS), power transformation of precipitation (PTR), empirical quantile mapping (EQM), gamma quantile mapping (GQM) and gamma-pareto quantile mapping (GPQM) were applied to adjust the biases of historical monthly precipitation outputs from five General Circulation Models (GCMs) dynamically downscaled by two Regional Climate Models (RCMs) for a total of seven different GCM-RCM pairs over Costa Rica. High-resolution gridded precipitation observations were used for the control period 1951–1980 and validated over the period 1981–1995. Results show that considerable biases exist between uncorrected GCM-RCM outputs and observations, which largely depend on GCM-RCM pair, seasonality, climatic region and spatial resolution. After the application of bias correction, substantial biases reductions and comparable performances among most BC methods were observed for most GCM-RCM pairs; with EQM and DT marginally outperforming the remaining methods. Consequently, EQM and DT were selectively applied to correct the biases of precipitation projections from each individual GCM-RCM pair for a near-future (2011–2040), mid-future (2041–2070) and far-future (2071–2100) period under Representative Concentration Pathways (RCPs) 2.6, 4.5 and 8.5 using the control period 1961–1990. Results from the bias-corrected future ensemble-mean anticipate a marked decreasing trend in precipitation from near to far-future periods during the dry season (December, January, February (DJF) and March, April, May (MAM)) for RCP4.5 and 8.5; with pronounced drier conditions for those climatic regions draining towards the Pacific Ocean. In contrast, mostly wetter conditions are expected during the dry season under RCP2.6, particularly for the Caribbean region. In most of the country, the greatest decrease in precipitation is projected at the beginning of the rainy season (June, July, August (JJA)) for the far-future period under RCP8.5, except for the Caribbean region where mostly wetter conditions are anticipated. Regardless of future period, slight increases in precipitation with higher radiative forcing are expected for SON excluding the Caribbean region, where precipitation is likely to increase with increasing radiative forcing and future period. This study demonstrates that bias correction should be considered before direct application of GCM-RCM precipitation projections over complex territories such as Costa Rica.

In this paper, the HBV-TEC hydrological model is presented. The development of this model version, aims to provide researchers and scholars with a stable and robust implementation of the HBV hydrological model based on the R programming language. To evaluate its performance, the HBV-TEC model was applied to three subcatchments of the Aguacaliente river catchment, an experimental catchment located in the province of Cartago, Costa Rica. Results suggest a satisfactory performance of the model for two subcatchments and an unsatisfactory performance of the remaining subcatchment; most of which could be attributed to unsufficient meteorological data along with a highly heterogeneous spatial rainfall-distribution

Precipitation climatologies for the period 1961–1990 were generated for all climatic regions of Costa Rica using an irregular rain-gauge observational network comprised by 416 rain-gauge stations. Two sub-networks were defined: a high temporal resolution sub-network (HTR), including stations having at least 20 years of continuous records during the study period (157 in total); and a high spatial resolution sub-network (HSR), which includes all HTR-stations plus those stations with less than 20 years of continuous records (416 in total). Results from the kriging variance reduction efficiency (KRE) objective function between the two sub-networks, show that ordinary kriging (OK) is unable to fully explain the spatio-temporal variability of precipitation within most climatic regions if only stations from the HTR sub-network are used. Results also suggests that in most cases, it is beneficial to increase the density of the rain-gauge observational network at the expense of temporal fidelity, by including more stations even though their records may not represent the same time step. Thereafter, precipitation climatologies were generated using seven deterministic (IDW, TS2, TS2PARA, TS2LINEAR, TPS, MQS and NN) and two geostatistical (OK and KED) interpolation methods. Performance of the various interpolation methods was evaluated using cross validation technique, selecting the mean absolute error (MAE) and the root-mean square error (RMSE) as agreement metrics. Results suggest that IDW is marginally superior to OK and KED for most climatic regions. The remaining deterministic methods however, considerably deviate from IDW, which suggests that these methods are incapable of properly capturing the true-nature of spatial precipitation patterns over the considered climatic regions. The final generated IDW climatology was then validated against the Global Precipitation Climatology Centre (GPCC), Climate Research Unit (CRU) and WorldClim datasets, in which overall spatial and temporal coherence is considered satisfactory, giving assurance about the use this new climatology in the development of local climate impact studies.

This study aimed to assess the impacts of climate change on streamflow characteristics of five tropical catchments located in Costa Rica. An ensemble of five General Circulation Models (GCMs), namely HadGEM2-ES, CanESM2, EC-EARTH, MIROC5, MPI-ESM-LR dynamically downscaled by two Regional Climate Models (RCMs), specifically HadRM3P and RCA4, was selected to provide an overview of the impacts of different climate change scenarios under Representative Concentration Pathways (RCPs) 2.6, 4.5 and 8.5 using the 1961–1990 baseline period. The GR2M hydrological model was used to reproduce the historical monthly surface runoff patterns of each catchment. Following calibration and validation of the GRM2 model, the projected impact of climate change on streamflow was simulated for a near-future (2011–2040), mid-future (2041–2070) and far-future (2071–2100) for each catchment using the bias-corrected GCM-RCM multimodel ensemble-mean (MEM). Results anticipate wetter conditions for all catchments in the near-future and mid-future periods under RCPs 2.6 and 4.5, whereas dryer conditions are expected for the far-future period under RCP 8.5. Projected temperature trends indicate consistently warmer conditions with increasing radiative forcing and future periods. Streamflow changes across all catchments however are dominated by variations in projected precipitation. Wetter conditions for the near-future and mid-future horizons under RCPs 2.6 and 4.5 would result in higher runoff volumes, particularly during the late wet season (LWS). Conversely, dryer conditions for the far-future period under RCP8.5 would result in considerably lower runoff volumes during the early wet season (EWS) and the Mid-Summer Drought (MSD). In consequence, projected seasonal changes on streamflow across all catchments may result in more frequent flooding, droughts, and water supply shortage compared to historical hydrological regimes.

The hydraulic model EPANET was applied and calibrated for the water distribution system (WDS) of La Sirena, Colombia. The Parameter ESTimator (PEST) was used for parameter optimization and sensitivity analysis. Observation data included levels at water storage tanks and pressures at monitoring nodes. Adjustable parameters were grouped into different classes according to two different scenarios identified as constrained and unconstrained. These scenarios were established to evaluate the effect of parameter space size and compensating errors over the calibration process. Results from the unconstrained scenario, where 723 adjustable parameters were declared, showed that considerable compensating errors are introduced into the optimization process if all parameters were open to adjustment. The constrained scenario on the other hand, represented a more properly discretized scheme as parameters were grouped into classes of similar characteristics and insensitive parameters were fixed. This had a profound impact on the parameter space as adjustable parameters were reduced to 24. The constrained solution, even when it is valid only for the system's normal operating conditions, clearly demonstrates that Parallel PEST (PPEST) has the potential to be used in the calibration of WDS models. Nevertheless, further investigation is needed to determine PPEST's performance in complex WDS models.