An Analytical Model For Threshold Crossing Rates Of Concentration Fluctuations In Dispersing Plumes

A simple analytical model is developed for the meanupcrossing rate of plume concentration fluctuations assuming that thisprocess can be well approximated by a lognormal process. The resultingexpression requires only the specification of the in-plume fluctuationintensity and in-plume Taylor micro-time scale and, hence, does notexplicitly involve the joint probability density function of theconcentration and its derivative. The analytical model provides agood fit to some field measurements of the mean upcrossing rate ina dispersing plume.     

A model of impulsive loop flares revisited

A critical examination of the components of the recent impulsive loop flare model of Takakura is made. It is found that his analysis of the stability of the electron distribution resulting from anomalous heat conduction is in error and electron plasma waves would not be excited. Rather, in the regions where the electron/proton temperature ratioT _e/T _i 10, electrostatic ion-cyclotron waves would be excited and in the regions whereT _e 10, ion-acoustic waves would be excited. Ratios ofT _e/T _i 10 occur only in the late time development behind the conduction fronts. Since the anomalous resistivity due to electrostatic ion-cyclotron waves is fortuitously about 70% of the one used by Takakura, the general development will follow closely the one calculated by him. Because the anomalous resistivity due to ion-acoustic waves is about 95 times the one used by Takakura, the development in the parts of the loop whereT _e/T _i 10 for late times would be altered considerably.Also Guest Worker at NOAA Space Environment Laboratory, Boulder, Colorado, U.S.A.     

Study on uncertainty of inversion and resolution limit in the presence of noise

This study analyzed the uncertainty of inversion and the resolution limit in the presence of noise by means of statistical experiments. The exhaustive method is adopted to obtain the global optimal solution in each experiment. We found that even with small level of noise, solutions fluctuate in a large range for the thin bed. The distribution of solutions in the presence of noise is closely related to the spread of the cost function in the absence of noise. As a result, the area of a certain neighborhood around the true solution on the spread of the cost function in the absence of noise is used to evaluate the uncertainty of inversion and the resolution limit in the presence of noise. In the case that the SNR (signal-to-noise ratio) is 5 in this study, solutions focus around the true solution with a very small uncertainty only when the bed thickness is greater than the reciprocal of the double predominant frequency of the convoluting wavelet.     

An inter-comparison of tidal solutions computed with a range of unstructured grid models of the Irish and Celtic Sea Regions

Three finite element codes, namely TELEMAC, ADCIRC and QUODDY, are used to compute the spatial distributions of the M₂, M₄ and M₆ components of the tide in the sea region off the west coast of Britain. This region is chosen because there is an accurate topographic dataset in the area and detailed open boundary M₂ tidal forcing for driving the model. In addition, accurate solutions (based upon comparisons with extensive observations) using uniform grid finite difference models forced with these open boundary data exist for comparison purposes. By using boundary forcing, bottom topography and bottom drag coefficients identical to those used in an earlier finite difference model, there is no danger of comparing finite element solutions for “untuned unoptimised solutions” with those from a “tuned optimised solution”. In addition, by placing the open boundary in all finite element calculations at the same location as that used in a previous finite difference model and using the same M₂ tidal boundary forcing and water depths, a like with like comparison of solutions derived with the various finite element models was possible. In addition, this open boundary was well removed from the shallow water region, namely the eastern Irish Sea where the higher harmonics were generated. Since these are not included in the open boundary, forcing their generation was determined by physical processes within the models. Consequently, an inter-comparison of these higher harmonics generated by the various finite element codes gives some indication of the degree of variability in the solution particularly in coastal regions from one finite element model to another. Initial calculations using high-resolution near-shore topography in the eastern Irish Sea and including “wetting and drying” showed that M₂ tidal amplitudes and phases in the region computed with TELEMAC were in good agreement with observations. The ADCIRC code gave amplitudes about 30 cm lower and phases about 8° higher. For the M₄ tide, in the eastern Irish Sea amplitudes computed with TELEMAC were about 4 cm higher than ADCIRC on average, with phase differences of order 5°. For the M₆ component, amplitudes and phases showed significant small-scale variability in the eastern Irish Sea, and no clear bias between the models could be found. Although setting a minimum water depth of 5 m in the near-shore region, hence removing wetting and drying, reduced the small-scale variability in the models, the differences in M₂ and M₄ tide between models remained. For M₆, a significant reduction in variability occurred in the eastern Irish Sea when a minimum 5-m water depth was specified. In this case, TELEMAC gave amplitudes that were 1 cm higher and phases 30° lower than ADCIRC on average. For QUODDY in the eastern Irish Sea, average M₂ tidal amplitudes were about 10 cm higher and phase 8° higher than those computed with TELEMAC. For M₄, amplitudes were approximately 2 cm higher with phases of order 15° higher in the northern part of the region and 15° lower in the southern part. For M₆ in the north of the region, amplitudes were 2 cm higher and about 2 cm lower in the south. Very rapid M₆ tidal-phase changes occurred in the near-shore regions. The lessons learned from this model inter-comparison study are summarised in the final section of the paper. In addition, the problems of performing a detailed model–model inter-comparison are discussed, as are the enormous difficulties of conducting a true model skill assessment that would require detailed measurements of tidal boundary forcing, near-shore topography and precise knowledge of bed types and bed forms. Such data are at present not available.     

Tectonics of southwestern Mexico,isotopic evidence,nuclear Central America,Late Cretaceous break up

