On solution of the inverse problem for confined aquifer flow via maximum likelihood

Joint estimation of transmissivity (T) and storativity (S) in a confined aquifer is done via maximum likelihood (ML). The differential equation of groundwater flow is discretized by the finite-element method, leading to equation _t+x _t=u _t. Elements of matrices and , as well as estimated covariance matrix of noise termu _t, are functions of T and S. By minimizing the negative loglikelihood function corresponding to discretized groundwater flow equation with respect to T and S, ML estimators are obtained. The ML approach is found to yield accurate estimates of T and S (within 9 and 10% of their actual values, respectively) and showed quadratic convergence in Newton's search technique. Prediction of aquifer response, using ML estimators, results in estimated piezometric heads accurate to ±0.5 m from their actual, exact values. Statistical properties of ML estimators are derived and some basic results for statistical inference are given.     

A geostatistically based approach to the identification of aquifer transmissivities in Yolo Basin, California

A geostatistically based approach is developed for the identification of aquifer transmissivities in Yolo Basin, California. The approach combines weighted least-squares with universal kriging and cokriging techniques in an overall scheme that (1)considers a prioriknown information on aquifer transmissivity and specific capacities of wells, (2)considers uncertainties in water level and transmissivity data, and (3)estimates the reliability of the generated transmissivity values. Minimization of a global least-squares function that incorporates calibration and plausibility criteria leads to a transmissivity map that shows a good agreement with pumping-test results.     

Nest microhabitat selection in larks from a European semi-arid shrub-steppe: the role of sunlight and predation

We analysed nest site selection in two shrub-steppe larks, the lesser short-toed larkCalandrella rufescensand the Thekla larkGalerida theklae. Both species nest on the ground next to a bush preferably facing north or north-east. Nestling growth and fledging weight are greater when the nest faces these directions. Both species prefer to nest next to camephytes with little foliage, and the lesser short-toed lark chooses sites with high visibility. We conclude that sunlight and nest predation are fundamental selective pressures determining nest microhabitat selection in larks in semi-arid areas.     

Model prediction capacity of ephemeral gully evolution in conservation tillage systems

Ephemeral gully (EG) erosion has an important impact on agricultural soil losses and increases field surface hydrology connectivity and transport of pollutants to nearby water bodies. Watershed models including an EG component are scarce and not yet properly evaluated. The objective of this study is to evaluate the capacity of one such tool, AnnAGNPS, to simulate the evolution of two EG formed in a conservation tillage system. The dataset for model testing included runoff measurements and EG morphological characteristics during 3 years. Model evaluation focused on EG evolution of volume, width, and length model outputs, and included calibration and testing phases and a global sensitivity analysis (GSA). While the model did not fully reproduce width and length, the model efficiency to simulate EG volume was satisfactory for both calibration and testing phases, supporting the watershed management objectives of the model. GSA revealed that the most sensitive factors were EG depth, critical shear stress, headcut detachment exponent coefficient b, and headcut detachment leading coefficient a. For EG outputs the model was additive, showing low sensitivity to interactions between the inputs. Prediction of EG spatial evolution on conservation tillage systems requires improved development of gully erosion components, since many of the processes were developed originally for traditional tillage practices or larger channel systems. Our results identify the need for future research when EG form within conservation tillage systems, in particular to study gully headcut, soil erodibility, and width functions specific to these practices.     

Paleosecular variation and absolute geomagnetic paleointensity records retrieved from the Early Cretaceous Posadas Formation (Misiones,Argentina)

The Early Cretaceous may be considered a key period for understanding the evolution of the Earth’s magnetic field. Some still unsolved problems are related to the mode of paleosecular variation (PSV) of the Earth’s magnetic field before and during the Cretaceous Normal Superchron. We report here a detailed rock-magnetic, paleomagnetic and paleointensity investigation from 28 lava flows (331 standard paleomagnetic cores) collected in the Argentinean part of the Parana Flood Basalts (Formation Posadas) in order to contribute to the study of PSV during the early Cretaceous and to obtain precise Cretaceous paleomagnetic pole positions for stable South America. The average paleofield direction is precisely determined from 26 sites, which show small within-site dispersion and high directional stability. Five sites show evidences for the self-reversal of thermoremanent magnetization. 23 sites yielded normal polarity magnetization and only 3 are reversely magnetized. Moving windows averages were used to analyze the sequential variation of virtual geomagnetic pole’s (VGP) axial positions. Interestingly, the axial average VGP path traces an almost complete cycle around the geographical pole and passes near the location of all previously published Paraná Magmatic Province poles. Both paleomagnetic poles and average VGP paths are significantly different from the pole position suggested by fixed hotspot reconstructions, which may be due to true polar wander or the hotspot motion itself. Only 15 samples from 5 individual basaltic lava flows, yielded acceptable paleointensity estimates. The site mean paleointensities range from 25.2 ± 2.2 to 44.0 ± 2.2 μT. The virtual dipole moments (VDMs) range from 4.8 to 9.9 × 10²² Am². This correspond to a mean value of 7.7 ± 2.1 × 10²² Am² which is 96% of the present day geomagnetic field strength. These intensities agree with the relatively high values already reported for Early Cretaceous, which are consistent with some inferences from computer simulations previously published.     

