Consistent seasonal snow cover depth and duration variability over the Western Himalayas (WH)

Precipitation in solid form, i.e., snow, during winter season over the Western Himalayas (WH) leads to the build-up of seasonal snow cover. Seasonal snow cover build-up (snow cover depth and duration) largely depends on atmospheric variables such as temperature, precipitation, radiation, wind, etc. Integrated (combined) influence of atmospheric variables on seasonal snow cover gets reflected in terms of spatial and temporal variability in seasonal snow cover build-up pattern. Hence spatial and temporal variability of seasonal snow cover build-up can serve as a good indicator of climate change in high altitude mountainous regions like the WH. Consistent seasonal snow cover depth and duration, delay days and early melt days of consistent seasonal snow cover at 11 stations spread across different mountain ranges over the WH were analyzed. Mean, maximum and percentiles (25th, 50th, 75th, 90th and 95th) of consistent seasonal snow cover depth and duration show decline over the WH in the recent past 2–3 decades. Consistent seasonal snow cover is found to melt early and snow cover build-up pattern is found to show changes over the WH. Decline in consistent seasonal snow cover depth, duration and changing snow cover build-up pattern over the WH in recent decades indicate that WH has undergone considerable climate change and winter weather patterns are changing in the WH.     

Role of land state in a high resolution mesoscale model for simulating the Uttarakhand heavy rainfall event over India

In 2013, Indian summer monsoon witnessed a very heavy rainfall event (>30 cm/day) over Uttarakhand in north India, claiming more than 5000 lives and property damage worth approximately 40 billion USD. This event was associated with the interaction of two synoptic systems, i.e., intensified subtropical westerly trough over north India and north-westward moving monsoon depression formed over the Bay of Bengal. The event had occurred over highly variable terrain and land surface characteristics. Although global models predicted the large scale event, they failed to predict realistic location, timing, amount, intensity and distribution of rainfall over the region. The goal of this study is to assess the impact of land state conditions in simulating this severe event using a high resolution mesoscale model. The land conditions such as multi-layer soil moisture and soil temperature fields were generated from High Resolution Land Data Assimilation (HRLDAS) modelling system. Two experiments were conducted namely, (1) CNTL (Control, without land data assimilation) and (2) LDAS, with land data assimilation (i.e., with HRLDAS-based soil moisture and temperature fields) using Weather Research and Forecasting (WRF) modelling system. Initial soil moisture correlation and root mean square error for LDAS is 0.73 and 0.05, whereas for CNTL it is 0.63 and 0.053 respectively, with a stronger heat low in LDAS. The differences in wind and moisture transport in LDAS favoured increased moisture transport from Arabian Sea through a convectively unstable region embedded within two low pressure centers over Arabian Sea and Bay of Bengal. The improvement in rainfall is significantly correlated to the persistent generation of potential vorticity (PV) in LDAS. Further, PV tendency analysis confirmed that the increased generation of PV is due to the enhanced horizontal PV advection component rather than the diabatic heating terms due to modified flow fields. These results suggest that, two different synoptic systems merged by the strong interaction of moving PV columns resulted in the strengthening and further amplification of the system over the region in LDAS. This study highlights the importance of better representation of the land surface fields for improved prediction of localized anomalous weather event over India.     

Characteristics of the Asian–Pacific Oscillation in Boreal Summer Simulated by BCC_CSM with Different Horizontal Resolutions

The summer Asian–Pacific Oscillation(APO) is a major teleconnection pattern that reflects the zonal thermal contrast between East Asia and the North Pacific in the upper troposphere. The performance of Beijing Climate Center Climate System Models(BCC CSMs) with different horizontal resolutions, i.e., BCC CSM1.1 and BCC CSM1.1(m), in reproducing APO interannual variability, APO-related precipitation anomalies, and associated atmospheric circulation anomalies, is evaluated.The results show that BCC CSM1.1(m) can successfully capture the interannual variability of the summer APO index. It is also more capable in reproducing the APO's spatial pattern, compared to BCC CSM1.1, due to its higher horizontal resolution. Associated with a positive APO index, the northward-shifted and intensified South Asian high, strengthened extratropical westerly jet, and tropical easterly jet in the upper troposphere, as well as the southwesterly monsoonal flow over North Africa and the Indian Ocean in the lower troposphere, are realistically represented by BCC CSM1.1(m), leading to an improvement in reproducing the increased precipitation over tropical North Africa, South Asia, and East Asia, as well as the decreased precipitation over subtropical North Africa, Japan, and North America. In contrast, these features are less consistent with observations when simulated by BCC CSM1.1. Regression analysis further indicates that surface temperature anomalies over the North Pacific and the southern and western flanks of the Tibetan Plateau are reasonably reproduced by BCC CSM1.1(m), which contributes to the substantial improvement in the simulation of the characteristics of summer APO compared to that of BCC CSM1.1.     

