Insights into the Microwave Instruments Onboard the Fengyun 3DSatellite: Data Quality and Assimilation in theMet Office NWP System

正1Met Office, Exeter EX1 3PB, UK2National Satellite Meteorological Center, China Meteorological Administration, Beijing 100081, China     

Defining the shelf edge and the three‐dimensional shelf edge to slope facies variability in shelf‐edge deltas

Shelf‐edge deltas play a critical role in shelf‐margin accretion and deepwater sediment delivery, yet much remains to be understood about the detailed linkage between shelf edge and slope sedimentation. The shelf edge separates the flat‐lying shelf from steeper slope regions, and is observable in seismic data and continuous outcrops; however, it is commonly obscured in non‐continuous outcrops. Defining this zone is essential because it segregates areas dominated by shelf currents from those governed by gravity‐driven processes. Understanding this linkage is paramount for predicting and characterizing associated deepwater reservoirs. In the Tanqua Karoo Basin, the Permian Kookfontein Formation shelf‐slope clinothems are well‐exposed for 21 km along depositional strike and dip. Two independent methods identified the shelf‐edge position, indicating that it is defined by: (i) a transition from predominantly shelf‐current to gravitational deposits; (ii) an increase in soft‐sediment deformation; (iii) a significant gradient increase; and (iv) clinothem thickening. A quantitative approach was used to assess the impact of process‐regime variability along the shelf edge on downslope sedimentation. Facies proportions were quantified from sedimentary logs and photographic panels, and integrated with mapped key surfaces to construct a stratigraphic grid. Spatial variability in facies proportions highlights two types of shelf‐edge depositional zones within the same shelf‐edge delta. Where deposition occurred in fluvial‐dominated zones, the slope is sand rich, channelized with channels widening downslope, and rich in collapse features. Where deltaic deposits indicate considerable tidal reworking, the deposits are thin and pinch‐out close to the shelf edge, and the slope is sand poor and lacks channelization. Amplification of tidal energy, and decrease in fluvial drive on the shelf, coincides with a decrease in mouth bar and shelf‐edge collapse, and a lack of channelization on the slope. This analysis suggests that process‐regime variability along the shelf edge exercised significant control on shelf‐edge progradation, slope channelization and deepwater sediment delivery.     

Simulating Magnetized Jets

A suitable model for the macroscopic behavior of accretion disk-jet systems is provided by the equations of MagnetoHydroDynamics (MHD). These equations allow us to perform scale-encompassing numerical simulations of multidimensional nonlinear magnetized plasma flows. For that purpose, we continue the development and exploitation of the Versatile Advection Code (VAC) along with its recent extension which employs dynamically controlled grid adaptation. In the adaptive mesh refinement AMRVAC code, modules for simulating any-dimensional special relativistic hydro- and magnetohydrodynamic problems are currently operational. Here, we review recent 3D MHD simulations of fundamental plasma instabilities, relevant when dealing with cospatial shear flow and twisted magnetic fields. Such magnetized jet flows can be susceptible to a wide variety of hydro (e.g. Kelvin-Helmholtz) or magnetohydrodynamic (e.g. current driven kink) instabilities. Recent MHD computations of 3D jet flows have revealed how such mutually interacting instabilities can in fact aid in maintaining jet coherency. Another breakthrough from computational magnetofluid modeling is the demonstration of continuous, collimated, transmagnetosonic jet launching from magnetized accretion disks. Summarizing, MHD simulations are rapidly gaining realism and significantly advance our understanding of nonlinear astrophysical magnetofluid dynamics.     

On the interest of using field primary production data to calibrate phytoplankton rate processes in ecosystem models

In many ecosystem models based on empirical formulations, parameters generally are calibrated in order to achieve the best fit between measured and simulated chlorophyll a standing stocks. An accurate calibration of rate processes as primary production rarely is taken into account. In this paper, we test the usefulness of calibration of phytoplankton photosynthetic processes in an ecosystem model using field primary production data. We used 18 months of photosynthetic process data from the Baie des Veys ecosystem (Normandy, France). Five empirical formulations of photosynthesis–irradiance curve models amongst the most widely used were tested. In each formulation, the variability of photosynthetic parameters (i.e. the light-saturated rate of photosynthesis (P_max^B) and the initial slope of the photosynthesis–light curve (α^B)) was considered depending on environmental factors (temperature and nutrient availability). The fit of the five equations as well as the calibration of parameters on field measurements (i.e. the light-saturated rate of photosynthesis (P_ref^B), the initial slope of the photosynthesis–light curve (α_ref^B), the half-saturation constant for nitrogen (K_N) and silicates uptake (K_Si), and the coefficient in the exponential thermal effect (K_T)) was performed using the whole available data set of P vs. E curves (n = 143, P vs. E curves). Then, the Smith formulation allowing the best simulation of the Baie des Veys primary production and corresponding parameters were introduced in an ecosystem box model. This formulation led directly to a satisfactory representation of the Baie des Veys phytoplankton dynamics without additional calibration. Results obtained were compared with a more classical approach in which ecosystem models were calibrated using published values of parameters. This comparison showed that for the two years studied, annual primary production estimated through the ecosystem model was 13% and 26% higher with our approach than with the more classical approach. This work emphasizes the importance of accurately representing rate processes in ecosystem models in order to adequately simulate production as well as standing stocks.     

Estimating boundary currents from satellite altimetry: A case study for the east coast of India