Revisit the calibration errors on experimental slant total electron content (TEC) determined with GPS

The calibration errors on experimental slant total electron content (TEC) determined with global positioning system (GPS) observations is revisited. Instead of the analysis of the calibration errors on the carrier phase leveled to code ionospheric observable, we focus on the accuracy analysis of the undifferenced ambiguity-fixed carrier phase ionospheric observable determined from a global distribution of permanent receivers. The results achieved are: (1) using data from an entire month within the last solar cycle maximum, the undifferenced ambiguity-fixed carrier phase ionospheric observable is found to be over one order of magnitude more accurate than the carrier phase leveled to code ionospheric observable and the raw code ionospheric observable. The observation error of the undifferenced ambiguity-fixed carrier phase ionospheric observable ranges from 0.05 to 0.11 total electron content unit (TECU) while that of the carrier phase leveled to code and the raw code ionospheric observable is from 0.65 to 1.65 and 3.14 to 7.48 TECU, respectively. (2) The time-varying receiver differential code bias (DCB), which presents clear day boundary discontinuity and intra-day variability pattern, contributes the most part of the observation error. This contribution is assessed by the short-term stability of the between-receiver DCB, which ranges from 0.06 to 0.17 TECU in a single day. (3) The remaining part of the observation errors presents a sidereal time cycle pattern, indicating the effects of the multipath. Further, the magnitude of the remaining part implies that the code multipath effects are much reduced. (4) The intra-day variation of the between-receiver DCB of the collocated stations suggests that estimating DCBs as a daily constant can have a mis-modeling error of at least several tenths of 1 TECU.     

Crustal xenoliths from Calbuco Volcano,Andean Southern Volcanic Zone: implications for crustal composition and magma-crust interaction

Crustal xenoliths in the 1961 andesite flow of Calbuco Volcano, in the southern Southern Volcanic Zone (SSVZ) of the Andes, consist predominantly of pyroxene granulites and hornblende gabbronorites. The granulites contain plagioclase+pyroxene+magnetite±amphibole, and have pronounced granoblastic textures. Small amounts of relict amphibole surrounded by pyroxene-plagioclase-magnetite-glass symplectites are found in some specimens. These and similar textures in the gabbronorites are interpreted as evidence of dehydration melting. Mineral and bulk rock geochemical data indicate that the granulites are derived from an incompatible trace element depleted basaltic protolith that underwent two stages of metamorphism: a moderate pressure, high temperature stage accompanied by melting and melt extraction from some samples, followed by thermal metamorphism after entrainment in the Calbuco andesite lavas. High _Nd ^T values (+4.0 to +8.6), Nd-isotope model ages of 1.7–2.0 Ga, and trace element characteristics like chondrite normalized La/Yb< and La/Nb1 indicate that the protoliths were oceanic basalts. Similar oceanic metabasalts of greenschist to amphibolite facies are found in the Paleozoic metamorphic belt that underlies the Chilean coastal ranges. Mineral and bulk rock compositions of the gabbronorite xenoliths indicate that they are cognate, crystallizing from the basaltic andesite magma at Calbuco. Crystallization pressures for the gabbros based on total Al contents in amphibole are 6–8 kbar. These pressures point to middle to lower crustal storage of the Calbuco magma. Neither granulite nor gabbro xenoliths have the appropriate geochemical characteristics to be contaminants of Calbuco andesites, although an ancient sedimentary contaminant is indicated by the lava compositions. The presence of oceanic metabasaltic xenoliths, together with the sedimentary isotopic imprint, suggests that the lower crust beneath the volcano is analogous to the coastal metamorphic belt, which is an accretionary complex of intercalated basalts and sediments that formed along the Paleozoic Gondwanan margin. If this is the case, the geochemical composition of the lower and middle crust beneath the SSVZ is significantly different from that of most recent SSVZ volcanic rocks.