An algorithm to retrieve precipitation with synthetic aperture radar

This paper presents a new type of rainfall retrieval algorithm, called the model-oriented statistical and Volterra integration. It is a combination of the model-oriented statistical (MOS) and Volterra integral equation (VIE) approaches. The steps involved in this new algorithm can be briefly illustrated as follows. Firstly, information such as the start point and width of the rain is obtained through pre-analysis of the data received by synthetic aperture radar (SAR). Secondly, the VIE retrieval algorithm is employed over a short distance to obtain information on the shape of the rain. Finally, the rain rate can be calculated by using the MOS retrieval algorithm. Simulation results show that the proposed algorithm is effective and simple, and can lead to time savings of nearly 50% compared with MOS. An example of application of SAR data is also discussed, involving the retrieval of precipitation information over the South China Sea.     

Estimation of Above Ground Biomass for Central Indian Deciduous Forests Using ALOS PALSAR L-Band Data

Reliable and accurate estimates of tropical forest above ground biomass (AGB) are important to reduce uncertainties in carbon budgeting. In the present study we estimated AGB of central Indian deciduous forests of Madhya Pradesh (M.P.) state, India, using Advanced Land Observing Satellite – Phased Array type L-band Synthetic Aperture Radar (ALOS-PALSAR) L-band data of year 2010 in conjunction with field based AGB estimates using empirical models. Digital numbers of gridded 1?×?1° dual polarization (HH & HV) PALSAR mosaics for the study area were converted to normalized radar cross section (sigma naught - σ⁰). A total of 415 sampling plots (0.1 ha) data collected over the study area during 2009–10 was used in the present study. Plot-level AGB estimates using volume equations representative to the study area were computed using field inventory data. The plot-level AGB estimates were empirically modeled with the PALSAR backscatter information in HH, HV and their ratios from different forest types of the study area. The HV backscatter information showed better relation with field based AGB estimates with a coefficient of determination (R²) of 0.509 which was used to estimate spatial AGB of the study area. Results suggested a total AGB of 367.4 Mt for forests of M.P. state. Further, validation of the model was carried out using observed vs. predicted AGB estimates, which suggested a root mean square error (RMSE) of ±19.32 t/ha. The model reported robust and defensible relation for observed vs. predicted AGB values of the study area.     

Using a regional numerical weather prediction model for GNSS positioning over Brazil

The global navigation satellite system (GNSS) can provide centimeter positioning accuracy at low costs. However, in order to obtain the desired high accuracy, it is necessary to use high-quality atmospheric models. We focus on the troposphere, which is an important topic of research in Brazil where the tropospheric characteristics are unique, both spatially and temporally. There are dry regions, which lie mainly in the central part of the country. However, the most interesting area for the investigation of tropospheric models is the wet region which is located in the Amazon forest. This region substantially affects the variability of humidity over other regions of Brazil. It provides a large quantity of water vapor through the humidity convergence zone, especially for the southeast region. The interconnection and large fluxes of water vapor can generate serious deficiencies in tropospheric modeling. The CPTEC/INPE (Center for Weather Forecasting and Climate Studies/Brazilian Institute for Space Research) has been providing since July 2012 a numerical weather prediction (NWP) model for South America, known as Eta. It has yield excellent results in weather prediction but has not been used in GNSS positioning. This NWP model was evaluated in precise point positioning (PPP) and network-based positioning. Concerning PPP, the best positioning results were obtained for the station SAGA, located in Amazon region. Using the NWP model, the 3D RMS are less than 10 cm for all 24 h of data, whereas the values reach approximately 60 cm for the Hopfield model. For network-based positioning, the best results were obtained mainly when the tropospheric characteristics are critical, in which case an improvement of up to 7.2 % was obtained in 3D RMS using NWP models.     

Generating GPS satellite fractional cycle bias for ambiguity-fixed precise point positioning

