Validation of Three Daily Satellite Rainfall Products in a Humid Tropic Watershed, Brantas, Indonesia: Implications to Land Characteristics and Hydrological Modelling

A total of three different satellite products, CHIRPS, GPM, and PERSIANN, with different spatial resolutions, were examined for their ability to estimate rainfall data at a pixel level, using 30-year-long observations from six locations. Quantitative and qualitative accuracy indicators, as well as R2 and NSE from hydrological estimates, were used as the performance measures. The results show that all of the satellite estimates are unsatisfactory, giving the NRMSE ranging from 6 to 30% at a daily level, with CC only 0.21–0.36. Limited number of gauges, coarse spatial data resolution, and physical terrain complexity were found to be linked with low accuracy. Accuracy was slightly better in dry seasons or low rain rate classes. The errors increased exponentially with the increase in rain rates. CHIPRS and PERSIANN tend to slightly underestimate at lower rain rates, but do show a consistently better performance, with an NRMSE of 6–12%. CHRIPS and PERSIANN also exhibit better estimates of monthly flow data and water balance components, namely runoff, groundwater, and water yield. GPM has a better ability for rainfall event detections, especially during high rainfall events or extremes (>40 mm/day). The errors of the satellite products are generally linked to slope, wind, elevation, and evapotranspiration. Hydrologic simulations using SWAT modelling and the three satellite rainfall products show that CHIRPS slightly has the daily best performance, with R2 of 0.59 and 0.62, and NSE = 0.54, and the monthly aggregated improved at a monthly level. The water balance components generated at an annual level, using three satellite products, show that CHIRPS outperformed with a ration closer to one, though with a tendency to overestimate up to 3–4× times the data generated from the rainfall gauges. The findings of this study are beneficial in supporting efforts for improving satellite rainfall products and water resource implications.

Time-dependent solutions for daily-periodic slope flows driven by surface energy budget

<p>Diurnal wind systems generated from daytime heating and nighttime cooling of valleys and slopes are a very common feature over mountainous terrains. Despite their frequent occurrence and relevance for a variety of applications, ranging from pollutant transport to convection initiation, slope winds are far from being fully understood and still provide an open research topic.</p><p>A well-known steady-state analytical model <span>is the one</span> developed by Prandtl (1942). Then, a first time-dependent analytical model was proposed by F. Defant (1949) and later extended by Zardi and Serafin (2015). These models provide slope-normal profiles of temperature and along-slope wind velocity as a response to a sinusoidal forcing representing the surface temperature. The resulting profiles exhibit sinusoidal oscillations at every distance from the surface, although with different phase lags under different regimes, determined by different combinations of slope angle and stability of the unperturbed ambient atmosphere. As a consequence, they can not explain the observed differences between daytime upslope and nighttime downslope winds in terms of magnitude and height of the peak of wind velocity, as well as the different timing characterising nighttime, daytime, and the two reversal phases.</p><p>In the present work, the solutions derived in Zardi and Serafin (2015) are extended to include a more realistic daily-periodic surface forcing taking into account the daily evolution of the surface temperature computed on the basis of a surface energy <span>budget</span>. Incoming solar radiation is represented by means of a Fourier series expansion derived from well-established relationships taking into account latitude, day of the year, slope angle, exposition and other astronomical and atmospheric factors. Based on <span>these</span> expansions, suitable harmonic solutions are derived for the heat flux into the ground and sensible heat flux in the atmosphere, and hence for the daily evolution of slope-normal profiles of along-slope wind velocity and potential temperature.</p><p>References:</p><ul><li><span>Prandtl L. 1942. Führer durch die Strömungslehre, Chapter 5. Vieweg und Sohn: Braunschweig, Germany. </span>[English translation: Prandtl L. 1952. Mountain and valley winds in stratified air, in Essentials of Fluid Dynamics: 422–425. Hafner Publishing Company: New York, NY]</li> <li><span>Defant F. 1949. Zur Theorie der Hangwinde, nebst Bemerkungen zur Theorie der Berg- und Talwinde. </span>Arch. Meteorol. Geophys. Bioklimatol. A 1: 421–450</li> <li>Zardi D., Serafin S. 2015. An analytic solution for time‐periodic thermally driven slope flows. Q. J. R. Meteorol. Soc., 141, 1968–1974, https://doi.org/10.1002/qj.2485</li> </ul>

