Search

Satellite observations from Clouds and the Earth's Radiant Energy System (CERES) radiometers have produced over two decades of world-class data documenting time–space variations in Earth's top-of-atmosphere (TOA) radiation budget. In addition to energy exchanges among Earth and space, climate studies require accurate information on radiant energy exchanges at the surface and within the atmosphere. The CERES Cloud Radiative Swath (CRS) data product extends the standard Single Scanner Footprint (SSF) data product by calculating a suite of radiative fluxes from the surface to TOA at the instantaneous CERES footprint scale using the NASA Langley Fu–Liou radiative transfer model. Here, we describe the CRS flux algorithm and evaluate its performance against a network of ground-based measurements and CERES TOA observations. CRS all-sky downwelling broadband fluxes show significant improvements in surface validation statistics relative to the parameterized fluxes on the SSF product, including a ∼30%–40% (∼20%) reduction in SW↓ (LW↓) root-mean-square error (RMSΔ), improved correlation coefficients, and the lowest SW↓ bias over most surface types. RMSΔ and correlation statistics improve over five different surface types under both overcast and clear-sky conditions. The global mean computed TOA outgoing LW radiation (OLR) remains within <1% (2–3 W m−2) of CERES observations, while the global mean reflected SW radiation (RSW) is excessive by ∼3.5% (∼9 W m−2) owing to cloudy-sky computation errors. As we highlight using data from two remote field campaigns, the CRS data product provides many benefits for studies requiring advanced surface radiative fluxes. [ABSTRACT FROM AUTHOR]

Aerosol optical depths (AODs) used for the Edition 4.1 Clouds and the Earth's Radiant Energy System (CERES) Synoptic 1 ∘ (SYN1deg) product are evaluated. AODs are derived from Moderate Resolution Imaging Spectroradiometer (MODIS) observations and assimilated by an aerosol transport model (the Model for Atmospheric Transport and Chemistry – MATCH). As a consequence, clear-sky AODs closely match with those derived from MODIS instruments. AODs under all-sky conditions are larger than AODs under clear-sky conditions, which is supported by ground-based AErosol RObotic NETwork (AERONET) observations. When all-sky MATCH AODs are compared with Modern-Era Retrospective analysis for Research and Applications (Version 2; MERRA-2) AODs, MATCH AODs are generally larger than MERRA-2 AODs, especially over convective regions (e.g., the Amazon, central Africa, and eastern Asia). This variation is largely due to the differing methods of assimilating the MODIS AOD data product and the use of quality flags in our assimilation. Including AODs with larger retrieval uncertainty makes AODs over the convective regions larger. When AODs are used for clear-sky irradiance computations and computed downward shortwave irradiances are compared with ground-based observations, the computed instantaneous irradiances are 1 %–2 % larger than observed irradiances. The comparison of top-of-atmosphere clear-sky irradiances with those derived from CERES observations suggests that AODs used for surface radiation observation sites are 0.01–0.03 larger, which is within the uncertainty of instantaneous MODIS AODs. However, the comparison with AERONET AODs suggests that AODs used for computations over desert sites are 0.08 larger. The cause of positive biases in downward shortwave irradiance and in AOD for the desert sites is possibly due to the dust particle size and distribution, as defined by the MATCH transport model, and the transfer of that information into the radiative transfer model. [ABSTRACT FROM AUTHOR]

