A 19-Month Record of Marine Aerosol–Cloud–Radiation Properties Derived from DOE ARM Mobile Facility Deployment at the Azores. Part I: Cloud Fraction and Single-Layered MBL Cloud Properties

2014 ◽  
Vol 27 (10) ◽  
pp. 3665-3682 ◽  
Author(s):  
Xiquan Dong ◽  
Baike Xi ◽  
Aaron Kennedy ◽  
Patrick Minnis ◽  
Robert Wood
Keyword(s):  
Boundary Layer ◽  
High Pressure ◽  
Liquid Water ◽  
Mixed Boundary ◽  
Cloud Condensation Nuclei ◽  
Low Pressure ◽  
Cloud Fraction ◽  
Cloud Base ◽  
Microphysical Properties ◽  
Cloud Properties

Abstract A 19-month record of total and single-layered low (<3 km), middle (3–6 km), and high (>6 km) cloud fractions (CFs) and the single-layered marine boundary layer (MBL) cloud macrophysical and microphysical properties was generated from ground-based measurements at the Atmospheric Radiation Measurement Program (ARM) Azores site between June 2009 and December 2010. This is the most comprehensive dataset of marine cloud fraction and MBL cloud properties. The annual means of total CF and single-layered low, middle, and high CFs derived from ARM radar and lidar observations are 0.702, 0.271, 0.01, and 0.106, respectively. Greater total and single-layered high (>6 km) CFs occurred during the winter, whereas single-layered low (<3 km) CFs were more prominent during summer. Diurnal cycles for both total and low CFs were stronger during summer than during winter. The CFs are bimodally distributed in the vertical with a lower peak at ~1 km and a higher peak between 8 and 11 km during all seasons, except summer when only the low peak occurs. Persistent high pressure and dry conditions produce more single-layered MBL clouds and fewer total clouds during summer, whereas the low pressure and moist air masses during winter generate more total and multilayered clouds, and deep frontal clouds associated with midlatitude cyclones. The seasonal variations of cloud heights and thickness are also associated with the seasonal synoptic patterns. The MBL cloud layer is low, warm, and thin with large liquid water path (LWP) and liquid water content (LWC) during summer, whereas during winter it is higher, colder, and thicker with reduced LWP and LWC. The cloud LWP and LWC values are greater at night than during daytime. The monthly mean daytime cloud droplet effective radius re values are nearly constant, while the daytime droplet number concentration Nd basically follows the LWC variation. There is a strong correlation between cloud condensation nuclei (CCN) concentration NCCN and Nd during January–May, probably due to the frequent low pressure systems because upward motion brings more surface CCN to cloud base (well-mixed boundary layer). During summer and autumn, the correlation between Nd and NCCN is not as strong as that during January–May because downward motion from high pressure systems is predominant. Compared to the compiled aircraft in situ measurements during the Atlantic Stratocumulus Transition Experiment (ASTEX), the cloud microphysical retrievals in this study agree well with historical aircraft data. Different air mass sources over the ARM Azores site have significant impacts on the cloud microphysical properties and surface CCN as demonstrated by great variability in NCCN and cloud microphysical properties during some months.

scholarly journals Observed correlations between aerosol and cloud properties in an Indian Ocean trade cumulus regime

2015 ◽  
Vol 15 (20) ◽  
pp. 29347-29402 ◽  
Author(s):  
K. Pistone ◽  
P. S. Praveen ◽  
R. M. Thomas ◽  
V. Ramanathan ◽  
E. Wilcox ◽  
...  
Keyword(s):  
Boundary Layer ◽  
Water Vapor ◽  
Indian Ocean ◽  
Liquid Water ◽  
Radiative Forcing ◽  
Cloud Base ◽  
Future Research ◽  
Air Mass ◽  
Vapor Flux ◽  
Cloud Properties

