ALBERT

Description: A novel approach for the joint retrieval of aerosol optical depth (AOD) and aerosol type, using Meteosat Second Generation – Spinning Enhanced Visible and Infrared Imagers (MSG/SEVIRI) observations in two solar channels, is presented. The retrieval is based on a Time Series (TS) technique, which makes use of the two visible bands at 0.6 μm and 0.8 μm in three orderly scan times (15 min interval between two scans) to retrieve the AOD over land. Using the radiative transfer equation for plane-parallel atmosphere, two coupled differential equations for the upward and downward fluxes are derived. The boundary conditions for the upward and downward fluxes at the top and at the bottom of the atmosphere are used in these equations to provide an analytic solution for the AOD. To derive these fluxes, the aerosol single scattering albedo (SSA) and asymmetry factor are required to provide a solution. These are provided from a set of six pre-defined aerosol types with the SSA and asymmetry factor. We assume one aerosol type for a grid of 1°×1° and the surface reflectance changes little between two subsequent observations. A k-ratio approach is used in the inversion to find the best solution of atmospheric properties and surface reflectance. The k-ratio approach assumes that the surface reflectance is little influenced by aerosol scattering at 1.6 μm and therefore the ratio of surface reflectances in the solar band for two subsequent observations can be well-approximated by the ratio of the reflectances at 1.6 μm. A further assumption is that the surface reflectance varies only slightly over a period of 30 min. The algorithm makes use of numerical minimisation routines to obtain the optimal solution of atmospheric properties and surface reflectance by selection of the most suitable aerosol type from pre-defined sets. A detailed analysis of the retrieval results shows that it is suitable for AOD retrieval over land from SEVIRI data. Six AErosol RObotic NETwork (AERONET) sites with different surface types are used for detailed analysis and 42 other AERONET sites are used for validation. From 445 collocations representing stable and homogeneous aerosol type, we find that 〉75% of the MSG-retrieved AOD at 0.6 and 0.8 μm values compare favourably with AERONET observed AOD values, within an error envelope of ± 0.05 ± 0.15 τ and a high correlation coefficient (R〉0.86). The AOD datasets derived using the TS method with SEVIRI data is also compared with collocated AOD products derived from NASA TERRA and AQUA MODIS (The Moderate-resolution Imaging Spectroradiometer) data using the Dark Dense Vegetation (DDV) method and the Deep Blue algorithms. Using the TS method, the AOD could be retrieved for more pixels than with the NASA Deep Blue algorithm. This method is potentially also useful for surface reflectance retrieval using SEVIRI observations. The current paper focuses on AOD retrieval and analysis, and the analysis and validation of reflectance will be given in a following paper.

Description: In order to quantify the aerosol impact on climate, a range of aerosol parameters are required. In this paper, two-year of ground-based observations of aerosol optical properties from an urban site in Beijing of China are assessed. The aerosol absorption coefficient (σa), scattering coefficient (σs), as well as single scattering albedo (ω) are analyzed to aid in characterizing Beijing's urban aerosol. Two-year averages (and standard deviations) for σa at 532 nm, σs at 525 nm and ω at 525 nm are 56±49 Mm−1, 288±281 Mm−1 and 0.80±0.09, respectively. Meanwhile, there is a distinct diurnal variation for σa, with its minimum occurring at approximately 14:00 to 15:00 and maximum at midnight. σs peaks in the late morning and the minimum occurs in the evening. σs in summer is higher than that in winter. ω is also higher in summer than that in winter, except before 07:00 a.m., and peaks in the early afternoon. Both σa and σs show strong dependence on local wind in all four seasons. When the wind blows from the north with low speed (0–2 m/s), the values of σa are high, and in contrast, very low with wind speeds higher than 4 m/s. When the wind blows from south with low speed (0–4 m/s), σa is intermediate. The patterns of the wind dependence of σa indicates that σa is mainly dominated by local emissions. σs displays a similar dependence on wind speed and direction to σa, except in summer. In summer, the σs value is highest when wind is from southeast with speed of 0–6 m/s. This indicates that the particle pollution resulting from regional transport is only significant in the summer season. ω also shows wind dependence to some extent though not as strong as σa or σs. Overall, the wind dependence results provide valuable information about the locations of aerosol pollution sources and suggest that the air pollution in summer is a regional problem but in other seasons it is mainly affected by local urban emissions.