Spectral and broadband radiances and irradiances (fluxes) were measured from surface, airborne, and spaceborne platforms in the Chesapeake Lighthouse and Aircraft Measurements for Satellites (CLAMS) campaign. The radiation data obtained on the 4 clear days over ocean during CLAMS are analyzed here with the Coupled Ocean Atmosphere Radiative Transfer (COART) model. The model is successively compared with observations of broadband fluxes and albedos near the ocean surface from the Clouds and the Earth's Radiant Energy System (CERES) Ocean Validation Experiment (COVE) sea platform and a low-level OV-10 aircraft, of near-surface spectral albedos from COVE and OV-10, of broadband radiances at multiple angles and inferred top-of-atmosphere (TOA) fluxes froth CERES, and of spectral radiances at multiple angles from Airborne Multiangle Imaging Spectroradiometer (MISR), or 'AirMISR.' at 20-km altitude. The radiation measurements from different platforms are shown to be consistent with each other and with model results. The discrepancies between the model and observations at the surface are less than l0 W [m.sup.2] for downwelling and 2 W [m.sup.2] for upwelling fluxes. The model-observation discrepancies for shortwave ocean albedo are less than 8%; some discrepancies in spectral albedo are larger but less than 20%. The discrepancies between low-altitude aircraft and surface measurements are somewhat larger than those between the model and the surface measurements: the former are due to the effects of differences in height, aircraft pitch and roll, and the noise of spatial and temporal variations of atmospheric and oceanic properties. The discrepancy between the model and the CERES observations for the upwelling radiance is 5.9% for all angles; this is reduced to 4.9% if observations within 15[degrees] of the sun-glint angle are excluded. The measurements and model agree on the principal impacts that ocean optical properties have on upwelling radiation at low levels in the atmosphere. Wind-driven surface roughness significantly affects the upwelling radiances measured by aircraft and satellites at small sun-glint angles, especially in the near-infrared channel of MISR. Intercomparisons of various measurements and the model show that most of the radiation observations in CLAMS are robust, and that the coupled radiative transfer model used here accurately treats scattering and absorption processes in both the air and the water.

The diurnal cycle of outgoing longwave radiation (OLR) from the earth is analyzed by decomposing satellite observations into a set of empirical orthogonal functions (EOFs). The observations are from the Earth Radiation Budget Experiment (ERBE) scanning radiometer aboard the Earth Radiation Budget Satellite, which had a precessing orbit with 57[degrees] inclination. The diurnal cycles of land and ocean differ considerably. The first EOF for land accounts for 73% to 85% of the variance, whereas the first EOF for ocean accounts for only 16% to 20% of the variance, depending on season. The diurnal cycle for land is surprisingly symmetric about local noon for the first EOF, which is approximately a half-sine during day and flat at night. The second EOF describes lead-lag effects due to surface heating and cloud formation. For the ocean, the first EOF and second EOF are similar to that of land, except for spring, when the first ocean EOF is a semidiurnal cycle and the second ocean EOF is the half-sine. The first EOF for land has a daytime peak of about 50 W [m.sup.-2], whereas the first ocean EOF peaks at about 25 W [m.sup.-2]. The geographical and seasonal patterns of OLR diurnal cycle provide insights into the interaction of radiation with the atmosphere and surface and are useful for validating and upgrading circulation models.

Data from the Earth Radiation Budget Experiment scanning radiometer aboard the NOAA-9 operational meteorological satellite are used to investigate the spatial variability of outgoing longwave radiation (OLR). Daily and monthly radiation maps at 2.5 [degrees] latitude-longitude scale are used as a basis for the study. The regions of greatest variability are in the tropics and subtropics. Storm tracks such as the South Pacific convergence zone appear as regions of high OLR variability. Spatial spectra in longitude show two regimes of OLR. At large scales (wavenumbers less than 6), the spatial spectrum is flat. For wavenumbers greater than 10, the spectra decrease as wavenumber to the -3 power. The spatial spectrum of daily anomalies from the mean is a strong function of latitude and season, with interesting features. Correlations of daily anomalies from the monthly mean decrease exponentially in latitude but have a damped-wave structure in longitude. The spatial variability of the daily maps, as measured by degree variance, have 10 times the power at degree 24 than the monthly maps, but at scales between 1 and 10, the degree variance is practically the same for daily as for monthly.

Aerosol optical depths (AOD) used for the Edition 4.1 Clouds and the Earth's Radiant Energy System (CERES) Synoptic (SYN1deg) are evaluated. AODs are derived from Moderate Resolution Imaging Spectroradiometer (MODIS) observations and assimilated by an aerosol transport model (MATCH). As a consequence, clear-sky AODs closely match with those derived from MODIS instruments. AODs under all-sky conditions are larger than AODs under clear-sky conditions, which is supported by ground-based AERONET observations. When all-sky MATCH AODs are compared with Modern-Era Retrospective Analysis for Research and Applications (MERRA2) AODs, MATCH AODs are generally larger than MERRA2 AODS especially over convective regions (e.g. Amazon, central Africa, and eastern Asia). The difference is largely caused by MODIS AODs used for assimilation. Including AODs with larger retrieval uncertainty makes AODs over the convective regions larger. When AODs are used for clear-sky irradiance computations and computed downward shortwave irradiances are compared with ground based observations, the computed instantaneous irradiances are 1% to 2% larger than observed irradiances. The comparison of top-of-atmosphere clear-sky irradiances with those derived from CERES observations suggests that AODs used for surface radiation observation sites are larger by 0.01 to 0.03, which is within the uncertainty of instantaneous MODIS AODs. However, the comparison with AERONET AOD suggests AODs used for computations over desert sites are 0.08 larger. The cause of positive biases of downward shortwave irradiance and AODs for the desert sites are unknown. [ABSTRACT FROM AUTHOR]