Abstract. There are many contributing factors which determine the micro- and macrophysical properties of clouds, including atmospheric structure, dominant meteorological conditions, and aerosol concentration, all of which may be coupled to one another. In the quest to determine aerosol effects on clouds, these potential relationships must be understood, as changes in atmospheric conditions due to aerosol may change the expected magnitude of indirect effects by altering cloud properties in unexpected ways. Here we describe several observed correlations between aerosol conditions and cloud and atmospheric properties in the Indian Ocean winter monsoon season. In the CARDEX (Cloud, Aerosol, Radiative forcing, Dynamics EXperiment) field campaign conducted in February and March 2012 in the northern Indian Ocean, continuous measurements of atmospheric precipitable water vapor and the liquid water path (LWP) of trade cumulus clouds were made, concurrent with measurements of water vapor flux, cloud and aerosol vertical profiles, meteorological data, and surface and total-column aerosol. Here we present evidence of a positive correlation between aerosol and cloud LWP which becomes clear after the data are filtered to control for the natural meteorological variability in the region. We then use the aircraft and ground observatory measurements to explore the mechanisms behind the observed aerosol–LWP correlation. We determine that increased boundary-layer humidity lowering the cloud base is responsible for the observed increase in cloud liquid water. Large-scale analysis indicates that high pollution cases originate with a highly-polluted boundary layer air mass approaching the observatory from a northwesterly direction. This polluted mass exhibits higher temperatures and humidity than the clean case, the former of which may be attributable to heating due to aerosol absorption of solar radiation over the subcontinent. While high temperature conditions dispersed along with the high-aerosol anomaly, the high humidity condition was observed to instead develop along with the polluted air mass. We then explore potential causal mechanisms of the observed correlations, though future research will be needed to more fully describe the aerosol–humidity relationship.

scholarly journals Simultaneous and synergistic profiling of cloud and drizzle properties using ground-based observations

2017 ◽  
Author(s):  
Stephanie P. Rusli ◽  
David P. Donovan ◽  
Herman W. J. Russchenberg
Keyword(s):  
Liquid Water ◽  
Optimal Estimation ◽  
Liquid Water Content ◽  
Cloud Base ◽  
Unified Framework ◽  
Forward Models ◽  
Eddy Simulation ◽  
Brightness Temperatures ◽  
Microphysical Properties ◽  
Cloud Properties

Abstract. Despite the importance of radar reflectivity (Z) measurements in the retrieval of liquid water cloud properties, it remains non-trivial to interpret Z due to the possible presence of drizzle droplets within the clouds. So far, there has been no published work that utilizes Z to identify the presence of drizzle above the cloud base in an optimized and a physically-consistent manner. In this work, we develop a retrieval technique that exploits the synergy of different remote sensing systems to carry out this task and to subsequently profile the microphysical properties of the cloud and drizzle in a unified framework. This is accomplished by using ground-based measurements of Z, lidar attenuated backscatter below as well as above the cloud base, and microwave brightness temperatures. Fast physical forward models coupled to cloud and drizzle structure parametrization are used in an optimal estimation type framework in order to retrieve the best-estimate for the cloud and drizzle property profiles. The cloud retrieval is first evaluated using synthetic signals generated from large-eddy simulation output to verify the forward models used in the retrieval procedure and the vertical parametrization of the liquid water content. From this exercise it is found that, on average, the cloud properties can be retrieved within 5 % of the mean truth. The full cloud-drizzle retrieval method is then applied to a selected ACCEPT campaign dataset collected in Cabauw, The Netherlands. An assessment of the retrieval products is performed using three independent methods from the literature, each was specifically developed to retrieve only the cloud properties, the drizzle properties below the cloud base, or the drizzle fraction within the cloud, respectively. One-to-one comparisons, taking into account the uncertainties or limitations of each retrieval, show that our results are generally consistent with what is derived using the three independent methods.

scholarly journals Simultaneous and synergistic profiling of cloud and drizzle properties using ground-based observations

2017 ◽  
Vol 10 (12) ◽  
pp. 4777-4803 ◽  
Author(s):  
Stephanie P. Rusli ◽  
David P. Donovan ◽  
Herman W. J. Russchenberg
Keyword(s):  
Liquid Water ◽  
Optimal Estimation ◽  
Liquid Water Content ◽  
Cloud Base ◽  
Unified Framework ◽  
Forward Models ◽  
Eddy Simulation ◽  
Brightness Temperatures ◽  
Microphysical Properties ◽  
Cloud Properties