Description: North China Plain (NCP) is one of the most densely populated regions in China and has experienced enormous economic growth in the past decades. Its regional trace gas pollution has also become one of the top environmental concerns in China. Measurements of surface trace gases, including O3, NOx, SO2 and CO were carried out within the HaChi (Haze in China) project at Wuqing Meteorology Station, located between 2 mega-cities (Beijing and Tianjin) in the NCP, from 9 July 2009 to 21 January 2010. Detailed statistical analyses were made in order to provide information on the levels of the measured air pollutants and their characteristics. Gaseous air pollutant concentrations were also studied together with meteorological data and satellite data to help us better understand the causes of the observed variations in the trace gases during the field campaign. In comparison to measurements from other rural and background stations in the NCP, relatively high concentrations were detected in Wuqing, presumably due to regional mixing and transport of pollutants. Local meteorology had deterministic impacts on air pollution levels, which have to be accounted for when evaluating other effects on pollutant concentrations. Trace gas concentrations showed strong dependence on wind, providing information on regional pollution characteristics. O3 mixing ratio also showed clear dependencies on temperature and relative humidity.

Description: Agricultural biomass burning (ABB) in Central and East China occurs every year from May to October and peaks in June. The biomass burning event in June 2007 was very strong. During the period from 26 May to 16 June 2007, ABB occurred mainly in Anhui, Henan, Jiangsu and Shandong provinces. A comprehensive set of aerosol optical depth (AOD) data, produced by a merger of AOD product data from the Moderate Resolution Imaging Spectroradiometer (MODIS) and the Multiangle Imaging Spectroradiometer (MIRS), is used to study the spatial and temporal distribution of agricultural biomass aerosols in Central and East China combining with ground observations from both AErosol RObotic NETwork (AERONET) and China Aerosol Remote Sensing NETwork (CARSNET) measurements. We compared merged AOD data with single-sensor single-algorithm AOD data (MODIS Dark Target AOD data, MODIS Deep Blue AOD data, SRAP-MODIS AOD data and MISR AOD data). In this comparison, we found merged AOD products can improve the quality of AOD products from single-sensor single-algorithm data sets by expanding the spatial coverage of the study area and keeping the statistical confidence in AOD parameters. There existed high correlation (0.8479) between the merged AOD data and AERONET measurements. Our merged AOD data make use of synergetic information conveyed in all of the available satellite data. The merged AOD data were used for the analysis of the biomass burning event from 26 May to 16 June 2007 together with meteorological data. The merged AOD products and the ground observations from China suggest that biomass burning in Central and East China has had great impact on AOD over China. Influenced by this ABB, the highest AOD value in Beijing on 12 June 2007 reached 5.71.

Description: The eastern China seas are one of the largest marginal seas in the world, where high primary productivity and phytoplankton blooms are often observed. However, to date, little is known about the spatial and temporal variability of phytoplankton blooms in these areas due to the difficulty of the monitoring of bloom events by field measurement. In this study, 14-yr time series of satellite ocean color data from the Sea-Viewing Wide Field-of-view Sensor (SeaWiFS) and the Moderate Resolution Imaging Spectroradiometer (MODIS) onboard the Aqua satellite have been used to investigate the seasonal and inter-annual variability and long-term changes of phytoplankton blooms in the eastern China seas. We validated and calibrated the satellite-derive chlorophyll concentration in the eastern China seas based on extensive data sets from two large cruises. Overestimation of satellite-derive chlorophyll concentration caused by high water turbidity was found to be less than 10 μg L−1. This level can be used as a safe threshold for the identification of a phytoplankton bloom in a marginal sea with turbid waters. Annually, blooms mostly occur in the Changjiang Estuary and along the coasts of Zhejiang. The coasts of the northern Yellow Sea and Bohai Sea also have high-frequency blooms. The blooms have significant seasonal variation, with most of the blooms occurring in the spring (April–June) and summer (July–September). This study revealed a doubling in bloom intensity in the Yellow Sea and Bohai Sea during the past 14 yr (1998–2011), yet surprisingly, there has been no decadal increase or decrease of bloom intensity in despite of significant inter-annual variation in the Changjiang Estuary. The time series in situ datasets show that both the nitrate and phosphate concentrations increase more than twofold from 1998 to 2005. This might be the reason for the doubling of bloom intensity in the Yellow Sea and Bohai Sea. In addition, the ENSO and PDO can affect the inter-annual variation of bloom intensity in the eastern China seas.