Effects of water mass imbalance and hydrometeor transport on the enthalpy flux and water phase on diabatic heating rate in computing the regional energy and water budget of the atmosphere over ocean are investigated. Equations of energy and water budget of the atmospheric column that explicitly consider the velocity of liquid and ice cloud particles, and rain and snow are formulated by separating water variables from dry air. Differences of energy budget equations formulated in this study from those used in earlier studies are that 1) diabatic heating rate depends on water phase, 2) diabatic heating due to net condensation of nonprecipitating hydrometeors is included, and 3) hydrometeors can be advected with a different velocity from the dry-air velocity. Convergence of water vapor associated with phase change and horizontal transport of hydrometeors is to increase diabatic heating in the atmospheric column where hydrometeors are formed and exported and to reduce energy where hydrometeors are imported and evaporated. The process can improve the regional energy and water mass balance when energy data products are integrated. Effects of enthalpy transport associated with water mass transport through the surface are cooling to the atmosphere and warming to the ocean when the enthalpy is averaged over the global ocean. There is no net effect to the atmosphere and ocean columns combined. While precipitation phase changes the regional diabatic heating rate up to 15 W m−2, the dependence of the global mean value on the temperature threshold of melting snow to form rain is less than 1 W m−2. [ABSTRACT FROM AUTHOR]

Radiative heating rates computed with cloud properties derived from passive and active sensors are investigated. Zonal monthly radiative heating rate anomalies computed using both active and passive sensors show that larger variability in longwave cooling exists near the tropical tropopause and near the top of the boundary layer between ~50°N to ~50°S. Aerosol variability contributes to increases in shortwave heating rate variability. When zonal monthly mean cloud effects on the radiative heating rate computed with both active and passive sensors and those computed with passive sensor only are compared, the latter shows cooling and heating peaks corresponding to cloud top and base height ranges used for separating cloud types. The difference of these two sets of cloud radiative effect on heating rates in the middle to upper troposphere is larger than the radiative heating rate uncertainty estimated based on the difference of two active sensor radiative heating rate profile data products. In addition, radiative heating rate contribution to generation of eddy available potential energy is also investigated. Although radiation contribution to generation of eddy available potential energy averaged over a year and the entire globe is small, radiation increases the eddy available potential energy in the northern hemisphere during summer. Two key elements that longwave radiation contribute to the generation of eddy potential energy are (1) longitudinal temperature gradient in the atmosphere associated with land and ocean surface temperatures contrasts and absorption of longwave radiation emitted by the surface and (2) cooling near the cloud top of stratocumulus clouds. Key Points: Large longwave radiative heating rate variability exists near the tropical tropopause and near the top of the boundary layerPassive sensor derived radiative heating rate has a larger bias in the middle to upper troposphereRadiation significantly increases the generation of eddy potential energy in the northern hemisphere during summer [ABSTRACT FROM AUTHOR]

Tower-based data combined with high-resolution satellite products have been used to produce surface albedo at various spatial scales over land. Because tower-based albedo data are available at only a few sites, surface albedos using these combined data are spatially limited. Moreover, tower-based albedo data are not representative of highly heterogeneous regions. To produce areal-averaged and spectrally-resolved surface albedo for regions with various degrees of surface heterogeneity, we have developed a transmission-based retrieval and demonstrated its feasibility for relatively homogeneous land surfaces. Here, we demonstrate its feasibility for a highly heterogeneous coastal region. We use the atmospheric transmission measured during a 19-month period (June 2009-December 2010) by a ground-based Multi-Filter Rotating Shadowband Radiometer (MFRSR) at five wavelengths (0.415, 0.5, 0.615, 0.673 and 0.87 μm) at the Department of Energy's Atmospheric Radiation Measurement (ARM) Mobile Facility (AMF) site located on Graciosa Island. We compare the MFRSR-retrieved areal-averaged surface albedo with albedo derived from Moderate Resolution Imaging Spectroradiometer (MODIS) observations, and also a composite-based albedo. We demonstrate that these three methods produce similar spectral signatures of surface albedo; however, the MFRSR-retrieved albedo, is higher on average (≤0.04) than the MODIS-based areal-averaged surface albedo and the largest difference occurs in winter. [ABSTRACT FROM AUTHOR]

