Cell lineage consists of cell division timing, cell migration and cell fate, and is highly conserved during development of nematode species. An outstanding question is how differentiated cells are genetically and physically regulated in order to migrate to their precise destination among individuals. Here, we first generated a reference embryo using time-lapse 3 dimensional images of 222 wild-type C. elegans embryos at about 1.5-minute interval. This was achieved by automatic tracing and quantitative analysis of cellular phenotypes from 4- to 24-cell stage, including cell cycle duration, division orientation and migration trajectory. We next characterized cell division timing and cell kinematic state, which suggests that eight groups of cells can be clustered based on invariant and distinct division sequence. Cells may still be moving while others start to divide, indicating strong robustness against motional noise in developing embryo. We then devised a system-level phenotyping method for detecting mutant defect in global growth rate, cell cycle duration, division orientation and cell arrangement. A total of 758 genes were selected for perturbation by RNA interference followed by automatic phenotyping, which suggests a cryptic genetic architecture coordinating early morphogenesis spatially and temporally. The high-quality wild-type reference supports a conceptual close-packing model for cell arrangement during 4- to 8-cell stage, implying fundamental mechanical laws regulating the topological structure of early C. elegans embryo. Also, we observed a series of remarkable morphogenesis phenomena such as induced defect or recovery from defect in mutant embryo. To facilitate use of this quantification system, we built a software named STAR 1.0 for visualizing the wild-type reference and mutant phenotype. It also allows automatic phenotyping of new mutant embryo. Taken together, we not only provide a statistical wild-type reference with defined variability, but also shed light on both genetic and physical mechanisms coordinating early embryonic morphogenesis of C. elegans. The statistical reference permits a sensitive approach for mutant phenotype analysis, with which we phenotype a total of 1818 mutant embryos by depletion of 758 genes.