Description: The air-sea exchanges of CO2 in the world's 165 estuaries and 87 continental shelves are evaluated. Generally and in all seasons, upper estuaries with salinities of less than two are strong sources of CO2 (39 ± 56 mol C m−2 yr−1, negative flux indicates that the water is losing CO2 to the atmosphere); mid-estuaries with salinities of between 2 and 25 are moderate sources (17.5 ± 34 mol C m−2 yr−1) and lower estuaries with salinities of more than 25 are weak sources (8.4 ± 14 mol C m−2 yr−1). With respect to latitude, estuaries between 23.5 and 50° N have the largest flux per unit area (63 ± 101 mmol C m−2 d−1); these are followed by mid-latitude estuaries (23.5–0° S: 44 ± 29 mmol C m−2 d−1; 0–23.5° N: 39 ± 55 mmol C m−2 d−1), and then regions north of 50° N (36 ± 91 mmol C m−2 d−1). Estuaries south of 50° S have the smallest flux per unit area (9.5 ± 12 molC m−2 d−1). Mixing with low-pCO2 shelf waters, water temperature, residence time and the complexity of the biogeochemistry are major factors that govern the pCO2 in estuaries but wind speed, seldom discussed, is critical to controlling the air-water exchanges of CO2. The total annual release of CO2 from the world's estuaries is now estimated to be 0.10 PgC yr−1, which is much lower than published values mainly because of the contribution of a considerable amount of heretofore unpublished or new data from Asia and the Arctic. The Asian data, although indicating high in pCO2, are low in sea-to-air fluxes because the wind speeds are lower than previously determined values, which rely heavily on data from Europe and North America, where pCO2 is lower but wind speeds are much higher, such that the CO2 fluxes are higher than in Asia. Newly emerged CO2 flux data in the Arctic reveal that estuaries there mostly absorb, rather than release CO2. Most continental shelves, and especially those at high latitude, are under-saturated in terms of CO2 and absorb CO2 from the atmosphere in all seasons. Shelves between 0° and 23.5° S are on average a weak source and have a small flux per unit area of CO2 to the atmosphere. Water temperature, the spreading of river plumes, upwelling, and biological production seem to be the main factors in determining pCO2 in the shelves. Wind speed, again, is critical because at high latitudes, the winds tend to be strong. Since the surface water pCO2 values are low, the air-to-sea fluxes are high in regions above 50° N and below 50° S. At low latitudes, the winds tend to be weak, so the sea-to-air CO2 flux is small. Overall, the world's continental shelves absorb 0.4 PgC yr−1 from the atmosphere.

Description: Suspended sediment yields in the Upper Yangtze River and its four headwater tributaries (i.e. Jinsha, Min, Jialing and Wu) have declined significantly during recent decades. Compared with 1956–1970, mean annual suspended sediment yield during 2001–2011 was reduced by 84% in the Upper Yangtze River at Yichang, by 34% in the Jinsha at Pingshan, by 84% in the Jialing at Beibei, by 75% in the Wu at Wulong, and by 48% in the Min at Gaochang. Linking these observed decadal changes of runoff discharge and suspended sediment load to dam construction and multiple environmental rehabilitation projects (e.g. soil-water conservation, reforestation) during the past decades, it can be concluded that the construction of large dams on the main stem and major tributaries of the Upper Yangtze River has played a principal role in the reduction of fluvial suspended sediment yields, while the environment rehabilitation projects may make limited contributions to the changes in suspended sediment yields, except for the Jialing River.