Characterization of the morning transition from downslope to upslope winds and its connection with the nocturnal inversion breakup at the foot of a gentle slope.

<p>Thermally driven winds observed in complex terrain are characterized by a daily cycle dominated by two main phases: a diurnal phase in which winds blow upslope (anabatic), and a nocturnal one in which they revert their direction and blow downslope (katabatic). This alternating pattern also implies two transition phases, following sunrise and sunset respectively. </p><p>Here we study the upslope component of the slope wind with a focus on the morning transition from the katabatic to the anabatic flow based on data from the MATERHORN experiment, performed in Salt Lake Desert (Utah) between Fall 2012 and Spring 2013 (Fernando et al, 2015). </p><p>First of all, a criterion for the selection of purely thermally driven slope wind days is proposed and adopted to select five case studies, taken from both the spring and the autumn periods. Then, the analysis allowed the investigation of the driving mechanisms through the connection with the patterns of erosion of the nocturnal inversion in the valley bed at the foot of the slope under analysis. Three main patterns of erosion of the inversion in the particular topography of a gentle and isolated slope are identified: a) erosion due to upward growth of a convective boundary layer, b) erosion due to descent of the inversion top, and c) erosion due to a mix of the two previous mechanisms. The three patterns are then linked to the initiation of the transition by two different and competing mechanisms: mixing from above (top-down dilution) and surface heating from below. Finally, an analytical model for the description of slope circulation (Zardi and Serafin, 2015) has been used to diagnose the time of the transition.</p><p>Zardi, D. and S. Serafin, 2015: An analytic solution for daily-periodic thermally-driven slope flow. Quart. J. Roy. Meteor. Soc., 141, 1968–1974. </p><p>Fernando, H. J. S., Pardyjak, E. R., Di Sabatino, S., Chow, F. K., De Wekker, S. F. J., Hoch, S. W, Zsedrovits, T., 2015, The MATERHORN: Unraveling the intricacies of mountain weather. <em>Bulletin Of The American Meteorological Society</em>, 96, 1945-1967. </p>

Characterization of the morning transition from downslope to upslope winds and its connection with the nocturnal inversion breakup at the foot of a gentle slope

<p>Thermally driven winds observed in complex terrain are characterized by a daily cycle dominated by two main phases: a diurnal phase in which winds blow upslope (anabatic), and a nocturnal one in which they revert their direction and blow down slope (katabatic). This alternating pattern also implies two transition phases, following sunrise and sunset respectively. </p><p>Here we study the up-slope component of the slope wind with a focus on the morning transition based on from the MATERHORN experiment, performed in Salt Lake Desert (Utah) between Fall 2012 and Spring 2013. </p><p>The analysis develops along three main paths of investigation. The first one is the selection of the suitable conditions for the study of the diurnal component and the characterization of the morning transition. The second one focuses on the deep analysis of the erosion of the nocturnal inversion at the foot of the slope in order to investigate the physical mechanisms driving it. And the third one consists in the comparison between the experimental data and the results of an analytical model (Zardi and Serafin, 2015). The study of the morning transition in the selected case studies allowed its characterization in terms of the relation with the solar radiation cycle, in terms of its seasonality and in terms of its propagation along the slope and along the vertical direction. Most of the results of this investigation are related to the identification of the main mechanisms of erosion of the nocturnal inversion at the foot of the slope and to its role to the beginning of the transition itself. Finally, it is shown how the above model can fairly reproduce the cycle between anabatic and katabatic flow and their intensity.</p><p>Zardi, D. and S. Serafin, 2015: An analytic solution for daily-periodic thermally-driven slope flow. Quart. J. Roy. Meteor. Soc., 141, 1968–1974.</p>

