Methods and systems are disclosed to realistically simulate facial wrinkle aging of a person using a neutral state (natural look) image and one or more expression images (e.g., smile, frown, pout, wink) that induce wrinkles. The neutral and expression images are processed to simulate wrinkle aging by registering the wrinkles that are visible in the expression image onto the neutral image, thereby generating a wrinkle-aged simulated image. Advantageously, a person's own wrinkle histological data is utilized, hence providing an accurate and realistic wrinkle aging simulation. Similarly, the neutral image is processed to eliminate all visible wrinkles thereby generating a wrinkle de-aged simulation image. Additionally, blending of a neutral image with an aged or de-aged simulation image is disclosed, where the degree of blending is based on statistical modeling of skin condition with age and/or expected outcome of a particular type of treatment. The methods and systems disclosed have wide applicability, including, for example, areas such as dermatology, cosmetics and computer animation, among others.