Description: Groundwater overdraft has caused fast water level decline in the North China Plain (NCP) since the 1980s. Although many hydrological models have been developed for the NCP in the past few decades, most of them deal only with the groundwater component or only at local scales. In the present study, a coupled surface water–groundwater model using the MIKE SHE code has been developed for the entire alluvial plain of the NCP. All the major processes in the land phase of the hydrological cycle are considered in the integrated modeling approach. The most important parameters of the model are first identified by a sensitivity analysis process and then calibrated for the period 2000–2005. The calibrated model is validated for the period 2006–2008 against daily observations of groundwater heads. The simulation results compare well with the observations where acceptable values of root mean square error (RMSE) (most values lie below 4 m) and correlation coefficient (R) (0.36–0.97) are obtained. The simulated evapotranspiration (ET) is then compared with the remote sensing (RS)-based ET data to further validate the model simulation. The comparison result with a R2 value of 0.93 between the monthly averaged values of simulated actual evapotranspiration (AET) and RS AET for the entire NCP shows a good performance of the model. The water balance results indicate that more than 70% of water leaving the flow system is attributed to the ET component, of which about 0.25% is taken from the saturated zone (SZ); about 29% comes from pumping, including irrigation pumping and non-irrigation pumping (net pumping). Sustainable water management analysis of the NCP is conducted using the simulation results obtained from the integrated model. An effective approach to improve water use efficiency in the NCP is by reducing the actual ET, e.g. by introducing water-saving technologies and changes in cropping.

Description: Multiple-point geostatistical simulation (MPS) has recently become popular in stochastic hydrogeology, primarily because of its capability to derive multivariate distributions from a training image (TI). However, its application in three-dimensional (3-D) simulations has been constrained by the difficulty of constructing a 3-D TI. The object-based unconditional simulation program TiGenerator may be a useful tool in this regard; yet the applicability of such parametric training images has not been documented in detail. Another issue in MPS is the integration of multiple geophysical data. The proper way to retrieve and incorporate information from high-resolution geophysical data is still under discussion. In this study, MPS simulation was applied to different scenarios regarding the TI and soft conditioning. By comparing their output from simulations of groundwater flow and probabilistic capture zone, TI from both sources (directly converted from high-resolution geophysical data and generated by TiGenerator) yields comparable results, even for the probabilistic capture zones, which are highly sensitive to the geological architecture. This study also suggests that soft conditioning in MPS is a convenient and efficient way of integrating secondary data such as 3-D airborne electromagnetic data (SkyTEM), but over-conditioning has to be avoided.

Description: Uncertainty of groundwater model predictions has in the past mostly been related to uncertainty in the hydraulic parameters, whereas uncertainty in the geological structure has not been considered to the same extent. Recent developments in theoretical methods for quantifying geological uncertainty have made it possible to consider this factor in groundwater modeling. In this study we have applied the multiple-point geostatistical method (MPS) integrated in the Stanford Geostatistical Modeling Software (SGeMS) for exploring the impact of geological uncertainty on groundwater flow patterns for a site in Denmark. Realizations from the geostatistical model were used as input to a groundwater model developed from Modular three-dimensional finite-difference ground-water model (MODFLOW) within the Groundwater Modeling System (GMS) modeling environment. The uncertainty analysis was carried out in three scenarios involving simulation of groundwater head distribution and travel time. The first scenario implied 100 stochastic geological models all assigning the same hydraulic parameters for the same geological units. In the second scenario the same 100 geological models were subjected to model optimization, where the hydraulic parameters for each of them were estimated by calibration against observations of hydraulic head and stream discharge. In the third scenario each geological model was run with 216 randomized sets of parameters. The analysis documented that the uncertainty on the conceptual geological model was as significant as the uncertainty related to the embedded hydraulic parameters.