Abstract. Despite the importance of radar reflectivity (Z) measurements in the retrieval of liquid water cloud properties, it remains nontrivial to interpret Z due to the possible presence of drizzle droplets within the clouds. So far, there has been no published work that utilizes Z to identify the presence of drizzle above the cloud base in an optimized and a physically consistent manner. In this work, we develop a retrieval technique that exploits the synergy of different remote sensing systems to carry out this task and to subsequently profile the microphysical properties of the cloud and drizzle in a unified framework. This is accomplished by using ground-based measurements of Z, lidar attenuated backscatter below as well as above the cloud base, and microwave brightness temperatures. Fast physical forward models coupled to cloud and drizzle structure parameterization are used in an optimal-estimation-type framework in order to retrieve the best estimate for the cloud and drizzle property profiles. The cloud retrieval is first evaluated using synthetic signals generated from large-eddy simulation (LES) output to verify the forward models used in the retrieval procedure and the vertical parameterization of the liquid water content (LWC). From this exercise it is found that, on average, the cloud properties can be retrieved within 5 % of the mean truth. The full cloud–drizzle retrieval method is then applied to a selected ACCEPT (Analysis of the Composition of Clouds with Extended Polarization Techniques) campaign dataset collected in Cabauw, the Netherlands. An assessment of the retrieval products is performed using three independent methods from the literature; each was specifically developed to retrieve only the cloud properties, the drizzle properties below the cloud base, or the drizzle fraction within the cloud. One-to-one comparisons, taking into account the uncertainties or limitations of each retrieval, show that our results are consistent with what is derived using the three independent methods.

scholarly journals Arctic Mixed-Phase Stratiform Cloud Properties from Multiple Years of Surface-Based Measurements at Two High-Latitude Locations

2009 ◽  
Vol 66 (9) ◽  
pp. 2874-2887 ◽  
Author(s):  
Gijs de Boer ◽  
Edwin W. Eloranta ◽  
Matthew D. Shupe
Keyword(s):  
Liquid Water ◽  
Microwave Radiometer ◽  
Single Layer ◽  
Mixed Phase ◽  
Cloud Base ◽  
Microphysical Properties ◽  
Cloud Properties ◽  
Arctic Change ◽  
Mixed Phase Clouds ◽  
Water Contents

Abstract Macro- and microphysical properties of single-layer stratiform mixed-phase clouds are derived from multiple years of lidar, radar, and radiosonde observations. Measurements were made as part of the Mixed-Phase Arctic Clouds Experiment (MPACE) and the Study of Environmental Arctic Change (SEARCH) in Barrow, Alaska, and Eureka, Nunavut, Canada, respectively. Single-layer mixed-phase clouds occurred between 4% and 26% of the total time observed, varying with season and location. They had mean cloud-base heights between ∼700 and 2100 m and thicknesses between ∼200 and 700 m. Seasonal mean cloud optical depths ranged from 2.2 up. The clouds existed at temperatures of ∼242–271 K and occurred under different wind conditions, depending on season. Utilizing retrievals from a combination of lidar, radar, and microwave radiometer, mean cloud microphysical properties were derived, with mean liquid effective diameters estimated from 16 to 49 μm, mean liquid number densities on the order of 104–105 L−1, and mean water contents estimated between 0.07 and 0.28 g m−3. Ice precipitation was shown to have mean ice effective diameters of 50–125 μm, mean ice number densities on the order of 10 L−1, and mean water contents estimated between 0.012 and 0.031 g m−3. Mean cloud liquid water paths ranged from 25 to 100 g m−2. All results are compared to previous studies, and potential retrieval errors are discussed. Additionally, seasonal variation in macro- and microphysical properties was highlighted. Finally, fraction of liquid water to ice mass was shown to decrease with decreasing temperature.

scholarly journals A comparison of ship and satellite measurements of cloud properties with global climate model simulations in the southeast Pacific stratus deck

2010 ◽  
Vol 10 (14) ◽  
pp. 6527-6536 ◽  
Author(s):  
M. A. Brunke ◽  
S. P. de Szoeke ◽  
P. Zuidema ◽  
X. Zeng
Keyword(s):  
Diurnal Cycle ◽  
Climate Model ◽  
Global Climate ◽  
Global Climate Model ◽  
Cloud Fraction ◽  
Cloud Base ◽  
Liquid Water Path ◽  
Cloud Properties ◽  
Southeast Pacific ◽  
Cloud Thickness