With the development of precise point positioning (PPP), the School of Geodesy and Geomatics (SGG) at Wuhan University is now routinely producing GPS satellite fractional cycle bias (FCB) products with open access for worldwide PPP users to conduct ambiguity-fixed PPP solution. We provide a brief theoretical background of PPP and present the strategies and models to compute the FCB products. The practical realization of the two-step (wide-lane and narrow-lane) FCB estimation scheme is described in detail. With GPS measurements taken in various situations, i.e., static, dynamic, and on low earth orbit (LEO) satellites, the quality of FCB estimation and the effectiveness of PPP ambiguity resolution (AR) are evaluated. The comparison with CNES FCBs indicated that our FCBs had a good consistency with the CNES ones. For wide-lane FCB, almost all the differences of the two products were within ±0.05 cycles. For narrow-lane FCB, 87.8 % of the differences were located between ±0.05 cycles, and 97.4 % of them were located between ±0.075 cycles. The experimental results showed that, compared with conventional ambiguity-float PPP, the averaged position RMS of static PPP can be improved from (3.6, 1.4, 3.6) to (2.0, 1.0, 2.7) centimeters for ambiguity-fixed PPP. The average accuracy improvement in the east, north, and up components reached 44.4, 28.6, and 25.0 %, respectively. A kinematic, ambiguity-fixed PPP test with observation of 80 min achieved a position accuracy of better than 5 cm at the one-sigma level in all three coordinate components. Compared with the results of ambiguity-float, kinematic PPP, the positioning biases of ambiguity-fixed PPP were improved by about 78.2, 20.8, and 65.1 % in east, north, and up. The RMS of LEO PPP test was improved by about 23.0, 37.0, and 43.0 % for GRACE-A and GRACE-B in radial, tangential, and normal directions when AR was applied to the same data set. These results demonstrated that the SGG FCB products can be produced with high quality for users anywhere around the world to carry out ambiguity-fixed PPP solutions.     

Weak GPS signal tracking using FFT discriminator in open loop receiver

An open loop tracking architecture, which tracks GPS signals under weak and challenging conditions, is analyzed. The in-phase and quadrature-phase integration pair is regarded as a single tone complex signal. An FFT-based method is used as a frequency discriminator to estimate the Doppler frequency residual of the single tone signal. Another FFT-based method applies complex squaring to eliminate the effect of the navigation data bits polarities. The performance of the FFT-based discriminators is assessed in three criteria. Those criteria are the signal strength and dynamic range that can be tracked and the accuracy of the estimated Doppler frequency. In addition, the performance of the discriminators is analyzed to provide the theoretical and simulated peak detection probability. The results indicate that the FFT discriminator can track signals about 5 dB weaker than the signals that can be tracked by the complex squared FFT discriminator. In a quasi-static environment, the Doppler frequency residual can be assumed to be around zero, which can enable the FFT-based discriminators to track signals with approximately 2 dB less power. Moreover, the performance of the FFT-based discriminators is compared with the performance of two other frequency discriminators, namely the fast–slow and power-based. The comparison results indicate that these two frequency discriminators give higher frequency estimation accuracy, but they have a narrower dynamic range.     

Is the long-term variation of the estimated GPS differential code biases associated with ionospheric variability?

The global positioning system (GPS) differential code biases (DCB) provided by the International GNSS Service (IGS) show solar-cycle-like variation during 2002–2013. This study is to examine whether this variation of the GPS DCBs is associated with ionospheric variability. The GPS observations from low earth orbit (LEO) satellites including CHAMP, GRACE and Jason-1 are used to address this issue. The GPS DCBs estimated from the LEO-based observations at different orbit altitudes show a similar tendency as the IGS DCBs. However, this solar-cycle-like dependency is eliminated when the DCBs of 13 continuously operating GPS satellites are constrained to zero-mean. Our results thus revealed that ionospheric variation is not responsible for the long-term variation of the GPS DCBs. Instead, it is attributed to the GPS satellite replacement with different satellite types and the zero-mean condition imposed on all satellite DCBs.     

Testing block subdivision algorithms on block designs

Integrated land use–transportation models predict future transportation demand taking into account how households and firms arrange themselves partly as a function of the transportation system. Recent integrated models require parcels as inputs and produce household and employment predictions at the parcel scale. Block subdivision algorithms automatically generate parcel patterns within blocks. Evaluating block subdivision algorithms is done by way of generating parcels and comparing them to those in a parcel database. Three block subdivision algorithms are evaluated on how closely they reproduce parcels of different block types found in a parcel database from Montreal, Canada. While the authors who developed each of the algorithms have evaluated them, they have used their own metrics and block types to evaluate their own algorithms. This makes it difficult to compare their strengths and weaknesses. The contribution of this paper is in resolving this difficulty with the aim of finding a better algorithm suited to subdividing each block type. The proposed hypothesis is that given the different approaches that block subdivision algorithms take, it’s likely that different algorithms are better adapted to subdividing different block types. To test this, a standardized block type classification is used that consists of mutually exclusive and comprehensive categories. A statistical method is used for finding a better algorithm and the probability it will perform well for a given block type. Results suggest the oriented bounding box algorithm performs better for warped non-uniform sites, as well as gridiron and fragmented uniform sites. It also produces more similar parcel areas and widths. The Generalized Parcel Divider 1 algorithm performs better for gridiron non-uniform sites. The Straight Skeleton algorithm performs better for loop and lollipop networks as well as fragmented non-uniform and warped uniform sites. It also produces more similar parcel shapes and patterns.