Investigating the dynamics of thermally driven up-slope flows: analysis of data from measurements in the Inn Valley (Austria) and comparison with a simple model

<p><span>Diurnal wind systems typically develop in mountainous areas following the daytime heating and nighttime cooling of sloping surfaces. While down-slope winds have been extensively treated in the literature, up-slope winds have received much less attention. In particular, the physical mechanisms associated with the development of these winds, as well as the search for appropriate parameterization of turbulent fluxes of mass, momentum, and heat over slopes in numerical weather prediction models are still open research topics.</span></p><p><span>Here we present some preliminary results from the analysis of turbulence data (sonic wind speed, temperature, humidity, and turbulent fluxes) collected at two slope stations which are part of the i-Box initiative. The i-Box project (Rotach et al. 2017) aims at studying turbulent exchange processes in complex terrain areas. The experimental setup is composed of six stations disseminated in the surroundings of the alpine city of Innsbruck, in the Inn Valley. The two stations adopted for the present study are located at different points on the valley sidewalls, one with a slope angle of 27° (labelled NF27) and one with a slope angle of 10° (NF10). Both stations are located over slopes covered by alpine meadow and at an altitude of about 1000 m MSL (400 m above the valley floor). The station NF27 has two measurement points, 1.5 and 6.8 m AGL, while the station NF10 has one measurement point, at 6.2 m AGL.</span></p><p><span>The analysis shows that criteria proposed in the literature for the selection of valley-wind days may not apply for the identification of slope-wind days. Furthermore, from the analysis of second order moments, scaling relationships are derived for up-slope flow conditions. In addition, measurements representing the evolution of the up-slope flow structure from the early morning to the mid-afternoon are compared with an existing, simplified, analytical model, which provides the evolution of the vertical profiles of temperature and along-slope wind velocity as generated by a sinusoidal forcing representing the daily cycle of surface temperature. An improvement of the existing model, where the surface energy budget is considered as the boundary condition for the surface temperature, is also tested.</span></p>

Temporal and spatial variability in surface roughness and accumulation rate around 88° S from repeat airborne geophysical surveys

We use repeat high-resolution airborne geophysical data consisting of laser altimetry, snow, and Ku-band radar and optical imagery acquired in 2014, 2016, and 2017 to analyze the spatial and temporal variability in surface roughness, slope, wind deposition, and snow accumulation at 88∘ S, an elevation bias validation site for ICESat-2 and potential validation site for CryoSat-2. We find significant small-scale variability (<10 km) in snow accumulation based on the snow radar subsurface stratigraphy, indicating areas of strong wind redistribution are prevalent at 88∘ S. In general, highs in snow accumulation rate correspond with topographic lows, resulting in a negative correlation coefficient of r2=-0.32 between accumulation rate and MSWD (mean slope in the mean wind direction). This relationship is strongest in areas where the dominant wind direction is parallel to the survey profile, which is expected as the geophysical surveys only capture a two-dimensional cross section of snow redistribution. Variability in snow accumulation appears to correlate with variability in MSWD. The correlation coefficient between the standard deviations of accumulation rate and MSWD is r2=0.48, indicating a stronger link between the standard deviations than the actual parameters. Our analysis shows that there is no simple relationship between surface slope, wind direction, and snow accumulation rates for the overall survey area. We find high variability in surface roughness derived from laser altimetry measurements on length scales smaller than 10 km, sometimes with very distinct and sharp transitions. Some areas also show significant temporal variability over the course of the 3 survey years. Ultimately, there is no statistically significant slope-independent relationship between surface roughness and accumulation rates within our survey area. The observed correspondence between the small-scale temporal and spatial variability in surface roughness and backscatter, as evidenced by Ku-band radar signal strength retrievals, will make it difficult to develop elevation bias corrections for radar altimeter retrieval algorithms.

Temporal and spatial variability in surface roughness and accumulation rate around 88° S from repeat airborne geophysical surveys