NASA's Clouds and the Earth's Radiant Energy System (CERES) project integrates CERES, Moderate Resolution Imaging Spectroradiometer (MODIS), and geostationary satellite observations to provide top-of-atmosphere (TOA) irradiances derived from broadband radiance observations by CERES instruments. It also uses snow cover and sea ice extent retrieved from microwave instruments as well as thermodynamic variables from reanalysis. In addition, these variables are used for surface and atmospheric irradiance computations. The CERES project provides TOA, surface, and atmospheric irradiances in various spatial and temporal resolutions. These data sets are for climate research and evaluation of climate models. Long-term observations are required to understand how the Earth system responds to radiative forcing. A simple model is used to estimate the time to detect trends in TOA reflected shortwave and emitted longwave irradiances. [ABSTRACT FROM AUTHOR]

The Clouds and the Earth's Radiant Energy System Synoptic (SYN1deg), edition 3, product provides climate-quality global 3-hourly 1° × 1°gridded top of atmosphere, in-atmosphere, and surface radiant fluxes. The in-atmosphere surface fluxes are computed hourly using a radiative transfer code based upon inputs from Terra and Aqua Moderate Resolution Imaging Spectroradiometer (MODIS), 3-hourly geostationary (GEO) data, and meteorological assimilation data from the Goddard Earth Observing System. The GEO visible and infrared imager calibration is tied to MODIS to ensure uniform MODIS-like cloud properties across all satellite cloud datasets. Computed surface radiant fluxes are compared to surface observations at 85 globally distributed land (37) and ocean buoy (48) sites as well as several other publicly available global surface radiant flux data products. Computed monthly mean downward fluxes from SYN1deg have a bias (standard deviation) of 3.0 W m−2 (5.7%) for shortwave and −4.0 W m−2 (2.9%) for longwave compared to surface observations. The standard deviation between surface downward shortwave flux calculations and observations at the 3-hourly time scale is reduced when the diurnal cycle of cloud changes is explicitly accounted for. The improvement is smaller for surface downward longwave flux owing to an additional sensitivity to boundary layer temperature/humidity, which has a weaker diurnal cycle compared to clouds. [ABSTRACT FROM AUTHOR]

The energy budgets over land and oceans are still afflicted with considerable uncertainties, despite their key importance for terrestrial and maritime climates. We evaluate these budgets as represented in 43 CMIP5 climate models with direct observations from both surface and space and identify substantial biases, particularly in the surface fluxes of downward solar and thermal radiation. These flux biases in the various models are then linearly related to their respective land and ocean means to infer best estimates for present day downward solar and thermal radiation over land and oceans. Over land, where most direct observations are available to constrain the surface fluxes, we obtain 184 and 306 Wm for solar and thermal downward radiation, respectively. Over oceans, with weaker observational constraints, corresponding estimates are around 185 and 356 Wm. Considering additionally surface albedo and emissivity, we infer a surface absorbed solar and net thermal radiation of 136 and −66 Wm over land, and 170 and −53 Wm over oceans, respectively. The surface net radiation is thus estimated at 70 Wm over land and 117 Wm over oceans, which may impose additional constraints on the poorly known sensible/latent heat flux magnitudes, estimated here near 32/38 Wm over land, and 16/100 Wm over oceans. Estimated uncertainties are on the order of 10 and 5 Wm for most surface and TOA fluxes, respectively. By combining these surface budgets with satellite-determined TOA budgets we quantify the atmospheric energy budgets as residuals (including ocean to land transports), and revisit the global mean energy balance. [ABSTRACT FROM AUTHOR]

