Progressive Transformation of a Flux Rope to an ICME

The solar wind conditions at one astronomical unit (AU) can be strongly disturbed by interplanetary coronal mass ejections (ICMEs). A subset, called magnetic clouds (MCs), is formed by twisted flux ropes that transport an important amount of magnetic flux and helicity, which is released in CMEs. At 1 AU from the Sun, the magnetic structure of MCs is generally modeled by neglecting their expansion during the spacecraft crossing. However, in some cases, MCs present a significant expansion. We present here an analysis of the huge and significantly expanding MC observed by the Wind spacecraft during 9 – 10 November 2004. This MC was embedded in an ICME. After determining an approximate orientation for the flux rope using the minimum variance method, we obtain a precise orientation of the cloud axis by relating its front and rear magnetic discontinuities using a direct method. This method takes into account the conservation of the azimuthal magnetic flux between the inbound and outbound branches and is valid for a finite impact parameter (i.e., not necessarily a small distance between the spacecraft trajectory and the cloud axis). The MC is also studied using dynamic models with isotropic expansion. We have found (6.2±1.5)×10²⁰ Mx for the axial flux and (78±18)×10²⁰ Mx for the azimuthal flux. Moreover, using the direct method, we find that the ICME is formed by a flux rope (MC) followed by an extended coherent magnetic region. These observations are interpreted by considering the existence of a previously larger flux rope, which partially reconnected with its environment in the front. We estimate that the reconnection process started close to the Sun. These findings imply that the ejected flux rope is progressively peeled by reconnection and transformed to the observed ICME (with a remnant flux rope in the front part).     

A Stochastic Approach to Estimate Distribution of Built-Up Area in Regions with Thick Tree Cover

Buildings and other human-made constructions have been accepted as an indicator of human habitation and are identified as built-up area. Identification of built-up area in a region and its subsequent measurement is a key step in many fields of studies like urban planning, environmental studies, and population demography. Remote sensing techniques utilising medium resolution images (e.g. LISS III, Landsat) are extensively used for the extraction of the built-up area as high-resolution images are expensive, and its processing is difficult. Extraction of built land use from medium resolution images poses a challenge in regions like Western-Ghats, North-East regions of India, and countries in tropical region, due to the thick evergreen tree cover. The spectral signature of individual houses with a small footprint are easily overpowered by the overlapping tree canopy in a medium resolution image when the buildings are not clustered. Kerala is a typical case for this scenario. The research presented here proposes a stochastic-dasymetric process to aid in the built-up area recognition process by taking Kerala as a case study. The method utilises a set of ancillary information to derive a probability surface. The ancillary information used here includes distance from road junctions, distance from road network, population density, built-up space visible in the LISS III image, the population of the region, and the household size. The methodology employs logistic regression and Monte Carlo simulation in two sub processes. The algorithm estimates the built-up area expected in the region and distributes the estimated built-up area among pixels according to the probability estimated from the ancillary information. The output of the algorithm has two components. The first component is an example scenario of the built-up area distribution. The second component is a probability surface, where the value of each pixel denotes the probability of that pixel to have a significant built-up area within it. The algorithm is validated for regions in Kerala and found to be significant. The model correctly predicted the built-up pixel count count over a validation grid of 900 m in 95.2% of the cases. The algorithm is implemented using Python and ArcGIS.     

Comparison of Vegetation Indices from Two Ground Based Sensors

The best and commonly used ground-based sensor to monitor crop growth, ASD FieldSpecPro Spectroradiometer (Analytical Spectral Devices, Boulder, CO, USA) is a passive sensor, which can be used under adequate light condition. However, now-a-days active sensors such as GreenSeeker? (GS) handheld crop response (Trimble Agriculture division, USA) are used for monitoring crop growth and are flexible in terms of timeliness and illumination conditions besides being cheaper than the ASD. Before its wide use, the suitability and accuracy of GS should be assessed by comparing the NDVI measured by this instrument with that by ASD, under diverse wheat growing conditions of India. Keeping this in view, the present experiment was undertaken with the following objectives: (1) to find out the temporal variation of NDVI measured both by ASD and GS treatments, (2) to find out relationship between the NDVI measured through ASD and GS and, (3) to evaluate the suitability of GS for NDVI measurements. It was observed that the numerical value of NDVI as measured by GS was always significantly (P < 0.05) lower than that measured by ASD for all the experiments under study. The NDVI-ASD and NDVI-GS were significantly positively correlated (P < 0.01) with the correlation coefficients being +0.94, +0.88 and +0.87 for irrigation and nitrogen experiment, irrigation and cultivars experiment, and tillage, residue and nitrogen experiments, respectively. Further, the regression equation developed between the NDVI-ASD and NDVI-GS: [NDVI-GS = 1.070 × (NDVI-ASD ? 0.292] can be successfully used to compute the NDVI of ASD from that computed by GS.     

Estimation of Partial Pressure of Carbon Dioxide and Air-Sea Fluxes in Hooghly Estuary Based on In Situ and Satellite Observations

An empirical model is developed and used with remotely sensed predictors: sea surface temperature (SST) and chlorophyll-a concentration (Chl-a), to compute surface water partial pressure of carbon dioxide (pCO_2w) and air-sea fluxes of CO₂ in the Hooghly estuary and its adjacent coastal oceans. In situ observations used here were based on measurements carried out in this region during winter and summer periods in 2008. The estimated pCO_2w compares well with the in situ observations at root mean square error ±18 μatm. In winter, estimated pCO_2w ranges between 320 and 500 μatm with large values (>400 μatm) on the south-western and south-eastern flanks of the coastal domain and lower values (340–375 μatm) on the main-channel. In summer, it remained spatially uniform at 450 μatm. Extrapolation of the results over the study region based on the Moderate Imaging Specroradiometer (MODIS) measured SST and Chl-a suggests that the region is a strong source of atmospheric CO₂ during the summer with net release of 0.095 Tg C year^?1 (equivalent to mean flux of 90 molC m^?2 year^?1) and is a weak source during the winter with net release of 0.006 Tg C yr^?1 (0.5 molC m^?2 year^?1) from the geographical extent of 6000 Km² area.     

重力场数学拟合模型解的影响因素探究

为了提高地球重力场模型不适定方程求解的精度,该文采用谱分析方法从级数展开阶数、数据采样率及数据缺失量3个方面探索影响数学拟合效果的根本因素:从常用的三角级数及勒让德级数模型出发,引出重力场拟合模型球谐函数模型,观察在改变级数展开阶数、数据采样率及数据缺失量等情况下所对应设计矩阵谱结构的变化,并从微观上研究影响误差分配的有关因素及最小奇异值对误差的决定性作用,为探求重力场模型解不准的原因及实现更高精度的全球重力场模型的建立提供参考。     

Nonlinear least-squares method via an isomorphic geometrical setup

The resolution of a nonlinear parametric adjustment model is addressed through an isomorphic geometrical setup with tensor structure and notation, represented by a u-dimensional “model surface” embedded in a flat n-dimensional “observational space”. Then observations correspond to the observational-space coordinates of the pointQ, theu initial parameters correspond to the model-surface coordinates of the “initial” pointP, and theu adjusted parameters correspond to the model-surface coordinates of the “least-squares” point . The least-squares criterion results in a minimum-distance property implying that the vector Q must be orthogonal to the model surface. The geometrical setup leads to the solution of modified normal equations, characterized by a positive-definite matrix. The latter contains second-order and, optionally, thirdorder partial derivatives of the observables with respect to the parameters. This approach significantly shortens the convergence process as compared to the standard (linearized) method.