Abstract. Here, liquid water path (LWP), cloud fraction, cloud top height, and cloud base height retrieved by a suite of A-train satellite instruments (the CPR aboard CloudSat, CALIOP aboard CALIPSO, and MODIS aboard Aqua) are compared to ship observations from research cruises made in 2001 and 2003–2007 into the stratus/stratocumulus deck over the southeast Pacific Ocean. It is found that CloudSat radar-only LWP is generally too high over this region and the CloudSat/CALIPSO cloud bases are too low. This results in a relationship (LWP~h9) between CloudSat LWP and CALIPSO cloud thickness (h) that is very different from the adiabatic relationship (LWP~h2) from in situ observations. Such biases can be reduced if LWPs suspected to be contaminated by precipitation are eliminated, as determined by the maximum radar reflectivity Zmax>−15 dBZ in the apparent lower half of the cloud, and if cloud bases are determined based upon the adiabatically-determined cloud thickness (h~LWP1/2). Furthermore, comparing results from a global model (CAM3.1) to ship observations reveals that, while the simulated LWP is quite reasonable, the model cloud is too thick and too low, allowing the model to have LWPs that are almost independent of h. This model can also obtain a reasonable diurnal cycle in LWP and cloud fraction at a location roughly in the centre of this region (20° S, 85° W) but has an opposite diurnal cycle to those observed aboard ship at a location closer to the coast (20° S, 75° W). The diurnal cycle at the latter location is slightly improved in the newest version of the model (CAM4). However, the simulated clouds remain too thick and too low, as cloud bases are usually at or near the surface.

scholarly journals Investigation of Regional and Seasonal Variations in Marine Boundary Layer Cloud Properties from MODIS Observations

2008 ◽  
Vol 21 (19) ◽  
pp. 4955-4973 ◽  
Author(s):  
Michael P. Jensen ◽  
Andrew M. Vogelmann ◽  
William D. Collins ◽  
Guang J. Zhang ◽  
Edward P. Luke
Keyword(s):  
Boundary Layer ◽  
General Circulation ◽  
Climate Models ◽  
Marine Boundary Layer ◽  
General Circulation Models ◽  
Test Bed ◽  
Liquid Water Path ◽  
Microphysical Properties ◽  
Cloud Properties ◽  
Peak Versus

Abstract To aid in understanding the role that marine boundary layer (MBL) clouds play in climate and assist in improving their representations in general circulation models (GCMs), their long-term microphysical and macroscale characteristics are quantified using observations from the Moderate Resolution Imaging Spectroradiometer (MODIS) instrument aboard the National Aeronautics and Space Administration’s (NASA’s) Terra satellite. Six years of MODIS pixel-level cloud products are used from oceanic study regions off the west coasts of California, Peru, the Canary Islands, Angola, and Australia where these cloud types are common. Characterizations are given for their organization (macroscale structure), the associated microphysical properties, and the seasonal dependencies of their variations for scales consistent with the size of a GCM grid box (300 km × 300 km). MBL mesoscale structure is quantified using effective cloud diameter CD, which is introduced here as a simplified measure of bulk cloud organization; it is straightforward to compute and provides descriptive information beyond that offered by cloud fraction. The interrelationships of these characteristics are explored while considering the influences of the MBL state, such as the occurrence of drizzle. Several commonalities emerge for the five study regions. MBL clouds contain the best natural examples of plane-parallel clouds, but overcast clouds occur in only about 25% of the scenes, which emphasizes the importance of representing broken MBL cloud fields in climate models (that are subgrid scale). During the peak months of cloud occurrence, mesoscale organization (larger CD) increases such that the fractions of scenes characterized as “overcast” and “clumped” increase at the expense of the “scattered” scenes. Cloud liquid water path and visible optical depth usually trend strongly with CD, with the largest values occurring for scenes that are drizzling. However, considerable interregional differences exist in these trends, suggesting that different regression functionalities exist for each region. For peak versus off-peak months, the fraction of drizzling scenes (as a function of CD) are similar for California and Angola, which suggests that a single probability distribution function might be used for their drizzle occurrence in climate models. The patterns are strikingly opposite for Peru and Australia; thus, the contrasts among regions may offer a test bed for model simulations of MBL drizzle occurrence.

Reconciling and Validating the Cloud Thickness and Liquid Water Path Tendencies Proposed by R. Wood and J. J. van der Dussen et al.

