Precise orbit determination and prediction is essential for navigation satellites to provide a precise navigation result to users. Recently, orbit determination accuracy of GPS and Galileo satellites has reached to centimeter level. Japanese QZSS (Quasi Zenith Satellite System) also tries to provide centimeter level orbit determination result, and JAXA (Japan aerospace exploration agency) has developed precise orbit and clock determination tool called MADOCA (Multi-GNSS Advanced Demonstration tool for Orbit and Clock Analysis). To determine a precise orbit, disturbances acting on the navigation satellites have to be accurately calculated. Gravitational disturbances (e.g., geo- potential, third body, ocean tide, and so on) have been accurately modeled by using many early studies since the behaviors of these disturbances are similar for all satellites. On the other hands, non- gravitational disturbances (e.g., solar radiation pressure and thermal re-radiation pressure) have large uncertainties since their behavior depends on position, attitude, surface geometry, and optical properties of each satellite. This research focuses on accurate non-gravitational modeling to improve orbit determination accuracy.