In the paleogeographic reconstruction of Mexico and northern Central America, an ever-increasing amount of evidence shows that the entire region is a collage of suspect terranes transported from abroad, whose timing and sense of motion are now beginning to be understood. Among these, the Chortis block (nuclear Central America) and the Baja California Peninsula have been proposed as pieces of continent separated from the Pacific coast of southwestern Mexico, that have moved either southeastward by the Farallon plate or northwestward by the Kula plate. Previous studies mainly confined to the northern margin of the Chortis block, confirmed a left-lateral displacement of 130 km in Neogene time. Further studies made northwestward along the Mexican coast provided a better understanding of magmatic and metamorphic processes in the area, and suggested times of detachment increased to 30 Ma, 40 Ma, and 66 Ma. The pre-detachment westernmost position of the block has changed, depending on the model chosen, from Puerto Vallarta and beyond, to the current position. Here we show that the isotopic mineral ages from coastal granites along the coast from Puerto Vallarta, Jalisco (80 Ma) to Puerto Angel, Oaxaca (11 Ma) record systematic decrease of cooling ages from NW to SE. This pattern is interpreted to result from the progressive uplift of rocks exposed at the present-day coast in that direction, such uplift occurred in response to the development of the Middle America Trench at the newly formed continental margin when the Chortis block was sliding at an average rate of 1.5 cm/year in a sinistral sense to its present position. Our results also constrain the position of the Kula-Farallon spreading axis north of Puerto Vallarta. These observations led us to conclude that several indicators point to this time and region for the onset of strike-slip drifting of the Chortis block toward its current position. Here, we also present several view points in terms of other possilble interpretations to different tectonic, geologic and isotopic data sets published recently by different authors.     

Integration of reflection seismic and sediment grain-size data from Lake Khubsugul (Northern Mongolia): a reply to Prokopenko and Kendall

Prokopenko and Kendall (J Paleolimnol doi:, 2008) criticise the work presented in Fedotov et al. (J Paleolimnol 39:335–348, 2008), and instead propose an alternative interpretation for the grain-size evolution recorded in the KDP-01 core, retrieved from the central part of Lake Khubsugul. Their interpretation is based (i) on a seismic-stratigraphic re-interpretation of sparker seismic profile khub012 (which they copied from Fedotov et al. (EOS Trans 87:246–250, 2006)), (ii) on the presupposition that changes in lake level are the dominant control on facies distribution in Lake Khubsugul, and (iii) on the invalidation of our age-depth model. In this reply to their comment, we demonstrate that they interpreted seismic artefacts and geometries caused by changes in profile orientation as true stratigraphic features and that the lake-level reconstruction they derive from this interpretation is therefore incorrect. We also demonstrate that their grain-size predictions, which they consider to be predominantly driven by changes in lake level, are inconsistent with the measured sulphate concentration, which is a demonstrated proxy of lake level in Lake Khubsugul, and with the measured grain-size record. Finally, we point out that even if there would be a problem with the age-depth model, this problem would not affect the part of the sedimentary sequence discussed in Fedotov et al. (J Paleolimnol 39:335–348, 2008).     

Artificial light at night: implications for early life stages development in four temperate freshwater fish species

Flora and fauna have both evolved under a natural cycle of light and dark. But especially in urban areas, the night is now increasingly disturbed by artificial light. Many traits and behaviours in fish are triggered by a circadian clock, for example hatching and swim bladder inflation, which predominantly take place at dusk or night. As lighting becomes brighter and extends farther into rural areas, the distinction between day and night becomes increasingly blurred. Therefore, the loss of diurnal trigger by artificial light at night was hypothesized having deleterious effects on these traits and impact fish reproduction. To assess these effects, eggs of four native freshwater fishes, Eurasian perch Perca fluviatilis, roach Rutilus rutilus, bleak Alburnus alburnus and chub Leuciscus cephalus, were incubated under two different light conditions: a photoperiod of 14 h light:10 h darkness (LD) and continuous illumination (LL). The time to hatch and swim bladder inflation was recorded. The species showed inconsistent reactions to the light treatments. In roach and bleak, the time to 50% hatch was longer in LL, whereas continuous lighting had an accelerating effect in chub. Incubation in LL elongated the hatching period in perch and roach and, in perch, the onset of darkness seemed to trigger hatching. The swim bladder inflation was significantly promoted by continuous light in chub and bleak but was not affected in roach. In conclusion, nocturnal artificial illumination could have an effect on hatching and initial swim bladder filling by masking the day–night-change and thereby diminish the trigger effect. However, the reactions were species specific and the increase in variation indicated a lack of diurnal triggering, whilst a general deleterious effect of artificial light at night has not been identified on early life stages.     

The role of prior model calibration on predictions with ensemble Kalman filter

This paper, based on a real world case study (Limmat aquifer, Switzerland), compares inverse groundwater flow models calibrated with specified numbers of monitoring head locations. These models are updated in real time with the ensemble Kalman filter (EnKF) and the prediction improvement is assessed in relation to the amount of monitoring locations used for calibration and updating. The prediction errors of the models calibrated in transient state are smaller if the amount of monitoring locations used for the calibration is larger. For highly dynamic groundwater flow systems a transient calibration is recommended as a model calibrated in steady state can lead to worse results than a noncalibrated model with a well-chosen uniform conductivity. The model predictions can be improved further with the assimilation of new measurement data from on-line sensors with the EnKF. Within all the studied models the reduction of 1-day hydraulic head prediction error (in terms of mean absolute error [MAE]) with EnKF lies between 31% (assimilation of head data from 5 locations) and 72% (assimilation of head data from 85 locations). The largest prediction improvements are expected for models that were calibrated with only a limited amount of historical information. It is worthwhile to update the model even with few monitoring locations as it seems that the error reduction with EnKF decreases exponentially with the amount of monitoring locations used. These results prove the feasibility of data assimilation with EnKF also for a real world case and show that improved predictions of groundwater levels can be obtained.