The NOAA-9 Earth Radiation Budget Experiment (ERBE) scanner measured broadband shortwave, longwave, and total radiances from February 1985 through January 1987. These scanner radiances are reprocessed using the more recent Clouds and the Earth's Radiant Energy System (CERES) unfiltering algorithm. The scene information, including cloud properties, required for reprocessing is derived using Advanced Very High Resolution Radiometer (AVHRR) data on board NOAA-9, while no imager data were used in the original ERBE unfiltering. The reprocessing increases the NOAA-9 ERBE scanner unfiltered longwave radiances by 1.4%-2.0% during daytime and 0.2%-0.3% during nighttime relative to those derived from the ERBE unfiltering algorithm. Similarly, the scanner unfiltered shortwave radiances increase by ~1% for clear ocean and land and decrease for all-sky ocean, land, and snow/ice by ~1%. The resulting NOAA-9 ERBE scanner unfiltered radiances are then compared with NOAA-9 nonscanner irradiances by integrating the ERBE scanner radiance over the nonscanner field of view. The comparison indicates that the integrated scanner radiances are larger by 0.9% for shortwave and 0.7% smaller for longwave. A sensitivity study shows that the one-standard-deviation uncertainties in the agreement are ±2.5%, ±1.2%, and ±1.8% for the shortwave, nighttime longwave, and daytime longwave irradiances, respectively. The NOAA-9 and ERBS nonscanner irradiances are also compared using 2 years of data. The comparison indicates that the NOAA-9 nonscanner shortwave, nighttime longwave, and daytime longwave irradiances are 0.3% larger, 0.6% smaller, and 0.4% larger, respectively. The longer observational record provided by the ERBS nonscanner plays a critical role in tying the CERES-like NOAA-9 ERBE scanner dataset from the mid-1980s to the present-day CERES scanner data record. [ABSTRACT FROM AUTHOR]

NASA's Clouds and the Earth's Radiant Energy System (CERES) project is responsible for operation and data processing of observations from scanning radiometers on board the Tropical Rainfall Measuring Mission (TRMM), Terra, Aqua, and Suomi National Polar-Orbiting Partnership (NPP) sateliites. The clouds and radiative swath (CRS) CERES data product contains irradiances computed using a radiative transfer model for nearly all CERES footprints in addition to top-of-atmosphere (TOA) irradiances derived from observed radiances by CERES instruments. This paper describes a method to constrain computed irradiances by CERES-derived TOA irradiances using Lagrangian multipliers. Radiative transfer model inputs include profiles of atmospheric temperature, humidity, aerosols and ozone, surface temperature and albedo, and up to two sets of cloud properties for a CERES footprint. Those inputs are adjusted depending on predefined uncertainties to match computed TOA and CERES-derived TOA irradiance. Because CERES instantaneous irradiances for an individual footprint also include uncertainties, primarily due to the conversion of radiance to irradiance using anisotropic directional models, the degree of the constraint depends on CERES-derived TOA irradiance as well. As a result of adjustment, TOA computed-minus-observed standard deviations are reduced from 8 to 4 W m-2 for longwave irradiance and from 15 to 6 W m-2 for shortwave irradiance. While agreement of computed TOA with CERES-derived irradiances improves, comparisons with surface observations show that model constrainment to the TOA does not reduce computation bias error at the surface. After constrainment, shortwave down at the surface has an increased bias (standard deviation) of 1% (0.5%) and longwave increases by 0.2% (0.1%). Clear-sky changes are negligible. [ABSTRACT FROM AUTHOR]

Identification of clear-sky snow and ice is an important step in the production of cryosphere radiation budget products, which are used in the derivation of long-term data series for climate research. In this paper, a new method of clear-sky snow/ice identification for Moderate Resolution Imaging Spectroradiometer (MODIS) is presented. The algorithm's goal is to enhance the identification of snow and ice within the Clouds and the Earth's Radiant Energy System (CERES) data after application of the standard CERES scene identification scheme. The input of the algorithm uses spectral radiances from five MODIS bands and surface skin temperature available in the CERES Single Scanner Footprint (SSF) product. The algorithm produces a cryosphere rating from an aggregated test: a higher rating corresponds to a more certain identification of the clear-sky snow/ice-covered scene. Empirical analysis of regions of interest representing distinctive targets such as snow, ice, ice and water clouds, open waters, and snow-free land selected from a number of MODIS images shows that the cryosphere rating of snow/ice targets falls into 95% confidence intervals lying above the same confidence intervals of all other targets. This enables recognition of clear-sky cryosphere by using a single threshold applied to the rating, which makes this technique different from traditional branching techniques based on multiple thresholds. Limited tests show that the established threshold clearly separates the cryosphere rating values computed for the cryosphere from those computed for noncryosphere scenes, whereas individual tests applied consequently cannot reliably identify the cryosphere for complex scenes. [ABSTRACT FROM AUTHOR]