2015 ◽  
Vol 72 (5) ◽  
pp. 2033-2040 ◽  
Author(s):  
Mohamed S. Ghonima ◽  
Joel R. Norris ◽  
Thijs Heus ◽  
Jan Kleissl
Keyword(s):  
Boundary Layer ◽  
Liquid Water ◽  
Mixed Boundary ◽  
Liquid Water Path ◽  
Water Path ◽  
Large Eddy ◽  
Heat And Moisture ◽  
Cloud Thickness ◽  
Detailed Derivation ◽  
Good Agreement

Abstract A detailed derivation of stratocumulus cloud thickness and liquid water path tendencies as a function of the well-mixed boundary layer mass, heat, and moisture budget equations is presented. The derivation corrects an error in the cloud thickness tendency equation derived by R. Wood to make it consistent with the liquid water path tendency equation derived by J. J. van der Dussen et al. The validity of the tendency equations is then tested against the output of large-eddy simulations of a typical stratocumulus-topped boundary layer case and is found to be in good agreement.

scholarly journals Climatology of stratocumulus cloud morphologies: microphysical properties and radiative effects

2014 ◽  
Vol 14 (13) ◽  
pp. 6695-6716 ◽  
Author(s):  
A. Muhlbauer ◽  
I. L. McCoy ◽  
R. Wood
Keyword(s):  
Physical Properties ◽  
Global Scale ◽  
Cloud Fraction ◽  
Cloud Base ◽  
Trade Wind ◽  
Radiative Effects ◽  
Regional Variability ◽  
Microphysical Properties ◽  
Low Clouds ◽  
Low Cloud

Abstract. An artificial neural network cloud classification scheme is combined with A-train observations to characterize the physical properties and radiative effects of marine low clouds based on their morphology and type of mesoscale cellular convection (MCC) on a global scale. The cloud morphological categories are (i) organized closed MCC, (ii) organized open MCC and (iii) cellular but disorganized MCC. Global distributions of the frequency of occurrence of MCC types show clear regional signatures. Organized closed and open MCCs are most frequently found in subtropical regions and in midlatitude storm tracks of both hemispheres. Cellular but disorganized MCC are the predominant type of marine low clouds in regions with warmer sea surface temperature such as in the tropics and trade wind zones. All MCC types exhibit a pronounced seasonal cycle. The physical properties of MCCs such as cloud fraction, radar reflectivity, drizzle rates and cloud top heights as well as the radiative effects of MCCs are found highly variable and a function of the type of MCC. On a global scale, the cloud fraction is largest for closed MCC with mean cloud fractions of about 90%, whereas cloud fractions of open and cellular but disorganized MCC are only about 51% and 40%, respectively. Probability density functions (PDFs) of cloud fractions are heavily skewed and exhibit modest regional variability. PDFs of column maximum radar reflectivities and inferred cloud base drizzle rates indicate fundamental differences in the cloud and precipitation characteristics of different MCC types. Similarly, the radiative effects of MCCs differ substantially from each other in terms of shortwave reflectance and transmissivity. These differences highlight the importance of low-cloud morphologies and their associated cloudiness on the shortwave cloud forcing.

scholarly journals Version 4 CALIPSO IIR ice and liquid water cloud microphysical properties, Part I: the retrieval algorithms

2020 ◽  
Author(s):  
Anne Garnier ◽  
Jacques Pelon ◽  
Nicolas Pascal ◽  
Mark A. Vaughan ◽  
Philippe Dubuisson ◽  
...  
Keyword(s):  
Liquid Water ◽  
Cloud Base ◽  
Transfer Model ◽  
Effective Diameter ◽  
Ice Clouds ◽  
Clear Sky ◽  
Brightness Temperatures ◽  
Microphysical Properties ◽  
Data Products ◽  
Level 2