We present a methodology to derive surface geostrophic current from a newly released altimetric sea-level data set. TOPEX/Poseidon data were first completely reprocessed from Geophysical Data Records using new algorithms accommodating marginal seas and coastal conditions. The methodology applied to the reprocessed data essentially consists of a smoothing of the raw along-track coastal altimetric data at scales at which the geostrophic equilibrium holds. This was reduced to a computational procedure using a set of objective criteria. We have applied the method to the East India Coastal Current (EICC) at the western boundary of the Bay of Bengal. This paper first examines the quality of the new data set, which compares well with tide-gauge data; the current we derived is consistent with independent estimates. Our methodology reveals the full spectrum of the along-shore current, ranging from intra-seasonal to inter-annual time scales, from the deep ocean to the shelf-break area where the EICC exists. The algorithm can be applied to any coastal region where an order of the Rossby radius can be defined, and it therefore opens up bright prospects for mapping the variability of other boundary-current systems in the world ocean from altimetry.     

Exploring Geological and Topographical Controls on Low Flows with Hydrogeological Models

This study investigates how catchment properties influence low-flow dynamics. With 496 synthetic models composed of a bedrock and an alluvial aquifer, we systematically assess the impact of the hydraulic conductivity of both lithologies, of the hillslope and of the river slope on catchment dynamics. The physically based hydrogeological simulator HydroGeoSphere is employed, which allows obtaining a range of low-flow indicators. The hydraulic conductivity of the bedrock K_bedrock, a proxy for transmissivity, is the only catchment property exerting an overall control on low flows and explains 60% of the variance of Q95/Q50. The difference in dynamics of catchments with same K_bedrock depends on hillslope gradients and the alluvial aquifer properties. The buffering capacity of the bedrock is mainly related to K_bedrock and the hillslope gradient. We thus propose the dimensionless bedrock productivity index (BPI) that combines these characteristics with the mean net precipitation. For bedrock only models, the BPI explains 82% of the variance of the ratio of Q95 to mean net precipitation. The alluvial aquifer can significantly influence low flows when the bedrock productivity is limited. Although our synthetic catchment setup is simple, it is far more complex than the available analytical approaches or conceptual hydrological models. The direct application of the results to existing catchments requires nevertheless careful consideration of the local geological topographic and climatic conditions. This study provides quantitative insight into the complex interrelations between geology, topography and low-flow dynamics and challenges previous studies which neglect or oversimplify geological characteristics in the assessment of low flows.     

Gauging extreme floods on YouTube: application of LSPIV to home movies for the post‐event determination of stream discharges

Movies taken by witnesses of extreme flood events are increasingly available on video sharing websites. They potentially provide highly valuable information on flow velocities and hydraulic processes that can help improve the post‐flood determination of discharges in streams and flooded areas. We investigated the troubles and potential of applying the now mature large‐scale particle image velocimetry (LSPIV) technique to such flood movies that are recorded under non‐ideal conditions. Processing was performed using user‐friendly, free software only, such as Fudaa‐LSPIV. Typical issues related to the image processing and to the hydrological analysis are illustrated using a selected example of a pulsed flash‐flood flow filmed in a mountainous torrent. Simple corrections for lens distortion (fisheye) and limited incoherent camera movement (shake) were successfully applied, and the related errors were reduced to a few percents. Testing the different image resolution levels offered by YouTube showed that the difference in time‐averaged longitudinal velocity was less than 5% compared with full resolution. A limited number of GRPs, typically 10, is required, but they must be adequately distributed around the area of interest. The indirect determination of the water level is the main source of uncertainty in the results, usually much more than errors because of the longitudinal slope and waviness of the free‐surface of the flow. The image‐based method yielded direct discharge estimates of the base flow between pulses, of the pulse waves, and of the time‐averaged flow over a movie sequence including a series of five pulses. A comparison with traditional indirect determination methods showed that the critical‐depth method may produce significantly biassed results for such a fast, unsteady flow, while the slope‐area method seems to be more robust but would overestimate the time‐averaged flow rate if applied to the high‐water marks of a pulsed flow. Copyright © 2015 John Wiley & Sons, Ltd.     

MODIS Derived Equilibrium Line Altitude Estimates for Purogangri Ice Cap,Tibetan Plateau,and their Relation to Climatic Predictors (2001–2012)

The variation of the equilibrium line altitude can be used as an indicator for glacier mass balance variability. Snow lines at the end of the ablation period are suitable proxies for the annual equilibrium line altitude on glaciers. We investigate snow lines at Purogangri ice cap on the central Plateau in order to study the interannual variability of glacier mass balance. Datasets of the daily Moderate Resolution Imaging Spectroradiometer snow product MOD10A1 were used to infer transient snow line variability during 2001–2012 and to derive regional‐scale, annual equilibrium line altitude. The Moderate Resolution Imaging Spectroradiometer snow albedo embedded within the snow product was compared with high‐resolution Landsat imagery. An albedo threshold was established to differentiate between ice and snow and the 13th percentile of the altitudes of snow‐covered pixels was chosen to represent the snow line altitude. The second maximum of the snow line altitudes in the ablation period was taken as a proxy for the annual equilibrium line altitude. A linear correlation analysis was carried out (1) between interannual variability of the equilibrium line altitude at Purogangri ice cap and various climate elements derived from the High Asia Reanalysis, and (2) between interannual variability of the equilibrium line altitude and the circulation indices North Atlantic Oscillation and Indian Summer Monsoon. Results suggest that air temperature and meridional wind speed above ground in July, as well as the lower tropospheric zonal wind in June and August play a crucial role in the development of the annual equilibrium line altitude.