Spectral and broadband radiances and irradiances (fluxes) were measured from surface, airborne, and spaceborne platforms in the Chesapeake Lighthouse and Aircraft Measurements for Satellites (CLAMS) campaign. The radiation data obtained on the 4 clear days over ocean during CLAMS are analyzed here with the Coupled Ocean–Atmosphere Radiative Transfer (COART) model. The model is successively compared with observations of broadband fluxes and albedos near the ocean surface from the Clouds and the Earth's Radiant Energy System (CERES) Ocean Validation Experiment (COVE) sea platform and a low-level OV-10 aircraft, of near-surface spectral albedos from COVE and OV-10, of broadband radiances at multiple angles and inferred top-of-atmosphere (TOA) fluxes from CERES, and of spectral radiances at multiple angles from Airborne Multiangle Imaging Spectroradiometer (MISR), or “AirMISR,” at 20-km altidude. The radiation measurements from different platforms are shown to be consistent with each other and with model results. The discrepancies between the model and observations at the surface are less than 10 W m-2 for downwelling and 2 W m-2 for upwelling fluxes. The model–observation discrepancies for shortwave ocean albedo are less than 8%; some discrepancies in spectral albedo are larger but less than 20%. The discrepancies between low-altitude aircraft and surface measurements are somewhat larger than those between the model and the surface measurements; the former are due to the effects of differences in height, aircraft pitch and roll, and the noise of spatial and temporal variations of atmospheric and oceanic properties. The discrepancy between the model and the CERES observations for the upwelling radiance is 5.9% for all angles; this is reduced to 4.9% if observations within 15° of the sun-glint angle are excluded. The measurements and model agree on the principal impacts that ocean optical properties have on upwelling radiation at low levels in the atmosphere. Wind-driven surface roughness significantly affects the upwelling radiances measured by aircraft and satellites at small sun-glint angles, especially in the near-infrared channel of MISR. Intercomparisons of various measurements and the model show that most of the radiation observations in CLAMS are robust, and that the coupled radiative transfer model used here accurately treats scattering and absorption processes in both the air and the water. [ABSTRACT FROM AUTHOR]

Investigates the spatial variability of outgoing longwave radiation using data from the Earth Radiation Budget Experiment scanning radiometer aboard the NOAA-9 operational meteorological satellite. Variances; Longitudinal variability; Correlations; Degree variances.

The Clouds and the Earth's Radiant Energy System (CERES) Energy Balanced and Filled (EBAF) Edition 4.1 data product provides global surface irradiances. Uncertainties in the global and regional monthly and annual mean all-sky net shortwave, longwave, and shortwave plus longwave (total) irradiances are estimated using ground-based observations. Error covariance is derived from surface irradiance sensitivity to surface, atmospheric, cloud and aerosol property perturbations. Uncertainties in global annual mean net shortwave, longwave, and total irradiances at the surface are, respectively, 5.7 Wm−2, 6.7 Wm−2, and 9.7 Wm−2. In addition, the uncertainty in surface downward irradiance monthly anomalies and their trends are estimated based on the difference derived from EBAF surface irradiances and observations. The uncertainty in the decadal trend suggests that when differences of decadal global mean downward shortwave and longwave irradiances are, respectively, greater than 0.45 Wm−2 and 0.52 Wm−2, the difference is larger than 1σ uncertainties. However, surface irradiance observation sites are located predominately over tropical oceans and the northern hemisphere mid-latitude. As a consequence, the effect of a discontinuity introduced by using multiple geostationary satellites in deriving cloud properties is likely to be excluded from these trend and decadal change uncertainty estimates. Nevertheless, the monthly anomaly timeseries of radiative cooling in the atmosphere (multiplied by −1) agrees reasonably well with the anomaly time series of diabatic heating derived from global mean precipitation and sensible heat flux with a correlation coefficient of 0.46. [ABSTRACT FROM AUTHOR]