Shape optimization of airfoil in WIG craft has been performed by considering the ground effect. The WIG craft should satisfy various aerodynamic characteristics such as lift, lift to drag ratio, and static height stability. However, they show a strong trade-off phenomenon so that it is difficult to satisfy aerodynamic properties simultaneously. Optimization is carried out through the multi-objective genetic algorithm. A multi-objective optimization means that each objective is considered separately instead of weighting. Due to the trade-off, pareto sets and non-dominated solutions can be obtained instead of the unique solution. NACA0015 airfoil is considered as a baseline model, shapes of airfoil are parameterized and rebuilt with four-Bezier curves. There are eighteen design variables and three objective functions. The range of design variables and their resolutions are two primary keys for the successful optimization. By two preliminary optimizations, the variation can be reduced effectively. After thirty evolutions, the non-dominated pareto individuals of twenty seven are obtained. Pareto sets are all the set of possible and excellent solution across the design space. At any selections of the pareto set, these are no better solutions in all design space.