Abstract. Following the release of the Version 4 Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) data products from the Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) mission, a new version 4 (V4) of the CALIPSO Imaging Infrared Radiometer (IIR) Level 2 data products has been developed. The IIR Level 2 data products include cloud effective emissivities and cloud microphysical properties such as effective diameter and ice or liquid water path estimates. Dedicated retrievals for water clouds were added in V4, taking advantage of the high sensitivity of the IIR retrieval technique to small particle sizes. This paper (Part I) describes the improvements in the V4 algorithms compared to those used in the version 3 (V3) release, while results will be presented in a companion (Part II) paper. To reduce biases at very small emissivities that were made evident in V3, the radiative transfer model used to compute clear sky brightness temperatures over oceans has been updated and tuned for the simulations using MERRA-2 data to match IIR observations in clear sky conditions. Furthermore, the clear-sky mask has been refined compared to V3 by taking advantage of additional information now available in the V4 CALIOP 5-km layer products used as an input to the IIR algorithm. After sea surface emissivity adjustments, observed and computed brightness temperatures differ by less than ± 0.2 K at night for the three IIR channels centered at 08.65, 10.6, and 12.05 µm, and inter-channel biases are reduced from several tens of Kelvin in V3 to less than 0.1 K in V4. We have also aimed at improving retrievals in ice clouds having large optical depths by refining the determination of the radiative temperature needed for emissivity computation. The initial V3 estimate, namely the cloud centroid temperature derived from CALIOP, is corrected using a parameterized function of temperature difference between cloud base and top altitudes, cloud absorption optical depth, and the CALIOP multiple scattering correction factor. As shown in Part II, this improvement reduces the low biases at large optical depths that were seen in V3, and increases the number of retrievals in dense ice clouds. As in V3, the IIR microphysical retrievals use the concept of microphysical indices applied to the pairs of IIR channels at 12.05 μm and 10.6 μm and at 12.05 μm and 08.65 μm. The V4 algorithm uses ice look-up tables (LUTs) built using two ice crystal models from the recent TAMUice 2016 database, namely the single hexagonal column model and the 8-element column aggregate model, from which bulk properties are synthesized using a gamma size distribution. Four sets of effective diameters derived from a second approach are also reported in V4. Here, the LUTs are analytical functions relating microphysical index applied to IIR channels 12.05 µm and 10.6 µm and effective diameter as derived from in situ measurements at tropical and mid-latitudes during the TC4 and SPARTICUS field experiments.

scholarly journals Vertical Velocity Statistics in Fair-Weather Cumuli at the ARM TWP Nauru Climate Research Facility

2010 ◽  
Vol 23 (24) ◽  
pp. 6590-6604 ◽  
Author(s):  
Pavlos Kollias ◽  
Bruce Albrecht
Keyword(s):  
Boundary Layer ◽  
Numerical Models ◽  
Radar Data ◽  
Cloud Fraction ◽  
Cloud Base ◽  
Fair Weather ◽  
Research Facility ◽  
Climate Research ◽  
Velocity Statistics ◽  
Air Motion

Abstract Fair-weather cumuli are fundamental in regulating the vertical structure of water vapor and entropy in the lowest 2–3 km of the earth’s atmosphere over vast areas of the oceans. In this study, a long record of profiling cloud radar observations at the Atmospheric Radiation Measurement Program (ARM) Climate Research Facility (ACRF) at Nauru Island is used to investigate cloud vertical air motion statistics over an 8-yr observing period. Appropriate processing of the observed low radar reflectivities provides radar volume samples that contain only small cloud droplets; thus, the Doppler velocities are used as air motion tracers. The technique is applied to shallow boundary layer clouds (less than 1000 m thick) during the 1999–2007 period when radar data are available. Using the boundary layer winds from the soundings obtained at the Nauru ACRF, the fair-weather cumuli fields are classified in easterly and westerly boundary layer wind regimes. This distinction is necessary to separate marine-forced (westerlies) from land-forced (easterlies) shallow clouds because of a well-studied island effect at the Nauru ACRF. The two regimes exhibit large diurnal differences in cloud fraction and cloud dynamics as manifested by the analysis of the hourly averaged vertical air motion statistics. The fair-weather cumuli fields associated with easterlies exhibit a strong diurnal cycle in cloud fraction and updraft strength and fraction, indicating a strong influence of land-forced clouds. In contrast over the fair-weather cumuli with oceanic origin, land-forced clouds are characterized by uniform diurnal cloudiness and persistent updrafts at the cloud-base level. This study provides a unique observational dataset appropriate for testing fair-weather cumulus mass flux and turbulence parameterizations in numerical models.

Sign in / Sign up

Export Citation Format

Share Document