We use repeat high-resolution airborne geophysical data consisting of laser altimetry, snow and Ku-band radar and optical imagery acquired in 2014, 2016 and 2017 to analyze the spatial and temporal variability in surface roughness, slope, wind deposition, and snow accumulation at 88° S as this is a bias validation site for ICESat-2 and may be a potential validation site for CryoSat-2. We find significant small–scale variability (

Evolution of Late Wintertime Thermally Driven Local Flows and Temperature Fields over Far-Western Nepal

Atmospheric processes over the Himalayan complex terrain are yet to be studied extensively. Only a few significant researches are reported from this region and the Far-Western Region (FWR) of Nepal still remains untouched. Thus, the present study was conceived to understand the meteorological flow characteristics and thermal environment over the region and associated areas during the late wintertime with the application of the state-of-the-art-of Weather Research and Forecasting (WRF) Modeling System. The study revealed that the northern mountainous region developed strong down slope wind during the night and morning times, which sweeps out the southern plain area of Nepal and may reach just beyond the border. The wind over the plain was very shallow whose depth was just about 100 m. The down slope winds over the southern slope of the Daijee and Nandhaur mountain ranges were significantly enhanced by the subsidence of the southerly wind that prevails above 1 km height above the mean sea level. Close to the noon time a very gentle southerly valley wind from the southern plain replaced the nighttime down slope. Very shallow but strong surface inversion builds up over the plain that breaks up in the late morning. The depth of the mixed layer and the valley wind may reach up to 1km in the afternoon. The thermal environment over the FWR of Nepal was fairly hot that may remain around 35°C in the afternoon around the Mahendranagar area whereas the temperature during the nighttime may go as low as 23°C. The study revealed that, contrary to the general perception, temperature over plain areas of Nepal was significantly higher than further southern areas belonging to India. The meteorological flow fields over the FWR of Nepal executed diurnal periodicity with little day-to-day variation during the late wintertime.Journal of Institute of Science and TechnologyVolume 21, Issue 1, August 2016, page: 35-47

Quantifying horizontal and vertical tracer mass fluxes in an idealized valley during daytime

The transport and mixing of pollution during the daytime evolution of a valley boundary layer is studied in an idealized way. The goal is to quantify horizontal and vertical tracer mass fluxes between four different valley volumes: the convective boundary layer, the slope wind layer, the stable core, and the atmosphere above the valley. For this purpose, large eddy simulations (LES) are conducted with the Weather Research and Forecasting (WRF) model for a quasi-two-dimensional valley. The valley geometry consists of two slopes with constant slope angle and is homogeneous in the along-valley direction. The surface sensible heat flux is horizontally homogeneous and prescribed by a sine function. The initial sounding is characterized by an atmosphere at rest and a constant Brunt–Väisälä frequency. Various experiments are conducted for different combinations of surface heating amplitudes and initial stability conditions. A passive tracer is released with an arbitrary but constant rate at the valley floor and resulting tracer mass fluxes are evaluated between the aforementioned volumes.As a result of the surface heating, a convective boundary layer is established in the lower part of the valley with a stable layer on top – the so-called stable core. The height of the slope wind layer, as well as the wind speed within, decreases with height due to the vertically increasing stability. Hence, the mass flux within the slope wind layer decreases with height as well. Due to mass continuity, this along-slope mass flux convergence leads to a partial redirection of the flow from the slope wind layer towards the valley centre and the formation of a horizontal intrusion above the convective boundary layer. This intrusion is associated with a transport of tracer mass from the slope wind layer towards the valley centre. A strong static stability and/or weak forcing lead to large tracer mass fluxes associated with this phenomenon. The total export of tracer mass out of the valley atmosphere increases with decreasing stability and increasing forcing. The effects of initial stability and forcing can be combined to a single parameter, the breakup parameter B. An analytical function is presented that describes the exponential decrease of the percentage of exported tracer mass with increasing B. This study is limited by the idealization of the terrain shape, stratification, and forcing, but quantifies transport processes for a large range of forcing amplitudes and atmospheric stability.